It has recently been shown that ancestors of New Guineans and Bougainville Islanders have inherited a proportion of their ancestry from Denisovans, an archaic hominin group from Siberia. However, only a sparse sampling of populations from Southeast Asia and Oceania were analyzed. Here, we quantify Denisova admixture in 33 additional populations from Asia and Oceania. Aboriginal Australians, Near Oceanians, Polynesians, Fijians, east Indonesians, and Mamanwa (a "Negrito" group from the Philippines) have all inherited genetic material from Denisovans, but mainland East Asians, western Indonesians, Jehai (a Negrito group from Malaysia), and Onge (a Negrito group from the Andaman Islands) have not. These results indicate that Denisova gene flow occurred into the common ancestors of New Guineans, Australians, and Mamanwa but not into the ancestors of the Jehai and Onge and suggest that relatives of present-day East Asians were not in Southeast Asia when the Denisova gene flow occurred. Our finding that descendants of the earliest inhabitants of Southeast Asia do not all harbor Denisova admixture is inconsistent with a history in which the Denisova interbreeding occurred in mainland Asia and then spread over Southeast Asia, leading to all its earliest modern human inhabitants. Instead, the data can be most parsimoniously explained if the Denisova gene flow occurred in Southeast Asia itself. Thus, archaic Denisovans must have lived over an extraordinarily broad geographic and ecological range, from Siberia to tropical Asia.
By investigating genealogical relationships, we estimated the phylogenetic history and biogeography in the megophryid genus Leptobrachium (sensu lato, including Vibrissaphora) from southern China, Indochina, Thailand and the Sundaland. The genealogical relationships among the 30 named and unnamed taxa were estimated using 2009 bp of sequences from the mitochondrial DNA genes 12S rRNA, tRNA(val), and 16S rRNA using maximum parsimony, maximum likelihood, and Bayesian inference methods. The genus Leptobrachium was a well-supported monophyletic group that contained two major clades. One clade had three subclades primarily from disjunct regions including Borneo, Peninsular Malaysia and Java, and Thailand. The Bornean subclade included one species each from the Philippines and Sumatra. The other major clade consisted of two subclades, one from Indochina and the other from southern China (Vibrissaphora). Divergence times estimated an old evolutionary history of each subclade, one that could not be explained by the geohistory of Southeast Asian major landmasses.
The extent to which response to environmental change is mediated by species-specific ecology is an important aspect of the population histories of tropical taxa. During the Pleistocene glacial cycles and associated sea level fluctuations, the Sunda region in Southeast Asia experienced concurrent changes in landmass area and the ratio of forest to open habitat, providing an ideal setting to test the expectation that habitat associations played an important role in determining species' response to the opportunity for geographic expansion. We used mitochondrial control region sequences and six microsatellite loci to compare the phylogeographic structure and demographic histories of four broadly sympatric species of Old World fruit bats in the genus, Cynopterus. Two forest-associated species and two open-habitat generalists were sampled along a latitudinal transect in Singapore, peninsular Malaysia, and southern Thailand. Contrary to expectations based on habitat associations, the geographic scale of population structure was not concordant across ecologically similar species. We found evidence for long and relatively stable demographic history in one forest and one open-habitat species, and inferred non-coincident demographic expansions in the second forest and open-habitat species. Thus, while these results indicate that Pleistocene climate change did not have a single effect on population structure across species, a correlation between habitat association and response to environmental change was supported in only two of four species. We conclude that interactions between multiple factors, including historical and contemporary environmental change, species-specific ecology and interspecific interactions, have shaped the recent evolutionary histories of Cynopterus fruit bats in Southeast Asia.
Biogeographic connections between Australia and other continents are still poorly understood although the plate tectonics of the Indo-Pacific region is now well described. Eupetes macrocerus is an enigmatic taxon distributed in a small area on the Malay Peninsula and on Sumatra and Borneo. It has generally been associated with Ptilorrhoa in New Guinea on the other side of Wallace's Line, but a relationship with the West African Picathartes has also been suggested. Using three nuclear markers, we demonstrate that Eupetes is the sister taxon of the South African genus Chaetops, and their sister taxon in turn being Picathartes, with a divergence in the Eocene. Thus, this clade is distributed in remote corners of Africa and Asia, which makes the biogeographic history of these birds very intriguing. The most parsimonious explanation would be that they represent a relictual basal group in the Passerida clade established after a long-distance dispersal from the Australo-Papuan region to Africa. Many earlier taxonomic arrangements may have been based on assumptions about relationships with similar-looking forms in the same, or adjacent, biogeographic regions, and revisions with molecular data may uncover such cases of neglect of ancient relictual patterns reflecting past connections between the continents.
Long columns of migrating larval sciarid armyworms were discovered in central and northern Japan, specifically Kanagawa, Gunma, Miyagi and Akita prefectures, as well as Hokkaido. This is the first examination of armyworms in East Asia. In Europe, armyworms have been identified as Sciara militaris, belonging to the family Sciaridae (sciarid flies or black fungus gnats), by rearing them to adulthood. In Japan, we were unable to obtain live samples for rearing; therefore, DNA barcodes were obtained from the samples of armyworms collected in the Gunma and Miyagi prefectures. The DNA barcodes were compared with those obtained from the following samples: pupae of S. militaris from UK, adults of Sciara kitakamiensis, Sciara humeralis, Sciara hemerobioides, Sciara thoracica, Sciara helvola and Sciara melanostyla from Japan, and adults of one undescribed Sciara species from Malaysia. Neighbour-joining, maximum parsimony, and maximum likelihood analyses revealed that the armyworms discovered in Japan are S. kitakamiensis. Although adults of this species have been recorded in several locations in Japan, this is the first report of migrating larval armyworms. DNA barcodes were effectively used to link different life stages of this species. The average intraspecific and interspecific pairwise genetic distances of the genus Sciara were 0.3% and 12.6%, respectively. The present study illustrates that DNA barcodes are an effective means of identifying sciarid flies in Japan.
Nasopharyngeal carcinoma (NPC) is one of the most common cancers in Malaysia, mainly occurring among the Chinese population. To detect common genetic alterations in NPC, we screened seven cases of NPC using the comparative genomic hybridization (CGH) technique. Before proceeding to the CGH technique, the tumors were first confirmed to consist of 75% tumor cells or more. In brief, the technique consists of binding tumor DNA with normal DNA and human Cot-1 DNA, which is then hybridized to normal metaphase spreads. The slides were then counterstained with 4,6 diamino-2-phenylindole (DAPI II) for detection. Analyses were performed using CGH software (Cytovision). We found genetic alterations in all seven NPC samples. The common chromosomal gains (57%, four cases) were found on chromosome arms 1q, 4p, 5, 7q, 11, 14p, 15q, 18p, and 21p, and common chromosomal losses (43%, three cases) were found on chromosome arm 16p. Our results showed chromosomal alterations in all seven NPC cases in the Malaysian population. This result provides the platform for further investigations to locate tumor suppressor genes and oncogenes at specific chromosomal regions in Malaysian NPC patients.
Oil palm (Elaeis guineensis Jacq.) requires 19 years per cycle of phenotypic selection. The use of molecular markers may reduce the generation interval and the cost of oil-palm breeding. Our objectives were to compare, by simulation, the response to phenotypic selection, marker-assisted recurrent selection (MARS), and genomewide selection with small population sizes in oil palm, and assess the efficiency of each method in terms of years and cost per unit gain. Markers significantly associated with the trait were used to calculate the marker scores in MARS, whereas all markers were used (without significance tests) to calculate the marker scores in genomewide selection. Responses to phenotypic selection and genomewide selection were consistently greater than the response to MARS. With population sizes of N = 50 or 70, responses to genomewide selection were 4-25% larger than the corresponding responses to phenotypic selection, depending on the heritability and number of quantitative trait loci. Cost per unit gain was 26-57% lower with genomewide selection than with phenotypic selection when markers cost US $1.50 per data point, and 35-65% lower when markers cost $0.15 per data point. With population sizes of N = 50 or 70, time per unit gain was 11-23 years with genomewide selection and 14-25 years with phenotypic selection. We conclude that for a realistic yet relatively small population size of N = 50 in oil palm, genomewide selection is superior to MARS and phenotypic selection in terms of gain per unit cost and time. Our results should be generally applicable to other tree species that are characterized by long generation intervals, high costs of maintaining breeding plantations, and small population sizes in selection programs.
The major histocompatibility complex (MHC) plays a crucial role in the immune system, and in some species, it is a target by which individuals choose mates to optimize the fitness of their offspring, potentially mediated by olfactory cues. Under the genetic compatibility hypothesis, individuals are predicted to choose mates with compatible MHC alleles, to increase the fitness of their offspring. Studies of MHC-based mate choice in wild mammals are under-represented currently, and few investigate more than one class of MHC genes. We investigated mate choice based on the compatibility of MHC class I and II genes in a wild population of European badgers (Meles meles). We also investigated mate choice based on microsatellite-derived pairwise relatedness, to attempt to distinguish MHC-specific effects from genomewide effects. We found MHC-assortative mating, based on MHC class II, but not class I genes. Parent pairs had smaller MHC class II DRB amino acid distances and smaller functional distances than expected from random pairings. When we separated the analyses into within-group and neighbouring-group parent pairs, only neighbouring-group pairs showed MHC-assortative mating, due to similarity at MHC class II loci. Our randomizations showed no evidence of genomewide-based inbreeding, based on 35 microsatellite loci; MHC class II similarity was therefore the apparent target of mate choice. We propose that MHC-assortative mate choice may be a local adaptation to endemic pathogens, and this assortative mate choice may have contributed to the low MHC genetic diversity in this population.
The ratio of genetic diversity on X chromosomes relative to autosomes in organisms with XX/XY sex chromosomes could provide fundamental insight into the process of genome evolution. Here we report this ratio for 24 cynomolgus monkeys (Macaca fascicularis) originating in Indonesia, Malaysia, and the Philippines. The average X/A diversity ratios in these samples was 0.34 and 0.20 in the Indonesian-Malaysian and Philippine populations, respectively, considerably lower than the null expectation of 0.75. A Philippine population supposed to derive from an ancestral population by founding events showed a significantly lower ratio than the parental population, suggesting a demographic effect for the reduction. Taking sex-specific mutation rate bias and demographic effect into account, expected X/A diversity ratios generated by computer simulations roughly agreed with the observed data in the intergenic regions. In contrast, silent sites in genic regions on X chromosomes showed strong reduction in genetic diversity and the observed X/A diversity ratio in the genic regions cannot be explained by mutation rate bias and demography, indicating that natural selection also reduces the level of polymorphism near genes. Whole-genome analysis of a female cynomolgus monkey also supported the notion of stronger reduction of genetic diversity near genes on the X chromosome.
HIV-1 epidemics among men who have sex with men (MSM) continue to expand in developed and developing countries. Although HIV infection in MSM is amongst the highest of the key affected populations in many countries in Southeast Asia, comprehensive molecular epidemiological study of HIV-1 among MSM remains inadequate in the region including in Malaysia. Here, we reported the phylodynamic profiles of HIV-1 genotypes circulating among MSM population in Kuala Lumpur, Malaysia. A total of n = 459 newly-diagnosed treatment-naïve consenting subjects were recruited between March 2006 and August 2012, of whom 87 (18.9%) were self-reported MSM. Transmitted drug resistance mutations were absent in these isolates. Cumulatively, phylogenetic reconstructions of the pro-rt gene (HXB2∶2253-3275) showed that HIV-1 subtype B and CRF01_AE were predominant and contributed to approximately 80% of the total HIV-1 infection among MSM. In addition to numerous unique transmission lineages within these genotypes, twelve monophyletic transmission clusters of different sizes (2-7 MSM sequences, supported by posterior probability value of 1) were identified in Malaysia. Bayesian coalescent analysis estimated that the divergence times for these clusters were mainly dated between 1995 and 2005 with four major transmission clusters radiating at least 12 years ago suggesting that active spread of multiple sub-epidemic clusters occurred during this period. The changes in effective population size of subtype B showed an exponential growth within 5 years between 1988 and 1993, while CRF01_AE lineage exhibited similar expansion between 1993 and 2003. Our study provides the first insight of the phylodynamic profile of HIV-1 subtype B and CRF01_AE circulating among MSM population in Kuala Lumpur, Malaysia, unravelling the importance of understanding transmission behaviours as well as evolutionary history of HIV-1 in assessing the risk of outbreak or epidemic expansion.
Multivariate analyses were performed using 13 morphological traits and 13 molecular markers (10 SSRs and three ISSRs) to assess the phylogenetic relationship among tungro resistant genotypes. For morphological traits, the genotypes were grouped into six clusters, according to D(2) statistic and Canonical vector analysis. Plant height, days to flowering, days to maturity, panicle length, number of spikelet per panicle, number of unfilled grain per panicle and yield were important contributors to genetic divergence in 14 rice genotypes. Based on Nei's genetic distance for molecular studies, seven clusters were formed among the tungro resistant and susceptible genotypes. Mantel's test revealed a significant correlation (r = 0.834*) between the morphological and molecular data. To develop high yielding tungro resistant varieties based on both morphological and molecular analyses, crosses could be made with susceptible (BR10 and BR11) genotypes with low yielding but highly resistant genotypes, Sonahidemota, Kumragoir, Nakuchimota, Khaiyamota, Khairymota and Kachamota. The chi-square analysis for seven alleles (RM11, RM17, RM20, RM23, RM80, RM108 and RM531) of SSR and five loci (RY1, MR1, MR2, MR4 and GF5) of three ISSR markers in F2 population of cross, BR11×Sonahidemota, showed a good fit to the expected segregation ratio (1:2:1) for a single gene model.
In mycobacteria, conjugation differs from the canonical Hfr model, but is still poorly understood. Here, we quantified this evolutionary processe in a natural mycobacterial population, taking advantage of a large clinical strain collection of the emerging pathogen Mycobacterium abscessus (MAB).
Climatic differences across a taxon's range may be associated with specific bioenergetic demands and may result in genetics-based metabolic adaptation, particularly in aquatic ectothermic organisms that rely on heat exchange with the environment to regulate key physiological processes. Extending down the east coast of Australia, the Great Dividing Range (GDR) has a strong influence on climate and the evolutionary history of freshwater fish species. Despite the GDR acting as a strong contemporary barrier to fish movement, many species, and species with shared ancestries, are found on both sides of the GDR, indicative of historical dispersal events. We sequenced complete mitogenomes from the four extant species of the freshwater cod genus Maccullochella, two of which occur on the semi-arid, inland side of the GDR, and two on the mesic coastal side. We constructed a dated phylogeny and explored the relative influences of purifying and positive selection in the evolution of mitogenome divergence among species. Results supported mid- to late-Pleistocene divergence of Maccullochella across the GDR (220-710 thousand years ago), bringing forward previously reported dates. Against a background of pervasive purifying selection, we detected potentially functionally relevant fixed amino acid differences across the GDR. Although many amino acid differences between inland and coastal species may have become fixed under relaxed purifying selection in coastal environments rather than positive selection, there was evidence of episodic positive selection acting on specific codons in the Mary River coastal lineage, which has consistently experienced the warmest and least extreme climate in the genus.
Oil palm is the most productive oil-bearing crop. Although it is planted on only 5% of the total world vegetable oil acreage, palm oil accounts for 33% of vegetable oil and 45% of edible oil worldwide, but increased cultivation competes with dwindling rainforest reserves. We report the 1.8-gigabase (Gb) genome sequence of the African oil palm Elaeis guineensis, the predominant source of worldwide oil production. A total of 1.535 Gb of assembled sequence and transcriptome data from 30 tissue types were used to predict at least 34,802 genes, including oil biosynthesis genes and homologues of WRINKLED1 (WRI1), and other transcriptional regulators, which are highly expressed in the kernel. We also report the draft sequence of the South American oil palm Elaeis oleifera, which has the same number of chromosomes (2n = 32) and produces fertile interspecific hybrids with E. guineensis but seems to have diverged in the New World. Segmental duplications of chromosome arms define the palaeotetraploid origin of palm trees. The oil palm sequence enables the discovery of genes for important traits as well as somaclonal epigenetic alterations that restrict the use of clones in commercial plantings, and should therefore help to achieve sustainability for biofuels and edible oils, reducing the rainforest footprint of this tropical plantation crop.
We analyse molecular and phenotypic evolution in a group of taxonomically problematic Indomalayan pitvipers, the Trimeresurus sumatranus group. Mitochondrial DNA sequencing provides a well-resolved phylogeny, with each species representing a distinct lineage. Multivariate morphological analysis reveals a high level of phenotypic differentiation, which is congruent between the sexes but does not reflect phylogenetic history. An adaptive explanation for the observed pattern of differentiation is supported by independent contrasts analysis, which shows significant correlations between current ecology and the characters that most account for the variation between taxa, including those that are presently used to identify the species. Reduced precipitation and altitude, and increased temperature, are correlated with higher numbers of scales on the head, body and tail. It is hypothesized that scale number plays an important role in heat and water exchange by influencing the area of exposed of interstitial skin, and that colour pattern variation reflects selection pressures involving camouflage and thermoregulation. Ecological convergence in traits used for classification is found to have important implications for species identification where taxa are distributed over varying environments.
We investigate the evolution of host association in a cryptic complex of mutualistic Crematogaster (Decacrema) ants that inhabits and defends Macaranga trees in Southeast Asia. Previous phylogenetic studies based on limited samplings of Decacrema present conflicting reconstructions of the evolutionary history of the association, inferring both cospeciation and the predominance of host shifts. We use cytochrome oxidase I (COI) to reconstruct phylogenetic relationships in a comprehensive sampling of the Decacrema inhabitants of Macaranga. Using a published Macaranga phylogeny, we test whether the ants and plants have cospeciated. The COI phylogeny reveals 10 well-supported lineages and an absence of cospeciation. Host shifts, however, have been constrained by stem traits that are themselves correlated with Macaranga phylogeny. Earlier lineages of Decacrema exclusively inhabit waxy stems, a basal state in the Pachystemon clade within Macaranga, whereas younger species of Pachystemon, characterized by nonwaxy stems, are inhabited only by younger lineages of Decacrema. Despite the absence of cospeciation, the correlated succession of stem texture in both phylogenies suggests that Decacrema and Pachystemon have diversified in association, or codiversified. Subsequent to the colonization of the Pachystemon clade, Decacrema expanded onto a second clade within Macaranga, inducing the development of myrmecophytism in the Pruinosae group. Confinement to the aseasonal wet climate zone of western Malesia suggests myrmecophytic Macaranga are no older than the wet forest community in Southeast Asia, estimated to be about 20 million years old (early Miocene). Our calculation of COI divergence rates from several published arthropod studies that relied on tenable calibrations indicates a generally conserved rate of approximately 1.5% per million years. Applying this rate to a rate-smoothed Bayesian chronogram of the ants, the Decacrema from Macaranga are inferred to be at least 12 million years old (mid-Miocene). However, using the extremes of rate variation in COI produces an age as recent as 6 million years. Our inferred timeline based on 1.5% per million years concurs with independent biogeographical events in the region reconstructed from palynological data, thus suggesting that the evolutionary histories of Decacrema and their Pachystemon hosts have been contemporaneous since the mid-Miocene. The evolution of myrmecophytism enabled Macaranga to radiate into enemy-free space, while the ants' diversification has been shaped by stem traits, host specialization, and geographic factors. We discuss the possibility that the ancient and exclusive association between Decacrema and Macaranga was facilitated by an impoverished diversity of myrmecophytes and phytoecious (obligately plant inhabiting) ants in the region.
Benthic species, though ecologically important, are vulnerable to genetic loss and population size reduction due to impacts from fishing trawls. An assessment of genetic diversity and population structure is therefore needed to assist in a resource management program. To address this issue, the two-spined yellowtail stargazer (Uranoscopus cognatus) was collected within selected locations in the Indo-West Pacific (IWP). The partial mitochondrial DNA cytochrome c oxidase subunit 1 and the nuclear DNA recombination activating gene 1 were sequenced. Genetic diversity analyses revealed that the populations were moderately to highly diversified (haplotype diversity, H = 0.490-0.900, nucleotide diversity, π = 0.0010-0.0034) except sampling station (ST) 1 and 14. The low diversity level, however was apparent only in the matrilineal marker (H = 0.118-0.216; π = 0.0004-0.0008), possibly due to stochastic factors or anthropogenic stressors. Population structure analyses revealed a retention of ancestral polymorphism that was likely due to incomplete lineage sorting in U. cognatus, and prolonged vicariance by the Indo-Pacific Barrier has partitioned them into separate stock units. Population segregation was also shown by the phenotypic divergence in allopatric populations, regarding the premaxillary protrusion, which is possibly associated with the mechanism for upper jaw movement in biomechanical feeding approaches. The moderate genetic diversity estimated for each region, in addition to past population expansion events, indicated that U. cognatus within the IWP was still healthy and abundant (except in ST1 and 14), and two stock units were identified to be subjected to a specific resource management program.
The family Poritidae formerly included 6 genera: Alveopora, Goniopora, Machadoporites, Porites, Poritipora, and Stylaraea. Morphologically, the genera can be differentiated based on the number of tentacles, the number of septa and their arrangement, the length of the polyp column, and the diameter of the corallites. However, the phylogenetic relationships within and between the genera are unknown or contentious. On the one hand, Alveopora has been transferred to the Acroporidae recently because it was shown to be more closely related to this family than to the Poritidae by previous molecular studies. On the other hand, Goniopora is morphologically similar to 2 recently described genera, Machadoporites and Poritipora, particularly with regard to the number of septa (approximately 24), but they have not yet been investigated at the molecular level. In this study, we analyzed 93 samples from all 5 poritid genera and Alveopora using 2 genetic markers (the barcoding region of the mitochondrial COI and the ITS region of the nuclear rDNA) to investigate their phylogenetic relationships and to revise their taxonomy. The reconstructed molecular trees confirmed that Alveopora is genetically distant from all poritid genera but closely related to the family Acroporidae, whereas the other genera are genetically closely related. The molecular trees also revealed that Machadoporites and Poritipora were indistinguishable from Goniopora. However, Goniopora stutchburyi was genetically isolated from the other congeneric species and formed a sister group to Goniopora together with Porites and Stylaraea, thus suggesting that 24 septa could be an ancestral feature in the Poritidae. Based on these data, we move G. stutchburyi into a new genus, Bernardpora gen. nov., whereas Machadoporites and Poritipora are merged with Goniopora.
Transferrin is a protein super-family involved in iron transport, a central process in cellular homeostasis. Throughout the evolution of vertebrates, transferrin members have diversified into distinct subfamilies including serotransferrin, ovotransferrin, lactoferrin, melanotransferrin, the inhibitor of carbonic anhydrase, pacifastin, and the major yolk protein in sea urchin. Previous phylogenetic analyses have established the branching order of the diverse transferrin subfamilies but were mostly focused on the transferrin repertoire present in mammals. Here, we conduct a comprehensive phylogenetic analysis of transferrin protein sequences in sequenced vertebrates, placing a special focus on the less-studied nonmammalian vertebrates. Our analyses uncover a novel transferrin clade present across fish, sauropsid, and amphibian genomes but strikingly absent from mammals. Our reconstructed scenario implies that this novel class emerged through a duplication event at the vertebrate ancestor, and that it was subsequently lost in the lineage leading to mammals. We detect footprints of accelerated evolution following the duplication event, which suggest positive selection and early functional divergence of this novel clade. Interestingly, the loss of this novel class of transferrin in mammals coincided with the divergence by duplication of lactoferrin and serotransferrin in this lineage. Altogether, our results provide novel insights on the evolution of iron-binding proteins in the various vertebrate groups.