RNA Sequencing and Co-expressed Long Non-coding RNA in Modern and Wild Wheats

Genome-wide characterization of drought-responsive long non-coding RNAs in sorghum (Sorghum bicolor)

Drought stress strongly affects crop yield. Although knowledge of long non-coding RNAs (lncRNAs) has been updated continuously and rapidly, information about lncRNAs in drought resistance regulation is extremely limited in sorghum. Here, lncRNA-sequencing was performed with seedlings of a sorghum cultivar (Jinza29) under three water control treatments to investigate the mechanism of lncRNAs responsible for drought resistance in sorghum. A total of 377 differentially expressed lncRNAs (DElncRNAs) were identified. We also predicted 4322 and 2827 transcripts as potential cis-target and trans-target genes for drought-responsive lncRNAs, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that those target genes exhibited marked enrichment into "oxidoreductase activity", "signal transducer activity", "DNA repair", "photosynthesis", "glutathione metabolism", and "phenylpropanoid biosynthesis" and other terms associated with abiotic stress resistance. Moreover, several lncRNAs were estimated to modulate the expression of other genes related to stress response and photosynthetic carbon metabolism. Additionally, we found 107 DElncRNAs that might be candidate target mimics for 56 miRNAs. LncRNAs play important roles in drought adaptation of sorghum through interacting with protein-encoding genes. The obtained results provided novel insights into the biological characteristics of lncRNAs and offered potential regulatory factors for genetically enhancing drought resistance in sorghum.

Exploring the emerging role of long non-coding RNAs (lncRNAs) in plant biology: Functions, mechanisms of action, and future directions

Long non-coding RNAs (lncRNAs) are a class of RNA transcripts that surpass 200 nucleotides in length and lack discernible coding potential. LncRNAs that have been functionally characterized have pivotal functions in several plant processes, including the regulation of flowering, and development of lateral roots. It also plays a crucial role in the plant's response to abiotic stressors and exhibits vital activities in environmental adaptation. The progress in NGS (next-generation sequencing) and functional genomics technology has facilitated the discovery of lncRNA in plant species. This review is a brief explanation of lncRNA genomics, its molecular role, and the mechanism of action in plants. The review also addresses the challenges encountered in this field and highlights promising molecular and computational methodologies that can aid in the comparative and functional analysis of lncRNAs.

Advanced Biotechnological Interventions in Mitigating Drought Stress in Plants

This comprehensive article critically analyzes the advanced biotechnological strategies to mitigate plant drought stress. It encompasses an in-depth exploration of the latest developments in plant genomics, proteomics, and metabolomics, shedding light on the complex molecular mechanisms that plants employ to combat drought stress. The study also emphasizes the significant advancements in genetic engineering techniques, particularly CRISPR-Cas9 genome editing, which have revolutionized the creation of drought-resistant crop varieties. Furthermore, the article explores microbial biotechnology's pivotal role, such as plant growth-promoting rhizobacteria (PGPR) and mycorrhizae, in enhancing plant resilience against drought conditions. The integration of these cutting-edge biotechnological interventions with traditional breeding methods is presented as a holistic approach for fortifying crops against drought stress. This integration addresses immediate agricultural needs and contributes significantly to sustainable agriculture, ensuring food security in the face of escalating climate change challenges.

Transcriptomic analysis of ncRNAs and mRNAs interactions during drought stress in switchgrass

Switchgrass (Panicum virgatum L.) plays a pivotal role as a bioenergy feedstock in the production of cellulosic ethanol and contributes significantly to enhancing ecological grasslands and soil quality. The utilization of non-coding RNAs (ncRNAs) has gained momentum in deciphering the intricate genetic responses to abiotic stress in various plant species. Nevertheless, the current research landscape lacks a comprehensive exploration of the responses of diverse ncRNAs, including long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), and microRNAs (miRNAs), to drought stress in switchgrass. In this study, we employed whole transcriptome sequencing to comprehensively characterize the expression profiles of both mRNA and ncRNAs during episodes of drought stress in switchgrass. Our analysis identified a total of 12,511 mRNAs, 59 miRNAs, 38 circRNAs, and 368 lncRNAs that exhibited significant differential expression between normal and drought-treated switchgrass leaves. Notably, the majority of up-regulated mRNAs displayed pronounced enrichment within the starch and sucrose metabolism pathway, as validated through KEGG analysis. Co-expression analysis illuminated that differentially expressed (DE) lncRNAs conceivably regulated 1308 protein-coding genes in trans and 7110 protein-coding genes in cis. Furthermore, both cis- and trans-target mRNAs of DE lncRNAs exhibited enrichment in four common KEGG pathways. The intricate interplay between lncRNAs and circRNAs with miRNAs via miRNA response elements was explored within the competitive endogenous RNA (ceRNA) network framework. As a result, we constructed elaborate regulatory networks, including lncRNA-novel_miRNA480-mRNA, lncRNA-novel_miRNA304-mRNA, lncRNA/circRNA-novel_miRNA122-PvSS4, and lncRNA/circRNA-novel_miRNA14-PvSS4, and subsequently validated the functionality of the target gene, starch synthase 4 (PvSS4). Furthermore, through the overexpression of PvSS4, we ascertained its capacity to enhance drought tolerance in yeast. However, it is noteworthy that PvSS4 did not exhibit any discernible impact under salt stress conditions. These findings, as presented herein, not only contribute substantively to our understanding of ceRNA networks but also offer a basis for further investigations into their potential functions in response to drought stress in switchgrass.

Long non-coding RNA-mediated epigenetic response for abiotic stress tolerance in plants

Plants perceive environmental fluctuations as stress and confront several stresses throughout their life cycle individually or in combination. Plants have evolved their sensing and signaling mechanisms to perceive and respond to a variety of stresses. Epigenetic regulation plays a critical role in the regulation of genes, spatiotemporal expression of genes under stress conditions and imparts a stress memory to encounter future stress responses. It is quintessential to integrate our understanding of genetics and epigenetics to maintain plant fitness, achieve desired genetic gains with no trade-offs, and durable long-term stress tolerance. The long non-coding RNA >200 nts having no coding potential (or very low) play several roles in epigenetic memory, contributing to the regulation of gene expression and the maintenance of cellular identity which include chromatin remodeling, imprinting (dosage compensation), stable silencing, facilitating nuclear organization, regulation of enhancer-promoter interactions, response to environmental signals and epigenetic switching. The lncRNAs are involved in a myriad of stress responses by activation or repression of target genes and hence are potential candidates for deploying in climate-resilient breeding programs. This review puts forward the significant roles of long non-coding RNA as an epigenetic response during abiotic stresses in plants and the prospects of deploying lncRNAs for designing climate-resilient plants.

Transcriptomic Complexity of Culm Growth and Development in Different Types of Moso Bamboo

Moso bamboo is capable of both sexual and asexual reproduction during natural growth, resulting in four distinct types of culms: the bamboo shoot-culm, the seedling stem, the leptomorph rhizome, and a long-ignored culm-the outward-rhizome. Sometimes, when the outward rhizomes break through the soil, they continue to grow longitudinally and develop into a new individual. However, the roles of alternative transcription start sites (aTSS) or termination sites (aTTS) as well as alternative splicing (AS) have not been comprehensively studied for their development. To re-annotate the moso bamboo genome and identify genome-wide aTSS, aTTS, and AS in growing culms, we utilized single-molecule long-read sequencing technology. In total, 169,433 non-redundant isoforms and 14,840 new gene loci were identified. Among 1311 lncRNAs, most of which showed a positive correlation with their target mRNAs, one-third of these IncRNAs were preferentially expressed in winter bamboo shoots. In addition, the predominant AS type observed in moso bamboo was intron retention, while aTSS and aTTS events occurred more frequently than AS. Notably, most genes with AS events were also accompanied by aTSS and aTTS events. Outward rhizome growth in moso bamboo was associated with a significant increase in intron retention, possibly due to changes in the growth environment. As different types of moso bamboo culms grow and develop, a significant number of isoforms undergo changes in their conserved domains due to the regulation of aTSS, aTTS, and AS. As a result, these isoforms may play different roles than their original functions. These isoforms then performed different functions from their original roles, contributing to the transcriptomic complexity of moso bamboo. Overall, this study provided a comprehensive overview of the transcriptomic changes underlying different types of moso bamboo culm growth and development.

Comprehensive investigation of long non-coding RNAs in an endophytic fungus Calcarisporium arbuscula NRRL 3705

