Asymmetric hybridization in Rhododendron agastum: a hybrid taxon comprising mainly F1s in Yunnan, China

Recent advancements in the physiological, genetic, and genomic research on Rhododendrons for trait improvement

High species diversity, hybridization potential, broad geographical dispersal range and ornamental characteristics (i.e., attractive size, shape, structure, flowers, and evergreen) have fetched a good international market for Rhododendron. However, most species are restricted to specific geographic areas due to their habitat specificity in acidic soil and cold climates, resulting many species being classified under threat categories of the IUCN. In this review, advances in research on Rhododendron for improvement to floral display quality and stress resistance have been described. The low genetic barrier among species has created opportunities for extensive hybridization and ploidy alteration for introducing quality and adaptive traits during the development of new varieties. Recent technological advances have supported investigations into the mechanism of flower development, as well as cold tolerance and pathogen resistance mechanisms in the Rhododendron. However, most of the species have limited adaptability to drought, line-tolerance, pathogen resistance, and high-temperature conditions and this resistance ability present in few species largely remains unexplored. Additionally, the available genetic diversity and genomic information on species, and possibilities for their application in molecular breeding have been summarized. Overall, genomic resource data are scarce in the majority of the members of this genus. Finally, various research gaps such as genetic mapping of quality traits, understanding the molecular mechanism of quality-related traits and genomic assortment in Rhododendron members have been discussed in the future perspective section.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04006-6.

Pervasive hybridization during evolutionary radiation of Rhododendron subgenus Hymenanthes in mountains of southwest China

Radiations are especially important for generating species biodiversity in mountainous ecosystems. The contribution of hybridization to such radiations has rarely been examined. Here, we use extensive genomic data to test whether hybridization was involved in evolutionary radiation within Rhododendron subgenus Hymenanthes, whose members show strong geographic isolation in the mountains of southwest China. We sequenced genomes for 143 species of this subgenus and 93 species of four other subgenera, and found that Hymenanthes was monophyletic and radiated during the late Oligocene to middle Miocene. Widespread hybridization events were inferred within and between the identified clades and subclades. This suggests that hybridization occurred both early and late during diversification of subgenus Hymenanthes, although the extent to which hybridization, speciation through mixing-isolation-mixing or hybrid speciation, accelerated the diversification needs further exploration. Cycles of isolation and contact in such and other montane ecosystems may have together promoted species radiation through hybridization between diverging populations and species. Similar radiation processes may apply to other montane floras in this region and elsewhere.

Resolution, conflict and rate shifts: insights from a densely sampled plastome phylogeny for Rhododendron (Ericaceae)

BACKGROUND AND AIMS

Rhododendron is a species-rich and taxonomically challenging genus due to recent adaptive radiation and frequent hybridization. A well-resolved phylogenetic tree would help to understand the diverse history of Rhododendron in the Himalaya-Hengduan Mountains where the genus is most diverse.

METHODS

We reconstructed the phylogeny based on plastid genomes with broad taxon sampling, covering 161 species representing all eight subgenera and all 12 sections, including ~45 % of the Rhododendron species native to the Himalaya-Hengduan Mountains. We compared this phylogeny with nuclear phylogenies to elucidate reticulate evolutionary events and clarify relationships at all levels within the genus. We also estimated the timing and diversification history of Rhododendron, especially the two species-rich subgenera Rhododendron and Hymenanthes that comprise >90 % of Rhododendron species in the Himalaya-Hengduan Mountains.

KEY RESULTS

The full plastid dataset produced a well-resolved and supported phylogeny of Rhododendron. We identified 13 clades that were almost always monophyletic across all published phylogenies. The conflicts between nuclear and plastid phylogenies suggested strongly that reticulation events may have occurred in the deep lineage history of the genus. Within Rhododendron, subgenus Therorhodion diverged first at 56 Mya, then a burst of diversification occurred from 23.8 to 17.6 Mya, generating ten lineages among the component 12 clades of core Rhododendron. Diversification in subgenus Rhododendron accelerated c. 16.6 Mya and then became fairly continuous. Conversely, Hymenanthes diversification was slow at first, then accelerated very rapidly around 5 Mya. In the Himalaya-Hengduan Mountains, subgenus Rhododendron contained one major clade adapted to high altitudes and another to low altitudes, whereas most clades in Hymenanthes contained both low- and high-altitude species, indicating greater ecological plasticity during its diversification.

CONCLUSIONS

The 13 clades proposed here may help to identify specific ancient hybridization events. This study will help to establish a stable and reliable taxonomic framework for Rhododendron, and provides insight into what drove its diversification and ecological adaption. Denser sampling of taxa, examining both organelle and nuclear genomes, is needed to better understand the divergence and diversification history of Rhododendron.

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.

Species divergence with gene flow and hybrid speciation on the Qinghai-Tibet Plateau

The Qinghai-Tibet Plateau (QTP) sensu lato (sl), comprising the platform, the Himalaya and the Hengduan Mountains, is characterized by a large number of endemic plant species. This evolutionary cradle may have arisen from explosive species diversification because of geographic isolation. However, gene flow has been widely detected during the speciation processes of all groups examined, suggesting that natural selection may have also played an important role during species divergence in this region. In addition, natural hybrids have been recovered in almost all species-rich genera. This suggests that numerous species in this region are still 'on the speciation pathway to complete reproductive isolation (RI)'. Such hybrids could directly develop into new species through hybrid polyploidization and homoploid hybrid speciation (HHS). HHS may take place more easily than previously thought through alternate inheritance of alleles of parents at multiple RI loci. Therefore, isolation, selection and hybridization could together have promoted species diversification of numerous plant genera on the QTP sl. We emphasize the need for identification and functional analysis of alleles of major genes for speciation, and especially encourage investigations of parallel adaptive divergence causing RI across different lineages within similar but specific habitats in this region.

Cryptic Species Diversification of the Pedicularis siphonantha Complex (Orobanchaceae) in the Mountains of Southwest China Since the Pliocene

Morphological approaches often fail to delimit species in recently derived species complexes. This can be exacerbated in historical collections which may have lost key features in specimen preparation and preservation. Here, we examine the Pedicularis siphonantha complex, endemic to the Mountains of Southwest China. This complex is characterized by its red/purple/pink and long-tubular corolla, and twisted, beaked galea. However, herbarium specimens are often difficult to identify to species. Molecular approaches using nrITS or nuclear ribosomal internal transcribed spacer (nrITS) + plastid DNA (ptDNA) have been successfully used for species identification in Pedicularis. To resolve taxonomic confusion in the Pedicularis siphonantha complex, we reconstructed phylogenies of the complex using nrITS and four plastid DNA loci (matK, rbcL, trnH-psbA, and trnL-F). To recover as much of the phylogenetic history as possible, we sampled individuals at the population level. Topological incongruence between the nrITS and ptDNA datasets was recovered in clades including two widely distributed species, Pedicularis milliana and Pedicularis tenuituba. Based on morphological, geographical, and genetic evidence, we suggest that hybridization/introgression has occurred between P. milliana and Pedicularis sigmoidea/Pedicularis sp. 1 in the Yulong Snow Mountain of Lijiang, northwest Yunnan, and between P. tenuituba and Pedicularis leptosiphon in Ninglang, northwest Yunnan. After removing conflicting DNA regions in Pedicularis dolichosiphon (nrITS) and P. milliana (ptDNA), the concatenated nrITS and ptDNA phylogenies distinguish 11 species in the P. siphonantha complex, including two undescribed species, from the Jiaozi and Yulong Snow Mountains, respectively. Phylogeographical analyses indicate that the P. siponantha complex originated from south of the Hengduan Mountains, expanding north to the Himalayas and the Yunnan-Guizhou Plateau. Moreover, the uplift of the Qinghai-Tibet Plateau and climate oscillations may have driven further diversification in the complex.

