Transcriptomic analysis of a tertiary relict plant, extreme xerophyte Reaumuria soongorica to identify genes related to drought adaptation

Key anti-freeze genes and pathways of Lanzhou lily (Lilium davidii, var. unicolor) during the seedling stage

Temperature is one of the most important environmental factors for plant growth, as low-temperature freezing damage seriously affects the yield and distribution of plants. The Lanzhou lily (Lilium davidii, var. unicolor) is a famous ornamental plant with high ornamental value. Using an Illumina HiSeq transcriptome sequencing platform, sequencing was conducted on Lanzhou lilies exposed to two different temperature conditions: a normal temperature treatment at 20°C (A) and a cold treatment at -4°C (C). After being treated for 24 hours, a total of 5848 differentially expressed genes (DEGs) were identified, including 3478 significantly up regulated genes and 2370 significantly down regulated genes, accounting for 10.27% of the total number of DEGs. Quantitative real-time PCR (QRT-PCR) analysis showed that the expression trends of 10 randomly selected DEGs coincided with the results of high-throughput sequencing. In addition, genes responding to low-temperature stress were analyzed using the interaction regulatory network method. The anti-freeze pathway of Lanzhou lily was found to involve the photosynthetic and metabolic pathways, and the key freezing resistance genes were the OLEO3 gene, 9 CBF family genes, and C2H2 transcription factor c117817_g1 (ZFP). This lays the foundation for revealing the underlying mechanism of the molecular anti-freeze mechanism in Lanzhou lily.

RADseq-based population genomic analysis and environmental adaptation of rare and endangered recretohalophyte Reaumuria trigyna

Genetic diversity reflects the survival potential, history, and population dynamics of an organism. It underlies the adaptive potential of populations and their response to environmental change. Reaumuria trigyna is an endemic species in the Eastern Alxa and West Ordos desert regions in China. The species has been considered a good candidate to explore the unique survival strategies of plants that inhabit this area. In this study, we performed population genomic analyses based on restriction-site associated DNA sequencing to understand the genetic diversity, population genetic structure, and differentiation of the species. Analyses of 92,719 high-quality single-nucleotide polymorphisms (SNPs) indicated that overall genetic diversity of R. trigyna was low (HO = 0.249 and HE = 0.208). No significant genetic differentiation was observed among the investigated populations. However, a subtle population genetic structure was detected. We suggest that this might be explained by adaptive diversification reinforced by the geographical isolation of populations. Overall, 3513 outlier SNPs were located in 243 gene-coding sequences in the R. trigyna transcriptome. Potential sites under diversifying selection occurred in genes (e.g., AP2/EREBP, E3 ubiquitin-protein ligase, FLS, and 4CL) related to phytohormone regulation and synthesis of secondary metabolites which have roles in adaptation of species. Our genetic analyses provide scientific criteria for evaluating the evolutionary capacity of R. trigyna and the discovery of unique adaptions. Our findings extend knowledge of refugia, environmental adaption, and evolution of germplasm resources that survive in the Ordos area.

Physiological and transcriptome analyses highlight multiple pathways involved in drought stress in Medicago falcata

Medicago falcata is one of the leguminous forage crops, which grows well in arid and semiarid region. To fully investigate the mechanism of drought resistance response in M. falcata, we challenged the M. falcata plants with 30% PEG-6000, and performed physiological and transcriptome analyses. It was found that, the activities of antioxidant enzymes (eg. SOD, POD, and CAT) and soluble sugar content were all increased in the PEG-treated group, as compared to the control group. Transcriptome results showed that a total of 706 genes were differentially expressed in the PEG-treated plants in comparison with the control. Gene enrichment analyses on differentially expressed genes revealed that a number of genes in various pathway were significantly enriched, including the phenylpropanoid biosynthesis (ko00940) and glycolysis/gluconeogenesis (ko00010), indicating the involvement of these key pathways in drought response. Furthermore, the expression levels of seven differentially expressed genes were verified to be involved in drought response in M. falcata by qPCR. Taken together, these results will provide valuable information related to drought response in M. falcata and lay a foundation for molecular studies and genetic breeding of legume crops in future research.

Transcriptome profiling reveals that foliar water uptake occurs with C3 and crassulacean acid metabolism facultative photosynthesis in Tamarix ramosissima under extreme drought

Tamarix ramosissima is a typical desert plant species that is widely distributed in the desert areas of Northwest China. It plays a significant role in sand fixation and soil water conservation. In particular, how it uses water to survive in the desert plays an important role in plant growth and ecosystem function. Previous studies have revealed that T. ramosissima can alleviate drought by absorbing water from its leaves under extreme drought conditions. To date, there is no clear molecular regulation mechanism to explain foliar water uptake (FWU). In the present study, we correlated diurnal meteorological data, sap flow and photosynthetic parameters to determine the physical and biological characteristics of FWU. Our results suggested that the lesser the groundwater, the easier it is for T. ramosissima to absorb water via the leaves. Gene ontology annotation and Kyoto Encyclopaedia of Genes and Genomes pathway analysis of the transcriptome profile of plants subjected to high humidity suggested that FWU was highly correlated to carbohydrate metabolism, energy transfer, pyruvate metabolism, hormone signal transduction and plant-pathogen interaction. Interestingly, as a C₃ plant, genes such as PEPC, PPDK, MDH and RuBP, which are involved in crassulacean acid metabolism (CAM) photosynthesis, were highly upregulated and accompanied by FWU. Therefore, we proposed that in the case of sufficient water supply, C₃ photosynthesis is used in T. ramosissima, whereas in cases of extreme drought, starch is degraded to provide CO₂ for CAM photosynthesis to make full use of the water obtained via FWU and the water that was transported or stored to assimilating branches and stems. This study may provide not only an important theoretical foundation for FWU and conversion from C₃ plants to CAM plants but also for engineering improved photosynthesis in high-yield drought-tolerant plants and mitigation of climate change-driven drought.

Rapidly evolving genetic features for desert adaptations in Stipagrostis pennata

BACKGROUND

Stipagrostis pennata is distributed in the mobile and semi-mobile sand dunes which can adapt well to extreme environments such as drought and high temperature. It is a pioneer plant species with potential for stabilizing sand dunes and ecological restoration. It can settle on moving sand dunes earlier than other desert plants. It can effectively improve the stability of sand dunes and help more plants settle down and increase plant diversity. However, despite its important ecological value, the genetic resources available for this species are limited.

RESULTS

We used single-molecule real-time sequencing technology to obtain the complete full-length transcriptome of Stipagrostis pennata, including 90,204 unigenes with an average length of 2624 bp. In addition, the 5436 transcription factors identified in these unigenes are rich in stress resistance genes, such as MYB-related, C3H, bHLH, GRAS and HSF, etc., which may play a role in adapting to desert drought and strong wind stress. Intron retention events are abundant alternative splicing events. Stipagrostis pennata has experienced stronger positive selection, accelerating the fixation of advantageous variants. Thirty-eight genes, such as CPP/TSO1-like gene, have evolved rapidly and may play a role in material transportation, flowering and seed formation.

CONCLUSIONS

The present study captures the complete full-length transcriptome of Stipagrostis pennata and reveals its rapid evolution. The desert adaptation in Stipagrostis pennata is reflected in the regulation of gene expression and the adaptability of gene function. Our findings provide a wealth of knowledge for the evolutionary adaptability of desert grass species.

An updated phylogenetic and biogeographic analysis based on genome skimming data reveals convergent evolution of shrubby habit in Clematis in the Pliocene and Pleistocene

