Long non-coding RNAs and their biological roles in plants

Unraveling the lncRNA-miRNA-mRNA Regulatory Network Involved in Poplar Coma Development through High-Throughput Sequencing

Poplar coma, the fluff-like appendages of seeds originating from the differentiated surface cells of the placenta and funicle, aids in the long-distance dispersal of seeds in the spring. However, it also poses hazards to human safety and causes pollution in the surrounding environment. Unraveling the regulatory mechanisms governing the initiation and development of coma is essential for addressing this issue comprehensively. In this study, strand-specific RNA-seq was conducted at three distinct stages of coma development, revealing 1888 lncRNAs and 52,810 mRNAs. The expression profiles of lncRNAs and mRNAs during coma development were analyzed. Subsequently, potential target genes of lncRNAs were predicted through co-localization and co-expression analyses. Integrating various types of sequencing data, lncRNA-miRNA-TF regulatory networks related to the initiation of coma were constructed. Utilizing identified differentially expressed genes encoding kinesin and actin, lncRNA-miRNA-mRNA regulatory networks associated with the construction and arrangement of the coma cytoskeleton were established. Additionally, relying on differentially expressed genes encoding cellulose synthase, sucrose synthase, and expansin, lncRNA-miRNA-mRNA regulatory networks related to coma cell wall synthesis and remodeling were developed. This study not only enhances the comprehension of lncRNA but also provides novel insights into the molecular mechanisms governing the initiation and development of poplar coma.

Long noncoding RNA from Betula platyphylla, BplncSIR1, confers salt tolerance by regulating BpNAC2 to mediate reactive oxygen species scavenging and stomatal movement

Long noncoding RNAs (lncRNAs) play an important role in abiotic stress tolerance. However, their function in conferring abiotic stress tolerance is still unclear. Herein, we characterized the function of a salt-responsive nuclear lncRNA (BplncSIR1) from Betula platyphylla (birch). Birch plants overexpressing and knocking out for BplncSIR1 were generated. BplncSIR1 was found to improve salt tolerance by inducing antioxidant activity and stomatal closure, and also accelerate plant growth. Chromatin isolation by RNA purification (ChIRP) combined with RNA sequencing indicated that BplncSIR1 binds to the promoter of BpNAC2 (encoding NAC domain-containing protein 2) to activate its expression. Plants overexpressing and knocking out for BpNAC2 were generated. Consistent with that of BplncSIR1, overexpression of BpNAC2 also accelerated plant growth and conferred salt tolerance. In addition, BpNAC2 binds to different cis-acting elements, such as G-box and 'CCAAT' sequences, to regulate the genes involved in salt tolerance, resulting in reduced ROS accumulation and decreased water loss rate by stomatal closure. Taken together, BplncSIR1 serves as the regulator of BpNAC2 to induce its expression in response to salt stress, and activated BpNAC2 accelerates plant growth and improves salt tolerance. Therefore, BplncSIR1 might be a candidate gene for molecular breeding to cultivate plants with both a high growth rate and improved salt tolerance.

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.

Beyond the genome: lncRNAs as regulators of the PI3K/AKT pathway in lung cancer

Lung cancer is a prevalent and devastating disease, representing a significant global health burden. Despite advancements in therapeutic strategies, the molecular mechanisms underlying its pathogenesis remain incompletely understood. Lung cancer typically displays the deregulated activity of the phosphoinositide 3-kinase/protein kinase B (PI3K/AKT) pathway, which is vital for cell proliferation, survival, and metastasis. Emerging evidence suggests that long non-coding RNA (lncRNAs) can modulate the PI3K/AKT pathway, offering new insights into lung cancer biology and potential therapeutic opportunities. These lncRNA act as either oncogenes, promoting pathway activation, or tumour suppressors, attenuating pathway signalling. The dysregulation of lncRNA is associated with various cellular processes, including apoptosis, cell cycle control, epithelial-mesenchymal transition (EMT), and angiogenesis, ultimately influencing lung cancer growth and metastasis. The development of novel therapeutic strategies, such as small interfering RNAs (siRNAs), antisense oligonucleotides, and CRISPR/Cas9-mediated gene editing, holds promise for restoring lncRNAs dysregulation and re-establishing the equilibrium of the PI3K/AKT pathway. The emerging role of lncRNAs as regulators of the PI3K/AKT pathway sheds new light on the complex molecular landscape of lung cancer. Understanding the interplay between lncRNA and the PI3K/AKT pathway could lead to the identification of novel biomarkers for prognosis and therapeutic targets for precision medicine. The potential of lncRNAs-based therapeutics may pave the way for more effective and personalized treatment approaches in lung cancer and potentially other malignancies with dysregulated PI3K/AKT signalling. This review aims to explore the emerging role of lncRNAs as key regulators of the PI3K/AKT pathway in lung cancer.

Engineering traits through CRISPR/cas genome editing in woody species to improve forest diversity and yield

Dangers confronting forest ecosystems are many and the strength of these biological systems is deteriorating, thus substantially affecting tree physiology, phenology, and growth. The establishment of genetically engineered trees into degraded woodlands, which would be adaptive to changing climate, could help in subsiding ecological threats and bring new prospects. This should not be resisted due to the apprehension of transgene dispersal in forests. Consequently, it is important to have a deep insight into the genetic structure and phenotypic limits of the reproductive capability of tree stands/population(s) to endure tolerance and survival. Importantly, for a better understanding of genes and their functional mechanisms, gene editing (GeEd) technology is an excellent molecular tool to unravel adaptation progressions. Therefore, GeEd could be harnessed for resolving the allelic interactions for the creation of gene diversity, and transgene dispersal may be alleviated among the population or species in different bioclimatic zones around the globe. This review highlights the potential of the CRISPR/Cas tools in genomic, transcriptomic, and epigenomic-based assorted and programmable alterations of genes in trees that might be able to fix the trait-specific gene function. Also, we have discussed the application of diverse forms of GeEd to genetically improve several traits, such as wood density, phytochemical constituents, biotic and abiotic stress tolerance, and photosynthetic efficiency in trees. We believe that the technology encourages fundamental research in the forestry sector besides addressing key aspects, which might fasten tree breeding and germplasm improvement programs worldwide.

LncRNA109897-JrCCR4-JrTLP1b forms a positive feedback loop to regulate walnut resistance against anthracnose caused by Colletotrichum gloeosporioides

Walnut anthracnose induced by Colletotrichum gloeosporioides is a disastrous disease that severely restricts the development of the walnut industry in China. Long non-coding RNAs (lncRNAs) are involved in adaptive responses to disease, but their roles in the regulation of walnut anthracnose resistance response are not well defined. In this study, transcriptome analysis demonstrated that a C. gloeosporioides-induced lncRNA, lncRNA109897, located upstream from the target gene JrCCR4, upregulated the expression of JrCCR4. JrCCR4 interacted with JrTLP1b and promoted its transcriptional activity. In turn, JrTLP1b induced the transcription of lncRNA109897 to promote its expression. Meanwhile, transient expression in walnut leaves and stable transformation of Arabidopsis thaliana further proved that lncRNA, JrCCR4, and JrTLP1b improve the resistance of C. gloeosporioides. Collectively, these findings provide insights into the mechanism by which the lncRNA109897-JrCCR4-JrTLP1b transcriptional cascade regulates the resistance of walnut to anthracnose.

Insights to Gossypium defense response against Verticillium dahliae: the Cotton Cancer

