Recent allopolyploidization and transcriptomic asymmetry in the mangrove shrub Acanthus tetraploideus

BACKGROUND

Mangrove species are vital to the ecosystems of tropical and subtropical coastlines worldwide. Despite the underexplored role of polyploidization in these species, deciphering its impact on gene expression is essential for understanding the connection between polyploidization and species diversification. Our initial investigation, integrating multiple nuclear loci with morphological and cytological data, indicates that the tetraploid Acanthus tetraploideus likely originated from allopolyploidization events involving the diploid species A. ilicifolius and A. ebracteatus. Expanding on these insights, this study utilizes genome-wide evidence to confirm the divergence patterns among extant Acanthus mangrove diploids and to investigate the origin and transcriptome asymmetry of the tetraploid A. tetraploideus.

RESULTS

Phylogenetic analysis and molecular dating revealed a closer evolutionary relationship between A. ebracteatus and A. volubilis than between A. ebracteatus and A. ilicifolius, diverged approximately 6.92 Mya and 9.59 Mya, respectively. Analysis of individual whole transcriptomes revealed that homeologous sequences in A. tetraploideus were preferentially clustered with A. ilicifolius and A. ebracteatus, rather than A. volubilis, in a roughly 1:1 ratio. The high similarity in nucleotide sequences and homologous polymorphisms between the tetraploid A. tetraploideus and its two parental diploids, A. ebracteatus and A. ilicifolius, supports the hypothesis of a recent allopolyploid origin for A. tetraploideus. Estimation of homeolog expression revealed a general attenuation of homeolog expression divergence in A. tetraploideus compared to the in silico parental mix, with 22.87% and 67.66% of genes exhibiting biased homeolog expression, respectively. Further investigation identified remarkable retention of parental expression dominance in the tetraploid, suggesting that parental genetic legacy substantially influences the reconfiguration of homeolog expression in the derived tetraploid. Meanwhile, the observation of numerous novel expression patterns between the two homeolog sets suggests that the transcriptome shock (i.e., the transcriptomic changes induced by interspecific hybridization) associated with allopolyploidization and subsequent post-polyploid evolutionary processes also significantly impact transcriptome asymmetry in A. tetraploideus. While no strong evidence directly links transcriptomic changes to specific adaptive traits, the patterns in unbiased and novelly biased genes in A. tetraploideus suggest adaptations to stable polyploidy. Unbiased genes involved in fundamental cellular processes and novelly biased genes related to chromosome dynamics and cell cycle regulation may stabilize polyploid genomes, supporting the species' establishment and long-term success. These findings underscore the role of transcriptomic stability in polyploid adaptation.

CONCLUSIONS

Our study sheds light on the evolutionary origins and the intricate transcriptional reconfiguration of the tetraploid A. tetraploideus. These insights significantly enhance our comprehension of the pivotal role that polyploidization plays in speciation and adaptative evolution of mangrove species.

New estimates and synthesis of chromosome numbers, ploidy levels and genome size variation in Allium sect. Codonoprasum: advancing our understanding of the unresolved diversification and evolution of this section

BACKGROUND

The genus Allium is known for its high chromosomal variability, but most chromosome counts are based on a few individuals and genome size (GS) reports are limited in certain taxonomic groups. This is evident in the Allium sect. Codonoprasum, a species-rich (> 150 species) and taxonomically complex section with weak morphological differences between taxa, the presence of polyploidy and frequent misidentification of taxa. Consequently, a significant proportion of older karyological reports may be unreliable and GS data are lacking for the majority of species within the section. This study, using chromosome counting and flow cytometry (FCM), provides the first comprehensive and detailed insight into variation in chromosome number, polyploid frequency and distribution, and GS in section members, marking a step towards understanding the unresolved diversification and evolution of this group.

RESULTS

We analysed 1578 individuals from 316 populations of 25 taxa and reported DNA ploidy levels and their GS, with calibration from chromosome counts in 22 taxa. Five taxa had multiple ploidy levels. First estimates of GS were obtained for 16 taxa. A comprehensive review of chromosome number and DNA-ploidy levels in 129 taxa of the section revealed that all taxa have x = 8, except A. rupestre with two polyploid series (x = 8, descending dysploidy x = 7), unique for this section. Diploid taxa dominated (72.1%), while di- & polyploid (12.4%) and exclusively polyploid (15.5%) taxa were less common. Ploidy diversity showed that diploid taxa dominated in the eastern Mediterranean and decreased towards the west and north, whereas only polyploid cytotypes of di- & polyploid taxa or exclusively polyploid taxa dominated in northern and northwestern Europe. A 4.1-fold variation in GS was observed across 33 taxa analysed so far (2C = 22.3-92.1 pg), mainly due to polyploidy, with GS downsizing observed in taxa with multiple ploidy levels. Intra-sectional GS variation suggests evolutionary relationships, and intraspecific GS variation within some taxa may indicate taxonomic heterogeneity and/or historical migration patterns.

CONCLUSIONS

Our study showed advantages of FCM as an effective tool for detecting ploidy levels and determining GS within the section. GS could be an additional character in understanding evolution and phylogenetic relationships within the section.

Across a phylogeographic break in the Qinling Mountains-Huaihe River Line: Quaternary evolutionary history of a medicinal and edible homologous plant (Allium macrostemon) in China

Biogeographic barriers to gene flow are central to studies of plant phylogeography. There are many physical and geographic barriers in China, but few studies have used molecular ecological evidence to investigate the natural geographic isolation barrier of the Qinling Mountains-Huaihe River Line (QHL). Allium macrostemon is a precious Chinese perennial herb belonging to the Amaryllidaceae family. It is used as a food and medicine, with a variety of health and healing properties. Five SSR markers, three chloroplast DNA (cpDNA) markers (psbA-trnH, rps16 and trnL-F), one nuclear ribosomal DNA (nrDNA) marker (ITS), and simplified genome GBS sequencing were used to analyse the genetic diversity and structure of A. macrostemon. Combining SSR, cpDNA, nrDNA ITS data and GBS analysis results, we divided A. macrostemon populations into northern and southern groups, with the southern group further divided into southwestern and central-southeastern groups. Niche simulation results reveal that the distribution area of A. macrostemon will reach its maximum in the future. These data indicate that the regional separation of A. macrostemon has been maintained by the combined influence of a geographical barrier and Quaternary climate, and that the back-and-forth fluctuations of QHL and Quaternary climate have played an important role in this process. QHL acts as a north-south dividing line in phylogeography and population genetic structure, promoting physical geographic isolation. This study provides a theoretical basis for the conservation, development, and utilization of A. macrostemon resources. It further provides a reference for understanding the systematic geographical pattern of the large-scale spatial distribution of plants in China and enriches our understanding of Quaternary plant evolution in areas with complex terrain.

