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

An intronless tau class glutathione transferase detoxifies several herbicides in flufenacet-resistant ryegrass

Resistance to preemergence herbicides, e.g. inhibitors of the biosynthesis of very-long-chain fatty acids (VLCFAs), is evolving in response to increased use of these compounds. Grass weeds such as ryegrasses (Lolium spp.) have accumulated resistance to various herbicide modes of action. Here, an RNA-seq analysis was conducted using 3 ryegrass populations resistant to the VLCFA biosynthesis inhibitor flufenacet to investigate this phenomenon. Besides various transcripts, including putative long noncoding RNAs (lncRNAs), a single putatively functional tau class glutathione transferase (GST) was constitutively differentially expressed. It was further induced by herbicide application. This GST was expressed as a recombinant protein in Escherichia coli along with other GSTs and detoxified flufenacet rapidly in vitro. Detoxification rates of other herbicides tested in vitro were in accordance with cross-resistance patterns previously determined in vivo. A genome-wide GST analysis revealed that the candidate GST was located in a cluster of 3 intronless GSTs. Their intronless nature possibly results from the retroposition of cellular mRNAs followed by tandem duplication and may affect gene expression. The large number of GSTs (≥195) in the genome of rigid ryegrass (Lolium rigidum) compared with other plant organisms is likely a key factor in the ability of this weed to evolve resistance to different herbicide chemistries. However, in the case of flufenacet resistance, a single upregulated GST with high affinity for the substrate flufenacet possibly contributes overproportionally to rapid herbicide detoxification in planta. The regulation of this gene and the role of differentially expressed transcripts, including various putative lncRNAs, require further investigation.

A multispecies study reveals the diversity and potential regulatory role of long noncoding RNAs in cucurbits

Plant long noncoding RNAs (lncRNAs) exhibit features such as tissue-specific expression, spatiotemporal regulation, and stress responsiveness. Although diverse studies support the regulatory role of lncRNAs in model plants, our knowledge about lncRNAs in crops is limited. We employ a custom pipeline on a dataset of over 1000 RNA-seq samples across nine representative species of the family Cucurbitaceae to predict 91 209 nonredundant lncRNAs. The lncRNAs were characterized according to three confidence levels and classified by their genomic context into intergenic, natural antisense, intronic, and sense-overlapping. Compared with protein-coding genes, lncRNAs were, on average, expressed at low levels and displayed significantly higher specificity when considering tissue, developmental stages, and stress responsiveness. The evolutionary analysis indicates higher positional conservation than sequence conservation, probably linked to the conserved modular motifs within syntenic lncRNAs. Moreover, a positive correlation between the expression of intergenic/natural antisense lncRNAs and their closest/parental gene was observed. For those intergenic, the correlation decreases with the distance to the neighboring gene, supporting that their potential cis-regulatory effect is within a short-range. Furthermore, the analysis of developmental studies showed that a conserved NAT-lncRNA family is differentially expressed in a coordinated way with their cognate sense protein-coding genes. These genes code for proteins associated with phloem development, thus providing insights about the potential involvement of some of the identified lncRNAs in a developmental process. We expect that this extensive inventory will constitute a valuable resource for further research lines focused on elucidating the regulatory mechanisms mediated by lncRNAs in cucurbits.

A resource of identified and annotated lincRNAs expressed during somatic embryogenesis development in Norway spruce

Long non-coding RNAs (lncRNAs) have emerged as important regulators of many biological processes, although their regulatory roles remain poorly characterized in woody plants, especially in gymnosperms. A major challenge of working with lncRNAs is to assign functional annotations, since they have a low coding potential and low cross-species conservation. We utilised an existing RNA-Sequencing resource and performed short RNA sequencing of somatic embryogenesis developmental stages in Norway spruce (Picea abies L. Karst). We implemented a pipeline to identify lncRNAs located within the intergenic space (lincRNAs) and generated a co-expression network including protein coding, lincRNA and miRNA genes. To assign putative functional annotation, we employed a guilt-by-association approach using the co-expression network and integrated these results with annotation assigned using semantic similarity and co-expression. Moreover, we evaluated the relationship between lincRNAs and miRNAs, and identified which lincRNAs are conserved in other species. We identified lincRNAs with clear evidence of differential expression during somatic embryogenesis and used network connectivity to identify those with the greatest regulatory potential. This work provides the most comprehensive view of lincRNAs in Norway spruce and is the first study to perform global identification of lincRNAs during somatic embryogenesis in conifers. The data have been integrated into the expression visualisation tools at the PlantGenIE.org web resource to enable easy access to the community. This will facilitate the use of the data to address novel questions about the role of lincRNAs in the regulation of embryogenesis and facilitate future comparative genomics studies.

Underground communication: Long non-coding RNA signaling in the plant rhizosphere

Long non-coding RNAs (lncRNAs) have emerged as integral gene-expression regulators underlying plant growth, development, and adaptation. To adapt to the heterogeneous and dynamic rhizosphere, plants use interconnected regulatory mechanisms to optimally fine-tune gene-expression-governing interactions with soil biota, as well as nutrient acquisition and heavy metal tolerance. Recently, high-throughput sequencing has enabled the identification of plant lncRNAs responsive to rhizosphere biotic and abiotic cues. Here, we examine lncRNA biogenesis, classification, and mode of action, highlighting the functions of lncRNAs in mediating plant adaptation to diverse rhizosphere factors. We then discuss studies that reveal the significance and target genes of lncRNAs during developmental plasticity and stress responses at the rhizobium interface. A comprehensive understanding of specific lncRNAs, their regulatory targets, and the intricacies of their functional interaction networks will provide crucial insights into how these transcriptomic switches fine-tune responses to shifting rhizosphere signals. Looking ahead, we foresee that single-cell dissection of cell-type-specific lncRNA regulatory dynamics will enhance our understanding of the precise developmental modulation mechanisms that enable plant rhizosphere adaptation. Overcoming future challenges through multi-omics and genetic approaches will more fully reveal the integral roles of lncRNAs in governing plant adaptation to the belowground environment.

Exploring the enigma: history, present, and future of long non-coding RNAs in cancer

Long noncoding RNAs (lncRNAs), which are more than 200 nucleotides in length and do not encode proteins, play crucial roles in governing gene expression at both the transcriptional and posttranscriptional levels. These molecules demonstrate specific expression patterns in various tissues and developmental stages, suggesting their involvement in numerous developmental processes and diseases, notably cancer. Despite their widespread acknowledgment and the growing enthusiasm surrounding their potential as diagnostic and prognostic biomarkers, the precise mechanisms through which lncRNAs function remain inadequately understood. A few lncRNAs have been studied in depth, providing valuable insights into their biological activities and suggesting emerging functional themes and mechanistic models. However, the extent to which the mammalian genome is transcribed into functional noncoding transcripts is still a matter of debate. This review synthesizes our current understanding of lncRNA biogenesis, their genomic contexts, and their multifaceted roles in tumorigenesis, highlighting their potential in cancer-targeted therapy. By exploring historical perspectives alongside recent breakthroughs, we aim to illuminate the diverse roles of lncRNA and reflect on the broader implications of their study for understanding genome evolution and function, as well as for advancing clinical applications.

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.

Validation of the Interaction Between Stress-Responsive LncRNAs and RNA-Binding Proteins In Vitro

Plant long non-coding RNAs (lncRNAs) have been implicated in many biological processes, including responses to abiotic stresses, yet their detailed functions and especially their modes of action are still underexplored. lncRNAs often interact with proteins to participate in multiple levels of gene regulation. Therefore, identifying the RNA-binding proteins and validating their interaction with lncRNAs will be instrumental in revealing the functions of lncRNAs. Here, we describe two major methods to determine the interaction between lncRNA and proteins in vitro, RNA pull-down, and RNA EMSA.

Population-level annotation of lncRNAs in Arabidopsis reveals extensive expression variation associated with transposable element-like silencing

