Rapidly evolving genes and stress adaptation of two desert poplars, Populus euphratica and P. pruinosa

Populus euphratica GRP2 Interacts with Target mRNAs to Negatively Regulate Salt Tolerance by Interfering with Photosynthesis, Na+, and ROS Homeostasis

The transcription of glycine-rich RNA-binding protein 2 (PeGRP2) transiently increased in the roots and shoots of Populus euphratica (a salt-resistant poplar) upon initial salt exposure and tended to decrease after long-term NaCl stress (100 mM, 12 days). PeGRP2 overexpression in the hybrid Populus tremula × P. alba '717-1B4' (P. × canescens) increased its salt sensitivity, which was reflected in the plant's growth and photosynthesis. PeGRP2 contains a conserved RNA recognition motif domain at the N-terminus, and RNA affinity purification (RAP) sequencing was developed to enrich the target mRNAs that physically interacted with PeGRP2 in P. × canescens. RAP sequencing combined with RT-qPCR revealed that NaCl decreased the transcripts of PeGRP2-interacting mRNAs encoding photosynthetic proteins, antioxidative enzymes, ATPases, and Na+/H+ antiporters in this transgenic poplar. Specifically, PeGRP2 negatively affected the stability of the target mRNAs encoding the photosynthetic proteins PETC and RBCMT; antioxidant enzymes SOD[Mn], CDSP32, and CYB1-2; ATPases AHA11, ACA8, and ACA9; and the Na+/H+ antiporter NHA1. This resulted in (i) a greater reduction in Fv/Fm, YII, ETR, and Pn; (ii) less pronounced activation of antioxidative enzymes; and (iii) a reduced ability to maintain Na+ homeostasis in the transgenic poplars during long-term salt stress, leading to their lowered ability to tolerate salinity stress.

Spatio-temporal transcriptome dynamics coordinate rapid transition of core crop functions in 'lactating' pigeon

Pigeons (Columba livia) are among a select few avian species that have developed a specialized reproductive mode wherein the parents produce a 'milk' in their crop to feed newborn squabs. Nonetheless, the transcriptomic dynamics and role in the rapid transition of core crop functions during 'lactation' remain largely unexplored. Here, we generated a de novo pigeon genome assembly to construct a high resolution spatio-temporal transcriptomic landscape of the crop epithelium across the entire breeding stage. This multi-omics analysis identified a set of 'lactation'-related genes involved in lipid and protein metabolism, which contribute to the rapid functional transitions in the crop. Analysis of in situ high-throughput chromatin conformation capture (Hi-C) sequencing revealed extensive reorganization of promoter-enhancer interactions linked to the dynamic expression of these 'lactation'-related genes between stages. Moreover, their expression is spatially localized in specific epithelial layers, and can be correlated with phenotypic changes in the crop. These results illustrate the preferential de novo synthesis of 'milk' lipids and proteins in the crop, and provides candidate enhancer loci for further investigation of the regulatory elements controlling pigeon 'lactation'.

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

