Evolution of the thermostability of actin-depolymerizing factors enhances the adaptation of pollen germination to high temperature.

Double fertilization in many flowering plants (angiosperms) often occurs during the hot summer season, but the mechanisms that enable angiosperms to adapt specifically to high temperatures are largely unknown. The actin cytoskeleton is essential for pollen germination and the polarized growth of pollen tubes, yet how this process responds to high temperatures remains unclear. Here, we reveal that the high thermal stability of 11 Arabidopsis (Arabidopsis thaliana) actin-depolymerizing factors (ADFs) is significantly different: ADFs that specifically accumulate in tip-growing cells (pollen and root hairs) exhibit high thermal stability. Through ancestral protein reconstruction, we found that subclass II ADFs (expressed specifically in pollen) have undergone a dynamic wave-like evolution of the retention, loss, and regeneration of thermostable sites. Additionally, the sites of AtADF7 with high thermal stability are conserved in ADFs specific to angiosperm pollen. Moreover, the high thermal stability of ADFs is required to regulate actin dynamics and turnover at high temperatures to promote pollen germination. Collectively, these findings suggest strategies for the adaptation of sexual reproduction to high temperature in angiosperms at the cell biology level.

Genome evolution of the psammophyte Pugionium for desert adaptation and further speciation.

Deserts exert strong selection pressures on plants, but the underlying genomic drivers of ecological adaptation and subsequent speciation remain largely unknown. Here, we generated de novo genome assemblies and conducted population genomic analyses of the psammophytic genus Pugionium (Brassicaceae). Our results indicated that this bispecific genus had undergone an allopolyploid event, and the two parental genomes were derived from two ancestral lineages with different chromosome numbers and structures. The postpolyploid expansion of gene families related to abiotic stress responses and lignin biosynthesis facilitated environmental adaptations of the genus to desert habitats. Population genomic analyses of both species further revealed their recent divergence with continuous gene flow, and the most divergent regions were found to be centered on three highly structurally reshuffled chromosomes. Genes under selection in these regions, which were mainly located in one of the two subgenomes, contributed greatly to the interspecific divergence in microhabitat adaptation.

Whole-genome resequencing and transcriptome analyses of four generation mutants to reveal spur-type and skin-color related genes in apple (Malus domestica Borkh. Cv. Red delicious).

BACKGROUND: Bud sport is a kind of somatic mutation that usually occurred in apple. 'Red Delicious' is considered to be a special plant material of bud sport, whereas the genetic basis of plant mutants is still unknown. In this study, we used whole-genome resequencing and transcriptome sequencing to identify genes related to spur-type and skin-color in the 'Red Delicious' (G0) and its four generation mutants including 'Starking Red' (G1), 'Starkrimson' (G2), 'Campbell Redchief' (G3) and 'Vallee Spur' (G4). RESULTS: The number of single nucleotide polymorphisms (SNPs), insertions and deletions (InDels) and structural variations (SVs) were decreased in four generation mutants compared to G0, and the number of unique SNPs and InDels were over 9-fold and 4-fold higher in G1 versus (vs.) G2 and G2 vs. G3, respectively. Chromosomes 2, 5, 11 and 15 carried the most SNPs, InDels and SVs, while chromosomes 1 and 6 carried the least. Meanwhile, we identified 4,356 variation genes by whole-genome resequencing and transcriptome, and obtained 13 and 16 differentially expressed genes (DEGs) related to spur-type and skin-color by gene expression levels. Among them, DELLA and 4CL7 were the potential genes that regulate the difference of spur-type and skin-color characters, respectively. CONCLUSIONS: Our study identified potential genes associated with spur-type and skin-color differences in 'Red Delicious' and its four generation mutants, which provides a theoretical foundation for the mechanism of the apple bud sport.

The C2 H2 -type zinc finger transcription factor OSIC1 positively regulates stomatal closure under osmotic stress in poplar.

