Tissue-Specific Transcriptome Analysis Reveals Multiple Responses to Salt Stress in Populus euphratica Seedlings

The Euphrates Poplar Responses to Abiotic Stress and Its Unique Traits in Dry Regions of China (Xinjiang and Inner Mongolia): What Should We Know?

At the moment, drought, salinity, and low-temperature stress are ubiquitous environmental issues. In arid regions including Xinjiang and Inner Mongolia and other areas worldwide, the area of tree plantations appears to be rising, triggering tree growth. Water is a vital resource in the agricultural systems of countries impacted by aridity and salinity. Worldwide efforts to reduce quantitative yield losses on Populus euphratica by adapting tree plant production to unfavorable environmental conditions have been made in response to the responsiveness of the increasing control of water stress. Although there has been much advancement in identifying the genes that resist abiotic stresses, little is known about how plants such as P. euphratica deal with numerous abiotic stresses. P. euphratica is a varied riparian plant that can tolerate drought, salinity, low temperatures, and climate change, and has a variety of water stress adaptability abilities. To conduct this review, we gathered all available information throughout the Web of Science, the Chinese National Knowledge Infrastructure, and the National Center for Biotechnology Information on the impact of abiotic stress on the molecular mechanism and evolution of gene families at the transcription level. The data demonstrated that P. euphratica might gradually adapt its stomatal aperture, photosynthesis, antioxidant activities, xylem architecture, and hydraulic conductivity to endure extreme drought and salt stress. Our analyses will give readers an understanding of how to manage a gene family in desert trees and the influence of abiotic stresses on the productivity of tree plants. They will also give readers the knowledge necessary to improve biotechnology-based tree plant stress tolerance for sustaining yield and quality trees in China's arid regions.

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

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

SMRT and Illumina sequencing provide insights into mechanisms of lignin and terpenoids biosynthesis in Pinus massoniana Lamb

Wood and oleoresin are important industrial raw materials with high economic value; however, their molecular formation and biosynthesis mechanisms in different tissues of Pinus massoniana remain unexplored. Therefore, we used single-molecule real-time sequencing technology (SMRT) and Illumina RNA sequencing to establish a transcriptome dataset and explore the expression pattern of genes related to secondary metabolites involved in wood formation and oleoresin biosynthesis in six different P. massoniana tissues. In total, 63.58 Gb of polymerase reads were obtained, including 41,407 isoforms with an average length of 1822 bp. We identified 3939 and 8785 isoforms and 161 and 481 transcription factors with tissue expression specificity and in the reproductive and vegetative organs, respectively. Eighty isoforms were annotated as cellulose synthases and 224 isoforms involved in lignin biosynthesis were enriched. Additionally, we identified 217 isoforms involved in the terpenoid biosynthesis pathway, with needles having the most tissue-specific genes for terpenoid biosynthesis. Some isoforms related to lignin biosynthesis were highly expressed in the xylem, according to the results of transcriptome sequencing and real-time quantitative reverse-transcription polymerase chain reaction. Our research confirmed the advantages of SMRT sequencing and provided valuable information for the transcriptional annotation of P. massoniana, which will be beneficial for producing better raw wood and oleoresin materials.

Conjoint Analysis of Genome-Wide lncRNA and mRNA Expression during the Salicylic Acid Response in Populus × euramericana

Long noncoding RNAs (lncRNAs) participate in a wide range of biological processes, but lncRNAs in plants remain largely unknown; in particular, we lack a systematic identification of plant lncRNAs involved in hormone responses. To explore the molecular mechanism of the response of poplar to salicylic acid (SA), the changes in protective enzymes, which are closely related to plant resistance induced by exogenous SA, were studied, and the expression of mRNA and lncRNA were determined by high-throughput RNA sequencing. The results showed that the activities of phenylalanine ammonia lyase (PAL) and polyphenol oxidase (PPO), in the leaves of Populus × euramericana, were significantly increased by exogenous SA application. High-throughput RNA sequencing showed that 26,366 genes and 5690 lncRNAs were detected under the different treatment conditions: SA and H₂O application. Among these, 606 genes and 49 lncRNAs were differentially expressed. According to target prediction, lncRNAs and target genes involved in light response, stress response, plant disease resistance, and growth and development, were differentially expressed in SA-treated leaves. Interaction analysis showed that lncRNA-mRNA interactions, following exogenous SA, were involved in the response of poplar leaves to the external environment. Our study provides a comprehensive view of Populus × euramericana lncRNAs and offers insights into the potential functions and regulatory interactions of SA-responsive lncRNAs, thus forming the foundation for future functional analysis of SA-responsive lncRNAs in Populus × euramericana.