Long noncoding RNAs (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, approximately 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 (P. euphratica and P. pruinosa) than between salt-tolerant and salt-sensitive poplars. Among these lncRNAs, the antisense lncRNA lncERF024 was induced by salt and differentiated expression between salt-sensitive and salt-tolerant poplars. 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 novel insight into how the diversification of lncRNA expression contributes to plant adaptation traits and showed that lncERF024 may be involved in the regulation both of gene expression and protein function conferring salt tolerance in Populus.

Comparative transcriptome analyses of three Gentiana species provides signals for the molecular footprints of selection effects and the phylogenetic relationships

Species in Gentiana section Cruciata are important alpine plants with a center of diversity and speciation in Qinghai-Tibet Plateau (QTP), and some of these species are sympatrically distributed in northeastern QTP. Studies on genome features and natural selection signatures of sympatric species in section Crucata have been impeded by a lack of genomic resources. Here, we showed transcript characterizations and molecular footprints of selection effects on G. straminea, G. dahurica and G. officinalis based on the comparative transcriptome. A total of 62.97 Gb clean reads were obtained with unigene numbers per species ranging from 141,819 to 236,408 after assembly. We found that these three species had similar distribution of functional categories in different databases, and key enzyme-encoding genes involved in the iridoids biosynthesis were also obtained. The selective pressure analyses indicated that most paired orthologs between these three species were subject to negative selection, and only a low proportion of the orthologs that underwent positive selection were detected. We found that some positive selected genes were involved in "catalytic activity", "metabolic process", "response to stimulus" and "response to stress". Besides, large numbers of SSR primer pairs with transferabilities were successfully designed based on the available transcriptome datasets of three Gentiana species. The phylogenetic relationships reconstructed based on 352 single-copy nuclear genes provided a rough phylogenetic framework for this genus and confirmed the monophyly of section Cruciata. Our study not only provides insights for the natural selection effects on sympatric Gentiana species, but also enhances future genetic breeding or evolutionary studies on Qinjiao species.

How voles adapt to subterranean lifestyle: Insights from RNA-seq

Life under the earth surface is highly challenging and associated with a number of morphological, physiological and behavioral modifications. Subterranean niche protects animals from predators, fluctuations in environmental parameters, but is characterized by high levels of carbon dioxide and low levels of oxygen and implies high energy requirements associated with burrowing. Moreover, it lacks most of the sensory inputs available above ground. The current study describes results from RNA-seq analysis of four subterranean voles from subfamily Arvicolinae: Prometheomys schaposchnikowi, Ellobius lutescens, Terricola subterraneus, and Lasiopodomys mandarinus. Original RNA-seq data were obtained for eight species, for nine species, SRA data were downloaded from the NCBI SRA database. Additionally assembled transcriptomes of Mynomes ochrogaster and Cricetulus griseus were included in the analysis. We searched for the selection signatures and parallel amino acid substitutions in a total of 19 species. Even within this limited data set, we found significant changes of dN/dS ratio by free-ratio model analysis for subterranean Arvicolinae. Parallel substitutions were detected in genes RAD23B and PYCR2. These genes are associated with DNA repair processes and response to oxidative stress. Similar substitutions were discovered in the RAD23 genes for highly specialized subterranean Heterocephalus glaber and Fukomys damarensis. The most pronounced signatures of adaptive evolution related to subterranean niche within species of Arvicolinae subfamily were detected for Ellobius lutescens. Our results suggest that genomic adaptations can occur very quickly so far as the amount of selection signatures was found to be compliant with the degree of specialization to the subterranean niche and independent from the evolutionary age of the taxon. We found that the number of genomic signatures of selection does not depend on the age of the taxon, but is positively correlated with the degree of specialization to the subterranean niche.

Gene family expansions in Antarctic winged midge as a strategy for adaptation to cold environments

Parochlus steinenii is the only flying insect native to Antarctica. To elucidate the molecular mechanisms underlying its adaptation to cold environments, we conducted comparative genomic analyses of P. steinenii and closely related lineages. In an analysis of gene family evolution, 68 rapidly evolving gene families, involved in the innate immune system, unfolded protein response, DNA packaging, protein folding, and unsaturated fatty acid biosynthesis were detected. Some gene families were P. steinenii-specific and showed phylogenetic instability. Acyl-CoA delta desaturase and heat shock cognate protein 70 (Hsc70) were representative gene families, showing signatures of positive selection with multiple gene duplication events. Acyl-CoA delta desaturases may play pivotal roles in membrane fluidity, and expanded Hsc70 genes may function as chaperones or thermal sensors in cold environments. These findings suggest that multiple gene family expansions contributed to the adaptation of P. steinenii to cold environments.

Transcriptome-Wide Analysis Revealed the Potential of the High-Affinity Potassium Transporter (HKT) Gene Family in Rice Salinity Tolerance via Ion Homeostasis

The high-affinity potassium transporter (HKT) genes are key ions transporters, regulating the plant response to salt stress via sodium (Na⁺) and potassium (K⁺) homeostasis. The main goal of this research was to find and understand the HKT genes in rice and their potential biological activities in response to brassinosteroids (BRs), jasmonic acid (JA), seawater, and NaCl stress. The in silico analyses of seven OsHKT genes involved their evolutionary tree, gene structures, conserved motifs, and chemical properties, highlighting the key aspects of OsHKT genes. The Gene Ontology (GO) analysis of HKT genes revealed their roles in growth and stress responses. Promoter analysis showed that the majority of the HKT genes participate in abiotic stress responses. Tissue-specific expression analysis showed higher transcriptional activity of OsHKT genes in roots and leaves. Under NaCl, BR, and JA application, OsHKT1 was expressed differentially in roots and shoots. Similarly, the induced expression pattern of OsHKT1 was recorded in the seawater resistant (SWR) cultivar. Additionally, the Na⁺ to K⁺ ratio under different concentrations of NaCl stress has been evaluated. Our data highlighted the important role of the OsHKT gene family in regulating the JA and BR mediated rice salinity tolerance and could be useful for rice future breeding programs.

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.

Large-scale phylogenomics reveals ancient introgression in Asian Hepatica and new insights into the origin of the insular endemic Hepatica maxima

Hepatica maxima is native to Ulleungdo, which is one of the oceanic islands in Korea, and it likely originated via anagenetic speciation from the Korean mainland species H. asiatica. However, the relationships among the Asian lineages remain unresolved. Phylogenomics based on plant genomes can provide new insights into the evolutionary history of plants. We first generated plastid, mitochondrial and transcriptome sequences of the insular endemic species H. maxima. Using the genomic data for H. maxima, we obtained a phylogenomic dataset consisting of 76 plastid, 37 mitochondrial and 413 nuclear genes from Asian Hepatica and two outgroups. Coalescent- and concatenation-based methods revealed cytonuclear and organellar discordance in the lineage. The presence of gynodioecy with cytoplasmic male sterility in Asian Hepatica suggests that the discordance is correlated with potential disruption of linkage disequilibrium between the organellar genomes. Species network analyses revealed a deep history of hybridization and introgression in Asian Hepatica. We discovered that ancient and recent introgression events occurred throughout the evolutionary history of the insular endemic species H. maxima. The introgression may serve as an important source of genetic variation to facilitate adaptation to the Ulleungdo environment.

Comparative Transcriptomics Reveals Gene Families Associated with Predatory Behavior in Photuris femme fatale Fireflies

Identifying the basis of phenotypic variation is a key objective of genetics. This work has been mostly limited to model systems with a plethora of genetic manipulation and functional characterization tools. With the development of high-throughput sequencing and new computational tools, it is possible to identify candidate genes related to phenotypic variation in non-model organisms. Fireflies are excellent for studying phenotypic variation because of their diverse and well-characterized behaviors. Most adult fireflies emit a single mating flash pattern and do not eat. In contrast, adult females of many species in the genus Photuris employ multiple flash patterns and prey upon mate-seeking males of other firefly species. To investigate the genetic basis for this variation, we used comparative transcriptomics to identify positively selected genes between a predatory firefly, Photuris sp., and a non-predatory relative, Photuris frontalis, controlling for genes generally under selection in fireflies by comparing to a Photinus firefly. Nine gene families were identified under positive selection in the predatory versus non-predatory Photuris comparison, including genes involved in digestion, detoxification, vision, reproduction, and neural processes. These results generate intriguing hypotheses about the genetic basis for insect behavior and highlight the utility of comparative transcriptomic tools to investigate complex behaviors in non-model systems.

Molecular footprints of selection effects and whole genome duplication (WGD) events in three blueberry species: detected by transcriptome dataset

BACKGROUND

Both selection effects and whole genome duplication played very important roles in plant speciation and evolution, and to decipher the corresponding molecular footprint has always been a central task of geneticists. Vaccinium is species rich genus that comprised of about 450 species, and blueberry is one of the most important species of Vaccinium genus, which is gaining popularity because of high healthful value. In this article, we aimed to decipher the molecular footprints of natural selection on the single copy genes and WGD events occur in the evolutionary history of blueberry species.

RESULTS

We identified 30,143, 29,922 and 28,891 putative protein coding sequences from 45,535, 42,914 and 43,630 unigenes assembled from the leaves' transcriptome assembly of 19 rabbiteye (T1), 13 southern highbush (T2) and 22 northern highbush (T3) blueberry cultivars. A total of 17, 21 and 27 single copy orthologs were found to undergone positive selection in T1 versus T2, T1 versus T3, and T2 versus T3, respectively, and these orthologs were enriched in metabolic pathways including "Terpenoid backbone biosynthesis", "Valine, leucine and isoleucine biosynthesis", "Butanoate metabolism", "C5-Branched dibasic acid metabolism" "Pantothenate and CoA biosynthesis". We also detected significant molecular footprints of a recent (about 9.04 MYA), medium (about 43.44 MYA) and an ancient (about 116.39 MYA) WGD events that occurred in the evolutionary history of three blueberry species.

CONCLUSION

Some important functional genes revealed positive selection effect in blueberry. At least three rounds of WGD events were detected in the evolutionary history of blueberry species. Our work provides insights about the genetic mechanism of adaptive evolution in blueberry and species radiation of Vaccinium in short geological scale time.

Change in Spatial Distribution Patterns and Regeneration of Populus euphratica under Different Surface Soil Salinity Conditions

Ecological conservation and restoration have increasingly captured attention worldwide due to the degradation of ecosystems. As one of the most ecologically fragile areas, the Tarim River Basin, of Xinjiang, China, encountered serious decline of desert riparian forests. The Chinese government has implemented the "Ecological Water Conveyance Project" (EWCP) since 2000, protecting and restoring the dominant natural species of the desert riparian forests, i.e., Populus euphratica Oliv. The regenerative effect after the water conveyance was noteworthy. For the purpose of clarifying the mechanism of P. euphratica forest regeneration to find a better prescription for the ecological restoration works in the Tarim River Basin, we investigated the relationship between the distribution of P. euphratica and soil salinity. Experimentally evaluated the effects of surface soil salinity on P. euphratica seed germination and the influence of river flooding on the salinity of surface soils. The results showed that (1) P. euphratica trees mainly spread along the river channel within 2 km; with increasing vertical distance to the channel, the number of trees declined significantly; (2) where the salinity of the surface soil is high, there are less living P. euphratica trees; (3) the germination of P. euphratica seeds decreases with increased soil conductivity; when the soil conductivity was higher than 7 ms/cm, the germination of P. euphratica seeds was severely curtailed. (4) Flooding regimes were a pre-condition of P. euphratica restoration; they had profound effects on improving the germination of the seeds via ameliorating water conditions and reducing salinity. Our results point out that the most efficient ecological prescription for restoring and protecting desert riparian forests is to induce flooding twice yearly during June to August with 10- to 15-day durations each time. Such a plan (especially in the Tarim River Basin) should prioritize the protection of seedlings.

Comparative transcriptomic analysis reveals genetic divergence and domestication genes in Diospyros

BACKGROUND

Persimmon (Diospyros kaki) is the most economically cultivated species belonging to the genus Diospyros. However, little is known about the interspecific diversity and mechanism of domestication, partly due to the lack of genomic information that is available for closely related species of D. kaki (DK). Here, we performed transcriptome sequencing on nine samples, including DK, a variety of DK and seven closely related species, to evaluate the interspecific genetic divergence and to identify candidate genes involved in persimmon domestication.

RESULTS

We obtained a total of 483,421 unigenes with N50 at 1490 bp in the nine Diospyros samples and identified 2603 orthogroups that were shared among all the samples using OrthoMCL analysis. A phylogenetic tree was established based on the tandem 2603 one-to-one single copy gene alignments, showing that DK was closely related to D. kaki var. silvestris (DKV) and that it clustered with the clade of D. deyangnsis (DD) and was farthest from the D. cathayensis (DC) species. The nonsynonymous substitutions (Ka), via synonymous substitution (Ks) ratios, was directly proportional to the genetic relationship of the different species. The higher the Ka/Ks ratios, the longer the distance was. Moreover, 31 positively selected genes (PSGs) involved in carbohydrate metabolism and phenolic metabolism were identified and isolated, and nearly all PSGs except the MATE gene had a high expression in the DK or DKV species. It was hypothesized that these genes might contribute to the domestication of the DK species. Finally, we developed the expressed sequence tag-simple sequence repeat (EST-SSR) and identified 2 unique amplicons DKSSR10 and DKSSR39: the former was absent in the DC species but was present in the other species, the latter had a long amplification product in the DJ species.

CONCLUSION

This study presents the first transcriptome resources for the closely related species of persimmon and reveals interspecific genetic divergence. It is speculated that DK is derived from the hybridization of DD and DO species. Furthermore, our analysis suggests candidate PSGs that may be crucial for the adaptation, domestication, and speciation of persimmon relatives and suggests that DKSSR10 and DKSSSR39 could potentially serve as species-specific molecular markers.

Salicylic Acid Alleviated Salt Damage of Populus euphratica: A Physiological and Transcriptomic Analysis

Populus euphratica Oliv. is a model tree for studying abiotic stress, especially salt stress response. Salt stress is one of the most extensive abiotic stresses, which has an adverse effect on plant growth and development. Salicylic acid (SA) is an important signaling molecule that plays an important role in modulating the plant responses to abiotic stresses. To answer whether the endogenous SA can be induced by salt stress, and whether SA effectively alleviates the negative effects of salt on poplar growth is the main purpose of the study. To elucidate the effects of SA and salt stress on the growth of P. euphratica, we examined the morphological and physiological changes of P. euphratica under 300 mM NaCl after treatment with different concentrations of SA. A pretreatment of P. euphratica with 0.4 mM SA for 3 days effectively improved the growth status of plants under subsequent salt stress. These results indicate that appropriate concentrations of exogenous SA can effectively counteract the negative effect of salt stress on growth and development. Subsequently, transcripts involved in salt stress response via SA signaling were captured by RNA sequencing. The results indicated that numerous specific genes encoding mitogen-activated protein kinase, calcium-dependent protein kinase, and antioxidant enzymes were upregulated. Potassium transporters and Na+/H+ antiporters, which maintain K+/Na+ balance, were also upregulated after SA pretreatment. The transcriptome changes show that the ion transport and antioxidant enzymes were the early enhanced systems in response of P. euphratica to salt via SA, expanding our knowledge about SA function in salt stress defense in P. euphratica. This provides a solid foundation for future study of functional genes controlling effective components in metabolic pathways of trees.

Genome-wide identification and expression analysis of HKT transcription factor under salt stress in nine plant species

High-affinity K+ (HKT) gene family is regulated the transport of Na+ and maintain the balance between Na+ and K+ in the process of plant growth, development and response to abiotic stress. Despite this fact, systemic and comprehensive studies on HKT in multiply plants remains unknown. A total of 29 HKT genes distributed on nine species were identified. Phylogenetic tree analysis results indicated that HKT genes were divided into five homology subfamilies. Combining structural analysis with protein contains five highly conservative motifs, HKT family has similar gene structures and special gene characteristics. Finally, the expression patterns of HKT showed two different dramatic changes in different organs and tissues under different salt stress in multiply plants. This study has many implications for research into the comparative genomics analysis of HTK gene family, which revealed regulation mechanism of HKT genes, is valuable for understanding development and response to abiotic stress in plant.

Transcriptomic Analysis of Seed Germination Under Salt Stress in Two Desert Sister Species (Populus euphratica and P. pruinosa)

As a major abiotic stress, soil salinity limits seed germination and plant growth, development and production. Seed germination is highly related not only to the seedlings survival rate but also subsequent vegetative growth. Populus euphratica and P. pruinosa are closely related species that show a distinguished adaptability to salinity stress. In this study, we performed an integrative transcriptome analyses of three seed germination phases from P. euphratica and P. pruinosa under salt stress. A two-dimensional data set of this study provides a comprehensive view of the dynamic biochemical processes that underpin seed germination and salt tolerance. Our analysis identified 12831 differentially expressed genes (DEGs) for seed germination processes and 8071 DEGs for salt tolerance in the two species. Furthermore, we identified the expression profiles and main pathways in each growth phase. For seed germination, a large number of DEGs, including those involved in energy production and hormonal regulation pathways, were transiently and specifically induced in the late phase. In the comparison of salt tolerance between the two species, the flavonoid and brassinosteroid pathways were significantly enriched. More specifically, in the flavonoid pathway, FLS and F3'5'H exhibited significant differential expression. In the brassinosteroid pathway, the expression levels of DWF4, BR6OX2 and ROT3 were notably higher in P. pruinosa than in P. euphratica. Our results describe transcript dynamics and highlight secondary metabolite pathways involved in the response to salt stress during the seed germination of two desert poplars.

The draft genome sequence of a desert tree Populus pruinosa

Populus pruinosa is a large tree that grows in deserts and shows distinct differences in both morphology and adaptation compared to its sister species, P. euphratica. Here we present a draft genome sequence for P. pruinosa and examine genomic variations between the 2 species. A total of 60 Gb of clean reads from whole-genome sequencing of a P. pruinosa individual were generated using the Illumina HiSeq2000 platform. The assembled genome is 479.3 Mb in length, with an N50 contig size of 14.0 kb and a scaffold size of 698.5 kb; 45.47% of the genome is composed of repetitive elements. We predicted 35 131 protein-coding genes, of which 88.06% were functionally annotated. Gene family clustering revealed 224 unique and 640 expanded gene families in the P. pruinosa genome. Further evolutionary analysis identified numerous genes with elevated values for pairwise genetic differentiation between P. pruinosa and P. euphratica. We provide the genome sequence and gene annotation for P. pruinosa. A large number of genetic variations were recovered by comparison of the genomes between P. pruinosa and P. euphratica. These variations will provide a valuable resource for studying the genetic bases for the phenotypic and adaptive divergence of the 2 sister species.

A comparative transcriptomic analysis reveals the core genetic components of salt and osmotic stress responses in Braya humilis

Braya humilis is a member of the Euclidieae tribe within the family Brassicaceae. This species exhibits a broad range of adaptations to different climatic zones and latitudes as it has a distribution that ranges from northern Asia to the arctic-alpine regions of northern North America. In China, B. humilis is mainly found on the Qinghai-Tibetan Plateau (QTP) and in adjacent arid regions. In this study, we sequenced a sample from an arid region adjacent to the QTP using the Illumina platform generating a total of 46,485 highly accurate unigenes, of which 78.41% were annotated by BLASTing versus public protein databases. The B. humilis transcriptome is characterized by a high level of sequence conservation compared with its close relative, Arabidopsis thaliana. We also used reciprocal blast to identify shared orthologous genes between B. humilis and four other sequenced Brassicaceae species (i.e. A. thaliana, A. lyrata, Capsella rubella, and Thellungiella parvula). To enable precise characterization of orthologous genes, the early-diverging basal angiosperm Amborella trichopoda was also included. A total of 6,689 orthologous genes were identified before stricter criteria for the determination of e-values, amino acid hit lengths, and identity values was applied to further reduce this list. This led to a final list of 381 core orthologous genes for B. humilis; 39 out of these genes are involved in salt and osmotic stress responses and estimations of nonsynonymous/synonymous substitution ratios for this species and A. thaliana orthologs show that these genes are under purifying selection in B. humilis. Expression of six genes was detected in B. humilis seedlings under salt and osmotic stress treatments. Comparable expression patterns to their counterparts in Arabidopsis suggest that these orthologous genes are both sequence and functional conservation. The results of this study demonstrate that the environmental adaptations of B. humilis are mainly the results of preexisting genetic components. Future work will be required to characterize the expression patterns of these orthologous genes in natural populations and will provide further insights into the adaptive mechanisms underlying the wide range of B. humilis adaptations.

Comparative Transcriptome Analysis Reveals Adaptive Evolution of Notopterygium incisum and Notopterygium franchetii, Two High-Alpine Herbal Species Endemic to China

The extreme conditions (e.g., cold, low oxygen, and strong ultraviolet radiation) of the high mountains provide an ideal natural laboratory for studies on speciation and the adaptive evolution of organisms. Up to now, few genome/transcriptome-based studies have been carried out on how plants adapt to conditions at extremely high altitudes. Notopterygium incisum and Notopterygiumfranchetii (Notopterygium, Apiaceae) are two endangered high-alpine herbal plants endemic to China. To explore the molecular genetic mechanisms of adaptation to high altitudes, we performed high-throughput RNA sequencing (RNA-seq) to characterize the transcriptomes of the two species. In total, more than 130 million sequence reads, 81,446 and 63,153 unigenes with total lengths of 86,924,837 and 62,615,693 bp, were generated for the two herbal species, respectively. OrthoMCL analysis identified 6375 single-copy orthologous genes between N. incisum and N. franchetii. In total, 381 positively-selected candidate genes were identified for both plants by using estimations of the non-synonymous to synonymous substitution rate. At least 18 of these genes potentially participate in RNA splicing, DNA repair, glutathione metabolism and the plant–pathogen interaction pathway, which were further enriched in various functional gene categories possibly responsible for environment adaptation in high mountains. Meanwhile, we detected various transcription factors that regulated the material and energy metabolism in N. incisum and N. franchetii, which probably play vital roles in the tolerance to stress in surroundings. In addition, 60 primer pairs based on orthologous microsatellite-containing sequences between the both Notopterygium species were determined. Finally, 17 polymorphic microsatellite markers (SSR) were successfully characterized for the two endangered species. Based on these candidate orthologous and SSR markers, we detected that the adaptive evolution and species divergence of N. incisum and N. franchetii were significantly associated with the extremely heterogeneous environments and climatic oscillations in high-altitude areas. This work provides important insights into the molecular mechanisms of adaptation to high-altitudes in alpine herbal plants.

Comparative Transcriptome and Chloroplast Genome Analyses of Two Related Dipteronia Species

Dipteronia (order Sapindales) is an endangered genus endemic to China and has two living species, D.sinensis and D. dyeriana. The plants are closely related to the genus Acer, which is also classified in the order Sapindales. Evolutionary studies on Dipteronia have been hindered by the paucity of information on their genomes and plastids. Here, we used next generation sequencing to characterize the transcriptomes and complete chloroplast genomes of both Dipteronia species. A comparison of the transcriptomes of both species identified a total of 7814 orthologs. Estimation of selection pressures using Ka/Ks ratios showed that only 30 of 5435 orthologous pairs had a ratio significantly >1, i.e., showing positive selection. However, 4041 orthologs had a Ka/Ks < 0.5 (p < 0.05), suggesting that most genes had likely undergone purifying selection. Based on orthologous unigenes, 314 single copy nuclear genes (SCNGs) were identified. Through a combination of de novo and reference guided assembly, plastid genomes were obtained; that of D. sinensis was 157,080 bp and that of D. dyeriana was 157,071 bp. Both plastid genomes encoded 87 protein coding genes, 40 tRNAs, and 8 rRNAs; no significant differences were detected in the size, gene content, and organization of the two plastomes. We used the whole chloroplast genomes to determine the phylogeny of D. sinensis and D. dyeriana and confirmed that the two species were highly divergent. Overall, our study provides comprehensive transcriptomic and chloroplast genomic resources, which will be valuable for future evolutionary studies of Dipteronia.

Genome-wide transcriptional and physiological responses to drought stress in leaves and roots of two willow genotypes

BACKGROUND

Drought is a major environmental stress that can have severe impacts on plant productivity and survival. Understanding molecular mechanisms of drought responses is crucial in order to breed for drought adapted plant cultivars. The aim of the present study was to investigate phenotypic and transcriptional drought responses in two willow genotypes (520 and 592) originating from an experimental cross between S. viminalis × (S. viminalis × S. schwerinii). Willows are woody perennials in the Salicaceae plant family that are grown as bioenergy crops worldwide.

METHODS

An experiment was conducted where plants were exposed to drought and different eco-physiological parameters were assessed. RNA-seq data was furthermore generated with the Illumina technology from root tips and leaves from plants grown in drought and well-watered (WW) conditions. The RNA-seq data was assembled de novo with the Trinity assembler to create a reference gene set to which the reads were mapped in order to obtain differentially expressed genes (DEGs) between the drought and WW conditions. To investigate molecular mechanisms involved in the drought response, GO enrichment analyses were conducted. Candidate genes with a putative function in the drought response were also identified.

RESULTS

A total of 52,599 gene models were obtained and after filtering on gene expression (FPKM ≥ 1), 35,733 gene models remained, of which 24,421 contained open reading frames. A total of 5,112 unique DEGs were identified between drought and WW conditions, of which the majority were found in the root tips. Phenotypically, genotype 592 displayed less growth reduction in response to drought compared to genotype 520. At the transcriptional level, genotype 520 displayed a greater response in the leaves as more DEGs were found in genotype 520 compared to genotype 592. In contrast, the transcriptional responses in the root tips were rather similar between the two genotypes. A core set of candidate genes encoding proteins with a putative function in drought response was identified, for example MYBs and bZIPs as well as chlorophyll a/b binding proteins.

DISCUSSION

We found substantial differences in drought responses between the genotypes, both at the phenotypic and transcriptional levels. In addition to the genotypic variation in several traits, we also found indications for genotypic variation in trait plasticity, which could play a role in drought adaptation. Furthermore, the two genotypes displayed overall similar transcriptional responses in the root tips, but more variation in the leaves. It is thus possible that the observed phenotypic differences could be a result of transcriptional differences mostly at the leaf level.

CONCLUSIONS

This study has contributed to a better general understanding of drought responses in woody plants, specifically in willows, and has implications for breeding research towards more drought adapted plants.

On the salty side of life: molecular, physiological and anatomical adaptation and acclimation of trees to extreme habitats

Saline and sodic soils that cannot be used for agriculture occur worldwide. Cultivating stress-tolerant trees to obtain biomass from salinized areas has been suggested. Various tree species of economic importance for fruit, fibre and timber production exhibit high salinity tolerance. Little is known about the mechanisms enabling tree crops to cope with high salinity for extended periods. Here, the molecular, physiological and anatomical adjustments underlying salt tolerance in glycophytic and halophytic model tree species, such as Populus euphratica in terrestrial habitats, and mangrove species along coastlines are reviewed. Key mechanisms that have been identified as mediating salt tolerance are discussed at scales from the genetic to the morphological level, including leaf succulence and structural adjustments of wood anatomy. The genetic and transcriptomic bases for physiological salt acclimation are salt sensing and signalling networks that activate target genes; the target genes keep reactive oxygen species under control, maintain the ion balance and restore water status. Evolutionary adaptation includes gene duplication in these pathways. Strategies for and limitations to tree improvement, particularly transgenic approaches for increasing salt tolerance by transforming trees with single and multiple candidate genes, are discussed.

Genome-wide Identification of WRKY Genes in the Desert Poplar Populus euphratica and Adaptive Evolution of the Genes in Response to Salt Stress

WRKY transcription factors play important roles in plant development and responses to various stresses in plants. However, little is known about the evolution of the WRKY genes in the desert poplar species Populus euphratica, which is highly tolerant of salt stress. In this study, we identified 107 PeWRKY genes from the P. euphratica genome and examined their evolutionary relationships with the WRKY genes of the salt-sensitive congener Populus trichocarpa. Ten PeWRKY genes are specific to P. euphratica, and five of these showed altered expression under salt stress. Furthermore, we found that two pairs of orthologs between the two species showed evidence of positive evolution, with dN/dS ratios>1 (nonsynonymous/synonymous substitutions), and both of them altered their expression in response to salinity stress. These findings suggested that both the development of new genes and positive evolution in some orthologs of the WRKY gene family may have played an important role in the acquisition of high salt tolerance by P. euphratica.

Transcriptome sequencing of three Ranunculus species (Ranunculaceae) reveals candidate genes in adaptation from terrestrial to aquatic habitats

Adaptation to aquatic habitats is a formidable challenge for terrestrial angiosperms that has long intrigued scientists. As part of a suite of work to explore the molecular mechanism of adaptation to aquatic habitats, we here sequenced the transcriptome of the submerged aquatic plant Ranunculus bungei, and two terrestrial relatives R. cantoniensis and R. brotherusii, followed by comparative evolutionary analyses to determine candidate genes for adaption to aquatic habitats. We obtained 126,037, 140,218 and 114,753 contigs for R. bungei, R. cantoniensis and R. brotherusii respectively. Bidirectional Best Hit method and OrthoMCL method identified 11,362 and 8,174 1:1:1 orthologous genes (one ortholog is represented in each of the three species) respectively. Non-synonymous/synonymous (dN/dS) analyses were performed with a maximum likelihood method and an approximate method for the three species-pairs. In total, 14 genes of R. bungei potentially involved in the adaptive transition from terrestrial to aquatic habitats were identified. Some of the homologs to these genes in model plants are involved in vacuole protein formation, regulating 'water transport process' and 'microtubule cytoskeleton organization'. Our study opens the door to understand the molecular mechanism of plant adaptation from terrestrial to aquatic habitats.

Transcriptome assembly, profiling and differential gene expression analysis of the halophyte Suaeda fruticosa provides insights into salt tolerance

BACKGROUND

Improvement of crop production is needed to feed the growing world population as the amount and quality of agricultural land decreases and soil salinity increases. This has stimulated research on salt tolerance in plants. Most crops tolerate a limited amount of salt to survive and produce biomass, while halophytes (salt-tolerant plants) have the ability to grow with saline water utilizing specific biochemical mechanisms. However, little is known about the genes involved in salt tolerance. We have characterized the transcriptome of Suaeda fruticosa, a halophyte that has the ability to sequester salts in its leaves. Suaeda fruticosa is an annual shrub in the family Chenopodiaceae found in coastal and inland regions of Pakistan and Mediterranean shores. This plant is an obligate halophyte that grows optimally from 200-400 mM NaCl and can grow at up to 1000 mM NaCl. High throughput sequencing technology was performed to provide understanding of genes involved in the salt tolerance mechanism. De novo assembly of the transcriptome and analysis has allowed identification of differentially expressed and unique genes present in this non-conventional crop.

RESULTS

Twelve sequencing libraries prepared from control (0 mM NaCl treated) and optimum (300 mM NaCl treated) plants were sequenced using Illumina Hiseq 2000 to investigate differential gene expression between shoots and roots of Suaeda fruticosa. The transcriptome was assembled de novo using Velvet and Oases k-45 and clustered using CDHIT-EST. There are 54,526 unigenes; among these 475 genes are downregulated and 44 are upregulated when samples from plants grown under optimal salt are compared with those grown without salt. BLAST analysis identified the differentially expressed genes, which were categorized in gene ontology terms and their pathways.

CONCLUSIONS

This work has identified potential genes involved in salt tolerance in Suaeda fruticosa, and has provided an outline of tools to use for de novo transcriptome analysis. The assemblies that were used provide coverage of a considerable proportion of the transcriptome, which allows analysis of differential gene expression and identification of genes that may be involved in salt tolerance. The transcriptome may serve as a reference sequence for study of other succulent halophytes.

Transcriptome Profiles of Populus euphratica upon Heat Shock stress

Heat stress, which strongly affects plant performance and often results in reduced vegetative growth and yields depression, has become an increasingly serious global problem. Populus euphratica Oliv. which has been considered as a tree model for the study of plant response to abiotic stresses, could be resistant to an extremely wide environmental temperature range (–40 °C to 45 °C). Previous study is mainly focused on its gene regulation upon drought and salt stress. However, little is known about gene regulation at the global transcriptome level upon heat stress. To understand the gene network controlling heat stress in P. euphratica, a transcriptome sequencing using Illumina Hiseq 2000 was performed to generate a 10 gigabases depth for each sample in the tissue of leaf. 119,573 unigeneswere generated with an average length of 474 bp. Approximately 49,605 (41.49%) unigenes exhibited significantly different expressions between two libraries. Among these unigenes, 11,165 (9.34%) were upregulated and 38,440 (32.15%) were down regulated. Heat shock proteins classified as molecular chaperones showed a significant percentage (1.13%) in the up regulated group. Heat responsive genes, such as polyubiquitins, were over expressed in heat treated sample. GO enrichment analysis revealed that the Go terms for differentially expressed unigenes were significantly enriched in hormone-mediated signal, biological process regulation and metabolic process regulation. Our data revealed a global transcriptome picture of P. euphratica in response to heat shock. The identified potential heat stress-related transcripts can be used to infer the gene regulation networks underlying the molecular mechanisms of heat response in P. euphratica.

Response to drought and salt stress in leaves of poplar (Populus alba × Populus glandulosa): expression profiling by oligonucleotide microarray analysis

Drought and salt stresses are major environmental constraints on forest productivity. To identify genes responsible for stress tolerance, we conducted a genome-wide analysis in poplar (Populus alba × Populus glandulosa) leaves exposed to drought and salt (NaCl) stresses. We investigated gene expression at the mRNA level using oligonucleotide microarrays containing 44,718 genes from Populus trichocarpa. A total of 1604 and 1042 genes were up-regulated (≥2-fold; P value < 0.05) by drought and salt stresses, respectively, and 765 genes were up-regulated by both stresses. In addition, 2742 and 1685 genes were down-regulated by drought and salt stresses, respectively, and 1564 genes were down-regulated by both stresses. The large number of genes regulated by both stresses suggests that crosstalk occurs between the drought and salt stress responses. Most up-regulated genes were involved in functions such as subcellular localization, signal transduction, metabolism, and transcription. Among the up-regulated genes, we identified 47 signaling proteins, 65 transcription factors, and 43 abiotic stress-related genes. Several genes were modulated by only one of the two stresses. About 25% of the genes significantly regulated by these stresses are of unknown function, suggesting that poplar may provide an opportunity to discover novel stress-related genes.

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

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.

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.