Diploid hybrid origin ofOstryopsis intermedia(Betulaceae) in the Qinghai-Tibet Plateau triggered by Quaternary climate change

Documenting homoploid hybrid speciation

Homoploid hybrid speciation is challenging to document because hybridization can lead to outcomes other than speciation. Thus, some authors have argued that establishment of homoploid hybrid speciation should include evidence that reproductive barriers isolating the hybrid neo-species from its parental species were derived from hybridization. While this criterion is difficult to satisfy, several recent papers have successfully employed a common pipeline to identify candidate genes underlying such barriers and (in one case) to validate their function. We describe this pipeline, its application to several plant and animal species and what we have learned about homoploid hybrid speciation as a consequence. We argue that - given the ubiquity of admixture and the polygenic basis of reproductive isolation - homoploid hybrid speciation could be much more common and more protracted than suggested by earlier conceptual arguments and theoretical studies.

Cryptic diversity and rampant hybridization in annual gentians on the Qinghai-Tibet Plateau revealed by population genomic analysis

Understanding the evolutionary and ecological processes involved in population differentiation and speciation provides critical insights into biodiversity formation. In this study, we employed 29,865 single nucleotide polymorphisms (SNPs) and complete plastomes to examine genomic divergence and hybridization in Gentiana aristata, which is endemic to the Qinghai-Tibet Plateau (QTP) region. Genetic clustering revealed that G. aristata is characterized by geographic genetic structures with five clusters (West, East, Central, South and North). The West cluster has a specific morphological character (i.e., blue corolla) and higher values of FST compared to the remaining clusters, likely the result of the geological barrier formed by the Yangtze River. The West cluster diverged from the other clusters in the Early Pliocene; these remaining clusters diverged from one another in the Early Quaternary. Phylogenetic reconstructions based on SNPs and plastid data revealed substantial cyto-nuclear conflicts. Genetic clustering and D-statistics demonstrated rampant hybridization between the Central and North clusters, along the Bayankala Mountains, which form the geological barrier between the Central and North clusters. Species distribution modeling demonstrated the range of G. aristata expanded since the Last Interglacial period. Our findings provide genetic and morphological evidence of cryptic diversity in G. aristata, and identified rampant hybridization between genetic clusters along a geological barrier. These findings suggest that geological barriers and climatic fluctuations have an important role in triggering diversification as well as hybridization, indicating that cryptic diversity and hybridization are essential factors in biodiversity formation within the QTP region.

Insights into the phylogenetic relationships and species boundaries of the Myricaria squamosa complex (Tamaricaceae) based on the complete chloroplast genome

Myricaria plants are widely distributed in Eurasia and are helpful for windbreak and embankment protection. Current molecular evidence has led to controversy regarding species boundaries within the Myricaria genus and interspecific phylogenetic relationships between three specific species-M. bracteata, M. paniculata and M. squamosa-which have remained unresolved. This study treated these three unresolved taxa as a species complex, named the M. squamosa complex. The genome skimming approach was used to determine 35 complete plastome sequences and nuclear ribosomal DNA sequences for the said complex and other closely related species, followed by de novo assembly. Comparative analyses were conducted across Myricaria to identify the genome size, gene content, repeat type and number, SSR (simple sequence repeat) abundance, and codon usage bias of chloroplast genomes. Tree-based species delimitation results indicated that M. bracteata, M. paniculata and M. squamosa could not be distinguished and formed two monophyletic lineages (P1 and P2) that were clustered together. Compared to plastome-based species delimitation, the standard nuclear DNA barcode had the lowest species resolution, and the standard chloroplast DNA barcode and group-specific barcodes delimitated a maximum of four out of the five species. Plastid phylogenomics analyses indicated that the monophyletic M. squamosa complex is comprised of two evolutionarily significant units: one in the western Tarim Basin and the other in the eastern Qinghai-Tibet Plateau. This finding contradicts previous species discrimination and promotes the urgent need for taxonomic revision of the threatened genus Myricaria. Dense sampling and plastid genomes will be essential in this effort. The super-barcodes and specific barcode candidates outlined in this study will aid in further studies of evolutionary history.

Multiple paternally inherited chloroplast capture events associated with Taxus speciation in the Hengduan Mountains

Mountainous regions provide a multitude of habitats and opportunities for complex speciation scenarios. Hybridization leading to chloroplast capture, which can be revealed by incongruent phylogenetic trees, is one possible outcome. Four allopatric Taxus lineages (three species and an undescribed lineage) from the Hengduan Mountains, southwest China, exhibit conflicting phylogenetic relationships between nuclear and chloroplast phylogenies. Here, we use multi-omic data at the population level to investigate their historical speciation processes. Population genomic analysis based on ddRAD-seq data revealed limited contemporary inter-specific gene flow involving only populations located close to another species. In a historical context, chloroplast and nuclear data (transcriptome) consistently showed conflicting phylogenetic relationships for T. florinii and the Emei type lineage. ILS and chloroplast recombination were excluded as possible causes, and transcriptome and ddRAD-seq data revealed an absence of the mosaic nuclear genomes that characterize hybrid origin scenarios. Therefore, T. florinii appears to have originated when a lineage of T. florinii captured the T. chinensis plastid type, whereas plastid introgression in the opposite direction generated the Emei Type. All four species have distinct ecological niche based on community investigations and ecological niche analyses. We propose that the origins of both species represent very rare examples of chloroplast capture events despite the paternal cpDNA inheritance of gymnosperms. Specifically, allopatrically and/or ecologically diverged parental species experienced a rare secondary contact, subsequent hybridization and reciprocal chloroplast capture, generating two new lineages, each of which acquired a unique ecological niche. These events might have been triggered by orogenic activities of the Hengduan Mountains and an intensification of the Asian monsoon in the late Miocene, and may represent a scenario more common in these mountains than presently known.

