Genome-Wide Identification and Characterization of Long Non-Coding RNA in Wheat Roots in Response to Ca2+ Channel Blocker

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 emerging role of long non-coding RNAs (lncRNAs) in plant biology: Functions, mechanisms of action, and future directions

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

Leveraging a meta-learning approach to advance the accuracy of Nav blocking peptides prediction

The voltage-gated sodium (Nav) channel is a crucial molecular component responsible for initiating and propagating action potentials. While the α subunit, forming the channel pore, plays a central role in this function, the complete physiological function of Nav channels relies on crucial interactions between the α subunit and auxiliary proteins, known as protein-protein interactions (PPI). Nav blocking peptides (NaBPs) have been recognized as a promising and alternative therapeutic agent for pain and itch. Although traditional experimental methods can precisely determine the effect and activity of NaBPs, they remain time-consuming and costly. Hence, machine learning (ML)-based methods that are capable of accurately contributing in silico prediction of NaBPs are highly desirable. In this study, we develop an innovative meta-learning-based NaBP prediction method (MetaNaBP). MetaNaBP generates new feature representations by employing a wide range of sequence-based feature descriptors that cover multiple perspectives, in combination with powerful ML algorithms. Then, these feature representations were optimized to identify informative features using a two-step feature selection method. Finally, the selected informative features were applied to develop the final meta-predictor. To the best of our knowledge, MetaNaBP is the first meta-predictor for NaBP prediction. Experimental results demonstrated that MetaNaBP achieved an accuracy of 0.948 and a Matthews correlation coefficient of 0.898 over the independent test dataset, which were 5.79% and 11.76% higher than the existing method. In addition, the discriminative power of our feature representations surpassed that of conventional feature descriptors over both the training and independent test datasets. We anticipate that MetaNaBP will be exploited for the large-scale prediction and analysis of NaBPs to narrow down the potential NaBPs.

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.

Epigenetic mechanisms of lncRNA in response to thermal stress during embryogenesis of allotetraploid Cyprinus carpio

Embryogenesis and epigenetic mechanisms of lncRNA may play an important role in the formation of temperature tolerance in allotetraploid Cyprinus carpio. To investigate the response of lncRNA to thermal stress during embryogenesis of C. carpio, transcriptome sequencing was performed on 81 embryo or larva samples from different early development stages and temperatures. We identified 45,097 lncRNAs and analyzed transcriptome variation during embryogenesis. Stage-specific and temperature-specific DE lncRNAs and DEGs were screened. GO and KEGG analysis identified numerous pathways involved in thermal stress. Temperature-specific regulation of cis-/trans-/antisense lncRNAs was analyzed. Interaction network analysis identified 6 hub lncRNAs and many hub genes, such as cdk1 and hsf1. Decreased expression of many essential genes regulated by lncRNAs may lead to the death of embryos at 33 °C. Our findings provide new insights into the regulation of lncRNA in thermal stress response during embryogenesis and contribute to the understanding of environmental adaptation of allotetraploid species.

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.

Kinase MxMPK4-1 and calmodulin-binding protein MxIQM3 enhance apple root acidification during Fe deficiency

Iron (Fe) deficiency is a long-standing issue in plant mineral nutrition. Ca2+ signals and the mitogen-activated protein kinase (MAPK) cascade are frequently activated in parallel to perceive external cues, but their interplay under Fe deficiency stress remains largely unclear. Here, the kinase MxMPK4-1, which is induced during the response to Fe deficiency stress in apple rootstock Malus xiaojinensis, cooperates with IQ-motif containing protein3 (MxIQM3). MxIQM3 gene expression, protein abundance, and phosphorylation level increased under Fe deficiency stress. The overexpression of MxIQM3 in apple calli and rootstocks mitigated the Fe deficiency phenotype and improved stress tolerance, whereas RNA interference or silencing of MxIQM3 in apple calli and rootstocks, respectively, worsened the phenotype and reduced tolerance to Fe deficiency. MxMPK4-1 interacted with MxIQM3 and subsequently phosphorylated MxIQM3 at Ser393, and co-expression of MxMPK4-1 and MxIQM3 in apple calli and rootstocks enhanced Fe deficiency responses. Furthermore, MxIQM3 interacted with the central-loop region of the plasma membrane (PM) H+-ATPase MxHA2. Phospho-mimicking mutation of MxIQM3 at Ser393 inhibited binding to MxHA2, but phospho-abolishing mutation promoted interaction with both the central-loop and C terminus of MxHA2, demonstrating phosphorylation of MxIQM3 caused dissociation from MxHA2 and therefore increased H+ secretion. Moreover, Ca2+/MxCAM7 (Calmodulin7) regulated the MxMPK4-1-MxIQM3 module in response to Fe deficiency stress. Overall, our results demonstrate that MxMPK4-1-MxIQM3 forms a functional complex and positively regulates PM H+-ATPase activity in Fe deficiency responses, revealing a versatile mechanism of Ca2+/MxCAM7 signaling and MAPK cascade under Fe deficiency stress.

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.

