Strong reproductive isolation despite occasional hybridization between a widely distributed and a narrow endemic Rhododendron species

Synthesis and Scope of the Role of Postmating Prezygotic Isolation in Speciation

How barriers to gene flow arise and are maintained are key questions in evolutionary biology. Speciation research has mainly focused on barriers that occur either before mating or after zygote formation. In comparison, postmating prezygotic (PMPZ) isolation-a barrier that acts after gamete release but before zygote formation-is less frequently investigated but may hold a unique role in generating biodiversity. Here we discuss the distinctive features of PMPZ isolation, including the primary drivers and molecular mechanisms underpinning PMPZ isolation. We then present the first comprehensive survey of PMPZ isolation research, revealing that it is a widespread form of prezygotic isolation across eukaryotes. The survey also exposes obstacles in studying PMPZ isolation, in part attributable to the challenges involved in directly measuring PMPZ isolation and uncovering its causal mechanisms. Finally, we identify outstanding knowledge gaps and provide recommendations for improving future research on PMPZ isolation. This will allow us to better understand the nature of this often-neglected reproductive barrier and its contribution to speciation.

Partial and asymmetrical reproductive isolation between two sympatric tropical shrub species: Cnidoscolus aconitifolius and C. souzae (Euphorbiaceae)

Reproductive isolation is conferred by several barriers that occur at different stages of reproduction. Comprehensive reviews on the topic have identified that barriers occurring prior to zygote formation are often stronger than those that occur afterward. However, the overrepresentation of temperate perennial herbs in the current literature precludes any generalization of this pattern to plants that present other life forms and patterns of distribution. Here, we assessed reproductive isolation barriers and their absolute contribution to reproductive isolation and asymmetry in Cnidoscolus aconitifolius and C. souzae, two closely related tropical shrub species that co-occur on the Yucatan peninsula. The reproductive barriers assessed were phenological mismatch, pollinator differentiation, pollen-pistil incompatibility (three pre-zygotic barriers), fruit set failure, and seed unviability (post-zygotic barriers). Reproductive isolation between the study species was found to be complete in the direction C. aconitifolius to C. souzae, but only partial in the opposite direction. One post-zygotic barrier was the strongest example. Most barriers, particularly the pre-zygotic examples, were asymmetrical and predicted the direction of heterospecific pollen flow and hybrid formation from C. souzae to C. aconitifolius. Both parental species, as well as the hybrids, were diploid and had a chromosome number 2n = 36. More studies with tropical woody perennials are required to fully determine whether this group of plants consistently shows stronger post-zygotic barriers.

Root microbiota alters response to root rot in Rhododendron delavayi Franch

Root microbiota have a significant effect on plant health. However, the role of root microbiota in the resistance of Rhododendron against root rot is not known. In this study, we employed amplicon 16S and ITS sequencing to investigate the bacterial and fungal communities associated with four distinct niches (bulk soil, rhizosphere, rhizoplane, and endosphere) of both healthy and diseased Rhododendron plants in the Baili Rhododendron nature reserve in China. The amplicon data analysis identified 182 bacterial genera and 141 fungal genera that were impacted by root rot across all niches. Specifically, the rhizoplane appeared to exert a selective gating effect, resulting in a reduction in the complexity of bacterial communities, but not fungal communities, in wild Rhododendron delavayi Franch roots. Nevertheless, the stress induced by root rot led to alterations in the root microbiota and compromised the gating function of the rhizoplane, thereby significantly increasing the complexity of the bacterial community within the plant root. In the root tissue following root rot outbreak, the relative abundance of the pathogenic species Pezicula brunnea and Diaporthe helianthi was enriched by as much as 6.13% and 1.71%, respectively. These findings provide novel insights into the contribution of enrichment of root-associated microbiota to wild plant hosts under the disease stress of root rot. The root rot-causing pathogenic fungi may interact with beneficial bacteria and induce plants to send out "cry for help" signals, which may encourage the specific assembly of microbiota. In the Rhododendron delavayi Franch root microbiota, we found 23 potentially beneficial microbes. Notably, certain beneficial bacteria, such as Sporolactobacillus and Stenotrophomonas, were found to accumulate in the rhizoplane and endosphere under root rot disease stress. Overall, our results lend support to our hypothesis that Rhododendron recruits protective microbes as a strategy to suppress root rot outbreaks. Future endeavors in isolating beneficial microbes capable of mitigating root rot have the potential to enhance plant resilience against root diseases.