Plastome phylogenomics of Allaeanthus, Broussonetia and Malaisia (Dorstenieae, Moraceae) and the origin of B. × kazinoki

Species of Broussonetia have been essential in the development of papermaking technology. In Japan and Korea, a hybrid between B. monoica and B. papyrifera (= B. × kazinoki) known as kōzo and daknamu is still the major source of raw materials for making traditional paper washi and hanji, respectively. Despite their cultural and practical significance, however, the origin and taxonomy of kōzo and daknamu remain controversial. Additionally, the long-held generic concept of Broussonetia s.l., which included Sect. Allaeanthus and Sect. Broussonetia, was challenged as phylogenetic analyses showed Malaisia is sister to the latter section. To re-examine the taxonomic proposition that recognizes Allaeanthus, Broussonetia, and Malaisia (i.e., Broussonetia alliance), plastome and nuclear ribosomal DNA (nrDNA) sequences of six species of the alliance were assembled. Characterized by the canonical quadripartite structure, genome alignments and contents of the six plastomes (160,121-162,594 bp) are highly conserved, except for the pseudogenization and/or loss of the rpl22 gene. Relationships of the Broussonetia alliance are identical between plastome and nrDNA trees, supporting the maintenance of Malaisia and the resurrection of Allaeanthus. The phylogenomic relationships also indicate that the monoecy in B. monoica is a derived state, possibly resulting from hybridization between the dioecious B. kaempferi (♀) and B. papyrifera (♂). Based on the hypervariable ndhF-rpl32 intergenic spacer selected by sliding window analysis, phylogeographic analysis indicates that B. monoica is the sole maternal parent of B. × kazinoki and that daknamu carries multiple haplotypes, while only one haplotype was detected in kōzo. Because hybridizations between B. monoica and B. papyrifera are unidirectional and have occurred rarely in nature, our data suggest that daknamu might have originated via deliberate hybrid breeding selected for making hanji in Korea. On the contrary, kōzo appears to have a single origin and the possibility of a Korean origin cannot be ruled out.

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.

Natural Hybrid Origin of the Controversial "Species" Clematis × pinnata (Ranunculaceae) Based on Multidisciplinary Evidence

Interspecific hybridization is common and has often been viewed as a driving force of plant diversity. However, it raises taxonomic problems and thus impacts biodiversity estimation and biological conservation. Although previous molecular phylogenetic studies suggested that interspecific hybridization may be rather common in Clematis, and artificial hybridization has been widely applied to produce new Clematis cultivars for nearly two centuries, the issue of natural hybridization of Clematis has never been addressed in detail. In this study, we tested the hybrid origin of a mesophytic and cold-adapted vine species, Clematis pinnata, which is a rare and taxonomically controversial taxon endemic to northern China. Using field investigations, flow cytometry (FCM), phylogenomic analysis, morphological statistics, and niche modeling, we tested hybrid origin and species status of C. pinnata. The FCM results showed that all the tested species were homoploid (2n = 16). Phylonet and HyDe analyses based on transcriptome data showed the hybrid origins of C. × pinnata from either C. brevicaudata × C. heracleifolia or C. brevicaudata × C. tubulosa. The plastome phylogeny depicted that C. × pinnata in different sampling sites originated by different hybridization events. Morphological analysis showed intermediacy of C. × pinnata between its putative parental species in many qualitative and quantitative characters. Niche modeling results suggested that C. × pinnata had not been adapted to a novel ecological niche independent of its putative parents. These findings demonstrated that plants of C. × pinnata did not formed a self-evolved clade and should not be treated as a species. The present study also suggests that interspecific hybridization is a common mechanism in Clematis to generate diversity and variation, and it may play an important role in the evolution and diversification of this genus. Our study implies that morphological diversity caused by natural hybridization may overstate the real species diversity in Clematis.

Incongruent phylogenies and their implications for the study of diversification, taxonomy, and genome size evolution of Rhododendron

PREMISE

Classification of taxa depends on the quality of inferred phylogenies. Rhododendron, a highly species-rich genus (>1156 species) of woody plants, has a highly debated infrageneric classification, due to its huge diversity, homoplasy in key characters, and incongruence among data sets. We provide a broad coverage of representative species to resolve Rhododendron infrageneric phylogeny and highlight the areas of incongruence. We further investigate the effect of polyploidy and genome size evolution on diversification of Rhododendron.

METHODS

We generated two plastid and two nuclear loci for 260 Rhododendron species. We analyzed the loci separately as well as concatenated, utilizing both likelihood and Bayesian methods. We tested incongruence both among the data sets and with previous studies. We estimated genome sizes for 125 species through flow cytometry.

RESULTS

Our results suggest stronger support for larger subgenera; however, the smaller subgenera pose several problems; for example, R. tomentosum (former genus Ledum) occupies incongruent positions based on different DNA regions. The main shift to higher diversification in the genus occurs in the Himalayan/Southeast Asian clade of R. subg. Hymenanthes. We found that polyploidy occurs in almost all subgenera but most frequently within R. subg. Rhododendron sections Rhododendron and Schistanthe.

CONCLUSIONS

We endorse the recognition of five major clades at the subgeneric level, but a number of species cannot be confidently assigned to these clades due to incongruency. With regard to genome size evolution, results support previous reports that genome sizes of tropical plants are lower than those of colder and temperate regions and that genome downsizing promotes diversification.

Testing genome skimming for species discrimination in the large and taxonomically difficult genus Rhododendron

Standard plant DNA barcodes based on 2-3 plastid regions, and nrDNA ITS show variable levels of resolution, and fail to discriminate among species in many plant groups. Genome skimming to recover complete plastid genome sequences and nrDNA arrays has been proposed as a solution to address these resolution limitations. However, few studies have empirically tested what gains are achieved in practice. Of particular interest is whether adding substantially more plastid and nrDNA characters will lead to an increase in discriminatory power, or whether the resolution limitations of standard plant barcodes are fundamentally due to plastid genomes and nrDNA not tracking species boundaries. To address this, we used genome skimming to recover near-complete plastid genomes and nuclear ribosomal DNA from Rhododendron species and compared discrimination success with standard plant barcodes. We sampled 218 individuals representing 145 species of this species-rich and taxonomically difficult genus, focusing on the global biodiversity hotspots of the Himalaya-Hengduan Mountains. Only 33% of species were distinguished using ITS+matK+rbcL+trnH-psbA. In contrast, 55% of species were distinguished using plastid genome and nrDNA sequences. The vast majority of this increase is due to the additional plastid characters. Thus, despite previous studies showing an asymptote in discrimination success beyond 3-4 plastid regions, these results show that a demonstrable increase in discriminatory power is possible with extensive plastid genome data. However, despite these gains, many species remain unresolved, and these results also reinforce the need to access multiple unlinked nuclear loci to obtain transformative gains in species discrimination in plants.