Population structure and phylogeography of three closely related tree peonies

Paeonia decomposita, Paeonia rotundiloba, and Paeonia rockii are three closely related species of Sect. Moutan is distributed in the montane area of the Eastern Hengduan Mountain region. Understanding the population history of these three tree peony species could contribute to unraveling the evolutionary patterns of undergrowth species in this hotspot area. We used one nuclear DNA marker (internal transcribed spacer region, ITS) and two chloroplast DNA markers (matK, ycf1) to reconstruct the phylogeographic pattern of the populations. In total, 228 individuals from 17 populations of the three species were analyzed in this study. Three nuclear clades (Clade I - Clade III) and four maternal clades (Clade A - Clade D) were reconstructed. Molecular dating suggested that young lineages diverged during the late Pliocene and early Pleistocene, younger than the uplift of the Hengduan Mountains but older than the last glacial maximum (LGM). Significant population and phylogeographic structures were detected at both markers. Furthermore, the populations of these tree peonies were overall at equilibrium during the climatic oscillations of the Pleistocene. The simulated palaeoranges of the three species during the LGM period mostly overlapped, which could have led to cross-breeding events. We propose an evolutionary scenario in which mountain orogenesis around the Hengduan Mountain area triggered parapatric isolation between maternal lineages of tree peonies. Subsequent climatic fluctuations drove migration and range recontact of these populations along the valleys. This detailed evolutionary history provides new insights into the phylogeographic pattern of species from mountain-valley systems.

Plastome phylogenomics provide new perspective into the phylogeny and evolution of Betulaceae (Fagales)

BACKGROUND

Betulaceae is a relatively small but morphologically diverse family, with many species having important economic and ecological values. Although plastome structure of Betulaceae has been reported sporadically, a comprehensive exploration for plastome evolution is still lacking. Besides, previous phylogenies had been constructed based on limited gene fragments, generating unrobust phylogenetic framework and hindering further studies on divergence ages, biogeography and character evolution. Here, 109 plastomes (sixteen newly assembled and 93 previously published) were subject to comparative genomic and phylogenomic analyses to reconstruct a robust phylogeny and trace the diversification history of Betulaceae.

RESULTS

All Betulaceae plastomes were highly conserved in genome size, gene order, and structure, although specific variations such as gene loss and IR boundary shifts were revealed. Ten divergent hotspots, including five coding regions (P_i > 0.02) and five noncoding regions (P_i > 0.035), were identified as candidate DNA barcodes for phylogenetic analysis and species delimitation. Phylogenomic analyses yielded high-resolution topology that supported reciprocal monophyly between Betula and Alnus within Betuloideae, and successive divergence of Corylus, Ostryopsis, and Carpinus-Ostrya within Coryloideae. Incomplete lineage sorting and hybridization may be responsible for the mutual paraphyly between Ostrya and Carpinus. Betulaceae ancestors originated from East Asia during the upper Cretaceous; dispersals and subsequent vicariance accompanied by historical environment changes contributed to its diversification and intercontinental disjunction. Ancestral state reconstruction indicated the acquisition of many taxonomic characters was actually the results of parallel or reversal evolution.

CONCLUSIONS

Our research represents the most comprehensive taxon-sampled and plastome-level phylogenetic inference for Betulaceae to date. The results clearly document global patterns of plastome structural evolution, and established a well-supported phylogeny of Betulaceae. The robust phylogenetic framework not only provides new insights into the intergeneric relationships, but also contributes to a perspective on the diversification history and evolution of the family.

Genomic divergence between two sister Ostrya species through linked selection and recombination

Studying the evolution of genomic divergence between lineages is a topical issue in evolutionary biology. However, the evolutionary forces that shape the heterogeneous divergence of the genomic landscape are still poorly understood. Here, two wind-pollinated sister-species (Ostrya japonica and O. chinensis) are used to explore what these potential forces might be. A total of 40 individuals from 16 populations across their main distribution areas in China were sampled for genome-wide resequencing. Population demography analyses revealed that these two sister-species diverged at 3.06-4.43 Mya. Both population contraction and increased gene flow were detected during glacial periods, suggesting secondary contact at those times. All three parameters (D XY, π, and ρ) decreased in those regions showing high levels of differentiation (F ST). These findings indicate that linked selection and recombination played a key role in the genomic heterogeneous differentiation between the two Ostrya species. Genotype-environment association analyses showed that precipitation was the most important ecological factor for speciation. Such environmentally related genes and positive selection genes may have contributed to local adaptation and the maintenance of species boundaries.

Response of distribution patterns of two closely related species in Taxus genus to climate change since last inter-glacial