A phenylalanine ammonia lyase from Fritillaria unibracteata promotes drought tolerance by regulating lignin biosynthesis and SA signaling pathway

Drought is one of the key threatening environmental factors for plant and agriculture. Phenylalanine ammonia lyase (PAL) is a key enzyme involved in plant defense against abiotic stress, however, the role of PAL in drought tolerance remains elusive. Here, a PAL member (FuPAL1) containing noncanonical Ala-Ser-Gly triad was isolated from Fritillaria unibracteata, one important alpine pharmaceutical plant. FuPAL1, mainly distributed in cytosol, was more conserved than FuCOMT and FuCHI at both nucleotide and amino acid levels. FuPAL1 was overexpressed in Escherichia coli and the purified recombinant FuPAL1 protein showed catalytic preference on L-Phe than L-Tyr. Homology modeling and site-mutation of FuPAL1 exhibited FuPAL1 took part in the ammonization process by forming MIO-like group, and Phe141, Ser208, Ileu218 and Glu490 played key roles in substrate binding and (or) catalysis. HPLC analysis showed that lignin and salicylic acid levels increased but total flavonoid levels decreased in FuPAL1 transgenic Arabidopsis compared to wild-type plants. Moreover, FuPAL1 transgenic Arabidopsis significantly enhanced its drought tolerance, which suggested that FuPAL1 mediated tolerance to drought by inducing the biosynthesis and accumulation of salicylic acid and lignin. Taken together, our results confirmed that the FuPAL1 played an important role in drought tolerance, and FuPAL1 might be a valuable target for genetic improvement of drought resistance in future.

Comparative Transcriptomic Analysis of Root and Leaf Transcript Profiles Reveals the Coordinated Mechanisms in Response to Salinity Stress in Common Vetch

Owing to its strong environmental suitability to adverse abiotic stress conditions, common vetch (Vicia sativa) is grown worldwide for both forage and green manure purposes and is an important protein source for human consumption and livestock feed. The germination of common vetch seeds and growth of seedlings are severely affected by salinity stress, and the response of common vetch to salinity stress at the molecular level is still poorly understood. In this study, we report the first comparative transcriptomic analysis of the leaves and roots of common vetch under salinity stress. A total of 6361 differentially expressed genes were identified in leaves and roots. In the roots, the stress response was dominated by genes involved in peroxidase activity. However, the genes in leaves focused mainly on Ca²⁺ transport. Overexpression of six salinity-inducible transcription factors in yeast further confirmed their biological functions in the salinity stress response. Our study provides the most comprehensive transcriptomic analysis of common vetch leaf and root responses to salinity stress. Our findings broaden the knowledge of the common and distinct intrinsic molecular mechanisms within the leaves and roots of common vetch and could help to develop common vetch cultivars with high salinity tolerance.

Regulatory short RNAs: A decade's tale for manipulating salt tolerance in plants

Salt stress is a globally increasing environmental detriment to crop growth and productivity. Exposure to salt stress evokes a complex medley of cellular signals, which rapidly reprogram transcriptional and metabolic networks to shape plant phenotype. To date, genetic engineering approaches were used with success to enhance salt tolerance; however, their performance is yet to be evaluated under realistic field conditions. Regulatory short non-coding RNAs (rsRNAs) are emerging as next-generation candidates for engineering salt tolerance in crops. In view of this, the present review provides a comprehensive analysis of a decade's worth of functional studies on non-coding RNAs involved in salt tolerance. Further, we have integrated this knowledge of rsRNA-mediated regulation with the current paradigm of salt tolerance to highlight two regulatory complexes (RCs) for regulating salt tolerance in plants. Finally, a knowledge-driven roadmap is proposed to judiciously utilize RC component(s) for enhancing salt tolerance in crops.

