Genomic insights into biased allele loss and increased gene numbers after genome duplication in autotetraploid Cyclocarya paliurus

BACKGROUND

Autopolyploidy is a valuable model for studying whole-genome duplication (WGD) without hybridization, yet little is known about the genomic structural and functional changes that occur in autopolyploids after WGD. Cyclocarya paliurus (Juglandaceae) is a natural diploid-autotetraploid species. We generated an allele-aware autotetraploid genome, a chimeric chromosome-level diploid genome, and whole-genome resequencing data for 106 autotetraploid individuals at an average depth of 60 × per individual, along with 12 diploid individuals at an average depth of 90 × per individual.

RESULTS

Autotetraploid C. paliurus had 64 chromosomes clustered into 16 homologous groups, and the majority of homologous chromosomes demonstrated similar chromosome length, gene numbers, and expression. The regions of synteny, structural variation and nonalignment to the diploid genome accounted for 81.3%, 8.8% and 9.9% of the autotetraploid genome, respectively. Our analyses identified 20,626 genes (69.18%) with four alleles and 9191 genes (30.82%) with one, two, or three alleles, suggesting post-polyploid allelic loss. Genes with allelic loss were found to occur more often in proximity to or within structural variations and exhibited a marked overlap with transposable elements. Additionally, such genes showed a reduced tendency to interact with other genes. We also found 102 genes with more than four copies in the autotetraploid genome, and their expression levels were significantly higher than their diploid counterparts. These genes were enriched in enzymes involved in stress response and plant defense, potentially contributing to the evolutionary success of autotetraploids. Our population genomic analyses suggested a single origin of autotetraploids and recent divergence (~ 0.57 Mya) from diploids, with minimal interploidy admixture.

CONCLUSIONS

Our results indicate the potential for genomic and functional reorganization, which may contribute to evolutionary success in autotetraploid C. paliurus.

Genomic Analysis of Plastid-Nuclear Interactions and Differential Evolution Rates in Coevolved Genes across Juglandaceae Species

The interaction between the nuclear and chloroplast genomes in plants is crucial for preserving essential cellular functions in the face of varying rates of mutation, levels of selection, and modes of transmission. Despite this, identifying nuclear genes that coevolve with chloroplast genomes at a genome-wide level has remained a challenge. In this study, we conducted an evolutionary rate covariation analysis to identify candidate nuclear genes coevolving with chloroplast genomes in Juglandaceae. Our analysis was based on 4,894 orthologous nuclear genes and 76 genes across seven chloroplast partitions in nine Juglandaceae species. Our results indicated that 1,369 (27.97%) of the nuclear genes demonstrated signatures of coevolution, with the Ycf1/2 partition yielding the largest number of hits (765) and the ClpP1 partition yielding the fewest (13). These hits were found to be significantly enriched in biological processes related to leaf development, photoperiodism, and response to abiotic stress. Among the seven partitions, AccD, ClpP1, MatK, and RNA polymerase partitions and their respective hits exhibited a narrow range, characterized by dN/dS values below 1. In contrast, the Ribosomal, Photosynthesis, Ycf1/2 partitions and their corresponding hits, displayed a broader range of dN/dS values, with certain values exceeding 1. Our findings highlight the differences in the number of candidate nuclear genes coevolving with the seven chloroplast partitions in Juglandaceae species and the correlation between the evolution rates of these genes and their corresponding chloroplast partitions.

Genomic Insights into Adaptation to Karst Limestone and Incipient Speciation in East Asian Platycarya spp. (Juglandaceae)

When challenged by similar environmental conditions, phylogenetically distant taxa often independently evolve similar traits (convergent evolution). Meanwhile, adaptation to extreme habitats might lead to divergence between taxa that are otherwise closely related. These processes have long existed in the conceptual sphere, yet molecular evidence, especially for woody perennials, is scarce. The karst endemic Platycarya longipes and its only congeneric species, Platycarya strobilacea, which is widely distributed in the mountains in East Asia, provide an ideal model for examining the molecular basis of both convergent evolution and speciation. Using chromosome-level genome assemblies of both species, and whole-genome resequencing data from 207 individuals spanning their entire distribution range, we demonstrate that P. longipes and P. strobilacea form two species-specific clades, which diverged around 2.09 million years ago. We find an excess of genomic regions exhibiting extreme interspecific differentiation, potentially due to long-term selection in P. longipes, likely contributing to the incipient speciation of the genus Platycarya. Interestingly, our results unveil underlying karst adaptation in both copies of the calcium influx channel gene TPC1 in P. longipes. TPC1 has previously been identified as a selective target in certain karst-endemic herbs, indicating a convergent adaptation to high calcium stress among karst-endemic species. Our study reveals the genic convergence of TPC1 among karst endemics and the driving forces underneath the incipient speciation of the two Platycarya lineages.