Climate change affects the species spatio-temporal distribution deeply. However, how climate affects the spatio-temporal distribution pattern of related species on the large scale remains largely unclear. Here, we selected two closely related species in Taxus genus Taxus chinensis and Taxus mairei to explore their distribution pattern. Four environmental variables were employed to simulate the distribution patterns using the optimized Maxent model. The results showed that the highly suitable area of T. chinensis and T. mairei in current period was 1.616 × 105 km2 and 3.093 × 105 km2, respectively. The distribution area of T. chinensis was smaller than that of T. mairei in different periods. Comparison of different periods shown that the distribution area of the two species was almost in stasis from LIG to the future periods. Temperature and precipitation were the main climate factors that determined the potential distribution of the two species. The centroids of T. chinensis and T. mairei were in Sichuan and Hunan provinces in current period, respectively. In the future, the centroid migration direction of the two species would shift towards northeast. Our results revealed that the average elevation distribution of T. chinensis was higher than that of T. mairei. This study sheds new insights into the habitat preference and limiting environment factors of the two related species and provides a valuable reference for the conservation of these two threatened species.

Molecular signatures of parallel adaptive divergence causing reproductive isolation and speciation across two genera

Parallel evolution of reproductive isolation (PERI) provides strong evidence for natural selection playing a fundamental role in the origin of species. However, PERI has been rarely demonstrated for well established species drawn from different genera. In particular, parallel molecular signatures for the same genes in response to similar habitat divergence in such different lineages is lacking. Here, based on whole-genome sequencing data, we first explore the speciation process in two sister species of Carpinus (Betulaceae) in response to divergence for temperature and soil-iron concentration in habitats they occupy in northern and southwestern China, respectively. We then determine whether parallel molecular mutations occur during speciation in this pair of species and also in another sister-species pair of the related genus, Ostryopsis, which occupy similarly divergent habitats in China. We show that gene flow occurred during the origin of both pairs of sister species since approximately 9.8 or approximately 2 million years ago, implying strong natural selection during divergence. Also, in both species pairs we detected concurrent positive selection in a gene (LHY) for flowering time and in two paralogous genes (FRO4 and FRO7) of a gene family known to be important for iron tolerance. These changes were in addition to changes in other major genes related to these two traits. The different alleles of these particular candidate genes possessed by the sister species of Carpinus were functionally tested and indicated likely to alter flowering time and iron tolerance as previously demonstrated in the pair of Ostryopsis sister species. Allelic changes in these genes may have effectively resulted in high levels of prezygotic reproductive isolation to evolve between sister species of each pair. Our results show that PERI can occur in different genera at different timescales and involve similar signatures of molecular evolution at genes or paralogues of the same gene family, causing reproductive isolation as a consequence of adaptation to similarly divergent habitats.

•

PERI provides strong evidence for natural selection playing a fundamental role in the origin of species

•

PERI is rarely demonstrated for well-established species drawn from different genera

•

We detected PERI across two genera (Carpinus and Ostryopsis) in the family Betulaceae

•

PERI can occur in different genera at different timescales and involve molecular signatures at similar pathways

The phylogeographic history of Megistostegium (Malvaceae) in the dry, spiny thickets of southwestern Madagascar using RAD-seq data and ecological niche modeling

The spiny thicket of southwestern Madagascar represents an extreme and ancient landscape with extraordinary levels of biodiversity and endemism. Few hypotheses exist for explaining speciation in the region and few plant studies have explored hypotheses for species diversification. Here, we investigate three species in the endemic genus Megistostegium (Malvaceae) to evaluate phylogeographic structure and explore the roles of climate, soil, and paleoclimate oscillations on population divergence and speciation throughout the region. We combine phylogenetic and phylogeographic inference of RADseq data with ecological niche modeling across space and time. Population structure is concurrent with major rivers in the region and we identify a new, potentially important biogeographic break coincident with several landscape features. Our data further suggests that niches occupied by species and populations differ substantially across their distribution. Paleodistribution modeling provide evidence that past climatic change could be responsible for the current distribution, population structure, and maintenance of species in Megistostegium.

Comparative phylogeography study reveals introgression and incomplete lineage sorting during rapid diversification of Rhodiola

BACKGROUND AND AIMS

Many plant taxa in the Qinghai-Tibetan Plateau (QTP) and the Hengduan Mountains (HM) radiated rapidly during the Quaternary but with frequent secondary contact between diverging populations. Incomplete lineage sorting and introgressive hybridization might be involved during the rapid radiation, but their effects on phylogeography have not been fully determined.

METHODS

We investigated the chloroplast DNA (cpDNA)/internal transcribed spacer (ITS) sequence variations of 611 samples of Rhodiola bupleuroides, R. discolor, R. fastigiata and R. chrysanthemifolia from the QTP and HM to compare the phylogeographic patterns between the four species with different evolutionary histories, geographic ranges and reproductive modes.

KEY RESULTS

The divergence times of these species were consistent with the last peak of in situ speciation in the HM. While closely related species exhibited different phylogeographic patterns, they shared several ribotypes and haplotypes in sympatric populations, suggesting introgressive hybridization. A significant phylogenetic discordance between ribotypes and haplotypes was detected in three species, implying incomplete lineage sorting. Rhodiola discolor houses an extraordinary richness of cpDNA haplotypes, and this finding may be attributed to adaptive radiation.