Molecular dissection of maize seedling salt tolerance using a genome-wide association analysis method

Salt stress is a major devastating abiotic factor that affects the yield and quality of maize. However, knowledge of the molecular mechanisms of the responses to salt stress in maize is limited. To elucidate the genetic basis of salt tolerance traits, a genome-wide association study was performed on 348 maize inbred lines under normal and salt stress conditions using 557 894 single nucleotide polymorphisms (SNPs). The phenotypic data for 27 traits revealed coefficients of variation of >25%. In total, 149 significant SNPs explaining 6.6%-11.2% of the phenotypic variation for each SNP were identified. Of the 104 identified quantitative trait loci (QTLs), 83 were related to salt tolerance and 21 to normal traits. Additionally, 13 QTLs were associated with two to five traits. Eleven and six QTLs controlling salt tolerance traits and normal root growth, respectively, co-localized with QTL intervals reported previously. Based on functional annotations, 13 candidate genes were predicted. Expression levels analysis of 12 candidate genes revealed that they were all responsive to salt stress. The CRISPR/Cas9 technology targeting three sites was applied in maize, and its editing efficiency reached 70%. By comparing the biomass of three CRISPR/Cas9 mutants of ZmCLCg and one zmpmp3 EMS mutant with their wild-type plants under salt stress, the salt tolerance function of candidate genes ZmCLCg and ZmPMP3 were confirmed. Chloride content analysis revealed that ZmCLCg regulated chloride transport under sodium chloride stress. These results help to explain genetic variations in salt tolerance and provide novel loci for generating salt-tolerant maize lines.

Tissue tolerance mechanisms conferring salinity tolerance in a halophytic perennial species Nitraria sibirica Pall

Plant salt tolerance relies on a coordinated functioning of different tissues and organs. Salinity tissue tolerance is one of the key traits that confer plant adaptation to saline environment. This trait implies maintenance low cytosolic Na+/K+ ratio in metabolically active cellular compartments. In this study, we used Nitraria sibirica Pall., a perennial woody halophyte species, to understand the mechanistic basis of its salinity tissue tolerance. The results showed that the growth of seedlings was stimulated by 100-200 mM NaCl treatment. The ions distribution analysis showed that the leaves act as an Na+ sink, while the plant roots possess superior K+ retention. The excessive Na+ absorbed from the soil was mainly transported to the shoot and was eventuallysequestrated into mesophyll vacuoles in the leaves. As a result, N. sibirica could keep the optimal balance of K+/Na+ at a tissue- and cell-specific level under saline condition. To enable this, N. sibirica increased both vacuolar H+-ATPase and H+-PPase enzymes activities and up-regulated the expressions of NsVHA, NsVP1 and NsNHX1 genes. Vacuolar Na+ sequestration in the leaf mesophyll, mediated by NsVHA, NsVP1 and NsNHX1, reduced the Na+ concentration in cytosol and inhibited further K+ loss. Meanwhile, N. sibirica enhanced the Two Pore K+ expression at the transcriptional level to promote K+ efflux from vacuole into cytoplasm, assisting in maintaining cytosolic K+ homeostasis. It is concluded that the tissue tolerance traits such as vacuolar Na+ sequestration and intracellular K+ homeostasis are critical to confer adaptation of N. sibirica to soil salinity.

A gene that underwent adaptive evolution, LAC2 (LACCASE), in Populus euphratica improves drought tolerance by improving water transport capacity