The soil-borne pathogen Verticillium dahliae, also referred as "The Cotton Cancer," is responsible for causing Verticillium wilt in cotton crops, a destructive disease with a global impact. To infect cotton plants, the pathogen employs multiple virulence mechanisms such as releasing enzymes that degrade cell walls, activating genes that contribute to virulence, and using protein effectors. Conversely, cotton plants have developed numerous defense mechanisms to combat the impact of V. dahliae. These include strengthening the cell wall by producing lignin and depositing callose, discharging reactive oxygen species, and amassing hormones related to defense. Despite the efforts to develop resistant cultivars, there is still no permanent solution to Verticillium wilt due to a limited understanding of the underlying molecular mechanisms that drive both resistance and pathogenesis is currently prevalent. To address this challenge, cutting-edge technologies such as clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9), host-induced gene silencing (HIGS), and gene delivery via nano-carriers could be employed as effective alternatives to control the disease. This article intends to present an overview of V. dahliae virulence mechanisms and discuss the different cotton defense mechanisms against Verticillium wilt, including morphophysiological and biochemical responses and signaling pathways including jasmonic acid (JA), salicylic acid (SA), ethylene (ET), and strigolactones (SLs). Additionally, the article highlights the significance of microRNAs (miRNAs), circular RNAs (circRNAs), and long non-coding RNAs (lncRNAs) in gene expression regulation, as well as the different methods employed to identify and functionally validate genes to achieve resistance against this disease. Gaining a more profound understanding of these mechanisms could potentially result in the creation of more efficient strategies for combating Verticillium wilt in cotton crops.

The coordinated regulation of early meiotic stages is dominated by non-coding RNAs and stage-specific transcription in wheat

Reproductive success hinges on precisely coordinated meiosis, yet our understanding of how structural rearrangements of chromatin and phase transitions during meiosis are transcriptionally regulated is limited. In crop plants, detailed analysis of the meiotic transcriptome could identify regulatory genes and epigenetic regulators that can be targeted to increase recombination rates and broaden genetic variation, as well as provide a resource for comparison among eukaryotes of different taxa to answer outstanding questions about meiosis. We conducted a meiotic stage-specific analysis of messenger RNA (mRNA), small non-coding RNA (sncRNA), and long intervening/intergenic non-coding RNA (lincRNA) in wheat (Triticum aestivum L.) and revealed novel mechanisms of meiotic transcriptional regulation and meiosis-specific transcripts. Amidst general repression of mRNA expression, significant enrichment of ncRNAs was identified during prophase I relative to vegetative cells. The core meiotic transcriptome was comprised of 9309 meiosis-specific transcripts, 48 134 previously unannotated meiotic transcripts, and many known and novel ncRNAs differentially expressed at specific stages. The abundant meiotic sncRNAs controlled the reprogramming of central metabolic pathways by targeting genes involved in photosynthesis, glycolysis, hormone biosynthesis, and cellular homeostasis, and lincRNAs enhanced the expression of nearby genes. Alternative splicing was not evident in this polyploid species, but isoforms were switched at phase transitions. The novel, stage-specific regulatory controls uncovered here challenge the conventional understanding of this crucial biological process and provide a new resource of requisite knowledge for those aiming to directly modulate meiosis to improve crop plants. The wheat meiosis transcriptome dataset can be queried for genes of interest using an eFP browser located at https://bar.utoronto.ca/efp_wheat/cgi-bin/efpWeb.cgi?dataSource=Wheat_Meiosis.

Small RNAs and Karma methylation in Elaeis guineensis mother palms are linked to high clonal mantling

The mantled phenotype is an abnormal somaclonal variant arising from the oil palm cloning process and severe phenotypes lead to oil yield losses. Hypomethylation of the Karma retrotransposon within the B-type MADS-box EgDEF1 gene has been associated with this phenotype. While abnormal Karma-EgDEF1 hypomethylation was detected in mantled clones, we examined the methylation state of Karma in ortets that gave rise to high mantling rates in their clones. Small RNAs (sRNAs) were proposed to play a role in Karma hypomethylation as part of the RNA-directed DNA methylation process, hence differential expression analysis of sRNAs between the ortet groups was conducted. While no sRNA was differentially expressed at the Karma-EgDEF1 region, three sRNA clusters were differentially regulated in high-mantling ortets. The first two down-regulated clusters were possibly derived from long non-coding RNAs while the third up-regulated cluster was derived from the intron of a DnaJ chaperone gene. Several predicted mRNA targets for the first two sRNA clusters conversely displayed increased expression in high-mantling relative to low-mantling ortets. These predicted mRNA targets may be associated with defense or pathogenesis response. In addition, several differentially methylated regions (DMRs) were identified in Karma and its surrounding regions, mainly comprising subtle CHH hypomethylation in high-mantling ortets. Four of the 12 DMRs were located in a region corresponding to hypomethylated areas at the 3'end of Karma previously reported in mantled clones. Further investigations on these sRNAs and DMRs may indicate the predisposition of certain ortets towards mantled somaclonal variation.

The long noncoding RNA LINC15957 regulates anthocyanin accumulation in radish

Radish (Raphanus sativus L.) is an important root vegetable crop belonging to the Brassicaceae family. Anthocyanin rich radish varieties are popular among consumers because of their bright color and high nutritional value. However, the underlying molecular mechanism responsible for skin and flesh induce anthocyanin biosynthesis in transient overexpression, gene silencing and transcriptome sequencing were used to verify its function in radish anthocyanin accumulation, radish remains unclear. Here, we identified a long noncoding RNA LINC15957, overexpression of LINC15957 was significantly increased anthocyanin accumulation in radish leaves, and the expression levels of structural genes related to anthocyanin biosynthesis were also significantly increased. Anthocyanin accumulation and expression levels of anthocyanin biosynthesis genes were significantly reduced in silenced LINC15957 flesh when compared with control. By the transcriptome sequencing of the overexpressed LINC15957 plants and the control, 5,772 differentially expressed genes were identified. A total of 3,849 differentially expressed transcription factors were identified, of which MYB, bHLH, WD40, bZIP, ERF, WRKY and MATE were detected and differentially expressed in the overexpressed LINC15957 plants. KEGG enrichment analysis revealed the genes were significant enriched in tyrosine, L-Phenylalanine, tryptophan, phenylpropanol, and flavonoid biosynthesis. RT-qPCR analysis showed that 8 differentially expressed genes (DEGs) were differentially expressed in LINC15957-overexpressed plants. These results suggested that LINC15957 involved in regulate anthocyanin accumulation and provide abundant data to investigate the genes regulate anthocyanin biosynthesis in radish.

Regulatory miPEP Open Reading Frames Contained in the Primary Transcripts of microRNAs

This review aims to consider retrospectively the available data on the coding properties of pri-microRNAs and the regulatory functions of their open reading frames (ORFs) and the encoded peptides (miPEPs). Studies identifying miPEPs and analyzing the fine molecular mechanisms of their functional activities are reviewed together with a brief description of the methods to identify pri-miRNA ORFs and the encoded protein products. Generally, miPEPs have been identified in many plant species of several families and in a few animal species. Importantly, molecular mechanisms of the miPEP action are often quite different between flowering plants and metazoan species. Requirement for the additional studies in these directions is highlighted by alternative findings concerning negative or positive regulation of pri-miRNA/miRNA expression by miPEPs in plants and animals. Additionally, the question of how miPEPs are distributed in non-flowering plant taxa is very important for understanding the evolutionary origin of such micropeptides. Evidently, further extensive studies are needed to explore the functions of miPEPs and the corresponding ORFs and to understand the full set of their roles in eukaryotic organisms. Thus, we address the most recent integrative views of different genomic, physiological, and molecular aspects concerning the expression of miPEPs and their possible fine functions.

Astrocytic MicroRNAs and Transcription Factors in Alzheimer's Disease and Therapeutic Interventions

Astrocytes are important for maintaining cholesterol metabolism, glutamate uptake, and neurotransmission. Indeed, inflammatory processes and neurodegeneration contribute to the altered morphology, gene expression, and function of astrocytes. Astrocytes, in collaboration with numerous microRNAs, regulate brain cholesterol levels as well as glutamatergic and inflammatory signaling, all of which contribute to general brain homeostasis. Neural electrical activity, synaptic plasticity processes, learning, and memory are dependent on the astrocyte-neuron crosstalk. Here, we review the involvement of astrocytic microRNAs that potentially regulate cholesterol metabolism, glutamate uptake, and inflammation in Alzheimer's disease (AD). The interaction between astrocytic microRNAs and long non-coding RNA and transcription factors specific to astrocytes also contributes to the pathogenesis of AD. Thus, astrocytic microRNAs arise as a promising target, as AD conditions are a worldwide public health problem. This review examines novel therapeutic strategies to target astrocyte dysfunction in AD, such as lipid nanodiscs, engineered G protein-coupled receptors, extracellular vesicles, and nanoparticles.