Long non-coding RNAs (lncRNAs) play an important role in eukaryotic cells. However, there is no report of lncRNAs in endophytic fungi Calcarisporium arbuscula. Here, in Calcarisporium arbuscula NRRL 3705, an endophytic fungus predominantly producing mycotoxins aurovertins, the genome-wide identification of lncRNAs was carried out based on RNA-Seq. Totally, 1332 lncRNAs were identified, including 1082 long intergenic noncoding RNAs, 64 long intronic noncoding RNAs and 186 long noncoding natural antisense transcripts. The average length of lncRNA and mRNA were 254 and 1102 bp, respectively. LncRNAs were shorter, with fewer exons and lower expression levels. Moreover, there were 39 up-regulated lncRNAs and 10 down-regulated lncRNAs in the ΔaurA mutant, which lacks the aurovertin biosynthetic enzyme AurA. Interestingly, expression of genes related to the metabolism of linoleic acid and methane were significantly down regulated in the ΔaurA mutant. This study enriches the endophytic fungal lncRNA database and provide a basis for further research.

Genome-Wide Analysis of Long Non-Coding RNAs Related to UV-B Radiation in the Antarctic Moss Pohlia nutans

Antarctic organisms are consistently suffering from multiple environmental pressures, especially the strong UV radiation caused by the loss of the ozone layer. The mosses and lichens dominate the vegetation of the Antarctic continent, which grow and propagate in these harsh environments. However, the molecular mechanisms and related regulatory networks of these Antarctic plants against UV-B radiation are largely unknown. Here, we used an integrated multi-omics approach to study the regulatory mechanism of long non-coding RNAs (lncRNAs) of an Antarctic moss (Pohlia nutans) in response to UV-B radiation. We identified a total of 5729 lncRNA sequences by transcriptome sequencing, including 1459 differentially expressed lncRNAs (DELs). Through functional annotation, we found that the target genes of DELs were significantly enriched in plant-pathogen interaction and the flavonoid synthesis pathway. In addition, a total of 451 metabolites were detected by metabonomic analysis, and 97 differentially change metabolites (DCMs) were found. Flavonoids account for 20% of the total significantly up-regulated metabolites. In addition, the comprehensive transcriptome and metabolome analyses revealed the co-expression pattern of DELs and DCMs of flavonoids. Our results provide insights into the regulatory network of lncRNA under UV-B radiation and the adaptation of Antarctic moss to the polar environments.

A Survey of the Transcriptomic Resources in Durum Wheat: Stress Responses, Data Integration and Exploitation

Durum wheat (Triticum turgidum subsp. durum (Desf.) Husn.) is an allotetraploid cereal crop of worldwide importance, given its use for making pasta, couscous, and bulgur. Under climate change scenarios, abiotic (e.g., high and low temperatures, salinity, drought) and biotic (mainly exemplified by fungal pathogens) stresses represent a significant limit for durum cultivation because they can severely affect yield and grain quality. The advent of next-generation sequencing technologies has brought a huge development in transcriptomic resources with many relevant datasets now available for durum wheat, at various anatomical levels, also focusing on phenological phases and environmental conditions. In this review, we cover all the transcriptomic resources generated on durum wheat to date and focus on the corresponding scientific insights gained into abiotic and biotic stress responses. We describe relevant databases, tools and approaches, including connections with other "omics" that could assist data integration for candidate gene discovery for bio-agronomical traits. The biological knowledge summarized here will ultimately help in accelerating durum wheat breeding.

Long Non-Coding RNAs of Plants in Response to Abiotic Stresses and Their Regulating Roles in Promoting Environmental Adaption

Abiotic stresses triggered by climate change and human activity cause substantial agricultural and environmental problems which hamper plant growth. Plants have evolved sophisticated mechanisms in response to abiotic stresses, such as stress perception, epigenetic modification, and regulation of transcription and translation. Over the past decade, a large body of literature has revealed the various regulatory roles of long non-coding RNAs (lncRNAs) in the plant response to abiotic stresses and their irreplaceable functions in environmental adaptation. LncRNAs are recognized as a class of ncRNAs that are longer than 200 nucleotides, influencing a variety of biological processes. In this review, we mainly focused on the recent progress of plant lncRNAs, outlining their features, evolution, and functions of plant lncRNAs in response to drought, low or high temperature, salt, and heavy metal stress. The approaches to characterize the function of lncRNAs and the mechanisms of how they regulate plant responses to abiotic stresses were further reviewed. Moreover, we discuss the accumulating discoveries regarding the biological functions of lncRNAs on plant stress memory as well. The present review provides updated information and directions for us to characterize the potential functions of lncRNAs in abiotic stresses in the future.

Genome-wide analysis of long non-coding RNAs in sugar beet (Beta vulgaris L.) under drought stress

Drought stress is one of the most severe abiotic stresses that restrict global crop production. Long non-coding RNAs (lncRNAs) have been proved to play a key role in response to drought stress. However, genome-wide identification and characterization of drought-responsive lncRNAs in sugar beet is still lacking. Thus, the present study focused on analyzing lncRNAs in sugar beet under drought stress. We identified 32017 reliable lncRNAs in sugar beet by strand-specific high-throughput sequencing. A total of 386 differentially expressed lncRNAs (DElncRNAs) were found under drought stress. The most significantly upregulated and downregulated lncRNAs were TCONS_00055787 (upregulated by more than 6000 fold) and TCONS_00038334 (downregulated by more than 18000 fold), respectively. Quantitative real-time PCR results exhibited a high concordance with RNA sequencing data, which conformed that the expression patterns of lncRNAs based on RNA sequencing were highly reliable. In addition, we predicted 2353 and 9041 transcripts that were estimated to be the cis- and trans-target genes of the drought-responsive lncRNAs. As revealed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the target genes of DElncRNAs were significantly enriched in organelle subcompartment, thylakoid, endopeptidase activity, catalytic activity, developmental process, lipid metabolic process, RNA polymerase activity, transferase activity, flavonoid biosynthesis and several other terms associated with abiotic stress tolerance. Moreover, 42 DElncRNAs were predicted as potential miRNA target mimics. LncRNAs have important effects on plant adaptation to drought conditions through the interaction with protein-encoding genes. The present study leads to greater insights into lncRNA biology and offers candidate regulators for improving the drought tolerance of sugar beet cultivars at the genetic level.

Genome-wide analysis of long non-coding RNAs under diel light exhibits role in floral development and the circadian clock in Arabidopsis thaliana

The circadian clock is regulated by signaling networks that enhance a plant's ability to coordinate internal events with the external environment. In this study, we examine the rhythmic expression of long non-coding RNAs (lncRNAs) using multiple transcriptomes of Arabidopsis thaliana in the diel light cycle and integrated this information to have a better understanding of the functions of lncRNAs in regulating the circadian clock. We identified 968, 1050, and 998 lncRNAs at 8 h light, 16 h light and 8 h dark conditions, respectively. Among these, 423, 486, and 417 lncRNAs were uniquely present at 8 h light, 16 h light, and 8 h dark, respectively, whereas 334 lncRNAs were common under the three conditions. The specificity of identified lncRNAs under different light conditions was verified using qRT-PCR. The identified lncRNAs were less GC-rich and expressed at a significantly lower level than the mRNAs of protein-coding genes. In addition, we identified enriched motifs in lncRNA transcribing regions that were associated with light-responsive genes (SORLREP and SORLIP), flower development (AGAMOUS), and circadian clock (CCA1) under all three light conditions. We identified 10 and 12 different lncRNAs targeting different miRNAs with perfect and interrupted complementarity (endogenous target mimic). These predicted lncRNA-interacting miRNAs govern the function of a set of genes involved in the developmental process, reproductive structure development, gene silencing and transcription regulation. We demonstrated that the lncRNA transcribing regions were enriched for epigenetic marks such as H3.3, H3K4me2, H3K4me3, H4K16ac, H3K36ac, H3K56ac and depleted for heterochromatic (H3K9me2 and H3K27me1) and repressive (H3K27me3) histone modifications. Further, we found that hypermethylated genomic regions negatively correlated with lncRNA transcribing regions. Overall, our study showed that lncRNAs expressed corresponding to the diel light cycle are implicated in regulating the circadian rhythm and governing the developmental stage-specific growth.

An intersectional analysis of LncRNAs and mRNAs reveals the potential therapeutic targets of Bi Zhong Xiao Decoction in collagen-induced arthritis rats

Background

Bi Zhong Xiao decoction (BZXD), a traditional Chinese herbal formula, has been used clinically for many years to treat rheumatoid arthritis (RA). Both clinical and experimental studies have revealed that BZXD is effective in treating RA, but the mechanism remains unclear. In this study, we aimed to explore the mechanism of efficacy of BZXD through transcriptomic analysis of lncRNA and mRNA.

Methods