Long noncoding RNAs (lncRNAs) are understudied and underannotated in plants. In mammals, lncRNA loci are nearly as ubiquitous as protein-coding genes, and their expression is highly variable between individuals of the same species. Using Arabidopsis thaliana as a model, we aimed to elucidate the true scope of lncRNA transcription across plants from different regions and study its natural variation. We used transcriptome deep sequencing data sets spanning hundreds of natural accessions and several developmental stages to create a population-wide annotation of lncRNAs, revealing thousands of previously unannotated lncRNA loci. While lncRNA transcription is ubiquitous in the genome, most loci appear to be actively silenced and their expression is extremely variable between natural accessions. This high expression variability is largely caused by the high variability of repressive chromatin levels at lncRNA loci. High variability was particularly common for intergenic lncRNAs (lincRNAs), where pieces of transposable elements (TEs) present in 50% of these lincRNA loci are associated with increased silencing and variation, and such lncRNAs tend to be targeted by the TE silencing machinery. We created a population-wide lncRNA annotation in Arabidopsis and improve our understanding of plant lncRNA genome biology, raising fundamental questions about what causes transcription and silencing across the genome.

RNA Polymerases IV and V Are Involved in Olive Fruit Development

Transcription is carried out in most eukaryotes by three multimeric complexes (RNA polymerases I, II and III). However, plants contain two additional RNA polymerases (IV and V), which have evolved from RNA polymerase II. RNA polymerases II, IV and V contain both common and specific subunits that may specialise some of their functions. In this study, we conducted a search for the genes that putatively code for the specific subunits of RNA polymerases IV and V, as well as those corresponding to RNA polymerase II in olive trees. Based on the homology with the genes of Arabidopsis thaliana, we identified 13 genes that putatively code for the specific subunits of polymerases IV and V, and 16 genes that code for the corresponding specific subunits of polymerase II in olives. The transcriptomic analysis by RNA-Seq revealed that the expression of the RNA polymerases IV and V genes was induced during the initial stages of fruit development. Given that RNA polymerases IV and V are involved in the transcription of long non-coding RNAs, we investigated their expression and observed relevant changes in the expression of this type of RNAs. Particularly, the expression of the intergenic and intronic long non-coding RNAs tended to increase in the early steps of fruit development, suggesting their potential role in this process. The positive correlation between the expression of RNA polymerases IV and V subunits and the expression of non-coding RNAs supports the hypothesis that RNA polymerases IV and V may play a role in fruit development through the synthesis of this type of RNAs.

Plant long non-coding RNAs: identification and analysis to unveil their physiological functions

Eukaryotic genomes encode thousands of RNA molecules; however, only a minimal fraction is translated into proteins. Among the non-coding elements, long non-coding RNAs (lncRNAs) play important roles in diverse biological processes. LncRNAs are associated mainly with the regulation of the expression of the genome; nonetheless, their study has just scratched the surface. This is somewhat due to the lack of widespread conservation at the sequence level, in addition to their relatively low and highly tissue-specific expression patterns, which makes their exploration challenging, especially in plant genomes where only a few of these molecules have been described completely. Recently published high-quality genomes of crop plants, along with new computational tools, are considered promising resources for studying these molecules in plants. This review briefly summarizes the characteristics of plant lncRNAs, their presence and conservation, the different protocols to find these elements, and the limitations of these protocols. Likewise, it describes their roles in different plant physiological phenomena. We believe that the study of lncRNAs can help to design strategies to reduce the negative effect of biotic and abiotic stresses on the yield of crop plants and, in the future, help create fruits and vegetables with improved nutritional content, higher amounts of compounds with positive effects on human health, better organoleptic characteristics, and fruits with a longer postharvest shelf life.

Transcriptome analysis reveals a lncRNA-miRNA-mRNA regulatory network in OsRpp30-mediated disease resistance in rice

BACKGROUND

Long non-coding RNAs (lncRNAs) play critical roles in various biological processes in plants. Extensive studies utilizing high-throughput RNA sequencing have revealed that many lncRNAs are involved in plant disease resistance. Oryza sativa RNase P protein 30 (OsRpp30) has been identified as a positive regulator of rice immunity against fungal and bacterial pathogens. Nevertheless, the specific functions of lncRNAs in relation to OsRpp30-mediated disease resistance in rice remain elusive.

RESULTS

We conducted a comprehensive analysis of lncRNAs, miRNAs, and mRNAs expression patterns in wild type (WT), OsRpp30 overexpression (OsRpp30-OE), and OsRpp30 knockout (OsRpp30-KO) rice plants. In total, we identified 91 differentially expressed lncRNAs (DElncRNAs), 1671 differentially expressed mRNAs (DEmRNAs), and 41 differentially expressed miRNAs (DEmiRNAs) across the different rice lines. To gain further insights, we investigated the interaction between DElncRNAs and DEmRNAs, leading to the discovery of 10 trans- and 27 cis-targeting pairs specific to the OsRpp30-OE and OsRpp30-KO samples. In addition, we constructed a competing endogenous RNA (ceRNA) network comprising differentially expressed lncRNAs, miRNAs, and mRNAs to elucidate their intricate interplay in rice disease resistance. The ceRNA network analysis uncovered a set of gene targets regulated by lncRNAs and miRNAs, which were found to be involved in pathogen recognition, hormone pathways, transcription factor activation, and other biological processes related to plant immunity.

CONCLUSIONS

Our study provides a comprehensive expression profiling of lncRNAs, miRNAs, and mRNAs in a collection of defense mutants in rice. To decipher the putative functional significance of lncRNAs, we constructed trans- and cis-targeting networks involving differentially expressed lncRNAs and mRNAs, as well as a ceRNA network incorporating differentially expressed lncRNAs, miRNAs, and mRNAs. Together, the findings from this study provide compelling evidence supporting the pivotal roles of lncRNAs in OsRpp30-mediated disease resistance in rice.

The plant noncoding transcriptome: a versatile environmental sensor

Plant noncoding RNA transcripts have gained increasing attention in recent years due to growing evidence that they can regulate developmental plasticity. In this review article, we comprehensively analyze the relationship between noncoding RNA transcripts in plants and their response to environmental cues. We first provide an overview of the various noncoding transcript types, including long and small RNAs, and how the environment modulates their performance. We then highlight the importance of noncoding RNA secondary structure for their molecular and biological functions. Finally, we discuss recent studies that have unveiled the functional significance of specific long noncoding transcripts and their molecular partners within ribonucleoprotein complexes during development and in response to biotic and abiotic stress. Overall, this review sheds light on the fascinating and complex relationship between dynamic noncoding transcription and plant environmental responses, and highlights the need for further research to uncover the underlying molecular mechanisms and exploit the potential of noncoding transcripts for crop resilience in the context of global warming.

Evaluation of efficacy of non-coding RNA in abiotic stress management of field crops: Current status and future prospective

Abiotic stresses are responsible for the major losses in crop yield all over the world. Stresses generate harmful ROS which can impair cellular processes in plants. Therefore, plants have evolved antioxidant systems in defence against the stress-induced damages. The frequency of occurrence of abiotic stressors has increased several-fold due to the climate change experienced in recent times and projected for the future. This had particularly aggravated the risk of yield losses and threatened global food security. Non-coding RNAs are the part of eukaryotic genome that does not code for any proteins. However, they have been recently found to have a crucial role in the responses of plants to both abiotic and biotic stresses. There are different types of ncRNAs, for example, miRNAs and lncRNAs, which have the potential to regulate the expression of stress-related genes at the levels of transcription, post-transcription, and translation of proteins. The lncRNAs are also able to impart their epigenetic effects on the target genes through the alteration of the status of histone modification and organization of the chromatins. The current review attempts to deliver a comprehensive account of the role of ncRNAs in the regulation of plants' abiotic stress responses through ROS homeostasis. The potential applications ncRNAs in amelioration of abiotic stresses in field crops also have been evaluated.

