Using multi-way admixture mapping to elucidate TB susceptibility in the South African Coloured population

Background

The admixed South African Coloured population is ideally suited to the discovery of tuberculosis susceptibility genetic variants and their probable ethnic origins, but previous attempts at finding such variants using genome-wide admixture mapping were hampered by the inaccuracy of local ancestry inference. In this study, we infer local ancestry using the novel algorithm implemented in RFMix, with the emphasis on identifying regions of excess San or Bantu ancestry, which we hypothesize may harbour TB susceptibility genes.

Results

Using simulated data, we demonstrate reasonable accuracy of local ancestry inference by RFMix, with a tendency towards miss-calling San ancestry as Bantu. Regions with either excess San ancestry or excess African (San or Bantu) ancestry are less likely to be affected by this bias, and we therefore proceeded to identify such regions, found in cases but not in controls (642 cases and 91 controls). A number of promising regions were found (overall p-values of 7.19×10^-5 for San ancestry and <2.00×10^-16 for African ancestry), including chromosomes 15q15 and 17q22, which are close to genomic regions previously implicated in TB. Promising immune-related susceptibility genes such as the GADD45A, OSM and B7-H5 genes are also harboured in the identified regions.

Conclusion

Admixture mapping is feasible in the South African Coloured population and a number of novel TB susceptibility genomic regions were uncovered.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1021) contains supplementary material, which is available to authorized users.     

Morphological variation of genetically confirmed Alouatta Pigra × A. palliata hybrids from a natural hybrid zone in Tabasco,Mexico

While hybridization has been reported for a large number of primate taxa, there is a general lack of data on hybrid morphology for wild individuals with known genetic ancestry. A confirmed hybrid zone for the closely related Neotropical primates Alouatta palliata and A. pigra has provided a unique opportunity to study primate hybrid morphological variation. Here we used molecular evidence based on mitochondrial, Y‐chromosome, and autosomal data to assess hybrid ancestry. We conducted univariate and multivariate statistical comparisons of morphometric data collected from individuals both outside and within the hybrid zone in Tabasco, Mexico. Our results show that of all the hybrids detected (N = 128), only 12% of them were approximately genetically intermediate, and none of them were first generation hybrids. Univariate pairwise comparisons among parental individuals, multigenerational backcrossed hybrids, and intermediate hybrids showed that overall, multigenerational backcrossed hybrids resemble the parental species with which they share most of their alleles. Conversely, intermediates were highly variable. Similarly, principal component analysis depicts an overlap between the parental species and their backcrosses when considering overall morphological differences. Finally, discriminant function analysis of the morphological variables was overall unreliable for classifying individuals into their assigned genotypic classes. Taken together, our results suggest that primate natural hybridization studies should incorporate molecular methods for determining ancestry, because morphology may not always be a reliable indicator of hybrid status. Hybrid zones could comprise a large number of multigenerational backcrossed hybrids that are indistinguishable from the parental species. The implications for studying hybridization in the primate fossil record are discussed. Am J Phys Anthropol, 2013. © 2012 Wiley Periodicals, Inc.     

In situ hybridization confirms jumping nucleolus organizing regions in Allium

In situ hybridization with a ¹²⁵I-rDNA clone from Vicia faba was performed against Allium cepa and three strains of top onion, which represent hybrids between A. cepa and A. fistulosum. In principle, the labelling patterns correspond to the patterns of the silver-stained nucleolus organizing regions (NORs) in the same species. This strongly supports the inference drawn from the Ag-NOR patterns that NORs can jump between terminal heterochromatin blocks of different Allium chromosomes in the parental species A. cepa as well as in their interspecific hybrids.     

Heterogeneity and concordance in locus‐specific differentiation and introgression between species of towhees