Species delimitation and hybridization history of a hazel species complex

BACKGROUND AND AIMS

Hybridization increases species adaptation and biodiversity but also obscures species boundaries. In this study, species delimitation and hybridization history were examined within one Chinese hazel species complex (Corylus chinensis-Corylus fargesii). Two species including four varieties have already been described for this complex, with overlapping distributions.

METHODS

A total of 322 trees from 44 populations of these four varieties across their ranges were sampled for morphological and molecular analyses. Climatic datasets based on 108 geographical locations were used to evaluate their niche differentiations. Flowering phenology was also observed for two co-occurring species or varieties.

KEY RESULTS

Four statistically different phenotypic clusters were revealed, but these clusters were highly inconsistent with the traditional taxonomic groups. All the clusters showed statistically distinct niches, with complete or partial geographical isolation. Only two clusters displayed a distributional overlap, but they had distinct flowering phenologies at the site where they co-occurred. Population-level evidence based on the genotypes of ten simple sequence repeat loci supported four phenotypic clusters. In addition, one cluster was shown to have an admixed genetic composition derived from the other three clusters through repeated historical hybridizations.

CONCLUSIONS

Based on our new evidence, it is better to treat the four clusters identified here as four independent species. One of them was shown to have an admixed genetic composition derived from the other three through repeated historical hybridizations. This study highlights the importance of applying integrative and statistical methods to infer species delimitations and hybridization history. Such a protocol should be adopted widely for future taxonomic studies.

Untargeted metabolite profiling of petal blight in field-grown Rhododendron agastum using GC-TOF-MS and UHPLC-QTOF-MS/MS

Petal blight caused by fungi is among the most destructive diseases of Rhododendron, especially Rhododendron agastum. Nonetheless, the metabolite changes that occur during petal blight are unknown. We used untargeted gas chromatography time-of-flight mass spectrometry (GC-TOF-MS) and ultra-high performance liquid chromatography quadrupole time-of-flight tandem mass spectrometry (UHPLC-QTOF-MS/MS) to compare the metabolite profiles of healthy and petal blight R. agastum flowers. Using GC-TOF-MS, 571 peaks were extracted, of which 189 metabolites were tentatively identified. On the other hand, 364 and 277 metabolites were tentatively identified in the positive and negative ionization modes of the UHPLC-QTOF-MS/MS, respectively. Principal component analysis (PCA) and orthogonal projections to latent structures-discriminant analysis (OPLS-DA) were able to clearly discriminate between healthy and petal blight flowers. Differentially abundant metabolites were primarily enriched in the biosynthesis of specialized metabolites. 17 accumulated specialized metabolites in petal blight flowers have been reported to have antifungal activity, and literature indicates that 9 of them are unique to plants. 3 metabolites (chlorogenic acid, medicarpin, and apigenin) are reportedly involved in resistance to blight caused by pathogens. We therefore speculate that the accumulation of chlorogenic acid, medicarpin, and apigenin may be involved in the resistance to petal blight. Our results suggest that these metabolites may be used as candidate biocontrol agents for the control fungal petal blight in Rhododendron.

Natural hybridization between two butterfly bushes in Tibet: dominance of F1 hybrids promotes strong reproductive isolation

BACKGROUND

It has been recognized that a certain amount of habitat disturbance is a facilitating factor for the occurrence of natural hybridization, yet to date we are unaware of any studies exploring hybridization and reproductive barriers in those plants preferentially occupying disturbed habitats. Buddleja plants (also called butterfly bush) generally do grow in disturbed habitats, and several species with hybrid origin have been proposed, based solely on morphological evidence.

RESULTS

In the present study, we test the hypothesis that B. × wardii is of natural hybridization origin in two sympatric populations of three taxa including B. × wardii and its parents (B. alternifolia and B. crispa) plus 4 referenced parental populations, using four nuclear genes and three chloroplast intergenic spacers, as well as with 10 morphological characters. Our results suggest that at both sites B. × wardii is likely to be a hybrid between B. alternifolia and B. crispa, and moreover, we confirm that most of the hybrids examined are F1s. That these plants are F1s is further supported by morphology, as no transgressive characters were detected. B. crispa was found to be the maternal parent in the Bahe (BH) population, from cpDNA evidence. However, in the Taji (TJ) population, the direction of hybridization was difficult to establish due to the shared cpDNA haplotypes between B. alternifolia and B. crispa, however we still predicted a similar unidirectional hybridization pattern due to results from cross-specific pollination treatments which supported the "SI × SC rule".

CONCLUSIONS

The presence of mainly F1 hybrids can successfully impede gene flow and thus maintain species boundaries in parental species in a typical distribution of Buddleja, i.e. in disturbed habitats.

Species Identification in the Rhododendron vernicosum-R. decorum Species Complex (Ericaceae)

Delimitating species boundaries is the primary aim of biological classification and could be critical for evaluating the evolving process of species and conserving biodiversity. Rhododendron is an iconic group with an extraordinary diversity in southwest China. However, it remains unknown whether the recorded species therein comprise independently evolving lineages or artificially delimitated morphological entities. In this study, we carried out species delimitation of four Rhododendron species in the R. vernicosum-R. decorum species complex based on morphological analyses and population genetic data from nuclear simple sequence repeats (SSR) markers. We randomly selected a total of 105 specimens of different individuals identified as four species across their distributional ranges to examine the statistically distinct phenotypic clusters based on 19 morphological traits. Similarly, we genotyped 55 individuals of four species from 21 populations using 15 SSR markers. The morphological analyses sorted R. decorum and the other three species into two different phenotypic clusters. The genetic clusters were consistent with the morphological clusters. However, we also recovered the third genetic cluster, comprising six R. vernicosum populations and containing the admixed genetic compositions of the other two distinct genetic clusters. This hybrid group was morphologically similar to the typical R. vernicosum (including the samples from its type specimen locality and both R. verruciferum and R. gonggashanense) but with more genetic ancestry from R. decorum. Based on our findings, we identify two distinct species and one putative hybrid group due to introgression in the R. vernicosum-R. decorum species complex. We propose to merge R. verruciferum and R. gonggashanense into R. vernicosum based on genetic compositions and our morphological analyses. The hybrid group inferred from our findings, however, needs further investigations.

Rhamnella intermedia (Rhamnaceae), a new evergreen species from southwest Guangxi