Convergent evolution, often viewed as the inevitable outcome of natural selection, has received special attention since the time of Darwin. Clematis is well known for its climbing habit, but it has some shrubby species, known as sect. Fruticella s.l. The shrubby Clematis species are distributed in the dry habitats of Central Asia and adjacent areas showing possible convergent evolution. In this study, we assembled the complete plastome and nuclear ribosomal DNA (nrDNA) sequences of 56 Clematis species, representing most sections and covering most of the shrubby species, to reconstruct their evolutionary histories. Using both maximum likelihood and Bayesian methods, the plastome and nrDNA datasets generated similar, but not identical, phylogenetic relationships, which are better resolved than in previous studies. Then, molecular dating, historical range reconstruction, and character optimization analyses were conducted based on this updated phylogenetic framework. All the morphological characters widely used for taxonomy were shown to have evolved multiple times. Molecular dating inferred that Clematis diverged from its sister in the mid Miocene, and all six major clades of Clematis originated during the late Miocene, with a species radiation during the Pliocene to Pleistocene. The results clearly showed that the shrubby habit evolved independently in four lineages of Clematis in Asia. We also revealed that the shrubby lineages have emerged since the very beginning of Pliocene. Asian monsoon variation in the Pliocene and glacial period fluctuation in the Pleistocene may be the driving forces for the origin and diversification of the shrubby Clematis in Central Asia and adjacent dry areas.

Fungal Endophytic Community and Diversity Associated with Desert Shrubs Driven by Plant Identity and Organ Differentiation in Extremely Arid Desert Ecosystem

Despite desert ecosystem being crucial to our understanding of natural geography, species evolution and global climate change, there is limited information on the dynamics of their composition and the diversity of endophytic fungi communities driven by plant identity and organ differentiation. Here, an extensive investigation of endophytic fungal microbiome in root, stem, and leaf organs associated with five xerophyte shrubs in an extremely arid desert, Northwest China, were examined. The fungal community dominated by Dothideomycetes and Pleosporales. Shrub species strongly drive the niche-based processes of endophytic fungi across the root, stem and leaf compartments. The diversity and composition of endophytic fungi in stem showed higher variability among plant species than leaf and root. The fungal communities in root libraries were more diverse and exhibited a remarkable differentiation of community composition. We further demonstrated the significant host preferences and tissue specificity of desert endophytic fungi, and unique specific taxa were also observed. The co-occurrence network revealed the coexistence of fungal endophytes in arid desert, and the root fungal network harbored the highest interspecies connectivity. Members of Pleosporales were the most common keystone species in the root fungal network. This is the first report of mycobiota in both plant species and organ differentiation in an extremely arid desert ecosystem.

Leaf size variations in a dominant desert shrub, Reaumuria soongarica, adapted to heterogeneous environments

The climate in arid Central Asia (ACA) has changed rapidly in recent decades, but the ecological consequences of this are far from clear. To predict the impacts of climate change on ecosystem functioning, greater attention should be given to the relationships between leaf functional traits and environmental heterogeneity. As a dominant constructive shrub widely distributed in ACA, Reaumuria soongarica provided us with an ideal model to understand how leaf functional traits of desert ecosystems responded to the heterogeneous environments of ACA. Here, to determine the influences of genetic and ecological factors, we characterized species-wide variations in leaf traits among 30 wild populations of R. soongarica and 16 populations grown in a common garden. We found that the leaf length, width, and leaf length to width ratio (L/W) of the northern lineage were significantly larger than those of other genetic lineages, and principal component analysis based on the in situ environmental factors distinguished the northern lineage from the other lineages studied. With increasing latitude, leaf length, width, and L/W in the wild populations increased significantly. Leaf length and L/W were negatively correlated with altitude, and first increased and then decreased with increasing mean annual temperature (MAT) and mean annual precipitation (MAP). Stepwise regression analyses further indicated that leaf length variation was mainly affected by latitude. However, leaf width was uncorrelated with altitude, MAT, or MAP. The common garden trial showed that leaf width variation among the eastern populations was caused by both local adaptation and phenotypic plasticity. Our findings suggest that R. soongarica preferentially changes leaf length to adjust leaf size to cope with environmental change. We also reveal phenotypic evidence for ecological speciation of R. soongarica. These results will help us better understand and predict the consequences of climate change for desert ecosystem functioning.

miRNA-mRNA Integrated Analysis Reveals Roles for miRNAs in a Typical Halophyte, Reaumuria soongorica, during Seed Germination under Salt Stress

MicroRNAs (miRNAs) are endogenous small RNAs that play a crucial role in plant growth, development, and environmental stress responses. Reaumuria soongorica is a typical halophyte that is widely distributed in saline–alkali desert regions. Under salt stress, R. soongorica can complete germination, a critical biological process in the life cycle of seed plants. To identify miRNAs and predict target mRNAs involved in seed germination during salt stress, nine small-RNA libraries were constructed and analyzed from R. soongorica seeds treated with various concentrations of NaCl. We also obtained transcriptome data under the same treatment conditions. Further analysis identified 88 conserved miRNAs representing 25 defined families and discovered 13 novel miRNAs from nine libraries. A co-expression analysis was performed on the same samples to identify putative miRNA–mRNA interactions that were responsive to salt stress. A comparative analysis of expression during germination under 273 (threshold) and 43 mM (optimal) NaCl treatments identified 13 differentially expressed miRNAs and 23 corresponding target mRNAs, while a comparison between 43 mM NaCl and non-salt-stress conditions uncovered one differentially expressed miRNA and one corresponding target mRNA. These results provide basic data for further study of molecular mechanisms involved in the germination of salt-stressed R. soongorica seeds, and also provide a reference for the improvement of salt tolerance during plant germination.

Crassulacean Acid Metabolism Abiotic Stress-Responsive Transcription Factors: a Potential Genetic Engineering Approach for Improving Crop Tolerance to Abiotic Stress

This perspective paper explores the utilization of abiotic stress-responsive transcription factors (TFs) from crassulacean acid metabolism (CAM) plants to improve abiotic stress tolerance in crop plants. CAM is a specialized type of photosynthetic adaptation that enhances water-use efficiency (WUE) by shifting CO2 uptake to all or part of the nighttime when evaporative water losses are minimal. Recent studies have shown that TF-based genetic engineering could be a useful approach for improving plant abiotic stress tolerance because of the role of TFs as master regulators of clusters of stress-responsive genes. Here, we explore the use of abiotic stress-responsive TFs from CAM plants to improve abiotic stress tolerance and WUE in crops by controlling the expression of gene cohorts that mediate drought-responsive adaptations. Recent research has revealed several TF families including AP2/ERF, MYB, WRKY, NAC, NF-Y, and bZIP that might regulate water-deficit stress responses and CAM in the inducible CAM plant Mesembryanthemum crystallinum under water-deficit stress-induced CAM and in the obligate CAM plant Kalanchoe fedtschenkoi. Overexpression of genes from these families in Arabidopsis thaliana can improve abiotic stress tolerance in A. thaliana in some instances. Therefore, we propose that TF-based genetic engineering with a small number of CAM abiotic stress-responsive TFs will be a promising strategy for improving abiotic stress tolerance and WUE in crop plants in a projected hotter and drier landscape in the 21st-century and beyond.

Integrative omics analysis in Pandanus odorifer (Forssk.) Kuntze reveals the role of Asparagine synthetase in salinity tolerance

Pandanus odorifer (Forssk) Kuntze grows naturally along the coastal regions and withstands salt-sprays as well as strong winds. A combination of omics approaches and enzyme activity studies was employed to comprehend the mechanistic basis of high salinity tolerance in P. odorifer. The young seedlings of P. odorifer were exposed to 1 M salt stress for up to three weeks and analyzed using RNAsequencing (RNAseq) and LC-MS. Integrative omics analysis revealed high expression of the Asparagine synthetase (AS) (EC 6.3.5.4) (8.95 fold) and remarkable levels of Asparagine (Asn) (28.5 fold). This indicated that salt stress promoted Asn accumulation in P. odorifer. To understand this further, the Asn biosynthesis pathway was traced out in P. odorifer. It was noticed that seven genes involved in Asn bisynthetic pathway namely glutamine synthetase (GS) (EC 6.3.1.2) glutamate synthase (GOGAT) (EC 1.4.1.14), aspartate kinase (EC 2.7.2.4), pyruvate kinase (EC 2.7.1.40), aspartate aminotransferase (AspAT) (EC 2.6.1.1), phosphoenolpyruvate carboxylase (PEPC) (EC 4.1.1.31) and AS were up-regulated under salt stress. AS transcripts were most abundant thereby showed its highest activity and thus were generating maximal Asn under salt stress. Also, an up-regulated Na+/H+ antiporter (NHX1) facilitated compartmentalization of Na+ into vacuoles, suggesting P. odorifer as salt accumulator species.