The combination method of ultra-high performance liquid chromatography-mass spectrometry/mass spectrometry was used to assess the quality of BZXD. The efficacy of BZXD in treating collagen-induced arthritis (CIA) was evaluated by clinical assessment, weight changes, hematoxylin–eosin and safranin o-fast green staining, and Micro-CT. Arraystar rat lncRNA-mRNA chip technology was used to determine the lncRNA and mRNA expression profiles of the Control, CIA and BZXD groups, and to screen gene expression profiles related to the curative effect of BZXD. A lncRNA-mRNA co-expression network was constructed for the therapeutic efficacy genes. Through GO function and KEGG pathway enrichment analysis, the biological functions and signaling pathways of therapeutic efficacy genes were determined. Based on fold change and functional annotation, key differentially expressed lncRNAs and mRNAs were selected for reverse transcription-quantitative polymerase chain reaction (RT-qPCR) validation. The functions of lncRNAs targeting mRNAs were verified in vitro.

Results

We demonstrated that BZXD could effectively reverse bone erosion. After BZXD treatment, up to 33 lncRNAs and 107 mRNAs differentially expressed genes were reversely regulated by BZXD. These differentially expressed lncRNAs are mainly involved in the biological process of the immune response and are closely related to the ECM-receptor interaction, MAPK signaling pathway, Focal adhesion, Ras signaling pathway, Antigen processing and presentation, and Chemokine signaling pathway. We identified four lncRNAs (uc.361−, ENSRNOT00000092834, ENSRNOT00000089244, ENSRNOT00000084631) and three mRNAs (Acvr2a, Cbx2, Morc4) as potential therapeutic targets for BZXD and their microarray data consistent with the RT-qPCR. In vitro experiments confirmed that silencing the lncRNAs ENSRNOT00000092834 and ENSRNOT00000084631 reversed the expression of target mRNAs.

Conclusions

This study elucidates the possible mechanism of BZXD reversing bone erosion in CIA rats from the perspective of lncRNA and mRNA. To provide a basis and direction for further exploration of the mechanism of BZXD in treating RA.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13020-022-00670-z.

Genomic sequencing of Thinopyrum elongatum chromosome arm 7EL, carrying fusarium head blight resistance, and characterization of its impact on the transcriptome of the introgressed line CS-7EL

Background

The tall wheatgrass species Thinopyrum elongatum carries a strong fusarium head blight (FHB) resistance locus located on the long arm of chromosome 7 (7EL) as well as resistance to leaf and stem rusts, all diseases with a significant impact on wheat production. Towards understanding the contribution of Th. elongatum 7EL to improvement of disease resistance in wheat, the genomic sequence of the 7EL fragment present in the wheat Chinese Spring (CS) telosomic addition line CS-7EL was determined and the contribution and impact of 7EL on the rachis transcriptome during FHB infection was compared between CS and CS-7EL.

Results

We assembled the Th. elongatum 7EL chromosome arm using a reference-guided approach. Combining this assembly with the available reference sequence for CS hexaploid wheat provided a reliable reference for interrogating the transcriptomic differences in response to infection conferred by the 7EL fragment. Comparison of the transcriptomes of rachis tissues from CS and CS-7EL showed expression of Th. elongatum transcripts as well as modulation of wheat transcript expression profiles in the CS-7EL line. Expression profiles at 4 days after infection with Fusarium graminearum, the causal agent of FHB, showed an increased in expression of genes associated with an effective defense response, in particular glucan endo-1,3-beta-glucosidases and chitinases, in the FHB-resistant line CS-7EL while there was a larger increase in differential expression for genes associated with the level of fungal infection in the FHB-susceptible line CS. One hundred and seven 7EL transcripts were expressed in the smallest 7EL region defined to carry FHB resistance.

Conclusion

7EL contributed to CS-7EL transcriptome by direct expression and through alteration of wheat transcript profiles. FHB resistance in CS-7EL was associated with transcriptome changes suggesting a more effective defense response. A list of candidate genes for the FHB resistance locus on 7EL has been established.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08433-8.

Global identification and integrated analysis of heat-responsive long non-coding RNAs in contrasting rice cultivars

Identified 2743 rice lncRNAs LncRNAs in response to heat stress Function prediction of HRLs Network among HRLs, genes and miRNAs co-localization of HRLs with QTLs Significant motifs in HRL sequences Long non-coding RNAs (lncRNAs) play vital roles in plant responses to environmental challenges. A better understanding of the gene regulation mediated by lncRNAs and their systematic identification would provide great benefits for modern agriculture. In this study, we performed strand-specific RNA sequencing for two rice varieties, heat-tolerant ZS97B and heat-susceptible SYD2 under heat stress. In total, 2743 putative lncRNAs were identified, and their expression profiles in response to heat treatments were established. We identified 231 differentially expressed lncRNAs (DELs) under heat stress, including 31 DELs common to both varieties and 103 and 97 specific to ZS97B and SYD2, respectively, all defined as heat-responsive lncRNAs (HRLs). The target-coding genes of HRLs were predicted, and GO and KEGG annotations of HRL targets revealed functions in which HRLs might be involved. The interaction network between HRLs, target genes and relevant miRNAs was constructed. The HRLs and their targets were compared with publicly available QTLs for rice seedling growth under heat stimulus. Ten HRLs and twelve target genes were linked with five heat stress-relevant QTLs. Sequence analysis revealed several motifs significantly enriched within the 231 HRL sequences. Our findings provide a valuable resource for further characterization of lncRNAs in terms of heat response and plant heat tolerance improvement.

Novel insights into water-deficit-responsive mRNAs and lncRNAs during fiber development in Gossypium hirsutum

BACKGROUND

The fiber yield and quality of cotton are greatly and periodically affected by water deficit. However, the molecular mechanism of the water deficit response in cotton fiber cells has not been fully elucidated.

RESULTS

In this study, water deficit caused a significant reduction in fiber length, strength, and elongation rate but a dramatic increase in micronaire value. To explore genome-wide transcriptional changes, fibers from cotton plants subjected to water deficit (WD) and normal irrigation (NI) during fiber development were analyzed by transcriptome sequencing. Analysis showed that 3427 mRNAs and 1021 long noncoding RNAs (lncRNAs) from fibers were differentially expressed between WD and NI plants. The maximum number of differentially expressed genes (DEGs) and lncRNAs (DERs) was identified in fibers at the secondary cell wall biosynthesis stage, suggesting that this is a critical period in response to water deficit. Twelve genes in cotton fiber were differentially and persistently expressed at ≥ five time points, suggesting that these genes are involved in both fiber development and the water-deficit response and could potentially be used in breeding to improve cotton resistance to drought stress. A total of 540 DEGs were predicted to be potentially regulated by DERs by analysis of coexpression and genomic colocation, accounting for approximately 15.76% of all DEGs. Four DERs, potentially acting as target mimics for microRNAs (miRNAs), indirectly regulated their corresponding DEGs in response to water deficit.

CONCLUSIONS

This work provides a comprehensive transcriptome analysis of fiber cells and a set of protein-coding genes and lncRNAs implicated in the cotton response to water deficit, significantly affecting fiber quality during the fiber development stage.

Non-Coding RNAs in Response to Drought Stress

Drought stress causes changes in the morphological, physiological, biochemical and molecular characteristics of plants. The response to drought in different plants may vary from avoidance, tolerance and escape to recovery from stress. This response is genetically programmed and regulated in a very complex yet synchronized manner. The crucial genetic regulations mediated by non-coding RNAs (ncRNAs) have emerged as game-changers in modulating the plant responses to drought and other abiotic stresses. The ncRNAs interact with their targets to form potentially subtle regulatory networks that control multiple genes to determine the overall response of plants. Many long and small drought-responsive ncRNAs have been identified and characterized in different plant varieties. The miRNA-based research is better documented, while lncRNA and transposon-derived RNAs are relatively new, and their cellular role is beginning to be understood. In this review, we have compiled the information on the categorization of non-coding RNAs based on their biogenesis and function. We also discuss the available literature on the role of long and small non-coding RNAs in mitigating drought stress in plants.

Epigenetic regulation of salinity stress responses in cereals

Cereals are important crops and are exposed to various types of environmental stresses that affect the overall growth and yield. Among the various abiotic stresses, salt stress is a major environmental factor that influences the genetic, physiological, and biochemical responses of cereal crops. Epigenetic regulation which includes DNA methylation, histone modification, and chromatin remodelling plays an important role in salt stress tolerance. Recent studies in rice genomics have highlighted that the epigenetic changes are heritable and therefore can be considered as molecular signatures. An epigenetic mechanism under salinity induces phenotypic responses involving modulations in gene expression. Association between histone modification and altered DNA methylation patterns and differential gene expression has been evidenced for salt sensitivity in rice and other cereal crops. In addition, epigenetics also creates stress memory that helps the plant to better combat future stress exposure. In the present review, we have discussed epigenetic influences in stress tolerance, adaptation, and evolution processes. Understanding the epigenetic regulation of salinity could help for designing salt-tolerant varieties leading to improved crop productivity.