A Feedback Loop of LINC00665 and the Wnt Signaling Pathway Expedites Osteosarcoma Cell Proliferation, Invasion, and Epithelial-Mesenchymal Transition

OBJECTIVES

Osteosarcoma (OS) is a malignant tumor with frequent occurrence among teenagers. Long non-coding RNAs (lncRNAs) play pro-cancer roles in many tumors. The purpose of this study was to figure out the functional role of a novel lncRNA long intergenic non-protein coding RNA 665 (LINC00665) in OS by observing the OS cell behaviors.

METHODS

Quantitative reverse transcription polymerase chain reaction (RT-qPCR) was used to analyze LINC00665 expression in OS cells. Cell function assays assessed the impacts of LINC00665 on OS cell phenotype. Immunofluorescence and western blot analyzed the function of LINC00665 on epithelial-mesenchymal transition (EMT) in OS. Moreover, mechanistic assays analyzed the downstream mechanism of LINC00665 in OS cells.

RESULTS

LINC00665 was significantly up-regulated in OS cells. LINC00665 silence facilitated OS cell proliferation, migration, invasion, and EMT while inhibiting cell apoptosis. Mechanically, LINC00665 acted as a competing endogenous RNA (ceRNA) to sponge miR-1249-5p and thereby modulated Wnt family member 2B (WNT2B) to activate Wnt pathway. Wnt pathway activated LINC00665 expression transcriptionally.

CONCLUSIONS

Our study uncovered the cancer-promoting role of LINC00665 in OS, and the feedback loop of LINC00665/miR-1249-5p/WNT2B/Wnt might be a potential target for OS treatment.

Regulatory non-coding RNAs: Emerging roles during plant cell reprogramming and in vitro regeneration

Plant regeneration is a well-known capacity of plants occurring either in vivo or in vitro. This potential is the basis for plant micropropagation and genetic transformation as well as a useful system to analyse different aspects of plant development. Recent studies have proven that RNA species with no protein-coding capacity are key regulators of cellular function and essential for cell reprogramming. In this review, the current knowledge on the role of several ncRNAs in plant regeneration processes is summarized, with a focus on cell fate reprogramming. Moreover, the involvement/impact of microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and small-interfering RNAs (siRNAs) in the regulatory networks of cell dedifferentiation, proliferation and differentiation is also analysed. A deeper understanding of plant ncRNAs in somatic cell reprogramming will allow a better modulation of in vitro regeneration processes such as organogenesis and somatic embryogenesis.

PpTCP18 is upregulated by lncRNA5 and controls branch number in peach (Prunus persica) through positive feedback regulation of strigolactone biosynthesis

Branch number is an important agronomic trait in peach (Prunus persica) trees because plant architecture affects fruit yield and quality. Although breeders can select varieties with different tree architecture, the biological mechanisms underlying architecture remain largely unclear. In this study, a pillar peach ('Zhaoshouhong') and a standard peach ('Okubo') were compared. 'Zhaoshouhong' was found to have significantly fewer secondary branches than 'Okubo'. Treatment with the synthetic strigolactone (SL) GR24 decreased branch number. Transcriptome analysis indicated that PpTCP18 (a homologous gene of Arabidopsis thaliana BRC1) expression was negatively correlated with strigolactone synthesis gene expression, indicating that PpTCP18 may play an important role in peach branching. Yeast one-hybrid, electrophoretic mobility shift, dual-luciferase assays and PpTCP18-knockdown in peach leaf buds indicated that PpTCP18 could increase expression of PpLBO1, PpMAX1, and PpMAX4. Furthermore, transgenic Arabidopsis plants overexpressing PpTCP18 clearly exhibited reduced primary rosette-leaf branches. Moreover, lncRNA sequencing and transient expression analysis revealed that lncRNA5 targeted PpTCP18, significantly increasing PpTCP18 expression. These results provide insights into the mRNA and lncRNA network in the peach SL signaling pathway and indicate that PpTCP18, a transcription factor downstream of SL signaling, is involved in positive feedback regulation of SL biosynthesis. This role of PpTCP18 may represent a novel mechanism in peach branching regulation. Our study improves current understanding of the mechanisms underlying peach branching and provides theoretical support for genetic improvement of peach tree architecture.

Population transcriptomic analysis identifies the comprehensive lncRNAs landscape of spike in wheat (Triticum aestivum L.)

BACKGROUND

Long noncoding RNAs (lncRNAs) are emerging as the important regulators involving in growth and development as well as stress response in plants. However, current lncRNA studies were mainly performed at the individual level and the significance of it is not well understood in wheat.

RESULTS

In this study, the lncRNA landscape of wheat spike was characterized through analysing a total of 186 spike RNA-seq datasets from 93 wheat genotypes. A total of 35,913 lncRNAs as well as 1,619 lncRNA-mRNA pairs comprised of 443 lncRNAs and 464 mRNAs were obtained. Compared to coding genes, these lncRNAs displayed rather low conservation among wheat and other gramineous species. Based on re-sequencing data, the genetic variations of these lncRNA were investigated and obvious genetic bottleneck were found on them during wheat domestication process. Furthermore, 122 lncRNAs were found to act as ceRNA to regulate endogenous competition. Finally, association and co-localization analysis of the candidate lncRNA-mRNA pairs identified 170 lncRNAs and 167 target mRNAs significantly associated with spike-related traits, including lncRNA.127690.1/TraesCS2A02G518500.1 (PMEI) and lncRNA.104854.1/TraesCS6A02G050300.1 (ATG5) associated with heading date and spike length, respectively.

CONCLUSIONS

This study reported the lncRNA landscape of wheat spike through the population transcriptome analysis, which not only contribute to better understand the wheat evolution from the perspective of lncRNA, but also lay the foundation for revealing roles of lncRNA playing in spike development.

Genome-wide analysis suggests the potential role of lncRNAs during seed development and seed size/weight determination in chickpea

The integrated transcriptome data analyses suggested the plausible roles of lncRNAs during seed development in chickpea. The candidate lncRNAs associated with QTLs and those involved in miRNA-mediated seed size/weight determination in chickpea have been identified. Long non-coding RNAs (lncRNAs) are important regulators of various biological processes. Here, we identified lncRNAs at seven successive stages of seed development in small-seeded and large-seeded chickpea cultivars. In total, 4751 lncRNAs implicated in diverse biological processes were identified. Most of lncRNAs were conserved between the two cultivars, whereas only a few of them were conserved in other plants, suggesting their species-specificity. A large number of lncRNAs differentially expressed between the two chickpea cultivars associated with seed development-related processes were identified. The lncRNAs acting as precursors of miRNAs and those mimicking target protein-coding genes of miRNAs involved in seed size/weight determination, including HAIKU1, BIG SEEDS1, and SHB1, were also revealed. Further, lncRNAs located within seed size/weight associated quantitative trait loci were also detected. Overall, we present a comprehensive resource and identified candidate lncRNAs that may play important roles during seed development and seed size/weight determination in chickpea.

WRKY63 transcriptional activation of COOLAIR and COLDAIR regulates vernalization-induced flowering

Arabidopsis (Arabidopsis thaliana) FLOWERING LOCUS C (FLC) acts as a key flowering regulator by repressing the expression of the floral integrator FLOWERING LOCUS T (FT). Prolonged exposure to cold (vernalization) induces flowering by reducing FLC expression. The long noncoding RNAs (lncRNAs) COOLAIR and COLDAIR, which are transcribed from the 3' end and the first intron of FLC, respectively, are important for FLC repression under vernalization. However, the molecular mechanism of how COOLAIR and COLDAIR are transcriptionally activated remains elusive. In this study, we found that the group-III WRKY transcription factor WRKY63 can directly activate FLC. wrky63 mutant plants display an early flowering phenotype and are insensitive to vernalization. Interestingly, we found that WRKY63 can activate the expression of COOLAIR and COLDAIR by binding to their promoters.WRKY63 therefore acts as a dual regulator that activates FLC directly under non-vernalization conditions but represses FLC indirectly during vernalization through inducing COOLAIR and COLDAIR. Furthermore, genome-wide occupancy profile analyses indicated that the binding of WRKY63 to vernalization-induced genes increases after vernalization. In addition, WRKY63 binding is associated with decreased levels of the repressive marker Histone H3 Lysine 27 trimethylation (H3K27me3). Collectively, our results indicate that WRKY63 is an important flowering regulator involved in vernalization-induced transcriptional regulation.

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.