The diagnostic and prognostic value of pseudogene SIGLEC17P in lung adenocarcinoma and a preliminary functional study

Among malignant tumors, lung adenocarcinoma (LUAD) is the leading cause of death worldwide. This study explored the diagnostic, prognostic value, and preliminary functional verification of sialic acid binding Ig like lectin 17, pseudogene (SIGLEC17P) in LUAD. Prognostic lncRNAs for LUAD were identified by The Cancer Genome Atlas and quantitative real-time PCR (qRT-PCR) was used to detect the expression of SIGLEC17P in LUAD and paracarcinoma tissues. Subsequently, lentiviral vectors were used to overexpress SIGLEC17P in A549 and H1299 cells. The effects of SIGLEC17P overexpression on the proliferation, migration, and invasiveness of LUAD cells (A549 and H1299) were evaluated by Cell Counting Kit-8, wound healing, and transwell migration assays, respectively. Bioinformatics analyses were performed to reveal the potential pathways in which SIGLEC17P is involved in LUAD. qRT-PCR results revealed low SIGLEC17P expression in LUAD tissues and a significant association with the N stage, T stage, and tumor node metastasis stage. Furthermore, the receiver operating characteristic curve demonstrated a reliable diagnostic value. The proliferation, migration, and invasion of LUAD cells were inhibited by overexpression of SIGLEC17P. Bioinformatics analyses suggested that SIGLEC17P might exert antioncogenic effects in LUAD through the mir-20-3p/ADH1B or mir-4476-5p/DPYSL axis. In summary, our results revealed that SIGLEC17P acts as a prognostic biomarker, independent prognostic factor, and potential therapeutic target for patients with LUAD.

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.

Functional analysis of long non-coding RNAs involved in alkaline stress responses in wheat

Saline-alkali soil is a major environmental problem affecting crop productivity. One of the most effective approaches to combat it is to breed stress-tolerant plants through genetic engineering. Shanrong No. 4 (SR4) is an alkaline-tolerant cultivar of bread wheat (Triticum aestivum) derived from asymmetric somatic hybridization between the common wheat cultivar Jinan 177 (JN177) and tall wheatgrass. In this study, we aimed to explore the structure and function of alkalinity stress-responsive long non-coding RNAs (lncRNAs) in wheat. Sequencing was employed to identify the lncRNAs associated with stress tolerance and their corresponding targets. Approximately 19 000 novel lncRNA sequences were detected in SR4 and JN177. Upon exposure to alkaline stress, SR4 differentially expressed 5691 lncRNAs, whilst JN177 differentially expressed 5932. We selected five of them (L0760, L6247, L0208, L2098, and L3065) and generated seedlings of transiently knocked down strains using the virus-induced gene-silencing method. Knockdown of L0760 and L2098 caused the plants to exhibit sensitivity to alkaline stress, whereas knockdown of L6247, L0208, and L3065 increased the ability of plants to tolerate alkaline stress. We constructed lncRNA-miRNA-target-mRNA networks and alkali-response-related lncRNA-target-mRNA association networks to analyse the functions of lncRNAs. Collectively, our results demonstrate that lncRNAs may perform different roles under alkaline stress conditions.

Integrated Analysis of miRNAs and Gene Expression Profiles Reveals Potential Biomarkers for Osteoarthritis

Purpose: Currently, the early diagnosis and treatment of osteoarthritis (OA) remain a challenge. In the present study, we attempted to explore potential biomarkers for the diagnosis and treatment of OA.