Allelopathic potential and allelochemical composition in different soil layers of Rhododendron delavayi forest, southwest China

Forests habituated by Rhododendron delavayi often lack understory vegetation, which could possibly be a consequence of allelopathy. It is a phenomenon by virtue of which certain plant species produce allelochemicals that affect the growth and behavior of surrounding plants. To elucidate the allelopathic potential and allelochemicals present in the different layers of a R. delavayi forest, extracts obtained from three layers of the forest were used for seed germination bioassays and gas chromatography-mass spectrometry (GC–MS) analysis. Aqueous extracts of the litter and humus layers significantly inhibited the seed germination of R. delavayi, Festuca arundinacea, and Lolium perenne, with the litter layer causing the strongest inhibitory effect. A total of 26 allelochemicals were identified in the litter, humus, and soil layers by GC–MS analysis. The primary allelochemicals in the soil and humus layer were organic acids, while the main allelochemicals in the litter layer were phenolic acids. The redundancy analysis revealed the significance of total nitrogen (TN) and relative water content (RWC) in explaining the distribution of the allelochemicals. The results indicated that the litter layer exerted the maximum allelopathic effect due to presence of maximum amount of allelochemical especially the phenolic acids.

Unusual patterns of hybridization involving two alpine Salvia species: Absence of both F1 and backcrossed hybrids

Natural hybridization plays an important role in speciation; however, we still know little about the mechanisms underlying the early stages of hybrid speciation. Hybrid zones are commonly dominated by F₁s, or backcrosses, which impedes further speciation. In the present study, morphological traits and double digest restriction-site associated DNA sequencing (ddRAD-seq) data have been used to confirm natural hybridization between Salvia flava and S. castanea, the first case of identification of natural hybridization using combined phenotypic and molecular evidence in the East Asian clade of Salvia. We further examined several reproductive barriers in both pre-zygotic and post-zygotic reproductive stages to clarify the causes and consequences of the hybridization pattern. Our results revealed that reproductive isolation between the two species was strong despite the occurrence of hybridization. Interestingly, we found that most of the hybrids were likely to be F₂s. This is a very unusual pattern of hybridization, and has rarely been reported before. The prevalence of geitonogamy within these self-compatible hybrids due to short distance foraging by pollinators might explain the origin of this unusual pattern. F₂s can self-breed and develop further, therefore, we might be witnessing the early stages of hybrid speciation. Our study provides a new case for understanding the diversification of plants on the Qinghai-Tibet Plateau.

Microhabitat and Pollinator Differentiation Drive Reproductive Isolation between Two Sympatric Salvia Species (Lamiaceae)

Evaluation of multiple barriers contributing to reproductive isolation between sympatric plant species is key to understanding the mechanism of their coexistence; however, such investigations in biodiversity hotspots are still rare. In this study, we investigated and compared geography, microhabitat, phenology, flora, and pollinators, in addition to pollen–pistil interactions, seed production, and seed germination of the closely related sympatric Salvia digitaloides and S. flava on Yulong Snow Mountain, Southwestern Yunnan, China. The geographic distribution of these species overlapped, but their adaptation to physical and chemical properties of soil microhabitats differed. They shared the same flowering time but differed in flower size, style length, nectar volume, sugar concentration, and flower longevity. Both species shared bumblebees as effective pollinators, but flower constancy for the two species was relatively strong. Pollen tube growth, seed production, and seed germination were lower in interspecific than in intraspecific crosses. Our study suggested that microhabitat and pollinator isolation acted as the most important isolating barriers in maintaining the coexistence of the two Salvia species. Our study also highlighted that post-pollination barriers play an important role in preventing the gene flow between these two Salvia species.