CONCLUSION

In addition to geographic isolation and climate oscillations during the Quaternary, both introgressive hybridization and incomplete lineage sorting play important roles in species that experienced rapid diversification in the QTP and HM.

Can Ecological Farming Systems Positively Affect Household Income from Agriculture? A Case Study of the Suburban Area of Hanoi, Vietnam

From 1995, Vietnam started building safety standards for agricultural products divided into three labels, “safe”, “VietGAP”, and “organic” to deal with environmental problems and health issues resulting from the intensification of farming methods and increasing awareness of food safety. This research examines the reasons for government support of agricultural production according to those standards instead of conventional methods. To this end, we characterised the current farming systems, identified factors affecting the generation of cash income from agriculture, and investigated the influence of the farming system on farm income in the suburban area of Hanoi. Data collection involved a survey based on a structured questionnaire, key informants’ interviews, and focus group discussions conducted in 2015 with a sample size of 312 respondents. Three forms of farming systems: conventional, safe, and organic were identified. Literature suggested that ecological farming systems that are safe and organic are supported by stakeholders in Vietnam to maintain agricultural sustainability. A multiple linear regression model was applied to identify factors that affect cash incomes from agricultural activities. This research found a direct correlation between the ecological factor of farming systems and respondents’ cash income from agricultural activities. Additionally, the amount of family labour and respondents’ higher-level education and farm experiences appeared to positively influence the cash incomes from agriculture. On the other hand, farm size and membership in a farmers’ association were identified as the major factors that negatively affect agricultural cash income in the study area.

Species divergence and phylogeography of Corylus heterophylla Fisch complex (Betulaceae): Inferred from molecular, climatic and morphological data

Historical geo-climatic changes have shaped the geographical distributions and genetic diversity of numerous plant taxa in East Asia, which promote species divergence and ultimately speciation. Here, we integrated multiple approaches, including molecular phylogeography, ecological niche modeling, and morphological traits to examine the nucleotide diversity and interspecific divergence within Corylus heterophylla complex (C. heterophylla, C. kweichowensis, and C. yunnanensis). These three sibling taxa harbored similar high levels of nucleotide diversity at the species level. The molecular data (SCNG and cpDNA) unanimously supported the division of C. heterophylla complex into two major clades, with C. yunnanensis diverged earlier from the complex, whereas C. heterophylla and C. kweichowensis could hardly be separated. The split between the two clades (c. 12.89 Ma) coincided with the formation of Sichuan Basin in the middle Miocene, while the divergence among and within the five subclades (YUN1-YUN3, HK1-HK2) occurred from the late Miocene to the Pleistocene. C. heterophylla of northern China experienced glacial contraction and interglacial expansion during the Quaternary, whereas C. kweichowensis and C. yunnanensis of southern China presented population expansion even during the last glacial maximum. Despite of high levels of genetic admixture between C. heterophylla and C. kweichowensis, significant ecological and morphological discrepancy as well as incomplete geographic isolation indicated that adaptive evolution triggered by divergent selection may have played important roles in incipient ecological speciation.

When tropical and subtropical congeners met: Multiple ancient hybridization events within Eriobotrya in the Yunnan-Guizhou Plateau, a tropical-subtropical transition area in China

Global climate changes during the Miocene may have created ample opportunities for hybridization between members of tropical and subtropical biomes at the boundary between these zones. Yet, very few studies have explored this possibility. The Yunnan-Guizhou Plateau (YGP) in Southwest China is a biodiversity hotspot for vascular plants, located in a transitional area between the floristic regions of tropical Southeast Asia and subtropical East Asia. The genus Eriobotrya (Rosaceae) comprises both tropical and subtropical taxa, with 12 species recorded in the YGP, making it a suitable basis for testing the hypothesis of between-biome hybridization. Therefore, we surveyed the evolutionary history of Eriobotrya by examining three chloroplast regions and five nuclear genes for 817 individuals (47 populations) of 23 Eriobotrya species (including 19 populations of 12 species in the YGP), plus genome re-sequencing of 33 representative samples. We concluded that: (1) phylogenetic positions for 16 species exhibited strong cytonuclear conflicts, most probably due to ancient hybridization; (2) the YGP is a hotspot for hybridization, with 11 species showing clear evidence of chloroplast capture; and (3) Eriobotrya probably originated in tropical Asia during the Eocene. From the Miocene onwards, the intensification of the Eastern Asia monsoon and global cooling may have shifted the tropical-subtropical boundary and caused secondary contact between species, thus providing ample opportunity for hybridization and diversification of Eriobotrya, especially in the YGP. Our study highlights the significant role that paleoclimate changes probably played in driving hybridization and generating rich species diversity in climate transition zones.

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.

Recurrent hybridization underlies the evolution of novelty in Gentiana (Gentianaceae) in the Qinghai-Tibetan Plateau