Methods: The differentially expressed genes (DEGs) were identified based on three mRNA datasets of synovial tissues for OA patients and normal controls downloaded from the Gene Expression Omnibus (GEO) database. Furthermore, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were used for evaluating gene function related categories. Then, miRNA sequencing was performed for differentially expressed miRNAs’ identification. Finally, weighted gene co-expression network analysis (WGCNA) was performed for genes detected by the three mRNA datasets and a competing endogenous RNA (ceRNA) network with DEGs and differentially expressed microRNAs (miRNAs) was constructed for central genes identification. In addition, the relationship between central gene expression and immune infiltration was analyzed, and the candidate agents for OA were predicted based on the Connectivity Map database. Quantitative RT-PCR (qRT-PCR), Western blotting analysis, and immunofluorescent staining were performed to validate the expression levels of differentially expressed miRNAs and differentially expressed target genes in normal and OA tissues and chondrocytes. MiRNA–mRNA network was also validated in chondrocytes in vitro.

Results: A total of 259 DEGs and 26 differentially expressed miRNAs were identified, among which 94 miRNA–mRNA interactions were predicted. The brown module in WGCNA was most closely correlated with the clinical traits of OA. After overlapping the brown module genes with miRNA–mRNA pairs, 27 miRNA–mRNA pairs were obtained. A ceRNA network was constructed with 5505 lncRNA–miRNA–mRNA interactions. B-cell translocation gene 2(BTG2), Abelson-related gene (ABL2), and vascular endothelial growth factor A (VEGFA) were identified to be the central genes with good predictive performance, which were significantly correlated with immune cell infiltration in OA, reflected by declined activated dendritic cells (aDCs), and elevated contents of B cells, macrophages, neutrophils, and T helper cells. Anisomycin, MG-132, thapsigargin, and lycorine were predicted to be the potential candidate agents for OA intervention. In vitro, the expression levels of differentially expressed miRNAs and biomarkers identified in the present study were consistent with the results obtained in normal or OA knee cartilage tissues and chondrocytes. Furthermore, BTG2 was identified to be negatively regulated by miR-125a-5p.

Conclusion: BTG2, ABL2, and VEGFA can be regarded as potential predictive and treatment biomarkers for OA, which might guide the clinical therapy of OA.

Detection of Root Physiological Parameters and Potassium and Calcium Currents in the Rhizoplane of the Apple Rootstock Superior Line 12-2 With Improved Apple Replant Disease Resistance

The cultivation of resistant rootstocks is one of the more effective ways to mitigate apple replant disease (ARD). We performed an ion current test, a pot experiment, and a pathogen infection test on the apple rootstocks 12-2 (self-named), T337, and M26. The ion current test showed that exposure to ARD soil extract for 30 min had a significant effect on K+ ion currents at the meristem, elongation, and mature zones of the M26 rhizoplane and on Ca2+ currents in the meristem and elongation zones. ARD also had a significant effect on Ca2+ currents in the meristem, elongation, and mature zones of the T337 rhizoplane. Exposure to ARD soil extract for 5 min had a significant effect on K+ currents in the meristem, elongation, and mature zones of 12-2 and on the Ca2+ currents in the elongation and mature zones. Compared to a 5-min exposure, a 30-min exposure to ARD extract had a less pronounced effect on K+ and Ca2+ currents in the 12-2 rhizoplane. The pot experiment showed that ARD soil had no significant effect on any root architectural or physiological parameters of 12-2. By contrast, ARD soil significantly reduced some root growth indices and the dry and fresh weights of T337 and M26 compared with controls on sterilized soil. ARD also had a significant effect on root metabolic activity, root antioxidant enzyme activity (except superoxide dismutase for T337), and malondialdehyde content of T337 and M26. Pathogen infection tests showed that Fusarium proliferatum MR5 significantly affected the root structure and reduced the root metabolic activity of T337 and M26. It also reduced their root antioxidant enzyme activities (except catalase for T337) and significantly increased the root malondialdehyde content, reactive oxygen levels, and proline and soluble sugar contents. By contrast, MR5 had no such effects on 12-2. Based on these results, 12-2 has the potential to serve as an important ARD-resistant rootstock.

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

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

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

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

Genome-wide identification and characterization of long non-coding RNAs conferring resistance to Colletotrichum gloeosporioides in walnut (Juglans regia)

BACKGROUND

Walnut anthracnose caused by Colletotrichum gloeosporioides (Penz.) Penz. and Sacc. is an important walnut production problem in China. Although the long non-coding RNAs (lncRNAs) are important for plant disease resistance, the molecular mechanisms underlying resistance to C. gloeosporioides in walnut remain poorly understood.