Transcriptome changes in the phenylpropanoid pathway in senescing leaves of Toona sinensis

Toona sinensis is a deciduous tree native to eastern and southeastern Asia that has important culinary and cultural values. To expand current knowledge of the transcriptome and functional genomics in this species, a de novo transcriptome sequence analysis of young and mature leaf tissues of T. sinensis was performed using the Illumina platform. Over 8.1 Gb of data were generated, assembled into 64,541 unigenes, and annotated with known biological functions. Proteins involved in primary metabolite biosynthesis were identified based on similarities to known proteins, including some related to biosynthesis of carbohydrates, amino acids, lipids, and energy. Analysis of unigenes differentially expressed between young and mature leaves (transcriptomic libraries 'YL' and 'ML', respectively) showed that the KEGG pathways of phenylpropanoid, naringenin, lignin, cutin, suberin, and wax biosynthesis were significantly enriched in mature leaves. These results not only expand knowledge of transcriptome characteristics for this valuable species, but also provide a useful transcriptomic dataset to accelerate the researches on its metabolic mechanisms and functional genomics. This study can also further the understanding of unique aromatic metabolism and Chinese medicinal properties of T. sinensis.

Yerba mate (Ilex paraguariensis, A. St.-Hil.) de novo transcriptome assembly based on tissue specific genomic expression profiles

Background

The most common infusion in southern Latin-American countries is prepared with dried leaves of Ilex paraguariensis A. St.-Hil., an aboriginal ancestral beverage known for its high polyphenols concentration currently consumed in > 90% of homes in Argentina, in Paraguay and Uruguay. The economy of entire provinces heavily relies on the production, collection and manufacture of Ilex paraguariensis, the fifth plant species with highest antioxidant activity. Polyphenols are associated to relevant health benefits including strong antioxidant properties. Despite its regional relevance and potential biotechnological applications, little is known about functional genomics and genetics underlying phenotypic variation of relevant traits. By generating tissue specific transcriptomic profiles, we aimed to comprehensively annotate genes in the Ilex paraguariensis phenylpropanoid pathway and to evaluate differential expression profiles.

Results

In this study we generated a reliable transcriptome assembly based on a collection of 15 RNA-Seq libraries from different tissues of Ilex paraguariensis. A total of 554 million RNA-Seq reads were assembled into 193,897 transcripts, where 24,612 annotated full-length transcripts had complete ORF. We assessed the transcriptome assembly quality, completeness and accuracy using BUSCO and TransRate; consistency was also evaluated by experimentally validating 11 predicted genes by PCR and sequencing. Functional annotation against KEGG Pathway database identified 1395 unigenes involved in biosynthesis of secondary metabolites, 531 annotated transcripts corresponded to the phenylpropanoid pathway. The top 30 differentially expressed genes among tissue revealed genes involved in photosynthesis and stress response. These significant differences were then validated by qRT-PCR.

Conclusions

Our study is the first to provide data from whole genome gene expression profiles in different Ilex paraguariensis tissues, experimentally validating in-silico predicted genes key to the phenylpropanoid (antioxidant) pathway. Our results provide essential genomic data of potential use in breeding programs for polyphenol content. Further studies are necessary to assess if the observed expression variation in the phenylpropanoid pathway annotated genes is related to variations in leaves’ polyphenol content at the population scale. These results set the current reference for Ilex paraguariensis genomic studies and provide a substantial contribution to research and biotechnological applications of phenylpropanoid secondary metabolites.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5240-6) contains supplementary material, which is available to authorized users.

Sodium-Related Adaptations to Drought: New Insights From the Xerophyte Plant Zygophyllum xanthoxylum

Understanding the unusual physiological mechanisms that enable drought tolerance in xerophytes will be of considerable benefit because of the potential to identify novel and key genetic elements for future crop improvements. These plants are interesting because they are well-adapted for life in arid zones; Zygophyllum xanthoxylum, for example, is a typical xerophytic shrub that inhabits central Asian deserts, accumulating substantial levels of sodium (Na+) in its succulent leaves while growing in soils that contain very low levels of this ion. The physiological importance of this unusual trait to drought adaptations remains poorly understood, however. Thus, 2-week-old Z. xanthoxylum plants were treated with 50 mM NaCl (Na) for 7 days in this study in order to investigate their drought tolerance, leaf osmotic potential (Ψs) related parameters, anatomical characteristics, and transpiration traits. The results demonstrated that NaCl treatment significantly enhanced both the survivability and durability of Z. xanthoxylum plants under extreme drought conditions. The bulk of the Na+ ions encapsulated in plants was overwhelmingly allocated to leaves rather than roots or stems under drought conditions; thus, compared to the control, significantly more Na+ compared to other solutes such as K+, Ca2+, Cl-, sugars, and proline accumulated in the leaves of NaCl-treated plants and led to a marked decrease (31%) in leaf Ψs. In addition, the accumulation of Na+ ions also resulted in mesophyll cell enlargement and leaf succulence, enabling the additional storage of water; Na+ ions also reduced the rate of water loss by decreasing stomatal density and down-regulating stomatal aperture size. The results of this study demonstrate that Z. xanthoxylum has evolved a notable ability to utilize Na+ ions to lower Ψs, swell its leaves, and decrease stomatal aperture sizes, in order to enable the additional uptake and storage of water and mitigate losses. These distinctive drought adaption characteristics mean that the xerophytic plant Z. xanthoxylum presents a fascinating case study for the potential identification of important and novel genetic elements that could improve crops. This report provides insights on the eco-physiological role of sodium accumulation in xerophytes adapted to extremely arid habitats.

Transcriptomes Divergence of Ricotia lunaria Between the Two Micro-Climatic Divergent Slopes at "Evolution Canyon" I, Israel

As one of the hotspot regions for sympatric speciation studies, Evolution Canyon (EC) became an ideal place for its high level of microclimatic divergence interslopes. In this study, to highlight the genetic mechanisms of sympatric speciation, phenotypic variation on flowering time and transcriptomic divergence were investigated between two ecotypes of Ricotia lunaria, which inhabit the opposite temperate and tropical slopes of EC I (Lower Nahal Oren, Mount Carmel, Israel) separated by 100 m at the bottom of the slopes. Growth chamber results showed that flowering time of the ecotype from south-facing slope population # 3 (SFS 3) was significantly 3 months ahead of the north-facing slope population # 5 (NFS 5). At the same floral development stage, transcriptome analysis showed that 1,064 unigenes were differentially expressed between the two ecotypes, which enriched in the four main pathways involved in abiotic and/or biotic stresses responses, including flavonoid biosynthesis, α-linolenic acid metabolism, plant-pathogen interaction and linoleic acid metabolism. Furthermore, based on Ka/Ks analysis, nine genes were suggested to be involved in the ecological divergence between the two ecotypes, whose homologs functioned in RNA editing, ABA signaling, photoprotective response, chloroplasts protein-conducting channel, and carbohydrate metabolism in Arabidopsis thaliana. Among them, four genes, namely, SPDS1, FCLY, Tic21 and BGLU25, also showed adaptive divergence between R. lunaria and A. thaliana, suggesting that these genes could play an important role in plant speciation, at least in Brassicaceae. Based on results of both the phenotype of flowering time and comparative transcriptome, we hypothesize that, after long-time local adaptations to their interslope microclimatic environments, the molecular functions of these nine genes could have been diverged between the two ecotypes. They might differentially regulate the expression of the downstream genes and pathways that are involved in the interslope abiotic stresses, which could further diverge the flowering time between the two ecotypes, and finally induce the reproductive isolation establishment by natural selection overruling interslope gene flow, promoting sympatric speciation.