Salt and drought impair plant osmotic homeostasis and greatly limit plant growth and development. Plants decrease stomatal aperture to reduce water loss and maintain osmotic homeostasis, leading to improved stress tolerance. Herein, we identified the C2 H2 transcription factor gene OSMOTIC STRESS INDUCED C2 H2 1 (OSIC1) from Populus alba var. pyramidalis to be induced by salt, drought, polyethylene glycol 6000 (PEG6000) and abscisic acid (ABA). Overexpression of OSIC1 conferred transgenic poplar more tolerance to high salinity, drought and PEG6000 treatment by reducing stomatal aperture, while its mutant generated by the CRISPR/Cas9 system showed the opposite phenotype. Furthermore, OSIC1 directly up-regulates PalCuAOζ in vitro and in vivo, encoding a copper-containing polyamine oxidase, to enhance H2 O2 accumulation in guard cells and thus modulates stomatal closure when stresses occur. Additionally, ABA-, drought- and salt-induced PalMPK3 phosphorylates OSIC1 to increase its transcriptional activity to PalCuAOζ. This regulation of OSIC1 at the transcriptional and protein levels guarantees rapid stomatal closure when poplar responds to osmotic stress. Our results revealed a novel transcriptional regulatory mechanism of H2 O2 production in guard cells mediated by the OSIC1-PalCuAOζ module. These findings deepen our understanding of how perennial woody plants, like poplar, respond to osmotic stress caused by salt and drought and provide potential targets for breeding.

A CLE-BAM-CIK signalling module controls root protophloem differentiation in Arabidopsis.

Exogenous application of CLAVATA3 (CLV3)/EMBRYO SURROUNDING REGION (CLE) peptides suppresses protophloem differentiation and leads to the consumption of the proximal root meristem. However, the exact CLE peptides and the corresponding receptor complex regulating protophloem differentiation have not yet been clarified. Through expression pattern and phylogenetic analyses, CLE25/26/45 were identified as candidate peptides. Further genetic analyses, physiological assays and specific protophloem marker observations indicated that CLE25/26/45, BARELY ANY MERISTEM1/3 (BAM1/3) and CLV3 INSENSITIVE KINASEs (CIKs) are involved in regulating protophloem differentiation. The cle25 26 45 and cik2 3 4 5 6 mutation can greatly rescue the root defects of brevis radix (brx) and octopus (ops) mutants. The protophloem differentiation and proximal root meristem consumption of clv1 bam1 3 and cik2 3 4 5 6 were insensitive to CLE25/26/45 treatments. cle25 26 45, clv1 bam1 3 and cik2 3 4 5 6 displayed similar premature protophloem. In addition, CLE25/26/45 induced the interactions between BAMs and CIKs in vivo. Furthermore, CLE25/26/45 enhanced the phosphorylation levels of CIKs, which were greatly impaired in clv1 bam1 3 mutant. Our work clarifies that the CLE25/26/45-BAM1/3-CIK2/3/4/5/6 signalling module genetically acts downstream of BRX and OPS to suppress protophloem differentiation.

Effects of biodegradable film mulching on bacterial diversity in soils.

The spread of biodegradable plastic films (BDFs) not only increases grain yield but also reduces environmental pollution from plastic film to a large extent. Soil microbes are considered to be involved in biodegradation processes. However, the study of microbe diversity in soil mulched with biodegradable plastic film remains limited. Here, we compared the diversity of microbes between soils with biodegradable film and nonbiodegradable film (NBDF) mulch. The results showed that BDFs affected total C, P and NH4+-N, especially organism C content, as well as microbe species richness (ACE; Chao1) and diversity (Simpson index; Shannon index). In terms of dominant phyla and genera, BDFs and NBDF can influence the abundance of disparate species. Furthermore, BDFs could also contribute to improving the richness of the important functional bacterial groups in soil, e.g., Pedomicrobium and Comamonas, both of which are involved in the degradation of plastic residues in soil. Finally, we found that BDFs improved the transformation of nitrogen by significantly increasing the abundances of Nitrobacter and Nitrospira. Our results highlight the impact of BDF mulch on the abundance of functional bacteria in the soil.

Ancient polymorphisms and divergence hitchhiking contribute to genomic islands of divergence within a poplar species complex.

How genome divergence eventually leads to speciation is a topic of prime evolutionary interest. Genomic islands of elevated divergence are frequently reported between diverging lineages, and their size is expected to increase with time and gene flow under the speciation-with-gene-flow model. However, such islands can also result from divergent sorting of ancient polymorphisms, recent ecological selection regardless of gene flow, and/or recurrent background selection and selective sweeps in low-recombination regions. It is challenging to disentangle these nonexclusive alternatives, but here we attempt to do this in an analysis of what drove genomic divergence between four lineages comprising a species complex of desert poplar trees. Within this complex we found that two morphologically delimited species, Populus euphratica and Populus pruinosa, were paraphyletic while the four lineages exhibited contrasting levels of gene flow and divergence times, providing a good system for testing hypotheses on the origin of divergence islands. We show that the size and number of genomic islands that distinguish lineages are not associated with either rate of recent gene flow or time of divergence. Instead, they are most likely derived from divergent sorting of ancient polymorphisms and divergence hitchhiking. We found that highly diverged genes under lineage-specific selection and putatively involved in ecological and morphological divergence occur both within and outside these islands. Our results highlight the need to incorporate demography, absolute divergence measurement, and gene flow rate to explain the formation of genomic islands and to identify potential genomic regions involved in speciation.