Genome-wide identification and characterization of Puccinia striiformis-responsive lncRNAs in Triticum aestivum

Wheat stripe rust (yellow rust) caused by Puccinia striiformis f. sp. tritici (Pst) is a serious biotic stress factor limiting wheat production worldwide. Emerging evidence demonstrates that long non-coding RNAs (lncRNAs) participate in various developmental processes in plants via post-transcription regulation. In this study, RNA sequencing (RNA-seq) was performed on a pair of near-isogenic lines-rust resistance line FLW29 and rust susceptible line PBW343-which differed only in the rust susceptibility trait. A total of 6,807 lncRNA transcripts were identified using bioinformatics analyses, among which 10 lncRNAs were found to be differentially expressed between resistance and susceptible lines. In order to find the target genes of the identified lncRNAs, their interactions with wheat microRNA (miRNAs) were predicted. A total of 199 lncRNAs showed interactions with 65 miRNAs, which further target 757 distinct mRNA transcripts. Moreover, detailed functional annotations of the target genes were used to identify the candidate genes, pathways, domains, families, and transcription factors that may be related to stripe rust resistance response in wheat plants. The NAC domain protein, disease resistance proteins RPP13 and RPM1, At1g58400, monodehydroascorbate reductase, NBS-LRR-like protein, rust resistance kinase Lr10-like, LRR receptor, serine/threonine-protein kinase, and cysteine proteinase are among the identified targets that are crucial for wheat stripe rust resistance. Semiquantitative PCR analysis of some of the differentially expressed lncRNAs revealed variations in expression profiles of two lncRNAs between the Pst-resistant and Pst-susceptible genotypes at least under one condition. Additionally, simple sequence repeats (SSRs) were also identified from wheat lncRNA sequences, which may be very useful for conducting targeted gene mapping studies of stripe rust resistance in wheat. These findings improved our understanding of the molecular mechanism responsible for the stripe rust disease that can be further utilized to develop wheat varieties with durable resistance to this disease.

Long non-coding RNA and microRNA landscape of two major domesticated cotton species

Allotetraploid cotton plants Gossypium hirsutum and Gossypium barbadense have been widely cultivated for their natural, renewable textile fibres. Even though ncRNAs in domesticated cotton species have been extensively studied, systematic identification and annotation of lncRNAs and miRNAs expressed in various tissues and developmental stages under various biological contexts are limited. This influences the comprehension of their functions and future research on these cotton species. Here, we report high confidence lncRNAs and miRNA collection from G. hirsutum accession and G. barbadense accession using large-scale RNA-seq and small RNA-seq datasets incorporated into a user-friendly database, CoNCRAtlas. This database provides a wide range and depth of lncRNA and miRNA annotation based on the systematic integration of extensive annotations such as expression patterns derived from transcriptome data analysis in thousands of samples, as well as multi-omics annotations. We assume this comprehensive resource will accelerate evolutionary and functional studies in ncRNAs and inform future breeding programs for cotton improvement. CoNCRAtlas is accessible at http://www.nipgr.ac.in/CoNCRAtlas/.

Genome-wide identification and characterization of long non-coding RNA in barley roots in response to Piriformospora indica colonization

Currently, there is very limited information about long noncoding RNAs (lncRNAs) found in barley. It remains unclear whether barley lncRNAs are responsive to Piriformospora indica (P. indica) colonization.We found that barley roots exhibited fast development and that large roots branched after P. indica colonization. Genome-wide high-throughput RNA-seq and bioinformatic analysis showed that 4356 and 5154 differentially expressed LncRNAs (DELs) were found in response to P. indica at 3 and 7 days after colonization (dai), respectively, and 2456 DELs were found at 7 dai compared to 3 dai. Based on the coexpression correlation of lncRNAmRNA, we found that 98.6% of lncRNAs were positively correlated with 3430 mRNAs at 3 dai and 7 dai. Further GO analysis showed that 30 lncRNAs might be involved in the regulation of gene transcription; 23 lncRNAs might participate in cell cycle regulation. Moreover, the metabolite analysis indicated that chlorophyll a, sucrose, protein, gibberellin, and auxin were in accordance with the results of the transcriptome, and the respective lncRNAs were positively correlated with these target RNAs. Gene silencing suggested that lncRNA TCONS_00262342 is probably a key regulator of GA₃ synthesis pathway, which participates in P. indica and barley interactions. We concluded that acting as a molecular material basis and resource, lncRNAs respond to P. indica colonization by regulating metabolite content in barley and coordinate the complex regulatory process of higher life by constructing highly positive correlations with their target mRNAs.

Roles of long non-coding RNAs in plant immunity

Robust plant immune systems are fine-tuned by both protein-coding genes and non-coding RNAs. Long non-coding RNAs (lncRNAs) refer to RNAs with a length of more than 200 nt and usually do not have protein-coding function and do not belong to any other well-known non-coding RNA types. The non-protein-coding, low expression, and non-conservative characteristics of lncRNAs restrict their recognition. Although studies of lncRNAs in plants are in the early stage, emerging studies have shown that plants employ lncRNAs to regulate plant immunity. Moreover, in response to stresses, numerous lncRNAs are differentially expressed, which manifests the actions of low-expressed lncRNAs and makes plant-microbe/insect interactions a convenient system to study the functions of lncRNAs. Here, we summarize the current advances in plant lncRNAs, discuss their regulatory effects in different stages of plant immunity, and highlight their roles in diverse plant-microbe/insect interactions. These insights will not only strengthen our understanding of the roles and actions of lncRNAs in plant-microbe/insect interactions but also provide novel insight into plant immune responses and a basis for further research in this field.

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.

De novo full-length transcriptome analysis of two ecotypes of Phragmites australis (swamp reed and dune reed) provides new insights into the transcriptomic complexity of dune reed and its long-term adaptation to desert environments

BACKGROUND

The extremely harsh environment of the desert is changing dramatically every moment, and the rapid adaptive stress response in the short term requires enormous energy expenditure to mobilize widespread regulatory networks, which is all the more detrimental to the survival of the desert plants themselves. The dune reed, which has adapted to desert environments with complex and variable ecological factors, is an ideal type of plant for studying the molecular mechanisms by which Gramineae plants respond to combinatorial stress of the desert in their natural state. But so far, the data on the genetic resources of reeds is still scarce, therefore most of their research has focused on ecological and physiological studies.

RESULTS

In this study, we obtained the first De novo non-redundant Full-Length Non-Chimeric (FLNC) transcriptome databases for swamp reeds (SR), dune reeds (DR) and the All of Phragmites australis (merged of iso-seq data from SR and DR), using PacBio Iso-Seq technology and combining tools such as Iso-Seq3 and Cogent. We then identified and described long non-coding RNAs (LncRNA), transcription factor (TF) and alternative splicing (AS) events in reeds based on a transcriptome database. Meanwhile, we have identified and developed for the first time a large number of candidates expressed sequence tag-SSR (EST-SSRs) markers in reeds based on UniTransModels. In addition, through differential gene expression analysis of wild-type and homogenous cultures, we found a large number of transcription factors that may be associated with desert stress tolerance in the dune reed, and revealed that members of the Lhc family have an important role in the long-term adaptation of dune reeds to desert environments.

CONCLUSIONS

Our results provide a positive and usable genetic resource for Phragmites australis with a widespread adaptability and resistance, and provide a genetic database for subsequent reeds genome annotation and functional genomic studies.

Regulatory roles of noncoding RNAs in callus induction and plant cell dedifferentiation