Integrative mRNA and Long Noncoding RNA Analysis Reveals the Regulatory Network of Floral Bud Induction in Longan (Dimocarpus longan Lour.)

Longan (Dimocarpus longan Lour.) is a tropical/subtropical fruit tree of significant economic importance. Floral induction is an essential process for longan flowering and plays decisive effects on the longan yield. Due to the instability of flowering, it is necessary to understand the molecular mechanisms of floral induction in longan. In this study, mRNA and long noncoding RNA (lncRNA) transcriptome sequencing were performed using the apical buds of fruiting branches as materials. A total of 7,221 differential expressions of mRNAs (DEmRNAs) and 3,238 differential expressions of lncRNAs (DElncRNAs) were identified, respectively. KEGG enrichment analysis of DEmRNAs highlighted the importance of starch and sucrose metabolic, circadian rhythms, and plant hormone signal transduction pathways during floral induction. Combining the analysis of weighted gene co-expression network (WGCNA) and expression pattern of DEmRNAs in the three pathways, specific transcriptional characteristics at each stage during floral induction and regulatory network involving co-expressed genes were investigated. The results showed that sucrose metabolism and auxin signal transduction may be crucial for the growth and maturity of autumn shoots in September and October (B1-B2 stage); starch and sucrose metabolic, circadian rhythms, and plant hormone signal transduction pathways participated in the regulation of floral bud physiological differentiation together in November and December (B3-B4 stage) and the crosstalk among three pathways was also found. Hub genes in the co-expression network and key DEmRNAs in three pathways were identified. The circadian rhythm genes FKF1 and GI were found to activate SOC1gene through the photoperiod core factor COL genes, and they were co-expressed with auxin, gibberellin, abscisic acid, ethylene signaling genes, and sucrose biosynthesis genes at B4 stage. A total of 12 hub-DElncRNAs had potential for positively affecting their distant target genes in three putative key pathways, predominantly in a co-transcriptional manner. A hypothetical model of regulatory pathways and key genes and lncRNAs during floral bud induction in longan was proposed finally. Our studies will provide valuable clues and information to help elucidate the potential molecular mechanisms of floral initiation in longan and woody fruit trees.

How do plants remember drought?

Plants develop both short-term and transgenerational memory of drought stress through epigenetic regulation of transcription for a better response to subsequent exposure. Recurrent spells of droughts are more common than a single drought, with intermittent moist recovery intervals. While the detrimental effects of the first drought on plant structure and physiology are unavoidable, if survived, plants can memorize the first drought to present a more robust response to the following droughts. This includes a partial stomatal opening in the watered recovery interval, higher levels of osmoprotectants and ABA, and attenuation of photosynthesis in the subsequent exposure. Short-term drought memory is regulated by ABA and other phytohormone signaling with transcriptional memory behavior in various genes. High levels of methylated histones are deposited at the drought-tolerance genes. During the recovery interval, the RNA polymerase is stalled to be activated by a pause-breaking factor in the subsequent drought. Drought leads to DNA demethylation near drought-response genes, with genetic control of the process. Progenies of the drought-exposed plants can better adapt to drought owing to the inheritance of particular methylation patterns. However, a prolonged watered recovery interval leads to loss of drought memory, mediated by certain demethylases and chromatin accessibility factors. Small RNAs act as critical regulators of drought memory by altering transcript levels of drought-responsive target genes. Further studies in the future will throw more light on the genetic control of drought memory and the interplay of genetic and epigenetic factors in its inheritance. Plants from extreme environments can give queues to understanding robust memory responses at the ecosystem level.

Whole-Transcriptome Analysis Reveals Autophagy Is Involved in Early Senescence of zj-es Mutant Rice

Senescence is a necessary stage of plant growth and development, and the early senescence of rice will lead to yield reduction and quality decline. However, the mechanisms of rice senescence remain obscure. In this study, we characterized an early-senescence rice mutant, designated zj-es (ZheJing-early senescence), which was derived from the japonica rice cultivar Zhejing22. The mutant zj-es exhibited obvious early-senescence phenotype, such as collapsed chloroplast, lesions in leaves, declined fertility, plant dwarf, and decreased agronomic traits. The ZJ-ES gene was mapped in a 458 kb-interval between the molecular markers RM5992 and RM5813 on Chromosome 3, and analysis suggested that ZJ-ES is a novel gene controlling rice early senescence. Subsequently, whole-transcriptome RNA sequencing was performed on zj-es and its wild-type rice to dissect the underlying molecular mechanism for early senescence. Totally, 10,085 differentially expressed mRNAs (DEmRNAs), 1,253 differentially expressed lncRNAs (DElncRNAs), and 614 differentially expressed miRNAs (DEmiRNAs) were identified, respectively, in different comparison groups. Based on the weighted gene co-expression network analysis (WGCNA), the co-expression turquoise module was found to be the key for the occurrence of rice early senescence. Furthermore, analysis on the competing endogenous RNA (CeRNA) network revealed that 14 lncRNAs possibly regulated 16 co-expressed mRNAs through 8 miRNAs, and enrichment analysis showed that most of the DEmRNAs and the targets of DElncRNAs and DEmiRNAs were involved in reactive oxygen species (ROS)-triggered autophagy-related pathways. Further analysis showed that, in zj-es, ROS-related enzyme activities were markedly changed, ROS were largely accumulated, autophagosomes were obviously observed, cell death was significantly detected, and lesions were notably appeared in leaves. Totally, combining our results here and the remaining research, we infer that ROS-triggered autophagy induces the programmed cell death (PCD) and its coupled early senescence in zj-es mutant rice.

Revealing the novel complexity of plant long non-coding RNA by strand-specific and whole transcriptome sequencing for evolutionarily representative plant species

Background

Previous studies on plant long noncoding RNAs (lncRNAs) lacked consistency and suffered from many factors like heterogeneous data sources and experimental protocols, different plant tissues, inconsistent bioinformatics pipelines, etc. For example, the sequencing of RNAs with poly(A) tails excluded a large portion of lncRNAs without poly(A), and use of regular RNA-sequencing technique did not distinguish transcripts’ direction for lncRNAs. The current study was designed to systematically discover and analyze lncRNAs across eight evolutionarily representative plant species, using strand-specific (directional) and whole transcriptome sequencing (RiboMinus) technique.

Results

A total of 39,945 lncRNAs (25,350 lincRNAs and 14,595 lncNATs) were identified, which showed molecular features of lncRNAs that are consistent across divergent plant species but different from those of mRNA. Further, transposable elements (TEs) were found to play key roles in the origination of lncRNA, as significantly large number of lncRNAs were found to contain TEs in gene body and promoter region, and transcription of many lncRNAs was driven by TE promoters. The lncRNA sequences were divergent even in closely related species, and most plant lncRNAs were genus/species-specific, amid rapid turnover in evolution. Evaluated with PhastCons scores, plant lncRNAs showed similar conservation level to that of intergenic sequences, suggesting that most lincRNAs were young and with short evolutionary age. INDUCED BY PHOSPHATE STARVATION (IPS) was found so far to be the only plant lncRNA group with conserved motifs, which may play important roles in the adaptation of terrestrial life during migration from aquatic to terrestrial. Most highly and specially expressed lncRNAs formed co-expression network with coding genes, and their functions were believed to be closely related to their co-expression genes.

Conclusion

The study revealed novel features and complexity of lncRNAs in plants through systematic analysis, providing important insights into the origination and evolution of plant lncRNAs.

Supplementary Information

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