Genomic Insights into Adaptation to Karst Limestone and Incipient Speciation in East Asian Platycarya spp. (Juglandaceae)

When challenged by similar environmental conditions, phylogenetically distant taxa often independently evolve similar traits (convergent evolution). Meanwhile, adaptation to extreme habitats might lead to divergence between taxa that are otherwise closely related. These processes have long existed in the conceptual sphere, yet molecular evidence, especially for woody perennials, is scarce. The karst endemic Platycarya longipes, and its only congeneric species, P. strobilacea, which is widely distributed in the mountains in East Asia, provide an ideal model for examining the molecular basis of both convergent evolution and speciation. Using chromosome-level genome assemblies of both species, and whole genome resequencing data from 207 individuals spanning their entire distribution range, we demonstrate that P. longipes and P. strobilacea form two species-specific clades, which diverged around 2.09 million years ago. We find an excess of genomic regions exhibiting extreme interspecific differentiation, potentially due to long-term selection in P. longipes, likely contributing to the incipient speciation of the genus Platycarya. Interestingly, our results unveil underlying karst adaptation in both copies of the calcium influx channel gene TPC1 in P. longipes. TPC1 has previously been identified as a selective target in certain karst-endemic herbs, indicating a convergent adaptation to high calcium stress among karst-endemic species. Our study reveals the genic convergence of TPC1 among karst endemics, and the driving forces underneath the incipient speciation of the two Platycarya lineages.

ploidyfrost: Reference-free estimation of ploidy level from whole genome sequencing data based on de Bruijn graphs

Polyploidy is ubiquitous and its consequences are complex and variable. A change of ploidy level generally influences genetic diversity and results in morphological, physiological and ecological differences between cells or organisms with different ploidy levels. To avoid cumbersome experiments and take advantage of the less biased information provided by the vast amounts of genome sequencing data, computational tools for ploidy estimation are urgently needed. Until now, although a few such tools have been developed, many aspects of this estimation, such as the requirement of a reference genome, the lack of informative results and objective inferences, and the influence of false positives from errors and repeats, need further improvement. We have developed ploidyfrost, a de Bruijn graph-based method, to estimate ploidy levels from whole genome sequencing data sets without a reference genome. ploidyfrost provides a visual representation of allele frequency distribution generated using the ggplot2 package as well as quantitative results using the Gaussian mixture model. In addition, it takes advantage of colouring information encoded in coloured de Bruijn graphs to analyse multiple samples simultaneously and to flexibly filter putative false positives. We evaluated the performance of ploidyfrost by analysing highly heterozygous or repetitive samples of Cyclocarya paliurus and a complex allooctoploid sample of Fragaria × ananassa. Moreover, we demonstrated that the accuracy of analysis results can be improved by constraining a threshold such as Cramér's V coefficient on variant features, which may significantly reduce the side effects of sequencing errors and annoying repeats on the graphical structure constructed.

Population-genomic analyses reveal bottlenecks and asymmetric introgression from Persian into iron walnut during domestication

Background

Persian walnut, Juglans regia, occurs naturally from Greece to western China, while its closest relative, the iron walnut, Juglans sigillata, is endemic in southwest China; both species are cultivated for their nuts and wood. Here, we infer their demographic histories and the time and direction of possible hybridization and introgression between them.

Results

We use whole-genome resequencing data, different population-genetic approaches (PSMC and GONE), and isolation-with-migration models (IMa3) on individuals from Europe, Iran, Kazakhstan, Pakistan, and China. IMa3 analyses indicate that the two species diverged from each other by 0.85 million years ago, with unidirectional gene flow from eastern J. regia and its ancestor into J. sigillata, including the shell-thickness gene. Within J. regia, a western group, located from Europe to Iran, and an eastern group with individuals from northern China, experienced dramatically declining population sizes about 80 generations ago (roughly 2400 to 4000 years), followed by an expansion at about 40 generations, while J. sigillata had a constant population size from about 100 to 20 generations ago, followed by a rapid decline.