Evolution of RGF/GLV/CLEL Peptide Hormones and Their Roles in Land Plant Growth and Regulation.

Rooting is a key innovation during plant terrestrialization. RGFs/GLVs/CLELs are a family of secreted peptides, playing key roles in root stem cell niche maintenance and pattern formation. The origin of this peptide family is not well characterized. RGFs and their receptor genes, RGIs, were investigated comprehensively using phylogenetic and genetic analyses. We identified 203 RGF genes from 24 plant species, representing a variety of land plant lineages. We found that the RGF genes originate from land plants and expand via multiple duplication events. The lineage-specific RGF duplicates are retained due to their regulatory divergence, while a majority of RGFs experienced strong purifying selection in most land plants. Functional analysis indicated that RGFs and their receptor genes, RGIs, isolated from liverwort, tomato, and maize possess similar biological functions with their counterparts from Arabidopsis in root development. RGFs and RGIs are likely coevolved in land plants. Our studies shed light on the origin and functional conservation of this important peptide family in plant root development.

Plant Actin-Depolymerizing Factors Possess Opposing Biochemical Properties Arising from Key Amino Acid Changes throughout Evolution.

Functional divergence in paralogs is an important genetic source of evolutionary innovation. Actin-depolymerizing factors (ADFs) are among the most important actin binding proteins and are involved in generating and remodeling actin cytoskeletal architecture via their conserved F-actin severing or depolymerizing activity. In plants, ADFs coevolved with actin, but their biochemical properties are diverse. Unfortunately, the biochemical function of most plant ADFs and the potential mechanisms of their functional divergence remain unclear. Here, in vitro biochemical analyses demonstrated that all 11 ADF genes in Arabidopsis thaliana exhibit opposing biochemical properties. Subclass III ADFs evolved F-actin bundling (B-type) function from conserved F-actin depolymerizing (D-type) function, and subclass I ADFs have enhanced D-type function. By tracking historical mutation sites on ancestral proteins, several fundamental amino acid residues affecting the biochemical functions of these proteins were identified in Arabidopsis and various plants, suggesting that the biochemical divergence of ADFs has been conserved during the evolution of angiosperm plants. Importantly, N-terminal extensions on subclass III ADFs that arose from intron-sliding events are indispensable for the alteration of D-type to B-type function. We conclude that the evolution of these N-terminal extensions and several conserved mutations produced the diverse biochemical functions of plant ADFs from a putative ancestor.

Genome sequence and genetic transformation of a widely distributed and cultivated poplar.

Populus alba is widely distributed and cultivated in Europe and Asia. This species has been used for diverse studies. In this study, we assembled a de novo genome sequence of P. alba var. pyramidalis (= P. bolleana) and confirmed its high transformation efficiency and short transformation time by experiments. Through a process of hybrid genome assembly, a total of 464 M of the genome was assembled. Annotation analyses predicted 37 901 protein-coding genes. This genome is highly collinear to that of P. trichocarpa, with most genes having orthologs in the two species. We found a marked expansion of gene families related to histone and the hormone auxin but loss of disease resistance genes in P. alba if compared with the closely related P. trichocarpa. The genome sequence presented here represents a valuable resource for further molecular functional analyses of this species as a new tree model, poplar breeding practices and comparative genomic analyses across different poplars.

Interspecific delimitation and relationships among four Ostrya species based on plastomes.

BACKGROUND: The genus Ostrya (Betulaceae) contains eight species and four of them are distributed in China. However, studies based on limited informative sites of several chloroplast markers failed to resolve interspecific delimitation and relationships among the four Chinese species. In this study, we aimed to use the whole chloroplast genomes to address these two issues. RESULTS: We assembled and annotated 33 complete chloroplast genomes (plastomes) of the four Chinese species, representing 17 populations across most of their geographical distributions. Each species contained samples of several individuals that cover most of geographic distributions of the species. All plastomes are highly conserved in genome structure and gene order, with a total length of 158-159 kb and 122 genes. Phylogenetic analyses of whole plastomes, non-coding regions and protein-coding genes produced almost the same topological relationships. In contrast to the well-delimitated species boundary inferred from the nuclear ITS sequence variations, three of the four species are non-monophyletic in the plastome trees, which is consistent with previous studies based on a few chloroplast markers. CONCLUSIONS: The high incongruence between the ITS and plastome trees may suggest the widespread occurrences of hybrid introgression and incomplete lineage sorting during the divergence of these species. In addition, the plastomes with more informative sites compared with a few chloroplast markers still failed to resolve the phylogenetic relationships of the four species, and further studies involving population genomic data may be needed to better understand their evolutionary histories.