Demographic patterns of two related desert shrubs with overlapping distributions in response to past climate changes

Numerous studies have revealed that past geological events and climatic fluctuations had profoundly affected the genetic structure and demographic patterns of species. However, related species with overlapping ranges may have responded to such environmental changes in different ways. In this study, we compared the genetic structure and population dynamics of two typical desert shrubs with overlapping distributions in northern China, Nitraria tangutorum and Nitraria sphaerocarpa, based on chloroplast DNA (cpDNA) variations and species distribution models. We sequenced two cpDNA fragments (trnH-trnA and atpH-atpI) in 633 individuals sampled from 52 natural populations. Twenty-four chlorotypes, including eight rare chlorotypes, were identified, and a single dominant haplotype (H4) widely occurred in the entire geographical ranges of the two species. There were also a few distinctive chlorotypes fixed in different geographical regions. Population structure analyses suggested that the two species had significantly different levels of total genetic diversity and interpopulation differentiation, which was highly likely correlated with the special habitat preferences of the two species. A clear phylogeographic structure was identified to exist among populations of N. sphaerocarpa, but not exist for N. tangutorum. The neutral tests, together with the distribution of pairwise differences revealed that N. tangutorum experienced a sudden demographic expansion, and its expansion approximately occurred between 21 and 7 Kya before present, while a rapid range expansion was not identified for N. sphaerocarpa. The ecological niche modeling (ENM) analysis indicated that the potential ranges of two species apparently fluctuated during the past and present periods, with obvious contraction in the Last Glacial Maximum (LGM) and recolonization in the present, respectively, comparing to the Last Interglacial (LIG). These findings suggest that the two species extensively occurred in the Northwest of China before the Quaternary, and the current populations of them originated from a few separated glacial refugia following their habitat fragmentation in the Quarternary. Our results provide new insights on the impact of past geological and climatic fluctuations on the population dynamics of desert plants in northwestern China, and further enforce the hypothesis that there were several independent glacial refugia for these species during the Quaternary glaciations.

Phylogenomic analyses and chromosome ploidy identification reveal multiple cryptic species in Allium sikkimense complex (Amaryllidaceae)

Polyploidization is a process that typically leads to instantaneous reproductive isolation and has, therefore, been considered as one of the major evolutionary forces in the species-rich Hengduan Mountains (HM), yet this topic remains poorly studied in the region. Allium sikkimense and its relatives (about eight species) compose a natural diploid-polyploid complex with the highest diversity in the HM and adjacent areas. A combination of nuclear ribosomal DNA (nrDNA), plastome, transcriptome, and ploidy identification through chromosome counting and flow cytometry is employed to reconstruct the phylogenetic relationships in this complex and to investigate the frequency and the evolutionary significance of polyploidy in the complex. The plastome failed to resolve the phylogenetic relationships of the different species in the A. sikkimense complex, and the phylogenetic tree based on nrDNA also has limited resolution. However, our study reveals a well-resolved phylogenetic framework for species in the A. sikkimense complex using more than 1,000 orthologous genes from the transcriptome data. Previously recognized morphospecies A. sikkimense are non-monophyletic and comprise at least two independently evolved lineages (i.e., cryptic species), each forming a clade with different diploid species in this complex. The embedded pattern of octoploid A. jichouense and tetraploid A. sp. nov. within different polyploid samples of A. sikkimense supports a possible scenario of budding speciation (via niche divergence). Furthermore, our results reveal that co-occurring species in the A. sikkimense complex usually have different ploidy levels, suggesting that polyploidy is an important process for reproductive isolation of sympatric Allium species. Phylogenetic network analyses suggested that the phylogenetic relationships of the A. sikkimense complex, allowing for reticulation events, always fit the dataset better than a simple bifurcating tree. In addition, the included or exserted filaments, which have long been used to delimit species, are highly unreliable taxonomically due to their extensive parallel and convergent evolution.

Genetic Structure of the Liriope muscari Polyploid Complex and the Possibility of Its Genetic Disturbance in Japan

Anthropogenic activities, such as the movement of plants through greening, can result in genetic disturbance that can interfere with local adaptation in wild populations. Although research is underway to prevent genetic disturbance associated with greening, genetic disturbance of intraspecific polyploidy, which is estimated to be present in 24% of vascular plants, has not been well studied. Liriope muscari is a polyploid complex with known diploid (2n = 36), tetraploid (2n = 72), and hexaploid (2n = 108) forms. The plants of this species tolerate dry and hot conditions and are therefore frequently used for greening and gardening. However, the distribution of this polyploid in Japan, its genetic structure, and genetic disturbance are not known. In this study, we investigated the polyploidy distribution and genetic structure in naturally distributed L. muscari in Japan using chloroplast DNA (cpDNA) haplotypes and nuclear DNA (nDNA). Commercially produced individuals were also studied and compared with natural populations to assess any genetic disturbance of the ploidy complex in this species. Chromosome counts, cpDNA, and nDNA results showed three genetically and cytologically distinct groups in Japan: first, a tetraploid group in mainland Japan; second, a hexaploid group in the Ryukyu Islands; and third, a diploid and tetraploid group in the Ryukyu Islands. Significant isolation by distance was also detected within the three groups (p = 0.001). Genetic disturbance due to greening and gardening should be avoided among the three groups. Genetic disturbance can be reduced by using individuals derived from natural populations that are close to the sites used for greening and gardening. For commercially produced individuals, genetic disturbance is unlikely in the Kanto region, an area of high usage, while genetic disturbance is thought possible in the Ryukyu Islands.

Phylogeography and Population Genetics Analyses Reveal Evolutionary History of the Desert Resource Plant Lycium ruthenicum (Solanaceae)

Climactic oscillations during the Quaternary played a significant role in the formation of genetic diversity and historical demography of numerous plant species in northwestern China. In this study, we used 11 simple sequence repeats derived from expressed sequence tag (EST-SSR), two chloroplast DNA (cpDNA) fragments, and ecological niche modeling (ENM) to investigate the population structure and the phylogeographic history of Lycium ruthenicum, a plant species adapted to the climate in northwestern China. We identified 20 chloroplast haplotypes of which two were dominant and widely distributed in almost all populations. The species has high haplotype diversity and low nucleotide diversity based on the cpDNA data. The EST-SSR results showed a high percentage of total genetic variation within populations. Both the cpDNA and EST-SSR results indicated no significant differentiation among populations. By combining the evidence from ENM and demographic analysis, we confirmed that both the last interglacial (LIG) and late-glacial maximum (LGM) climatic fluctuations, aridification might have substantially narrowed the distribution range of this desert species, the southern parts of the Junggar Basin, the Tarim Basin, and the eastern Pamir Plateau were the potential glacial refugia for L. ruthenicum during the late middle Pleistocene to late Pleistocene Period. During the early Holocene, the warm, and humid climate promoted its demographic expansion in northwestern China. This work may provide new insights into the mechanism of formation of plant diversity in this arid region.

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.