Expression Patterns and Regulation of Non-Coding RNAs during Synthesis of Cellulose in Eucalyptus grandis Hill

Cellulose, an essential structural component in the plant cell wall and a renewable biomass resource, plays a significant role in nature. Eucalyptus’s excellent timber tree species (including Eucalyptus grandis Hill) provide many raw materials for the paper and wood industries. The synthesis of cellulose is a very complex process involving multiple genes and regulated by various biological networks. However, research on regulating associated genes and non-coding RNAs during cellulose synthesis in E. grandis remains lacking. In this study, the wood anatomical characteristics and chemical indexes of E. grandis were analyzed by taking three different parts (diameter at breast height (DBH), middle and upper part of the trunk) from the main stem of E. grandis as raw materials. The role of non-coding RNAs (Long non-coding RNA, lncRNA; Micro RNA, miRNA; Circle RNA, circRNA) on regulating candidate genes was presented, and the network map of ceRNA (Competing endogenous RNA) regulation during wood cellulose biosynthesis of E. grandis was constructed. The transcriptome sequencing of nine samples obtained from the trunk of the immature xylem in E. grandis at DBH, middle and upper parts had a 95.81 G clean reading, 57,480 transcripts, 7365 lncRNAs, and 5180 circRNAs. Each sample had 172–306 known miRNAs and 1644–3508 new miRNAs. A total of 190 DE-lncRNAs (Differentially expressed long non-coding RNAs), 174 DE-miRNAs (Differentially expressed micro RNAs), and 270 DE-circRNAs (Differentially expressed circle RNAs) were obtained by comparing transcript expression levels. Four lncRNAs and nine miRNAs were screened out, and the ceRNA regulatory network was constructed. LncRNA1 and lncRNA4 regulated the genes responsible for cellulose synthesis in E. grandis, which were overexpressed in 84K (Populus Alba × Populus glandulosa) poplar. The cellulose and lignin content in lncRNA4-oe were significantly higher than wild type 84K poplar and lncRNA1-oe. The average plant height, middle and basal part of the stem diameter in lncRNA4-oe were significantly higher than the wild type. However, there was no significant difference between the growth of lncRNA1-oe and the wild type. Further studies are warranted to explore the molecular regulatory mechanism of cellulose biosynthesis in Eucalyptus species.

Identification and characterization of early Fusarium wilt responsive mRNAs and long non-coding RNAs in banana root using high-throughput sequencing

Fusarium wilt disease, caused by Fusarium oxysporum f.sp. cubense (Foc), has been recognized as the most devastating disease to banana. The regulatory role of long non-coding RNAs (lncRNAs) in plant defense has been verified in many plant species. However, the understanding of their role during early FocTR4 (Foc tropical race 4) infection stage is very limited. In this study, lncRNA sequencing was used to reveal banana root transcriptome profile changes during early FocTR4 infection stages. Quantitative real time PCR (qRT-PCR) was performed to confirm the expression of eight differentially expressed (DE) lncRNAs (DELs) and their predicted target genes (DETs), and three DE genes (DEGs). Totally, 12,109 lncRNAs, 36,519 mRNAs and 2642 novel genes were obtained, of which 1398 (including 78 DELs, 1220 DE known genes and 100 DE novel genes) were identified as FocTR4 responsive DE transcripts. Gene function analysis revealed that most DEGs were involved in biosynthesis of secondary metabolites, plant–pathogen interaction, plant hormone signal transduction, phenylalanine metabolism, phenylpropanoid biosynthesis, alpha-linolenic acid metabolism and so on. Coincidently, many DETs have been identified as DEGs in previous transcriptome studies. Moreover, many DETs were found to be involved in ribosome, oxidative phosphorylation, lipoic acid metabolism, ubiquitin mediated proteolysis, N-glycan biosynthesis, protein processing in endoplasmic reticulum and DNA damage response pathways. QRT-PCR result showed the expression patterns of the selected transcripts were mostly consistent with our lncRNA sequencing data. Our present study showed the regulatory role of lncRNAs on known biotic and abiotic stress responsive genes and some new-found FocTR4 responsive genes, which can provide new insights into FocTR4-induced changes in the banana root transcriptome during the early pathogen infection stage.

Recognizing the hidden half in wheat: root system attributes associated with drought tolerance

Improving drought tolerance in wheat is crucial for maintaining productivity and food security. Roots are responsible for the uptake of water from soil, and a number of root traits are associated with drought tolerance. Studies have revealed many quantitative trait loci and genes controlling root development in plants. However, the genetic dissection of root traits in response to drought in wheat is still unclear. Here, we review crop root traits associated with drought, key genes governing root development in plants, and quantitative trait loci and genes regulating root system architecture under water-limited conditions in wheat. Deep roots, optimal root length density and xylem diameter, and increased root surface area are traits contributing to drought tolerance. In view of the diverse environments in which wheat is grown, the balance among root and shoot traits, as well as individual and population performance, are discussed. The known functions of key genes provide information for the genetic dissection of root development of wheat in a wide range of conditions, and will be beneficial for molecular marker development, marker-assisted selection, and genetic improvement in breeding for drought tolerance.

Long non-coding RNA exploration for mesenchymal stem cell characterisation

BACKGROUND

The development of RNA sequencing (RNAseq) and the corresponding emergence of public datasets have created new avenues of transcriptional marker search. The long non-coding RNAs (lncRNAs) constitute an emerging class of transcripts with a potential for high tissue specificity and function. Therefore, we tested the biomarker potential of lncRNAs on Mesenchymal Stem Cells (MSCs), a complex type of adult multipotent stem cells of diverse tissue origins, that is frequently used in clinics but which is lacking extensive characterization.

RESULTS

We developed a dedicated bioinformatics pipeline for the purpose of building a cell-specific catalogue of unannotated lncRNAs. The pipeline performs ab initio transcript identification, pseudoalignment and uses new methodologies such as a specific k-mer approach for naive quantification of expression in numerous RNAseq data. We next applied it on MSCs, and our pipeline was able to highlight novel lncRNAs with high cell specificity. Furthermore, with original and efficient approaches for functional prediction, we demonstrated that each candidate represents one specific state of MSCs biology.

CONCLUSIONS

We showed that our approach can be employed to harness lncRNAs as cell markers. More specifically, our results suggest different candidates as potential actors in MSCs biology and propose promising directions for future experimental investigations.

Interplay between miRNAs and lncRNAs: Mode of action and biological roles in plant development and stress adaptation

Plants employ sophisticated mechanisms to control developmental processes and to cope with environmental changes at transcriptional and post-transcriptional levels. MicroRNAs (miRNAs) and long noncoding RNAs (lncRNAs), two classes of endogenous noncoding RNAs, are key regulators of gene expression in plants. Recent studies have identified the interplay between miRNAs and lncRNAs as a novel regulatory layer of gene expression in plants. On one hand, miRNAs target lncRNAs for the production of phased small interfering RNAs (phasiRNAs). On the other hand, lncRNAs serve as origin of miRNAs or regulate the accumulation or activity of miRNAs at transcription and post-transcriptional levels. Theses lncRNA-miRNA interplays are crucial for plant development, physiology and responses to biotic and abiotic stresses. In this review, we summarize recent advances in the biological roles, interaction mechanisms and computational predication methods of the interplay between miRNAs and lncRNAs in plants.

From Trash to Luxury: The Potential Role of Plant LncRNA in DNA Methylation During Abiotic Stress

Remarkable progress has been made in elucidating important roles of plant non-coding RNAs. Among these RNAs, long noncoding RNAs (lncRNAs) have gained widespread attention, especially their role in plant environmental stress responses. LncRNAs act at different levels of gene expression regulation, and one of these mechanisms is by recruitment of DNA methyltransferases or demethylases to regulate the target gene transcription. In this mini-review, we highlight the function of lncRNAs, including their potential role in RNA-directed DNA Methylation (RdDM) silencing pathway and their potential function under abiotic stresses conditions. Moreover, we also present and discuss studies of lncRNAs in crops. Finally, we propose a path outlook for future research that may be important for plant breeding.

Long non-coding RNAs: emerging players regulating plant abiotic stress response and adaptation

BACKGROUND