A comprehensive review on the documented characteristics of four Reticulitermes termites (Rhinotermitidae, Blattodea) of China

Termites are known as social insects worldwide. Presently in China 473 species, 44 genera and 4 families of termites have been reported. Of them, 111 Reticulitermes species are widely spread in different zones of China. The dispersion flight season of these Chinese Reticulitermes species are usually started from February to June, but in some regions different species are distributed, sharing their boundaries and having overlapping flight seasons. These reasons become important sources of hybridization between two different heterospecific populations of termites. It was confirmed that the fertilized eggs and unfertilized eggs of some Reticulitermes termites have the capacity of cleavage. While the unfertilized eggs of R. aculabialis, R. chinensis and R. labralis cleaved normally and the only R. aculabialis unfertilized eggs develop in embryos. While, the R. flaviceps and R. chinensis were observed with their abnormal embryonic development, and not hatching of eggs parthenogenetically. They were reported more threatening to Chinese resources as they propagate with parthenogenesis, hybridization and sexual reproduction. Eggshell and macrophiles of eggs play important roles in species identification and control. Although, they are severe pests and cause a wide range of damages to wooden structures and products in homes, buildings, building materials, trees, crops, and forests in China's Mainland.

Natural hybridization among three Rhododendron species (Ericaceae) revealed by morphological and genomic evidence

BACKGROUND

Natural hybridization can influence the adaptive response to selection and accelerate species diversification. Understanding the composition and structure of hybrid zones may elucidate patterns of hybridization processes that are important to the formation and maintenance of species, especially for taxa that have experienced rapidly adaptive radiation. Here, we used morphological traits, ddRAD-seq and plastid DNA sequence data to investigate the structure of a Rhododendron hybrid zone and uncover the hybridization patterns among three sympatric and closely related species.

RESULTS

Our results show that the hybrid zone is complex, where bi-directional hybridization takes place among the three sympatric parental species: R. spinuliferum, R. scabrifolium, and R. spiciferum. Hybrids between R. spinuliferum and R. spiciferum (R. ×duclouxii) comprise multiple hybrid classes and a high proportion of F₁ generation hybrids, while a novel hybrid taxon between R. spinuliferum and R. scabrifolium dominated the F₂ generation, but no backcross individuals were detected. The hybrid zone showed basically coincident patterns of population structure between genomic and morphological data.

CONCLUSIONS

Natural hybridization exists among the three Rhododendron species in the hybrid zone, although patterns of hybrid formation vary between hybrid taxa, which may result in different evolutionary outcomes. This study represents a unique opportunity to dissect the ecological and evolutionary mechanisms associated with adaptive radiation of Rhododendron species in a biodiversity hotspot.

Chromosome-level genome assembly and population genetic analysis of a critically endangered rhododendron provide insights into its conservation

Rhododendrons are woody plants, famous throughout the world as having high horticultural value. However, many wild species are currently threatened with extinction. Here, we report for the first time a high-quality, chromosome-level genome of Rhododendron griersonianum, which has contributed to approximately 10% of all horticultural rhododendron varieties but which in its wild form has been evaluated as critically endangered. The final genome assembly, which has a contig N50 size of approximately 34 M and a total length of 677 M, is the highest-quality genome sequenced within the genus to date, in part due to its low heterozygosity (0.18%). Identified repeats constitute approximately 57% of the genome, and 38 280 protein-coding genes were predicted with high support. We further resequenced 31 individuals of R. griersonianum as well as 30 individuals of its widespread relative R. delavayi, and performed additional conservation genomic analysis. The results showed that R. griersonianum had lower genetic diversity (θ = 2.58e-3; π = 1.94e-3) when compared not only to R. delavayi (θ = 11.61e-3, π = 12.97e-3), but also to most other woody plants. Furthermore, three severe genetic bottlenecks were detected using both the Stairway plot and fastsimcoal2 analysis, which are thought to have occurred in the late Middle Pleistocene and the Last Glacial Maximum (LGM) period. After these bottlenecks, R. griersonianum recovered and maintained a constant effective population size (>25 000) until now. Intriguingly, R. griersonianum has accumulated significantly more deleterious mutations in the homozygous state than R. delavayi, and several deleterious mutations (e.g., in genes involved in the response to heat stress) are likely to have harmed the adaptation of this plant to its surroundings. This high-quality, chromosome-level genome and the population genomic analysis of the critically endangered R. griersonianum will provide an invaluable resource as well as insights for future study in this species to facilitate conservation and in the genus Rhododendron in general.