The maintenance or breakdown of reproductive isolation is an observable outcome of secondary contact between species. In cases where hybrids beyond the F1 are formed, the representation of each species' ancestry can vary dramatically among genomic regions. This genomic heterogeneity in ancestry and introgression can offer insight into evolutionary processes, particularly if introgression is compared in multiple hybrid zones. Similarly, considerable heterogeneity exists across the genome in the extent to which populations and species have diverged, reflecting the combined effects of different evolutionary processes on genetic variation. We studied hybridization across two hybrid zones of two phenotypically well‐differentiated bird species in Mexico (Pipilo maculatus and P. ocai), to investigate genomic heterogeneity in differentiation and introgression. Using genotyping‐by‐sequencing (GBS) and hierarchical Bayesian models, we genotyped 460 birds at over 41 000 single nucleotide polymorphism (SNP) loci. We identified loci exhibiting extreme introgression relative to the genome‐wide expectation using a Bayesian genomic cline model. We also estimated locus‐specific F_ST and identified loci with exceptionally high genetic divergence between the parental species. We found some concordance of locus‐specific introgression in the two independent hybrid zones (6–20% of extreme loci shared across zones), reflecting areas of the genome that experience similar gene flow when the species interact. Additionally, heterogeneity in introgression and divergence across the genome revealed another subset of loci under the influence of locally specific factors. These results are consistent with a history in which reproductive isolation has been influenced by a common set of loci in both hybrid zones, but where local environmental and stochastic factors also lead to genomic differentiation.     

A Continuous Correlated Beta Process Model for Genetic Ancestry in Admixed Populations

Admixture and recombination create populations and genomes with genetic ancestry from multiple source populations. Analyses of genetic ancestry in admixed populations are relevant for trait and disease mapping, studies of speciation, and conservation efforts. Consequently, many methods have been developed to infer genome-average ancestry and to deconvolute ancestry into continuous local ancestry blocks or tracts within individuals. Current methods for local ancestry inference perform well when admixture occurred recently or hybridization is ongoing, or when admixture occurred in the distant past such that local ancestry blocks have fixed in the admixed population. However, methods to infer local ancestry frequencies in isolated admixed populations still segregating for ancestry do not exist. In the current paper, I develop and test a continuous correlated beta process model to fill this analytical gap. The method explicitly models autocorrelations in ancestry frequencies at the population-level and uses discriminant analysis of SNP windows to take advantage of ancestry blocks within individuals. Analyses of simulated data sets show that the method is generally accurate such that ancestry frequency estimates exhibited low root-mean-square error and were highly correlated with the true values, particularly when large (±10 or ±20) SNP windows were used. Along these lines, the proposed method outperformed post hoc inference of ancestry frequencies from a traditional hidden Markov model (i.e., the linkage model in structure), particularly when admixture occurred more distantly in the past with little on-going gene flow or was followed by natural selection. The reliability and utility of the method was further assessed by analyzing genetic ancestry in an admixed human population (Uyghur) and three populations from a hybrid zone between Mus domesticus and M. musculus. Considerable variation in ancestry frequencies was detected within and among chromosomes in the Uyghur, with a large region of excess French ancestry harboring a gene with a known disease association. Similar variation was detected in the mouse hybrid zone, with notable constancy in regions of excess ancestry among admixed populations. By filling what has been an analytical gap, the proposed method should be a useful tool for many biologists. A computer program (popanc), written in C++, has been developed based on the proposed method and is available on-line at http://sourceforge.net/projects/popanc/.     

A consensus tree approach for reconstructing human evolutionary history and detecting population substructure

The random accumulation of variations in the human genome over time implicitly encodes a history of how human populations have arisen, dispersed, and intermixed since we emerged as a species. Reconstructing that history is a challenging computational and statistical problem but has important applications both to basic research and to the discovery of genotype-phenotype correlations. We present a novel approach to inferring human evolutionary history from genetic variation data. We use the idea of consensus trees, a technique generally used to reconcile species trees from divergent gene trees, adapting it to the problem of finding robust relationships within a set of intraspecies phylogenies derived from local regions of the genome. Validation on both simulated and real data shows the method to be effective in recapitulating known true structure of the data closely matching our best current understanding of human evolutionary history. Additional comparison with results of leading methods for the problem of population substructure assignment verifies that our method provides comparable accuracy in identifying meaningful population subgroups in addition to inferring relationships among them. The consensus tree approach thus provides a promising new model for the robust inference of substructure and ancestry from large-scale genetic variation data.     

PHYLOGENETIC INFERENCE OF NUPTIAL TRAIT EVOLUTION IN THE CONTEXT OF ASYMMETRICAL INTROGRESSION IN NORTH AMERICAN DARTERS (TELEOSTEI)