The immobile nature of plants means that they can be frequently confronted by various biotic and abiotic stresses during their lifecycle. Among the various abiotic stresses, water stress, temperature extremities, salinity, and heavy metal toxicity are the major abiotic stresses challenging overall plant growth. Plants have evolved complex molecular mechanisms to adapt under the given abiotic stresses. Long non-coding RNAs (lncRNAs)-a diverse class of RNAs that contain > 200 nucleotides(nt)-play an essential role in plant adaptation to various abiotic stresses.

RESULTS

LncRNAs play a significant role as 'biological regulators' for various developmental processes and biotic and abiotic stress responses in animals and plants at the transcription, post-transcription, and epigenetic level, targeting various stress-responsive mRNAs, regulatory gene(s) encoding transcription factors, and numerous microRNAs (miRNAs) that regulate the expression of different genes. However, the mechanistic role of lncRNAs at the molecular level, and possible target gene(s) contributing to plant abiotic stress response and adaptation, remain largely unknown. Here, we review various types of lncRNAs found in different plant species, with a focus on understanding the complex molecular mechanisms that contribute to abiotic stress tolerance in plants. We start by discussing the biogenesis, type and function, phylogenetic relationships, and sequence conservation of lncRNAs. Next, we review the role of lncRNAs controlling various abiotic stresses, including drought, heat, cold, heavy metal toxicity, and nutrient deficiency, with relevant examples from various plant species. Lastly, we briefly discuss the various lncRNA databases and the role of bioinformatics for predicting the structural and functional annotation of novel lncRNAs.

CONCLUSIONS

Understanding the intricate molecular mechanisms of stress-responsive lncRNAs is in its infancy. The availability of a comprehensive atlas of lncRNAs across whole genomes in crop plants, coupled with a comprehensive understanding of the complex molecular mechanisms that regulate various abiotic stress responses, will enable us to use lncRNAs as potential biomarkers for tailoring abiotic stress-tolerant plants in the future.

Uncloaking lncRNA-meditated gene expression as a potential regulator of CMS in cotton (Gossypium hirsutum L.)

Cytoplasmic male sterility is a well-proven mechanism for cotton hybrid production. Long non-coding RNAs belong to a class of transcriptional regulators that function in multiple biological processes. The cDNA libraries from the flower buds of the cotton CGMS, it's restorer (Rf) and maintainer lines were sequenced using high throughput NGS technique. A total of 1531 lncRNAs showed significant differential expression patterns between these three lines. Functional analysis of the co-expression network of lncRNA-mRNA using gene ontology vouchsafes that, lncRNAs play a crucial role in cytoplasmic male sterility and fertility restoration through pollen development, INO80 complex, development of anther wall tapetum, chromatin remodeling, and histone modification. Additionally, 94 lncRNAs were identified as putative precursors of 49 miRNAs. qRT-PCR affirms the concordance of expression pattern to RNA-seq data. These findings divulge the lncRNA driven miRNA-mediated regulation of gene expression profiling superintended for a better understanding of the CMS mechanisms of cotton.

High throughput deep sequencing elucidates the important role of lncRNAs in Foxtail millet response to herbicides

Long non-coding RNAs (lncRNAs) play an important function in plant growth and development as well as response to stresses. However, little information was known in foxtail millet; no study was reported on lncRNAs in plant response to herbicide treatment. In this study, by using deep sequencing and advanced bioinformatic analysis, a total of 2547 lncRNAs were identified, including 787 known and 1760 novel lncRNAs. These lncRNAs are distributed across all 9 chromosomes, and the majority were located in the intergenic region with 1-2 exons. These lncRNAs were differentially expressed between different genotypes under different herbicide treatments. lncRNAs regulate plant growth and development as well as response to herbicide treatments through targeting protein-coding genes that directly relate to chemical metabolism and defense system. Multiple potential target genes and lncRNA-mRNA-miRNA gene networks were discovered. These results elucidate the potential roles of lncRNAs in plant response to herbicides.

LncRNA improves cold resistance of winter wheat by interacting with miR398

One of the important functions of long non-coding RNA (lncRNA) is to be competing endogenous RNAs (ceRNAs). As miR398 is reported to respond to different stressors, it is necessary to explore its relationship with lncRNA in the cold resistance mechanism of winter wheat. Tae-miR398-precursor sequence was isolated from the winter wheat (Triticum aestivum). RLM-RACE verified that tae-miR398 cleaved its target CSD1. Quantitative detection at 5°C, -10°C and -25°C showed that the expression of tae-miR398 decreased in response to low temperatures, whereas CSD1 showed an opposite expression pattern. LncR9A, lncR117 and lncR616 were predicted and verified to interact with miR398. tae-miR398 and three lncRNAs were transferred into Arabidopsis thaliana respectively. The lncR9A were transferred into Brachypodium distachyom. Transgenic plants were cultivated at -8°C and assessed for the expression of malondialdehyde, chlorophyll, superoxide dismutase and miR398-lncRNA-target mRNA. The results demonstrate that tae-miR398 regulates low temperature tolerance by downregulating its target, CSD1. lncRNA regulates the expression of CSD1 indirectly by competitively binding miR398, which, in turn, affects the resistance of Dn1 to cold. miR398-regulation triggers a regulatory loop that is critical to cold stress tolerance in wheat. Our findings offer an improved strategy to crop plants with enhanced stress tolerance.

A genome-wide identification, characterization and functional analysis of salt-related long non-coding RNAs in non-model plant Pistacia vera L. using transcriptome high throughput sequencing

Long non-coding RNAs (lncRNAs) play crucial roles in regulating gene expression in response to plant stresses. Given the importance regulatory roles of lncRNAs, providing methods for predicting the function of these molecules, especially in non-model plants, is strongly demanded by researchers. Here, we constructed a reference sequence for lncRNAs in P. vera (Pistacia vera L.) with 53220 transcripts. In total, we identified 1909 and 2802 salt responsive lncRNAs in Ghazvini, a salt tolerant cultivar, after 6 and 24 h salt treatment, respectively and 1820 lncRNAs in Sarakhs, a salt sensitive cultivar, after 6 h salt treatment. Functional analysis of these lncRNAs by several hybrid methods, revealed that salt responsive NAT-related lncRNAs associated with transcription factors, CERK1, LEA, Laccase genes and several genes involved in the hormone signaling pathways. Moreover, gene ontology (GO) enrichment analysis of salt responsive target genes related to top five selected lncRNAs showed their involvement in the regulation of ATPase, cation transporter, kinase and UDP-glycosyltransferases genes. Quantitative real-time PCR (qRT-PCR) experiment results of lncRNAs, pre-miRNAs and mature miRNAs were in accordance with our RNA-seq analysis. In the present study, a comparative analysis of differentially expressed lncRNAs and microRNA precursors between salt tolerant and sensitive pistachio cultivars provides valuable knowledge on gene expression regulation under salt stress condition.

Insight into the Role of Epigenetic Processes in Abiotic and Biotic Stress Response in Wheat and Barley

Environmental stresses such as salinity, drought, heat, freezing, heavy metal and even pathogen infections seriously threaten the growth and yield of important cereal crops including wheat and barley. There is growing evidence indicating that plants employ sophisticated epigenetic mechanisms to fine-tune their responses to environmental stresses. Here, we provide an overview of recent developments in understanding the epigenetic processes and elements—such as DNA methylation, histone modification, chromatin remodeling, and non-coding RNAs—involved in plant responses to abiotic and biotic stresses in wheat and barley. Potentials of exploiting epigenetic variation for the improvement of wheat and barley are discussed.

Maize ZmBES1/BZR1-5 Decreases ABA Sensitivity and Confers Tolerance to Osmotic Stress in Transgenic Arabidopsis

The BRI1-EMS suppressor 1 (BES1)/brassinazole-resistant 1 (BZR1) transcription factors, key components in the brassinosteroid signaling pathway, play pivotal roles in plant growth and development. However, the function of BES1/BZR1 in crops during stress response remains poorly understood. In the present study, we characterized ZmBES1/BZR1-5 from maize, which was localized to the nucleus and was responsive to abscisic acid (ABA), salt and drought stresses. Heterologous expression of ZmBES1/BZR1-5 in transgenic Arabidopsis resulted in decreased ABA sensitivity, facilitated shoot growth and root development, and enhanced salt and drought tolerance with lower malondialdehyde (MDA) content and relative electrolyte leakage (REL) under osmotic stress. The RNA sequencing (RNA-seq) analysis revealed that 84 common differentially expressed genes (DEGs) were regulated by ZmBES1/BZR1-5 in transgenic Arabidopsis. Subsequently, gene ontology and KEGG pathway enrichment analyses showed that the DEGs were enriched in response to stress, secondary metabolism and metabolic pathways. Furthermore, 30 DEGs were assigned to stress response and possessed 2-15 E-box elements in their promoters, which could be potentially recognized and bound by ZmBES1/BZR1-5. Taken together, our results reveal that the ZmBES1/BZR1-5 transcription factor positively regulates salt and drought tolerance by binding to E-box to induce the expression of downstream stress-related genes. Therefore, our study contributes to the better understanding of BES1/BZR1 function in the stress response of plants.