Conclusions

Both J. regia and J. sigillata appear to have suffered sudden population declines during their domestication, suggesting that the bottleneck scenario of plant domestication may well apply in at least some perennial crop species. Introgression from introduced J. regia appears to have played a role in the domestication of J. sigillata.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02720-z.

Dead-End Hybridization in Walnut Trees Revealed by Large-Scale Genomic Sequence Data

Although hybridization plays a large role in speciation, some unknown fraction of hybrid individuals never reproduces, instead remaining as genetic dead-ends. We investigated a morphologically distinct and culturally important Chinese walnut, Juglans hopeiensis, suspected to have arisen from hybridization of Persian walnut (J. regia) with Asian butternuts (J. cathayensis, J. mandshurica, and hybrids between J. cathayensis and J. mandshurica). Based on 151 whole-genome sequences of the relevant taxa, we discovered that all J. hopeiensis individuals are first-generation hybrids, with the time for the onset of gene flow estimated as 370,000 years, implying both strong postzygotic barriers and the presence of J. regia in China by that time. Six inversion regions enriched for genes associated with pollen germination and pollen tube growth may be involved in the postzygotic barriers that prevent sexual reproduction in the hybrids. Despite its long-recurrent origination and distinct traits, J. hopeiensis does not appear on the way to speciation.

The rediscovery of Caryapoilanei (Juglandaceae) after 63 years reveals a new record from China

Despite having been first published in 1941, Caryapoilanei (A.Chev.) J.-F.Leroy is only known from three collections in Vietnam, Laos and Thailand. It has not been recollected since then and was long suspected to have become extinct through repeated deforestation events. Here, we report the rediscovery, and meanwhile the first new record in China, of this extremely rare gigantic hickory species at Yunnan province 63 years after its last collection in 1958. Three small patchy subpopulations were found with a total of about 50 adult trees having diameter at breast height (DBH) larger than 60 cm, together with some seedlings and saplings, but the fruit set was low. Based on new and fresh material, we present a revised morphological description of C.poilanei, and an updated distribution map for the species. In addition, we also provide a key for the hickories in China. Lastly, we suggest C.poilanei should be listed as a Critically Endangered (CR) species according to the latest IUCN Red List Categories and Criteria.

Do annual and perennial populations of an insect-pollinated plant species differ in mating system?

BACKGROUND AND AIMS

Theory predicts that outcrossing should be more prevalent among perennials than annuals, a pattern confirmed by comparative evidence from diverse angiosperm families. However, intraspecific comparisons between annual and perennial populations are few because such variation is uncommon among flowering plants. Here, we test the hypothesis that perennial populations outcross more than annual populations by investigating Incarvillea sinensis, a wide-ranging insect-pollinated herb native to China. The occurrence of both allopatric and sympatric populations allows us to examine the stability of mating system differences between life histories under varying ecological conditions.

METHODS

We estimated outcrossing rates and biparental inbreeding in 16 allopatric and five sympatric populations in which both life histories coexisted using 20 microsatellite loci. In each population we measured height, branch number, corolla size, tube length and herkogamy for ~30 individuals. In a sympatric population, we recorded daily flower number, pollinator visitation and the fruit and seed set of annual and perennial plants.

KEY RESULTS

As predicted, outcrossing rates (t) were considerably higher in perennial (mean = 0.76) than annual (mean = 0.09) populations. This difference in mating system was also maintained at sympatric sites where plants grew intermixed. In both allopatric and sympatric populations the degree of herkogamy was consistently larger in outcrossing than selfing plants. Perennials were more branched, with more and larger flowers than in annuals. In a sympatric population, annuals had a significantly higher fruit and seed set than perennials.

CONCLUSIONS

Genetically based differences in herkogamy between annuals and perennials appear to play a key role in governing outcrossing rates in populations, regardless of variation in local ecological conditions. The maintenance of mating system and life history trait differentiation between perennial and annual populations of I. sinensis probably results from correlated evolution in response to local environmental conditions.

Recent demographic histories of temperate deciduous trees inferred from microsatellite markers

Background