Plant regulatory noncoding RNAs (ncRNAs) have emerged as key modulators of gene expression during callus induction. Their further study may promote the design of innovative plant tissue culture protocols. The use of plants by humans has recently taken on a new and expanding insight due to the advent of genetic engineering technologies. In this context, callus cultures have shown remarkable potential for synthesizing valuable biomolecules, crop improvement, plant micropropagation, and biodiversity preservation. A crucial stage in callus production is the conversion of somatic cells into totipotent cells; compelling evidence indicates that stress factors, transcriptional regulators, and plant hormones can trigger this biological event. Besides, posttranscriptional regulators of gene expression might be essential participants in callus induction. However, research related to the analysis of noncoding RNAs (ncRNAs) that modulate callogenesis and plant cell dedifferentiation in vitro is still at an early stage. During the last decade, some relevant studies have enlightened the fact that different classes of ncRNAs, such as microRNAs (miRNAs), small interfering RNAs (siRNAs), and long noncoding RNAs (lncRNAs) are implicated in plant cell dedifferentiation through regulating the expression levels of diverse gene targets. Hence, understanding the molecular relevance of these ncRNAs in the aforesaid biological processes might represent a promising source of new biotechnological approaches for callus culture and plant improvement. In this current work, we review the experimental evidence regarding the prospective roles of ncRNAs in callus induction and plant cell dedifferentiation to promote this field of study.

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.

Reverse engineering environmental metatranscriptomes clarifies best practices for eukaryotic assembly

BACKGROUND

Diverse communities of microbial eukaryotes in the global ocean provide a variety of essential ecosystem services, from primary production and carbon flow through trophic transfer to cooperation via symbioses. Increasingly, these communities are being understood through the lens of omics tools, which enable high-throughput processing of diverse communities. Metatranscriptomics offers an understanding of near real-time gene expression in microbial eukaryotic communities, providing a window into community metabolic activity.

RESULTS

Here we present a workflow for eukaryotic metatranscriptome assembly, and validate the ability of the pipeline to recapitulate real and manufactured eukaryotic community-level expression data. We also include an open-source tool for simulating environmental metatranscriptomes for testing and validation purposes. We reanalyze previously published metatranscriptomic datasets using our metatranscriptome analysis approach.

CONCLUSION

We determined that a multi-assembler approach improves eukaryotic metatranscriptome assembly based on recapitulated taxonomic and functional annotations from an in-silico mock community. The systematic validation of metatranscriptome assembly and annotation methods provided here is a necessary step to assess the fidelity of our community composition measurements and functional content assignments from eukaryotic metatranscriptomes.

A long-read and short-read transcriptomics approach provides the first high-quality reference transcriptome and genome annotation for Pseudotsuga menziesii (Douglas-fir)

Douglas-fir (Pseudotsuga menziesii) is native to western North America. It grows in a wide range of environmental conditions and is an important timber tree. Although there are several studies on the gene expression responses of Douglas-fir to abiotic cues, the absence of high-quality transcriptome and genome data is a barrier to further investigation. Like for most conifers, the available transcriptome and genome reference dataset for Douglas-fir remains fragmented and requires refinement. We aimed to generate a highly accurate, and complete reference transcriptome and genome annotation. We deep-sequenced the transcriptome of Douglas-fir needles from seedlings that were grown under nonstress control conditions or a combination of heat and drought stress conditions using long-read (LR) and short-read (SR) sequencing platforms. We used 2 computational approaches, namely de novo and genome-guided LR transcriptome assembly. Using the LR de novo assembly, we identified 1.3X more high-quality transcripts, 1.85X more "complete" genes, and 2.7X more functionally annotated genes compared to the genome-guided assembly approach. We predicted 666 long noncoding RNAs and 12,778 unique protein-coding transcripts including 2,016 putative transcription factors. We leveraged the LR de novo assembled transcriptome with paired-end SR and a published single-end SR transcriptome to generate an improved genome annotation. This was conducted with BRAKER2 and refined based on functional annotation, repetitive content, and transcriptome alignment. This high-quality genome annotation has 51,419 unique gene models derived from 322,631 initial predictions. Overall, our informatics approach provides a new reference Douglas-fir transcriptome assembly and genome annotation with considerably improved completeness and functional annotation.

Wheat Long Noncoding RNAs from Organelle and Nuclear Genomes Carry Conserved microRNA Precursors Which May Together Comprise Intricate Networks in Insect Responses

Long noncoding RNAs (lncRNAs) are a diverse class of noncoding RNAs that are typically longer than 200 nucleotides but lack coding potentials. Advances in deep sequencing technologies enabled a better exploration of this type of noncoding transcripts. The poor sequence conservation, however, complicates the identification and annotation of lncRNAs at a large scale. Wheat is among the leading food staples worldwide whose production is threatened by both biotic and abiotic stressors. Here, we identified putative lncRNAs from durum wheat varieties that differ in stem solidness, a major source of defense against wheat stem sawfly, a devastating insect pest. We also analyzed and annotated lncRNAs from two bread wheat varieties, resistant and susceptible to another destructive pest, orange wheat blossom midge, with and without infestation. Several putative lncRNAs contained potential precursor sequences and/or target regions for microRNAs, another type of regulatory noncoding RNAs, which may indicate functional networks. Interestingly, in contrast to lncRNAs themselves, microRNAs with potential precursors within the lncRNA sequences appeared to be highly conserved at the sequence and family levels. We also observed a few putative lncRNAs that have perfect to near-perfect matches to organellar genomes, supporting the recent observations that organellar genomes may contribute to the noncoding transcript pool of the cell.

A long non-coding RNA PelncRNA1 is involved in Phyllostachys edulis response to UV-B stress

Phyllostachys edulis (moso bamboo) is China's most widespread bamboo species, with significant economic and ecological values. Long non-coding RNA (lncRNA) is a type of regulatory RNA that is longer than 200 nucleotides and incapable of encoding proteins, and is frequently involved in regulating biotic and abiotic stress and plant development. However, the biological functions of lncRNA in moso bamboo are unknown. In this study, a lncRNA (named PelncRNA1) differentially expressed following UV-B treatment was discovered in the whole transcriptome sequencing database of moso bamboo. The target genes were filtered and defined by correlation analysis of PelncRNA1 and gene expression pattern. The expression levels of PelncRNA1 and its target genes were verified using qRT-PCR. The results demonstrated that the expression levels of PelncRNA1 and its target genes increased during UV-B treatment. In Arabidopsis transgenic seedlings and moso bamboo protoplasts, PelncRNA1 was discovered to influence the expression of its target genes when overexpressed. In addition, transgenic Arabidopsis showed higher tolerance to UV-B stress. These results suggest that PelncRNA1 and its target genes are involved in the response of moso bamboo to UV-B stress. The novel findings would contribute to our understanding of how lncRNAs regulate the response to abiotic stresses in moso bamboo.

Role of long non-coding RNAs in rice reproductive development

Rice is a staple crop, feeding over half of the global population. The future demand of population growth and climate change requires substantial rice improvement. Recent advances in rice genomics have highlighted the vital role of the non-coding part of the genome. The protein-coding regions account for only a tiny portion of the eukaryotic genome, and most of the genomic regions transcribe copious amounts of non-coding RNAs. Of these, the long non-coding RNAs, including linear non-coding RNAs (lncRNAs) and circular non-coding RNAs (circRNAs), have been shown to play critical roles in various developmental processes by regulating the expression of genes and functions of proteins at transcriptional, post-transcriptional and post-translational levels. With the advances in next-generation sequencing technologies, a substantial number of long non-coding RNAs have been found to be expressed in plant reproductive organs in a cell- and tissue-specific manner suggesting their reproductive development-related functions. Accumulating evidence points towards the critical role of these non-coding RNAs in flowering, anther, and pollen development, ovule and seed development and photoperiod and temperature regulation of male fertility. In this mini review, we provide a brief overview of the role of the linear and circular long non-coding RNAs in rice reproductive development and control of fertility and crop yield.

LncDC: a machine learning-based tool for long non-coding RNA detection from RNA-Seq data