Long Non-coding RNA in Plants in the Era of Reference Sequences

The discovery of non-coding RNAs (ncRNAs), and the subsequent elucidation of their functional roles, was largely delayed due to the misidentification of non-protein-coding parts of DNA as "junk DNA," which forced ncRNAs into the shadows of their protein-coding counterparts. However, over the past decade, insight into the important regulatory roles of ncRNAs has led to rapid progress in their identification and characterization. Of the different types of ncRNAs, long non-coding RNAs (lncRNAs), has attracted considerable attention due to their mRNA-like structures and gene regulatory functions in plant stress responses. While RNA sequencing has been commonly used for mining lncRNAs, a lack of widespread conservation at the sequence level in addition to relatively low and highly tissue-specific expression patterns challenges high-throughput in silico identification approaches. The complex folding characteristics of lncRNA molecules also complicate target predictions, as the knowledge about the interaction interfaces between lncRNAs and potential targets is insufficient. Progress in characterizing lncRNAs and their targets from different species may hold the key to efficient identification of this class of ncRNAs from transcriptomic and potentially genomic resources. In wheat and barley, two of the most important crops, the knowledge about lncRNAs is very limited. However, recently published high-quality genomes of these crops are considered as promising resources for the identification of not only lncRNAs, but any class of molecules. Considering the increasing demand for food, these resources should be used efficiently to discover molecular mechanisms lying behind development and a/biotic stress responses. As our understanding of lncRNAs expands, interactions among ncRNA classes, as well as interactions with the coding sequences, will likely define novel functional networks that may be modulated for crop improvement.

Genome-Wide Identification of lncRNAs Involved in Fertility Transition in the Photo-Thermosensitive Genic Male Sterile Rice Line Wuxiang S

Long non-coding RNAs (lncRNAs) act as universal regulators of various biological processes, but no genome-wide screening of lncRNAs involved in the fertility transition of the photo-thermosensitive genic male sterile (PTGMS) rice line has been reported. Here, we performed strand-specific RNA sequencing at three developmental stages of a novel PTGMS line Wuxiang S (WXS). A total of 3,948 lncRNAs were identified; 622 of these were detected as differentially expressed lncRNAs (DE-lncRNAs) between male-sterile WXS (WXS-S) and male-fertile WXS (WXS-F). A large proportion of lncRNAs differentially expressed at the stage of pollen mother cells meiosis, suggested that it may be the most critical stage for fertility transition of WXS. Furthermore, functional annotation of the cis- and trans- targets of DE-lncRNAs showed that 150 targets corresponding to 141 DE-lncRNAs were identified as involved in anther and pollen development. Moreover, computational analysis predicted 97 lncRNAs as precursors for 72 miRNAs, and 94 DE-lncRNAs as potential endogenous target mimics (eTMs) for 150 miRNAs. Finally, using the dual luciferase reporter assays, we demonstrated that two lncRNAs act as eTMs to regulate the expression of the SPL and GRF genes by competing for the shared osa-miR156 and osa-miR396, respectively. These genomic characteristics, differential expression, and interaction of lncRNAs with miRNAs and mRNAs contribute to our understanding of the roles of lncRNAs during the fertility transition in PTGMS rice lines.

Transcriptomic and evolutionary analysis of the mechanisms by which P. argentatum, a rubber producing perennial, responds to drought

BACKGROUND

Guayule (Parthenium argentatum Gray) is a drought tolerant, rubber producing perennial shrub native to northern Mexico and the US Southwest. Hevea brasiliensis, currently the world's only source of natural rubber, is grown as a monoculture, leaving it vulnerable to both biotic and abiotic stressors. Isolation of rubber from guayule occurs by mechanical harvesting of the entire plant. It has been reported that environmental conditions leading up to harvest have a profound impact on rubber yield. The link between rubber biosynthesis and drought, a common environmental condition in guayule's native habitat, is currently unclear.

RESULTS

We took a transcriptomic and comparative genomic approach to determine how drought impacts rubber biosynthesis in guayule. We compared transcriptional profiles of stem tissue, the location of guayule rubber biosynthesis, collected from field-grown plants subjected to water-deficit (drought) and well-watered (control) conditions. Plants subjected to the imposed drought conditions displayed an increase in production of transcripts associated with defense responses and water homeostasis, and a decrease in transcripts associated with rubber biosynthesis. An evolutionary and comparative analysis of stress-response transcripts suggests that more anciently duplicated transcripts shared among the Asteraceae, rather than recently derived duplicates, are contributing to the drought response observed in guayule. In addition, we identified several deeply conserved long non-coding RNAs (lncRNAs) containing microRNA binding motifs. One lncRNA in particular, with origins at the base of Asteraceae, may be regulating the vegetative to reproductive transition observed in water-stressed guayule by acting as a miRNA sponge for miR166.

CONCLUSIONS

These data represent the first genomic analyses of how guayule responds to drought like conditions in agricultural production settings. We identified an inverse relationship between stress-responsive transcripts and those associated with precursor pathways to rubber biosynthesis suggesting a physiological trade-off between maintaining homeostasis and plant productivity. We also identify a number of regulators of abiotic responses, including transcription factors and lncRNAs, that are strong candidates for future projects aimed at modulating rubber biosynthesis under water-limiting conditions common to guayules' native production environment.

Transcriptome Analysis Implicates Involvement of Long Noncoding RNAs in Cytoplasmic Male Sterility and Fertility Restoration in Cotton

The cytoplasmic male sterility (CMS)/restorer-of-fertility system is an important tool to exploit heterosis during commercially hybrid seed production. The importance of long noncoding RNAs (lncRNAs) in plant development is recognized, but few analyses of lncRNAs during anther development of three-line hybrid cotton (CMS-D2 line A, maintainer line B, restorer-of-fertility line R) have been reported. Here, we performed transcriptome sequencing during anther development in three-line hybrid cotton. A total of 80,695 lncRNAs were identified, in which 43,347 and 44,739 lncRNAs were differentially expressed in A–B and A–R comparisons, respectively. These lncRNAs represent functional candidates involved in CMS and fertility restoration. GO analysis indicated that cellular hormone metabolic processes and oxidation–reduction reaction processes might be involved in CMS, and cellular component morphogenesis and small molecular biosynthetic processes might participate in fertility restoration. Additionally, 63 lncRNAs were identified as putative precursors of 35 miRNAs, and quantitative reverse transcription polymerase chain reaction (qRT-PCR) showed a similar expression pattern to RNA-seq data. Furthermore, construction of lncRNA regulatory networks indicated that several miRNA–lncRNA–mRNA networks might be involved in CMS and fertility restoration. Our findings provide systematic identification of lncRNAs during anther development and lays a solid foundation for the regulatory mechanisms and utilization in hybrid cotton breeding.

RNAseq analysis reveals drought-responsive molecular pathways with candidate genes and putative molecular markers in root tissue of wheat

Drought is one of the major impediments in wheat productivity. Traditional breeding and marker assisted QTL introgression had limited success. Available wheat genomic and RNA-seq data can decipher novel drought tolerance mechanisms with putative candidate gene and marker discovery. Drought is first sensed by root tissue but limited information is available about how roots respond to drought stress. In this view, two contrasting genotypes, namely, NI5439 41 (drought tolerant) and WL711 (drought susceptible) were used to generate ~78.2 GB data for the responses of wheat roots to drought. A total of 45139 DEGs, 13820 TF, 288 miRNAs, 640 pathways and 435829 putative markers were obtained. Study reveals use of such data in QTL to QTN refinement by analysis on two model drought-responsive QTLs on chromosome 3B in wheat roots possessing 18 differentially regulated genes with 190 sequence variants (173 SNPs and 17 InDels). Gene regulatory networks showed 69 hub-genes integrating ABA dependent and independent pathways controlling sensing of drought, root growth, uptake regulation, purine metabolism, thiamine metabolism and antibiotics pathways, stomatal closure and senescence. Eleven SSR markers were validated in a panel of 18 diverse wheat varieties. For effective future use of findings, web genomic resources were developed. We report RNA-Seq approach on wheat roots describing the drought response mechanisms under field drought conditions along with genomic resources, warranted in endeavour of wheat productivity.

Long noncoding RNAs and miRNAs regulating terpene and tartaric acid biosynthesis in rose-scented geranium