Rhamnella intermedia, a new evergreen species from southwest Guangxi, is described and illustrated in this study. This species is similar to R. brachycarpa by the size and ratio of length to width of dried fruit and seeds, by which it differs from R. rubrinervis and R. tonkinensis. However, it differs from R. brachycarpa by rarely mucronate seed apices, larger ratio of length to width of leaves, leaf apices acuminate to long acuminate, shorter leaf petioles, and longer fruiting pedicels. Principal component analysis based on phenotypic traits further recognised three separated groups. Rhamnella rubrinervis and R. tonkinensis were clustered into one group; the other two groups represented R. brachycarpa and two Guangxi populations, respectively. Furthermore, phylogenetic analysis of nuclear ITS sequence variations highly supported that the two Guangxi populations represented an independent evolutionary lineage and were closest to R. rubrinervis. Four fixed nucleotide sites were found and were different from R. rubrinervis. However, besides the differentiated traits in seeds and fruit, densely pilose young branches also separated them from R. rubrinervis. In addition, during our field investigations, none of the three closely related species were found at locations where this new species was distributed. Therefore, this new species, based on the two Guangxi populations, is named R. intermedia. The key to four closely related species is also presented.

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.

Characterization of Volatile Compounds in Four Different Rhododendron Flowers by GC×GC-QTOFMS

Volatile compounds in flowers of Rhododendron delavayi, R. agastum, R. annae, and R. irroratum were analyzed using comprehensive two-dimensional gas chromatography-mass spectrometry (GC×GC) coupled with high-resolution quadrupole time-of-flight mass spectrometry (QTOFMS). A significantly increased number of compounds was separated by GC×GC compared to conventional one-dimensional GC (1DGC), allowing more comprehensive understanding of the volatile composition of Rhododendron flowers. In total, 129 volatile compounds were detected and quantified. Among them, hexanal, limonene, benzeneacetaldehyde, 2-nonen-1-ol, phenylethyl alcohol, citronellal, isopulegol, 3,5-dimethoxytoluene, and pyridine are the main compounds with different content levels in all flower samples. 1,2,3-trimethoxy-5-methyl-benzene exhibits significantly higher content in R. irroratum compared to in the other three species, while isopulegol is only found in R. irroratum and R. agastum.

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.

Incomplete reproductive isolation between Rhododendron taxa enables hybrid formation and persistence

The evolutionary consequences of hybridization ultimately depend on the magnitude of reproductive isolation between hybrids and their parents. We evaluated the relative contributions of pre- and post-zygotic barriers to reproduction for hybrid formation, hybrid persistence and potential for reproductive isolation of hybrids formed between two Rhododendron species, R. spiciferum and R. spinuliferum. Our study established that incomplete reproductive isolation promotes hybrid formation and persistence and delays hybrid speciation. All pre-zygotic barriers to reproduction leading to hybrid formation are incomplete: parental species have overlapping flowering; they share the same pollinators; reciprocal assessments of pollen tube germination and growth do not differ among parents. The absence of post-zygotic barriers between parental taxa indicates that the persistence of hybrids is likely. Reproductive isolation was incomplete between hybrids and parents in all cases studied, although asymmetric differences in reproductive fitness were prevalent and possibly explain the genetic structure of natural hybrid swarms where hybridization is known to be bidirectional but asymmetric. Introgression, rather than speciation, is a probable evolutionary outcome of hybridization between the two Rhododendron taxa. Our study provides insights into understanding the evolutionary implications of natural hybridization in woody plants.

Maintenance of species boundaries in three sympatric Ligularia (Senecioneae, Asteraceae) species

The key process in speciation concerns the formation and maintenance of reproductive isolating barriers between diverging lineages. Although species boundaries are frequently investigated between two species across many taxa, reproductive isolating barriers among multiple species (>2) that would represent the most common phenomenon in nature, remain to be clarified. Here, we use double digest restriction-site associated DNA (ddRAD) sequencing to examine patterns of hybridization at a sympatric site where three Ligularia species grow together and verify whether those patterns contribute to the maintenance of boundaries among species. The results based on the RAD SNP datasets indicated hybridization Ligularia cyathiceps × L. duciformis and L. duciformis × L. yunnanensis were both restricted to F₁ s plus a few first-generation backcrosses and no gene introgression were identified, giving rise to strong reproductive isolation among hybridizing species. Moreover, hybrid swarm simulation, using HYBRIDLAB, indicated the RAD SNP datasets had sufficient discriminatory power for accurate hybrid detection. We conclude that parental species show strong reproductive isolation and they still maintain species boundaries, which may be the key mechanism to maintain species diversity of Ligularia in the eastern Qinghai-Tibetan Plateau and adjacent areas. Moreover, this study highlights the effectiveness of RAD sequencing in hybridization studies.

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.

Evidence for continual hybridization rather than hybrid speciation between Ligularia duciformis and L. paradoxa (Asteraceae)

Background

Hybrids possess phenotypic traits that are often intermediate between their parental taxa, which commonly serves as evidence of hybridization in morphological analyses. Natural hybridization has been shown to occur frequently in Ligularia (Asteraceae). In a previous study, Ligularia ×maoniushanensis was demonstrated as a natural hybrid species between L. duciformis and L. paradoxa based on morphological and reproductive traits.

Methods

We used three chloroplast (cpDNA) fragments (psbA-trnH, trnL-rpl32 and trnQ-5′rps16), the nuclear ribosomal internal transcribed spacer (nrITS), and co-dominant SSR and dominant ISSR markers to study natural hybridization between L. duciformis and L. paradoxa growing sympatrically in two locations. Parental taxa were inferred using network analyses of cpDNA and nrITS haplotypes. Admixture among individuals was examined using the Bayesian clustering programs STRUCTURE and NewHybrids based on the SSR and ISSR data; and potential introgression in the SSR loci was assessed using the INTROGRESS package.

Results

The putative parental species were clearly distinguished from other sympatric Ligularia species by nrITS data, and L. ×maoniushanensis individuals were confirmed to be the hybrid offspring of L. duciformis and L. paradoxa. Moreover, introgression was detected among several individuals morphologically identified as L. duciformis or L. paradoxa. Analyses of the cpDNA data revealed primarily unidirectional hybridization between L. duciformis and L. paradoxa, with L. paradoxa as the maternal parent in Mt. Maoniu, whereas bidirectional but asymmetrical hybridization was inferred to occur in Heihai Lake. The STRUCTURE analyses based on the SSR data detected two distinct clusters among the three taxa. The NewHybrids analyses showed that individuals circumscribed as L. ×maoniushanensis were dominated by early- and later-generation and backcrossing hybrids. The NewHybrids results based on the ISSR data were congruent with SSR results. In addition, introgression was detected in some SSR loci, and heterogeneity among loci was found in terms of detected patterns of introgression.

Conclusions

Our data provide strong evidence for hybridization and introgression between L. duciformis and L. paradoxa. Ligularia ×maoniushanensis was demonstrated to be of hybrid origin. Since no evident reproductive isolation was found between the two parental species, detected hybrids appear to be part of hybrid swarms resulting from frequent and ongoing gene flow, which might impede the formation of a new hybrid species.

The draft genome assembly of Rhododendron delavayi Franch. var. delavayi