Genome-wide patterns of copy number variation in the Chinese yak genome.

BACKGROUND: Copy number variation (CNV) represents an important source of genetic divergence that can produce drastic phenotypic differences and may therefore be subject to selection during domestication and environmental adaptation. To investigate the evolutionary dynamics of CNV in the yak genome, we used a read depth approach to detect CNV based on genome resequencing data from 14 wild and 65 domestic yaks and determined CNV regions related to domestication and adaptations to high-altitude. RESULTS: We identified 2,634 CNV regions (CNVRs) comprising a total of 153 megabases (5.7 % of the yak genome) and 3,879 overlapping annotated genes. Comparison between domestic and wild yak populations identified 121 potentially selected CNVRs, harboring genes related to neuronal development, reproduction, nutrition and energy metabolism. In addition, we found 85 CNVRs that are significantly different between domestic yak living in high- and low-altitude areas, including three genes related to hypoxia response and six related to immune defense. This analysis shows that genic CNVs may play an important role in phenotypic changes during yak domestication and adaptation to life at high-altitude. CONCLUSIONS: We present the first refined CNV map for yak along with comprehensive genomic analysis of yak CNV. Our results provide new insights into the genetic basis of yak domestication and adaptation to living in a high-altitude environment, as well as a valuable genetic resource that will facilitate future CNV association studies of important traits in yak and other bovid species.

The salinity tolerant poplar database (STPD): a comprehensive database for studying tree salt-tolerant adaption and poplar genomics.

BACKGROUND: Soil salinity is a significant factor that impairs plant growth and agricultural productivity, and numerous efforts are underway to enhance salt tolerance of economically important plants. Populus species are widely cultivated for diverse uses. Especially, they grow in different habitats, from salty soil to mesophytic environment, and are therefore used as a model genus for elucidating physiological and molecular mechanisms of stress tolerance in woody plants. DESCRIPTION: The Salinity Tolerant Poplar Database (STPD) is an integrative database for salt-tolerant poplar genome biology. Currently the STPD contains Populus euphratica genome and its related genetic resources. P. euphratica, with a preference of the salty habitats, has become a valuable genetic resource for the exploitation of tolerance characteristics in trees. This database contains curated data including genomic sequence, genes and gene functional information, non-coding RNA sequences, transposable elements, simple sequence repeats and single nucleotide polymorphisms information of P. euphratica, gene expression data between P. euphratica and Populus tomentosa, and whole-genome alignments between Populus trichocarpa, P. euphratica and Salix suchowensis. The STPD provides useful searching and data mining tools, including GBrowse genome browser, BLAST servers and genome alignments viewer, which can be used to browse genome regions, identify similar sequences and visualize genome alignments. Datasets within the STPD can also be downloaded to perform local searches. CONCLUSIONS: A new Salinity Tolerant Poplar Database has been developed to assist studies of salt tolerance in trees and poplar genomics. The database will be continuously updated to incorporate new genome-wide data of related poplar species. This database will serve as an infrastructure for researches on the molecular function of genes, comparative genomics, and evolution in closely related species as well as promote advances in molecular breeding within Populus. The STPD can be accessed at http://me.lzu.edu.cn/stpd/ .

Transcriptome differences between two sister desert poplar species under salt stress.