Long non-coding RNAs (lncRNAs) play an essential role in diverse biological processes and disease development. Accurate classification of lncRNAs and mRNAs is important for the identification of tissue- or disease-specific lncRNAs. Here, we present our tool LncDC (Long non-coding RNA detection) that is able to accurately predict lncRNAs with an XGBoost model using features extracted from RNA sequences, secondary structures, and translated proteins. Benchmarking experiments showed that LncDC consistently outperformed six state-of-the-art tools in distinguishing lncRNAs from mRNAs. Notably, the use of sequence and secondary structure (SASS) k-mer score features and flexible ORF features improved the classification capability of LncDC. We anticipate that LncDC will definitely promote the discovery of more and novel disease-specific lncRNAs. LncDC is implemented in Python and freely available at https://github.com/lim74/LncDC .

Genome-wide identification and characterization of lncRNAs in sunflower endosperm

BACKGROUND

Long non-coding RNAs (lncRNAs), as important regulators, play important roles in plant growth and development. The expression and epigenetic regulation of lncRNAs remain uncharacterized generally in plant seeds, especially in the transient endosperm of the dicotyledons.

RESULTS

In this study, we identified 11,840 candidate lncRNAs in 12 day-after-pollination sunflower endosperm by analyzing RNA-seq data. These lncRNAs were evenly distributed in all chromosomes and had specific features that were distinct from mRNAs including tissue-specificity expression, shorter and fewer exons. By GO analysis of protein coding genes showing strong correlation with the lncRNAs, we revealed that these lncRNAs potential function in many biological processes of seed development. Additionally, genome-wide DNA methylation analyses revealed that the level of DNA methylation at the transcription start sites was negatively correlated with gene expression levels in lncRNAs. Finally, 36 imprinted lncRNAs were identified including 32 maternally expressed lncRNAs and four paternally expressed lncRNAs. In CG and CHG context, DNA methylation levels of imprinted lncRNAs in the upstream and gene body regions were slightly lower in the endosperm than that in embryo tissues, which indicated that the maternal demethylation potentially induce the paternally bias expression of imprinted lncRNAs in sunflower endosperm.

CONCLUSION

Our findings not only identified and characterized lncRNAs on a genome-wide scale in the development of sunflower endosperm, but also provide novel insights into the parental effects and epigenetic regulation of lncRNAs in dicotyledonous seeds.

Non-coding RNA expression analysis revealed the molecular mechanism of flag leaf heterosis in inter-subspecific hybrid rice

Heterosis has been used widespread in agriculture, but its molecular mechanism is inadequately understood. Plants have a large number of non-coding RNAs (ncRNAs), among them, functional ncRNAs that have been studied widely containing long non-coding RNA (lncRNA) and circular RNA (circRNA) that play a role in varied biological processes, as well as microRNA (miRNA), which can not only regulate the post-transcriptional expression of target genes, but also target lncRNA and circRNA then participate the competing endogenous RNA (ceRNA) regulatory network. However, the influence of these three ncRNAs and their regulatory relationships on heterosis is unknown in rice. In this study, the expression profile of ncRNAs and the ncRNA regulatory network related to heterosis were comprehensively analyzed in inter-subspecific hybrid rice. A total of 867 miRNAs, 3,278 lncRNAs and 2,521 circRNAs were identified in the hybrid and its parents. Analysis of the global profiles of these three types of ncRNAs indicated that significant differences existed in the distribution and sequence characteristics of the corresponding genes. The numbers of miRNA and lncRNA in hybrid were higher than those in its parents. A total of 784 ncRNAs (169 miRNAs, 573 lncRNAs and 42 circRNAs) showed differentially expressed in the hybrid, and their target/host genes were vital in stress tolerance, growth and development in rice. These discoveries suggested that the expression plasticity of ncRNA has an important role of inter-subspecific hybrid rice heterosis. It is worth mentioning that miRNAs exhibited substantially more variations between hybrid and parents compared with observed variation for lncRNA and circRNA. Non-additive expression ncRNAs and allele-specific expression genes-related ncRNAs in hybrid were provided in this study, and multiple sets of ncRNA regulatory networks closely related to heterosis were obtained. Meanwhile, heterosis-related regulatory networks of ceRNA (lncRNA and circRNA) and miRNA were also demonstrated.

A Long Noncoding RNA Derived from lncRNA-mRNA Networks Modulates Seed Vigor

The discovery of long noncoding RNAs (lncRNAs) has filled a great gap in our understanding of posttranscriptional gene regulation in a variety of biological processes related to plant stress responses. However, systematic analyses of the lncRNAs expressed in rice seeds that germinate under cold stress have been elusive. In this study, we performed strand-specific whole transcriptome sequencing in germinated rice seeds under cold stress and normal temperature. A total of 6258 putative lncRNAs were identified and expressed in a stage-specific manner compared to mRNA. By investigating the targets of differentially expressed (DE) lncRNAs of LT-I (phase I of low temperature)/NT-I (phase I of normal temperature), it was shown that the auxin-activated signaling pathway was significantly enriched, and twenty-three protein-coding genes with most of the members of the SAUR family located in chromosome 9 were identified as the candidate target genes that may interact with five lncRNAs. A seed vigor-related lncRNA, SVR, which interplays with the members of the SAUR gene family in cis was eventually identified. The CRISPR/Cas 9 engineered mutations in SVR cause delay of germination. The findings provided new insights into the connection between lncRNAs and the auxin-activated signaling pathway in the regulation of rice seed vigor.

Global Identification of White Lupin lncRNAs Reveals Their Role in Cluster Roots under Phosphorus Deficiency

Phosphorus (P) deficiency heterogeneously affected plant nutritional status and physiological performance, ultimately leading to a severe yield reduction. A few putative long non-coding RNAs (lncRNAs) responding to P-starvation in the model crops Arabidopsis thaliana and Oryza sativa have been characterized. White lupin (Lupinus albus) is of prime importance, and is a legume with increasing agronomic value as a protein crop as it exhibits extreme tolerance to nutrient deficiency, particularly P deficiency. Despite its adapted nature to P deficiency, nothing is known about low P-induced lncRNAs in white lupin roots. To address this issue, we identified 39,840 mRNA and 2028 lncRNAs in the eight developmental stages of white lupin root (S0-S7 and lateral root, LR) grown under P deficiency. From these 2028 lncRNAs, 1564 were intergenic and 464 natural antisense intergenic transcript (NAT) lncRNAs. We further predicted six potential targets of miRNAs with twelve lncRNAs, which may regulate P-deficiency-related processes. Moreover, the weighted gene co-expression network analysis (WGCNA) revealed seven modules that were correlated with the expression pattern of lncRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed 606 GO terms and 27 different pathways including signal transduction, energy synthesis, detoxification, and Pi transport. In addition, we screened 13 putative lncRNAs that showed a distinct expression pattern in each root, indicating their role in the P deficiency regulatory network. Therefore, white lupin may be a reference legume to characterize P-deficiency-responsive novel lncRNAs, which would highlight the role of lncRNAs in the regulation of plant responses to P deficiency.

B, Iquebal MA, Bhattacharya RC, Rai A, Kumar D. Genome wide identification of lncRNAs and circRNAs having regulatory role in fruit shelf life in health crop cucumber (Cucumis sativus L.)