Rhododendron delavayi Franch. is globally famous as an ornamental plant. Its distribution in southwest China covers several different habitats and environments. However, not much research had been conducted on Rhododendron spp. at the molecular level, which hinders understanding of its evolution, speciation, and synthesis of secondary metabolites, as well as its wide adaptability to different environments. Here, we report the genome assembly and gene annotation of R. delavayi var. delavayi (the second genome sequenced in the Ericaceae), which will facilitate the study of the family. The genome assembly will have further applications in genome-assisted cultivar breeding. The final size of the assembled R. delavayi var. delavayi genome (695.09 Mb) was close to the 697.94 Mb, estimated by k-mer analysis. A total of 336.83 gigabases (Gb) of raw Illumina HiSeq 2000 reads were generated from 9 libraries (with insert sizes ranging from 170 bp to 40 kb), achieving a raw sequencing depth of ×482.6. After quality filtering, 246.06 Gb of clean reads were obtained, giving ×352.55 coverage depth. Assembly using Platanus gave a total scaffold length of 695.09 Mb, with a contig N50 of 61.8 kb and a scaffold N50 of 637.83 kb. Gene prediction resulted in the annotation of 32 938 protein-coding genes. The genome completeness was evaluated by CEGMA and BUSCO and reached 95.97% and 92.8%, respectively. The gene annotation completeness was also evaluated by CEGMA and BUSCO and reached 97.01% and 87.4%, respectively. Genome annotation revealed that 51.77% of the R. delavayi genome is composed of transposable elements, and 37.48% of long terminal repeat elements (LTRs). The de novo assembled genome of R. delavayi var. delavayi (hereinafter referred to as R. delavayi) is the second genomic resource of the family Ericaceae and will provide a valuable resource for research on future comparative genomic studies in Rhododendron species. The availability of the R. delavayi genome sequence will hopefully provide a tool for scientists to tackle open questions regarding molecular mechanisms underlying environmental interactions in the genus Rhododendron, more accurately understand the evolutionary processes and systematics of the genus, facilitate the identification of genes encoding pharmaceutically important compounds, and accelerate molecular breeding to release elite varieties.

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.

Population Genetic Differentiation and Taxonomy of Three Closely Related Species of Saxifraga (Saxifragaceae) from Southern Tibet and the Hengduan Mountains

The effects of rapid, recent uplift of the Hengduan Mountains on evolution and diversification of young floristic lineages still remain unclear. Here, we investigate diversification of three closely related Saxifraga species with a distribution restricted to the Hengduan Mountains (HM) and southern Tibet, and comment on their taxonomy based on molecular evidence. Three chloroplast DNA fragments (rbcL, trnL-F, trnS-G) and the nuclear ribosomal DNA internal transcribed spacer (ITS) were employed to study genetic structure across 104 individuals from 12 populations of Saxifraga umbellulata, S. pasumensis, and S. banmaensis. Chloroplast DNA (cpDNA) phylogenies revealed two well supported clades, corresponding to S. umbellulata and S. pasumensis plus S. banmaensis. Topology of the ITS phylogeny was largely congruent with that generated from cpDNA haplotypes, but with minor conflicts which might be caused by incomplete lineage sorting. Analyses of molecular variance of both cpDNA and ITS datasets revealed that most variation was held between S. pasumensis s.l. (with S. banmaensis) and S. umbellulata (92.31% for cpDNA; 69.78% for ITS), suggesting a high degree of genetic divergence between them. Molecular clock analysis based on ITS dataset suggested that the divergence between S. pasumensis s.l. and S. umbellulata can be dated to 8.50 Ma, probably a result of vicariant allopatric diversification associated with the uplift events of the HM. Vicariance associated with HM uplifts may also have been responsible for infraspecific differentiation in S. pasumensis. In contrast, infraspecific differentiation in S. umbellulata was most likely triggered by Quaternary glaciations. The much lower levels of gene diversity within populations of S. pasumensis compared with S. umbellulata could have resulted from both range contractions and human collection on account of its putative medicinal properties. Combining evidence from morphology, geographical distributions and molecular phylogenetic data, we recommend that S. banmaensis should be treated as a synonym of S. pasumensis which in turn, and based on the same sources of evidence, should be treated as a separate species rather than as a variety of S. umbellulata.

Why is population information crucial for taxonomy? A case study involving a hybrid swarm and related varieties

Hybridization has become a focal topic in evolutionary biology, and many taxonomists are aware that the process occurs more frequently than previously assumed. Nonetheless many species and varieties are still described without explicitly considering the possibility of hybridization, especially in countries that have relatively short scientific histories, but which often possess the highest species diversities. Furthermore, new taxa are often described based only on herbarium specimens, not taking into account information from wild populations, significantly decreasing the potential to detect morphologies arising from hybridization at this crucial descriptive stage. We used morphological data from a hybrid swarm involving two Rhododendron species to showcase possible character combinations in intermediates. Certain characters used to distinguish taxa were more variable within the same individual than between species, emphasizing the importance of population information for an adequate choice of characters. Most described varieties of the two species fell within the spectrum of hybrid morphology, suggesting that these taxa would be unlikely to have merited formal description if contemporary standards had been employed. In all investigated cases the hybrid nature of described varieties seems to have been detectable with adequate morphological data alone, if populations had been assessed. A post hoc assessment of taxa is often complicated, especially if certain types of information are not provided. To avoid accumulation of such invalid taxa, careful scrutiny should be employed for new descriptions. Hybrids (not hybrid species) described as taxa obscure valuable information about natural processes and impact negatively on further research that depends on taxonomic data.

Fitness of natural willow hybrids in a pioneer mosaic hybrid zone

Hybrid fitness is an important parameter to predict the evolutionary consequences of a hybridization event and to characterize hybrid zones. We studied fitness parameters of F1 and later-generation hybrids between the lowland species Salix purpurea and the alpine S. helvetica that have recently emerged during colonization of an alpine glacier forefield. Fruit production (number of capsules per catkin and fruit set) did not differ between hybrids and parents, but the number of seeds per capsule of F1 hybrids was slightly lower than that of later-generation hybrids and of the parents. Germination rates and seedling growth were tested on three substrates (pH 4.5, 7.0, and 8.0). Germination rates of seeds collected from F1 hybrids were lower on acid and neutral substrates, but equal at pH 8.0 compared to all other groups, while the seeds from later-generation hybrids performed as well as the parents on all three substrates. In seedling growth, the colonizer S. purpurea performed better than all other taxa on all three substrates, while hybrids resembled the subalpine species S. helvetica. Results suggest that endogenous selection acts against F1 hybrids, but favors fitter genotypes in later-generation hybrids. Exogenous selection via soil pH appears to be weak during seedling establishment. The pioneer vegetation on the glacier forefield may offer sufficient niche space for hybrid seedlings. Owing to the relatively high fitness of the hybrids and the scattered distribution of hybrids and parental individuals on the glacier forefield, this hybrid zone can be assigned to a mosaic model, probably facilitating gene flow and introgression between the parental species. As establishment of the hybrid zone appears to be linked to a colonization process, we propose to call it a pioneer mosaic hybrid zone.

Identification and cloning of class II and III chitinases from alkaline floral nectar of Rhododendron irroratum, Ericaceae