Introgressive hybridization and incomplete lineage sorting complicate the inference of phylogeny, and available species‐tree methods do not simultaneously account for these processes. Both hybridization and ancestral polymorphism have been invoked to explain divergent phylogenies inferred from different datasets for Stigmacerca, a clade of 11 North American darter species. Species of Stigmacerca are characterized by a mating system involving parental care with males guarding nesting territories and fertilized eggs. Males of four species of Stigmacerca develop egg‐mimic nuptial structures on their second dorsal fins during the breeding season. Previous phylogenies suggest contrasting scenarios for the evolution of this nuptial trait. Using a combination of coalescent‐based methods, we analyzed a dataset comprising a mitochondrial gene and 15 nuclear loci to estimate relationships and simultaneously test for introgressive hybridization. Our analyses identified several instances of interspecific gene flow involving both cytoplamsmic haplotypes and nuclear alleles. The new phylogeny was used to infer a single origin and recent loss of egg‐mimic structures in Stigmacerca and led to the discovery of a phylogenetically distinct species. Our results highlight the limited strategies available to account for introgressive hybridization in the inference of species relationships and the likely effects of this process on reconstructing trait evolution.     

A performance study of the impact of recombination on species tree analysis

Background

The most widely used state-of-the-art methods for reconstructing species phylogenies from genomic sequence data assume that sampled loci are identically and independently distributed. In principle, free recombination between loci and a lack of intra-locus recombination are necessary to satisfy this assumption. Few studies have quantified the practical impact of recombination on species tree inference methods, and even fewer have used genomic sequence data for this purpose. One prominent exception is the 2012 study of Lanier and Knowles. A main finding from the study was that species tree inference methods are relatively robust to intra-locus recombination, assuming free recombination between loci. The latter assumption means that the open question regarding the impact of recombination on species tree analysis is not fully resolved.

Results

The goal of this study is to further investigate this open question. Using simulations based upon the multi-species coalescent-with-recombination model as well as empirical datasets, we compared common pipeline-based techniques for inferring species phylogenies. The simulation conditions included a range of dataset sizes and several choices for recombination rate which was either uniform across loci or incorporated recombination hotspots. We found that pipelines which explicitly utilize inferred recombination breakpoints to delineate recombination-free intervals result in greater accuracy compared to widely used alternatives that preprocess sequences based upon linkage disequilibrium decay. Furthermore, the use of a relatively simple approach for recombination breakpoint inference does not degrade the accuracy of downstream species tree inference compared to more accurate alternatives.

Conclusions

Our findings clarify the impact of recombination upon current phylogenomic pipelines for species tree inference. Pipeline-based approaches which utilize inferred recombination breakpoints to densely sample loci across genomic sequences can tolerate intra-locus recombination and violations of the assumption of free recombination between loci.

     

Haplotype inference in crossbred populations without pedigree information

Background

Current methods for haplotype inference without pedigree information assume random mating populations. In animal and plant breeding, however, mating is often not random. A particular form of nonrandom mating occurs when parental individuals of opposite sex originate from distinct populations. In animal breeding this is called crossbreeding and hybridization in plant breeding. In these situations, association between marker and putative gene alleles might differ between the founding populations and origin of alleles should be accounted for in studies which estimate breeding values with marker data. The sequence of alleles from one parent constitutes one haplotype of an individual. Haplotypes thus reveal allele origin in data of crossbred individuals.

Results

We introduce a new method for haplotype inference without pedigree that allows nonrandom mating and that can use genotype data of the parental populations and of a crossbred population. The aim of the method is to estimate line origin of alleles. The method has a Bayesian set up with a Dirichlet Process as prior for the haplotypes in the two parental populations. The basic idea is that only a subset of the complete set of possible haplotypes is present in the population.

Conclusion

Line origin of approximately 95% of the alleles at heterozygous sites was assessed correctly in both simulated and real data. Comparing accuracy of haplotype frequencies inferred with the new algorithm to the accuracy of haplotype frequencies inferred with PHASE, an existing algorithm for haplotype inference, showed that the DP algorithm outperformed PHASE in situations of crossbreeding and that PHASE performed better in situations of random mating.     

Identification of recent hybridization between gray wolves and domesticated dogs by SNP genotyping

The ability to detect recent hybridization between dogs and wolves is important for conservation and legal actions, which often require accurate and rapid resolution of ancestry. The availability of a genetic test for dog–wolf hybrids would greatly support federal and legal enforcement efforts, particularly when the individual in question lacks prior ancestry information. We have developed a panel of 100 unlinked ancestry-informative SNP markers that can detect mixed ancestry within up to four generations of dog–wolf hybridization based on simulations of seven genealogical classes constructed following the rules of Mendelian inheritance. We establish 95 % confidence regions around the spatial clustering of each genealogical class using a tertiary plot of allele dosage and heterozygosity. The first- and second-backcrossed-generation hybrids were the most distinct from parental populations, with >90 % correctly assigned to genealogical class. In this article we provide a tool kit with population-level statistical quantification that can detect recent dog–wolf hybridization using a panel of dog–wolf ancestry-informative SNPs with divergent allele frequency distributions.     