A long noncoding RNA functions in high-light-induced anthocyanin accumulation in apple by activating ethylene synthesis

Anthocyanin production in apple (Malus domestica) fruit and their consequent coloration can be induced by high-light treatment. The hormone ethylene is also essential for this coloration, but the regulatory relationships that link ethylene and light with anthocyanin-associated coloration are not well defined. In this study, we observed that high-light treatment of apple fruit increased anthocyanin accumulation more than moderate-light treatment did and was the main contributor of induced ethylene production and activation of anthocyanin biosynthesis. A transcriptome study of light-treated apple fruit suggested that a long noncoding RNA (lncRNA), MdLNC610, the corresponding gene of which is physically located downstream from the 1-aminocyclopropane-1-carboxylate oxygenase (ACO) ethylene biosynthesis gene MdACO1, likely affects anthocyanin biosynthesis under high-light treatment. Expression and promoter β-glucuronidase reporter analyses further showed that MdLNC610 upregulates expression of MdACO1 and so likely participates in high-light-induced ethylene biosynthesis. Overexpression of MdACO1 and MdLNC610 in apple fruit and calli indicated that a major increase in MdLNC610 expression activates MdACO1 expression, thereby causing an increase in ethylene production and anthocyanin levels. These results suggest that MdLNC610 participates in the regulation of high-light-induced anthocyanin production by functioning as a positive regulator to promote MdACO1 expression and ethylene biosynthesis. Our study provides insights into the relationship between mRNA and lncRNA networks in the ethylene biosynthetic pathway and anthocyanin accumulation in apple fruit.

Deciphering the role of miRNA in reprogramming plant responses to drought stress

Drought is the most prevalent environmental stress that affects plants' growth, development, and crop productivity. However, plants have evolved adaptive mechanisms to respond to the harmful effects of drought. They reprogram their: transcriptome, proteome, and metabolome that alter their cellular and physiological processes and establish cellular homeostasis. One of the crucial regulatory processes that govern this reprogramming is post-transcriptional regulation by microRNAs (miRNAs). miRNAs are small non-coding RNAs, involved in the downregulation of the target mRNA via translation inhibition/mRNA degradation/miRNA-mediated mRNA decay/ribosome drop off/DNA methylation. Many drought-inducible miRNAs have been identified and characterized in plants. Their main targets are regulatory genes that influence growth, development, osmotic stress tolerance, antioxidant defense, phytohormone-mediated signaling, and delayed senescence during drought stress. Overexpression of drought-responsive miRNAs (Osa-miR535, miR160, miR408, Osa-miR393, Osa-miR319, and Gma-miR394) in certain plants has led to tolerance against drought stress indicating their vital role in stress mitigation. Similarly, knock down (miR166/miR398c) or deletion (miR169 and miR827) of miRNAs has also resulted in tolerance to drought stress. Likewise, engineered Arabidopsis plants with miR165, miR166 using short tandem target mimic strategy, exhibited drought tolerance. Since miRNAs regulate the expression of an array of drought-responsive genes, they can act as prospective targets for genetic manipulations to enhance drought tolerance in crops and achieve sustainable agriculture. Further investigations toward functional characterization of diverse miRNAs, and understanding stress-responses regulated by these miRNAs and their utilization in biotechnological applications is highly recommended.

Biogenesis, Functions, Interactions, and Resources of Non-Coding RNAs in Plants

Plant transcriptomes encompass a large number of functional non-coding RNAs (ncRNAs), only some of which have protein-coding capacity. Since their initial discovery, ncRNAs have been classified into two broad categories based on their biogenesis and mechanisms of action, housekeeping ncRNAs and regulatory ncRNAs. With advances in RNA sequencing technology and computational methods, bioinformatics resources continue to emerge and update rapidly, including workflow for in silico ncRNA analysis, up-to-date platforms, databases, and tools dedicated to ncRNA identification and functional annotation. In this review, we aim to describe the biogenesis, biological functions, and interactions with DNA, RNA, protein, and microorganism of five major regulatory ncRNAs (miRNA, siRNA, tsRNA, circRNA, lncRNA) in plants. Then, we systematically summarize tools for analysis and prediction of plant ncRNAs, as well as databases. Furthermore, we discuss the silico analysis process of these ncRNAs and present a protocol for step-by-step computational analysis of ncRNAs. In general, this review will help researchers better understand the world of ncRNAs at multiple levels.

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.

A SA-regulated lincRNA promotes Arabidopsis disease resistance by modulating pre-rRNA processing

Regulation of gene expression at translational level has been shown critical for plant defense against pathogen infection. Pre-rRNA processing is essential for ribosome biosynthesis and thus affects protein translation. It remains unknown if plants modulate pre-rRNA processing as a translation regulatory mechanism for disease resistance. In this study, we show a 5' snoRNA capped and 3' polyadenylated (SPA) lincRNA named SUNA1 promotes disease resistance involved in modulating pre-rRNA processing in Arabidopsis. SUNA1 expression is highly induced by Pst DC3000 infection, which is impaired in SA biosynthesis-defective mutant sid2 and signaling mutant npr1. Consistently, SA triggers SUNA1 expression dependent on NPR1. Functional analysis indicates that SUNA1 plays a positive role in Arabidopsis defense against Pst DC3000 relying on its snoRNA signature motifs. Potential mechanism study suggests that the nucleus-localized SUNA1 interacts with the nucleolar methyltransferase fibrillarin to modulate SA-controlled pre-rRNA processing, then enhancing the translational efficiency (TE) of some defense genes in Arabidopsis response to Pst DC3000 infection. NPR1 appears to have similar effects as SUNA1 on pre-rRNA processing and TE of defense genes. Together, these studies reveal one kind of undescribed antibacterial translation regulatory mechanism, in which SA-NPR1-SUNA1 signaling cascade controls pre-rRNA processing and TE of certain defense genes in Arabidopsis.

Long Intergenic Noncoding RNA (lincRNA) Discovery from Non-Strand-Specific RNA-Seq Data

Long noncoding RNAs (lncRNAs) are transcripts over 200 base pairs in length without discernible protein coding potential. Long intergenic noncoding RNAs (lincRNAs) constitute a subset of lncRNAs, which do not overlap protein coding genes. Here we describe a detailed pipeline for lincRNA discovery from publicly available non-stranded RNA-Seq datasets. The pipeline presented can be applied to any plant species for which RNA-Seq data and a reference genome sequence are available.

Transcriptional Association between mRNAs and Their Paired Natural Antisense Transcripts Following Fusarium oxysporum Inoculation in Brassica rapa L

Long noncoding RNAs (lncRNAs) play important roles in abiotic and biotic stress responses; however, studies on the mechanism of regulation of lncRNA expression are limited in plants. The present study examined the relationship between lncRNA expression level and two active histone modifications (H3K4me3 and H3K36me3) in Brassica rapa. Both histone marks were enriched in the chromatin regions encoding lncRNAs, especially around the transcription start site. The transcription level of long intergenic noncoding RNAs was positively associated with the level of H3K4me3 and H3K36me3, while this association was not observed in natural antisense RNAs (NATs) and intronic noncoding RNAs. As coordinate expression of mRNAs and paired NATs under biotic stress treatment has been identified, the transcriptional relationship between mRNAs and their paired NATs following Fusarium oxysporum f. sp. conglutinans (Foc) inoculation was examined. A positive association of expression levels between mRNAs and their paired NATs following Foc inoculation was observed. This association held for several defense-response-related genes and their NAT pairs. These results suggest that coordinate expression of mRNAs and paired NATs plays a role in the defense response against Foc.

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.

The parent-of-origin lncRNA MISSEN regulates rice endosperm development

The cereal endosperm is a major factor determining seed size and shape. However, the molecular mechanisms of endosperm development are not fully understood. Long noncoding RNAs (lncRNAs) function in various biological processes. Here we show a lncRNA, MISSEN, that plays an essential role in early endosperm development in rice (Oryza sativa). MISSEN is a parent-of-origin lncRNA expressed in endosperm, and negatively regulates endosperm development, leading to a prominent dent and bulge in the seed. Mechanistically, MISSEN functions through hijacking a helicase family protein (HeFP) to regulate tubulin function during endosperm nucleus division and endosperm cellularization, resulting in abnormal cytoskeletal polymerization. Finally, we revealed that the expression of MISSEN is inhibited by histone H3 lysine 27 trimethylation (H3K27me3) modification after pollination. Therefore, MISSEN is the first lncRNA identified as a regulator in endosperm development, highlighting the potential applications in rice breeding.