Plio-Pleistocene climatic change drives allopatric speciation and population divergence within the Scrophularia incisa complex (Scrophulariaceae) of desert and steppe subshrubs in Northwest China

Numerous temperate plants and animals on the Qinghai-Tibet Plateau (QTP) are hypothesized to have differentiated due to vicariant allopatric speciation associated with the geologic uplifts. However, this hypothesis has rarely been tested through a phylogeographic study of relative species in a broader geographic context, including the QTP, Tianshan Mountains, Mongolian Plateau, and surrounding regions. To understand the speciation and diversification process of plants across this wide area, phylogeographic analysis were examined from Scrophularia incisa and two other closely relative species comprising S. kiriloviana and S. dentata. Thirty-two populations of the three close relatives were genotyped using chloroplast DNA fragments and nuclear microsatellite loci to assess population structure and diversity, supplemented by phylogenetic dating, ancestral area reconstructions and species distribution modelings, as well as niche identity tests. Our chloroplast DNA (cpDNA) phylogeny showed that this monophyletic group of desert and steppe semi-shrub is derived from a Middle Pliocene ancestor of the Central Asia. Lineages in Central Asia vs. China diverged through climate/tectonic-induced vicariance during Middle Pliocene. Genetic and ENM data in conjunction with niche differentiation analyses support that the divergence of S. incisa, S. dentata and S. kiriloviana in China lineage proceeded through allopatric speciation, might triggered by early Pleistocene climate change of increase of aridification and enlargement of deserts, while subsequent climate-induced cycles of range contractions/expansions enhanced the geographical isolation and habit fragmentation of these taxa. These findings highlight the importance of the Plio-Pleistocene climate change in shaping genetic diversity and driving speciation in temperate steppes and deserts of Northwestern China.

Cytotype distribution and chloroplast phylogeography of the Actinidia chinensis complex

BACKGROUND

Plant phylogeographic studies of species in subtropical China have mainly focused on rare and endangered species, whereas few studies have been conducted on taxa with relatively wide distribution, especially polyploid species. We investigated the cytotype and haplotype distribution pattern of the Actinidia chinensis complex, a widespread geographically woody liana with variable ploidy in subtropical China comprising two varieties, with three chloroplast fragments DNA (ndhF-rpl132, rps16-trnQ and trnE-trnT). Macroevolutionary, microevolutionary and niche modeling tools were also combined to disentangle the origin and the demographic history of the species or cytotypes.

RESULTS

The ploidy levels of 3338 individuals from 128 populations sampled throughout the species distribution range were estimated with flow cytometry. The widespread cytotypes were diploids followed by tetraploids and hexaploids, whereas triploids and octoploids occurred in a few populations. Thirty-one chloroplast haplotypes were detected. The genetic diversity and genetic structure were found to be high between varieties (or ploidy races) chinensis and deliciosa. Our results revealed that these two varieties inhabit significantly different climatic niche spaces. Ecological niche models (ENMs) indicate that all varieties' ranges contracted during the Last Inter Glacial (LIG), and expanded eastward or northward during the Last Glacial Maximum (LGM).

CONCLUSIONS

Pliocene and Plio-Pleistocene climatic fluctuations and vicariance appear to have played key roles in shaping current population structure and historical demography in the A. chinensis complex. The polyploidization process also appears to have played an important role in the historical demography of the complex through improving their adaptability to environmental changes.

The effects of migration load, selfing, inbreeding depression, and the genetics of adaptation on autotetraploid versus diploid establishment in peripheral habitats

The distribution and abundance of polyploids has intrigued biologists since their discovery in the early 20th century. A pattern in nature that may give insight to processes that shape the distribution and abundance of polyploids is that polyploid populations are sometimes associated with peripheral habitats within the range of a species of mixed ploidy. Here, adaptation and competition of a diploid versus an autotetraploid population in a peripheral habitat are examined theoretically. It is shown that a nascent autotetraploid population adapts to and outcompetes a diploid population in the periphery when the rate of gamete dispersal is high, and when the mode of gene action is recessive for moderate to high rates of selfing. With additive or dominant modes of gene action, the conditions for an autotetraploid to outcompete a diploid in the periphery appear determined more by the rate of selfing and less by gamete dispersal. All of these results are based on empirical work that suggests inbreeding depression is higher in diploids versus autotetraploids. Generally, the results indicate that, although autotetraploids incur minority cytotype exclusion, diploids face burdens themselves. In the case of adaptation to a peripheral habitat, this burden is migration load from gamete and propagule dispersal.

Intricate Distribution Patterns of Six Cytotypes of Allium oleraceum at a Continental Scale: Niche Expansion and Innovation Followed by Niche Contraction With Increasing Ploidy Level

The establishment and success of polyploids are thought to often be facilitated by ecological niche differentiation from diploids. Unfortunately, most studies compared diploids and polyploids, ignoring variation in ploidy level in polyploids. To fill this gap, we performed a large-scale study of 11,163 samples from 1,283 populations of the polyploid perennial geophyte Allium oleraceum with reported mixed-ploidy populations, revealed distribution ranges of cytotypes, assessed their niches and explored the pattern of niche change with increasing ploidy level. Altogether, six ploidy levels (3x-8x) were identified. The most common were pentaploids (53.6%) followed by hexaploids (22.7%) and tetraploids (21.6%). Higher cytotype diversity was found at lower latitudes than at higher latitudes (>52° N), where only tetraploids and pentaploids occurred. We detected 17.4% of mixed-ploidy populations, usually as a combination of two, rarely of three, cytotypes. The majority of mixed-ploidy populations were found in zones of sympatry of the participating cytotypes, suggesting they have arisen through migration (secondary contact zone). Using coarse-grained variables (climate, soil), we found evidence of both niche expansion and innovation in tetraploids related to triploids, whereas higher ploidy levels showed almost zero niche expansion, but a trend of increased niche unfilling of tetraploids. Niche unfilling in higher ploidy levels was caused by a contraction of niche envelopes toward lower continentality of the climate and resulted in a gradual decrease of niche breadth and a gradual shift in niche optima. Field-recorded data indicated wide habitat breadth of tetraploids and pentaploids, but also a pattern of increasing synanthropy in higher ploidy levels. Wide niche breadth of tetra- and pentaploids might be related to their multiple origins from different environmental conditions, higher "age", and retained sexuality, which likely preserve their adaptive potential. In contrast, other cytotypes with narrower niches are mostly asexual, probably originating from a limited range of contrasting environments. Persistence of local ploidy mixtures could be enabled by the perenniality of A. oleraceum and its prevalence of vegetative reproduction, facilitating the establishment and decreasing exclusion of minority cytotype due to its reproductive costs. Vegetative reproduction might also significantly accelerate colonization of new areas, including recolonization of previously glaciated areas.