The Qinghai-Tibetan Plateau (QTP) and adjacent areas are centres of diversity for several alpine groups. Although it is known that the QTP acted as a source area for diversification of the alpine genus Gentiana, the evolutionary processes underlying diversity in this genus, especially the formation of narrow endemics, are still poorly understood. Hybridization has been proposed as a driver of plant endemism in the QTP but few cases have been documented with genetic data. Here, we describe a new endemic species in Gentiana section Cruciata as G. hoae sp. nov., and explore its evolutionary history with complete plastid genomes and nuclear ribosomal internal transcribed spacer sequence data. Genetic divergence within G. hoae ~3 million years ago was followed by postglacial expansion on the QTP, suggesting Pleistocene glaciations as a key factor shaping the population history of G. hoae. Furthermore, a mismatch between plastid and nuclear data suggest that G. hoae participated in historical hybridization, while population sequencing show this species continues to hybridize with the co-occurring congener G. straminea in three locations. Our results indicate that hybridization may be a common process in the evolution of Gentiana and may be widespread among recently diverged taxa of the QTP.

Hybrid speciation via inheritance of alternate alleles of parental isolating genes

It is increasingly realized that homoploid hybrid speciation (HHS), which involves no change in chromosome number, is an important mechanism of speciation. HHS will likely increase in frequency as ecological and geographical barriers between species are continuing to be disrupted by human activities. HHS requires the establishment of reproductive isolation between a hybrid and its parents, but the underlying genes and genetic mechanisms remain largely unknown. In this study, we reveal by integrated approaches that reproductive isolation originates in one homoploid hybrid plant species through the inheritance of alternate alleles at genes that determine parental premating isolation. The parent species of this hybrid species are reproductively isolated by differences in flowering time and survivorship on soils containing high concentrations of iron. We found that the hybrid species inherits alleles of parental isolating major genes related to flowering time from one parent and alleles of major genes related to iron tolerance from the other parent. In this way, it became reproductively isolated from one parent by the difference in flowering time and from the other by habitat adaptation (iron tolerance). These findings and further modeling results suggest that HHS may occur relatively easily via the inheritance of alternate parental premating isolating genes and barriers.

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.

Phylogenomic analyses based on genome-skimming data reveal cyto-nuclear discordance in the evolutionary history of Cotoneaster (Rosaceae)

As a consequence of hybridization, polyploidization, and apomixis, the genus Cotoneaster (Rosaceae) represents one of the most complicated and controversial lineages in Rosaceae, with ca. 370 species which have been classified into two subgenera and several sections, and is notorious for its taxonomic difficulty. The infrageneric relationships and taxonomy of Cotoneaster have remained poorly understood. Previous studies have focused mainly on natural hybridization involving only several species, and phylogeny based on very limited markers. In the present study, the sequences of complete chloroplast genomes and 204 low-copy nuclear genes of 72 accessions, representing 69 species as ingroups, were used to conduct the most comprehensive phylogenetic analysis so far for Cotoneaster. Based on the sequences of complete chloroplast genomes and many nuclear genes, our analyses yield two robust phylogenetic trees respectively. Chloroplast genome and nuclear data confidently resolved relationships of this genus into two major clades which largely supported current classification based on morphological evidence. However, conflicts between the chloroplast genome and low-copy nuclear phylogenies were observed in both the species level and clade level. Cyto-nuclear discordance in the phylogeny could be caused by frequent hybridization events and incomplete sorting lineage (ILS). In addition, our divergence-time analysis revealed an evolutionary radiation of the genus from late Miocene to date.

Asymmetric Reproductive Barriers and Gene Flow Promote the Rise of a Stable Hybrid Zone in the Mediterranean High Mountain

Hybrid zones have the potential to shed light on evolutionary processes driving adaptation and speciation. Secondary contact hybrid zones are particularly powerful natural systems for studying the interaction between divergent genomes to understand the mode and rate at which reproductive isolation accumulates during speciation. We have studied a total of 720 plants belonging to five populations from two Erysimum (Brassicaceae) species presenting a contact zone in the Sierra Nevada mountains (SE Spain). The plants were phenotyped in 2007 and 2017, and most of them were genotyped the first year using 10 microsatellite markers. Plants coming from natural populations were grown in a common garden to evaluate the reproductive barriers between both species by means of controlled crosses. All the plants used for the field and greenhouse study were characterized by measuring traits related to plant size and flower size. We estimated the genetic molecular variances, the genetic differentiation, and the genetic structure by means of the F-statistic and Bayesian inference. We also estimated the amount of recent gene flow between populations. We found a narrow unimodal hybrid zone where the hybrid genotypes appear to have been maintained by significant levels of a unidirectional gene flow coming from parental populations and from weak reproductive isolation between them. Hybrid plants exhibited intermediate or vigorous phenotypes depending on the analyzed trait. The phenotypic differences between the hybrid and the parental plants were highly coherent between the field and controlled cross experiments and through time. The highly coherent results obtained by combining field, experimental, and genetic data demonstrate the existence of a stable and narrow unimodal hybrid zone between Erysimum mediohispanicum and Erysimum nevadense at the high elevation of the Sierra Nevada mountains.

RAD-sequencing reveals patterns of diversification and hybridization, and the accumulation of reproductive isolation in a clade of partially sympatric, tropical island trees