Noncoding-RNA-Mediated Regulation in Response to Macronutrient Stress in Plants

Macronutrient elements including nitrogen (N), phosphorus (P), potassium (K), calcium (Ca), magnesium (Mg), and sulfur (S) are required in relatively large and steady amounts for plant growth and development. Deficient or excessive supply of macronutrients from external environments may trigger a series of plant responses at phenotypic and molecular levels during the entire life cycle. Among the intertwined molecular networks underlying plant responses to macronutrient stress, noncoding RNAs (ncRNAs), mainly microRNAs (miRNAs) and long ncRNAs (lncRNAs), may serve as pivotal regulators for the coordination between nutrient supply and plant demand, while the responsive ncRNA-target module and the interactive mechanism vary among elements and species. Towards a comprehensive identification and functional characterization of nutrient-responsive ncRNAs and their downstream molecules, high-throughput sequencing has produced massive omics data for comparative expression profiling as a first step. In this review, we highlight the recent findings of ncRNA-mediated regulation in response to macronutrient stress, with special emphasis on the large-scale sequencing efforts for screening out candidate nutrient-responsive ncRNAs in plants, and discuss potential improvements in theoretical study to provide better guidance for crop breeding practices.

Long noncoding RNA lncRNA354 functions as a competing endogenous RNA of miR160b to regulate ARF genes in response to salt stress in upland cotton

Long non-coding RNAs (lncRNAs) play important roles in response to biotic and abiotic stress through acting as competing endogenous RNAs (ceRNAs) to decoy mature miRNAs. However, whether this mechanism is involved in cotton salt stress response remains unknown. We report the characterization of an endogenous lncRNA, lncRNA354, whose expression was reduced in salt-treated cotton and was localized at the nucleus and cytoplasm. Using endogenous target mimic (eTM) analysis, we predicted that lncRNA354 had a potential binding site for miR160b. Transient expression in tobacco demonstrated that lncRNA354 was a miR160b eTM and attenuated miR160b suppression of its target genes, including auxin response factors (ARFs). Silencing or overexpressing lncRNA354 affected the expression of miR160b and target ARFs. Silencing lncRNA354 and targets GhARF17/18 resulted in taller cotton plants and enhanced the resistant to salt stress. Overexpression of lncRNA354 and targets GhARF17/18 in Arabidopsis led to dwarf plants, decreased root dry weight and reduced salt tolerance. Our results show that the lncRNA354-miR160b effect on GhARF17/18 expression may modulate auxin signalling and thus affect growth. These results also shed new light on a mechanism of lncRNA-associated responses to salt stress.

Global identification of long non-coding RNAs involved in the induction of spinach flowering

Background

Spinach is a beneficial annual vegetable species and sensitive to the bolting or early flowering, which causes a large reduction in quality and productivity. Indeed, bolting is an event induced by the coordinated effects of various environmental factors and endogenous genetic components. Although some key flowering responsive genes have been identified in spinach, non-coding RNA molecules like long non-coding RNAs (lncRNAs) were not investigated yet. Herein, we used bioinformatic approaches to analyze the transcriptome datasets from two different accessions Viroflay and Kashan at two vegetative and reproductive stages to reveal novel lncRNAs and the construction of the lncRNA-mRNA co-expression network. Additionally, correlations among gene expression modules and phenotypic traits were investigated; day to flowering was chosen as our interesting trait.

Results

In the present study, we identified a total of 1141 lncRNAs, of which 111 were differentially expressed between vegetative and reproductive stages. The GO and KEGG analyses carried out on the cis target gene of lncRNAs showed that the lncRNAs play an important role in the regulation of flowering spinach. Network analysis pinpointed several well-known flowering-related genes such as ELF, COL1, FLT, and FPF1 and also some putative TFs like MYB, WRKY, GATA, and MADS-box that are important regulators of flowering in spinach and could be potential targets for lncRNAs.

Conclusions

This study is the first report on identifying bolting and flowering-related lncRNAs based on transcriptome sequencing in spinach, which provides a useful resource for future functional genomics studies, genes expression researches, evaluating genes regulatory networks and molecular breeding programs in the regulation of the genetic mechanisms related to bolting in spinach.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07989-1.

RNA-Sequencing Based microRNA Expression Signature of Colorectal Cancer: The Impact of Oncogenic Targets Regulated by miR-490-3p

To elucidate novel aspects of the molecular pathogenesis of colorectal cancer (CRC), we have created a new microRNA (miRNA) expression signature based on RNA-sequencing. Analysis of the signature showed that 84 miRNAs were upregulated, and 70 were downregulated in CRC tissues. Interestingly, our signature indicated that both guide and passenger strands of some miRNAs were significantly dysregulated in CRC tissues. These findings support our earlier data demonstrating the involvement of miRNA passenger strands in cancer pathogenesis. Our study focused on downregulated miR-490-3p and investigated its tumor-suppressive function in CRC cells. We successfully identified a total of 38 putative oncogenic targets regulated by miR-490-3p in CRC cells. Among these targets, the expression of three genes (IRAK1: p = 0.0427, FUT1: p = 0.0468, and GPRIN2: p = 0.0080) significantly predicted 5-year overall survival of CRC patients. Moreover, we analyzed the direct regulation of IRAK1 by miR-490-3p, and its resultant oncogenic function in CRC cells. Thus, we have clarified a part of the molecular pathway of CRC based on the action of tumor-suppressive miR-490-3p. This new miRNA expression signature of CRC will be a useful tool for elucidating new molecular pathogenesis in this disease.

Comparative Transcriptome Analysis Provides Insights into Response of Ulva compressa to Fluctuating Salinity Conditions

Ulva compressa, a green tide-forming species, can adapt to hypo-salinity conditions, such as estuaries and brackish lakes. To understand the underlying molecular mechanisms of hypo-salinity stress tolerance, transcriptome-wide gene expression profiles in U. compressa were created using digital gene expression profiles. The RNA-seq data were analyzed based on the comparison of differently expressed genes involved in specific pathways under hypo-salinity and recovery conditions. The up-regulation of genes in photosynthesis and glycolysis pathways may contribute to the recovery of photosynthesis and energy metabolism, which could provide sufficient energy for the tolerance under long-term hyposaline stress. Multiple strategies, such as ion transportation and osmolytes metabolism, were performed to maintain the osmotic homeostasis. Additionally, several long noncoding RNA were differently expressed during the stress, which could play important roles in the osmotolerance. Our work will serve as an essential foundation for the understanding of the tolerance mechanism of U. compressa under the fluctuating salinity conditions.

Advances in Cereal Crop Genomics for Resilience under Climate Change

Adapting to climate change, providing sufficient human food and nutritional needs, and securing sufficient energy supplies will call for a radical transformation from the current conventional adaptation approaches to more broad-based and transformative alternatives. This entails diversifying the agricultural system and boosting productivity of major cereal crops through development of climate-resilient cultivars that can sustainably maintain higher yields under climate change conditions, expanding our focus to crop wild relatives, and better exploitation of underutilized crop species. This is facilitated by the recent developments in plant genomics, such as advances in genome sequencing, assembly, and annotation, as well as gene editing technologies, which have increased the availability of high-quality reference genomes for various model and non-model plant species. This has necessitated genomics-assisted breeding of crops, including underutilized species, consequently broadening genetic variation of the available germplasm; improving the discovery of novel alleles controlling important agronomic traits; and enhancing creation of new crop cultivars with improved tolerance to biotic and abiotic stresses and superior nutritive quality. Here, therefore, we summarize these recent developments in plant genomics and their application, with particular reference to cereal crops (including underutilized species). Particularly, we discuss genome sequencing approaches, quantitative trait loci (QTL) mapping and genome-wide association (GWAS) studies, directed mutagenesis, plant non-coding RNAs, precise gene editing technologies such as CRISPR-Cas9, and complementation of crop genotyping by crop phenotyping. We then conclude by providing an outlook that, as we step into the future, high-throughput phenotyping, pan-genomics, transposable elements analysis, and machine learning hold much promise for crop improvements related to climate resilience and nutritional superiority.