Cucumber is an extremely perishable vegetable; however, under room conditions, the fruits become unfit for consumption 2-3 days after harvesting. One natural variant, DC-48 with an extended shelf-life was identified, fruits of which can be stored up to 10-15 days under room temperature. The genes involved in this economically important trait are regulated by non-coding RNAs. The study aims to identify the long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) by taking two contrasting genotypes, DC-48 and DC-83, at two different fruit developmental stages. The upper epidermis of the fruits was collected at 5 days and 10 days after pollination (DAP) for high throughput RNA sequencing. The differential expression analysis was performed to identify differentially expressed (DE) lncRNAs and circRNAs along with the network analysis of lncRNA, miRNA, circRNA, and mRNA interactions. A total of 97 DElncRNAs were identified where 18 were common under both the developmental stages (8 down regulated and 10 upregulated). Based on the back-spliced reads, 238 circRNAs were found to be distributed uniformly throughout the cucumber genomes with the highest numbers (71) in chromosome 4. The majority of the circRNAs (49%) were exonic in origin followed by inter-genic (47%) and intronic (4%) origin. The genes related to fruit firmness, namely, polygalacturonase, expansin, pectate lyase, and xyloglucan glycosyltransferase were present in the target sites and co-localized networks indicating the role of the lncRNA and circRNAs in their regulation. Genes related to fruit ripening, namely, trehalose-6-phosphate synthase, squamosa promoter binding protein, WRKY domain transcription factors, MADS box proteins, abscisic stress ripening inhibitors, and different classes of heat shock proteins (HSPs) were also found to be regulated by the identified lncRNA and circRNAs. Besides, ethylene biosynthesis and chlorophyll metabolisms were also found to be regulated by DElncRNAs and circRNAs. A total of 17 transcripts were also successfully validated through RT PCR data. These results would help the breeders to identify the complex molecular network and regulatory role of the lncRNAs and circRNAs in determining the shelf-life of cucumbers.

DeepPlnc: Bi-modal deep learning for highly accurate plant lncRNA discovery

We present here a bi-modal CNN based deep-learning system, DeepPlnc, to identify plant lncRNAs with high accuracy while using sequence and structural properties. Unlike most of the existing software, it works accurately even in conditions with ambiguity of boundaries and incomplete sequences. It scored consistently high for performance metrics while breaching accuracy of >98% when tested across a large number of validated instances. During multiple benchmarkings it consistently outperformed all the compared tools and maintained a highly significant lead in the range of 2.5%- 4.6% from the second best performing tool (p-value << 0.01). DeepPlnc was used to annotate a de novo assembled transcriptome of a himalayan species where again it suggested its much better suitability for genome and transcriptome annotation purposes than the existing tools. DeepPlnc has been made freely available as a web-server and stand-alone program at https://scbb.ihbt.res.in/DeepPlnc/.

Full-Length Transcriptional Analysis of the Same Soybean Genotype With Compatible and Incompatible Reactions to Heterodera glycines Reveals Nematode Infection Activating Plant Defense Response

Full-length transcriptome sequencing with long reads is a powerful tool to analyze transcriptional and post-transcriptional events; however, it has not been applied on soybean (Glycine max). Here, a comparative full-length transcriptome analysis was performed on soybean genotype 09-138 infected with soybean cyst nematode (SCN, Heterodera glycines) race 4 (SCN4, incompatible reaction) and race 5 (SCN5, compatible reaction) using Oxford Nanopore Technology. Each of 9 full-length samples collected 8 days post inoculation with/without nematodes generated an average of 6.1 GB of clean data and a total of 65,038 transcript sequences. After redundant transcripts were removed, 1,117 novel genes and 41,096 novel transcripts were identified. By analyzing the sequence structure of the novel transcripts, a total of 28,759 complete open reading frame (ORF) sequences, 5,337 transcription factors, 288 long non-coding RNAs, and 40,090 novel transcripts with function annotation were predicted. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of differentially expressed genes (DEGs) revealed that growth hormone, auxin-activated signaling pathway and multidimensional cell growth, and phenylpropanoid biosynthesis pathway were enriched by infection with both nematode races. More DEGs associated with stress response elements, plant-hormone signaling transduction pathway, and plant-pathogen interaction pathway with more upregulation were found in the incompatible reaction with SCN4 infection, and more DEGs with more upregulation involved in cell wall modification and carbohydrate bioprocess were detected in the compatible reaction with SCN5 infection when compared with each other. Among them, overlapping DEGs with a quantitative difference was triggered. The combination of protein-protein interaction with DEGs for the first time indicated that nematode infection activated the interactions between transcription factor WRKY and VQ (valine-glutamine motif) to contribute to soybean defense. The knowledge of the SCN-soybean interaction mechanism as a model will present more understanding of other plant-nematode interactions.

Sequencing and de novo transcriptome assembly for discovering regulators of gene expression in Jack (Artocarpus heterophyllus)

Jack (Artocarpus heterophyllus) is a multipurpose fruit-tree species with minimal genomic resources. The study reports developing comprehensive transcriptome data containing 80,411 unigenes with an N50 value of 1265 bp. We predicted 64,215 CDSs from the unigenes and annotated and functionally categorized them into the biological process (23,230), molecular function (27,149), and cellular components (17,284). From 80,411 unigenes, we discovered 16,853 perfect SSRs with 192 distinct repeat motif types reiterating 4 to 22 times. Besides, we identified 2741 TFs from 69 TF families, 53 miRNAs from 19 conserved miRNA families, 25,953 potential lncRNAs, and placed three functional eTMs in different lncRNA-miRNA pairs. The regulatory networks involving genes, TFs, and miRNAs identified several regulatory and regulated nodes providing insight into miRNAs' gene associations and transcription factor-mediated regulation. The comparison of expression patterns of some selected miRNAs vis-à-vis their corresponding target genes showed an inverse relationship indicating the possible miRNA-mediated regulation of the genes.

Role of phasiRNAs in plant-pathogen interactions: molecular perspectives and bioinformatics tools

The genome of an organism is regulated in concert with the organized action of various genetic regulators at different hierarchical levels. Small non-coding RNAs are one of these regulators, among which microRNAs (miRNAs), a distinguished sRNA group with decisive functions in the development, growth and stress-responsive activities of both plants as well as animals, are keenly explored over a good number of years. Recent studies in plants revealed that apart from the silencing activity exhibited by miRNAs on their targets, miRNAs of specific size and structural features can direct the phasing pattern of their target loci to form phased secondary small interfering RNAs (phasiRNAs). These trigger-miRNAs were identified to target both coding and long non-coding RNAs that act as potent phasiRNA precursors or PHAS loci. The phasiRNAs produced thereby exhibit a role in enhancing further downstream regulation either on their own precursors or on those transcripts that are distinct from their genetic source of origin. Hence, these tiny regulators can stimulate an elaborative cascade of interacting RNA networks via cis and trans-regulatory mechanisms. Our review focuses on the comprehensive understanding of phasiRNAs and their trigger miRNAs, by giving much emphasis on their role in the regulation of plant defense responses, together with a summary of the computational tools available for the prediction of the same.

Single-molecule real-time sequencing of the full-length transcriptome of purple garlic (Allium sativum L. cv. Leduzipi) and identification of serine O-acetyltransferase family proteins involved in cysteine biosynthesis

BACKGROUND

Garlic (Allium sativum L.), whose bioactive components are mainly organosulfur compounds (OSCs), is a herbaceous perennial widely consumed as a green vegetable and a condiment. Yet, the metabolic enzymes involved in the biosynthesis of OSCs are not identified in garlic.

RESULTS

Here, a full-length transcriptome of purple garlic was generated via PacBio and Illumina sequencing, to characterize the garlic transcriptome and identify key proteins mediating the biosynthesis of OSCs. Overall, 22.56 Gb of clean data were generated, resulting in 454 698 circular consensus sequence (CCS) reads, of which 83.4% (379 206) were identified as being full-length non-chimeric reads - their further transcript clustering facilitated identification of 36 571 high-quality consensus reads. Once corrected, their genome-wide mapping revealed that 6140 reads were novel isoforms of known genes, and 2186 reads were novel isoforms from novel genes. We detected 1677 alternative splicing events, finding 2902 genes possessing either two or more poly(A) sites. Given the importance of serine O-acetyltransferase (SERAT) in cysteine biosynthesis, we investigated the five SERAT homologs in garlic. Phylogenetic analysis revealed a three-tier classification of SERAT proteins, each featuring a serine acetyltransferase domain (N-terminal) and one or two hexapeptide transferase motifs. Template-based modeling showed that garlic SERATs shared a common homo-trimeric structure with homologs from bacteria and other plants. The residues responsible for substrate recognition and catalysis were highly conserved, implying a similar reaction mechanism. In profiling the five SERAT genes' transcript levels, their expression pattern varied significantly among different tissues.