Drought severely limits plant development and growth; accordingly, plants have evolved strategies to prevent water loss and adapt to water deficit conditions. However, experimental cases that corroborate these evolutionary processes are limited. The LACCASEs (LACs) family is involved in various plant development and growth processes. Here, we performed an evolutionary analysis of LACs from Populus euphratica and characterized the functions of LACs in Arabidopsis and poplar. The results showed that in PeuLACs, multiple gene duplications led to apparent functional redundancy as the result of various selective pressures. Among them, PeuLAC2 underwent strong positive selection. Heterologous expression analyses showed that the overexpression of PeuLAC2 alters the xylem structure of plants, including thickening the secondary cell wall (SCW) and increasing the fiber cell length and stem tensile strength. Altogether, these changes improve the water transport capacity of plants. The analysis of the physiological experimental results showed that PeuLAC2-OE lines exhibited a stronger antioxidant response and greater drought tolerance than WT. Three genes screened by transcriptome analysis, NAC025, BG1, and UGT, that are associated with SCW synthesis and drought stress were all upregulated in the PeuLAC2-OE lines, implying that the overexpression of PeuLAC2 thickened the SCW, improved the water transport capacity of the plant, and further enhanced its drought tolerance. Our study highlights that genes that have undergone adaptive evolution may participate in the development of adaptive traits in P. euphratica and that PeuLAC2 could be a candidate gene for molecular genetic breeding in trees.

Full-Length Transcriptome Analysis of Four Different Tissues of Cephalotaxus oliveri

Cephalotaxus oliveri is a tertiary relict conifer endemic to China, regarded as a national second-level protected plant in China. This species has experienced severe changes in temperature and precipitation in the past millions of years, adapting well to harsh environments. In view of global climate change and its endangered conditions, it is crucial to study how it responds to changes in temperature and precipitation for its conservation work. In this study, single-molecule real-time (SMRT) sequencing and Illumina RNA sequencing were combined to generate the complete transcriptome of C. oliveri. Using the RNA-seq data to correct the SMRT sequencing data, the four tissues obtained 63,831 (root), 58,108 (stem), 33,013 (leaf) and 62,436 (male cone) full-length unigenes, with a N50 length of 2523, 3480, 3181, and 3267 bp, respectively. Additionally, 35,887, 11,306, 36,422, and 25,439 SSRs were detected for the male cone, leaf, root, and stem, respectively. The number of long non-coding RNAs predicted from the root was the largest (11,113), and the other tissues were 3408 (stem), 3193 (leaf), and 3107 (male cone), respectively. Functional annotation and enrichment analysis of tissue-specific expressed genes revealed the special roles in response to environmental stress and adaptability in the different four tissues. We also characterized the gene families and pathways related to abiotic factors. This work provides a comprehensive transcriptome resource for C. oliveri, and this resource will facilitate further studies on the functional genomics and adaptive evolution of C. oliveri.

A Tissue-Specific Landscape of Alternative Polyadenylation, lncRNAs, TFs, and Gene Co-expression Networks in Liriodendron chinense

Liriodendron chinense is an economically and ecologically important deciduous tree species. Although the reference genome has been revealed, alternative polyadenylation (APA), transcription factors (TFs), long non-coding RNAs (lncRNAs), and co-expression networks of tissue-specific genes remain incompletely annotated. In this study, we used the bracts, petals, sepals, stamens, pistils, leaves, and shoot apex of L. chinense as materials for hybrid sequencing. On the one hand, we improved the annotation of the genome. We detected 13,139 novel genes, 7,527 lncRNAs, 1,791 TFs, and 6,721 genes with APA sites. On the other hand, we found that tissue-specific genes play a significant role in maintaining tissue characteristics. In total, 2,040 tissue-specific genes were identified, among which 9.2% of tissue-specific genes were affected by APA, and 1,809 tissue-specific genes were represented in seven specific co-expression modules. We also found that bract-specific hub genes were associated plant defense, leaf-specific hub genes were involved in energy metabolism. Moreover, we also found that a stamen-specific hub TF Lchi25777 may be involved in the determination of stamen identity, and a shoot-apex-specific hub TF Lchi05072 may participate in maintaining meristem characteristic. Our study provides a landscape of APA, lncRNAs, TFs, and tissue-specific gene co-expression networks in L. chinense that will improve genome annotation, strengthen our understanding of transcriptome complexity, and drive further research into the regulatory mechanisms of tissue-specific genes.

TEnGExA: an R package based tool for tissue enrichment and gene expression analysis