Conservation Genomics of a Threatened Rhododendron: Contrasting Patterns of Population Structure Revealed From Neutral and Selected SNPs

Though it is well-acknowledged that next generation sequencing (NGS) technologies can provide further insights into plant conservation management than traditional molecular markers, studies employing NGS to address conservation genomics and subsequent conservation strategies for threatened plants are still rare. Rhododendron is the largest genus of woody plants in China, and many species are threatened, however, to date there has been no conservation genetic research using NGS in this genus. In the present study, we investigated the conservation genetics of R. cyanocarpum, a threatened species endemic to the Cangshan Mountains in Yunnan, China, using a double digest restriction-site-associated DNA-sequencing (ddRAD-seq) approach. Due to the availability of sufficient SNPs, we were able to distinguish between neutral and putatively selected SNPs and were able to further investigate the genetic diversity, population structure, and differentiation in R. cyanocarpum, as well as make an estimation of its demographic history. A total of 6,584 SNPs were obtained, of which 5,729 were neutral (detected using Tajima’s D). In terms of the 5,729 neutral SNPs, R. cyanocarpum had a higher genetic diversity (π = 0.0702 ± 0.0017, H_e = 0.0675 ± 0.0016) than other plant species assessed using Rad-seq methods, while population differentiation (F_st from 0.0314 to 0.0452) was weak. Interestingly, contrasting patterns of population structure were revealed from all neutral and selected SNPs, with distinct genetic clusters forming for all SNPs and neutral SNPs, but no distinct subgroups for selected ones. Moreover, we were able to detect changes in effective population size (N_e) of R. cyanocarpum from 150,000 years ago, including a bottleneck event ca. 60,000 years ago, followed by recovery of N_e over a short period, and a subsequent gradual decline in N_e to date. Implications for conserving R. cyanocarpum based on these main results are then discussed.

Comparative population genetic analyses suggest hybrid origin of R hododendron pubicostatum, an endangered plant species with extremely small populations endemic to Yunnan, China

Gene flow between sympatric congeneric plants is thought to be very common and may pose serious threats to endangered species. In the present study, we evaluate the genetic diversity and divergence of three sympatric Rhododendron species in Jiaozi Mountain using newly developed microsatellites through the Illumina MiSeq sequencing approach. Genetic diversity of all three Rhododendron species studied was moderate in comparison to genetic parameters previously reported from species of this genus. Interestingly, genetic structure analysis of the three species identified a possible hybrid origin of the threatened Rh. pubicostatum. This sympatry should be considered a unimodal hybrid zone, since R h. pubicostatum is predominant here. Unimodal hybrid zones are uncommon in Rhododendron, despite the fact that hybridization frequently occurs in the genus. Issues pertaining to the conservation of R h. pubicostatum resulting from admixture of genetic material from its parental species are discussed.

A Genomewide Scan for Genetic Structure and Demographic History of Two Closely Related Species, Rhododendron dauricum and R. mucronulatum (Rhododendron, Ericaceae)