Class II and III chitinases belonging to different glycoside hydrolase families were major nectarins in Rhododendron irroratum floral nectar which showed significant chitinolytic activity. Previous studies have demonstrated antimicrobial activity in plant floral nectar, but the molecular basis for the mechanism is still poorly understood. Two chitinases, class II (Rhchi2) and III (Rhchi3), were characterized from alkaline Rhododendron irroratum nectar by both SDS-PAGE and mass spectrometry. Rhchi2 (27 kDa) and Rhchi3 (29 kDa) are glycoside hydrolases (family 19 and 18) with theoretical pI of 8.19 and 7.04. The expression patterns of Rhchi2 and Rhchi3 were analyzed by semi-quantitative RT-PCR. Rhchi2 is expressed in flowers (corolla nectar pouches) and leaves while Rhchi3 is expressed in flowers. Chitinase in concentrated protein and fresh nectar samples was visualised by SDS-PAGE and chitinolytic activity in fresh nectar was determined spectrophotometrically via chitin-azure. Full length gene sequences were cloned with Tail-PCR and RACE. The amino acid sequence deduced from the coding region for these proteins showed high identity with known chitinases and predicted to be located in extracellular space. Fresh R. irroratum floral nectar showed significant chitinolytic activity. Our results demonstrate that class III chitinase (GH 18 family) also exists in floral nectar. The functional relationship between class II and III chitinases and the role of these pathogenesis-related proteins in antimicrobial activity in nectar is suggested.

Molecular evidence of hybrid zones of Cedrela (Meliaceae) in the Yungas of Northwestern Argentina

In the Yungas of Northwestern Argentina, three endangered species of Cedrela (C. angustifolia, C. saltensis, and C. balansae) follow altitudinal gradients of distribution with contact zones between them. We sampled 210 individuals from 20 populations that spanned most of Cedrela's geographical range in the Yungas, and used Amplified Fragment Length Polymorphism (AFLP) markers and DNA sequences of the nuclear Internal Transcribed Spacer (ITS) to investigate hybrid zones. Data analyses employed an array of complementary methods, including principal coordinate analyses, Bayesian clustering analyses, maximum likelihood tree-building, and network techniques. Both nuclear molecular systems - AFLP and ITS - provided insights into the evolutionary history of Cedrela in the Yungas in a congruent manner. We uncovered strong support for the occurrence of natural hybridization between C. balansae and C. saltensis. Additionally, we identified hybrid zones in areas of sympatry (at both the Calilegua National Park and the San Andrés farm) and in transition zones from 820 to 1100meters above sea level (localities of Pintascayo and Acambuco). There was no evidence for hybridization of either C. balansae or C. saltensis with C. angustifolia. The role of hybrid populations in conservation and use of genetic resources in the Yungas were discussed.

Genetic origin and composition of a natural hybrid poplar Populus × jrtyschensis from two distantly related species

BACKGROUND

The factors that contribute to and maintain hybrid zones between distinct species are highly variable, depending on hybrid origins, frequencies and fitness. In this study, we aimed to examine genetic origins, compositions and possible maintenance of Populus × jrtyschensis, an assumed natural hybrid between two distantly related species. This hybrid poplar occurs mainly on the floodplains along the river valleys between the overlapping distributions of the two putative parents.

RESULTS

We collected 566 individuals from 45 typical populations of P. × jrtyschensis, P. nigra and P. laurifolia. We genotyped them based on the sequence variations of one maternally inherited chloroplast DNA (cpDNA) fragment and genetic polymorphisms at 20 SSR loci. We further sequenced eight nuclear genes for 168 individuals from 31 populations. Two groups of cpDNA haplotypes characteristic of P. nigra and P. laurifolia respectively were both recovered for P. × jrtyschensis. Genetic structures and coalescent tests of two sets of nuclear population genetic data suggested that P. × jrtyschensis originated from hybridizations between the two assumed parental species. All examined populations of P. × jrtyschensis comprise mainly F1 hybrids from interspecific hybridizations between P. nigra and P. laurifolia. In the habitats of P. × jrtyschensis, there are lower concentrations of soil nitrogen than in the habitats occupied by the other two species.

CONCLUSIONS

Our extensive examination of the genetic composition of P. × jrtyschensis suggested that it is typical of F1-dominated hybrid zones. This finding plus the low concentration of soil nitrogen in the floodplain soils support the F1-dominated bounded hybrid superiority hypothesis of hybrid zone maintenance for this particular hybrid poplar.

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

Reproductive isolation (RI) plays an important role for speciation, but assessing reproductive barriers at all life-cycle stages remains challenging. In plants, most studies addressing the topic have been focusing on herbs with short generation times. The present study attempted to quantify several reproductive barriers between a hybridizing species pair of long-lived woody rhododendrons. Consistent with findings of previous studies, pre-zygotic reproductive barriers contributed more to total RI than post-zygotic reproductive barriers. Especially in the more widespread species geographic isolation was an important barrier, and pollinator constancy contributed exceptionally to RI in both species. Additionally to strong pre-zygotic reproductive barriers, post-zygotic reproductive barriers were considerable, and had asymmetric tendencies favoring one of the species as maternal parent. Overall, despite occasional hybridization, the present study provides evidence for strong RI between R. cyanocarpum and R. delavayi.

Evidence for asymmetrical hybridization despite pre- and post-pollination reproductive barriers between two Silene species

Interspecific hybridization is widespread among plants; nevertheless, pre- and post-zygotic isolating mechanisms may maintain species integrity for interfertile species in sympatry despite some gene flow. Interspecific hybridization and potential isolating barriers were evaluated between co-flowering Silene asclepiadea and Silene yunnanensis in an alpine community in southwest China. We investigated morphological and molecular (nuclear microsatellites and chloroplast gene sequence) variation in sympatric populations of S. asclepiadea and S. yunnanensis. Additionally, we analyzed pollinator behaviour and compared reproductive success between the putative hybrids and their parental species. Both the molecular and morphological data indicate that there were putative natural hybrids in the field, with S. asclepiadae the ovule parent and S. yunnanensis the pollen parent. Bumblebees were the primary visitors to S. asclepiadae and putative hybrids, while butterflies were the primary visitors to S. yunnanensis Pollen production and viability were significantly lower in putative hybrids than the parental species. The direction of hybridization is quite asymmetric from S. yunnanensis to S. asclepiadea Protandry combined with later peak flowering of S. yunnanensis, and pollinator preference may have contributed to the asymmetric pattern of hybridization, but putative hybrids were rare. Our results thus suggest that despite gene flow, S. asclepiadea and S. yunnanensis can maintain species boundaries, perhaps as a result of floral isolation and low fecundity of the hybrids.

Hybridization in closely related Rhododendron species: half of all species-differentiating markers experience serious transmission ratio distortion