RESULTS

The anthracnose-resistant F26 fruits from the B26 clone and the anthracnose-susceptible F423 fruits from the 4-23 clone of walnut were used as the test materials. Specifically, we performed a comparative transcriptome analysis of F26 and F423 fruit bracts to identify differentially expressed LncRNAs (DELs) at five time-points (tissues at 0 hpi, pathological tissues at 24 hpi, 48 hpi, 72 hpi, and distal uninoculated tissues at 120 hpi). Compared with F423, a total of 14,525 DELs were identified, including 10,645 upregulated lncRNAs and 3846 downregulated lncRNAs in F26. The number of upregulated lncRNAs in F26 compared to in F423 was significantly higher at the early stages of C. gloeosporioides infection. A total of 5 modules related to disease resistance were screened by WGCNA and the target genes of lncRNAs were obtained. Bioinformatic analysis showed that the target genes of upregulated lncRNAs were enriched in immune-related processes during the infection of C. gloeosporioides, such as activation of innate immune response, defense response to bacterium, incompatible interaction and immune system process, and enriched in plant hormone signal transduction, phenylpropanoid biosynthesis and other pathways. And 124 known target genes for 96 hub lncRNAs were predicted, including 10 known resistance genes. The expression of 5 lncRNAs and 5 target genes was confirmed by qPCR, which was consistent with the RNA-seq data.

CONCLUSIONS

The results of this study provide the basis for future functional characterizations of lncRNAs regarding the C. gloeosporioides resistance of walnut fruit bracts.

Genome-wide discovery and characterization of long noncoding RNAs in African oil palm (Elaeis guineensis Jacq.)

Long noncoding RNAs (lncRNAs) are an important class of genes and play important roles in a range of biological processes. However, few reports have described the identification of lncRNAs in oil palm. In this study, we applied strand specific RNA-seq with rRNA removal to identify 1,363 lncRNAs from the equally mixed tissues of oil palm spear leaf and six different developmental stages of mesocarp (8–24 weeks). Based on strand specific RNA-seq data and 18 released oil palm transcriptomes, we systematically characterized the expression patterns of lncRNA loci and their target genes. A total of 875 uniq target genes for natural antisense lncRNAs (NAT-lncRNA, 712), long intergenic noncoding RNAs (lincRNAs, 92), intronic-lncRNAs (33), and sense-lncRNAs (52) were predicted. A majority of lncRNA loci (77.8%–89.6%) had low expression in 18 transcriptomes, while only 89 lncRNA loci had medium to high expression in at least one transcriptome. Coexpression analysis between lncRNAs and their target genes indicated that 6% of lncRNAs had expression patterns positively correlated with those of target genes. Based on single nucleotide polymorphism (SNP) markers derived from our previous research, 6,882 SNPs were detected for lncRNAs and 28 SNPs belonging to 21 lncRNAs were associated with the variation of fatty acid contents. Moreover, seven lncRNAs showed expression patterns positively correlated expression pattern with those of genes in de novo fatty acid synthesis pathways. Our study identified a collection of lncRNAs for oil palm and provided clues for further research into lncRNAs that may regulate mesocarp development and lipid metabolism.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Genome-wide identification and characterization of novel non-coding RNA-derived SSRs in wheat

Expression of eukaryotic genes is largely regulated by non-coding RNAs (ncRNA). Sequence variations in the regulatory RNAs may have critical biological consequences including transcriptional and post-transcriptional gene regulation. ncRNA-derived markers thus can be proved useful in molecular breeding, QTL mapping and association studies for trait dissection. In present study, we identified a total of 661 SSRs dwelling in pre-miRNA (15), small nuclear RNA (25) and lncRNA (621). Of these, 46 were validated and 100% amplification success was observed in selected wheat genotypes. A set of 36 ncRNA-SSRs markers was utilized for genetic variability assessment in forty-eight Indian wheat genotypes (which includes bread wheat, durum wheat and relatives). Number of alleles ranged from 1 to 4 with an average of two alleles per SSR locus. Mean PIC, observed heterozygosity and Shannon information index were found to be 0.258, 0.37 and 0.476 which suggests ncRNA-SSRs show higher polymorphism compared to genic SSRs but lower polymorphism compared to genomic SSRs. Thirty-six ncRNA-SSRs showed transferability ranging from 42.1% to 100%. Average genetic dissimilarity among wheat genotypes was found to be 0.29 based on Jaccard's dissimilarity. This is the first report of ncRNA-SSRs in wheat which will be useful for molecular breeding and genetic improvement of wheat.