Understanding the processes of divergence and speciation is an important task for evolutionary research, and climate oscillations play a pivotal role. We estimated the genetic structure and demographic history of two closely related species of Rhododendron, R. dauricum, and R. mucronulatum, distributed in northeastern China using 664,406 single nucleotide polymorphic loci of specific-locus amplified fragment sequencing (SLAF-seq) and 4 chloroplast DNA (cpDNA) fragments, sampling 376 individuals from 39 populations of these two species across their geographic distributions. The geographical distribution of cpDNA haplotypes revealed that R. dauricum and R. mucronulatum have different spatial genetic structures and haplotype diversity. Analysis of molecular variance (AMOVA) results showed that these two species have significant genetic differentiation and that the phylogeny demonstrates that these two species clustered a monophyletic group based on SLAF data, respectively, but not in cpDNA data. The evidence of significant gene flow was also detected from R. mucronulatum to R. dauricum. A deep divergence between the two species was observed and occurred during the early Oligocene. The niche models showed that the two species have different demographic histories. Thus, our results imply that geography and climate changes played important roles in the evolutionary process of R. dauricum and R. mucronulatum, and although there was an interspecific gene flow, the divergence was maintained by natural selection.

Homoploid hybridization of plants in the Hengduan mountains region

The Hengduan Mountains Region (HMR) is a major global biodiversity hotspot. Complex tectonic and historical climatic conditions created opportunities for natural interspecific hybridization. Likewise, anthropogenic disturbance potentially raises the frequency of hybridization. Among species studies to date, the frequency of homoploid hybridization appears in the HMR. Of nine taxa in which natural hybridization has been detected, three groups are involved in homoploid hybrid speciation, and species pairs from the remaining six genera suggest that continuous gene flow occurs in hybrid zones. Reproductive isolation may greatly affect the dynamic and architecture of hybrid zones in the HMR. Asymmetrical hybridization and introgression can primarily be attributed to both prezygotic and postzygotic barriers. The frequent observation of such asymmetry may imply that reproductive barrier contributes to maintaining species boundaries in the alpine region. Ecological isolations with environmental disturbance may promote breeding barriers between parental species and hybrids. Hybrid zones may be an important phase for homoploid hybrid speciation. Hybrid zones potentially provided abundant genetic resources for the diversification of the HMR flora. The ecological and molecular mechanisms of control and mediation for natural hybridization will help biologists to understand the formation of biodiversity in the HMR. More researches from ecological and molecular aspects were required in future studies.

Conservation implications of asymmetric introgression and reproductive barriers in a rare primrose species

BACKGROUND

Primula is a large genus of flowering herbs well known for their heterostyly. Currently few natural hybrids are known and reproductive barriers in this genus in the wild have received little attention. However, there is instance of hybridization between rare and widely-spread species, and conservation implications of such situation is poorly understood. In the present study, we investigated hybridization patterns and reproductive barriers between a wide spread species, Primula poissonii and a rare species P. anisodora, of which only three populations are currently known.

RESULTS

Pollinator-mediated reproductive isolation was strong between parental species but not significant between hybrids and parental species. Hand pollination experiments showed significant reduction of both fruit- and seed-set for heterospecific pollination as compared with conspecific pollination for both parental species. Furthermore, hybrids had higher fruit- and seed-set when pollinated with P. anisodora pollen as opposed to P. poissonii pollen. Microsatellites identified backcrosses to P. anisodora in two of the three populations of P. anisodora, and additionally more individuals of P. anisodora showed introgression from P. poissonii than vice versa.

CONCLUSIONS

These results provide evidence for potential genetic swamping of the P. anisodora populations, which could pose a serious threat for this locally endemic species.

Molecular evidence for asymmetric hybridization in three closely related sympatric species