An increasing number of studies of hybridization in recent years have revealed that complete reproductive isolation between species is frequently not finalized in more or less closely related organisms. Most of these species do, however, seem to retain their phenotypical characteristics despite the implication of gene flow, highlighting the remaining gap in our knowledge of how much of an organism’s genome is permeable to gene flow, and which factors promote or prevent hybridization. We used AFLP markers to investigate the genetic composition of three populations involving two interfertile Rhododendron species: two sympatric populations, of which only one contained hybrids, and a further hybrid-dominated population. No fixed differences between the species were found, and only 5.8% of the markers showed some degree of species differentiation. Additionally, 45.5% of highly species-differentiating markers experienced significant transmission distortion in the hybrids, which was most pronounced in F1 hybrids, suggesting that factors conveying incompatibilities are still segregating within the species. Furthermore, the two hybrid populations showed stark contrasting composition of hybrids; one was an asymmetrically backcrossing hybrid swarm, while in the other, backcrosses were absent, thus preventing gene flow.

Species delimitation in plants using the Qinghai-Tibet Plateau endemic Orinus (Poaceae: Tridentinae) as an example

BACKGROUND AND AIMS

Accurate identification of species is essential for the majority of biological studies. However, defining species objectively and consistently remains a challenge, especially for plants distributed in remote regions where there is often a lack of sufficient previous specimens. In this study, multiple approaches and lines of evidence were used to determine species boundaries for plants occurring in the Qinghai-Tibet Plateau, using the genus Orinus (Poaceae) as a model system for an integrative approach to delimiting species.

METHODS

A total of 786 individuals from 102 populations of six previously recognized species were collected for niche, morphological and genetic analyses. Three plastid DNA regions (matK, rbcL and trnH-psbA) and one nuclear DNA region [internal transcribed space (ITS)] were sequenced.

KEY RESULTS

Whereas six species had been previously recognized, statistical analyses based on character variation, molecular data and niche differentiation identified only two well-delimited clusters, together with a third possibly originating from relatively recent hybridization between, or historical introgression from, the other two.

CONCLUSIONS

Based on a principle of integrative species delimitation to reconcile different sources of data, the results provide compelling evidence that the six previously recognized species of the genus Orinus that were examined should be reduced to two, with new circumscriptions, and a third, identified in this study, should be described as a new species. This empirical study highlights the value of applying genetic differentiation, morphometric statistics and ecological niche modelling in an integrative approach to re-circumscribing species boundaries. The results produce relatively objective, operational and unbiased taxonomic classifications of plants occurring in remote regions.

Variation in the frequency and extent of hybridization between Leucosceptrum japonicum and L. stellipilum (Lamiaceae) in the Central Japanese Mainland

Variations in the frequency and extent of hybridization among mixed populations located in the same contact zone provide natural laboratories for the study of extrinsic reproductive isolation maintaining species integrity. In this study, we examined the pattern of hybridization between L. japonicum and L. stellipilum among mixed populations in different localities of a contact zone. The genetic structures from three sympatric populations and six mixed populations in the hybrid zone, and five reference populations far from the contact zone, were characterized using 10 neutral nuclear microsatellite markers. Evidence from genetic distance-based clustering analysis, the frequency distribution of admixture proportion values, and the hybrid category assignment approaches indicated that the frequency and extent of hybridization varied considerably among populations in the contact zone between L. japonicum and L. stellipilum. One likely explanation is that variation in exogenous (ecological) selection among populations might contribute to differences in frequency and extent of hybridization. The present study will facilitate future research exploring the evolution of reproductive isolation between L. japonicum and L. stellipilum.

AFLPsim: an R package to simulate and detect dominant markers under selection in hybridizing populations

BACKGROUND

In spite of a large diversity of approaches to investigate loci under selection from a population genetic perspective, very few programs have been specifically designed to date to test selection in hybrids using dominant markers. In addition, simulators of dominant markers are very scarce and they do not usually take into account hybridization.

RESULTS

Here, we present a new, multifunctional, R package for dominant genetic markers, AFLPsim. This package can simulate dominant markers in hybridizing populations and implements genome scan methods for detecting outlier dominant loci in hybrids. In addition, it includes tools for further manipulating the results, plotting them and other tasks. We describe and tabulate the major functions implemented in AFLPsim. In addition, we provide some demonstration of its use and we perform a comparative study with other software. Finally, we conclude by briefly describing the input and output formats.

CONCLUSIONS

The R package AFLPsim application provides several useful tools in the context of hybridization studies. It can simulate dominant markers in hybridizing populations and predict their demographic evolution. In addition, we implement a new genome scan method for detecting outlier dominant loci in hybrids, which shows a rather high sensitivity and is very conservative in comparison with Gagnaire et al.'s, Bayescan and introgress. The application is downloadable at http://cran.r-project.org/web/packages/AFLPsim/.

Molecular evidence for natural hybridization between wild loquat (Eriobotrya japonica) and its relative E. prinoides

BACKGROUND

Interspecific hybridization has long been recognized as a pivotal process in plant evolution and speciation. It occurs fairly common in the genera of the subtribe Pyrinae. In Eriobotrya, a small tree genus of Pyrinae, E. prinoides var. daduheensis has been recognized as either a variety of E. prinoides, a natural hybrid between E. prinoides and E. japonica, or a variety of E. japonica. However, to date, there has been no convincing evidence on its status.

RESULTS

Four nuclear genes and two chloroplast regions were sequenced in 89 individuals of these three Eriobotrya taxa from two locations where they coexist. A few fixed nucleotide substitutions or gaps were found in each of the investigated nuclear and chloroplast loci between E. japonica and E. prinoides. Of the 35 individuals of E. prinoides var. daduheensis, 33 showed nucleotide additivity of E. japonica and E. prinoides in at least one nuclear gene, and 10 of them harboured nucleotide additivity at all the four nuclear genes. Most haplotypes of E. prinoides var. daduheensis were also shared with those of E. japonica and E. prinoides. In the two chloroplast regions, 28 and 7 individuals were identical with E. japonica and E. prinoides, respectively.

CONCLUSIONS

Our study provides compelling evidence for a hybrid status for E. prinoides var. daduheensis. Most hybrid individuals are later-generation hybrids. Both E. japonica and E. prinoides can serve as female parent. Differential adaptation might maintain the species boundary of E. prinoides and E. japonica in the face of hybridization and potential introgression.

Unidirectional hybridization and reproductive barriers between two heterostylous primrose species in north-west Yunnan, China

BACKGROUND AND AIMS

Heteromorphy in flowers has a profound effect on breeding patterns within a species, but little is known about how it affects reproductive barriers between species. The heterostylous genus Primula is very diverse in the Himalaya region, but hybrids there have been little researched. This study examines in detail a natural hybrid zone between P. beesiana and P. bulleyana.

METHODS

Chloroplast sequencing, AFLP (amplified fragment length polymorphism) markers and morphological comparisons were employed to characterize putative hybrids in the field, using synthetic F1s from hand pollination as controls. Pollinator visits to parent species and hybrids were observed in the field. Hand pollinations were conducted to compare pollen tube growth, seed production and seed viability for crosses involving different morphs, species and directions of crossing.

KEY RESULTS