A common pattern observed in temperate tree clades is that species are often morphologically distinct and partially interfertile but maintain species cohesion despite ongoing hybridization where ranges overlap. Although closely related species commonly occur in sympatry in tropical ecosystems, little is known about patterns of hybridization within a clade over time, and the implications of this hybridization for the maintenance of species boundaries. In this study, we focused on a clade of sympatric trees in the genus Diospyros in the Mascarene islands and investigated whether species are genetically distinct, whether they hybridize, and how patterns of hybridization are related to the time since divergence among species. We sampled multiple populations from each of 12 Mascarene Diospyros species, generated genome-wide single nucleotide polymorphism data using 2bRADseq, and conducted population genomic and phylogenomic analyses. We found that Mascarene Diospyros species diverged millions of years ago and are today largely genetically distinct from one another. Although hybridization was observed between closely related species belonging to the same subclade, more distantly related species showed little evidence of interspecific hybridization. Phylogenomic analyses also suggested that introgression has occurred during the evolutionary history of the clade. This suggests that, as diversification progressed, interspecific hybridization occurred among species, but became infrequent as lineages diverged from one another and evolved reproductive barriers. Species now coexist in partial sympatry, and experience limited hybridization between close relatives. Additional research is needed to better understand the role that introgression may have played in adaptation and diversification of Mascarene Diospyros, and its relevance for conservation.

Hybridization Between Yuccas From Baja California: Genomic and Environmental Patterns

Hybridization can occur when two geographically isolated species are reproductively compatible and have come into sympatry due to range shifts. Yucca and yucca moths exhibit obligate pollination mutualism; yucca moths are responsible for the gene flow mediated by pollen among yucca populations. In the Baja California Peninsula, there are two yucca sister species, Y. capensis and Y. valida, that have coevolved with the same pollinator, Tegeticula baja. Both yucca species are endemic to the peninsula, and their current distributions are allopatric. Based on their morphological characteristics, it has been suggested that some plants growing in the southern part of the Magdalena flatland, a spatially disjunct part of Yucca valida's range, have hybrid origins. We conducted genomic and climatic analyses of the two yucca species as well as the putative hybrid populations. We genotyped 3,423 single nucleotide polymorphisms in 120 individuals sampled from 35 localities. We applied Bayesian tests and geographic cline analyses to the genomic data. Using climatic information from the occurrence sites, we projected species distribution models in different periods to assess changes in the distributional range, and we performed a statistical test to define the niche divergence between the paternal species and the putative hybrid populations. Structure analysis revealed mixed ancestry in the genome of hybrid populations, and the Bayesian models supported a scenario of post-divergence gene flow between the yucca species. Our species distribution models reveal that the geographical ranges of the parental species overlapped mainly during the Last Glacial Maximum, which could facilitate genetic admixture between those species. Finally, we found that most of the assessed environmental axes between the parents and hybrid populations are divergent, indicating that the climatic niche of the hybrid populations is shifting from that of the populations' progenitors. Our results show that the populations in the southern part of the Magdalena flatland are the result of combination of the genetic components of two species. Hybrid individuals with this novel genomic combination arose in a different habitat than their parental species, and they exhibit ecological divergence, which contributes to reproductive isolation through spatial and temporal barriers.

Population subdivision and hybridization in a species complex of Gentiana in the Qinghai-Tibetan Plateau

BACKGROUND AND AIMS

Hosting several global biodiversity hotspots, the region of the Qinghai-Tibetan Plateau (QTP) is exceptionally species-rich and harbours a remarkable level of endemism. Yet, despite a growing number of studies, factors fostering divergence, speciation and ultimately diversity remain poorly understood for QTP alpine plants. This is particularly the case for the role of hybridization. Here, we explored the evolutionary history of three closely related Gentiana endemic species, and tested whether our results supported the mountain geo-biodiversity hypothesis (MGH).

METHODS

We genotyped 69 populations across the QTP with one chloroplast marker and 12 nuclear microsatellite loci. We performed phylogeographical analysis, Bayesian clustering, approximate Bayesian computation and principal components analysis to explore their genetic relationship and evolutionary history. In addition, we modelled their distribution under different climates.

KEY RESULTS

Each species was composed of two geographically distinct clades, corresponding to the south-eastern and north-western parts of their distribution. Thus Gentiana veitchiorum and G. lawrencei var. farreri, which diverged recently, appear to have shared at least refugia in the past, from which their range expanded later on. Indeed, climatic niche modelling showed that both species went through continuous expansion from the Last Interglacial Maximum to the present day. Moreover, we have evidence of hybridization in the northwest clade of G. lawrencei var. farreri, which probably occurred in the refugium located on the plateau platform. Furthermore, phylogenetic and population genetic analyses suggested that G. dolichocalyx should be a geographically limited distinct species with low genetic differentiation from G. lawrencei var. farreri.

CONCLUSIONS

Climatic fluctuations in the region of the QTP have played an important role in shaping the current genetic structure of G. lawrencei var. farreri and G. veitchiorum. We argue that a species pump effect did occur prior to the Last Interglacial Maximum, thus lending support to the MGH. However, our results do depart from expectations as suggested in the MGH for more recent distribution range and hybridization dynamics.

Speciation and subsequent secondary contact in two edaphic endemic primroses driven by Pleistocene climatic oscillation