Insights into phylogeny, age and evolution of Allium (Amaryllidaceae) based on the whole plastome sequences

BACKGROUND AND AIMS

The genus Allium L., one of the largest monocotyledonous genera and one that includes many economically important crops with nutritional and medicinal value, has been the focus of classification or phylogeny studies for centuries. Recent studies suggested that the genus can be divided into 15 subgenera and 72 sections, which were further classified into three evolutionary lineages. However, the phylogenetic relationships reconstructed by one or two loci showed weaker support, especially for the third evolutionary lineage, which might not show the species relationships very clearly and could hinder further adaptive and evolutionary study.

METHODS

In this study, a total of 39 complete chloroplast genomes of Allium (covering 12 Allium subgenera) were collected, and combining these with 125 species of plastomes from 19 other families of monocots, we reconstructed the phylogeny of the genus Allium, estimated the origin and divergence time of the three evolutionary lineages and investigated the adaptive evolution in this genus and related families.

RESULTS

Our phylogenetic analysis confirmed the monophyly and three evolutionary lineages of Allium, while new species relationships were detected within the third evolutionary lineage. The divergence time of the three evolutionary lineages was estimated to be in the early Eocene to the middle Miocene, and numerous positive selected genes (PSGs) and PSGs with high average Ka/Ks values were found in Allium species.

CONCLUSIONS

Our results detected a well-supported phylogenetic relationship of Allium. The PSGs and PSGs with high Ka/Ks values, as well as diversified morphologies, complicated chromosome characteristics and unique reproductive modes may play important roles in the adaptation and evolution of Allium species. This is the first study that conducted phylogenetic and evolutionary analyses on the genus Allium combined with the plastome and morphological and cytological data. We hope that this study can contribute to further analysis of Allium for other researchers.

Molecular Phylogeography and Evolutionary History of the Endemic Species Corydalis hendersonii (Papaveraceae) on the Tibetan Plateau Inferred From Chloroplast DNA and ITS Sequence Variation

In response to past climatic changes, the species with different habits or adaptive traits likely have experienced very different evolutionary histories, especially for species that restricted to high mountain areas. In order to trace how Quaternary climatic oscillations affected range distributions and intraspecific divergence of such alpine plants on the Tibetan Plateau, here, we investigated maternally inherited chloroplast DNA (cpDNA) markers and biparentally inherited nuclear ribosomal internal transcribed spacer (ITS) DNA variations and aimed to explore the phylogeographic history of the endemic alpine species Corydalis hendersonii Hemsl. (Papaveraceae). We sequenced four cpDNA fragments (trnS-trnG, trnT-trnL, atpH-atpI, and psbE-petL) and also the nuclear (ITS) region in 368 individuals from 30 populations across the species' range. The network and phylogenetic analysis based on cpDNA variations identified 15 chlorotypes that cluster into three distinct clades. However, our nuclear DNA results demonstrated that there were four genetic/geographical groups within C. hendersonii. Some common and highly divergent cpDNA and ITS haplotypes were distributed in the populations of central and northeastern Tibetan Plateau, and the highest nucleotide diversity and genetic differentiation were detected in the central region. Demographic tests further indicated that the populations of southwestern and western Tibet may have experienced recent range expansion, which most likely occurred during the last glacial maximum (LGM) and continued its expansion after the beginning of the Holocene. These two different groups of this species may have derived from potential refugia that existed in the central and/or northeastern regions of Tibet during recent interglacial periods. In addition, our AMOVA analyses detected high genetic differentiation along with the whole sampling range. Also, distinct phylogeographic structures were detected among populations of C. hendersonii based on both of cpDNA and ITS variation. These findings shed new light on the importance of climatic oscillations during Quaternary and complex local topography as causes of intraspecific diversification and demographic changes within cold-tolerant herbs in the Tibetan Plateau biodiversity hotspot.

Polyploidy promotes species diversification of Allium through ecological shifts

Despite the role of polyploidy in multiple evolutionary processes, its impact on plant diversification remains controversial. An increased polyploid frequency may facilitate speciation through shifts in ecology, morphology or both. Here we used Allium to evaluate: (1) the relationship between intraspecific polyploid frequency and species diversification rate; and (2) whether this process is associated with habitat and/or trait shifts. Using eight plastid and nuclear ribosomal markers, we built a phylogeny of 448 Allium species, representing 46% of the total. We quantified intraspecific ploidy diversity, heterogeneity in diversification rates and their relationship along the phylogeny using trait-dependent diversification models. Finally, we evaluated the association between polyploidisation and habitat or trait shifts. We detected high ploidy diversity in Allium and a polyploidy-related diversification rate shift with a probability of 95% in East Asia. Allium lineages with high polyploid frequencies had higher species diversification rates than those of diploids or lineages with lower polyploid frequencies. Shifts in speciation rates were strongly correlated with habitat shifts linked to particular soil conditions; 81.7% of edaphic variation could be explained by polyploidisation. Our study emphasises the role of intraspecific polyploid frequency combined with ecological drivers on Allium diversification, which may explain plant radiations more generally.

Population genetic structure, migration, and polyploidy origin of a medicinal species Gynostemma pentaphyllum (Cucurbitaceae)

Gynostemma pentaphyllum, a member of family Cucurbitaceae, is a perennial creeping herb used as a traditional medicinal plant in China. In this study, six polymorphic nSSR and four EST-SSR primers were used to genotype 1,020 individuals in 72 wild populations of G. pentaphyllum. The genetic diversity and population structure were investigated, and ecological niche modeling was performed to reveal the evolution and demographic history of its natural populations. The results show that G. pentaphyllum has a low level of genetic diversity and high level of variation among populations because of pervasive asexual propagation, genetic drift, and long-term habitat isolation. Results of the Mantel test demonstrate that the genetic distance and geographic distance are significantly correlated among G. pentaphyllum natural populations. The populations can be divided into two clusters on the basis of genetic structure. Asymmetrical patterns of historical gene flow were observed among the clusters. For the contemporary, almost all the bidirectional gene flow of the related pairs was symmetrical with slight differences. Recent bottlenecks were experienced by 34.72% of the studied populations. The geographic range of G. pentaphyllum continues to expand northward and eastward from Hengduan Mountains. The present distribution of G. pentaphyllum is a consequence of its complex evolution. Polyploidy in G. pentaphyllum is inferred to be polygenetic. Finally, G. pentaphyllum is a species in need of protection, so in situ and ex situ measures should be considered in the future.

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.