CONCLUSION

This study's findings deepen our knowledge of SERAT proteins, and provide timely genetic resources that could advance future exploration into garlic's genetic improvement and breeding. © 2021 Society of Chemical Industry.

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.

Genome-wide identification and characterization of Fusarium circinatum-responsive lncRNAs in Pinus radiata

Background

One of the most promising strategies of Pine Pitch Canker (PPC) management is the use of reproductive plant material resistant to the disease. Understanding the complexity of plant transcriptome that underlies the defence to the causal agent Fusarium circinatum, would greatly facilitate the development of an accurate breeding program. Long non-coding RNAs (lncRNAs) are emerging as important transcriptional regulators under biotic stresses in plants. However, to date, characterization of lncRNAs in conifer trees has not been reported. In this study, transcriptomic identification of lncRNAs was carried out using strand-specific paired-end RNA sequencing, from Pinus radiata samples inoculated with F. circinatum at an early stage of infection.

Results

Overall, 13,312 lncRNAs were predicted through a bioinformatics approach, including long intergenic non-coding RNAs (92.3%), antisense lncRNAs (3.3%) and intronic lncRNAs (2.9%). Compared with protein-coding RNAs, pine lncRNAs are shorter, have lower expression, lower GC content and harbour fewer and shorter exons. A total of 164 differentially expressed (DE) lncRNAs were identified in response to F. circinatum infection in the inoculated versus mock-inoculated P. radiata seedlings. The predicted cis-regulated target genes of these pathogen-responsive lncRNAs were related to defence mechanisms such as kinase activity, phytohormone regulation, and cell wall reinforcement. Co-expression network analysis of DE lncRNAs, DE protein-coding RNAs and lncRNA target genes also indicated a potential network regulating pectinesterase activity and cell wall remodelling.

Conclusions

This study presents the first comprehensive genome-wide analysis of P. radiata lncRNAs and provides the basis for future functional characterizations of lncRNAs in relation to pine defence responses against F. circinatum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08408-9.

Drought tolerance improvement in Solanum lycopersicum: an insight into "OMICS" approaches and genome editing

Solanum lycopersicum (tomato) is an internationally acclaimed vegetable crop that is grown worldwide. However, drought stress is one of the most critical challenges for tomato production, and it is a crucial task for agricultural biotechnology to produce drought-resistant cultivars. Although breeders have done a lot of work on the tomato to boost quality and quantity of production and enhance resistance to biotic and abiotic stresses, conventional tomato breeding approaches have been limited to improving drought tolerance because of the intricacy of drought traits. Many efforts have been made to better understand the mechanisms involved in adaptation and tolerance to drought stress in tomatoes throughout the years. "Omics" techniques, such as genomics, transcriptomics, proteomics, and metabolomics in combination with modern sequencing technologies, have tremendously aided the discovery of drought-responsive genes. In addition, the availability of biotechnological tools, such as plant transformation and the recently developed genome editing system for tomatoes, has opened up wider opportunities for validating the function of drought-responsive genes and the generation of drought-tolerant varieties. This review highlighted the recent progresses for tomatoes improvement against drought stress through "omics" and "multi-omics" technologies including genetic engineering. We have also discussed the roles of non-coding RNAs and genome editing techniques for drought stress tolerance improvement in tomatoes.

Differential microRNA expression, microRNA arm switching, and microRNA:long noncoding RNA interaction in response to salinity stress in soybean

BACKGROUND

Soybean is a major legume crop with high nutritional and environmental values suitable for sustainable agriculture. Noncoding RNAs (ncRNAs), including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), are important regulators of gene functions in eukaryotes. However, the interactions between these two types of ncRNAs in the context of plant physiology, especially in response to salinity stress, are poorly understood.

RESULTS

Here, we challenged a cultivated soybean accession (C08) and a wild one (W05) with salt treatment and obtained their small RNA transcriptomes at six time points from both root and leaf tissues. In addition to thoroughly analyzing the differentially expressed miRNAs, we also documented the first case of miRNA arm-switching (miR166m), the swapping of dominant miRNA arm expression, in soybean in different tissues. Two arms of miR166m target different genes related to salinity stress (chloroplastic beta-amylase 1 targeted by miR166m-5p and calcium-dependent protein kinase 1 targeted by miR166m-3p), suggesting arm-switching of miR166m play roles in soybean in response to salinity stress. Furthermore, two pairs of miRNA:lncRNA interacting partners (miR166i-5p and lncRNA Gmax_MSTRG.35921.1; and miR394a-3p and lncRNA Gmax_MSTRG.18616.1) were also discovered in reaction to salinity stress.

CONCLUSIONS

This study demonstrates how ncRNA involves in salinity stress responses in soybean by miRNA arm switching and miRNA:lncRNA interactions. The behaviors of ncRNAs revealed in this study will shed new light on molecular regulatory mechanisms of stress responses in plants, and hence provide potential new strategies for crop improvement.

Long Non-Coding RNAs as Emerging Regulators of Pathogen Response in Plants

Long non-coding RNAs (lncRNAs) are transcripts without protein-coding potential that contain more than 200 nucleotides that play important roles in plant survival in response to different stresses. They interact with molecules such as DNA, RNA, and protein, and play roles in the regulation of chromatin remodeling, RNA metabolism, and protein modification activities. These lncRNAs regulate the expression of their downstream targets through epigenetic changes, at the level of transcription and post-transcription. Emerging information from computational biology and functional characterization of some of them has revealed their diverse mechanisms of action and possible roles in biological processes such as flowering time, reproductive organ development, as well as biotic and abiotic stress responses. In this review, we have mainly focused on the role of lncRNAs in biotic stress response due to the limited availability of knowledge in this domain. We have discussed the available molecular mechanisms of certain known lncRNAs against specific pathogens. Further, considering that fungal, viral, and bacterial diseases are major factors in the global food crisis, we have highlighted the importance of lncRNAs against pathogen responses and the progress in plant research to develop a better understanding of their functions and molecular mechanisms.

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.

The evolutionary landscape and expression pattern of plant lincRNAs

Long intergenic non-coding RNAs (lincRNAs) are important regulators of cellular processes, including development and stress response. Many lincRNAs have been bioinformatically identified in plants, but their evolutionary dynamics and expression characteristics are still elusive. Here, we systematically identified thousands of lincRNAs in 26 plant species, including 6 non-flowering plants, investigated the conservation of the identified lincRNAs in different levels of plant lineages based on sequence and/or synteny homology and explored characteristics of the conserved lincRNAs during plant evolution and their co-expression relationship with protein-coding genes (PCGs). In addition to confirmation of the features well documented in literature for lincRNAs, such as species-specific, fewer exons, tissue-specific expression patterns and less abundantly expressed, we revealed that histone modification signals and/or binding sites of transcription factors were enriched in the conserved lincRNAs, implying their biological functionalities, as demonstrated by identifying conserved lincRNAs related to flower development in both the Brassicaceae and grass families and ancient lincRNAs potentially functioning in meristem development of non-flowering plants. Compared to PCGs, lincRNAs are more likely to be associated with transposable elements (TEs), but with different characteristics in different evolutionary lineages, for instance, the types of TEs and the variable level of association in lincRNAs with different conservativeness. Together, these results provide a comprehensive view on the evolutionary landscape of plant lincRNAs and shed new insights on the conservation and functionality of plant lincRNAs.

Epigenetics: a catalyst of plant immunity against pathogens