TESS3: fast inference of spatial population structure and genome scans for selection

Geography and landscape are important determinants of genetic variation in natural populations, and several ancestry estimation methods have been proposed to investigate population structure using genetic and geographic data simultaneously. Those approaches are often based on computer‐intensive stochastic simulations and do not scale with the dimensions of the data sets generated by high‐throughput sequencing technologies. There is a growing demand for faster algorithms able to analyse genomewide patterns of population genetic variation in their geographic context. In this study, we present TESS3 , a major update of the spatial ancestry estimation program TESS . By combining matrix factorization and spatial statistical methods, TESS3 provides estimates of ancestry coefficients with accuracy comparable to TESS and with run‐times much faster than the Bayesian version. In addition, the TESS3 program can be used to perform genome scans for selection, and separate adaptive from nonadaptive genetic variation using ancestral allele frequency differentiation tests. The main features of TESS3 are illustrated using simulated data and analysing genomic data from European lines of the plant species Arabidopsis thaliana.     

Quantitative evaluation of hybridization and the impact on biodiversity conservation

Anthropogenic hybridization is an increasing conservation threat worldwide. In South Africa, recent hybridization is threatening numerous ungulate taxa. For example, the genetic integrity of the near‐threatened bontebok (Damaliscus pygargus pygargus) is threatened by hybridization with the more common blesbok (D. p. phillipsi). Identifying nonadmixed parental and admixed individuals is challenging based on the morphological traits alone; however, molecular analyses may allow for accurate detection. Once hybrids are identified, population simulation software may assist in determining the optimal conservation management strategy, although quantitative evaluation of hybrid management is rarely performed. In this study, our objectives were to describe species‐wide and localized rates of hybridization in nearly 3,000 individuals based on 12 microsatellite loci, quantify the accuracy of hybrid assignment software (STRUCTURE and NEWHYBRIDS), and determine an optimal threshold of bontebok ancestry for management purposes. According to multiple methods, we identified 2,051 bontebok, 657 hybrids, and 29 blesbok. More than two‐thirds of locations contained at least some hybrid individuals, with populations varying in the degree of introgression. HYBRIDLAB was used to simulate four generations of coexistence between bontebok and blesbok, and to optimize a threshold of ancestry, where most hybrids will be detected and removed, and the fewest nonadmixed bontebok individuals misclassified as hybrids. Overall, a threshold Q‐value (admixture coefficient) of 0.90 would remove 94% of hybrid animals, while a threshold of 0.95 would remove 98% of hybrid animals but also 8% of nonadmixed bontebok. To this end, a threshold of 0.90 was identified as optimal and has since been implemented in formal policy by a provincial nature conservation agency. Due to widespread hybridization, effective conservation plans should be established and enforced to conserve native populations that are genetically unique.     

The allopolyploid origin of Narcissus obsoletus (Alliaceae): identification of parental genomes by karyotype characterization and genomic in situ hybridization

Karyological and genomic in situ hybridization (GISH) approaches provided evidence of the parentage and origin of the hybrid species Narcissus obsoletus. Here, we demonstrate that the putative parental species, N. serotinus L. and N. elegans (Haw.) Spach, recently proposed because of their intermediate morphological traits, have participated in the hybridization process forming this taxon. Karyotype characterization of parental genomes in populations from S Spain and N Morocco has revealed differences in chromosome length and karyotype asymmetry, highlighting their diploid nature. Multicolour GISH on metaphase plates of N. obsoletus, with N. serotinus and N. elegans DNA used as probes, showed differential fluorescent staining of 10 and 20 chromosomes from parental genomes, respectively. Both parental genomes were detected in the allopolyploid, albeit in a duplicated manner. Secondary hybridization between N. obsoletus and N. serotinus was also detected karyologically. Little karyological differentiation between different geographic regions was found in either N. serotinus or N. obsoletus. © 2009 The Linnean Society of London, Botanical Journal of the Linnean Society, 2009, 159 , 477–498.     

Are linear‐leaved Potamogeton hybrids really so rare? Molecular evidence for multiple hybridizations between P. acutifolius and P. compressus in central Europe