Natural hybridization is common in plants and results in different evolutionary consequences to hybridizing species. Pre- and post-zygotic reproductive isolating barriers can impede hybridization between closely related species to maintain their species integrity. In Northwest Yunnan, three Ligularia species (Ligularia cyathiceps, L. duciformis and L. yunnanensis) and four types of morphologically intermediate individuals were discovered growing together in an area subject to human disturbance. In this study, we used three low-copy nuclear loci to test the natural hybridization hypothesis and the hybridization direction was ascertained by three chloroplast DNA fragments. The results indicated there were two hybridization groups at the study site, L. cyathiceps × L. duciformis and L. duciformis × L. yunnanensis, and two types of morphologically intermediate individuals were produced by L. cyathiceps and L. duciformis, and another two types were produced by L. duciformis and L. yunnanensis, while no hybrids between L. cyathiceps and L. yunnanensis were observed. Both hybridizing groups showed bidirectional but asymmetric hybridization and the factors influencing the symmetry are discussed. Most hybrids produced by the two hybridization groups seemed to be F1 generation. Hybrids with different morphologies within the same hybridization group showed similar genetic components. The results suggest that although human disturbance may promote natural hybridization among the three Ligularia species bringing them together, hybrids are limited to F1s and therefore species boundaries might be maintained.

Characterization of sympatric Platanthera bifolia and Platanthera chlorantha (Orchidaceae) populations with intermediate plants

Platanthera bifolia and P. chlorantha are terrestrial and rewarding orchids with a wide Eurasian distribution. Although genetically closely related, they exhibit significant morphological, phenological and ecological differences that maintain reproductive isolation between the species. However, where both species co-occur, individuals with intermediate phenotypic traits, often considered as hybrids, are frequently observed. Here, we combined neutral genetic markers (AFLPs), morphometrics and floral scent analysis (GC-MS) to investigate two mixed Platanthera populations where morphologically intermediate plants were found. Self-pollination experiments revealed a low level of autogamy and artificial crossings combined with assessments of fruit set and seed viability, showed compatibility between the two species. The results of the genetic analyses showed that morphologically intermediate plants had similar genetic patterns as the P. bifolia group. These results are corroborated also by floral scent analyses, which confirmed a strong similarity in floral scent composition between intermediate morphotypes and P. bifolia. Therefore, this study provided a much more detailed picture of the genetic structure of a sympatric zone between two closely allied species and supports the hypothesis that intermediate morphotypes in sympatry could reflect an adaptive evolution in response to local pollinator-mediated selection.

Study on Reproductive Biology of Rhododendron longipedicellatum: A Newly Discovered and Special Threatened Plant Surviving in Limestone Habitat in Southeast Yunnan, China

Rhododendron longipedicellatum is a narrow endemic species and a subject of urgent demand in the domestic market and overseas. Its fascinating shapes, brilliantly gilvous flowers, and unusual flowering time endow this species with extremely high ornamental value. However, only five wild populations of R. longipedicellatum surviving in limestone habitat have been found through elaborate field investigation, and the number of the populations decreases further or is even confronted with risk of extinction due to the damage of human activities. To enhance the protection and utilization of R. longipedicellatum, this study systematically investigated several important aspects of reproductive biology, including floral syndrome, pollen viability and stigma receptivity, petal color reflectance, breeding system, and pollination biology. The results demonstrated that arched styles not only create obvious herkogamy that avoide self-pollination, but also effectively reduce rain damage to the intrinsic characteristics of the stigma surface secretions, promoting the female fitness of R. longipedicellatum in poor weather. Pollen viability maintained a high level over the flowering period. The reflectance spectrum of petals had two peaks at wavelengths of 360 and 580 nm. Tests of OCI, P/O and artificial pollination all indicated that R. longipedicellatum was self-compatible and that the breeding system was mixed mating. Geitonogamy mediated by Bombus braccatus was the primary pollination route in the natural environment, which suggested that the breeding system of R. longipedicellatum might be evolving from selfing to outcrossing. The pollination vector of R. longipedicellatum was very specific, in that only B. braccatus was confirmed to deliver pollen to the stigmas. Visitation frequency was influenced by the activity rhythms and resource requirements of the different castes (i.e., sex). B. braccatus workers were the most effective pollinators because of higher visitation frequency and more effective contribution to fruit production, whereas the presence of B. braccatus males might enhance pollen flow within the population to a certain extent. Finally, these findings not only provided a reliable theoretical basis for hybridization breeding of R. longipedicellatum as parents, but also laid a solid foundation for further molecular biology studies to more broadly reveal the mechanisms of its endangerment in the future.