Accurate inference of demographic histories for temperate tree species can aid our understanding of current climate change as a driver of evolution. Microsatellites are more suitable for inferring recent historical events due to their high mutation rates. However, most programs analyzing microsatellite data assume a strict stepwise mutation model (SMM), which could cause false detection of population shrinkage when microsatellite mutation does not follow SMM.

Results

This study aims to reconstruct the recent demographic histories of five cool-temperate tree species in Eastern Asia, Quercus mongolica, Q. liaotungensis, Juglans cathayensis, J. mandshurica and J. ailantifolia, by using 19 microsatellite markers with two methods considering generalized stepwise mutation model (GSM) (MIGRAINE and VarEff). Both programs revealed that all the five species experienced expansions after the Last Glacial Maximum (LGM). Within butternuts, J. cathayensis experienced a more serious bottleneck than the other species, and within oaks, Q. mongolica showed a moderate increase in population size and remained stable after the expansion. In addition, the point estimates of the multistep mutation proportion in the GSM model (p_GSM) for all five species were between 0.50 and 0.65, indicating that when inferring population demographic history of the cool-temperate forest species using microsatellite markers, it is better to assume a GSM rather than a SMM.

Conclusions

This study provides the first direct evidence that five cool-temperate tree species in East Asia have experienced expansions after the LGM with microsatellite data. Considering the mutation model of microsatellite has a vital influence on demographic inference, combining multiple programs such as MIGRAINE and VarEff can effectively reduce errors caused by inappropriate model selection and prior setting.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01805-w.

Phylogenomics Reveals an Ancient Hybrid Origin of the Persian Walnut

Persian walnut (Juglans regia) is cultivated worldwide for its high-quality wood and nuts, but its origin has remained mysterious because in phylogenies it occupies an unresolved position between American black walnuts and Asian butternuts. Equally unclear is the origin of the only American butternut, J. cinerea. We resequenced the whole genome of 80 individuals from 19 of the 22 species of Juglans and assembled the genome of its relatives Pterocarya stenoptera and Platycarya strobilacea. Using phylogenetic-network analysis of single-copy nuclear genes, genome-wide site pattern probabilities, and Approximate Bayesian Computation, we discovered that J. regia (and its landrace J. sigillata) arose as a hybrid between the American and the Asian lineages and that J. cinerea resulted from massive introgression from an immigrating Asian butternut into the genome of an American black walnut. Approximate Bayesian Computation modeling placed the hybrid origin in the late Pliocene, ∼3.45 My, with both parental lineages since having gone extinct in Europe.

Demographically idiosyncratic responses to climate change and rapid Pleistocene diversification of the walnut genus Juglans (Juglandaceae) revealed by whole-genome sequences

Whether species demography and diversification are driven primarily by extrinsic environmental changes such as climatic oscillations in the Quaternary or by intrinsic biological interactions like coevolution between antagonists is a matter of active debate. In fact, their relative importance can be assessed by tracking past population fluctuations over considerable time periods. We applied the pairwise sequentially Markovian coalescent approach on the genomes of 11 temperate Juglans species to estimate trajectories of changes in effective population size (N_e ) and used a Bayesian-coalescent based approach that simultaneously considers multiple genomes (G-PhoCS) to estimate divergence times between lineages. N_e curves of all study species converged 1.0 million yr ago, probably reflecting the time when the walnut genus last shared a common ancestor. This estimate was confirmed by the G-PhoCS estimates of divergence times. But all species did not react similarly to the dramatic climatic oscillations following early Pleistocene cooling, so the timing and amplitude of changes in N_e differed among species and even among conspecific lineages. The population histories of temperate walnut species were not driven by extrinsic environmental changes alone, and a key role was probably played by species-specific factors such as coevolutionary interactions with specialized pathogens.

Phylogeographic breaks within Asian butternuts indicate the existence of a phytogeographic divide in East Asia

East Asia has been hypothesized to be subdivided into two distinct northern and southern areas, separated by a band of dry climate that was far more severe in the early Tertiary but still exists today. However, this biogeographic hypothesis has rarely been tested using a molecular phylogeographic approach. We genotyped 70 populations throughout the distributional range of Asian butternuts (Juglans section Cardiocaryon) using eight chloroplast DNA regions, one single-copy nuclear gene, and 17 nuclear microsatellite loci, supplemented with paleodistribution modeling of the major genetic clades. The genetic data consistently identified two clades, one northern, comprising Juglans mandshurica and Juglans ailantifolia, and one southern, comprising Juglans cathayensis. The two clades diverged through climate-induced vicariance of an ancestral northern range during the mid-Miocene and remained mostly separate thereafter, with geographical isolation of the Japanese Islands and refugial isolation or secondary contacts in the late Pleistocene producing further subdivision within the northern clade. But beyond all that, we also discovered a role of environmental adaptation in maintaining and/or reinforcing the north-south divergence. Asian butternuts offer a strong case for the existence of a biogeographic divide between the northern and southern parts of East Asia during the Neogene and into the Pleistocene.