This study aimed to explore the noncoding RNAs, which have emerged as key regulatory molecules in biological processes, in rose-scented geranium. We analyzed RNA-seq data revealing 26 784 long noncoding RNAs (lncRNAs) and 871 miRNAs in rose-scented geranium. A total of 466 lncRNAs were annotated using different plant lncRNA public databases. Furthermore, 372 lncRNAs and 99 miRNAs were detected that target terpene and tartarate biosynthetic pathways. An interactome, comprising of lncRNAs, miRNAs, and mRNAs, was constructed that represents a noncoding RNA regulatory network of the target mRNAs. Real-time quantitative PCR expression validation was done for selected lncRNAs involved in the regulation of terpene and tartaric acid pathways. This study provides the first insights into the regulatory functioning of noncoding RNAs in rose-scented geranium.

Recent advances in the characterization of plant transcriptomes in response to drought, salinity, heat, and cold stress

Despite recent advancements in plant molecular biology and biotechnology, providing food security for an increasing world population remains a challenge. Drought (water scarcity), salinity, heat, and cold stress are considered major limiting factors that affect crop production both qualitatively and quantitatively. Therefore, the development of cost-effective and environmentally friendly strategies will be needed to resolve these agricultural problems. This will require a comprehensive understanding of transcriptomic alterations that occur in plants in response to varying levels of environmental stresses, singly and in combination. Here, we briefly discuss the current status and future challenges in plant research related to understanding transcriptional changes that occur in response to drought, salinity, heat, and cold stress.

Stress-responsive miRNAome of Glycine max (L.) Merrill: molecular insights and way forward

MAIN CONCLUSION

Analysis of stress-associated miRNAs of Glycine max (L.) Merrill reveals wider ramifications of small RNA-mediated (conserved and legume-specific miRNAs) gene regulatory foot prints in molecular adaptive responses. MicroRNAs (miRNAs) are indispensable components of gene regulatory mechanism of plants. Soybean is a crop of immense commercial potential grown worldwide for its edible oil and soy meal. Intensive research efforts, using the next generation sequencing and bioinformatics techniques, have led to the identification and characterization of numerous small RNAs, especially microRNAs (miRNAs), in soybean. Furthermore, studies have unequivocally demonstrated the significance of miRNAs during the developmental processes and various stresses in soybean. In this review, we summarize the current state of understanding of miRNA-based abiotic and biotic stress responses in soybean. In addition, the molecular insights gained from the stress-related soybean miRNAs have been compared to the miRNAs of other crops, especially legumes, and the core commonalities have been highlighted, though differences among them were not ignored. Nature of response of soybean-derived conserved miRNAs during various stresses was also analyzed to gain deeper insights regarding sRNAome-based defense responses. This review further provides way forward in legume small RNA transcriptomics based on the adaptive responses of soybean and other legume-derived miRNAs.

Revealing the dominant long noncoding RNAs responding to the infection with Colletotrichum gloeosporioides in Hevea brasiliensis

Background

Rubber tree (Hevea brasiliensis) acts as an important tropic economic crop and rubber tree anthracnose, mainly caused by Colletotrichum gloeosporioides, is one of the most common fungal disease, which leads to serious loss of rubber production. Therefore, the investigation on disease resistance is of great worldwide significance. In the past decades, substantial progress has been made on coding gene families related with plant disease resistance. However, in rubber tree, whether the disease resistance mechanism involves noncoding RNAs, especially long noncoding RNAs (lncRNAs), still remains poorly understood.

Results

Here, we modeled the development of H. brasiliensis leaf samples inoculated with C. gloeosporioides at divergent stages, explored to identify the expressed ncRNAs by RNA-seq, and investigated the dominant lncRNAs responding to the infection, through constructing a co-expressed network systematically. On the dominant lncRNAs, we explored the potential functional role of lncRNA11254 recruiting the transcription factor, and that lncRNA11041 and lncRNA11205 probably stimulate the accumulation of corresponding disease responsive miRNAs, and further modulate the expressions of target genes, accompanying with experimental examination.

Conclusions

Take together, computational analyses in silico and experimental evidences in our research collectively revealed the responsive roles of dominant lncRNAs to the pathogen. The results will provide new perspectives to unveil the plant disease resistance mechanisms, and will presumably provide a new theoretical basis and candidate prognostic markers for the optimization and innovation of genetic breeding for rubber tree.

Reviewers

This article was reviewed by Ryan McGinty and Roland Huber.

Electronic supplementary material

The online version of this article (10.1186/s13062-019-0235-z) contains supplementary material, which is available to authorized users.

Salinity-associated microRNAs and their potential roles in mediating salt tolerance in rice colonized by the endophytic root fungus Piriformospora indica

Piriformospora indica (P. indica), an endophytic root fungus, supports the growth and enhanced tolerance of plants to biotic and abiotic stresses. Several recent studies showed the significant role of small RNA (sRNA) molecules including microRNAs (miRNAs) in plant adaption to environmental stress, but little is known concerning the symbiosis-mediated salt stress tolerance regulated at miRNAs level. The overarching goal of this research is to elucidate the impact of miRNAs in regulating the P. indica-mediated salt tolerance in rice. Applying sRNA-seq analysis led to identify a set of 547 differentially abundant miRNAs in response to P. indica inoculation and salt stress. These included 206 rice-specific and 341 previously known miRNAs from other plant species. In silico analysis of miRNAs predictions of the differentially abundant miRNAs led to identifying of 193 putatively target genes, most of which were encoded either genes or transcription factors involved in nutrient uptake, sodium ion transporters, growth regulators, and auxin- responsive proteins. The rice-specific miRNAs targeted the transcription factors involved in the import of potassium ions into the root cells, the export of sodium ions, and plant growth and development. Interestingly, P. indica affected the differential abundance of miRNAs regulated genes and transcription factors linked to salt stress tolerance. Our data helps to understand the molecular basis of salt stress tolerance mediated by symbionts in plant and the potential impact of miRNAs for genetic improvement of rice varieties for tolerance to salt stress.

Chromosome-based survey sequencing reveals the genome organization of wild wheat progenitor Triticum dicoccoides

Wild emmer wheat (Triticum turgidum ssp. dicoccoides) is the progenitor of wheat. We performed chromosome-based survey sequencing of the 14 chromosomes, examining repetitive sequences, protein-coding genes, miRNA/target pairs and tRNA genes, as well as syntenic relationships with related grasses. We found considerable differences in the content and distribution of repetitive sequences between the A and B subgenomes. The gene contents of individual chromosomes varied widely, not necessarily correlating with chromosome size. We catalogued candidate agronomically important loci, along with new alleles and flanking sequences that can be used to design exome sequencing. Syntenic relationships and virtual gene orders revealed several small-scale evolutionary rearrangements, in addition to providing evidence for the 4AL-5AL-7BS translocation in wild emmer wheat. Chromosome-based sequence assemblies contained five novel miRNA families, among 59 families putatively encoded in the entire genome which provide insight into the domestication of wheat and an overview of the genome content and organization.

Genome-wide screening of long non-coding RNAs involved in rubber biosynthesis in Eucommia ulmoides

Increasing evidence indicates that long non-coding RNAs (lncRNAs) play pivotal roles in regulatory networks controlling plant and animal gene expression. However, lncRNA roles in regulating rubber biosynthesis in Eucommia ulmoides, an emerging source of natural rubber (Eu-rubber), are currently unknown. Here, we report on RNA deep-sequencing of E. ulmoides fruits at two developmental stages. Based on application of a stringent pipeline, 29,103 lncRNAs and 9,048 transcripts of uncertain coding potential (TUCPs) were identified. Two differentially expressed (DE) TUCPs appear to simultaneously regulate 12 protein-coding genes involved in Eu-rubber biosynthesis (GIEBs), as well as 95 DE genes. Functional categorization of these 95 DE genes indicated their involvement in subcellular microstructures and cellular processes, such as cell wall, cell division, and growth. These DE genes may participate in the differentiation and development of laticifers, where Eu-rubber is synthesized. A model is proposed in which "commanders" (DE TUCPs) direct the "builders" (DE genes) to construct a "storehouse" of materials needed for Eu-rubber synthesis, and the "workers" (GIEBs) to synthesize Eu-rubber. These findings provide insights into both cis- and trans-polyisoprene biosynthesis in plants, laying the foundation for additional studies of this crucial process.

Multiple Regression Analysis Reveals MicroRNA Regulatory Networks in Oryza sativa under Drought Stress