RNA-seq data analysis with rapidly advancing high-throughput sequencing technology, nowadays provides large number of transcripts or genes to perform downstream analysis including functional annotation and pathway analysis. However for the data from multiple tissues, downstream analysis with tissue-specific or tissue-enriched transcripts is highly preferable. However, there is still a need of tool for quickly performing tissue-enrichment and gene expression analysis irrespective of number of input genes or tissues at various fragments per kilobase of transcript per million fragments mapped (FPKM) thresholds. To fulfill this need, we presented a freely available R package and web-interface tool, TEnGExA, which allows tissue-enrichment analysis (TEA) for any number of genes or transcripts for any species provided only a read-count or FPKM-value matrix as input. Based on the different FPKM value and fold thresholds, TEnGExA classifies the user provided gene lists into tissue-enriched or tissue-specific transcripts along with other standard classes. By analyzing the published sample data from human, plant and microorganism, we signifies that TEnGExA can easily handle complex or large data from any species to provided tissue-enriched gene list for downstream analysis in quick time. In summary, TEnGExA is quick, easy to use and an efficient tool for TEA. The R package is freely available at https://github.com/ubagithub/TEnGExA/ and the GUI web interface is accessible at http://webtom.cabgrid.res.in/tissue_enrich/.

Transcriptome analysis of crude oil degrading Pseudomonas aeruginosa strains for identification of potential genes involved in crude oil degradation

In the microbial world, bacteria are the most effective agents in petroleum hydrocarbons (PHs) degradation, utilization/mineralization and they serve as essential degraders of crude oil contaminated environment. Some genes and traits are involved in the hydrocarbon utilization process for which transcriptome analyses are important to identify differentially expressed genes (DEGs) among different conditions, leading to a new understanding of genes or pathways associated with crude oil degradation. In this work, three crude oil utilizing Pseudomonas aeruginosa strains designated as N002, TP16 and J001 subjected to transcriptome analyses revealed a total of 81, 269 and 137 significant DEGs. Among them are 80 up-regulated genes and one downregulated gene of N002, 121 up- regulated and 148 down-regulated genes of TP16, 97 up-regulated and 40 down-regulated genes of J001 which are involved in various metabolic pathways. TP16 strain has shown more number of DEGs upon crude oil treatment in comparison to the other two strains. Through quantitative real time polymerase chain reaction (qRT-PCR), the selected DEGs of each strain from transcriptome data were substantiated. The results have shown that the up- regulated and down-regulated genes observed by qRT-PCR were consistent with transcriptome data. Taken together, our transcriptome results have revealed that TP16 is a potential P. aeruginosa strain for functional analysis of identified potential DEGs involved in crude oil degradation.

Improved genome assembly provides new insights into genome evolution in a desert poplar (Populus euphratica)

Populus euphratica is well adapted to extreme desert environments and is an important model species for elucidating the mechanisms of abiotic stress resistance in trees. The current assembly of P. euphratica genome is highly fragmented with many gaps and errors, thereby impeding downstream applications. Here, we report an improved chromosome-level reference genome of P. euphratica (v2.0) using single-molecule sequencing and chromosome conformation capture (Hi-C) technologies. Relative to the previous reference genome, our assembly represents a nearly 60-fold improvement in contiguity, with a scaffold N50 size of 28.59 Mb. Using this genome, we have found that extensive expansion of Gypsy elements in P. euphratica led to its rapid increase in genome size compared to any other Salicaceae species studied to date, and potentially contributed to adaptive divergence driven by insertions near genes involved in stress tolerance. We also detected a wide range of unique structural rearrangements in P. euphratica, including 2,549 translocations, 454 inversions, 121 tandem and 14 segmental duplications. Several key genes likely to be involved in tolerance to abiotic stress were identified within these regions. This high-quality genome represents a valuable resource for poplar breeding and genetic improvement in the future, as well as comparative genomic analysis with other Salicaceae species.

An AP2/ERF transcription factor ERF139 coordinates xylem cell expansion and secondary cell wall deposition