Seasonal variation in the mating system of a selfing annual with large floral displays

BACKGROUND AND AIMS

Flowering plants display considerable variation in mating system, specifically the relative frequency of cross- and self-fertilization. The majority of estimates of outcrossing rate do not account for temporal variation, particularly during the flowering season. Here, we investigated seasonal variation in mating and fertility in Incarvillea sinensis (Bignoniaceae), an annual with showy, insect-pollinated, 'one-day' flowers capable of delayed selfing. We examined the influence of several biotic and abiotic environmental factors on day-to-day variation in fruit set, seed set and patterns of mating.

METHODS

We recorded daily flower number and pollinator abundance in nine 3 × 3-m patches in a population at Mu Us Sand land, Inner Mongolia, China. From marked flowers we collected data on daily fruit and seed set and estimated outcrossing rate and biparental inbreeding using six microsatellite loci and 172 open-pollinated families throughout the flowering period.

KEY RESULTS

Flower density increased significantly over most of the 50-d flowering season, but was associated with a decline in levels of pollinator service by bees, particularly on windy days. Fruit and seed set declined over time, especially during the latter third of the flowering period. Multilocus estimates of outcrossing rate were obtained using two methods (the programs MLTR and BORICE) and both indicated high selfing rates of ∼80 %. There was evidence for a significant increase in levels of selfing as the flowering season progressed and pollinator visitation declined. Biparental inbreeding also declined significantly as the flowering season progressed.

CONCLUSIONS

Temporal variation in outcrossing rates may be a common feature of the mating biology of annual, insect-pollinated plants of harsh environments but our study is the first to examine seasonal mating-system dynamics in this context. Despite having large flowers and showy floral displays, I. sinensis attracted relatively few pollinators. Delayed selfing by corolla dragging largely explains the occurrence of mixed mating in I. sinensis, and this mode of self-fertilization probably functions to promote reproductive assurance when pollinator service is limited by windy environmental conditions and competition from co-occurring flowering plants.

Evolutionary history of a widespread tree species Acer mono in East Asia

East Asia has the most diverse temperate flora in the world primarily due to the lack of Pleistocene glaciation and the geographic heterogeneity. Although increasing phylogeography studies in this region provided more proofs in this issue, discrepancies and uncertainty still exist, especially in northern temperate deciduous broad-leaved and coniferous mixed forest region (II). And a widespread plant species could reduce the complexity to infer the relationship between diversity and physiographical pattern. Hence, we studied the evolution history of a widespread temperate tree, Acer mono, populations in region II and the influence of physiographic patterns on intraspecific genetic diversity. Analyses of chloroplast sequences and nuclear microsatellites indicated high levels of genetic diversity. The diversity distribution was spatially heterogeneous and a latitudinal cline existed in both markers. The spatial distribution pattern between genetic diversity within A. mono and the diversity at species level was generally consistent. Western subtropical evergreen broad-leaved forest subregion (IVb) had a unique ancient chloroplast clade (CP3) and a nuclear gene pool (GP5) with dominance indicating the critical role of this area in species diversification. Genetic data and ecological niche model results both suggested that populations in region II disappeared during the last glacial maximum (LGM) and recovered from south of Changbai Mt. and the Korean Peninsula. Two distribution centers were likely during the LGM, one in the north edge of warm temperate deciduous broad-leaved forest region (III) and another in the south edge of region III. This was reflected by the genetic pattern with two spatially independent genetic groups. This study highlights the key role of region III in sustaining genetic diversity in the northern range and connecting diversity between southern and northern range. We elucidated the diversity relationship between vegetation regions which could facilitate the understanding of biodiversity origin and maintenance in East Asia.

Development of polymorphic microsatellite markers for Incarvillea sinensis (Bignoniaceae)

PREMISE OF THE STUDY