Environmental variation shapes genetic variation in Bouteloua gracilis: Implications for restoration management of natural populations and cultivated varieties in the southwestern United States

With the increasing frequency of large-scale restoration efforts, the need to understand the adaptive genetic structure of natural plant populations and their relation to heavily utilized cultivars is critical. Bouteloua gracilis (blue grama) is a wind-dispersed, perennial grass consisting of several cytotypes (2n = 2×-6×) with a widespread distribution in western North America. The species is locally dominant and used regularly in restoration treatments. Using amplified fragment length polymorphism (AFLP) and cpDNA analyses, we assessed the genetic variability and adaptive genetic structure of blue grama within and among 44 sampling sites that are representative of the species' environmental and habitat diversity in the southwestern United States. Five cultivars were also included to investigate genetic diversity and differentiation in natural versus cultivated populations. Three main findings resulted from this study: (a) Ninety-four polymorphic AFLP markers distinguished two population clusters defined largely by samples on and off the Colorado Plateau; (b) substructure of samples on the Colorado Plateau was indicated by genetic divergence between boundary and interior regions, and was supported by cytotype distribution and cpDNA analysis; and (c) six AFLP markers were identified as "outliers," consistent with being under selection. These loci were significantly correlated to mean annual temperature, mean annual precipitation, precipitation of driest quarter, and precipitation of wettest quarter in natural populations, but not in cultivated samples. Marker × environment relationships were found to be largely influenced by cytotype and cultivar development. Our results demonstrate that blue grama is genetically variable, and exhibits genetic structure, which is shaped, in part, by environmental variability across the Colorado Plateau. Information from our study can be used to guide the selection of seed source populations for commercial development and long-term conservation management of B. gracilis, which could include genetic assessments of diversity and the adaptive potential of both natural and cultivated populations for wildland restoration.

Diverse range dynamics and dispersal routes of plants on the Tibetan Plateau during the late Quaternary

Phylogeographical studies have suggested that several plant species on the Tibetan Plateau (TP) underwent recolonization during the Quaternary and may have had distinct range dynamics in response to the last glacial. To further test this hypothesis and locate the possible historical dispersal routes, we selected 20 plant species from different parts of the TP and modeled their geographical distributions over four time periods using species distribution models (SDMs). Furthermore, we applied the least-cost path method together with SDMs and shared haplotypes to estimate their historical dispersal corridors. We identified three general scenarios of species distribution change during the late Quaternary: the 'contraction-expansion' scenario for species in the northeastern TP, the 'expansion-contraction' scenario for species in the southeast and the 'stable' scenario for widespread species. During the Quaternary, we identified that these species were likely to recolonize along the low-elevation valleys, huge mountain ranges and flat plateau platform (e.g. the Yarlung Zangbo Valley and the Himalaya). We inferred that Quaternary cyclic glaciations along with the various topographic and climatic conditions of the TP could have resulted in the diverse patterns of range shift and dispersal of Tibetan plant species. Finally, we believe that this study would provide valuable insights for the conservation of alpine species under future climate change.

The Quaternary evolutionary history, potential distribution dynamics, and conservation implications for a Qinghai-Tibet Plateau endemic herbaceous perennial, Anisodus tanguticus (Solanaceae)

Various hypotheses have been proposed about the Quaternary evolutionary history of plant species on the Qinghai-Tibet Plateau (QTP), yet only a handful of studies have considered both population genetics and ecological niche context. In this study, we proposed and compared climate refugia hypotheses based on the phylogeographic pattern of Anisodus tanguticus (three plastid DNA fragments and nuclear internal transcribed spacer regions from 32 populations) and present and past species distribution models (SDMs). We detected six plastid haplotypes in two well-differentiated lineages. Although all haplotypes could be found in its western (sampling) area, only haplotypes from one lineage occurred in its eastern area. Meanwhile, most genetic variations existed between populations (F ST = 0.822). The SDMs during the last glacial maximum and last interglacial periods showed range fragmentation in the western area and significant range contraction in the eastern area, respectively, in comparison with current potential distribution. This species may have undergone intraspecific divergence during the early Quaternary, which may have been caused by survival in different refugia during the earliest known glacial in the QTP, rather than geological isolation due to orogenesis events. Subsequently, climate oscillations during the Quaternary resulted in a dynamic distribution range for this species as well as the distribution pattern of its plastid haplotypes and nuclear genotypes. The interglacial periods may have had a greater effect on A. tanguticus than the glacial periods. Most importantly, neither genetic data nor SDM alone can fully reveal the climate refugia history of this species. We also discuss the conservation implications for this important Tibetan folk medicine plant in light of these findings and SDMs under future climate models. Together, our results underline the necessity to combine phylogeographic and SDM approaches in future investigations of the Quaternary evolutionary history of species in topographically complex areas, such as the QTP.

RAD-seq data point to a northern origin of the arctic-alpine genus Cassiope (Ericaceae)

Many arctic-alpine plants display a highly disjunct distribution between the Arctic/Boreal regions and the southern Asian mountains. Two main hypotheses have been proposed to explain the origin of this biogeographic pattern: (1) south-to-north migration in the late Pliocene/early Pleistocene, and (2) north-to-south migration during the Miocene. The genus Cassiope is disjunctly distributed between the Arctic/Boreal regions and the Himalayan-Hengduan Mountains (HHM) and was selected to test these hypotheses. We constructed a fossil-calibrated phylogeny of Ericaceae using two plastid regions to estimate the crown group age of Cassiope, and used sequence data from thousands of loci produced by restriction site associated DNA sequencing (RAD-seq) to reconstruct the phylogeny of Cassiope. We also performed Bayesian divergence time analysis and biogeographic analysis. The Cassiope crown group was estimated to have originated in the Miocene, which predates the onset of Northern hemisphere glaciation. All HHM species formed a clade together with one eastern Siberian species, and this clade was sister to all other Arctic/Boreal species. This topology implies a northern origin of Cassiope, which is confirmed by our biogeographic analysis. Our results thus suggest that the ancient north-to-south migration hypothesis is most consistent with the origin of Cassiope.

Revisiting the evolutionary events in Allium subgenus Cyathophora (Amaryllidaceae): Insights into the effect of the Hengduan Mountains Region (HMR) uplift and Quaternary climatic fluctuations to the environmental changes in the Qinghai-Tibet Plateau