Climatic perturbation during the Pleistocene era has played a major role in plant evolutionary history by altering species distribution range. However, the relative roles of climatic and geographic factors in the distribution dynamics remain poorly understood; in particular, the edaphic endemics. In this paper, we examine the evolutionary history of two ultramafic primroses, Primula hidakana and Primula takedana. These species are ecologically and morphologically distinct with disjunct distributions on Hokkaido Island, Japan. Primula hidakana is found on various rocks in southern Hokkaido and P. takedana in serpentine areas in northern Hokkaido. We performed population genetics analyses on nuclear and chloroplast data sets and tested alternative phylogenetic models of divergence using approximate Bayesian computation (ABC) analyses. Nuclear microsatellite loci clearly distinguished the two sister taxa. In contrast, chloroplast sequence variations were shared between P. takedana and P. hidakana. ABC analyses based on nuclear data supported a secondary contact scenario involving asymmetrical gene flow from P. hidakana to P. takedana. Paleodistribution modeling also supported the divergence model, and predicted their latitudinal range shifts leading to past secondary contact. Our findings highlight the importance of the distribution dynamics during the Pleistocene climatic oscillations in the evolution of serpentine plants, and demonstrate that tight species cohesion between serpentine and nonserpentine sister taxa has been maintained despite past interspecific gene flow across soil boundaries.

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.

Genetic data reveals a complex history of multiple admixture events in presently allopatric wild gingers (Asarum spp.) showing intertaxonomic clinal variation in calyx lobe length

Clinal variation is a major pattern of observed phenotypic diversity and identifying underlying demographic processes is a necessary step to understand the establishment of clinal variation. The wild ginger series Sakawanum (genus Asarum) comprises four taxa, which exhibit intertaxonomic clinal variation in calyx lobe length across two continental islands isolated by a sea strait. To test alternative hypotheses of the evolutionary history and to determine the implications for the formation of clinal variation, we conducted approximate Bayesian computation (ABC) analysis and ecological niche modeling (ENM). ABC analysis indicated that the scenario assuming multiple admixture events was strongly supported. This scenario assumed two admixture events occurred between morphologically distinct taxa, likely leading to the generation of intermediate taxa. One of the admixture events was estimated to have occurred during the last glacial maximum (LGM), during which the taxa were estimated to have formed a common refugia in southern areas by ENM analysis. Although four taxa are currently distributed allopatrically on different islands and trans-oceanic dispersal appears unlikely, the formation of a land bridge and the geographic range shift to refugia would have allowed secondary contact between previously isolated taxa. This study suggests that clinal variation can be shaped by demographic history including multiple admixtures due to climatic oscillations.

Geophysical upheavals and evolutionary diversification of plant species in the Himalaya

The Himalaya is one of the youngest and the loftiest mountain chains of the world; it is also referred to as the water tower of Asia. The Himalayan region harbors nearly 10,000 plant species constituting approximately 2.5% of the global angiosperm diversity of which over 4,000 are endemics. The present-day Himalayan flora consists of an admixture of immigrant taxa and diversified species over the last 40 million years. The interesting questions about the Himalayan flora discussed here are: how did the Himalaya achieve high endemic plant diversity starting with immigrant taxa and what were the main drivers of this diversity? This contribution aims to answer these questions and raise some more. We review and analyze existing information from diverse areas of earth and climate sciences, palaeobiology and phytogeography to evolve a bio-chronological record of plant species divergence and evolution in the Himalaya. From the analysis we infer the effects of major environmental upheavals on plant diversity in the region. The understanding developed in the following discussion is based on the idea that Himalaya experienced at least five phases of major geophysical upheavals, namely: (i) mega-collision between India and Eurasian plates, (ii) tectonic uplift in phases and progressive landform elevation, (iii) onset of southwest (SW) Indian monsoon, (iv) spurring of arid conditions in Central Asia, and (v) cyclic phases of cooling and warming in the Quaternary. The geophysical upheavals that were potentially disrupting for the ecosystem stability had a key role in providing impetus for biological diversification. The upheavals produced new geophysical environments, new ecological niches, imposed physical and physiological isolation barriers, acted as natural selection sieves and led to the formation of new species. This contribution aims to develop a comprehensive understanding of the plant biodiversity profile of the Himalaya in the context of complex, interconnected and dynamic relationship between earth system processes, climate and plant diversity.

Approximate Bayesian computation analysis of EST-associated microsatellites indicates that the broadleaved evergreen tree Castanopsis sieboldii survived the Last Glacial Maximum in multiple refugia in Japan

Climatic changes have played major roles in plants' evolutionary history. Glacial oscillations have been particularly important, but some of their effects on plants' populations are poorly understood, including the numbers and locations of refugia in Asian warm temperate zones. In the present study, we investigated the demographic history of the broadleaved evergreen tree species Castanopsis sieboldii (Fagaceae) during the last glacial period in Japan. We used approximate Bayesian computation (ABC) for model comparison and parameter estimation for the demographic modeling using 27 EST-associated microsatellites. We also performed the species distribution modeling (SDM). The results strongly support a demographic scenario that the Ryukyu Islands and the western parts in the main islands (Kyushu and western Shikoku) were derived from separate refugia and the eastern parts in the main islands and the Japan Sea groups were diverged from the western parts prior to the coldest stage of the Last Glacial Maximum (LGM). Our data indicate that multiple refugia survived at least one in the Ryukyu Islands, and the other three regions of the western and eastern parts and around the Japan Sea of the main islands of Japan during the LGM. The SDM analysis also suggests the potential habitats under LGM climate conditions were mainly located along the Pacific Ocean side of the coastal region. Our ABC-based study helps efforts resolve the demographic history of a dominant species in warm temperate broadleaved forests during and after the last glacial period, which provides a basic model for future phylogeographical studies using this approach.