Strengthening desert plant biotechnology research in the United Arab Emirates: a viewpoint

The biotechnology of desert plants is a vast subject. The main applications in this broad field of study comprises of plant tissue culture, genetic engineering, molecular markers and others. Biotechnology applications have the potential to address biodiversity conservation as well as agricultural, medicinal, and environmental issues. There is a need to increase our knowledge of the genetic diversity through the use of molecular genetics and biotechnological approaches in desert plants in the Arabian Gulf region including those in the United Arab Emirates (UAE). This article provides a prospective research for the study of UAE desert plant diversity through DNA fingerprinting as well as understanding the mechanisms of both abiotic stress resistance (including salinity, drought and heat stresses) and biotic stress resistance (including disease and insect resistance). Special attention is given to the desert halophytes and their utilization to alleviate the salinity stress, which is one of the major challenges in agriculture. In addition, symbioses with microorganisms are thought to be hypothesized as important components of desert plant survival under stressful environmental conditions. Thus, factors shaping the diversity and functionality of plant microbiomes in desert ecosystems are also emphasized in this article. It is important to establish a critical mass for biotechnology research and applications while strengthening the channels for collaboration among research/academic institutions in the area of desert plant biotechnology.

De novo transcriptomic analysis of leaf and fruit tissue of Cornus officinalis using Illumina platform

Cornus officinalis is one of the most widely used medicinal plants in China and other East Asian countries to cure diseases such as liver, kidney, cardiovascular diseases and frequent urination for thousands of years. It is a Level 3 protected species, and is one of the 42 national key protected wild species of animals and plants in China. However, the genetics and molecular biology of C. officinalis are poorly understood, which has hindered research on the molecular mechanism of its metabolism and utilization. Hence, enriching its genomic data and information is very important. In recent years, the fast-growing technology of next generation sequencing has provided an effective path to gain genomic information from nonmodel species. This study is the first to explore the leaf and fruit tissue transcriptome of C. officinalis using the Illumina HiSeq 4000 platform. A total of 57,954,134 and 60,971,652 clean reads from leaf and fruit were acquired, respectively (GenBank number SRP115440). The pooled reads from all two libraries were assembled into 56,392 unigenes with an average length 856 bp. Among these, 41,146 unigenes matched with sequences in the NCBI nonredundant protein database. The Gene Ontology database assigned 24,336 unigenes with biological process (83.26%), cellular components (53.58%), and molecular function (83.93%). In addition, 10,808 unigenes were assigned a KOG functional classification by the KOG database. Searching against the KEGG pathway database indicated that 18,435 unigenes were mapped to 371 KEGG pathways. Moreover, the edgeR database identified 4,585 significant differentially expressed genes (DEGs), of which 1,392 were up-regulated and 3,193 were down-regulated in fruit tissue compared with leaf tissue. Finally, we explored 581 transcription factors with 50 transcription factor gene families. Most DEGs and transcription factors were related to terpene biosynthesis and secondary metabolic regulation. This study not only represented the first de novo transcriptomic analysis of C. officinalis but also provided fundamental information on its genes and biosynthetic pathway. These findings will help us explore the molecular metabolism mechanism of terpene biosynthesis in C. officinalis.

The Combined Use of Proteomics and Transcriptomics Reveals a Complex Secondary Metabolite Network in Peperomia obtusifolia

Peperomia obtusifolia, an ornamental plant from the Piperaceae family, accumulates a series of secondary metabolites with interesting biological properties. From a biosynthesis standpoint, this species produces several benzopyrans derived from orsellinic acid, which is a polyketide typically found in fungi. Additionally, the chiral benzopyrans were reported as racemic and/or as diastereomeric mixtures, which raises questions about the level of enzymatic control in the cyclization step for the formation of the 3,4-dihydro-2H-pyran moiety. Therefore, this article describes the use of shotgun proteomic and transcriptome studies as well as phytochemical profiling for the characterization of the main biosynthesis pathways active in P. obtusifolia. This combined approach resulted in the identification of a series of proteins involved in its secondary metabolism, including tocopherol cyclase and prenyltransferases. The activity of these enzymes was supported by the phytochemical profiling performed in different organs of P. obtusifolia. However, the polyketide synthases possibly involved in the production of orsellinic acid could not be identified, suggesting that orsellinic acid may be produced by endophytes intimately associated with the plant.

Expression profiling of Chrysanthemum crassum under salinity stress and the initiation of morphological changes

Chrysanthemum crassum is a decaploid species of Chrysanthemum with high stress tolerance that allows survival under salinity stress while maintaining a relatively ideal growth rate. We previously recorded morphological changes after salt treatment, such as the expansion of leaf cells. To explore the underlying salinity tolerance mechanisms, we used an Illumina platform and obtained three sequencing libraries from samples collected after 0 h, 12 h and 24 h of salt treatment. Following de novo assembly, 154,944 transcripts were generated, and 97,833 (63.14%) transcripts were annotated, including 55 Gene Ontology (GO) terms and 128 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The expression profile of C. crassum was globally altered after salt treatment. We selected functional genes and pathways that may contribute to salinity tolerance and identified some factors involved in the salinity tolerance strategies of C. crassum, such as signal transduction, transcription factors and plant hormone regulation, enhancement of energy metabolism, functional proteins and osmolyte synthesis, reactive oxygen species (ROS) scavenging, photosystem protection and recovery, and cell wall protein modifications. Forty-six genes were selected for quantitative real-time polymerase chain reaction detection, and their expression patterns were shown to be consistent with the changes in their transcript abundance determined by RNA sequencing.

Making Plants Break a Sweat: the Structure, Function, and Evolution of Plant Salt Glands

Salt stress is a complex trait that poses a grand challenge in developing new crops better adapted to saline environments. Some plants, called recretohalophytes, that have naturally evolved to secrete excess salts through salt glands, offer an underexplored genetic resource for examining how plant development, anatomy, and physiology integrate to prevent excess salt from building up to toxic levels in plant tissue. In this review we examine the structure and evolution of salt glands, salt gland-specific gene expression, and the possibility that all salt glands have originated via evolutionary modifications of trichomes. Salt secretion via salt glands is found in more than 50 species in 14 angiosperm families distributed in caryophyllales, asterids, rosids, and grasses. The salt glands of these distantly related clades can be grouped into four structural classes. Although salt glands appear to have originated independently at least 12 times, they share convergently evolved features that facilitate salt compartmentalization and excretion. We review the structural diversity and evolution of salt glands, major transporters and proteins associated with salt transport and secretion in halophytes, salt gland relevant gene expression regulation, and the prospect for using new genomic and transcriptomic tools in combination with information from model organisms to better understand how salt glands contribute to salt tolerance. Finally, we consider the prospects for using this knowledge to engineer salt glands to increase salt tolerance in model species, and ultimately in crops.

Whole-transcriptome sequence analysis of differentially expressed genes in Phormium tenax under drought stress