The respective roles that the Hengduan Mountains Region (HMR) uplift around 4-3 Ma and Quaternary climatic oscillations played in causing the environmental changes in the Qinghai-Tibet Plateau (QTP) remain unknown. Here, we reconstruct the evolutionary history of two varieties of Allium cyathophorum and A. spicatum of subgenus Cyathophora, restricted to the HMR and the western QTP, respectively. Forty-five populations were surveyed for chloroplast and nuclear sequence variation to evaluate phylogenetic relationships, dates of divergence and ancestral area/inflorescence reconstructions. In addition, analyses were conducted on discernable micromorphologies, cytotypes and seed size variation. Our results indicated that two varieties of A. cyathophorum are separate species, i.e. A. farreri and A. cyathophorum, and the initial split of Cyathophora was triggered by the HMR uplift around 4-3 Ma. Subsequently, A. spicatum originated through the strengthened aridification in the western QTP induced vicariance of the ancestral populations in the HMR during the early Pleistocene. A self-sustaining allotetraploid species from A. farreri and A. cyathophorum was established during an interglacial period of penultimate glaciation of the QTP. Seed size variation also supports these by the colonization-competition tradeoff among small and large seeds. Our findings appear to suggest that the HMR uplift could have strengthened the development of the Asian monsoon regimes in this region and aridification in the western QTP, while the Quaternary climatic oscillations spurred the allopatric species' range shifts and created new open microhabitat for the alloploid species.

De Novo Transcriptome Assembly in Firmiana danxiaensis, a Tree Species Endemic to the Danxia Landform