Characterization of full-length transcriptome in Saccharum officinarum and molecular insights into tiller development

BACKGROUND

Although extensive breeding efforts are ongoing in sugarcane (Saccharum officinarum L.), the average yield is far below the theoretical potential. Tillering is an important component of sugarcane yield, however, the molecular mechanism underlying tiller development is still elusive. The limited genomic data in sugarcane, particularly due to its complex and large genome, has hindered in-depth molecular studies.

RESULTS

Herein, we generated full-length (FL) transcriptome from developing leaf and tiller bud samples based on PacBio Iso-Seq. In addition, we performed RNA-seq from tiller bud samples at three developmental stages (T0, T1 and T2) to uncover key genes and biological pathways involved in sugarcane tiller development. In total, 30,360 and 20,088 high-quality non-redundant isoforms were identified in leaf and tiller bud samples, respectively, representing 41,109 unique isoforms in sugarcane. Likewise, we identified 1063 and 1037 alternative splicing events identified in leaf and tiller bud samples, respectively. We predicted the presence of coding sequence for 40,343 isoforms, 98% of which was successfully annotated. Comparison with previous FL transcriptomes in sugarcane revealed 2963 unreported isoforms. In addition, we characterized 14,946 SSRs from 11,700 transcripts and 310 lncRNAs. By integrating RNA-seq with the FL transcriptome, 468 and 57 differentially expressed genes (DEG) were identified in T1vsT0 and T2vsT0, respectively. Strong up-regulation of several pyruvate phosphate dikinase and phosphoenolpyruvate carboxylase genes suggests enhanced carbon fixation and protein synthesis to facilitate tiller growth. Similarly, up-regulation of linoleate 9S-lipoxygenase and lipoxygenase genes in the linoleic acid metabolism pathway suggests high synthesis of key oxylipins involved in tiller growth and development.

CONCLUSIONS

Collectively, we have enriched the genomic data available in sugarcane and provided candidate genes for manipulating tiller formation and development, towards productivity enhancement in sugarcane.

Identification of the specific long-noncoding RNAs involved in night-break mediated flowering retardation in Chenopodium quinoa

Background

Night-break (NB) has been proven to repress flowering of short-day plants (SDPs). Long-noncoding RNAs (lncRNAs) play key roles in plant flowering. However, investigation of the relationship between lncRNAs and NB responses is still limited, especially in Chenopodium quinoa, an important short-day coarse cereal.

Results

In this study, we performed strand-specific RNA-seq of leaf samples collected from quinoa seedlings treated by SD and NB. A total of 4914 high-confidence lncRNAs were identified, out of which 91 lncRNAs showed specific responses to SD and NB. Based on the expression profiles, we identified 17 positive- and 7 negative-flowering lncRNAs. Co-expression network analysis indicated that 1653 mRNAs were the common targets of both types of flowering lncRNAs. By mapping these targets to the known flowering pathways in model plants, we found some pivotal flowering homologs, including 2 florigen encoding genes (FT (FLOWERING LOCUS T) and TSF (TWIN SISTER of FT) homologs), 3 circadian clock related genes (EARLY FLOWERING 3 (ELF3), LATE ELONGATED HYPOCOTYL (LHY) and ELONGATED HYPOCOTYL 5 (HY5) homologs), 2 photoreceptor genes (PHYTOCHROME A (PHYA) and CRYPTOCHROME1 (CRY1) homologs), 1 B-BOX type CONSTANS (CO) homolog and 1 RELATED TO ABI3/VP1 (RAV1) homolog, were specifically affected by NB and competed by the positive and negative-flowering lncRNAs. We speculated that these potential flowering lncRNAs may mediate quinoa NB responses by modifying the expression of the floral homologous genes.

Conclusions

Together, the findings in this study will deepen our understanding of the roles of lncRNAs in NB responses, and provide valuable information for functional characterization in future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07605-2.

Transcriptomic Profiling of Fe-Responsive lncRNAs and Their Regulatory Mechanism in Rice

Iron (Fe) deficiency directly affects crop growth and development, ultimately resulting in reduced crop yield and quality. Recently, long non-coding RNAs (lncRNAs) have been demonstrated to play critical regulatory roles in a multitude of pathways across numerous species. However, systematic screening of lncRNAs responding to Fe deficiency and their regulatory mechanism in plants has not been reported. In this work, 171 differently expressed lncRNAs (DE-lncRNAs) were identified based on analysis of strand-specific RNA-seq data from rice shoots and roots under Fe-deficient conditions. We also found several lncRNAs, which could generate miRNAs or act as endogenous target mimics to regulate expression of Fe-related genes. Analysis of interaction networks and gene ontology enrichment revealed that a number of DE-lncRNAs were associated with iron transport and photosynthesis, indicating a possible role of lncRNAs in regulation of Fe homeostasis. Moreover, we identified 76 potential lncRNA targets of OsbHLH156, a key regulator for transcriptional response to Fe deficiency. This study provides insight into the potential functions and regulatory mechanism of Fe-responsive lncRNAs and would be an initial and reference for any further studies regarding lncRNAs involved in Fe deficiency in plants.

Genome-wide identification and analysis of long non-coding RNAs involved in fatty acid biosynthesis in young soybean pods

Long non-coding RNAs (lncRNAs) are non-coding RNAs of more than 200 nucleotides. To date, the roles of lncRNAs in soybean fatty acid synthesis have not been fully studied. Here, the low-linolenic acid mutant 'MT72' and the wild-type control 'JN18' were used as materials. The lncRNAs in young pods at 30 and 40 days (d) after flowering were systematically identified and analyzed using transcriptome sequencing technology combined with bioinformatics tools. A total of 39,324 lncRNAs and 561 differentially expressed lncRNAs were identified. A lncRNAs-miRNAs-protein-coding genes (mRNAs) network was constructed, and 46 lncRNAs, 46 miRNAs and 137 mRNAs were found to be correlated. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of 12 targeted mRNAs in the competing endogenous RNA network showed that these lncRNAs may be involved in the biological processes of fatty acid transport, lipid synthesis and cell division. Finally, the expression levels of differentially expressed lncRNAs, miRNAs and mRNAs were verified using qRT-PCR. The expression patterns of most genes were consistent with the sequencing results. In conclusion, new information was provided for the study of fatty acid synthesis by lncRNAs in young soybean pods.

Genome-wide identification and characterization of long non-coding RNAs involved in flag leaf senescence of rice

KEY MESSAGE

This study showed the systematic identification of long non-coding RNAs (lncRNAs) involving in flag leaf senescence of rice, providing the possible lncRNA-mRNA regulatory relationships and lncRNA-miRNA-mRNA ceRNA networks during leaf senescence. LncRNAs have been reported to play crucial roles in diverse biological processes. However, no systematic identification of lncRNAs associated with leaf senescence in plants has been studied. In this study, a genome-wide high throughput sequencing analysis was performed using rice flag leaves developing from normal to senescence. A total of 3953 lncRNAs and 38757 mRNAs were identified, of which 343 lncRNAs and 9412 mRNAs were differentially expressed. Through weighted gene co-expression network analysis (WGCNA), 22 continuously down-expressed lncRNAs targeting 812 co-expressed mRNAs and 48 continuously up-expressed lncRNAs targeting 1209 co-expressed mRNAs were considered to be significantly associated with flag leaf senescence. Gene Ontology results suggested that the senescence-associated lncRNAs targeted mRNAs involving in many biological processes, including transcription, hormone response, oxidation-reduction process and substance metabolism. Additionally, 43 senescence-associated lncRNAs were predicted to target 111 co-expressed transcription factors. Interestingly, 8 down-expressed lncRNAs and 29 up-expressed lncRNAs were found to separately target 12 and 20 well-studied senescence-associated genes (SAGs). Furthermore, analysis on the competing endogenous RNA (CeRNA) network revealed that 6 down-expressed lncRNAs possibly regulated 51 co-expressed mRNAs through 15 miRNAs, and 14 up-expressed lncRNAs possibly regulated 117 co-expressed mRNAs through 21 miRNAs. Importantly, by expression validation, a conserved miR164-NAC regulatory pathway was found to be possibly involved in leaf senescence, where lncRNA MSTRG.62092.1 may serve as a ceRNA binding with miR164a and miR164e to regulate three transcription factors. And two key lncRNAs MSTRG.31014.21 and MSTRG.31014.36 also could regulate the abscisic-acid biosynthetic gene BGIOSGA025169 (OsNCED4) and BGIOSGA016313 (NAC family) through osa-miR5809. The possible regulation networks of lncRNAs involving in leaf senescence were discussed, and several candidate lncRNAs were recommended for prior transgenic analysis. These findings will extend the understanding on the regulatory roles of lncRNAs in leaf senescence, and lay a foundation for functional research on candidate lncRNAs.