Molecular data revealed all hybrid derivatives examined to be backcrosses of first or later generations towards P. bulleyana: all had the chloroplast DNA (cpDNA) of this species. Some individuals had morphological traits suggesting they were hybrids, but they were genetically similar to P. bulleyana; they might have been advanced generation backcrosses. Viable F1s could not be produced with P. bulleyana pollen on P. beesiana females, irrespective of the flower morphs used. Within-morph crosses for each species had very low (<10 %) seed viability, whereas crosses between pin P. bulleyana (female) and pin P. beesiana had a higher seed viability of 30 %. Thus genetic incompatibility mechanisms back up mechanical barriers to within-morph crosses in each species, but are not the same between the two species. The two species share their main pollinators, and pollinators were observed to fly between P. bulleyana and hybrids, suggesting that pollinator behaviour may not be an important isolating factor.

CONCLUSIONS

Hybridization is strongly asymmetric, with P. bulleyana the only possible mother and all detected hybrids being backcrosses in this direction. Partial ecological isolation and inhibition of heterospecific pollen, and possibly complete barriers to F1 formation on P. beesiana, may be enough to make F1 formation very rare in these species. Therefore, with no F1 detected, this hybrid zone may have a finite life span as successive generations become more similar to P. bulleyana.

Evolutionary diversifications of plants on the Qinghai-Tibetan Plateau

The Qinghai-Tibetan Plateau (QTP) is the highest and one of the most extensive plateaus in the world. Phylogenetic, phylogeographic, and ecological studies support plant diversifications on the QTP through multiple mechanisms such as allopatric speciation via geographic isolation, climatic oscillations and divergences, pollinator-mediated isolation, diploid hybridization and introgression, and allopolyploidy. These mechanisms have driven spectacular radiations and/or species diversifications in various groups of plants such as Pedicularis L., Saussurea DC., Rhododendron L., Primula L., Meconopsis Vig., Rhodiola L., and many lineages of gymnosperms. Nevertheless, much work is needed toward understanding the evolutionary mechanisms of plant diversifications on the QTP. Well-sampled biogeographic analyses of the QTP plants in the broad framework of the Northern Hemisphere as well as the Southern Hemisphere are still relatively few and should be encouraged in the next decade. This paper reviews recent evidence from phylogenetic and biogeographic studies in plants, in the context of rapid radiations, mechanisms of species diversifications on the QTP, and the biogeographic significance of the QTP in the broader context of both the Northern and Southern Hemisphere biogeography. Integrative multidimensional analyses of phylogeny, morphological innovations, geography, ecology, development, species interactions and diversifications, and geology are needed and should shed insights into the patterns of evolutionary assembly and radiations in this fascinating region.

Morphological and molecular evidence for natural hybridization in sympatric population of Roscoea humeana and R. cautleoides (Zingiberaceae)

The sympatric occurrence of some species in Roscoea is very common, but little information is available on natural hybridization. However, some intermediate individuals were found on the sympatric population of Roscoea humeana and R. cautleoides at Ganhaizi population in northwestern Yunnan Province, China. We suspected that these intermediate individuals were the hybrids of R. humeana and R. cautleoides from the previous evidence, but could not confirm them. In this study, morphometric analysis was followed by examination of HAT-RAPD polymorphisms to determine the occurrence of natural hybridization between sympatric R. humeana and R. cautleoides. The results showed that most morphological characters of the putative hybrids were found to be intermediate between those of R. humeana and R. cautleoides. Meanwhile, molecular analysis confirmed that the morphological intermediates were derived from hybridization between the two species. From the analysis of the NewHybrids, the hybridization individuals were mainly F₁s. These results indicated that interspecific hybridization between R. humeana and R. cautleoides indeed occurred in sympatric population.

Unusual patterns of hybridization involving a narrow endemic Rhododendron species (Ericaceae) in Yunnan, China

PREMISE OF THE STUDY

One potential threat to rare species is genetic swamping caused by hybridization, but few studies have quantified this threat. Rhododendron cyanocarpum is a narrow endemic species that occurs sympatrically with potentially interfertile congeners throughout its range within Yunnan, China. We searched the entire distribution of R. cyanocarpum for hybrids and examined the patterns of hybridization to assess potential threat from hybridization. •

METHODS

In a comprehensive field survey, we detected only one instance of hybridization involving R. cyanocarpum, with R. delavayi, at Huadianba near Dali. Material of both species and putative hybrids was examined using morphology, chloroplast DNA, nuclear ribosomal DNA, and Bayesian analysis of AFLP profiles. •

KEY RESULTS

Of 10 putative hybrids, two were F(1)(')s and at least seven were F(2)(')s. Four backcrosses to R. delavayi were detected among material with R. delavayi-like morphology within the hybrid zone. Backcrosses to R. cyanocarpum were not detected. Therefore F(2)(')s outnumbered all other classes within the hybrid zone, a situation not previously confirmed for plants and extremely rare generally. Hybridization was asymmetrical, with R. delavayi as the maternal parent in all but one of the hybrids detected. •

CONCLUSIONS

Although natural hybridization is common in Rhododendron, it is rare in R. cyanocarpum and is apparently not accompanied by backcrossing toward R. cyanocarpum. Hence, there is no immediate risk of genetic swamping, unless habitat disturbance increases and changes the patterns of hybridization. Our study is the first to report a plant hybrid zone dominated by F(2) hybrids. This pattern might contribute to species barrier maintenance.

Natural hybridization between Rhododendron delavayi and R. cyanocarpum (Ericaceae), from morphological, molecular and reproductive evidence

The natural hybridization that occurs between two sympatric species of Rhododendron subgenus Hymenanthes in Yunnan, China, was investigated. In field observations, it was noted that the putative hybrids between R. delavayi Franch. and R. cyanocarpum (Franch.) Franch. ex W.W. Sm. had intermediate morphologies. On the basis of morphology, chloroplast DNA (trnL-rpl32) and nuclear DNA (waxy), hybrids and parental species were identified. Hybridization occurred in both directions, but was asymmetrical, with R. delavayi as the major maternal parent. Reciprocal hand pollination treatments showed that either species, as pollen donor or pollen receiver, could produce fruits. It was noted that fruit set varied among treatments. The same pollinators (bumblebees) were shared in both parental species. From these results, we conclude that individuals with intermediate morphologies are indeed of hybrid origin from natural hybridization between R. cyanocarpum and R. delavayi. Furthermore, we presume the hybridization at the study site could have been initiated by habitat disturbance in the 1950s, and we may hence witness the early stages of hybrid swarm formation.

Testing hybridization hypotheses and evaluating the evolutionary potential of hybrids in mangrove plant species

Natural hybridization is of marked importance from global to local biological diversity. In mangroves, species ranges overlap extensively with one another and species share a long overlap of flowering time. Although hybridization has been suggested, patterns of hybridization and the evolutionary potential of hybrids are not yet fully understood. This study provides molecular evidence for the parental origins and status of hybrids in the dominant mangrove genus Rhizophora based on comparisons of chloroplast and nuclear phylogenies and estimations of genetic relatedness and structure from inter-simple sequence repeat (ISSR) markers. Phylogenetic analyses indicate that almost all species can act as maternal parents to hybrids and that hybridization can be bidirectional. Bayesian analyses indicate that hybrids are simple F(1) s, and no trace of backcrossing was detected within populations. Hybridization, for the most part, occurs almost only locally and dispersal of hybrid individuals is limited beyond the hybrid sites.