Phormium tenax is a kind of drought resistant garden plant with its rich and colorful leaves. To clarify the molecular mechanism of drought resistance in Phormium tenax, transcriptome was sequenced by the Illumina sequencing technology under normal and drought stress, respectively. A large number of contigs, transcripts and unigenes were obtained. Among them, only 30,814 unigenes were annotated by comparing with the protein databases. A total of 4,380 genes were differentially expressed, 2,698 of which were finally annotated under drought stress. Differentially expression analysis was also performed upon drought treatment. In KEGG pathway, the mechanism of drought resistance in Phormium tenax was explained from three aspects of metabolism and signaling of hormones, osmotic adjustment and reactive oxygen species metabolism. These results are helpful to understand the drought tolerance mechanism of Phormium tenax and will provide a precious genetic resource for drought-resistant vegetation breeding and research.

De Novo Assembly and Comparative Transcriptome Analysis Provide Insight into Lysine Biosynthesis in Toona sinensis Roem

Toona sinensis Roem is a popular leafy vegetable in Chinese cuisine and is also used as a traditional Chinese medicine. In this study, leaf samples were collected from the same plant on two development stages and then used for high-throughput Illumina RNA-sequencing (RNA-Seq). 125,884 transcripts and 54,628 unigenes were obtained through de novo assembly. A total of 25,570 could be annotated with known biological functions, which indicated that the T. sinensis leaves and shoots were undergoing multiple developmental processes especially for active metabolic processes. Analysis of differentially expressed unigenes between the two libraries showed that the lysine biosynthesis was an enriched KEGG pathway, and candidate genes involved in the lysine biosynthesis pathway in T. sinensis leaves and shoots were identified. Our results provide a primary analysis of the gene expression files of T. sinensis leaf and shoot on different development stages and afford a valuable resource for genetic and genomic research on plant lysine biosynthesis.

Genome-Scale Transcriptome Analysis of the Desert Shrub Artemisia sphaerocephala

BACKGROUND

Artemisia sphaerocephala, a semi-shrub belonging to the Artemisia genus of the Compositae family, is an important pioneer plant that inhabits moving and semi-stable sand dunes in the deserts and steppes of northwest and north-central China. It is very resilient in extreme environments. Additionally, its seeds have excellent nutritional value, and the abundant lipids and polysaccharides in the seeds make this plant a potential valuable source of bio-energy. However, partly due to the scarcity of genetic information, the genetic mechanisms controlling the traits and environmental adaptation capacity of A. sphaerocephala are unknown.

RESULTS

Here, we present the first in-depth transcriptomic analysis of A. sphaerocephala. To maximize the representation of conditional transcripts, mRNA was obtained from 17 samples, including living tissues of desert-growing A. sphaerocephala, seeds germinated in the laboratory, and calli subjected to no stress (control) and high and low temperature, high and low osmotic, and salt stresses. De novo transcriptome assembly performed using an Illumina HiSeq 2500 platform resulted in the generation of 68,373 unigenes. We analyzed the key genes involved in the unsaturated fatty acid synthesis pathway and identified 26 A. sphaerocephala fad2 genes, which is the largest fad2 gene family reported to date. Furthermore, a set of genes responsible for resistance to extreme temperatures, salt, drought and a combination of stresses was identified.

CONCLUSION

The present work provides abundant genomic information for functional dissection of the important traits of A. sphaerocephala and contributes to the current understanding of molecular adaptive mechanisms of A. sphaerocephala in the desert environment. Identification of the key genes in the unsaturated fatty acid synthesis pathway could increase understanding of the biological regulatory mechanisms of fatty acid composition traits in plants and facilitate genetic manipulation of the fatty acid composition of oil crops.

Different Sets of Post-Embryonic Development Genes Are Conserved or Lost in Two Caryophyllales Species (Reaumuria soongorica and Agriophyllum squarrosum)

Reaumuria soongorica and sand rice (Agriophyllum squarrosum) belong to the clade of Caryophyllales and are widely distributed in the desert regions of north China. Both plants have evolved many specific traits and adaptation strategies to cope with recurring environmental threats. However, the genetic basis that underpins their unique traits and adaptation remains unknown. In this study, the transcriptome data of R. soongorica and sand rice were compared with three other species with previously sequenced genomes (Arabidopsis thaliana, Oryza sativa, and Beta vulgaris). Four different gene sets were identified, namely, the genes conserved in both species, those lost in both species, those conserved in R. soongorica only, and those conserved in sand rice only. Gene ontology showed that post-embryonic development genes (PEDGs) were enriched in all gene sets, and different sets of PEDGs were conserved or lost in both the R. soongorica and sand rice genomes. Expression profiles of Arabidopsis orthologs further provided some clues to the function of the species-specific conserved PEDGs. Such orthologs included LEAFY PETIOLE, which could be a candidate gene involved in the development of branch priority in sand rice.

Different Sets of Post-Embryonic Development Genes Are Conserved or Lost in Two Caryophyllales Species (Reaumuria soongorica and Agriophyllum squarrosum)

Reaumuria soongorica and sand rice (Agriophyllum squarrosum) belong to the clade of Caryophyllales and are widely distributed in the desert regions of north China. Both plants have evolved many specific traits and adaptation strategies to cope with recurring environmental threats. However, the genetic basis that underpins their unique traits and adaptation remains unknown. In this study, the transcriptome data of R. soongorica and sand rice were compared with three other species with previously sequenced genomes (Arabidopsis thaliana, Oryza sativa, and Beta vulgaris). Four different gene sets were identified, namely, the genes conserved in both species, those lost in both species, those conserved in R. soongorica only, and those conserved in sand rice only. Gene ontology showed that post-embryonic development genes (PEDGs) were enriched in all gene sets, and different sets of PEDGs were conserved or lost in both the R. soongorica and sand rice genomes. Expression profiles of Arabidopsis orthologs further provided some clues to the function of the species-specific conserved PEDGs. Such orthologs included LEAFY PETIOLE, which could be a candidate gene involved in the development of branch priority in sand rice.

Transcriptome Analysis and Development of SSR Molecular Markers in Glycyrrhiza uralensis Fisch

Licorice is an important traditional Chinese medicine with clinical and industrial applications. Genetic resources of licorice are insufficient for analysis of molecular biology and genetic functions; as such, transcriptome sequencing must be conducted for functional characterization and development of molecular markers. In this study, transcriptome sequencing on the Illumina HiSeq 2500 sequencing platform generated a total of 5.41 Gb clean data. De novo assembly yielded a total of 46,641 unigenes. Comparison analysis using BLAST showed that the annotations of 29,614 unigenes were conserved. Further study revealed 773 genes related to biosynthesis of secondary metabolites of licorice, 40 genes involved in biosynthesis of the terpenoid backbone, and 16 genes associated with biosynthesis of glycyrrhizic acid. Analysis of unigenes larger than 1 Kb with a length of 11,702 nt presented 7,032 simple sequence repeats (SSR). Sixty-four of 69 randomly designed and synthesized SSR pairs were successfully amplified, 33 pairs of primers were polymorphism in in Glycyrrhiza uralensis Fisch., Glycyrrhiza inflata Bat., Glycyrrhiza glabra L. and Glycyrrhiza pallidiflora Maxim. This study not only presents the molecular biology data of licorice but also provides a basis for genetic diversity research and molecular marker-assisted breeding of licorice.

The role of East Asian monsoon system in shaping population divergence and dynamics of a constructive desert shrub Reaumuria soongarica

Both of the uplift of Qinghai-Tibet Plateau (QTP) and the development of East Asian monsoon system (EAMS) could have comprehensively impacted the formation and evolution of Arid Central Asia (ACA). To understand how desert plants endemic to ACA responded to these two factors, we profiled the historical population dynamics and distribution range shift of a constructive desert shrub Reaumuria soongarica (Tamaricaceae) based on species wide investigation of sequence variation of chloroplast DNA and nuclear ribosomal ITS. Phylogenetic analysis uncovered a deep divergence occurring at ca. 2.96 Mya between the western and eastern lineages of R. soongarica, and ecological niche modeling analysis strongly supported that the monsoonal climate could have fragmented its habitats in both glacial and interglacial periods and impelled its intraspecific divergence. Additionally, the population from the east monsoonal zone expanded rapidly, suggesting that the local monsoonal climate significantly impacted its population dynamics. The isolation by distance tests supported strong maternal gene flow along the direction of the East Asian winter monsoon, whose intensification induced the genetic admixture along the latitudinal populations of R. soongarica. Our results presented a new case that the development of EAMS had prominently impacted the intraspecific divergence and population dynamics of this desert plant.