Genome-wide analysis of long noncoding RNAs, 24-nt siRNAs, DNA methylation and H3K27me3 marks in Brassica rapa

Long noncoding RNAs (lncRNAs) are RNA fragments that generally do not code for a protein but are involved in epigenetic gene regulation. In this study, lncRNAs of Brassica rapa were classified into long intergenic noncoding RNAs, natural antisense RNAs, and intronic noncoding RNAs and their expression analyzed in relation to genome-wide 24-nt small interfering RNAs (siRNAs), DNA methylation, and histone H3 lysine 27 trimethylation marks (H3K27me3). More than 65% of the lncRNAs analyzed consisted of one exon, and more than 55% overlapped with inverted repeat regions (IRRs). Overlap of lncRNAs with IRRs or genomic regions encoding for 24-nt siRNAs resulted in increased DNA methylation levels when both were present. LncRNA did not overlap greatly with H3K27me3 marks, but the expression level of intronic noncoding RNAs that did coincide with H3K27me3 marks was higher than without H3K27me3 marks. The Brassica genus comprises important vegetables and oil seed crops grown across the world. B. rapa is a diploid (AA genome) thought to be one of the ancestral species of both B. juncea (AABB genome) and B. napus (AACC) through genome merging (allotetrapolyploidization). Complex genome restructuring and epigenetic alterations are thought to be involved in these allotetrapolyploidization events. Comparison of lncRNAs between B. rapa and B. nigra, B. oleracea, B. juncea, and B. napus showed the highest conservation with B. oleracea. This study presents a comprehensive analysis of the epigenome structure of B. rapa at multi-epigenetic levels (siRNAs, DNA methylation, H3K27me3, and lncRNAs) and identified a suite of candidate lncRNAs that may be epigenetically regulated in the Brassica genus.

Genome-Wide Identification and Characterization of Long Noncoding RNAs Involved in Chinese Wheat Mosaic Virus Infection of Nicotiana benthamiana

Simple Summary

Recent studies have shown that a large number of long noncoding RNAs (lncRNAs) can regulate various biological processes in animals and plants. However, the roles of long non-coding RNAs (lncRNAs) in the interaction between plants and viruses is unclear, particularly for the Chinese wheat mosaic virus (CWMV) interaction. In this study, we used a deep RNA sequencing strategy to profile lncRNAs involved in the response to CWMV infection in Nicotiana benthamiana and analyzed differentially expressed lncRNAs that responded to CWMV infection, using a bioinformatics method. We identified 1175 new lncRNAs in N. benthamiana infected with CWMV, with 65 lncRNAs showing differential expression. These lncRNAs were mainly enriched in plant hormone signal transduction and other pathways according to GO and KEGG pathway enrichment analyses. In addition, differential expression of XLOC_006393 after CWMV infection may be the precursor of NbmiR168c, which can respond to CWMV infection by modulating the expression of its target gene NbAGO1. We believe that our study makes a significant contribution to the literature because these results provide a valuable resource for studying lncRNAs involved in CWMV infection and improving the understanding of the molecular mechanism of CWMV infection.

Abstract

Recent studies have shown that a large number of long noncoding RNAs (lncRNAs) can regulate various biological processes in animals and plants. Although lncRNAs have been identified in many plants, they have not been reported in the model plant Nicotiana benthamiana. Particularly, the role of lncRNAs in plant virus infection remains unknown. In this study, we identified lncRNAs in N. benthamiana response to Chinese wheat mosaic virus (CWMV) infection by RNA sequencing. A total of 1175 lncRNAs, including 65 differentially expressed lncRNAs, were identified during CWMV infection. We then analyzed the functions of some of these differentially expressed lncRNAs. Interestingly, one differentially expressed lncRNA, XLOC_006393, was found to participate in CWMV infection as a precursor to microRNAs in N. benthamiana. These results suggest that lncRNAs play an important role in the regulatory network of N. benthamiana in response to CWMV infection.

Identification of long noncoding RNAs involved in resistance to downy mildew in Chinese cabbage

Brassica downy mildew, a severe disease caused by Hyaloperonospora brassicae, can cause enormous economic losses in Chinese cabbage (Brassica rapa L. ssp. pekinensis) production. Although some research has been reported recently concerning the underlying resistance to this disease, no studies have identified or characterized long noncoding RNAs involved in this defense response. In this study, using high-throughput RNA sequencing, we analyzed the disease-responding mRNAs and long noncoding RNAs in two resistant lines (T12-19 and 12-85) and one susceptible line (91-112). Clustering and Gene Ontology analysis of differentially expressed genes (DEGs) showed that more DEGs were involved in the defense response in the two resistant lines than in the susceptible line. Different expression patterns and proposed functions of differentially expressed long noncoding RNAs among T12-19, 12-85, and 91-112 indicated that each has a distinct disease response mechanism. There were significantly more cis- and trans-functional long noncoding RNAs in the resistant lines than in the susceptible line, and the genes regulated by these RNAs mostly participated in the disease defense response. Furthermore, we identified a candidate resistance-related long noncoding RNA, MSTRG.19915, which is a long noncoding natural antisense transcript of a MAPK gene, BrMAPK15. Via an agroinfiltration-mediated transient overexpression system and virus-induced gene silencing technology, BrMAPK15 was indicated to have a greater ability to defend against pathogens. MSTRG.19915-silenced seedlings showed enhanced resistance to downy mildew, probably because of the upregulated expression of BrMAPK15. This research identified and characterized long noncoding RNAs involved in resistance to downy mildew, laying a foundation for future in-depth studies of disease resistance mechanisms in Chinese cabbage.

Comprehensive transcriptome profiling reveals abundant long non-coding RNAs associated with development of the rice false smut fungus, Ustilaginoidea virens

Long non-coding RNAs (lncRNAs) play an important role in biological processes but regulation and function of lncRNAs remain largely unelucidated, especially in fungi. Ustilaginoidea virens is an economically important fungus causing a devastating disease of rice. By combining microscopic and RNA-seq analyses, we comprehensively characterized lncRNAs of this fungus in infection and developmental processes and defined four serial typical stages. RNA-seq analyses revealed 1724 lncRNAs in U. virens, including 1084 long intergenic non-coding RNAs (lincRNAs), 51 intronic RNAs (incRNAs), 566 natural antisense transcripts (lncNATs) and 23 sense transcripts. Gene Ontology enrichment of differentially expressed lincRNAs and lncNATs demonstrated that these were mainly involved in transport-related regulation. Functional studies of transport-related lncRNAs revealed that UvlncNAT-MFS, a cytoplasm localized lncNAT of a putative MFS transporter gene, UvMFS, could form an RNA duplex with UvMFS and was required for regulation of growth, conidiation and various stress responses. Our results were the first to elucidate the lncRNA profiles during infection and development of this important phytopathogen U. virens. The functional discovery of the novel lncRNA, UvlncNAT-MFS, revealed the potential of lncRNAs in regulation of life processes in fungi.