Differentiation of xylem elements involves cell expansion, secondary cell wall (SCW) deposition and programmed cell death. Transitions between these phases require strict spatiotemporal control. The function of Populus ERF139 (Potri.013G101100) in xylem differentiation was characterized in transgenic overexpression and dominant repressor lines of ERF139 in hybrid aspen (Populus tremula × tremuloides). Xylem properties, SCW chemistry and downstream targets were analyzed in both types of transgenic trees using microscopy techniques, Fourier transform-infrared spectroscopy, pyrolysis-GC/MS, wet chemistry methods and RNA sequencing. Opposite phenotypes were observed in the secondary xylem vessel sizes and SCW chemistry in the two different types of transgenic trees, supporting the function of ERF139 in suppressing the radial expansion of vessel elements and stimulating accumulation of guaiacyl-type lignin and possibly also xylan. Comparative transcriptomics identified genes related to SCW biosynthesis (LAC5, LBD15, MYB86) and salt and drought stress-responsive genes (ANAC002, ABA1) as potential direct targets of ERF139. The phenotypes of the transgenic trees and the stem expression profiles of ERF139 potential target genes support the role of ERF139 as a transcriptional regulator of xylem cell expansion and SCW formation, possibly in response to osmotic changes of the cells.

Conjoint Analysis of Genome-Wide lncRNA and mRNA Expression of Heteromorphic Leavesin Response to Environmental Heterogeneityin Populus euphratica

Heterophylly is the phenomenon of leaf forms varying along the longitudinal axis within a single plant. Populus euphratica, a heterophyllous woody plant, develops lanceolate leaves and dentate broad-ovate leaves on the bottom and top of the canopy, respectively, which are faced with different intensities of ambient solar radiation. However, the mechanism of the heteromorphic leaf response to the microenvironment in P. euphratica remains elusive. Here, we show that the dentate broad-ovate leaves have advantages in tolerating high light intensity, while lanceolate leaves are excellent at capturing light. Compared with lanceolate leaves, more trichomes, higher stomatal density, thicker lamina, and higher specific leaf weight were observed in dentate broad-ovate leaves. Furthermore, high-throughput RNA sequencing analysis revealed that the expression patterns of genes and long noncoding RNAs (lncRNAs) are different between the two heteromorphic leaves. A total of 36,492 genes and 1725 lncRNAs were detected, among which 586 genes and 54 lncRNAs were differentially expressed. Based on targets prediction, lncRNAs and target genes involved in light adaption, protein repair, stress response, and growth and development pathways were differentially expressed in heteromorphic leaves, 10 pairs of which were confirmed by quantitative real-time PCR. Additionally, the analysis of interactions indicated that lncRNA-mRNA interactions were involved in the response to the microenvironment of heteromorphic leaves. Taken together, these results suggest that the morphological features and joint regulation of lncRNA-mRNA in heteromorphic leaves may serve as survival strategies for P. euphratica, which could lead to optimal utilization of environmental factors.

Genome-Wide Identification of Long Noncoding RNAs and Their Responses to Salt Stress in Two Closely Related Poplars

Long noncoding RNAs (lncRNAs) are involved in various biological regulatory processes, but their roles in plants resistance to salt stress remain largely unknown. To systematically explore the characteristics of lncRNAs and their roles in plant salt responses, we conducted strand-specific RNA-sequencing of four tissue types with salt treatments in two closely related poplars (Populus euphratica and Populus alba var. pyramidalis), and a total of 10,646 and 10,531 lncRNAs were identified, respectively. These lncRNAs showed significantly lower values in terms of length, expression, and expression correction than with mRNA. We further found that about 40% and 60% of these identified lncRNAs responded to salt stress with tissue-specific expression patterns across the two poplars. Furthermore, lncRNAs showed weak evolutionary conservation in sequences and exhibited diverse regulatory styles; in particular, tissue- and species-specific responses to salt stress varied greatly in two poplars, for example, 322 lncRNAs were found highly expressed in P. euphratica but not in P. alba var. pyramidalis and 3,425 lncRNAs were identified to be species-specific in P. euphratica in response to salt stress. Moreover, tissue-specific expression of lncRNAs in two poplars were identified with predicted target genes included Aux/IAA, NAC, MYB, involved in regulating plant growth and the plant stress response. Taken together, the systematic analysis of lncRNAs between sister species enhances our understanding of the characteristics of lncRNAs and their roles in plant growth and salt response.