Transcriptome sequencing, annotation and polymorphism detection in the hop bush, Dodonaea viscosa

Background

The hop bush, Dodonaea viscosa, is a trans-oceanic species distributed oversix continents. It evolved in Australia where it is found over a wide range of habitat types and is an ecologically important species. Limited genomic resources are currently available for this species, thus our understanding of its evolutionary history and ecological adaptation is restricted. Here, we present a comprehensive transcriptome dataset for future genomic studies into this species.

Methods

We performed Illumina sequencing of cDNA prepared from leaf tissue collected from seven populations of D. viscosa ssp. angustissima and spatulata distributed along an environmental gradient in South Australia. Sequenced reads were assembled to provide a transcriptome resource. Contiguous sequences (contigs) were annotated using BLAST searches against the NCBI non-redundant database and gene ontology definitions were assigned. Single nucleotide polymorphisms were detected for the establishment of a genetic marker set. A comparison between the two subspecies was also carried out.

Results

Illumina sequencing returned 268,672,818 sequence reads, which were de novoassembled into 105,125 contigs. Contigs with significant BLAST alignments (E value < 1e^-5)numbered at 44,191, with 38,311 of these having their most significant hits to sequences from land plant species. Gene Ontology terms were assigned to 28,440 contigs and KEGG analysis identified 146 pathways that the gene products from 5,070 contigs are potentially involved in. The subspecies comparison identified 8,494 fixed SNP differences across 3,979 contiguous sequences, indicating a level of genetic differentiation between them. Across all samples, 248,235 SNPs were detected.

Conclusions

We have established a significant genomic data resource for D. viscosa,providing a comprehensive transcriptomic reference. Genetic differences among morphologically distinct subspecies were found. A wide range of putative gene regions were identified along with a large set of variable SNP markers, providing a basis for studies into the evolution and ecological adaptation of D. viscosa.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1987-1) contains supplementary material, which is available to authorized users.

Characterization of the Transcriptome of the Xerophyte Ammopiptanthus mongolicus Leaves under Drought Stress by 454 Pyrosequencing

BACKGROUND

Ammopiptanthus mongolicus (Maxim. Ex Kom.) Cheng f., an endangered ancient legume species, endemic to the Gobi desert in north-western China. As the only evergreen broadleaf shrub in this area, A. mongolicus plays an important role in the region's ecological-environmental stability. Despite the strong potential of A. mongolicus in providing new insights on drought tolerance, sequence information on the species in public databases remains scarce. To both learn about the role of gene expression in drought stress tolerance in A. mongolicus and to expand genomic resources for the species, transcriptome sequencing of stress-treated A. mongolicus plants was performed.

RESULTS

Using 454 pyrosequencing technology, 8,480 and 7,474 contigs were generated after de novo assembly of RNA sequences from leaves of untreated and drought-treated plants, respectively. After clustering using TGICL and CAP3 programs, a combined assembly of all reads produced a total of 11,357 putative unique transcripts (PUTs). Functional annotation and classification of the transcripts were conducted by aligning the 11,357 PUTs against the public protein databases and nucleotide database (Nt). Between control and drought-treated plants, 1,620 differentially expressed genes (DEGs) were identified, of which 1,106 were up-regulated and 514 were down-regulated. The differential expression of twenty candidate genes in metabolic pathways and transcription factors families related to stress-response were confirmed by quantitative real-time PCR. Representatives of several large gene families, such as WRKY and P5CS, were identified and verified in A. mongolicus for the first time.

CONCLUSIONS

The additional transcriptome resources, gene expression profiles, functional annotations, and candidate genes provide a more comprehensive understanding of the stress response pathways in xeric-adapted plant species such as A. mongolicus.

Deep sequencing-based characterization of transcriptome of trifoliate orange (Poncirus trifoliata (L.) Raf.) in response to cold stress

BACKGROUND

Trifoliate orange (Poncirus trifoliata (L.) Raf.) is extremely cold hardy after a full acclimation; however the underlying molecular mechanisms underlying this economically valuable trait remain poorly understood. In this study, global transcriptome profiles of trifoliate orange under cold conditions (4 °C) over a time course were generated by high-throughput sequencing.

RESULTS

More than 68 million high-quality reads were produced and assembled into a non-redundant data of 77,292 unigenes with an average length of 1112 bp (N50 = 1778 bp). Of these, 23,846 had significant sequence similarity to known genes and these were assigned to 61 gene ontology (GO) categories and 25 clusters of orthologous groups (COG) involved in 128 KEGG pathways. Sequences derived from cold-treated and control plants were mapped to the assembled transcriptome, resulting in the identification of 5549 differentially expressed genes (DEGs). These comprised 600 (462 up-regulated, 138 down-regulated), 2346 (1631 up-regulated, 715 down-regulated), and 5177 (2702 up-regulated, 2475 down-regulated) genes from the cold-treated samples at 6, 24 and 72 h, respectively. The accuracy of the RNA-seq derived transcript expression data was validated by analyzing the expression patterns of 17 DEGs by qPCR. Plant hormone signal transduction, plant-pathogen interaction, and secondary metabolism were the most significantly enriched GO categories amongst in the DEGs. A total of 60 transcription factors were shown to be cold responsive. In addition, a number of genes involved in the catabolism and signaling of hormones, such as abscisic acid, ethylene and gibberellin, were affected by the cold stress. Meanwhile, levels of putrescine progressively increased under cold, which was consistent with up-regulation of an arginine decarboxylase gene.

CONCLUSIONS

This dataset provides valuable information regarding the trifoliate orange transcriptome changes in response to cold stress and may help guide future identification and functional analysis of genes that are importnatn for enhancing cold hardiness.

De novo Transcriptome Assembly of Common Wild Rice (Oryza rufipogon Griff.) and Discovery of Drought-Response Genes in Root Tissue Based on Transcriptomic Data

Background

The perennial O. rufipogon (common wild rice), which is considered to be the ancestor of Asian cultivated rice species, contains many useful genetic resources, including drought resistance genes. However, few studies have identified the drought resistance and tissue-specific genes in common wild rice.

Results

In this study, transcriptome sequencing libraries were constructed, including drought-treated roots (DR) and control leaves (CL) and roots (CR). Using Illumina sequencing technology, we generated 16.75 million bases of high-quality sequence data for common wild rice and conducted de novo assembly and annotation of genes without prior genome information. These reads were assembled into 119,332 unigenes with an average length of 715 bp. A total of 88,813 distinct sequences (74.42% of unigenes) significantly matched known genes in the NCBI NT database. Differentially expressed gene (DEG) analysis showed that 3617 genes were up-regulated and 4171 genes were down-regulated in the CR library compared with the CL library. Among the DEGs, 535 genes were expressed in roots but not in shoots. A similar comparison between the DR and CR libraries showed that 1393 genes were up-regulated and 315 genes were down-regulated in the DR library compared with the CR library. Finally, 37 genes that were specifically expressed in roots were screened after comparing the DEGs identified in the above-described analyses.

Conclusion

This study provides a transcriptome sequence resource for common wild rice plants and establishes a digital gene expression profile of wild rice plants under drought conditions using the assembled transcriptome data as a reference. Several tissue-specific and drought-stress-related candidate genes were identified, representing a fully characterized transcriptome and providing a valuable resource for genetic and genomic studies in plants.

De Novo Assembly and Characterization of Fruit Transcriptome in Black Pepper (Piper nigrum)