Microsatellite markers were developed for Incarvillea sinensis var. sinensis (Bignoniaceae), an annual herb endemic to Inner Mongolia, to study the degree to which delayed self-fertilization is favored.

METHODS AND RESULTS

Eight polymorphic primer sets were isolated and characterized in two Inner Mongolia populations of I. sinensis var. sinensis with a relatively simple and fast subcloning method. Numbers of alleles per locus ranged from 2 to 7, with observed and expected heterozygosities ranging from 0 to 0.261 and from 0 to 0.778, respectively.

CONCLUSIONS

These markers will be useful for future studies of self-fertilization adaptability in I. sinensis var. sinensis.

Nuclear and chloroplast DNA phylogeography reveal two refuge areas with asymmetrical gene flow in a temperate walnut tree from East Asia

• Recently, there has been a debate about whether the temperate forests of East Asia merged or fragmented during glacial periods in the Pleistocene. Here, we tested these two opposing views through phylogeographical studies of the temperate-deciduous walnut tree, Juglans mandshurica (Juglandaceae) in northern and northeastern China, as well as Japan and Korea. • We assessed the genetic structure of 33 natural populations using 10 nuclear microsatellite loci and seven chloroplast DNA (cpDNA) fragments. • The cpDNA data showed the complete fixation of two different haplotype lineages in northeastern vs northern populations. This pronounced phylogeographic break was also indicated by nuclear microsatellite data, but there were disparities regarding individual populations. Among those populations fixed for haplotype A (the northeastern group), three were clustered in the northern group and four showed evidence of mixed ancestry based on microsatellite data. • Our results support the hypothesis that two independent refugia were maintained across the range of J. mandshurica in the north of China during the last glacial maximum, contrary to the inference that all temperate forests migrated to the south (25-30°N). The discordance between the patterns revealed by cpDNA and microsatellite data indicate that asymmetrical gene flow has occurred between the two refugia.

Variation in floral sex allocation and reproductive success in sequentially flowering inflorescence of Corydalis remota var. lineariloba (Fumariaceae)

In hermaphroditic plants, female reproductive success often varies among different positions within an inflorescence. However, few studies have evaluated the relative importance of underlying causes such as pollen limitation, resource limitation or architectural effect, and few have compared male allocation. During a 2-year investigation, we found that female reproductive success of an acropetally flowering species, Corydalis remota Fisch. ex Maxim. var. lineariloba Maxim. was significantly lower in the upper late developing flowers when compared with the lower early flowers. Supplementation with outcross pollen did not improve female reproductive success of the upper flowers, while removal of the lower developing fruits significantly increased female reproductive success of the upper flowers in both years, evidencing resource limitation of the upper flowers. Female production in upper flowers was greatly improved by simultaneous pollen supplementation of the upper flowers and removal of the lower fruits, suggesting that, when resources are abundant, pollen may limit the female reproductive success of the upper flowers. The less seed mass in the upper flowers didn't increase in all treatments due to architecture. In the upper flowers, ovule production was significantly lower and the pollen : ovule ratio was significantly higher. These results suggest that male-biased sex allocation in the upper flowers may lead to increased male reproductive success, whereas the lower flowers have higher female reproductive success.

Mating patterns and pollen dispersal in a heterodichogamous tree, Juglans mandshurica (Juglandaceae)

Mating patterns in heterodichogamous species are generally considered to be disassortative between flowering morphs, but this hypothesis has hitherto not been vigorously tested. Here, mating patterns and pollen dispersal were studied in Juglans mandshurica, a heterodichogamous wind-pollinated species that is widely distributed in northern and north-eastern China. Paternity analyses carried out on 11 microsatellite loci were used to estimate morph-specific rates of outcrossing and disassortative mating. Pollen dispersal and genetic structure were also investigated in the population under study. The mating pattern of J. mandshurica was highly outcrossing and disassortative. Pairwise values of intramorph relatedness were much higher than those of intermorph relatedness, and a low level of biparental inbreeding was detected. There was no significant difference in outcrossing and disassortative mating rates between the two morphs. The effective pollen dispersal distribution showed an excess of near-neighbor matings, and most offspring of individual trees were sired by one or two nearby trees. These results corroborate the previous suggestion that mating in heterodichogamous plant species is mainly disassortative between morphs, which not only prevents selfing but also effectively reduces intramorph inbreeding.