Weedy rice is a close relative of domesticated rice (Oryza sativa) that competes aggressively with the crop and limits rice productivity worldwide. Most genetic studies of weedy rice have focused on populations in regions where no reproductively compatible wild Oryza species occur (North America, Europe and northern Asia). Here, we examined the population genetics of weedy rice in Malaysia, where wild rice (O. rufipogon) can be found growing in close proximity to cultivated and weedy rice. Using 375 accessions and a combined analysis of 24 neutral SSR loci and two rice domestication genes (sh4, controlling seed shattering, and Bh4, controlling hull colour), we addressed the following questions: (i) What is the relationship of Malaysian weedy rice to domesticated and wild rice, and to weedy rice strains in the USA? (ii) To what extent does the presence of O. rufipogon influence the genetic and phenotypic diversity of Malaysian weeds? (iii) What do the distributions of sh4 and Bh4 alleles and associated phenotypes reveal about the origin and contemporary evolution of Malaysian weedy rice? Our results reveal the following: independent evolutionary origins for Malaysian weeds and US strains, despite their very close phenotypic resemblance; wild-to-weed gene flow in Malaysian weed populations, including apparent adaptive introgression of seed-shattering alleles; and a prominent role for modern Malaysian cultivars in the origin and recent proliferation of Malaysian weeds. These findings suggest that the genetic complexity and adaptability of weedy crop relatives can be profoundly influenced by proximity to reproductively compatible wild and domesticated populations.
Sundaland has a dynamic geographic history because its landmasses were periodically interconnected when sea levels fell during glacial periods. Superimposed on this geographic dynamism were environmental changes related to climatic oscillations. To investigate how tropical taxa responded to such changes, we studied the divergence and demographic history of two co-distributed rainforest passerine species, Arachnothera longirostra and Malacocincla malaccensis. We sampled birds primarily from Borneo and the Malay Peninsula, which straddle the now-submerged Sunda shelf, and analysed multilocus DNA data with a variety of coalescent and gene genealogy methods. Cross-shelf divergence in both species occurred well before the last glacial maximum, i.e., before the most recent land connection. However, post-divergence gene flow occurred, and it was more pronounced in A. longirostra (a highly vagile nectarivore/insectivore) than in M. malaccensis (an understory insectivore). Despite current habitat continuity on Borneo, the population of M. malaccensis in northeastern Borneo is substantially divergent from that on the rest of the island. The NE population experienced dramatic demographic fluctuations, probably because of competition with the other population, which expanded from western Borneo after the mid-Pleistocene. In contrast, the Borneo population of A. longirostra has little structure and appears to have experienced demographic expansion 16 kya, long after it had diverged from the Malay Peninsula population (630-690 kya). Malay Peninsula populations of both species have remained relatively stable. Overall, the most recent glacial period was not the chief determinant of the evolutionary dynamics of our study species, and in this respect, they are different from temperate species.
We examined mitochondrial DNA control region sequences of 73 Kinabatangan orangutans to test the hypothesis that the phylogeographical structure of the Bornean orangutan is influenced by riverine barriers. The Lower Kinabatangan Wildlife Sanctuary contains one of the most northern populations of orangutans (Pongo pygmaeus) on Borneo and is bisected by the Kinabatangan River, the longest river in Sabah. Orang-utan samples on either side of the river were strongly differentiated with a high Phi(ST) value of 0.404 (P < 0.001). Results also suggest an east-west gradient of genetic diversity and evidence for population expansion along the river, possibly reflecting a postglacial colonization of the Kinabatangan floodplain. We compared our data with previously published sequences of Bornean orangutans in the context of river catchment structure on the island and evaluated the general relevance of rivers as barriers to gene flow in this long-lived, solitary arboreal ape.
We investigate the geographical and historical context of diversification in a complex of mutualistic Crematogaster ants living in Macaranga trees in the equatorial rain forests of Southeast Asia. Using mitochondrial DNA from 433 ant colonies collected from 32 locations spanning Borneo, Malaya and Sumatra, we infer branching relationships, patterns of genetic diversity and population history. We reconstruct a time frame for the ants' diversification and demographic expansions, and identify areas that might have been refugia or centres of diversification. Seventeen operational lineages are identified, most of which can be distinguished by host preference and geographical range. The ants first diversified 16-20 Ma, not long after the onset of the everwet forests in Sundaland, and achieved most of their taxonomic diversity during the Pliocene. Pleistocene demographic expansions are inferred for several of the younger lineages. Phylogenetic relationships suggest a Bornean cradle and major axis of diversification. Taxonomic diversity tends to be associated with mountain ranges; in Borneo, it is greatest in the Crocker Range of Sabah and concentrated also in other parts of the northern northwest coast. Within-lineage genetic diversity in Malaya and Sumatra tends to also coincide with mountain ranges. A series of disjunct and restricted distributions spanning northern northwest Borneo and the major mountain ranges of Malaya and Sumatra, seen in three pairs of sister lineages, further suggests that these regions were rain-forest refuges during drier climatic phases of the Pleistocene. Results are discussed in the context of the history of Sundaland's rain forests.
Invasive species are one of the main sources of the ongoing global loss of biodiversity. Invasive ants are known as particularly damaging invaders and their introductions are often accompanied by population-level behavioural and genetic changes that may contribute to their success. Anoplolepis gracilipes is an invasive ant that has just recently received increased attention due to its negative impact on native ecosystems. We examined the behaviour and population structure of A. gracilipes in Sabah, Malaysia. A total of 475 individuals from 24 colonies were genotyped with eight microsatellite markers. Intracolonial relatedness was high, ranging from 0.37 to 1 (mean +/- SD: 0.82 +/- 0.04), while intercolonial relatedness was low (0.0 +/- 0.02, range -0.5-0.76). We compared five distinct sampling regions in Sabah and Brunei. A three-level hierarchical F-analysis revealed high genetic differentiation among colonies within the same region, but low genetic differentiation within colonies or across regions. Overall levels of heterozygosity were unusually high (mean H(O) = 0.95, mean H(E) = 0.71) with two loci being entirely heterozygous, indicating an unusual reproductive system in this species. Bioassays revealed a negative correlation between relatedness and aggression, suggesting kinship as one factor facilitating supercolony formation in this species. Furthermore, we genotyped one individual per nest from Sabah (22 nests), Sarawak (one nest), Brunei (three nests) and the Philippines (two nests) using two mitochondrial DNA markers. We found six haplotypes, two of which included 82.1% of all sequences. Our study shows that the sampled area in Sabah consists of a mosaic of differently interrelated nests in different stages of colony establishment. While some of the sampled colonies may belong to large supercolonies, others are more likely to represent recently introduced or dispersed propagules that are just beginning to expand.
The diamondback moth, Plutella xylostella, is renowned for developing resistance to insecticides and causing significant economic damage to Brassica vegetable crops throughout the world. Yet despite its economic importance, little is known about the population structure and movement patterns of this pest both at local and regional scales. In Australia, the movement patterns and insecticide resistance status of P. xylostella infesting canola, vegetables, forage brassicas and weeds have fundamental implications for the management of this pest. Here we use six polymorphic microsatellite loci to investigate population structure and gene flow in Australian populations of P. xylostella. Samples of P. xylostella from New Zealand, Malaysia, Indonesia and Kenya were also scored at these loci. We found no evidence of population structure within Australia, with most populations having low inbreeding coefficients and in Hardy-Weinberg equilibrium. In addition, a sample from the North Island of New Zealand was indistinguishable from the Australian samples. However, large genetic differences were found between the Australia/New Zealand samples and samples from Kenya, Malaysia and Indonesia. There was no relationship between genetic distance and geographic distance among Australian and New Zealand samples. Two of the loci were found to have null alleles, the frequency of which was increased in the populations outside the Australia/New Zealand region. We discuss these results with reference to insecticide resistance management strategies for P. xylostella in Australia.
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.
Macaranga (Euphorbiaceae) includes about 280 species with a palaeotropic distribution. The genus not only comprises some of the most prominent pioneer tree species in Southeast Asian lowland dipterocarp forests, it also exhibits a substantial radiation of ant-plants (myrmecophytes). Obligate ant-plant mutualisms are formed by about 30 Macaranga species and 13 ant species of the genera Crematogaster or Camponotus. To improve our understanding of the co-evolution of the ants and their host plants, we aim at reconstructing comparative organellar phylogeographies of both partners across their distributional range. Preliminary evidence indicated that chloroplast DNA introgression among closely related Macaranga species might occur. We therefore constructed a comprehensive chloroplast genealogy based on DNA sequence data from the noncoding ccmp2, ccmp6, and atpB-rbcL regions for 144 individuals from 41 Macaranga species, covering all major evolutionary lineages within the three sections that contain myrmecophytes. A total of 88 chloroplast haplotypes were identified, and grouped into a statistical parsimony network that clearly distinguished sections and well-defined subsectional groups. Within these groups, the arrangement of haplotypes followed geographical rather than taxonomical criteria. Thus, up to six chloroplast haplotypes were found within single species, and up to seven species shared a single haplotype. The spatial distribution of the chloroplast types revealed several dispersals between the Malay Peninsula and Borneo, and a deep split between Sabah and the remainder of Borneo. Our large-scale chloroplast genealogy highlights the complex history of migration, hybridization, and speciation in the myrmecophytes of the genus Macaranga. It will serve as a guideline for adequate sampling and data interpretation in phylogeographic studies of individual Macaranga species and species groups.
Behavioural observations suggest that orang-utans are semi-solitary animals with females being philopatric and males roaming more widely in search of receptive partners, leading to the prediction that females are more closely related than males at any given site. In contrast, our study presents evidence for male and female philopatry in the orang-utan. We examined patterns of relatedness and parentage in a wild orang-utan population in Borneo using noninvasively collected DNA samples from animals observed to defecate, and microsatellite markers to assess dispersal and mating strategies. Surprisingly, resident females were equally as related to other resident females (mean r(xy) = 0.303) as resident males were to other resident males (mean r(xy) = 0.305). Moreover, resident females were more related to each other and to the resident males than they were to nonresident females, and resident males were more related to each other (and resident females) than they were to nonresident males. We assigned genetic mothers to 12 individuals in the population, while sires could be identified for eight. Both flanged males and unflanged males achieved paternity, similar to findings reported for Sumatran orang-utans.
We investigated the genetic structure within and among Bornean orang-utans (Pongo pygmaeus) in forest fragments of the Lower Kinabatangan flood plain in Sabah, Malaysia. DNA was extracted from hair and faecal samples for 200 wild individuals collected during boat surveys on the Kinabatangan River. Fourteen microsatellite loci were used to characterize patterns of genetic diversity. We found that genetic diversity was high in the set of samples (mean H(E) = 0.74) and that genetic differentiation was significant between the samples (average F(ST) = 0.04, P < 0.001) with F(ST) values ranging from low (0.01) to moderately large (0.12) values. Pairwise F(ST) values were significantly higher across the Kinabatangan River than between samples from the same river side, thereby confirming the role of the river as a natural barrier to gene flow. The correlation between genetic and geographical distance was tested by means of a series of Mantel tests based on different measures of geographical distance. We used a Bayesian method to estimate immigration rates. The results indicate that migration is unlikely across the river but cannot be completely ruled out because of the limited F(ST) values. Assignment tests confirm the overall picture that gene flow is limited across the river. We found that migration between samples from the same side of the river had a high probability indicating that orang-utans used to move relatively freely between neighbouring areas. This strongly suggests that there is a need to maintain migration between isolated forest fragments. This could be done by restoring forest corridors alongside the river banks and between patches.
We examined differences in pollen dispersal efficiency between 2 years in terms of both spatial dispersal range and genetic relatedness of pollen in a tropical emergent tree, Dipterocarpus tempehes. The species was pollinated by the giant honeybee (Apis dorsata) in a year of intensive community-level mass-flowering or general flowering (1996), but by several species of moths in a year of less-intensive general flowering (1998). We carried out paternity analysis based on six DNA microsatellite markers on a total of 277 mature trees forming four spatially distinct subpopulations in a 70 ha area, and 147 and 188 2-year-old seedlings originating from seeds produced in 1996 and 1998 (cohorts 96 and 98, respectively). Outcrossing rates (0.93 and 0.96 for cohorts 96 and 98, respectively) did not differ between years. Mean dispersal distances (222 and 192 m) were not significantly different between the 2 years but marginally more biased to long distance in 1996. The mean relatedness among cross-pollinated seedlings sharing the same mothers in cohort 96 was lower than that in cohort 98. This can be attributed to the two facts that the proportion of intersubpopulations pollen flow among cross-pollination events was marginally higher in cohort 96 (44%) than in cohort 98 (33%), and that mature trees within the same subpopulations are genetically more related to each other than those between different subpopulations. We conclude that D. tempehes maintained effective pollen dispersal in terms of outcrossing rate and pollen dispersal distance in spite of the large difference in foraging characteristics between two types of pollinators. In terms of pollen relatedness, however, a slight difference was suggested between years in the level of biparental inbreeding.
Tropical rainforests in South-East Asia have been affected by climatic fluctuations during past glacial eras. To examine how the accompanying changes in land areas and temperature have affected the genetic properties of rainforest trees in the region, we investigated the phylogeographic patterns of a widespread dipterocarp species, Shorea leprosula. Two types of DNA markers were used: expressed sequence tag-based simple sequence repeats and chloroplast DNA (cpDNA) sequence variations. Both sets of markers revealed clear genetic differentiation between populations in Borneo and those in the Malay Peninsula and Sumatra (Malay/Sumatra). However, in the south-western part of Borneo, genetic admixture of the lineages was observed in the two marker types. Coalescent simulation based on cpDNA sequence variation suggested that the two lineages arose 0.28-0.09 million years before present and that following their divergence migration from Malay/Sumatra to Borneo strongly exceeded migration in the opposite direction. We conclude that the genetic structure of S. leprosula was largely formed during the middle Pleistocene and was subsequently modified by eastward migration across the subaerially exposed Sunda Shelf.
Selective logging and forest conversion to oil palm agriculture are rapidly altering tropical forests. However, functional responses of the soil microbiome to these land-use changes are poorly understood. Using 16S rRNA gene and shotgun metagenomic sequencing, we compared composition and functional attributes of soil biota between unlogged, once-logged and twice-logged rainforest, and areas converted to oil palm plantations in Sabah, Borneo. Although there was no significant effect of logging history, we found a significant difference between the taxonomic and functional composition of both primary and logged forests and oil palm. Oil palm had greater abundances of genes associated with DNA, RNA, protein metabolism and other core metabolic functions, but conversely, lower abundance of genes associated with secondary metabolism and cell-cell interactions, indicating less importance of antagonism or mutualism in the more oligotrophic oil palm environment. Overall, these results show a striking difference in taxonomic composition and functional gene diversity of soil microorganisms between oil palm and forest, but no significant difference between primary forest and forest areas with differing logging history. This reinforces the view that logged forest retains most features and functions of the original soil community. However, networks based on strong correlations between taxonomy and functions showed that network complexity is unexpectedly increased due to both logging and oil palm agriculture, which suggests a pervasive effect of both land-use changes on the interaction of soil microbes.
The Indo-Malayan bioregion has provided some of the most spectacular discoveries of new vertebrate species (e.g. saola, khanyou, bare-faced bulbul) over the last 25 years. Yet, very little is known about the processes that led to the current biodiversity in this region. We reconstructed the phylogeographic history of a group of closely related passerines, the Alophoixus bulbuls. These birds are continuously distributed in Indo-Malaya around the Thailand lowlands such that their distribution resembles a ring. Our analyses revealed a single colonization event of the mainland from Sundaland with sequential divergence of taxa from southwest to northeast characterized by significant gene flow between parapatric taxa, and reduced or ancient gene flow involving the two taxa at the extremities of the ring. We detected evidence of population expansion in two subspecies, including one that was involved in the closing of the ring. Hence, our analyses indicate that the diversification pattern of Alophoixus bulbuls fits a ring species model driven by geographic isolation. To our knowledge, the Alophoixus bulbuls represent the first case of a putative broken ring species complex in Indo-Malaya. We also discuss the implications of our results on our understanding of the biogeography in Indo-Malaya.
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.
Vivipary with precocious seedlings in mangrove plants was thought to be a hindrance to long-range dispersal. To examine the extent of seedling dispersal across oceans, we investigated the phylogeny and genetic structure among East Asiatic populations of Kandelia candel based on organelle DNAs. In total, three, 28 and seven haplotypes of the chloroplast DNA (cpDNA) atpB-rbcL spacer, cpDNA trnL-trnF spacer, and mitochondrial DNA (mtDNA) internal transcribed spacer (ITS) were identified, respectively, from 202 individuals. Three data sets suggested consistent phylogenies recovering two differentiated lineages corresponding to geographical regions, i.e. northern South-China-Sea + East-China-Sea region and southern South-China-Sea region (Sarawak). Phylogenetically, the Sarawak population was closely related to the Ranong population of western Peninsula Malaysia instead of other South-China-Sea populations, indicating its possible origin from the Indian Ocean Rim. No geographical subdivision was detected within the northern geographical region. An analysis of molecular variance (AMOVA) revealed low levels of genetic differentiation between and within mainland and island populations (phiCT = 0.015, phiSC = 0.037), indicating conspicuous long-distance seedling dispersal across oceans. Significant linkage disequilibrium excluded the possibility of recurrent homoplasious mutations as the major force causing phylogenetic discrepancy between mtDNA and the trnL-trnF spacer within the northern region. Instead, relative ages of alleles contributed to non-random chlorotype-mitotype associations and tree inconsistency. Widespread distribution and random associations (chi2 = 0.822, P = 0.189) of eight hypothetical ancestral cytotypes indicated the panmixis of populations of the northern geographical region as a whole. In contrast, rare and recently evolved alleles were restricted to marginal populations, revealing some preferential directional migration.
Pollen flow and population genetic structure among 30 potentially flowering individuals of Neobalanocarpus heimii, a tropical emergent tree, were investigated in a lowland tropical rainforest of Malaysia using microsatellite polymorphism. The 248 offspring in the vicinity of five reproductive trees of the 30 potentially flowering trees were used in paternity analysis for pollen-flow study. Four primer pairs, developed in different species of dipterocarps, were adopted to detect microsatellite polymorphism. Based upon microsatellite polymorphism, pollen flow and seed migration were detected. Pollen-flow events of more than 400 m were observed directly, based on paternity analysis in the study plot. The estimated average mating distance of the five reproductive trees was 524 m. This result suggests that reproduction of this species is mediated by a long-distance pollinator. The haplotypes of some offspring were not compatible with the nearest reproductive tree. Thus, the results suggest that some seeds are dispersed by a seed dispersal vector. Investigation of genetic structure showed significant and negative correlation of genetic relatedness and spatial distances between the 30 potentially flowering trees, but this correlation was weak. We suggest that long-distance gene flow and seed migration are responsible for the poorly developed genetic structure of this species.
Analyses of the spatial distribution pattern, spatial genetic structure and of genetic diversity were carried out in two tropical tree species with contrasting breeding systems and different ploidy levels using a 50-ha demographic plot in a lowland dipterocarp forest in Peninsular Malaysia. Shorea leprosula is a diploid and predominantly outcrossed species, whereas S. ovalis ssp. sericea is an autotetraploid species with apomictic mode of reproduction. Genetic diversity parameters estimated for S. leprosula using microsatellite were consistently higher than using allozyme. In comparisons with S. leprosula and other tropical tree species, S. ovalis ssp. sericea also displayed relatively high levels of genetic diversity. This might be explained by the lower pressure of genetic drift due to tetrasomic inheritance, and for autotetraploids each locus can accommodate up to four different alleles and this allows maintenance of more alleles at individual loci. The observed high levels of genetic diversity in S. ovalis ssp. sericea can also be due to a random retention of more heterogeneous individuals in the past, and the apomictic mode of reproduction might be an evolutionary strategy, which allows the species to maintain high levels of genetic diversity. The spatial distribution pattern analyses of both species showed significant levels of aggregation at small and medium but random distribution at the big diameter-class. The decrease in magnitude of spatial aggregation from small- to large-diameter classes might be due to compensatory mortality during recruitment and survival under competitive thinning process. Spatial genetic structure analyses for both species revealed significant spatial genetic structure for short distances in all the three diameter-classes. The magnitude of spatial genetic structure in both species was observed to be decreasing from smaller- to larger-diameter classes. The high spatial genetic structuring observed in S. ovalis ssp. sericea at the small-diameter class is due primarily to limited seed dispersal and apomictic mode of reproduction. The similar observation in S. leprosula, however, can be explained by limited seed and pollen dispersal, which supports further the fact that the species is pollinated by weak fliers, mainly of Thrips and Megalurothrips in the lowland dipterocarp forest.
Plasmodium knowlesi is a significant cause of human malaria transmitted as a zoonosis from macaque reservoir hosts in South-East Asia. Microsatellite genotyping has indicated that human infections in Malaysian Borneo are an admixture of two highly divergent sympatric parasite subpopulations that are, respectively, associated with long-tailed macaques (Cluster 1) and pig-tailed macaques (Cluster 2). Whole-genome sequences of clinical isolates subsequently confirmed the separate clusters, although fewer of the less common Cluster 2 type were sequenced. Here, to analyse population structure and genomic divergence in subpopulation samples of comparable depth, genome sequences were generated from 21 new clinical infections identified as Cluster 2 by microsatellite analysis, yielding a cumulative sample size for this subpopulation similar to that for Cluster 1. Profound heterogeneity in the level of intercluster divergence was distributed across the genome, with long contiguous chromosomal blocks having high or low divergence. Different mitochondrial genome clades were associated with the two major subpopulations, but limited exchange of haplotypes from one to the other was evident, as was also the case for the maternally inherited apicoplast genome. These findings indicate deep divergence of the two sympatric P. knowlesi subpopulations, with introgression likely to have occurred recently. There is no evidence yet of specific adaptation at any introgressed locus, but the recombinant mosaic types offer enhanced diversity on which selection may operate in a currently changing landscape and human environment. Loci responsible for maintaining genetic isolation of the sympatric subpopulations need to be identified in the chromosomal regions showing fixed differences.