Identification of lncRNAs involved in wheat tillering development in two pairs of near-isogenic lines

Emerging evidence demonstrates that lncRNAs participate in various developmental processes in plants via post-transcription regulation. However, few lncRNAs have been identified as regulators of tiller development in wheat (Triticum aestivum L.). In this study, high-throughput ribosomal depleted RNA sequencing was performed on the tillering nodes of two pairs of near-isogenic lines that differed only in the tillering trait. We identified 5399 lncRNA transcripts using bioinformational analyses. KEGG pathway analysis revealed 74 common differentially expressed lncRNAs substantially enriched in photosynthesis-related, phenylpropanoid biosynthesis, phosphatidylinositol signaling, brassinosteroid biosynthesis, zeatin biosynthesis, and carotenoid biosynthesis pathways. Detailed functional annotations of target genes were used to identify 27 tillering-associated lncRNAs. Among these, 10 were in photosynthesis-related pathways; 15 were in secondary metabolite pathways; and 8 were in plant hormone pathways, with 6 enriched in two kinds of pathways. These findings contribute to identifying tillering-associated lncRNAs in wheat and enable further investigation into the functions and roles of key candidate lncRNAs, and more experimental evidence was also needed if breeders wanted to utilize these candidate lncRNAs in wheat crop yield improvement in the future.

Influence of Cd exposure on H+ and Cd2+ fluxes in the leaf, stem and root of a novel aquatic hyperaccumulator - Microsorum pteropus

Microsorum pteropus has been proven to be a potential novel aquatic Cd hyperaccumulator. In this study, Non-invasive Micro-test Technology (NMT) was used to observe the ion fluxes of different M. pteropus tissues under Cd exposure. M. pteropus can hyperaccumulate more than 1000 mg/kg Cd in roots and leaves and approximately 600 mg/kg Cd in stems after seven days of exposure to 500 μM Cd, showing that this plant have a great capacity for Cd enrichment and resistance. The NMT test found H⁺ fluxes increased in all tissues after Cd exposure, with the largest increases being observed in stems, followed by the leaves and roots. Cd²⁺ fluxes showed different accumulation levels in different tissues, with low-level Cd exposure leading to influxes into roots and leaves, and high-level Cd exposure resulting in effluxes from roots. No significant influxes or effluxes were observed in leaves under high-level Cd exposure, or in stems under low- and high-levels of Cd exposure. However, transient high-level Cd exposure showed long-term Cd²⁺ influxes into roots and short-term Cd²⁺ effluxes out of stems and leaves. The roots of M. pteropus had greater regulation mechanisms for Cd enrichment and resistance, with influxes occurring following low-level exposure and effluxes occurring from high-level exposure. When exposed to Cd, M. pteropus stems showed less transportation and absorption. Low-level Cd exposure resulted in individual leaves directly absorbing Cd from hydroponic solutions. Different Cd enrichment and resistance mechanisms were exhibited by different M. pteropus tissues.

Identification and characterization of long non-coding RNAs regulating resistant starch biosynthesis in bread wheat (Triticum aestivum L.)

Resistant starch (RS) also known as healthy starch has shown several health benefits. Enhancing the RS through genetic modification approaches has huge commercial importance. Regulatory RNA like long non-coding RNA (lncRNA) plays an important role in gene regulation. In this study, we mined 63 transcriptome datasets of wheat belonging to 35 genotypes representing two seed developmental stages. Contrasting expression of a subset of lncRNAs in RS mutant lines compared to parent wheat variety 'C 306' signifies their probable role in RS biosynthesis. Further, lncRNA- TCONS_00130663 showed strong positive correlation (r² = 1) with LYPL gene and strong negative correlation with SBEIIb (r² = -0.94). We found TCONS_00130663 as positive regulator of LYPL gene through interaction with miR1128. Based on relative expression, in silico interaction and DSC analysis we hypothesize the dual role of TCONS_00130663 in RS type 2 and type 5. The study provides a useful resource for functional mechanism of lncRNAs.

Pentoxifylline alleviated cardiac injury via modulating the cardiac expression of lncRNA-00654-miR-133a-SOX5 mRNA in the rat model of ischemia-reperfusion