Many Firmiana species are locally endemic, providing an interesting system for studying adaptation and speciation. Among these species, F. danxiaensis is a tree species endemic to Mount Danxia in Guangdong, China, which is an area known for presenting the Danxia landform. How F. danxiaensis could have adapted to the stressful environment of rocky cliffs covered with barren soils in the Danxia landform is still unknown. In this study, we performed de novo assembly of the transcriptome of F. danxiaensis, obtaining 47,221 unigenes with an N50 value of 987 bp. Homology analysis showed that 32,318 of the unigenes presented hits in the NCBI non-redundant database, and 31,857 exhibited significant matches with the protein database of Theobroma cacao. Gene Ontology (GO) annotation showed that hundreds of unigenes participated in responses to various stresses or nutritional starvation, which may help us to understand the adaptation of F. danxiaensis to Danxia landform. Additionally, we found 263 genes related to responses to Cd, partially explaining the high accumulation of Cd observed in Firmiana species. The EuKaryotic Orthologous Groups (KOG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations revealed many genes playing roles in the biosynthesis of secondary metabolites and environmental adaptation, which may also contribute to the survivor and success of Firmiana species in extreme environments. Based on the obtained transcriptome, we further identified a Firmiana-specific whole-genome duplication event that occurred approximately 20 Mya, which may have provided raw materials for the diversification of Firmiana species.

Contrasting responses to Pleistocene climate changes: a case study of two sister species Allium cyathophorum and A. spicata (Amaryllidaceae) distributed in the eastern and western Qinghai-Tibet Plateau

It has been hypothesized that species occurring in the eastern and the western Qinghai-Tibet Plateau (QTP) responded differently to climate changes during the Pleistocene. Here, we test this hypothesis by phylogeographic analysis of two sister species, Allium cyathophorum and A. spicata. We sequenced two chloroplast DNA (cpDNA) fragments (accD-psaI and the rpl16 intron) of 150 individuals, and the nuclear (ITS) region of 114 individuals, from 19 populations throughout the distributional ranges of these species. The divergence between the two species was dated at 779 - 714 thousand years before the present and was likely initiated by the most major glaciation in the QTP. Analysis of chlorotype diversity showed that A. spicata, the species occurring in the western QTP, contains much lower genetic diversity (0.25) than A. cyathophorum (0.93), which is distributed in the eastern QTP. Moreover, multiple independent tests suggested that the A. spicata population had expanded recently, while no such expansion was detected in A. cyathophorum, indicating a contrasting pattern of responses to Pleistocene climate changes. These findings highlight the importance of geographical topography in determining how species responded to the climate changes that took place in the QTP during the Pleistocene.

Phylogeographical structure inferred from cpDNA sequence variation of Zygophyllum xanthoxylon across north-west China

Zygophyllum xanthoxylon, a desert species, displaying a broad east-west continuous distribution pattern in arid Northwestern China, can be considered as a model species to investigate the biogeographical history of this region. We sequenced two chloroplast DNA spacers (psbK-psbI and rpl32-trnL) in 226 individuals from 31 populations to explore the phylogeographical structure. Median-joining network was constructed and analysis of AMOVA, SMOVA, neutrality tests and distribution analysis were used to examine genetic structure and potential range expansion. Using species distribution modeling, the geographical distribution of Z. xanthoxylon was modeled during the present and at the Last Glacial Maximum (LGM). Among 26 haplotypes, one was widely distributed, but most was restricted to either the eastern or western region. The populations with the highest levels of haplotype diversity were found in the Tianshan Mountains and its surroundings in the west, and the Helan Mountains and Alxa Plateau in the east. AMOVA and SAMOVA showed that over all populations, the species lacks phylogeographical structure, which is speculated to be the result of its specific biology. Neutrality tests and mismatch distribution analysis support past range expansions of the species. Comparing the current distribution to those cold and dry conditions in LGM, Z. xanthoxylon had a shrunken and more fragmented range during LGM. Based on the evidences from phylogeographical patterns, distribution of genetic variability, and paleodistribution modeling, Z. xanthoxylon is speculated most likely to have originated from the east and migrated westward via the Hexi Corridor.

Himalayan-Tibetan plateau uplift drives divergence of polyploid poppies: Meconopsis Viguier (Papaveraceae)

Meconopsis Viguier (Papaveraceae) is an iconic genus of alpine forbs that includes medicinal and ornamental species. This study extends previous phylogenetic analyses of Meconopsis, using ITS sequences representing all the major Meconopsis clades. Phenotypic traits are also analysed for all described species. Our results show that Meconopsis evolved as a ≥ octaploid clade, with considerable interior structure reflecting further changes in ploidy levels as well as phenotypic differentiation. We support the exclusion of a few species as Cathcartia or Papaver, making Meconopsis a Tibetan region clade. Based on average rates of nucleotide substitution in angiosperm herbs we estimate that the Meconopsis clade diverged from the Meconella clade of Papaver approximately 16.6 Ma. This is soon after the 'hard' collision of the Indian continent with Asia caused uplift of the Himalaya and Hengduan ranges, greatly extended the Tibetan plateau, and initiated monsoonal climates. Eight major clades within Meconopsis are well supported and these correspond closely to previously recognised subgenus groups. However, the relationship among the clades is poorly resolved, probably because they diverged rapidly ∼15-11 Ma. Two of these clades are ∼dodecaploid but appear to have originated independently. The eight clades have distinct distributions, variously associated with the Himalaya, the eastern Plateau and Hengduan ranges. Some Meconopsis species were not monophyletic, suggesting that a combination of multilocus molecular and phenotypic traits is required when defining and revising species.

Evolutionary history of two divergent Dmrt1 genes reveals two rounds of polyploidy origins in gibel carp

Polyploidy lineages, despite very rare in vertebrates, have been proposed to play significant role in speciation and evolutionary success, but the occurrence history and consequences are still largely unknown. In this study, we used the conserved Dmrt1 to analyze polyploidy occurrence and evolutionary process in polyploid gibel carp. We identified two divergent Dmrt1 genes and respectively localized the two genes on three homologous chromosomes. Subsequently, the corresponding full-length cDNAs and genomic sequences of Dmrt1 genes were also characterized from the closely related species including Carassius auratus auratus and Cyprinus carpio, and their two Dmrt1 genes were respectively localized on two homologous chromosomes. Significantly, the evolutionary relationship analyses among cDNA and genomic DNA sequences of these Dmrt1 genes revealed two rounds of polyploidy origins in the gibel carp: an early polyploidy might result in an common tetraploid ancestor of Carassius auratus gibelio, Carassius auratus auratus and Cyprinus carpio before 18.49 million years ago (Mya), and an late polyploidy might occur from evolutionary branch of Carassius auratus at around 0.51 Mya, which lead to the occurrence of the hexaploid gibel carp. Therefore, this study provides clear genetic evidence for understanding occurrence time and historical process of polyploidy in polyploid vertebrates.

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.

Should I stay or should I go: biogeographic and evolutionary history of a polyploid complex (Chrysanthemum indicum complex) in response to Pleistocene climate change in China

Quaternary climatic oscillations greatly influenced the distribution and pattern of biodiversity in the Northern Hemisphere. Here we examine how such oscillations in South East Asia may have affected the demographic and evolutionary history of a polyploid plant complex associated with semi-dry habitats. We analyzed plastid and nuclear ribosomal DNA (rDNA) internal transcribed spacer (ITS) sequence variation within the Chrysanthemum indicum complex (Asteraceae), which comprises diploid and polyploid plants distributed throughout China. In total, 368 individuals from 47 populations across the geographical range of the complex were analyzed. We show that the relatively widespread tetraploid form of C. indicum expanded its range southward in the Pleistocene, possibly during the most recent or previous glacial period when conditions became drier and forests retreated in southern China. In marked contrast, diploid and other polyploid members of the complex failed to expand their ranges at these times or have since undergone range contractions in contrast to tetraploid C. indicum. We conclude that hybridization and gene flow between taxa occurred frequently during the evolutionary history of the complex, causing considerable sharing of chlorotypes and ITS types. Nevertheless, taxa within ploidy levels could be largely distinguished according to chlorotype and/or ITS type.

Evolutionary history of Purple cone spruce (Picea purpurea) in the Qinghai-Tibet Plateau: homoploid hybrid origin and Pleistocene expansion

Hybridization and introgression can play an important role in speciation. Here, we examine their roles in the origin and evolution of Picea purpurea, a diploid spruce species occurring on the Qinghai-Tibet Plateau (QTP). Phylogenetic relationships and ecological differences between this species and its relatives, P. schrenkiana, P. likiangensis and P. wilsonii, are unclear. To clarify them, we surveyed sequence variation within and between them for 11 nuclear loci, three chloroplast (cp) and two mitochondrial (mt) DNA fragments, and examined their ecological requirements using ecological niche modelling. Initial analyses based on 11 nuclear loci rejected a close relationship between P. schrenkiana and P. purpurea. BP&P tests and ecological niche modelling indicated substantial divergence between the remaining three species and supported the species status of P. purpurea, which contained many private alleles as expected for a well-established species. Sequence variation for cpDNA and mtDNA suggested a close relationship between P. purpurea and P. wilsonii, while variation at the nuclear se1364 gene suggested P. purpurea was more closely related to P. likiangensis. Analyses of genetic divergence, Bayesian clustering and model comparison using approximate Bayesian computation (ABC) of nuclear (nr) DNA variation all supported the hypothesis that P. purpurea originated by homoploid hybrid speciation from P. wilsonii and P. likiangensis. The ABC analysis dated the origin of P. purpurea at the Pleistocene, and the estimated hybrid parameter indicated that 69% of its nuclear composition was contributed by P. likiangensis and 31% by P. wilsonii. Our results further suggested that during or immediately following its formation, P. purpurea was subject to organelle DNA introgression from P. wilsonii such that it came to possess both mtDNA and cpDNA of P. wilsonii. The estimated parameters indicated that following its origin, P. purpurea underwent an expansion during/after the largest Pleistocene glaciation recorded for the QTP.

Identifying multiple origins of polyploid taxa: a multilocus study of the hybrid cloak fern (Astrolepis integerrima; Pteridaceae)

PREMISE OF THE STUDY

Molecular studies have shown that multiple origins of polyploid taxa are the rule rather than the exception. To understand the distribution and ecology of polyploid species and the evolutionary significance of polyploidy in general, it is important to delineate these independently derived lineages as accurately as possible. Although gene flow among polyploid lineages and backcrossing to their diploid parents often confound this process, such post origin gene flow is very infrequent in asexual polyploids. In this study, we estimate the number of independent origins of the apomictic allopolyploid fern Astrolepis integerrima, a morphologically heterogeneous species most common in the southwestern United States and Mexico, with outlying populations in the southeastern United States and the Caribbean.

METHODS

Plastid DNA sequence and AFLP data were obtained from 33 A. integerrima individuals. Phylogenetic analysis of the sequence data and multidimensional clustering of the AFLP data were used to identify independently derived lineages.

KEY RESULTS

Analysis of the two datasets identified 10 genetic groups within the 33 analyzed samples. These groups suggest a minimum of 10 origins of A. integerrima in the northern portion of its range, with both putative parents functioning as maternal donors, both supplying unreduced gametes, and both contributing a significant portion of their genetic diversity to the hybrids.

CONCLUSIONS

Our results highlight the extreme cryptic genetic diversity and systematic complexity that can underlie a single polyploid taxon.

Out of the Qinghai-Tibet Plateau: evidence for the origin and dispersal of Eurasian temperate plants from a phylogeographic study of Hippophaë rhamnoides (Elaeagnaceae)

Numerous temperate plants now distributed across Eurasia are hypothesized to have originated and migrated from the Qinghai-Tibet Plateau (QTP) and adjacent regions. However, this hypothesis has never been tested through a phylogeographic analysis of a widely distributed species. Here, we use Hippophaë rhamnoides as a model to test this hypothesis. We collected 635 individuals from 63 populations of the nine subspecies of H. rhamnoides. We sequenced two maternally inherited chloroplast (cp) DNA fragments and also the bi-paternally inherited nuclear ribosomal ITS. We recovered five major clades in phylogenetic trees constructed from cpDNA and internal transcribed spacer (ITS) sequence variation. Most sampled individuals of six subspecies that are distributed in northern China, central Asia and Asia Minor/Europe, respectively, comprised monophyletic clades (or subclades) nested within those found in the QTP. Two subspecies in the QTP were paraphyletic, while the placement of another subspecies from the Mongolian Plateau differed between the ITS and cpDNA phylogenetic trees. Our phylogeographic analyses supported an 'out-of-QTP' hypothesis for H. rhamnoides followed by allopatric divergence, hybridization and introgression. These findings highlight the complexity of intraspecific evolutions and the importance of the QTP as a center of origin for many temperate plants.

Genetic divergence, range expansion and possible homoploid hybrid speciation among pine species in Northeast China

Although homoploid hybrid speciation in plants is probably more common than previously realized, there are few well-documented cases of homoploid hybrid origin in conifers. We examined genetic divergence between two currently widespread pines in Northeast China, Pinus sylvestris var. mongolica and Pinus densiflora, and also whether two narrowly distributed pines in the same region, Pinus funebris and Pinus takahasii, might have originated from the two widespread species by homoploid hybrid speciation. Our results, based on population genetic analysis of chloroplast (cp), mitochondrial (mt) DNA, and nuclear gene sequence variation, showed that the two widespread species were divergent for both cp- and mtDNA variation, and also for haplotype variation at two of eight nuclear gene loci surveyed. Our analysis further indicated that P. sylvestris var. mongolica and P. densiflora remained allopatric during the most severe Quaternary glacial period that occurred in Northeast China, but subsequently exhibited rapid range expansions. P. funebris and P. takahasii, were found to contain a mixture of chlorotypes and nuclear haplotypes that distinguish P. sylvestris var. mongolica and P. densiflora, in support of the hypothesis that they possibly originated via homoploid hybrid speciation following secondary contact and hybridization between P. sylvestris var. mongolica and P. densiflora.

Direction and extent of organelle DNA introgression between two spruce species in the Qinghai-Tibetan Plateau

A recent model has shown that, during range expansion of one species in a territory already occupied by a related species, introgression should take place preferentially from the resident species towards the invading species and genome components experiencing low rates of gene flow should introgress more readily than those experiencing high rates of gene flow. Here, we use molecular markers from two organelle genomes with contrasted rates of gene flow to test these predictions by examining genetic exchanges between two morphologically distinct spruce Picea species growing in the Qinghai-Tibetan Plateau. The haplotypes from both mitochondrial (mt) DNA and chloroplast (cp) DNA cluster into two distinct lineages that differentiate allopatric populations of the two species. By contrast, in sympatry, the species share the same haplotypes, suggesting interspecific genetic exchanges. As predicted by the neutral model, all sympatric populations of the expanding species had received their maternally inherited mtDNA from the resident species, whereas for paternally inherited cpDNA introgression is more limited and not strictly unidirectional. Our results underscore cryptic introgressions of organelle DNAs in plants and the importance of considering rates of gene flow and range shifts to predict direction and extent of interspecific genetic exchanges.

Species delimitation and biogeography of two fir species (Abies) in central China: cytoplasmic DNA variation

It remains unclear how speciation history might contribute to species-specific variation and affect species delimitation. We examined concordance between cytoplasmic genetic variation and morphological taxonomy in two fir species, Abies chensiensis and A. fargesii, with overlapping distributions in central China. Range-wide genetic variation was investigated using mitochondrial (mt) and plastid (pt) DNA sequences, which contrast in their rates of gene flow. Four mtDNA haplotypes were recovered and showed no obvious species' bias in terms of relative frequency. In contrast, a high level of ptDNA variation was recorded in both species with 3 common ptDNA haplotypes shared between them and 21 rare ptDNA haplotypes specific to one or other species. We argue that the lack of concordance between morphological and molecular variation between the two fir species most likely reflects extensive ancestral polymorphism sharing for both forms of cytoplasmic DNA variation. It is feasible that a relatively fast mutation rate for ptDNA contributed to the production of many species-specific ptDNA haplotypes, which remained rare due to insufficient time passing for their spread and fixation in either species, despite high levels of intraspecific ptDNA gene flow. Our phylogeographic analyses further suggest that polymorphisms in both organelle genomes most likely originated during and following glacial intervals preceding the last glacial maximum, when species distributions became fragmented into several refugia and then expanded in range across central China.

Bridging global and microregional scales: ploidy distribution in Pilosella echioides (Asteraceae) in central Europe

BACKGROUND AND AIMS

A detailed knowledge of cytotype distribution can provide important insights into the evolutionary history of polyploid systems. This study aims to explore the spatial distribution of different cytotypes in Pilosella echioides at various spatial scales (from the whole distributional range to the population level) and to outline possible evolutionary scenarios for the observed geographic pattern.

METHODS

DNA-ploidy levels were estimated using DAPI flow cytometry in 4410 individuals of P. echioides from 46 populations spread over the entire distribution range in central Europe. Special attention was paid to the cytotype structure in the most ploidy-diverse population in south-west Moravia.

KEY RESULTS

Five different cytotypes (2x, 3x, 4x, 5x and 6x) were found, the last being recorded for the first time. Although ploidy-uniform (di- or tetraploid) sites clearly prevailed, nearly one-quarter of the populations investigated harboured more (up to all five) cytotypes. Whereas penta- and hexaploids constituted only a minority of the samples, a striking predominance of the triploid cytotype was observed in several populations.

CONCLUSIONS

The representative sampling confirmed previous data on cytotype distribution, i.e. the spatial aggregation of mixed-ploidy populations in south-west Moravia and Lower Austria and the predominance of ploidy-uniform populations in other parts of the area investigated. Recurrent origin of polyploids from diploid progenitors via unreduced gametes and their successful establishment are considered the key factors promoting intrapopulational ploidy mixture ('primary hybrid zones'). As an alternative to the generally accepted theory of cytotype co-existence based on the development of different means of inter-ploidy reproductive isolation, it is suggested that a long-term ploidy mixture can also be maintained in free-mating populations provided that the polyploids originate with a sufficient frequency. In addition, the prevalence (or subdominance) of the triploid cytotype in several mixed-ploidy populations represents the first evidence of such a phenomenon in plant systems with exclusively sexual reproduction.