BACKGROUND: Populus euphratica Oliv and P. pruinosa Schrenk (Salicaceae) both grow in dry desert areas with high summer temperatures. However, P. euphratica is distributed in dry deserts with deep underground water whereas P. pruinosa occurs in deserts in which there is underground water close to the surface. We therefore hypothesized that these two sister species may have evolved divergent regulatory and metabolic pathways during their interaction with different salt habitats and other stresses. To test this hypothesis, we compared transcriptomes from callus exposed to 24 h of salt stress and control callus samples from both species and identified differentially expressed genes (DEGs) and alternative splicing (AS) events that had occurred under salt stress. RESULTS: A total of 36,144 transcripts were identified and 1430 genes were found to be differentially expressed in at least one species in response to salt stress. Of these DEGs, 884 and 860 were identified in P. euphratica and P. pruinosa, respectively, while 314 DEGs were common to both species. On the basis of parametric analysis of gene set enrichment, GO enrichment in P. euphratica was found to be significantly different from that in P. pruinosa. Numerous genes involved in hormone biosynthesis, transporters and transcription factors showed clear differences between the two species in response to salt stress. We also identified 26,560 AS events which were mapped to 8380 poplar genomic loci from four libraries. GO enrichments for genes undergoing AS events in P. euphratica differed significantly from those in P. pruinosa. CONCLUSIONS: A number of salt-responsive genes in both P. euphratica and P. pruinosa were identified and candidate genes with potential roles in the salinity adaptation were proposed. Transcriptome comparisons of two sister desert poplar species under salt stress suggest that these two species may have developed different genetic pathways in order to adapt to different desert salt habitats. The DEGs that were found to be common to both species under salt stress may be especially important for future genetic improvement of cultivated poplars or other crops through transgenic approaches in order to increase tolerance of saline soil conditions.

Transcriptome dynamics of a desert poplar (Populus pruinosa) in response to continuous salinity stress.

KEY MESSAGE: Using RNA sequencing analysis, we identified 9,216 regulatory and salt-related genes with differential expression and temporal expression trends which provide a clear picture of transcriptomic dynamics in response to continuous salinity stress in a desert poplar, Populus pruinosa. Populus pruinosa Schrenk is native to the desert region of western China and extraordinarily well adapted to the local salt stress. Thus, it is an ideal model for studying plants' adaptation to salt stress, but its transcriptomic responses have not been previously characterized. Thus, we analyzed time- courses of these responses via a series of sequencings. In total, we generated 157.4 million 100 bp paired-end clean reads and identified 9,216 differentially expressed genes (DEGs) between salt-stressed calli and controls. Gene ontology classification analysis revealed that salt stress-related categories--including 'oxidation reduction', 'transcription factor activity', 'membrane' and 'ion channel activity'--were highly enriched among these DEGs. In addition, we grouped the 9,216 DEGs by their expression dynamics into four clusters, and the genes in each cluster showed enrichment for particular functional categories. We also found that most DEGs were activated within 24 h of the stress and their expression stabilized after 48 h. All these findings suggest that gene expression rapidly and coordinately changes during this species' adaptation to salt stress. In addition, the identified DEGs provide critical genetic resources for further functional analyses and indications of potential transgenic modifications for developing salt-tolerant poplars.

An alternative 3' splice site of PeuHKT1;3 improves the response to salt stress through enhancing affinity to K+ in Populus.

Alternative splicing (AS) serves as a crucial post-transcriptional regulator in plants that contributes to the resistance to salt stress. However, the underlying mechanism is largely unknown. In this research, we identified an important AS transcript in Populus euphratica, PeuHKT1:3a, generated by alternative 3' splice site splicing mode that resulted in the removal of 252 bases at the 5' end of the first exon in PeuHKT1:3. Protein sequence comparison showed that the site of AS occurred in PeuHKT1:3 is located at a crucial Ser residue within the first pore-loop domain, which leads to inefficient K+ transport in HKT I-type transporters. Expressing PeuHKT1;3a in an axt3 mutant yeast strain can effectively compensate for the lack of intracellular K+, whereas the expression of PeuHKT1;3 cannot yield the effect. Furthermore, in transgenic Arabidopsis and poplar plants, it was observed that lines expressing PeuHKT1;3a exhibited greater salt tolerance compared to those expressing the PeuHKT1;3 strain. Analysis of ion content and flux demonstrated that the transgenic PeuHKT1;3a line exhibited significantly higher K+ content compared to the PeuHKT1;3 line, while there was no significant difference in Na+ content. In conclusion, our findings revealed that AS can give rise to novel variants of HKT I-type proteins in P. euphratica with modified K+ selectivity to keep a higher K+/Na+ ratio to enhanced salt tolerance.

A chromosome-scale Rhubarb (Rheum tanguticum) genome assembly provides insights into the evolution of anthraquinone biosynthesis.