Pentoxifylline (PTX) protects from many cardiovascular complications. It plays a critical role in stem cell proliferation and differentiation. Here, the effect of PTX administration on cardiac ischemia and dysfunction was explored. PTX in 3 doses (20, 30, and 40 mg/kg), was administered in vivo 5 min before a 45 min occlusion of the left anterior descending artery, followed by a 120 min reperfusion in male Wistar rats. The left ventricular end-diastolic pressure and dP/dt_max were assessed. Blood and cardiac tissue samples were collected for measuring the levels of cardiac enzymes and the expression of lncRNA-00654-miR-133a-SOX5. Samples of left ventricles were collected and processed for light microscopic, immunohistochemical staining for c-kit (a marker for cardiac progenitor cells) and transmission electron microscopic examination. PTX administration showed improvements in cardiac function tests, enzymes, and myocytes. Microscopic features showed minimal cardiac edema, hemorrhage, cellular inflammatory infiltration and fibrosis in addition to increased c-kit + cells in cardiac tissue samples. Notably, this treatment also produced a dose-dependent decrease in lncRNA-00654 with an increase in SOX5 mRNA and miRNA-133a-3p expressions. In conclusion, PTX has the potential to alleviate cardiac injury and increase the number of c-kit + cells following ischemia-reperfusion in the rat model via modulation of lncRNA-00654 and miR-133a-SOX5 mRNA expressions.

Integrated Analysis of Long Non-coding RNAs (lncRNAs) and mRNAs Reveals the Regulatory Role of lncRNAs Associated With Salt Resistance in Camellia sinensis

Tea plant (Camellia sinensis), an important economic crop, is seriously affected by various abiotic stresses, including salt stress, which severely diminishes its widespread planting. However, little is known about the roles of long non-coding RNAs (lncRNAs) in transcriptional regulation under salt stress. In this study, high-throughput sequencing of tea shoots under salt-stress and control conditions was performed. Through sequencing analysis, 16,452 unique lncRNAs were identified, including 172 differentially expressed lncRNAs (DE-lncRNAs). The results of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of their cis- and trans-target genes showed that these DE-lncRNAs play important roles in many pathways such as the galactinol synthase (GOLS), calcium signaling pathway, and interact with transcription factors (TFs) under salt stress. The data from the gene-specific antisense oligodeoxynucleotide-mediated reduction in the lncRNA MSTRG.139242.1 and its predicted interacting gene, TEA027212.1 (Ca²⁺-ATPase 13), in tea leaves revealed that MSTRG.139242.1 may function in the response of tea plants to high salinity. In addition, 12 lncRNAs were predicted to be target mimics of 17 known mature miRNAs, such as miR156, that are related to the salt-stress response in C. sinensis. Our results provide new insights into lncRNAs as ubiquitous regulators in response to salt stress in tea plants.

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

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

Whole-transcriptome RNA sequencing reveals the global molecular responses and ceRNA regulatory network of mRNAs, lncRNAs, miRNAs and circRNAs in response to copper toxicity in Ziyang Xiangcheng (Citrus junos Sieb. Ex Tanaka)

BACKGROUND

Copper (Cu) toxicity has become a potential threat for citrus production, but little is known about related mechanisms. This study aims to uncover the global landscape of mRNAs, long non-coding RNAs (lncRNAs), circular RNAs (circRNAs) and microRNAs (miRNAs) in response to Cu toxicity so as to construct a regulatory network of competing endogenous RNAs (ceRNAs) and to provide valuable knowledge pertinent to Cu response in citrus.

RESULTS

Tolerance of four commonly used rootstocks to Cu toxicity was evaluated, and 'Ziyang Xiangcheng' (Citrus junos) was found to be the most tolerant genotype. Then the roots and leaves sampled from 'Ziyang Xiangcheng' with or without Cu treatment were used for whole-transcriptome sequencing. In total, 5734 and 222 mRNAs, 164 and 5 lncRNAs, 45 and 17 circRNAs, and 147 and 130 miRNAs were identified to be differentially expressed (DE) in Cu-treated roots and leaves, respectively, in comparison with the control. Gene ontology enrichment analysis showed that most of the DEmRNAs and targets of DElncRNAs and DEmiRNAs were annotated to the categories of 'oxidation-reduction', 'phosphorylation', 'membrane', and 'ion binding'. The ceRNA network was then constructed with the predicted pairs of DEmRNAs-DEmiRNAs and DElncRNAs-DEmiRNAs, which further revealed regulatory roles of these DERNAs in Cu toxicity.