Historical introgression among the species of Rodgersia (Saxifragaceae) in mountainous forests of southwest China

In the present study, we used genetic data and ecological niche modelling to explore possible historical introgressions among the species of Rodgersia (Saxifragaceae) in central-southwest China. Markedly differentiated chloroplast haplotypes were found in R. aesculifolia, R. sambucifolia and the Lijiang (LJ) population of R. pinnata, respectively, and differentiated chloroplast haplotypes within each of them showed the closest relationships with haplotypes from different species. ITS cloning did not reveal any shared ribotype between R. aesculifolia and the remaining species. Historical introgression between R. aesculifolia and R. sambucifolia (or R. pinnata) seems to be the most plausible explanation according to the geographical pattern and derivative status of putative introgressed chloroplast haplotypes, and also from morphological evidence. Introgressions were also found among R. sambucifolia, R. pinnata, and R. henricii from Yunnan. Frequent gene exchanges may have promoted the diversity of leaf shapes in this genus. Ecological niche modelling indicated that past secondary contact following range shifts during Pleistocene cold periods may have provided opportunities for ancient introgression between R. aesculifolia and adjacent species.

Rapid radiations of both kiwifruit hybrid lineages and their parents shed light on a two-layer mode of species diversification

Reticulate speciation caused by interspecific hybridization is now recognized as an important mechanism in the creation of biological diversity. However, depicting the patterns of phylogenetic networks for lineages that have undergone interspecific gene flow is challenging. Here we sequenced 25 taxa representing natural diversity in the genus Actinidia with an average mapping depth of 26× on the reference genome to reconstruct their reticulate history. We found evidence, including significant gene tree discordance, cytonuclear conflicts, and changes in genome-wide heterozygosity across taxa, collectively supporting extensive reticulation in the genus. Furthermore, at least two separate parental species pairs were involved in the repeated origin of the hybrid lineages, in some of which a further phase of syngameon was triggered. On the basis of the elucidated hybridization relationships, we obtained a highly resolved backbone phylogeny consisting of taxa exhibiting no evidence of hybrid origin. The backbone taxa have distinct demographic histories and are the product of recent rounds of rapid radiations via sorting of ancestral variation under variable climatic and ecological conditions. Our results suggest a mode for consecutive plant diversification through two layers of radiations, consisting of the rapid evolution of backbone lineages and the formation of hybrid swarms derived from these lineages.

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.

Transient hybridization, not homoploid hybrid speciation, between ancient and deeply divergent conifers

PREMISE OF THE STUDY

Homoploid hybrid speciation is receiving growing attention due the increasing recognition of its role in speciation. We investigate if individuals intermediate in morphology between the two species of the conifer genus Athrotaxis represent a homoploid hybrid species, A. laxifolia, or are spontaneous F1 hybrids.

METHODS

A total of 1055 individuals of Athrotaxis cupressoides and A. selaginoides, morphologically intermediate individuals, and two putative hybrid swarms were sampled across the range of the genus and genotyped with 13 microsatellites. We used simulations to test the power of our data to identify the pure species, F1s, F2s, and backcross generations.

KEY RESULTS

We found that Athrotaxis cupressoides and A. selaginoides are likely the most divergent congeneric conifers known, but the intermediates are F1 hybrids, sharing one allele each from A. cupressoides and A. selaginoides at six loci with completely species specific alleles. The hybrid swarms contain wide genetic variation with stronger affinities to the locally dominant species, A. selaginoides and A. selaginoides backcrosses outnumbering A. cupressoides backcrosses. In addition, we observed evidence for isolated advanced generation backcrosses within the range of the pure species.

CONCLUSIONS

We conclude that, even though they can be large and long-lived, Athrotaxis hybrid swarms are on a trajectory of decline and will eventually be reabsorbed by the parental species. However, this process may take millennia and fossil evidence suggests that such events have occurred repeatedly since the early Quaternary. Given this timeline, our study highlights the many obstacles to homoploid hybrid speciation.

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.

DNA barcoding evaluation and its taxonomic implications in the species-rich genus Primula L. in China

The genus Primula is extremely diverse in the east Himalaya-Hengduan Mountains (HHM) in China as a result of rapid radiation. In order to overcome the difficulty of morphological classification of this genus, we surveyed three plastid regions (rbcL, matK, and trnH-psbA) and two nuclear markers (ITS and ITS2) from 227 accessions representing 66 Primula species across 18 sections, to assess their discriminatory power as barcodes. We found that ITS alone or combined with plastid regions showed the best discrimination across different infrageneric ranks and at species level. We suggest rbcL + matK + ITS as the first choice at present to barcode Primula plants. Although the present barcoding combination performed poorly in many closely related species of Primula, it still provided many new insights into current Primula taxonomy, such as the underlying presence of cryptic species, and several potential improper taxonomic treatments. DNA barcoding is one useful technique in the integrative taxonomy of the genus Primula, but it still requires further efforts to improve its effectiveness in some taxonomically challenging groups.