Due to insufficient morphological differentiation, hybridization among linear‐leaved species of Potamogeton is rarely observed and recognized. Here, we applied molecular tools (sequencing of nrITS and rpl32‐trnL cpDNA intergenic spacer) to study several natural, central European populations of morphologically intermediate forms between two closely related species, P. acutifolius and P. compressus, to examine their taxonomical status and test whether they represented the variation range of the species or are results of interspecific hybridization. Both DNA regions provided distinguishing characters between putative parental species. The ITS sequences from all the morphologically intermediate specimens displayed an additive pattern combining the variation of both parental species, providing evidence for the hybrid origin of all these plants. This case study suggests that hybrids of linear‐leaved Potamogeton species are more common than usually believed, but have been difficult to identify without molecular tools.     

Enhanced Methods for Local Ancestry Assignment in Sequenced Admixed Individuals

Inferring the ancestry at each locus in the genome of recently admixed individuals (e.g., Latino Americans) plays a major role in medical and population genetic inferences, ranging from finding disease-risk loci, to inferring recombination rates, to mapping missing contigs in the human genome. Although many methods for local ancestry inference have been proposed, most are designed for use with genotyping arrays and fail to make use of the full spectrum of data available from sequencing. In addition, current haplotype-based approaches are very computationally demanding, requiring large computational time for moderately large sample sizes. Here we present new methods for local ancestry inference that leverage continent-specific variants (CSVs) to attain increased performance over existing approaches in sequenced admixed genomes. A key feature of our approach is that it incorporates the admixed genomes themselves jointly with public datasets, such as 1000 Genomes, to improve the accuracy of CSV calling. We use simulations to show that our approach attains accuracy similar to widely used computationally intensive haplotype-based approaches with large decreases in runtime. Most importantly, we show that our method recovers comparable local ancestries, as the 1000 Genomes consensus local ancestry calls in the real admixed individuals from the 1000 Genomes Project. We extend our approach to account for low-coverage sequencing and show that accurate local ancestry inference can be attained at low sequencing coverage. Finally, we generalize CSVs to sub-continental population-specific variants (sCSVs) and show that in some cases it is possible to determine the sub-continental ancestry for short chromosomal segments on the basis of sCSVs.     

Experimental Hybridization Reveals Biased Inheritance of the Internal Transcribed Spacer of the Nuclear Ribosomal DNA in Begonia ×taipeiensis

Begonia × taipeiensis C.-I Peng, a naturally occurring hybrid resulting from B. formosana × B. aptera, in Taiwan. To understand the inheritance of ribosomal DNA in unidirectional hybridization, experiments were conducted using B. formosana and B. aptera as ovule and pollen donors, respectively. The internal transcribed spacer region (ITS) of nuclear ribosomal DNA was amplified from the artificial hybrids, parental species, and natural hybrids. In contrast to the single type of ITS in the parental species, multiple sequences were cloned from both natural and artificial hybrids. A split decomposition network based on ITS nucleotide variation revealed that all but one (clone B14) of the hybrid sequences were “phylogenetically” closely related to B. formosana. Apparently, in such unidirectional hybridization, maternal DNA provided most of phylogenetic information. In the hybrid sequences, in addition to additive polymorphisms inherited from maternal (38.1%) and paternal (30.1%) plants, a novel nucleotide composition (31.8%) was also detected. The “new” characters are seen as noise in phylogenetic inference. They were probably obtained via intramolecular recombination, as gene conversion was not detected. The occurrence of genetic recombination appeared to be nonrandom, with a higher frequency in the ITS1 (3.14%) and ITS2 (3.42%) regions than in the 5.8S RNA gene (2.22%). Given the lack of sexual recombination in B. × taipeiensis and short time span, unequal crossing-over likely contributed to the heterogeneity of the ITS composition in the nuclear genome. Although the sterile hybrids have not attained their own lineage independent from the parental species, a high level of genetic diversity is transmitted asexually and maintained in these plants. Received 7 February 2001/ Accepted in revised form 31 May 2001     

AN EVALUATION OF THE HYBRID SPECIATION HYPOTHESIS FOR XIPHOPHORUS CLEMENCIAE BASED ON WHOLE GENOME SEQUENCES