CONCLUSIONS

A large number of mRNAs, lncRNAs, circRNAs, and miRNAs in 'Ziyang Xiangcheng' were altered in response to Cu toxicity, which may play crucial roles in mitigation of Cu toxicity through the ceRNA regulatory network in this Cu-tolerant rootstock.

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

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

LncRNA33732-respiratory burst oxidase module associated with WRKY1 in tomato- Phytophthora infestans interactions

Our previous studies indicated that tomato WRKY1 transcription factor acts as a positive regulator during tomato resistance to Phytophthora infestans. However, the molecular mechanism of WRKY1-mediated resistance regulation remains unclear. Here, we used a comparative transcriptome analysis between wild-type and WRKY1-overexpressing tomato plants to identify differentially expressed genes (DEGs) and long non-coding RNAs (DELs), and we examined long non-coding RNA (lncRNA)-gene networks. The promoter sequences of the upregulated DEGs and DELs were analyzed. Among 1073 DEGs and 199 DELs, 1 kb 5'-upstream regions of 59 DEGs and 22 DELs contain the W-box, the target sequence of the WRKY1. The results of promoter-β-glucuronidase (GUS) fusion and yeast one-hybrid assay showed that lncRNA33732 was activated by WRKY1 through sequence-specific interactions with the W-box element in its promoter. The overexpression and silencing analysis of lncRNA33732 in tomato showed that lncRNA33732 acts as a positive regulator and enhanced tomato resistance to P. infestans by induction of the expression of respiratory burst oxidase (RBOH) and increase in the accumulation of H₂ O₂ . When the expression of RBOH gene was inhibited in tomato plants, H₂ O₂ accumulation decreased and resistance were impaired. These findings suggest that lncRNA33732 activated by WRKY1 induces RBOH expression to increase H₂ O₂ accumulation in early defense reaction of tomato to P. infestans attack. Our results provide insights into the WRKY1-lncRNA33732-RBOH module involved in the regulation of H₂ O₂ accumulation and resistance to P. infestans, as well as provide candidates to enhance broad-spectrum resistance to pathogens in tomato.

Genome-wide identification of oil biosynthesis-related long non-coding RNAs in allopolyploid Brassica napus

Background

Long noncoding RNAs (lncRNAs) are transcripts longer than 200 bp that do not encode proteins but nonetheless have been shown to play important roles in various biological processes in plants. Brassica napus is an important seed oil crop worldwide and the target of many genetic improvement activities. To understand better the function of lncRNAs in regulating plant metabolic activities, we carried out a genome-wide lncRNA identification of lncRNAs in Brassica napus with a focus on lncRNAs involved in lipid metabolism. Twenty ribosomal RNA depleted strand specific RNA-seq (ssRNA-seq) datasets were generatred using RNAs isolated from B. napus seeds at four developmental stages. For comparison we also included 30 publically available RNA-seq datasets generated from poly(A) enriched mRNAs isolated from from various Brassica napus tissues in our analysis.

Results

A total of 8905 lncRNA loci were identified, including 7100 long intergenic noncoding RNA (lincRNA) loci and 1805 loci generating long noncoding natural antisense transcript (lncNAT). Many lncRNAs were identified only in the ssRNA-seq and poly(A) RNA-seq dataset, suggesting that B. napus has a large lncRNA repertoire and it is necessary to use libraries prepared from different tissues and developmental stages as well as different library preparation approaches to capture the whole spectrum of lncRNAs. Analysis of coexpression networks revealed that among the regulatory modules are networks containing lncRNAs and protein-coding genes related to oil biosynthesis indicating a possible role of lncRNAs in the control of lipid metabolism. One such example is that several lncRNAs are potential regulators of BnaC08g11970D that encodes oleosin1, a protein found in oil bodies and involved in seed lipid accumulation. We also observed that the expression levels of B. napus lncRNAs is positively correlated with their conservation levels.

Conclusions

We demonstrated that the B. napus genome has a large number of lncRNA and that these lncRNAs are expressed broadly across many developmental times and in different tissue types. We also provide evidence indicating that specific lncRNAs appear to be important regulators of lipid biosynthesis forming regulatory networks with transcripts involved in lipid biosynthesis. We also provide evidence that these lncRNAs are conserved in other species of the Brassicaceae family.

Electronic supplementary material

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

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

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