The plant immune system protects against pests and diseases. The recognition of stress-related molecular patterns triggers localised immune responses, which are often followed by longer-lasting systemic priming and/or up-regulation of defences. In some cases, this induced resistance (IR) can be transmitted to following generations. Such transgenerational IR is gradually reversed in the absence of stress at a rate that is proportional to the severity of disease experienced in previous generations. This review outlines the mechanisms by which epigenetic responses to pathogen infection shape the plant immune system across expanding time scales. We review the cis- and trans-acting mechanisms by which stress-inducible epigenetic changes at transposable elements (TEs) regulate genome-wide defence gene expression and draw particular attention to one regulatory model that is supported by recent evidence about the function of AGO1 and H2A.Z in transcriptional control of defence genes. Additionally, we explore how stress-induced mobilisation of epigenetically controlled TEs acts as a catalyst of Darwinian evolution by generating (epi)genetic diversity at environmentally responsive genes. This raises questions about the long-term evolutionary consequences of stress-induced diversification of the plant immune system in relation to the long-held dichotomy between Darwinian and Lamarckian evolution.

Genome-Wide Identification and Characterization of Long Noncoding RNAs in Populus × canescens Roots Treated With Different Nitrogen Fertilizers

Nitrate (NO₃ ^-) and ammonium (NH₄ ⁺) are the primary forms of inorganic nitrogen acquired by plant roots. LncRNAs, as key regulators of gene expression, are a class of non-coding RNAs larger than 200 bp. However, knowledge about the regulatory role of lncRNAs in response to different nitrogen forms remains limited, particularly in woody plants. Here, we performed strand-specific RNA-sequencing of P. × canescens roots under three different nitrogen fertilization treatments. In total, 324 lncRNAs and 6,112 mRNAs were identified as showing significantly differential expression between the NO₃ ^- and NH₄NO₃ treatments. Moreover, 333 lncRNAs and 6,007 mRNAs showed significantly differential expression between the NH₄ ⁺ and NH₄NO₃ treatments. Further analysis suggested that these lncRNAs and mRNAs have different response mechanisms for different nitrogen forms. In addition, functional annotation of cis and trans target mRNAs of differentially expressed lncRNAs indicated that 60 lncRNAs corresponding to 49 differentially expressed cis and trans target mRNAs were involved in plant nitrogen metabolism and amino acid biosynthesis and metabolism. Furthermore, 42 lncRNAs were identified as putative precursors of 63 miRNAs, and 28 differentially expressed lncRNAs were potential endogenous target mimics targeted by 96 miRNAs. Moreover, ceRNA regulation networks were constructed. MSTRG.6097.1, MSTRG.13550.1, MSTRG.2693.1, and MSTRG.12899.1, as hub lncRNAs in the ceRNA networks, are potential candidate lncRNAs for studying the regulatory mechanism in poplar roots under different nitrogen fertilization treatments. The results provide a basis for obtaining insight into the molecular mechanisms of lncRNA responses to different nitrogen forms in woody plants.

The long non-coding RNA landscape of Candida yeast pathogens

Long non-coding RNAs (lncRNAs) constitute a poorly studied class of transcripts with emerging roles in key cellular processes. Despite efforts to characterize lncRNAs across a wide range of species, these molecules remain largely unexplored in most eukaryotic microbes, including yeast pathogens of the Candida clade. Here, we analyze thousands of publicly available sequencing datasets to infer and characterize the lncRNA repertoires of five major Candida pathogens: Candida albicans, Candida tropicalis, Candida parapsilosis, Candida auris and Candida glabrata. Our results indicate that genomes of these species encode hundreds of lncRNAs that show levels of evolutionary constraint intermediate between those of intergenic genomic regions and protein-coding genes. Despite their low sequence conservation across the studied species, some lncRNAs are syntenic and are enriched in shared sequence motifs. We find co-expression of lncRNAs with certain protein-coding transcripts, hinting at potential functional associations. Finally, we identify lncRNAs that are differentially expressed during infection of human epithelial cells for four of the studied species. Our comprehensive bioinformatic analyses of Candida lncRNAs pave the way for future functional characterization of these transcripts.

Long non-coding RNAs (lncRNAs) play roles in key cellular processes, but remain largely unexplored in fungal pathogens such as Candida. Here, Hovhannisyan and Gabaldón analyze thousands of sequencing datasets to infer and characterize the lncRNA repertoires of five Candida species, paving the way for their future functional characterization.

Comparative Analysis of Coding and Non-Coding Features within Insect Tolerance Loci in Wheat with Their Homologs in Cereal Genomes

Food insecurity and malnutrition have reached critical levels with increased human population, climate fluctuations, water shortage; therefore, higher-yielding crops are in the spotlight of numerous studies. Abiotic factors affect the yield of staple food crops; among all, wheat stem sawfly (Cephus cinctus Norton) and orange wheat blossom midge (Sitodiplosis mosellana) are two of the most economically and agronomically harmful insect pests which cause yield loss in cereals, especially in wheat in North America. There is no effective strategy for suppressing this pest damage yet, and only the plants with intrinsic tolerance mechanisms such as solid stem phenotypes for WSS and antixenosis and/or antibiosis mechanisms for OWBM can limit damage. A major QTL and a causal gene for WSS resistance were previously identified in wheat, and 3 major QTLs and a causal gene for OWBM resistance. Here, we present a comparative analysis of coding and non-coding features of these loci of wheat across important cereal crops, barley, rye, oat, and rice. This research paves the way for our cloning and editing of additional WSS and OWBM tolerance gene(s), proteins, and metabolites.

Downregulation of lncRNA PpL-T31511 and Pp-miRn182 Promotes Hydrogen Cyanamide-Induced Endodormancy Release through the PP2C-H2O2 Pathway in Pear (Pyrus pyrifolia)

Bud endodormancy is an important, complex process subject to both genetic and epigenetic control, the mechanism of which is still unclear. The endogenous hormone abscisic acid (ABA) and its signaling pathway play important roles in the endodormancy process, in which the type 2C protein phosphatases (PP2Cs) is key to the ABA signal pathway. Due to its excellent effect on endodormancy release, hydrogen cyanamide (HC) treatment is considered an effective measure to study the mechanism of endodormancy release. In this study, RNA-Seq analysis was conducted on endodormant floral buds of pear (Pyrus pyrifolia) with HC treatment, and the HC-induced PP2C gene PpPP2C1 was identified. Next, software prediction, expression tests and transient assays revealed that lncRNA PpL-T31511-derived Pp-miRn182 targets PpPP2C1. The expression analysis showed that HC treatment upregulated the expression of PpPP2C1 and downregulated the expression of PpL-T31511 and Pp-miRn182. Moreover, HC treatment inhibited the accumulation of ABA signaling pathway-related genes and hydrogen peroxide (H₂O₂). Furthermore, overexpression of Pp-miRn182 reduced the inhibitory effect of PpPP2C1 on the H₂O₂ content. In summary, our study suggests that downregulation of PpL-T31511-derived Pp-miRn182 promotes HC-induced endodormancy release in pear plants through the PP2C-H₂O₂ pathway.

Advances in Non-Coding RNA Sequencing

Non-coding RNAs (ncRNAs) comprise a set of abundant and functionally diverse RNA molecules. Since the discovery of the first ncRNA in the 1960s, ncRNAs have been shown to be involved in nearly all steps of the central dogma of molecular biology. In recent years, the pace of discovery of novel ncRNAs and their cellular roles has been greatly accelerated by high-throughput sequencing. Advances in sequencing technology, library preparation protocols as well as computational biology helped to greatly expand our knowledge of which ncRNAs exist throughout the kingdoms of life. Moreover, RNA sequencing revealed crucial roles of many ncRNAs in human health and disease. In this review, we discuss the most recent methodological advancements in the rapidly evolving field of high-throughput sequencing and how it has greatly expanded our understanding of ncRNA biology across a large number of different organisms.