Rhubarb is the collective name for various perennial plants from the genus Rheum L. and the Polygonaceae family. They are one of the most ancient, commonly used, and important herbs in traditional Chinese medicine. Rhubarb is a major source of anthraquinones, but how they are synthesized remains largely unknown. Here, we generate a genome sequence assembly of one important medicinal rhubarb R. tanguticum at the chromosome level, with 2.76 Gb assembled into 11 chromosomes. The genome is shaped by two recent whole-genome duplication events and recent bursts of retrotransposons. Metabolic analyses show that the major anthraquinones are mainly synthesized in its roots. Transcriptomic analysis reveals a co-expression module with a high correlation to anthraquinone biosynthesis that includes key chalcone synthase genes. One CHS, four CYP450 and two BGL genes involved in secondary metabolism show significantly upregulated expression levels in roots compared with other tissues and clustered in the co-expression module, which implies that they may also act as candidate genes for anthraquinone biosynthesis. This study provides valuable insights into the genetic bases of anthraquinone biosynthesis that will facilitate improved breeding practices and agronomic properties for rhubarb in the future.

LncRNA expression analysis by comparative transcriptomics among closely related poplars and their regulatory roles in response to salt stress.

Long noncoding ribonucleic acids (lncRNAs) play crucial roles in regulating key biological processes; however, our knowledge of lncRNAs' roles in plant adaptive evolution is still limited. Here, we determined the divergence of conserved lncRNAs in closely related poplar species that were either tolerant or sensitive to salt stress by comparative transcriptome analysis. Among the 34,363 identified lncRNAs, ~3% were shared among poplar species with conserved sequences but diversified in their function, copy number, originating genomic region and expression patterns. Further cluster analysis revealed that the conserved lncRNAs showed more similar expression patterns within salt-tolerant poplars (Populus euphratica and P. pruinosa) than between salt-tolerant and salt-sensitive poplars. Among these lncRNAs, the antisense lncRNA lncERF024 was induced by salt and the differentiated expression between salt-sensitive and salt-tolerant poplars. The overexpression of lncERF024 in P. alba var. pyramidalis enhanced poplar tolerance to salt stress. Furthermore, RNA pull-down and RNA-seq analysis showed that numerous candidate genes or proteins associated with stress response and photosynthesis might be involved in salt resistance in PeulncERF024-OE poplars. Altogether, our study provided a novel insight into how the diversification of lncRNA expression contributes to plant adaptation traits and showed that lncERF024 may be involved in the regulation of both gene expression and protein function conferring salt tolerance in Populus.

Multi-omics data provide insight into the adaptation of the glasshouse plant Rheum nobile to the alpine subnival zone.

Subnival glasshouse plants provide a text-book example of high-altitude adaptation with reproductive organs enclosed in specialized semi-translucent bracts, monocarpic reproduction and continuous survival under stress. Here, we present genomic, transcriptomic and metabolomic analyses for one such plant, the Noble rhubarb (Rheum nobile). Comparative genomic analyses show that an expanded number of genes and retained genes from two recent whole-genome duplication events are both relevant to subnival adaptation of this species. Most photosynthesis genes are downregulated within bracts compared to within leaves, and indeed bracts exhibit a sharp reduction in photosynthetic pigments, indicating that the bracts no longer perform photosynthesis. Contrastingly, genes related to flavonol synthesis are upregulated, providing enhanced defense against UV irradiation damage. Additionally, anatomically abnormal mesophyll combined with the downregulation of genes related to mesophyll differentiation in bracts illustrates the innovation and specification of the glass-like bracts. We further detect substantial accumulation of antifreeze proteins (e.g. AFPs, LEAs) and various metabolites (e.g. Proline, Protective sugars, procyanidins) in over-wintering roots. These findings provide new insights into subnival adaptation and the evolution of glasshouse alpine plants.

Rapid radiation of Rheum (Polygonaceae) and parallel evolution of morphological traits.

In this study, we examined diversification history of Rheum and tested the hypothesis that morphological traits related to plant 'body-plans' evolved in parallel in this genus. We sequenced eight chloroplast DNA fragments (representing more than 8000 bps of the sequence for each species) of 34 species from the genus and 13 species from closely related genera. Phylogenetic analyses indicate that all species of Rheum form a monophyletic lineage sister to the two-species genus Oxyria, indicating that radiative diversifications have occurred in its evolutionary history. Our dating analyses further suggest that these radiations largely coincided with the extensive uplifts of the Qinghai-Tibetan Plateau (QTP). Ancestral state reconstruction and likelihood sensitivity tests strongly indicate that both decumbent and 'glasshouse-like' body-plan traits evolved in parallel in different clades. Our findings highlight the importance of the uplift of the QTP in promoting species diversification and the parallel evolution of morphological traits that are putatively adaptive during such an evolutionary history.