Black pepper is one of the most popular and oldest spices in the world and valued for its pungent constituent alkaloids. Pinerine is the main bioactive compound in pepper alkaloids, which perform unique physiological functions. However, the mechanisms of piperine synthesis are poorly understood. This study is the first to describe the fruit transcriptome of black pepper by sequencing on Illumina HiSeq 2000 platform. A total of 56,281,710 raw reads were obtained and assembled. From these raw reads, 44,061 unigenes with an average length of 1,345 nt were generated. During functional annotation, 40,537 unigenes were annotated in Gene Ontology categories, Kyoto Encyclopedia of Genes and Genomes pathways, Swiss-Prot database, and Nucleotide Collection (NR/NT) database. In addition, 8,196 simple sequence repeats (SSRs) were detected. In a detailed analysis of the transcriptome, housekeeping genes for quantitative polymerase chain reaction internal control, polymorphic SSRs, and lysine/ornithine metabolism-related genes were identified. These results validated the availability of our database. Our study could provide useful data for further research on piperine synthesis in black pepper.

RNA-seq Reveals Complicated Transcriptomic Responses to Drought Stress in a Nonmodel Tropic Plant, Bombax ceiba L

High-throughput transcriptome provides an unbiased approach for understanding the genetic basis and gene functions in response to different conditions. Here we sequenced RNA-seq libraries derived from a Bombax ceiba L. system under a controlled experiment. As a known medicinal and ornamental plant, B. ceiba grows mainly in hot-dry monsoon rainforests in Southeast Asia and Australia. Due to the specific growth environment, it has evolved a unique system that enables a physiologic response to drought stress. To date, few studies have characterized the genome-wide features of drought endurance in B. ceiba. In this study, we first attempted to characterize and identify the most differentially expressed genes and associated functional pathways under drought treatment and normal condition. Using RNA-seq technology, we generated the first transcriptome of B. ceiba and identified 59 differentially expressed genes with greater than 1,000-fold changes under two conditions. The set of upregulated genes implicates interplay among various pathways: plants growth, ubiquitin-mediated proteolysis, polysaccharides hydrolyzation, oxidative phosphorylation and photosynthesis, etc. In contrast, genes associated with stem growth, cell division, fruit ripening senescence, disease resistance, and proline synthesis are repressed. Notably, key genes of high RPKM levels in drought are AUX1, JAZ, and psbS, which are known to regulate the growth of plants, the resistance against abiotic stress, and the photosynthesis process. Furthermore, 16,656 microsatellite markers and 3,071 single-nucleotide polymorphisms (SNPs) were predicted by in silico methods. The identification and functional annotation of differentially expressed genes, microsatellites, and SNPs represent a major step forward and would serve as a valuable resource for understanding the complexity underlying drought endurance and adaptation in B. ceiba.

Transcriptomic Analysis of Drought Stress Responses in Ammopiptanthus mongolicus Leaves Using the RNA-Seq Technique

Ammopiptanthus mongolicus (Maxim. Ex Kom.) Cheng f., a relic tree of the Tertiary period, plays a critical role in maintaining desert ecosystems in the Mid-Asia region. Genome-scale gene expression profiling studies will provide deep insight into the molecular mechanism underlying the drought tolerance of A. mongolicus. In the present study, we investigated the transcriptional changes induced by drought treatment in A. mongolicus leaves by establishing a comprehensive transcriptome database and then performing a Digital Gene Expression (DGE) analysis using Solexa sequencing technology. A comprehensive transcriptome database was obtained by assembling the Illumina unigenes with expressed sequence tags (EST) available publicly, and other high throughput sequencing data. To analyze the dynamic and complicated gene regulation network during PEG6000-induced drought treatment in leaves of A. mongolicus, a time-course gene expression analysis was performed using tag-based DGE technology, which identified 437, 1,247 and 802 differentially expressed transcripts in 1, 24 and 72 h drought stress libraries, respectively. GO and KEGG analyses revealed hormone signal transduction and phenylpropanoid biosynthesis were enriched during drought treatment. A batch of drought-regulated transcription factor transcripts were identified, including the subsets of HD-ZIP, bZIP, WRKY, AP2/ERF and bHLH family members, which may play roles in drought response in A. mongolicus. The sequence collection assembled in the present study represents one of the most comprehensive transcriptome databases for A. mongolicus currently. The differentially expressed transcripts identified in our study provide a good start for identifying the key genes in stress response and performing functional analysis to reveal their roles in stress adaptation in planta.

Transcriptomic profiling of the salt-stress response in the halophyte Halogeton glomeratus

BACKGROUND

Halogeton glomeratus (H. glomeratus) is an extreme halophyte that is widely distributed in arid regions, including foothills, the Gobi desert of northwest China, and the marginal loess of Central Asia. However, research on the salt-tolerant mechanisms and genes of this species are limited because of a lack of genomic sequences. In the present study, the transcriptome of H. glomeratus was analyzed using next-generation sequencing technology to identify genes involved in salt tolerance and better understand mechanisms of salt response in the halophyte H. glomeratus.

RESULTS

Illumina RNA-sequencing was performed in five sequencing libraries that were prepared from samples treated with 200 mM NaCl for 6, 12, 24, and 72 h and a control sample to investigate changes in the H. glomeratus transcriptome in response to salt stress. The de novo assembly of five transcriptomes identified 50,267 transcripts. Among these transcripts, 31,496 (62.66%) were annotated, including 44 Gene Ontology (GO) terms and 128 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Compared with transcriptomes from the control and NaCl-treated samples, there were 2,223, 5,643, 7,510 and 10,908 genes that were differentially expressed after exposure to NaCl for 6, 12, 24, and 72 h, respectively. One hundred and eighteen salt-induced genes were common to at least two stages of salt stress, and 291 up-regulated genes were common to various stages of salt stress. Numerous genes that are related to ion transport, reactive oxygen species scavenging, energy metabolism, hormone-response pathways, and responses to biotic and abiotic stress appear to play a significant role in adaptation to salinity conditions in this species. The detection of expression patterns of 18 salt-induced genes by quantitative real-time polymerase chain reaction were basically consistent with their changes in transcript abundance determined by RNA sequencing.

CONCLUSIONS

Our findings provide a genomic sequence resource for functional genetic assignments of an extreme halophyte, H. glomeratus. We believe that the transcriptome datasets will help elucidate the genetic basis of this species' response to a salt environment and develop stress-tolerant crops based on favorable wild genetic resources.

Comparative transcriptome analysis of aboveground and underground tissues of Rhodiola algida, an important ethno-medicinal herb endemic to the Qinghai-Tibetan Plateau

Transcriptome sequencing is a powerful tool for the assessment of gene expression and the identification and characterization of molecular markers in non-model organisms. Rhodiola algida L. (Crassulaceae), endemic to the Qinghai-Tibetan Plateau, has long been used in traditional Chinese medicine to prevent altitude sickness and eliminate fatigue. Illumina-based high-throughput transcriptome sequencing of aboveground and underground tissues of R. algida respectively yielded 5.40 million and 5.18 million clean reads. A total of 82,664 unigenes averaging 577 bp in length were generated from the aboveground clean reads, with 86,237 unigenes of 502-bp mean length obtained from the underground tissues. Of 55,028 unigenes compared with sequences in UniProt databases, 20,413 unigenes had significant similarities with existing sequences in NR, NT, Swiss-Prot, GO, KEGG, and COG databases. Single nucleotide polymorphism (SNP) analysis identified 237,294 SNPs from 154,636 contigs of aboveground tissues and 197,540 SNPs from 144,963 underground-derived contigs. The information uncovered in this study should serve as a valuable resource for the characterization of important traits related to secondary metabolite formation and for the identification of associated molecular mechanisms.

De novo assembly of the desert tree Haloxylon ammodendron (C. A. Mey.) based on RNA-Seq data provides insight into drought response, gene discovery and marker identification