Drought is a major abiotic stress that reduces rice development and yield. miRNAs (microRNAs) are known to mediate posttranscriptional regulation under drought stress. Although the importance of individual miRNAs has been established, the crosstalks between miRNAs and mRNAs remain unearthed. Here we performed microarray analysis of miRNAs and matched mRNA expression profiles of drought-treated rice cultivar Nipponbare. Drought-responsive miRNA-mRNA regulations were identified by a combination of a partial least square (PLS) regression approach and sequence-based target prediction. A drought-induced network with 13 miRNAs and 58 target mRNAs was constructed, and four miRNA coregulatory modules were revealed. Functional analysis suggested that drought-response miRNA targets are enriched in hormone signaling, lipid and carbohydrate metabolism, and antioxidant defense. 13 candidate miRNAs and target genes were validated by RT-qPCR, hierarchical clustering, and ROC analysis. Two target genes (DWARF-3 and P0651G05.2) of miRNA coregulatory modules were further verified by RLM-5' RACE. Together, our integrative study of miRNA-mRNA interaction provided attractive candidates that will help elucidate the drought-response mechanisms in Oryza sativa.

Elevated carbon dioxide and drought modulate physiology and storage-root development in sweet potato by regulating microRNAs

Elevated CO₂ along with drought is a serious global threat to crop productivity. Therefore, understanding the molecular mechanisms plants use to protect these stresses is the key for plant growth and development. In this study, we mimicked natural stress conditions under a controlled Soil-Plant-Atmosphere-Research (SPAR) system and provided the evidence for how miRNAs regulate target genes under elevated CO₂ and drought conditions. Significant physiological and biomass data supported the effective utilization of source-sink (leaf to root) under elevated CO₂. Additionally, elevated CO₂ partially rescued the effect of drought on total biomass. We identified both known and novel miRNAs differentially expressed during drought, CO₂, and combined stress, along with putative targets. A total of 32 conserved miRNAs belonged to 23 miRNA families, and 25 novel miRNAs were identified by deep sequencing. Using the existing sweet potato genome database and stringent analyses, a total of 42 and 22 potential target genes were predicted for the conserved and novel miRNAs, respectively. These target genes are involved in drought response, hormone signaling, photosynthesis, carbon fixation, sucrose and starch metabolism, etc. Gene ontology and KEGG ontology functional enrichment revealed that these miRNAs might target transcription factors (MYB, TCP, NAC), hormone signaling regulators (ARF, AP2/ERF), cold and drought factors (corA), carbon metabolism (ATP synthase, fructose-1,6-bisphosphate), and photosynthesis (photosystem I and II complex units). Our study is the first report identifying targets of miRNAs under elevated CO₂ levels and could support the molecular mechanisms under elevated CO₂ in sweet potato and other crops in the future.

Comparative Transcriptome Profiling Analysis of Red- and White-Fleshed Strawberry (Fragaria�ananassa) Provides New Insight into the Regulation of the Anthocyanin Pathway

Anthocyanins are water-soluble pigments in plants. They confer both economic and healthy profits for humans. To gain a deeper insight into the regulation of anthocyanin biosynthesis in octoploid strawberry (Fragaria�ananassa; Fa), a widely consumed economically important fruit, we performed comparative transcriptomic analysis of red- and white-fleshed strawberry cultivars in two ripening stages. In total, 365,455 non-redundant transcripts were assembled from the RNA sequencing (RNAseq) data. Of this collection, 377 were annotated as putative anthocyanin-related transcripts. Differential expression analysis revealed that 57 anthocyanin biosynthesis transcripts were down-regulated, and 89 transcription factors (TFs) were either down- or up-regulated under anthocyanin deficiency. Additionally, amongst the 50,601 putative long non-coding RNAs (lncRNAs) identified here, 68 lncRNAs were differentially expressed and co-expressed with differentially expressed anthocyanin-related mRNAs; 2,070 co-expressing lncRNA-mRNA pairs were generated. Expression profile analysis revealed that it was the limited expression of FaF3'H (flavonoid 3'-hydroxylase) that blocked the cyanidin 3-glucoside accumulation in the two investigated strawberry cultivars. This was further supported by a transient overexpression experiment with FaMYB10. The down-regulated lncRNAs might participate in anthocyanin regulation by acting as targets for microRNAs (miRNAs). The level of competitive intensity in miRNA and lncRNA for the same mRNA targets was probably lower in the white-fleshed strawberries, which can release the repression effect of the mRNAs in red-fleshed strawberry as a result. This study for the first time presents lncRNAs related to anthocyanins in strawberries, provides new insights into the anthocyanin regulatory network and also lays the foundation for identifying new anthocyanin regulators in strawberry.

Transcriptional regulation of osmotic stress tolerance in wheat (Triticum aestivum L.)

The current review provides an updated, new insights into the regulation of transcription mediated underlying mechanisms of wheat plants to osmotic stress perturbations. Osmotic stress tolerance mechanisms being complex are governed by multiple factors at physiological, biochemical and at the molecular level, hence approaches like "OMICS" that can underpin mechanisms behind osmotic tolerance in wheat is of paramount importance. The transcription factors (TFs) are a class of molecular proteins, which are involved in regulation, modulation and orchestrating the responses of plants to a variety of environmental stresses. Recent reports have provided novel insights on the role of TFs in osmotic stress tolerance via direct molecular links. However, our knowledge on the regulatory role TFs during osmotic stress tolerance in wheat remains limited. The present review in its first part sheds light on the importance of studying the role of osmotic stress tolerance in wheat plants and second aims to decipher molecular mechanisms of TFs belonging to several classes, including DREB, NAC, MYB, WRKY and bHLH, which have been reported to engage in osmotic stress mediated gene expression in wheat and third part covers the systems biology approaches to understand the transcriptional regulation of osmotic stress and the role of long non-coding RNAs in response to osmotic stress with special emphasis on wheat. The current concept may lead to an understanding in molecular regulation and signalling interaction of TFs under osmotic stress to clarify challenges and problems for devising potential strategies to improve complex regulatory events involved in plant tolerance to osmotic stress adaptive pathways in wheat.

Genome-wide identification and characterization of mRNAs and lncRNAs involved in cold stress in the wild banana (Musa itinerans)

Cold stress seriously affects banana growth, yield and fruit quality. Long noncoding RNAs (lncRNAs) have been demonstrated as key regulators of biotic and abiotic stress in plants, but the identification and prediction of cold responsive mRNAs and lncRNAs in wild banana remains unexplored. In present study, a cold resistant wild banana line from China was used to profile the cold-responsive mRNAs and lncRNAs by RNA-seq under cold stress conditions, i.e. 13°C (critical growth temperature), 4°C (chilling temperature), 0°C (freezing temperature) and normal growing condition, i.e. 28°C (control group). A total of 12,462 lncRNAs were identified in cold-stressed wild banana. In mRNA, much more alternative splicing events occurred in wild banana under the cold stress conditions compared with that in the normal growing condition. The GO analysis of differential expression genes (DEGs) showed the biochemical processes and membrane related genes responded positively to the cold stress. The KEGG pathway enrichment analysis of the DEGs showed that the pathways of photosynthesis, photosynthesis–antenna proteins, circadian rhythm–plant, glutathione metabolism, starch and sucrose metabolism, cutin/suberine/biosynthesis were altered or affected by the cold stress conditions. Our analyses of the generated transcriptome and lncRNAs provide new insights into regulating expression of genes and lncRNAs that respond to cold stress in the wild banana.

A large-scale multiomics analysis of wheat stem solidness and the wheat stem sawfly feeding response, and syntenic associations in barley, Brachypodium, and rice

The wheat stem sawfly (WSS), Cephus cinctus Norton (Hymenoptera: Cephidae), is an important pest of wheat and other cereals, threatening the quality and quantity of grain production. WSS larvae feed and develop inside the stem where they are protected from the external environment; therefore, pest management strategies primarily rely on host plant resistance. A major locus on the long arm of wheat chromosome 3B underlies most of the variation in stem solidness; however, the impact of stem solidness on WSS feeding has not been completely characterized. Here, we used a multiomics approach to examine the response to WSS in both solid- and semi-solid-stemmed wheat varieties. The combined transcriptomic, proteomic, and metabolomic data revealed that two important molecular pathways, phenylpropanoid and phosphate pentose, are involved in plant defense against WSS. We also detected a general downregulation of several key defense transcripts, including those encoding secondary metabolites such as DIMBOA, tricetin, and lignin, which suggested that the WSS larva might interfere with plant defense. We comparatively analyzed the stem solidness genomic region known to be associated with WSS tolerance in wild emmer, durum, and bread wheats, and described syntenic regions in the close relatives barley, Brachypodium, and rice. Additionally, microRNAs identified from the same genomic region revealed potential regulatory pathways associated with the WSS response. We propose a model outlining the molecular responses of the WSS–wheat interactions. These findings provide insight into the link between stem solidness and WSS feeding at the molecular level.