Multiple factors contribute to reproductive isolation between two co-existing Habenaria species (Orchidaceae)

Reproductive isolation is a key feature that forms barriers to gene flow between distinct plants. In orchids, prezygotic reproductive isolation has been considered to be strong, because their associations with highly specific pollinators. In this study, the reproductive ecology and reproductive isolation of two sympatric Habenaria species, H. davidii and H. fordii, was investigated by floral phenology and morphology, hand-pollination experiments and visitor observation in southwest China. The two species were dependent on insects for pollination and completely self-compatible. A number of factors have been identified to limit gene flow between the two species and achieved full reproductive isolation. Ecogeographic isolation was a weak barrier. H. fordii and H. davidii had completely overlapped flowering periods, and floral morphology plays an important role in floral isolation. The two species shared the same hawkmoth pollinator, Cechenena lineosa, but the pollinaria of the two orchids were attached on different body parts of pollinators. Prezygotic isolation was not complete, but the interspecific pollination treatments of each species resulted in no seed sets, indicating that unlike many other orchid species, in which the postzygotic reproductive isolation is very weak or complete absence, the post-zygotic isolation strongly acted in the stage of seed production between two species. The results illustrate the reproductive isolation between two species involves multiple plant life-history stages and a variety of reproductive barriers can contribute to overall isolation.

Natural hybridization and reproductive isolation between two Primula species

Natural hybridization frequently occurs in plants and can facilitate gene flow between species, possibly resulting in species refusion. However, various reproductive barriers block the formation of hybrids and maintain species integrity. Here, we conducted a field survey to examine natural hybridization and reproductive isolation (RI) between sympatric populations of Primula secundiflora and P. poissonii using ten nuclear simple sequence repeat (SSR) loci. Although introgressive hybridization occurred, species boundaries between P. secundiflora and P. poissonii were maintained through nearly complete reproductive isolation. These interfertile species provide an excellent model for studying the RI mechanisms and evolutionary forces that maintain species boundaries.

Asymmetrical natural hybridization varies among hybrid swarms between two diploid Rhododendron species

Background and Aims

The extent to which hybridization leads to gene flow between plant species depends on the structure of hybrid populations. However, if this varies between locations, species barriers might prove permeable in some locations but not in others. To assess possible variation in hybrid population structure, the magnitude and direction of natural hybridization between two Chinese endemic species, Rhododendron spiciferum and Rhododendron spinuliferum , were evaluated.

Methods

Thirteen nuclear microsatellite markers were employed to characterize 566 individuals collected from 15 non-allopatric populations and nine allopatric parental populations. Chloroplast DNA (cpDNA) sequences were obtained from a subset of samples. Genetic structure and direction of gene flow was determined using a combination of STRUCTURE and NEWHYBRIDS analysis.

Key Results

Nuclear analysis revealed that parental taxa formed two genetically distinct clusters and hybrids shared the genetic background of both parents and did not form a separate genetic lineage. Overall, hybrid swarms were dominated by early- and later-generation hybrids, with a significantly higher proportion of hybrids (59·6 %) possessing >50 % R. spiciferum-like nuclear germplasm. The cpDNA analysis further indicated that a significantly greater proportion of hybrids (61·1 %) possessed the R. spiciferum cpDNA haplotype.

Conclusions

Gene flow between R. spiciferum and R. spinuliferum was found to be bidirectional in 14 of the 15 hybrid swarms and asymmetrical in six hybrid swarms. Asymmetrical gene flow was evident for only nuclear DNA (nDNA) in two populations, for only cpDNA in three populations, and for both nDNA and cpDNA in one population. Collectively, the variation in genetic structure found among the 15 hybrid swarms suggests that introgression rather than hybrid speciation is a more likely outcome of hybridization between these hybridizing taxa.