Once thought rare in animal taxa, hybridization has been increasingly recognized as an important and common force in animal evolution. In the past decade, a number of studies have suggested that hybridization has driven speciation in some animal groups. We investigate the signature of hybridization in the genome of a putative hybrid species, Xiphophorus clemenciae, through whole genome sequencing of this species and its hypothesized progenitors. Based on analysis of this data, we find that X. clemenciae is unlikely to have been derived from admixture between its proposed parental species. However, we find significant evidence for recent gene flow between Xiphophorus species. Although we detect genetic exchange in two pairs of species analyzed, the proportion of genomic regions that can be attributed to hybrid origin is small, suggesting that strong behavioral premating isolation prevents frequent hybridization in Xiphophorus. The direction of gene flow between species is potentially consistent with a role for sexual selection in mediating hybridization.     

MORPHOLOGICAL VARIABILITY IN GENETICALLY DEFINED CATEGORIES OF ANURAN HYBRIDS

Hybridization phenomena in anurans have traditionally been studied through morphological comparisons, under the assumption that various hybrids (e.g., F₁‘s, backcrosses) are predictably intermediate to parental species. We critically evaluate this assumption by examining morphology in genetically categorized hybrids between the treefrogs Hyla cinerea and H. gratiosa. A total of 202 frogs from a hybridizing population in Alabama were assayed for allozyme and mitochondrial DNA genotype and for a large suite of osteological characters. Discriminant analyses demonstrated distinct morphological separation between the genetically “pure” parental species. Morphometric analyses of genetically identified hybrids showed: 1) an extreme range of phenotypic expression within F₁ and backcross classes, and 2) no apparent directional parental bias on the F₁ phenotype. Had morphology alone been used as a guide, over 40 percent of the individuals with known hybrid ancestry would have been misclassified as “pure” parental species, and about 25 percent of the backcross individuals would not have been distinguished from F₁‘s. These results exemplify the utility of joint comparisons of morphology and genotypic constitution in studies of natural hybridization, and they emphasize the limitations inherent in describing hybrid classes solely by morphological criteria.     

How do species barriers decay? Concordance and local introgression in mosaic hybrid zones of mussels

The Mytilus complex of marine mussel species forms a mosaic of hybrid zones, found across temperate regions of the globe. This allows us to study ‘replicated’ instances of secondary contact between closely related species. Previous work on this complex has shown that local introgression is both widespread and highly heterogeneous, and has identified SNPs that are outliers of differentiation between lineages. Here, we developed an ancestry‐informative panel of such SNPs. We then compared their frequencies in newly sampled populations, including samples from within the hybrid zones, and parental populations at different distances from the contact. Results show that close to the hybrid zones, some outlier loci are near to fixation for the heterospecific allele, suggesting enhanced local introgression, or the local sweep of a shared ancestral allele. Conversely, genomic cline analyses, treating local parental populations as the reference, reveal a globally high concordance among loci, albeit with a few signals of asymmetric introgression. Enhanced local introgression at specific loci is consistent with the early transfer of adaptive variants after contact, possibly including asymmetric bi‐stable variants (Dobzhansky‐Muller incompatibilities), or haplotypes loaded with fewer deleterious mutations. Having escaped one barrier, however, these variants can be trapped or delayed at the next barrier, confining the introgression locally. These results shed light on the decay of species barriers during phases of contact.     

Morphological and molecular evidence of natural hybridization in Shorea (Dipterocarpaceae)

Shorea (Dipterocarpaceae) is a large genus in which many closely related species often grow together in Southeast Asian lowland tropical rain forests. Many Shorea species share common pollinators, and earlier studies suggested occurrence of interspecific hybridization and introgression. Here, we show morphological and molecular evidence of hybridization between Shorea species. In the census of all the trees of Shorea curtisii, Shorea leprosula, and Shorea parvifolia (>30 cm dbh) within the 164-ha area of Bukit Timah Nature Reserve in Singapore, we found 21 morphologically recognizable hybrid individuals. All of the putative hybrids could be distinguished obviously from the parental species on the basis of vegetative characters. Population genetic analysis of DNA sequences of two nuclear (GapC and PgiC) and chloroplast (trnL-trnF) regions demonstrated that each of the three species had several species-specific mutations. The nuclear sequences of the putative hybrids were heterozygote at all the species-specific sites between two parental species. Hybrid between S. curtisii and S. leprosula was found most, while S. curtisii × S. parvifolia and S. leprosula × S. parvifolia hybrids were also found. Almost no shared polymorphism between populations of the parental species suggests rarity of introgression. The study indicated that natural hybridization between sympatric Shorea species should not be uncommon, but all of the hybrid individuals were F₁, and the post-F₁ hybrids were considerably rare.