Background

Haloxylon ammodendron (C. A. Mey.) is widely distributed across a range of habitats, including gravel desert, clay desert, fixed and semi-fixed sand, and saline land in Asian and African deserts. To date, no genomic information or expressed sequence tag-simple sequence repeat (EST-SSR) marker has been reported for H. ammodendron plants.

Results

Using Illumina sequencing technology, we generated over two billion bases of high-quality sequence data on H. ammodendron and conducted de novo assembly and annotation of genes without prior genome information. These reads were assembled into 79,918 unigenes (mean length = 728 bp). Based on similarity searches comparing these unigenes with known proteins in the non-redundant (nr) protein database, 25,619 unigenes were functionally annotated with a cut-off E-value of 10^-5. In addition, DGE reads were mapped to the assembled transcriptome for gene expression analysis under drought stress. In total, 1,060 differentially expressed genes were identified. Among these genes, 356 genes were upregulated after drought treatment, and 704 genes were downregulated. We used the KEGG database to annotate these drought-induced genes; 207 unigenes were identified in the KEGG pathway annotation, and approximately 12.1% of the unigenes with known function fell into categories related to fatty acid metabolism, starch and sucrose metabolism, and nitrogen metabolism, suggesting that these pathways or processes may be involved in the drought response. Together, a total of 35 drought-inducible transcription factors were identified, including WRKY, MYB and bZIP family members.

Conclusions

Our study is the first to provide a transcriptome sequence resource for H. ammodendron plants and to determine its digital gene expression profile under drought conditions using the assembled transcriptome data for reference. These data provide a valuable resource for genetic and genomic studies of desert plants under abiotic conditions.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1111) contains supplementary material, which is available to authorized users.

RNA-seq analysis of Quercus pubescens Leaves: de novo transcriptome assembly, annotation and functional markers development

Quercus pubescens Willd., a species distributed from Spain to southwest Asia, ranks high for drought tolerance among European oaks. Q. pubescens performs a role of outstanding significance in most Mediterranean forest ecosystems, but few mechanistic studies have been conducted to explore its response to environmental constrains, due to the lack of genomic resources. In our study, we performed a deep transcriptomic sequencing in Q. pubescens leaves, including de novo assembly, functional annotation and the identification of new molecular markers. Our results are a pre-requisite for undertaking molecular functional studies, and may give support in population and association genetic studies.

254,265,700 clean reads were generated by the Illumina HiSeq 2000 platform, with an average length of 98 bp. De novo assembly, using CLC Genomics, produced 96,006 contigs, having a mean length of 618 bp. Sequence similarity analyses against seven public databases (Uniprot, NR, RefSeq and KOGs at NCBI, Pfam, InterPro and KEGG) resulted in 83,065 transcripts annotated with gene descriptions, conserved protein domains, or gene ontology terms. These annotations and local BLAST allowed identify genes specifically associated with mechanisms of drought avoidance. Finally, 14,202 microsatellite markers and 18,425 single nucleotide polymorphisms (SNPs) were, in silico, discovered in assembled and annotated sequences. We completed a successful global analysis of the Q. pubescens leaf transcriptome using RNA-seq. The assembled and annotated sequences together with newly discovered molecular markers provide genomic information for functional genomic studies in Q. pubescens, with special emphasis to response mechanisms to severe constrain of the Mediterranean climate. Our tools enable comparative genomics studies on other Quercus species taking advantage of large intra-specific ecophysiological differences.

Transcriptomic analysis of a psammophyte food crop, sand rice (Agriophyllum squarrosum) and identification of candidate genes essential for sand dune adaptation

Background

Sand rice (Agriophyllum squarrosum) is an annual desert plant adapted to mobile sand dunes in arid and semi-arid regions of Central Asia. The sand rice seeds have excellent nutrition value and have been historically consumed by local populations in the desert regions of northwest China. Sand rice is a potential food crop resilient to ongoing climate change; however, partly due to the scarcity of genetic information, this species has undergone only little agronomic modifications through classical breeding during recent years.

Results

We generated a deep transcriptomic sequencing of sand rice, which uncovers 67,741 unigenes. Phylogenetic analysis based on 221 single-copy genes showed close relationship between sand rice and the recently domesticated crop sugar beet. Transcriptomic comparisons also showed a high level of global sequence conservation between these two species. Conservation of sand rice and sugar beet orthologs assigned to response to salt stress gene ontology term suggests that sand rice is also a potential salt tolerant plant. Furthermore, sand rice is far more tolerant to high temperature. A set of genes likely relevant for resistance to heat stress, was functionally annotated according to expression levels, sequence annotation, and comparisons corresponding transcriptome profiling results in Arabidopsis.

Conclusions

The present work provides abundant genomic information for functional dissection of the important traits in sand rice. Future screening the genetic variation among different ecotypes and constructing a draft genome sequence will further facilitate agronomic trait improvement and final domestication of sand rice.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-872) contains supplementary material, which is available to authorized users.

Reference gene selection for quantitative real-time PCR normalization in Reaumuria soongorica

Despite its superiority for evaluating gene expression, real-time quantitative polymerase chain reaction (qPCR) results can be significantly biased by the use of inappropriate reference genes under different experimental conditions. Reaumuria soongorica is a dominant species of desert ecosystems in arid central Asia. Given the increasing interest in ecological engineering and potential genetic resources for arid agronomy, it is important to analyze gene function. However, systematic evaluation of stable reference genes should be performed prior to such analyses. In this study, the stabilities of 10 candidate reference genes were analyzed under 4 kinds of abiotic stresses (drought, salt, dark, and heat) within 4 accessions (HG010, HG020, XGG030, and XGG040) from 2 different habitats using 3 algorithms (geNorm, NormFinder, and BestKeeper). After validation of the ribulose-1,5-bisphosphate carboxylase/oxygenase large unite (rbcL) expression pattern, our data suggested that histone H2A (H2A) and eukaryotic initiation factor 4A-2 (EIF4A2) were the most stable reference genes, cyclophilin (CYCL) was moderate, and elongation factor 1α (EF1α) was the worst choice. This first systematic analysis for stably expressed genes will facilitate future functional analyses and deep mining of genetic resources in R. soongorica and other species of the Reaumuria genus.

Identification of differentially expressed genes in leaf of Reaumuria soongorica under PEG-induced drought stress by digital gene expression profiling

Reaumuria soongorica (Pall.) Maxim., a resurrection semi-shrub, is a typical constructive and dominant species in desert ecosystems in northwestern China. However, the gene expression characteristics of R. soongorica under drought stress have not been elucidated. Digital gene expression analysis was performed using Illumina technique to investigate differentially expressed genes (DEGs) between control and PEG-treated samples of R. soongorica. A total of 212,338 and 211,052 distinct tags were detected in the control and PEG-treated libraries, respectively. A total of 1,325 genes were identified as DEGs, 379 (28.6%) of which were up-regulated and 946 (71.4%) were down-regulated in response to drought stress. Functional annotation analysis identified numerous drought-inducible genes with various functions in response to drought stress. A number of regulatory proteins, functional proteins, and proteins induced by other stress factors in R. soongorica were identified. Alteration in the regulatory proteins (transcription factors and protein kinase) may be involved in signal transduction. Functional proteins, including flavonoid biosynthetic proteins, late embryogenesis abundant (LEA) proteins, small heat shock proteins (sHSP), and aquaporin and proline transporter may play protective roles in response to drought stress. Flavonoids, LEA proteins and sHSP function as reactive oxygen species scavenger or molecular chaperone. Aquaporin and proline transporters regulate the distribution of water and proline throughout the whole plant. The tolerance ability of R. soongorica may be gained through effective signal transduction and enhanced protection of functional proteins to reestablish cellular homeostasis. DEGs obtained in this study may provide useful insights to help further understand the drought-tolerant mechanism of R. soongorica.