Genome-wide transcriptional adaptation to salt stress in Populus

BACKGROUND

Adaptation to abiotic stresses is crucial for the survival of perennial plants in a natural environment. However, very little is known about the underlying mechanisms. Here, we adopted a liquid culture system to investigate plant adaptation to repeated salt stress in Populus trees.

RESULTS

We first evaluated phenotypic responses and found that plants exhibit better stress tolerance after pre-treatment of salt stress. Time-course RNA sequencing (RNA-seq) was then performed to profile changes in gene expression over 12 h of salt treatments. Analysis of differentially expressed genes (DEGs) indicated that significant transcriptional reprogramming and adaptation to repeated salt treatment occurred. Clustering analysis identified two modules of co-expressed genes that were potentially critical for repeated salt stress adaptation, and one key module for salt stress response in general. Gene Ontology (GO) enrichment analysis identified pathways including hormone signaling, cell wall biosynthesis and modification, negative regulation of growth, and epigenetic regulation to be highly enriched in these gene modules.

CONCLUSIONS

This study illustrates phenotypic and transcriptional adaptation of Populus trees to salt stress, revealing novel gene modules which are potentially critical for responding and adapting to salt stress.

Salt-Responsive Genes are Differentially Regulated at the Chromatin Levels Between Seedlings and Roots in Rice

The elucidation of epigenetic responses of salt-responsive genes facilitates understanding of the underlying mechanisms that confer salt tolerance in rice. However, it is still largely unknown how epigenetic mechanisms are associated with the expression of salt-responsive genes in rice and other crops. In this study, we reported tissue-specific gene expression and tissue-specific changes in chromatin modifications or signatures between seedlings and roots in response to salt treatment. Our study indicated that among six of individual mark examined (H3K4me3, H3K27me3, H4K12ac, H3K9ac, H3K27ac and H3K36me3), a positive association between salt-related changes in histone marks and the expression of differentially expressed genes (DEGs) was observed only for H3K9ac and H4K12ac in seedlings and H3K36me3 in roots. In contrast, chromatin states (CSs) with combinations of six histone modification marks played crucial roles in the differential expression of salt-responsive genes between seedlings and roots. Most importantly, CS7 containing the bivalent marks H3K4me3 and H3K27me3, with a mutual exclusion of functions with each other, displayed distinct functions in the expression of DEGs in both tissues. Specifically, H3K27me3 in CS7 mainly suppressed the expression of DEGs in roots, while H3K4me3 affected the expression of down- and up-regulated genes, possibly by antagonizing the repressive role of H3K27me3 in seedlings. Our findings indicate distinct impacts of the CSs on the differential expression of salt-responsive genes between seedlings and roots in rice, which provides an important background for understanding chromatin-based epigenetic mechanisms that might confer salt tolerance in plants.

Plant stress RNA-seq Nexus: a stress-specific transcriptome database in plant cells

Background

Abiotic and biotic stresses severely affect the growth and reproduction of plants and crops. Determining the critical molecular mechanisms and cellular processes in response to stresses will provide biological insight for addressing both climate change and food crises. RNA sequencing (RNA-Seq) is a revolutionary tool that has been used extensively in plant stress research. However, no existing large-scale RNA-Seq database has been designed to provide information on the stress-specific differentially expressed transcripts that occur across diverse plant species and various stresses.

Results

We have constructed a comprehensive database, the plant stress RNA-Seq nexus (PSRN), which includes 12 plant species, 26 plant-stress RNA-Seq datasets, and 937 samples. All samples are assigned to 133 stress-specific subsets, which are constructed into 254 subset pairs, a comparison between selected two subsets, for stress-specific differentially expressed transcript identification.

Conclusions

PSRN is an open resource for intuitive data exploration, providing expression profiles of coding-transcript/lncRNA and identifying which transcripts are differentially expressed between different stress-specific subsets, in order to support researchers generating new biological insights and hypotheses in molecular breeding or evolution. PSRN is freely available at http://syslab5.nchu.edu.tw/PSRN.

Electronic supplementary material

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