Genome-wide identification and functional prediction of novel and fungi-responsive lincRNAs in Triticum aestivum

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.

DNABERT-based explainable lncRNA identification in plant genome assemblies

Long non-coding ribonucleic acids (lncRNAs) have been shown to play an important role in plant gene regulation, involving both epigenetic and transcript regulation. LncRNAs are transcripts longer than 200 nucleotides that are not translated into functional proteins but can be translated into small peptides. Machine learning models have predominantly used transcriptome data with manually defined features to detect lncRNAs, however, they often underrepresent the abundance of lncRNAs and can be biased in their detection. Here we present a study using Natural Language Processing (NLP) models to identify plant lncRNAs from genomic sequences rather than transcriptomic data. The NLP models were trained to predict lncRNAs for seven model and crop species (Zea mays, Arabidopsis thaliana, Brassica napus, Brassica oleracea, Brassica rapa, Glycine max and Oryza sativa) using publicly available genomic references. We demonstrated that lncRNAs can be accurately predicted from genomic sequences with the highest accuracy of 83.4% for Z. mays and the lowest accuracy of 57.9% for B. rapa, revealing that genome assembly quality might affect the accuracy of lncRNA identification. Furthermore, we demonstrated the potential of using NLP models for cross-species prediction with an average of 63.1% accuracy using target species not previously seen by the model. As more species are incorporated into the training datasets, we expect the accuracy to increase, becoming a more reliable tool for uncovering novel lncRNAs. Finally, we show that the models can be interpreted using explainable artificial intelligence to identify motifs important to lncRNA prediction and that these motifs frequently flanked the lncRNA sequence.

Roles of long non-coding RNAs in plant immunity

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

Identification and validation of coding and non-coding RNAs involved in high-temperature-mediated seed dormancy in common wheat

INTRODUCTION

Seed dormancy (SD) significantly decreases under high temperature (HT) environment during seed maturation, resulting in pre-harvest sprouting (PHS) damage under prolonged rainfall and wet weather during wheat harvest. However, the molecular mechanism underlying HT-mediated SD remains elusiveSeed dormancy (SD) significantly decreases under high temperature (HT) environment during seed maturation, resulting in pre-harvest sprouting (PHS) damage under prolonged rainfall and wet weather during wheat harvest. However, the molecular mechanism underlying HT-mediated SD remains elusive.

METHODS

Here, the wheat landrace 'Waitoubai' with strong SD and PHS resistance was treated with HT from 21 to 35 days post anthesis (DPA). Then, the seeds under HT and normal temperature (NT) environments were collected at 21 DPA, 28 DPA, and 35 DPA and subjected to whole-transcriptome sequencing.

RESULTS

The phenotypic data showed that the seed germination percentage significantly increased, whereas SD decreased after HT treatment compared with NT, consistent with the results of previous studies. In total, 5128 mRNAs, 136 microRNAs (miRNAs), 273 long non-coding RNAs (lncRNAs), and 21 circularRNAs were found to be responsive to HT, and some of them were further verified through qRT-PCR. In particular, the known gibberellin (GA) biosynthesis gene TaGA20ox1 (TraesCS3D02G393900) was proved to be involved in HT-mediated dormancy by using the EMS-mutagenized wheat cultivar Jimai 22. Similarly, a novel gene TaCDPK21 (TraesCS7A02G267000) involved in the calcium signaling pathway was validated to be associated with HT-mediated dormancy by using the EMS mutant. Moreover, TaCDPK21 overexpression in Arabidopsis and functional complementarity tests supported the negative role of TaCDPK21 in SD. We also constructed a co-expression regulatory network based on differentially expressed mRNAs, miRNAs, and lncRNAs and found that a novel miR27319 was located at a key node of this regulatory network. Subsequently, using Arabidopsis and rice lines overexpressing miR27319 precursor or lacking miR27319 expression, we validated the positive role of miR27319 in SD and further preliminarily dissected the molecular mechanism of miR27319 underlying SD regulation through phytohormone abscisic acid and GA biosynthesis, catabolism, and signaling pathways.

DISCUSSION

These findings not only broaden our understanding of the complex regulatory network of HT-mediated dormancy but also provide new gene resources for improving wheat PHS resistance to minimize PHS damage by using the molecular pyramiding approach.

Overexpression of lncRNA77580 Regulates Drought and Salinity Stress Responses in Soybean

Emerging evidence indicates that long non-coding RNAs (lncRNAs) play important roles in diverse biological processes. However, the biological functions of most plant lncRNAs are still unknown. We previously discovered a soybean abiotic-stress-related lncRNA, lncRNA77580, and cloned the entire full-length sequence. Here, in order to fully identify the function of lncRNA77580 in soybean stress response, we created transgenic soybean lines overexpressing lncRNA77580. Compared with the wild type, overexpression of lncRNA77580 enhances the drought tolerance of soybean. However, the transgenic plants exhibit increased sensitivity to high salinity at the seedling stage. We found that lncRNA77580 modulates the transcription of different gene sets during salt and drought stress response. Under water deficit at the reproductive stage, lncRNA77580 overexpression increases the seed yield by increasing the seed number per plant. These results provide insight into the role of lncRNA77580 in soybean stress response.

Puccinia striiformis f. sp. tritici effectors in wheat immune responses

The obligate biotrophic fungus Puccinia striiformis f. sp. tritici, which causes yellow (stripe) rust disease, is among the leading biological agents resulting in tremendous yield losses on global wheat productions per annum. The combatting strategies include, but are not limited to, fungicide applications and the development of resistant cultivars. However, evolutionary pressure drives rapid changes, especially in its "effectorome" repertoire, thus allowing pathogens to evade and breach resistance. The extracellular and intracellular effectors, predominantly secreted proteins, are tactical arsenals aiming for many defense processes of plants. Hence, the identity of the effectors and the molecular mechanisms of the interactions between the effectors and the plant immune system have long been targeted in research. The obligate biotrophic nature of P. striiformis f. sp. tritici and the challenging nature of its host, the wheat, impede research on this topic. Next-generation sequencing and novel prediction algorithms in bioinformatics, which are accompanied by in vitro and in vivo validation approaches, offer a speedy pace for the discovery of new effectors and investigations of their biological functions. Here, we briefly review recent findings exploring the roles of P. striiformis f. sp. tritici effectors together with their cellular/subcellular localizations, host responses, and interactors. The current status and the challenges will be discussed. We hope that the overall work will provide a broader view of where we stand and a reference point to compare and evaluate new findings.

Harnessing genetic resistance to rusts in wheat and integrated rust management methods to develop more durable resistant cultivars

Wheat is one of the most important staple foods on earth. Leaf rust, stem rust and stripe rust, caused by Puccini triticina, Puccinia f. sp. graminis and Puccinia f. sp. striiformis, respectively, continue to threaten wheat production worldwide. Utilization of resistant cultivars is the most effective and chemical-free strategy to control rust diseases. Convectional and molecular biology techniques identified more than 200 resistance genes and their associated markers from common wheat and wheat wild relatives, which can be used by breeders in resistance breeding programmes. However, there is continuous emergence of new races of rust pathogens with novel degrees of virulence, thus rendering wheat resistance genes ineffective. An integration of genomic selection, genome editing, molecular breeding and marker-assisted selection, and phenotypic evaluations is required in developing high quality wheat varieties with resistance to multiple pathogens. Although host genotype resistance and application of fungicides are the most generally utilized approaches for controlling wheat rusts, effective agronomic methods are required to reduce disease management costs and increase wheat production sustainability. This review gives a critical overview of the current knowledge of rust resistance, particularly race-specific and non-race specific resistance, the role of pathogenesis-related proteins, non-coding RNAs, and transcription factors in rust resistance, and the molecular basis of interactions between wheat and rust pathogens. It will also discuss the new advances on how integrated rust management methods can assist in developing more durable resistant cultivars in these pathosystems.

Identification of long non-coding RNAs involved in floral scent of Rosa hybrida

Long non-coding RNAs (lncRNAs) were found to play important roles in transcriptional, post-transcriptional, and epigenetic gene regulation in various biological processes. However, lncRNAs and their regulatory roles remain poorly studied in horticultural plants. Rose is economically important not only for their wide use as garden and cut flowers but also as important sources of natural fragrance for perfume and cosmetics industry, but presently little was known about the regulatory mechanism of the floral scent production. In this paper, a RNA-Seq analysis with strand-specific libraries, was performed to rose flowers in different flowering stages. The scented variety 'Tianmidemeng' (Rosa hybrida) was used as plant material. A total of 13,957 lncRNAs were identified by mining the RNA-Seq data, including 10,887 annotated lncRNAs and 3070 novel lncRNAs. Among them, 10,075 lncRNAs were predicted to possess a total of 29,622 target genes, including 54 synthase genes and 24 transcription factors related to floral scent synthesis. 425 lncRNAs were differentially expressed during the flowering process, among which 19 were differentially expressed among all the three flowering stages. Using weighted correlation network analysis (WGCNA), we correlate the differentially-expressed lncRNAs to synthesis of individual floral scent compounds. Furthermore, regulatory function of one of candidate lncRNAs for floral scent synthesis was verified using VIGS method in the rose. In this study, we were able to show that lncRNAs may play important roles in floral scent production in the rose. This study also improves our understanding of how plants regulate their secondary metabolism by lncRNAs.

Regulatory non-coding RNA: The core defense mechanism against plant pathogens

Plant pathogens damage crops and threaten global food security. Plants have evolved complex defense networks against pathogens, using crosstalk among various signaling pathways. Key regulators conferring plant immunity through signaling pathways include protein-coding genes and non-coding RNAs (ncRNAs). The discovery of ncRNAs in plant transcriptomes was first considered "transcriptional noise". Recent reviews have highlighted the importance of non-coding RNAs. However, understanding interactions among different types of noncoding RNAs requires additional research. This review attempts to consider how long-ncRNAs, small-ncRNAs and circular RNAs interact in response to pathogenic diseases within different plant species. Developments within genomics and bioinformatics could lead to the further discovery of plant ncRNAs, knowledge of their biological roles, as well as an understanding of their importance in exploiting the recent molecular-based technologies for crop protection.

Genome-wide identification of long intergenic non-coding RNAs for Ralstonia solanacearum resistance in tomato (Solanum lycopersicum)

There is growing evidences indicating that long intergenic ncRNAs (lincRNAs) play key roles in plant development and stress responses. To research tomato lincRNA functions during the interaction between tomato and Ralstonia solanacearum, RNA-seq data of tomato plants inoculated with R. solanacearum was analyzed. In this study, 315 possible lincRNAs were identified from RNA-seq data. Then 23 differentially expressed lincRNAs between tomato plants inoculated with R. solanacearum and control were identified and a total of 171 possible target genes for these differentially expressed lincRNAs were predicted. Through GO and KEGG analysis, we found that lincRNA might be involved in jasmonic acid and ethylene signaling pathways to respond to tomato bacterial wilt infection. Furthermore, lincRNA may also be involved in regulating the expression of AGO protein. Subsequently, analysis of expression patterns between differentially expressed lincRNAs and adjacent mRNAs by qRT-PCR revealed that part of lincRNAs and their possible target genes exhibited positive correlation. Taken together, these results suggest that lincRNAs play potential roles in tomato against R. solanacearum infection and will provide fundamental information about the lincRNA-based plant defense mechanisms.

WGCNA Identifies a Comprehensive and Dynamic Gene Co-Expression Network That Associates with Smut Resistance in Sugarcane

Sugarcane smut is a major fungal disease caused by Sporisorium scitamineum, which seriously reduces the yield and quality of sugarcane. In this study, 36 transcriptome data were collected from two sugarcane genotypes, YT93-159 (resistant) and ROC22 (susceptible) upon S. scitamineum infection. Data analysis revealed 20,273 (12,659 up-regulated and 7614 down-regulated) and 11,897 (7806 up-regulated and 4091 down-regulated) differentially expressed genes (DEGs) in YT93-159 and ROC22, respectively. A co-expression network was then constructed by weighted gene co-expression network analysis (WGCNA), which identified 5010 DEGs in 15 co-expressed gene modules. Four of the 15 modules, namely, Skyblue, Salmon, Darkorange, and Grey60, were significantly associated with smut resistance. The GO and KEGG enrichment analyses indicated that the DEGs involving in these four modules could be enriched in stress-related metabolic pathways, such as MAPK and hormone signal transduction, plant-pathogen interaction, amino acid metabolism, glutathione metabolism, and flavonoid, and phenylpropanoid biosynthesis. In total, 38 hub genes, including six from the Skyblue module, four from the Salmon module, 12 from the Darkorange module, and 16 from the Grey60 module, were screened as candidate hub genes by calculating gene connectivity in the corresponding network. Only 30 hub genes were amplifiable with RT-qPCR, of which 27 were up-regulated upon S. scitamineum infection. The results were consistent with the trend of gene expression in RNA-Seq, suggesting their positive roles in smut resistance. Interestingly, the expression levels of AOX, Cyb5, and LAC were higher in ROC22 than in YT93-159, indicating these three genes may act as negative regulators in response to S. scitamineum infection. This study revealed the transcriptome dynamics in sugarcane challenged by S. scitamineum infection and provided gene targets for smut resistance breeding in sugarcane.

Understanding of the various aspects of gene regulatory networks related to crop improvement

The hierarchical relationship between transcription factors, associated proteins, and their target genes is defined by a gene regulatory network (GRN). GRNs allow us to understand how the genotype and environment of a plant are incorporated to control the downstream physiological responses. During plant growth or environmental acclimatization, GRNs are diverse and can be differently regulated across tissue types and organs. An overview of recent advances in the development of GRN that speed up basic and applied plant research is given here. Furthermore, the overview of genome and transcriptome involving GRN research along with the exciting advancement and application are discussed. In addition, different approaches to GRN predictions were elucidated. In this review, we also describe the role of GRN in crop improvement, crop plant manipulation, stress responses, speed breeding and identifying genetic variations/locus. Finally, the challenges and prospects of GRN in plant biology are discussed.

Genome-Wide Identification of Powdery Mildew Responsive Long Non-Coding RNAs in Cucurbita pepo

Cucurbita pepo L. is an essential economic vegetable crop worldwide, and its production is severely affected by powdery mildew (PM). However, our understanding of the molecular mechanism of PM resistance in C. pepo is very limited. Long non-coding RNAs (lncRNAs) play an important role in regulating plant responses to biotic stress. Here, we systematically identified 2,363 reliably expressed lncRNAs from the leaves of PM-susceptible (PS) and PM-resistant (PR) C. pepo. The C. pepo lncRNAs are shorter in length and expressed at a lower level than the protein-coding transcripts. Among the 2,363 lncRNAs, a total of 113 and 146 PM-responsive lncRNAs were identified in PS and PR, respectively. Six PM-responsive lncRNAs were predicted as potential precursors of microRNAs (miRNAs). In addition, 58 PM-responsive lncRNAs were predicted as targets of miRNAs and one PM-responsive lncRNA was predicted as an endogenous target mimic (eTM). Furthermore, a total of 5,200 potential cis target genes and 5,625 potential trans target genes were predicted for PM-responsive lncRNAs. Functional enrichment analysis showed that these potential target genes are involved in different biological processes, such as the plant-pathogen interaction pathway, MAPK signaling pathway, and plant hormone signal transduction pathway. Taken together, this study provides a comprehensive view of C. pepo lncRNAs and explores the putative functions of PM-responsive lncRNAs, thus laying the foundation for further study of the regulatory mechanisms of lncRNAs responding to PM.

The Characters of Non-Coding RNAs and Their Biological Roles in Plant Development and Abiotic Stress Response

Plant growth and development are greatly affected by the environment. Many genes have been identified to be involved in regulating plant development and adaption of abiotic stress. Apart from protein-coding genes, more and more evidence indicates that non-coding RNAs (ncRNAs), including small RNAs and long ncRNAs (lncRNAs), can target plant developmental and stress-responsive mRNAs, regulatory genes, DNA regulatory regions, and proteins to regulate the transcription of various genes at the transcriptional, posttranscriptional, and epigenetic level. Currently, the molecular regulatory mechanisms of sRNAs and lncRNAs controlling plant development and abiotic response are being deeply explored. In this review, we summarize the recent research progress of small RNAs and lncRNAs in plants, focusing on the signal factors, expression characters, targets functions, and interplay network of ncRNAs and their targets in plant development and abiotic stress responses. The complex molecular regulatory pathways among small RNAs, lncRNAs, and targets in plants are also discussed. Understanding molecular mechanisms and functional implications of ncRNAs in various abiotic stress responses and development will benefit us in regard to the use of ncRNAs as potential character-determining factors in molecular plant breeding.

Analysis of NIA and GSNOR family genes and nitric oxide homeostasis in response to wheat-leaf rust interaction

Nitric oxide (NO) modulates plant response to biotic and abiotic stresses by S-nitrosylation-mediated protein post-translational modification. Nitrate reductase (NR) and S-nitrosoglutathione reductase (GSNOR) enzymes are essential for NO synthesis and the maintenance of Nitric oxide/S-nitroso glutathione (NO/GSNO) homeostasis, respectively. S-nitrosoglutathione, formed by the S-nitrosylation reaction of NO with glutathione, plays a significant physiological role as the mobile reservoir of NO. The genome-wide analysis identified nine NR (NIA) and three GSNOR genes in the wheat genome. Phylogenic analysis revealed that the nine NIA genes +were clustered into four groups and the 3 GSNORs into two groups. qRT-PCR expression profiling of NIAs and GSNORs was done in Chinese spring (CS), a leaf rust susceptible wheat line showing compatible interaction, and Transfer (TR), leaf rust-resistant wheat line showing incompatible interaction, post-inoculation with leaf rust pathotype 77-5 (121-R-63). All the NIA genes showed upregulation during incompatible interaction in comparison with the compatible reaction. The GSNOR genes showed a variable pattern of expression: the TaGSNOR1 showed little change, whereas TaGSNOR2 showed higher expression during the incompatible response. TaGSNOR3 showed a rise of expression both in compatible and incompatible reactions. Before inoculation and after 72 h of pathogen inoculation, NO localization was studied in both compatible and incompatible reactions. The S-nitrosothiol accumulation, NR, and glutathione reductase activity showed a consistent increase in the incompatible interactions. The results demonstrate that both NR and GSNOR plays significant role in defence against the leaf rust pathogen in wheat by modulating NO homeostasis or signalling.

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

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

Identification and functional deciphering suggested the regulatory roles of long intergenic ncRNAs (lincRNAs) in increasing grafting pepper resistance to Phytophthora capsici

BACKGROUND

As a popular and valuable technique, grafting is widely used to protect against soil-borne diseases and nematodes in vegetable production. Growing evidences have revealed that long intergenic ncRNAs (lincRNAs) are strictly regulated and play essential roles in plants development and stress responses. Nevertheless, genome-wide identification and function deciphering of pepper lincRNAs, especially for their roles in improving grafting pepper resistance to Phytophthora capsici is largely unknown.

RESULTS

In this study, RNA-seq data of grafting and control pepper plants with or without P. capsici inoculation were used to identify lincRNAs. In total, 2,388 reliable lincRNAs were identified. They were relatively longer and contained few exons than protein-coding genes. Similar to coding genes, lincRNAs had higher densities in euchromatin regions; and longer chromosome transcribed more lincRNAs. Expression pattern profiling suggested that lincRNAs commonly had lower expression than mRNAs. Totally, 607 differentially expressed lincRNAs (DE-lincRANs) were identified, of which 172 were found between P. capsici resistance grafting pepper sample GR and susceptible sample LDS. The neighboring genes of DE-lincRNAs and miRNAs competitively sponged by DE-lincRNAs were identified. Subsequently, the expression level of DE-lincRNAs was further confirmed by qRT-PCR and regulation patterns between DE-lincRNAs and neighboring mRNAs were also validated. Function annotation revealed that DE-lincRNAs increased the resistance of grafting prepper to P. capsici by modulating the expression of disease-defense related genes through cis-regulating and/or lincRNA-miRNA-mRNA interaction networks.

CONCLUSIONS

This study identified pepper lincRNAs and suggested their potential roles in increasing the resistance level of grafting pepper to P. capsici.

Mining plant endogenous target mimics from miRNA-lncRNA interactions based on dual-path parallel ensemble pruning method

The interactions between microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) play important roles in biological activities. Specially, lncRNAs as endogenous target mimics (eTMs) can bind miRNAs to regulate the expressions of target messenger RNAs (mRNAs). A growing number of studies focus on animals, but the studies on plants are scarce and many functions of plant eTMs are unknown. This study proposes a novel ensemble pruning protocol for predicting plant miRNA-lncRNA interactions at first. It adaptively prunes the base models based on dual-path parallel ensemble method to meet the challenge of cross-species prediction. Then potential eTMs are mined from predicted results. The expression levels of RNAs are identified through biological experiment to construct the lncRNA-miRNA-mRNA regulatory network, and the functions of potential eTMs are inferred through enrichment analysis. Experiment results show that the proposed protocol outperforms existing methods and state-of-the-art predictors on various plant species. A total of 17 potential eTMs are verified by biological experiment to involve in 22 regulations, and 14 potential eTMs are inferred by Gene Ontology enrichment analysis to involve in 63 functions, which is significant for further research.

Fungi Inhabiting the Wheat Endosphere

Wheat production is influenced by changing environmental conditions, including climatic conditions, which results in the changing composition of microorganisms interacting with this cereal. The group of these microorganisms includes not only endophytic fungi associated with the wheat endosphere, both pathogenic and symbiotic, but also those with yet unrecognized functions and consequences for wheat. This paper reviews the literature in the context of the general characteristics of endophytic fungi inhabiting the internal tissues of wheat. In addition, the importance of epigenetic regulation in wheat-fungus interactions is recognized and the current state of knowledge is demonstrated. The possibilities of using symbiotic endophytic fungi in modern agronomy and wheat cultivation are also proposed. The fact that the current understanding of fungal endophytes in wheat is based on a rather small set of experimental conditions, including wheat genotypes, plant organs, plant tissues, plant development stage, or environmental conditions, is recognized. In addition, most of the research to date has been based on culture-dependent methods that exclude biotrophic and slow-growing species and favor the detection of fast-growing fungi. Additionally, only a few reports of studies on the entire wheat microbiome using high-throughput sequencing techniques exist. Conducting comprehensive research on the mycobiome of the endosphere of wheat, mainly in the context of the possibility of using this knowledge to improve the methods of wheat management, mainly the productivity and health of this cereal, is needed.

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.

Genome-wide profiling of long noncoding RNAs involved in wheat spike development

BACKGROUND

Long noncoding RNAs (lncRNAs) have been shown to play important roles in the regulation of plant growth and development. Recent transcriptomic analyses have revealed the gene expression profiling in wheat spike development, however, the possible regulatory roles of lncRNAs in wheat spike morphogenesis remain largely unclear.

RESULTS

Here, we analyzed the genome-wide profiling of lncRNAs during wheat spike development at six stages, and identified a total of 8,889 expressed lncRNAs, among which 2,753 were differentially expressed lncRNAs (DE lncRNAs) at various developmental stages. Three hundred fifteen differentially expressed cis- and trans-regulatory lncRNA-mRNA pairs comprised of 205 lncRNAs and 279 genes were predicted, which were found to be mainly involved in the stress responses, transcriptional and enzymatic regulations. Moreover, the 145 DE lncRNAs were predicted as putative precursors or target mimics of miRNAs. Finally, we identified the important lncRNAs that participate in spike development by potentially targeting stress response genes, TF genes or miRNAs.

CONCLUSIONS

This study outlines an overall view of lncRNAs and their possible regulatory networks during wheat spike development, which also provides an alternative resource for genetic manipulation of wheat spike architecture and thus yield.

Investigation of long non-coding RNAs as regulatory players of grapevine response to powdery and downy mildew infection

BACKGROUND

Long non-coding RNAs (lncRNAs) are regulatory transcripts of length > 200 nt. Owing to the rapidly progressing RNA-sequencing technologies, lncRNAs are emerging as considerable nodes in the plant antifungal defense networks. Therefore, we investigated their role in Vitis vinifera (grapevine) in response to obligate biotrophic fungal phytopathogens, Erysiphe necator (powdery mildew, PM) and Plasmopara viticola (downy mildew, DM), which impose huge agro-economic burden on grape-growers worldwide.

RESULTS

Using computational approach based on RNA-seq data, 71 PM- and 83 DM-responsive V. vinifera lncRNAs were identified and comprehensively examined for their putative functional roles in plant defense response. V. vinifera protein coding sequences (CDS) were also profiled based on expression levels, and 1037 PM-responsive and 670 DM-responsive CDS were identified. Next, co-expression analysis-based functional annotation revealed their association with gene ontology (GO) terms for 'response to stress', 'response to biotic stimulus', 'immune system process', etc. Further investigation based on analysis of domains, enzyme classification, pathways enrichment, transcription factors (TFs), interactions with microRNAs (miRNAs), and real-time quantitative PCR of lncRNAs and co-expressing CDS pairs suggested their involvement in modulation of basal and specific defense responses such as: Ca²⁺-dependent signaling, cell wall reinforcement, reactive oxygen species metabolism, pathogenesis related proteins accumulation, phytohormonal signal transduction, and secondary metabolism.

CONCLUSIONS

Overall, the identified lncRNAs provide insights into the underlying intricacy of grapevine transcriptional reprogramming/post-transcriptional regulation to delay or seize the living cell-dependent pathogen growth. Therefore, in addition to defense-responsive genes such as TFs, the identified lncRNAs can be further examined and leveraged to candidates for biotechnological improvement/breeding to enhance fungal stress resistance in this susceptible fruit crop of economic and nutritional importance.

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.

Genome-Wide Identification and Characterization of Potato Long Non-coding RNAs Associated With Phytophthora infestans Resistance

Long non-coding RNA (lncRNA) is a crucial regulatory mechanism in the plant response to biotic and abiotic stress. However, their roles in potato (Solanum tuberosum L.) resistance to Phytophthora infestans (P. infestans) largely remain unknown. In this study, we identify 2857 lncRNAs and 33,150 mRNAs of the potato from large-scale published RNA sequencing data. Characteristic analysis indicates a similar distribution pattern of lncRNAs and mRNAs on the potato chromosomes, and the mRNAs were longer and had more exons than lncRNAs. Identification of alternative splicing (AS) shows that there were a total of 2491 lncRNAs generated from AS and the highest frequency (46.49%) of alternative acceptors (AA). We performed R package TCseq to cluster 133 specific differentially expressed lncRNAs from resistance lines and found that the lncRNAs of cluster 2 were upregulated. The lncRNA targets were subject to KEGG pathway enrichment analysis, and the interactive network between lncRNAs and mRNAs was constructed by using GENIE3, a random forest machine learning algorithm. Transient overexpression of StLNC0004 in Nicotiana benthamiana significantly suppresses P. infestans growth compared with a control, and the expression of extensin (NbEXT), the ortholog of the StLNC0004 target gene, was significantly upregulated in the overexpression line. Together, these results suggest that lncRNAs play potential functional roles in the potato response to P. infestans infection.

PacBio full-length transcriptome of wild apple (Malus sieversii) provides insights into canker disease dynamic response

BACKGROUND

Valsa canker is a serious disease in the stem of Malus sieversii, caused by Valsa mali. However, little is known about the global response mechanism in M. sieversii to V. mali infection.

RESULTS

Phytohormone jasmonic acid (JA) and salicylic acid (SA) profiles and transcriptome analysis were used to elaborate on the dynamic response mechanism. We determined that the JA was initially produced to respond to the necrotrophic pathogen V. mali infection at the early response stage, then get synergistically transduced with SA to respond at the late response stage. Furthermore, we adopted Pacific Biosciences (PacBio) full-length sequencing to identify differentially expressed transcripts (DETs) during the canker response stage. We obtained 52,538 full-length transcripts, of which 8139 were DETs. Total 1336 lncRNAs, 23,737 alternative polyadenylation (APA) sites and 3780 putative transcription factors (TFs) were identified. Additionally, functional annotation analysis of DETs indicated that the wild apple response to the infection of V. mali involves plant-pathogen interaction, plant hormone signal transduction, flavonoid biosynthesis, and phenylpropanoid biosynthesis. The co-expression network of the differentially expressed TFs revealed 264 candidate TF transcripts. Among these candidates, the WRKY family was the most abundant. The MsWRKY7 and MsWRKY33 were highly correlated at the early response stage, and MsWRKY6, MsWRKY7, MsWRKY19, MsWRKY33, MsWRKY40, MsWRKY45, MsWRKY51, MsWRKY61, MsWRKY75 were highly correlated at the late stage.

CONCLUSIONS

The full-length transcriptomic analysis revealed a series of immune responsive events in M. sieversii in response to V. mali infection. The phytohormone signal pathway regulatory played an important role in the response stage. Additionally, the enriched disease resistance pathways and differentially expressed TFs dynamics collectively contributed to the immune response. This study provides valuable insights into a dynamic response in M. sieversii upon the necrotrophic pathogen V. mali infection, facilitates understanding of response mechanisms to canker disease for apple, and provides supports in the identification of potential resistance genes in M. sieversii.

Interactions and links among the noncoding RNAs in plants under stresses

The complex interplay among sRNAs, lncRNAs and circRNAs has been implicated in plants under biotic and abiotic stresses. Here, we review current advances in our understanding of ncRNA interactions and links, which have considerable potential for improving the agronomic traits and the environmental adaptability of plants. Plants can respond to biotic or abiotic stresses. To cope with various conditions, numerous intricate molecular regulatory mechanisms have evolved in plants. Noncoding RNAs (ncRNAs) can be divided into small noncoding RNAs (sRNAs), long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs). Emerging evidence has demonstrated that interplay among the ncRNAs acts as a novel layer in the regulatory mechanisms, which has attracted substantial interest. Links between sRNAs can affect plant immune responses and development in synergistic or antagonistic manners. Additionally, multiple interactions between lncRNAs and sRNAs are involved in crop breeding, disease resistance and high tolerance to environmental stresses. Here, we summarize current knowledge of the interactions and links among the ncRNAs in plant responses to stresses and the methods for identifying ncRNA interactions. Furthermore, challenges and prospects for further progress in elucidating ncRNA interactions and links are highlighted.

Screening and functional characterization of candidate resistance genes to powdery mildew from Dasypyrum villosum#4 in a wheat line Pm97033

KEY MESSAGE

Three genes designated DvLox, Pm21#4, and Pm21#4-H identified in a wheat-Dasypyrum villosum#4 T6V#4S·6DL translocation line Pm97033 conferred wheat for powdery mildew resistance. Powdery mildew (PM) caused by Blumeria graminis f. sp. tritici (Bgt) is one of the most devastating diseases in wheat. To date, only a few genes conferring resistance to wheat PM are cloned. Dasypyrum villosum is a wild relative of wheat, which provides Pm21 conferring wheat immunity to PM. In this study, we obtained many differentially expressed genes (DEGs) from a wheat-D. villosum#4 T6V#4S·6DL translocation line Pm97033 using RNA-sequencing. Among them, 7 DEGs associated with pathogen resistance were up-regulated in front of Bgt infection. Virus-induced gene silencing and transformation assays demonstrated that two of them, DvLox and Pm21#4 encoding a lipoxygenase and a encoding coiled-coil/nucleotide-binding site/leucine-rich repeat resistance protein, conferred wheat PM resistance. The transgenic wheat plants expressing DvLox enhanced PM resistance, and the transgenic wheat plants expressing Pm21#4 showed PM immunity. The Pm21#4-silenced Pm97033 plants by the cluster regularly interspaced short palindromic repeats-associated endonuclease (CRISPR/Cas9) system were susceptible to PM. Thus, Pm21#4 is a key gene contributing PM immune resistance in Pm97033. Constitutively expression of Pm21#4-H, which is silenced in Pm97033 and D. villosum#4, endowed a PM-susceptible wheat variety Fielder with PM immune resistance.

Genome-Wide Identification and Characterization of Fusarium graminearum-Responsive lncRNAs in Triticum aestivum

Although the war between wheat and Fusarium has been widely investigated for years, long noncoding RNAs (lncRNAs), which have been proven to regulate important processes in the development and stress responses of plants, are still poorly known in wheat against Fusarium. Herein, we systematically reveal the roles of wheat lncRNAs in the process of Fusarium graminearum infection by high-throughput RNA sequencing. Well over 4130 of the total 4276 differentially expressed lncRNAs were already specifically expressed at 12 h postinoculation (hpi), but only 89 of these were specifically expressed at 24 hpi, indicating that the initial stage was the crucial stage for lncRNA-mediated gene regulation of wheat defense against F. graminearum. Target analysis showed the lncRNAs participated in various biological stress processes and had exclusive regulation models at different infection stages. Further H₂O₂ accumulation and protein ubiquitination assays supported this idea. Moreover, two lncRNAs (XLOC_302848 and XLOC_321638) were identified as Fusarium seedling blight resistance candidates by lncRNA-target expression pattern validation, and two lncRNAs (XLOC_113815, XLOC_123624) were Fusarium head blight resistance potential regulators by cross-validating the RNAseq data with the refined meta-QTL of wheat FHB resistance. These findings extend our knowledge on wheat lncRNAs response to F. graminearum attack and provide new insights for the functional and molecular research of future interactions between wheat and Fusarium.

Systematic Identification and Analysis of Light-Responsive Circular RNA and Co-expression Networks in Lettuce (Lactuca sativa)

Circular RNA (circRNA) is a covalently-closed single-stranded RNA molecule that plays an important role in transcriptional regulation of gene expression in a variety of species. Light intensity is a pivotal environmental factor affecting plant growth and development. However, little is known regarding photoresponsive plant circRNAs. Here, we aimed to investigate the expression and function of circRNAs in lettuce leaves in response to different light intensity treatments. We performed RNA sequencing (RNA-Seq) on leaves of lettuce (Lactuca sativa) to determine circRNA expression profiles and reverse-transcription polymerase chain reaction (PCR) to validate the candidate circRNA molecules. We then combined bioinformatics approach to explore the function of the parental genes of circRNA, including network, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes pathway analysis. We identified 1650 circRNAs in lettuce, of which 1508 (86.40%) were derived from exons. Using real-time PCR, we characterized 10 validated differentially expressed circRNAs and their parental genes, all of which showed expression patterns consistent with RNA-Seq data. Interestingly, the expression of circRNA was, in some cases, inversely correlated with the expression of the parental gene. Furthermore, analysis of the circRNA-microRNA-mRNA network suggests that circRNAs may be involved in plant hormone signaling and chlorophyll metabolism during photoreactivity. These findings provide an essential reference basis for studying circRNAs' biological mechanisms in light-treated plants.

The Emerging Role of Long Non-Coding RNAs in Plant Defense Against Fungal Stress

Growing interest and recent evidence have identified long non-coding RNA (lncRNA) as the potential regulatory elements for eukaryotes. LncRNAs can activate various transcriptional and post-transcriptional events that impact cellular functions though multiple regulatory functions. Recently, a large number of lncRNAs have also been identified in higher plants, and an understanding of their functional role in plant resistance to infection is just emerging. Here, we focus on their identification in crop plant, and discuss their potential regulatory functions and lncRNA-miRNA-mRNA network in plant pathogen stress responses, referring to possible examples in a model plant. The knowledge gained from a deeper understanding of this colossal special group of plant lncRNAs will help in the biotechnological improvement of crops.

LncRNA39026 Enhances Tomato Resistance to Phytophthora infestans by Decoying miR168a and Inducing PR Gene Expression

Our previous study has indicated that a long noncoding RNA (lncRNA), lncRNA39026, can be responsive to Phytophthora infestans infection. However, the function and regulation mechanism of lncRNA39026 during tomato resistance to P. infestans are unknown. In this study, an lncRNA39026 sequence was cloned from tomato Zaofen No. 2, and this sequence contained an endogenous target mimicry for miR168a, which might suppress the expression of miR168a. LncRNA39026 was strongly downregulated at 3 h in the tomato plants infected with P. infestans, and its expression level displayed a negative correlation with the expression level of miR168a and a positive correlation with the expression levels of SlAGO1 genes (target gene of miR168a) upon P. infestans infection. Tomato plants in which lncRNA39026 was overexpressed displayed enhanced resistance to P. infestans, decreased level of miR168a, and increased level of SlAGO1, whereas the resistance was impaired, level of miR168a was increased, and level of SlAGO1 was decreased after lncRNA39026 silencing. In addition, lncRNA39026 could also induce the expression of pathogenesis-related (PR) genes, as shown by increased and decreased expression levels of PR genes in tomato plants with overexpressed and silenced lncRNA39026, respectively. The result demonstrated that lncRNA39026 might function to decoy miR168a and affect the expression of PR genes in tomato plants, increasing resistance to disease.

Genome-wide analysis of the transcriptional response to drought stress in root and leaf of common bean

Genes related to the response to drought stress in leaf and root tissue of drought-susceptible (DS) and tolerant (DT) genotypes were characterized by RNA-Seq. In total, 54,750 transcripts, representative of 28,590 genes, were identified; of these, 1,648 were of high-fidelity (merge of 12 libraries) and described for the first time in the Andean germplasm. From the 1,239 differentially expressed genes (DEGs), 458 were identified in DT, with a predominance of genes in categories of oxidative stress, response to stimulus and kinase activity. Most genes related to oxidation-reduction terms in roots were early triggered in DT (T75) compared to DS (T150) suggestive of a mechanism of tolerance by reducing the damage from ROS. Among the KEGG enriched by DEGs up-regulated in DT leaves, two related to the formation of Sulfur-containing compounds, which are known for their involvement in tolerance to abiotic stresses, were common to all treatments. Through qPCR, 88.64% of the DEGs were validated. A total of 151,283 variants were identified and functional effects estimated for 85,780. The raw data files were submitted to the NCBI database. A transcriptome map revealed new genes and isoforms under drought. These results supports a better understanding of the drought tolerance mechanisms in beans.

Genome-Wide Identification of lncRNAs During Rice Seed Development

Rice seed is a pivotal reproductive organ that directly determines yield and quality. Long non-coding RNAs (lncRNAs) have been recognized as key regulators in plant development, but the roles of lncRNAs in rice seed development remain unclear. In this study, we performed a paired-end RNA sequencing in samples of rice pistils and seeds at three and seven days after pollination (DAP) respectively. A total of 540 lncRNAs were obtained, among which 482 lncRNAs had significantly different expression patterns during seed development. Results from semi-qPCR conducted on 15 randomly selected differentially expressed lncRNAs suggested high reliability of the transcriptomic data. RNA interference of TCONS_00023703, which is predominantly transcribed in developing seeds, significantly reduced grain length and thousand-grain weight. These results expanded the dataset of lncRNA in rice and enhanced our understanding of the biological functions of lncRNAs in rice seed development.

Comparative transcriptome analysis uncovers regulatory roles of long non-coding RNAs involved in resistance to powdery mildew in melon

Background

Long non-coding RNAs (lncRNAs) are a class of non-coding RNAs with more than 200 nucleotides in length, which play vital roles in a wide range of biological processes. Powdery mildew disease (PM) has become a major threat to the production of melon. To investigate the potential roles of lncRNAs in resisting to PM in melon, it is necessary to identify lncRNAs and uncover their molecular functions. In this study, we compared the lncRNAs between a resistant and a susceptible melon in response to PM infection.

Results

It is reported that 11,612 lncRNAs were discovered, which were distributed across all 12 melon chromosomes, and > 85% were from intergenic regions. The melon lncRNAs have shorter transcript lengths and fewer exon numbers than protein-coding genes. In addition, a total of 407 and 611 lncRNAs were found to be differentially expressed after PM infection in PM-susceptible and PM-resistant melons, respectively. Furthermore, 1232 putative targets of differently expressed lncRNAs (DELs) were discovered and gene ontology enrichment (GO) analysis showed that these target genes were mainly enriched in stress-related terms. Consequently, co-expression patterns between LNC_018800 and CmWRKY21, LNC_018062 and MELO3C015771 (glutathione reductase coding gene), LNC_014937 and CmMLO5 were confirmed by qRT-PCR. Moreover, we also identified 24 lncRNAs that act as microRNA (miRNA) precursors, 43 lncRNAs as potential targets of 22 miRNA families and 13 lncRNAs as endogenous target mimics (eTMs) for 11 miRNAs.

Conclusion

This study shows the first characterization of lncRNAs involved in PM resistance in melon and provides a starting point for further investigation into the functions and regulatory mechanisms of lncRNAs in the resistance to PM.

Genome-Wide Identification of lncRNAs and Analysis of ceRNA Networks During Tomato Resistance to Phytophthora infestans

Our previous studies have revealed the function of long noncoding RNAs (lncRNAs) and microRNAs (miRNAs) in tomato in response to Phytophthora infestans infection. However, the interaction relationships between lncRNAs and miRNAs during tomato resistance to P. infestans infection are unknown. In this study, 9,011 lncRNAs were identified from tomato plants, including 115 upregulated and 81 downregulated lncRNAs. Among these, 148 were found to be differentially expressed and might affect the expression of 771 genes, which are composed of 887 matched lncRNA-mRNA pairs. In total, 88 lncRNAs were identified as endogenous RNAs (ceRNAs) and predicted to decoy 46 miRNAs. Degradome sequencing revealed that 11 miRNAs that were decoyed by 20 lncRNAs could target 30 genes. These lncRNAs, miRNAs, and target genes were predicted to form 10 regulatory modules. Among them, lncRNA42705/lncRNA08711, lncRNA39896, and lncRNA11265/lncRNA15816 might modulate MYB, HD-Zip, and NAC transcription factors by decoying miR159, miR166b, and miR164a-5p, respectively. Upon P. infestans infection, the expression levels of lncRNA42705 and lncRNA08711 displayed a negative correlation with the expression level of miR159 and a positive correlation with the expression levels of MYB genes. Tomato plants in which lncRNA42705 and lncRNA08711 were silenced displayed increased levels of miR159 and decreased levels of MYB, respectively. The result demonstrated that lncRNAs might function as ceRNAs to decoy miRNAs and affect their target genes in tomato plants, increasing resistance to disease.

Identification of lncRNAs and their regulatory relationships with target genes and corresponding miRNAs in melon response to powdery mildew fungi

Melon (Cucumis melo L.), an economically beneficial crop widely cultivated around the world, is vulnerable to powdery mildew (PM). However, the studies on molecular mechanism of melon response to PM fungi is still limited. Long non coding RNAs (lncRNAs) have emerged as new regulators in plants response to biotic stresses. We predicted and identified the intricate regulatory roles of lncRNAs in melon response to PM fungi. A total of 539 lncRNAs were identified from PM-resistant (MR-1) and susceptible melon (Top Mark), in which 254 were significantly altered after PM fungi infection. Multiple target genes of lncRNAs were found to be involved in the hydrolysis of chitin, callose deposition and cell wall thickening, plant-pathogen interaction and plant hormone signal transduction pathway. Additionally, a total of 42 lncRNAs possess the various functions with microRNAs (miRNAs), including lncRNAs that are targeted by miRNAs and function as miRNA precursors or miRNA sponges. These findings provide a comprehensive view of potentially functional lncRNAs, corresponding target genes and related lncRNA-miRNA pairs, which will greatly increase our knowledge of the mechanism underlying susceptibility and resistance to PM in melon.

Expressivity of the key genes associated with seed and pod development is highly regulated via lncRNAs and miRNAs in Pigeonpea

Non-coding RNA’s like miRNA, lncRNA, have gained immense importance as a significant regulatory factor in different physiological and developmental processes in plants. In an effort to understand the molecular role of these regulatory agents, in the present study, 3019 lncRNAs and 227 miRNAs were identified from different seed and pod developmental stages in Pigeonpea, a major grain legume of Southeast Asia and Africa. Target analysis revealed that 3768 mRNAs, including 83 TFs were targeted by lncRNAs; whereas 3060 mRNA, including 154 TFs, were targeted by miRNAs. The targeted transcription factors majorly belong to WRKY, MYB, bHLH, etc. families; whereas the targeted genes were associated with the embryo, seed, and flower development. Total 302 lncRNAs interact with miRNAs and formed endogenous target mimics (eTMs) which leads to sequestering of the miRNAs present in the cell. Expression analysis showed that notably, Cc_lncRNA-2830 expression is up-regulated and sequestrates miR160h in pod leading to higher expression of the miR160h target gene, Auxin responsive factor-18. A similar pattern was observed for SPIKE, Auxin signaling F-box-2, Bidirectional sugar transporter, and Starch synthetase-2 eTMs. All the identified target mRNAs code for transcription factor and genes are involved in the processes like cell division, plant growth and development, starch synthesis, sugar transportation and accumulation of storage proteins which are essential for seed and pod development. On a combinatorial basis, our study provides a lncRNA and miRNA based regulatory insight into the genes governing seed and pod development in Pigeonpea.

Maize transposable elements contribute to long non-coding RNAs that are regulatory hubs for abiotic stress response

Background

Several studies have mined short-read RNA sequencing datasets to identify long non-coding RNAs (lncRNAs), and others have focused on the function of individual lncRNAs in abiotic stress response. However, our understanding of the complement, function and origin of lncRNAs – and especially transposon derived lncRNAs (TE-lncRNAs) - in response to abiotic stress is still in its infancy.

Results

We utilized a dataset of 127 RNA sequencing samples that included total RNA datasets and PacBio fl-cDNA data to discover lncRNAs in maize. Overall, we identified 23,309 candidate lncRNAs from polyA+ and total RNA samples, with a strong discovery bias within total RNA. The majority (65%) of the 23,309 lncRNAs had sequence similarity to transposable elements (TEs). Most had similarity to long-terminal-repeat retrotransposons from the Copia and Gypsy superfamilies, reflecting a high proportion of these elements in the genome. However, DNA transposons were enriched for lncRNAs relative to their genomic representation by ~ 2-fold. By assessing the fraction of lncRNAs that respond to abiotic stresses like heat, cold, salt and drought, we identified 1077 differentially expressed lncRNA transcripts, including 509 TE-lncRNAs. In general, the expression of these lncRNAs was significantly correlated with their nearest gene. By inferring co-expression networks across our large dataset, we found that 39 lncRNAs are as major hubs in co-expression networks that respond to abiotic stress, and 18 appear to be derived from TEs.

Conclusions

Our results show that lncRNAs are enriched in total RNA samples, that most (65%) are derived from TEs, that at least 1077 are differentially expressed during abiotic stress, and that 39 are hubs in co-expression networks, including a small number that are evolutionary conserved. These results suggest that lncRNAs, including TE-lncRNAs, may play key regulatory roles in moderating abiotic responses.

Transcriptome-wide alternative splicing modulation during plant-pathogen interactions in wheat

Alternative splicing (AS) enhances the diversities of both transcripts and proteins in eukaryotes, which contribute to stress adaptation. To catalog wheat (Triticum aestivum L.) AS genes, we characterized 45 RNA-seq libraries from wheat seedlings infected by powdery mildew, Blumeria graminis f. sp. tritici (Bgt) or stripe rust fungus, Puccinia striiformis f. sp. tritici (Pst). We discovered that 11.2% and 10.4% of the multiexon genes had AS transcripts during Bgt and Pst infections, respectively. In response to fungal infection, wheat modulated AS not only in disease resistance proteins, but also in splicing related factors. Apart from the stress induced or activated splicing variants by pathogen, the differential expression profiles were fold increased through changing the ratio of full spliced transcripts versus intron retention (IR) transcripts. Comparing AS transcripts produced by the same gene in Bgt with Pst stress, the spliced terminal exons and the stranded introns are independent and different. This demonstrated that differential induction of specific splice variants were activated against two fungal pathogens. The specific induced AS genes in the Pst-resistant plants were enriched in improving the membrane permeability and protein modification ability, whereas gene expression involved in protein translation and transport were strengthened in Pst-susceptible plants.

Continuous salt stress-induced long non-coding RNAs and DNA methylation patterns in soybean roots

BACKGROUND

Environmental stimuli can activate a series of physiological and biochemical responses in plants accompanied by extensive transcriptional reprogramming. Long non-coding RNAs (lncRNAs), as versatile regulators, control gene expression in multiple ways and participate in the adaptation to biotic and abiotic stresses.

RESULTS

In this study, soybean seedlings were continuously cultured for 15 days with high salinity solutions started from seed germination. Strand-specific whole transcriptome sequencing and stringent bioinformatic analysis led to the identification of 3030 long intergenic non-coding RNAs (lincRNAs) and 275 natural antisense transcripts (lncNATs) in soybean roots. In contrast to mRNAs, newly identified lncRNAs exhibited less exons, similar AU content to UTRs, even distribution across the genome and low evolutionary conservation. Remarkably, more than 75% of discovered lncRNAs that were activated or up-regulated by continuous salt stress mainly targeted proteins with binding and catalytic activities. Furthermore, two DNA methylation maps with single-base resolution were generated by using reduced representation bisulfite sequencing, offering a genome-wide perspective and important clues for epigenetic regulation of stress-associated lncRNAs and protein-coding genes.

CONCLUSIONS

Taken together, our findings systematically demonstrated the characteristics of continuous salt stress-induced lncRNAs and extended the knowledge of corresponding methylation profiling, providing valuable evidence for a better understanding of how plants cope with long-term salt stress circumstances.

Identification of Long Non-Coding RNAs and the Regulatory Network Responsive to Arbuscular Mycorrhizal Fungi Colonization in Maize Roots

Recently, long noncoding RNAs (lncRNAs) have emerged as vital regulators of many biological processes in animals and plants. However, to our knowledge no investigations on plant lncRNAs which respond to arbuscular mycorrhizal (AM) fungi have been reported thus far. In this study, maize roots colonized with AM fungus were analyzed by strand-specific RNA-Seq to identify AM fungi-responsive lncRNAs and construct an associated regulatory network. A total of 1837 differentially expressed protein coding genes (DEGs) were identified from maize roots with Rhizophagus irregularis inoculation. Many AM fungi-responsive genes were homologs to MtPt4, STR, STR2, MtFatM, and enriched pathways such as fatty acid biosynthesis, response to phosphate starvation, and nitrogen metabolism are consistent with previous studies. In total, 5941 lncRNAs were identified, of which more than 3000 were new. Of those, 63 lncRNAs were differentially expressed. The putative target genes of differentially expressed lncRNAs (DELs) were mainly related to phosphate ion transmembrane transport, cellular response to potassium ion starvation, and lipid catabolic processes. Regulatory network analysis showed that DELs might be involved in the regulation of bidirectional nutrient exchange between plant and AM fungi as mimicry of microRNA targets. The results of this study can broaden our knowledge on the interaction between plant and AM fungi.

ceRNAs in plants: computational approaches and associated challenges for target mimic research

The competing endogenous RNA hypothesis has gained increasing attention as a potential global regulatory mechanism of microRNAs (miRNAs), and as a powerful tool to predict the function of many noncoding RNAs, including miRNAs themselves. Most studies have been focused on animals, although target mimic (TMs) discovery as well as important computational and experimental advances has been developed in plants over the past decade. Thus, our contribution summarizes recent progresses in computational approaches for research of miRNA:TM interactions. We divided this article in three main contributions. First, a general overview of research on TMs in plants is presented with practical descriptions of the available literature, tools, data, databases and computational reports. Second, we describe a common protocol for the computational and experimental analyses of TM. Third, we provide a bioinformatics approach for the prediction of TM motifs potentially cross-targeting both members within the same or from different miRNA families, based on the identification of consensus miRNA-binding sites from known TMs across sequenced genomes, transcriptomes and known miRNAs. This computational approach is promising because, in contrast to animals, miRNA families in plants are large with identical or similar members, several of which are also highly conserved. From the three consensus TM motifs found with our approach: MIM166, MIM171 and MIM159/319, the last one has found strong support on the recent experimental work by Reichel and Millar [Specificity of plant microRNA TMs: cross-targeting of mir159 and mir319. J Plant Physiol 2015;180:45-8]. Finally, we stress the discussion on the major computational and associated experimental challenges that have to be faced in future ceRNA studies.

Genome-wide identification and characterization of long non-coding RNAs responsive to Dickeya zeae in rice

Plant long non-coding RNA (lncRNA) is a type of newly emerging epigenetic regulator playing a critical role in plant growth, development, and biotic stress responses. However, it is unknown whether lncRNAs are involved in resistance responses between rice and Dickeya zeae, a bacterial agent causing rice foot rot disease. In this study, RNA-seq was performed to uncover the co-expression regulating networks mediated by D. zeae responsive lncRNAs and their candidate target genes. Of the 4709 lncRNAs identified, 2518 and 2191 were up- and down-regulated in response to D. zeae infection, respectively. Expression changes of 17 selected lncRNAs and their predicted targets with a potential role in defense response were investigated by qPCR. The expression levels of five lncRNAs were up-regulated while their cognate candidate target genes were down-regulated upon D. zeae infection. In addition, several lncRNAs were predicted to be target mimics of osa-miR396 and osa-miR156. These results suggest that lncRNAs might play a role in response to D. zeae infection by regulating the transcript levels of their targets and miRNAs in rice.

Integrated mRNA, sRNA, and degradome sequencing reveal oilseed rape complex responses to Sclerotinia sclerotiorum (Lib.) infection

Sclerotinia stem rot (SSR), caused by the fungal pathogen Sclerotinia sclerotiorum, is a devastating disease resulting in yield losses and decreases in seed quality in oilseed rape (Brassica napus) worldwide. However, the molecular mechanisms underlying the response of oilseed rape to S. sclerotiorum infection at the transcriptional and post-transcriptional levels are poorly understood. Here, we used an integrated omics approach (transcriptome, sRNAome, and degradome sequencing) on the Illumina platform to compare the RNA expression and post-transcriptional profiles of oilseed rape plants inoculated or not with S. sclerotiorum. In total, 7,065 differentially expressed genes (DEGs) compared with the mock-inoculated control at 48 hours post inoculation were identified. These DEGs were associated with protein kinases, signal transduction, transcription factors, hormones, pathogenesis-related proteins, secondary metabolism, and transport. In the sRNA-Seq analysis, 77 known and 176 novel miRNAs were identified; however, only 10 known and 41 novel miRNAs were differentially expressed between the samples inoculated or not with S. sclerotiorum. Degradome sequencing predicted 80 cleavage sites with 64 miRNAs. Integrated mRNA, sRNA and degradome sequencing analysis reveal oilseed rape complex responses to S. sclerotiorum infection. This study provides a global view of miRNA and mRNA expression profiles in oilseed rape following S. sclerotiorum infection.

Transcriptomic analyses of rice (Oryza sativa) genes and non-coding RNAs under nitrogen starvation using multiple omics technologies

BACKGROUND

Nitrogen (N) is a key macronutrient essential for plant growth, and its availability has a strong influence on crop development. The application of synthetic N fertilizers on crops has increased substantially in recent decades; however, the applied N is not fully utilized due to the low N use efficiency of crops. To overcome this limitation, it is important to understand the genome-wide responses and functions of key genes and potential regulatory factors in N metabolism.

RESULTS

Here, we characterized changes in the rice (Oryza sativa) transcriptome, including genes, newly identified putative long non-coding RNAs (lncRNAs), and microRNAs (miRNAs) and their target mRNAs in response to N starvation using four different transcriptome approaches. Analysis of rice genes involved in N metabolism and/or transport using strand-specific RNA-Seq identified 2588 novel putative lncRNA encoding loci. Analysis of previously published RNA-Seq datasets revealed a group of N starvation-responsive lncRNAs showing differential expression under other abiotic stress conditions. Poly A-primed sequencing (2P-Seq) revealed alternatively polyadenylated isoforms of N starvation-responsive lncRNAs and provided precise 3' end information on the transcript models of these lncRNAs. Analysis of small RNA-Seq data identified N starvation-responsive miRNAs and down-regulation of miR169 family members, causing de-repression of NF-YA, as confirmed by strand-specific RNA-Seq and qRT-PCR. Moreover, we profiled the N starvation-responsive down-regulation of root-specific miRNA, osa-miR444a.4-3p, and Degradome sequencing confirmed MADS25 as a novel target gene.

CONCLUSIONS

In this study, we used a combination of multiple RNA-Seq analyses to extensively profile the expression of genes, newly identified lncRNAs, and microRNAs in N-starved rice roots and shoots. Data generated in this study provide an in-depth understanding of the regulatory pathways modulated by N starvation-responsive miRNAs. The results of comprehensive, large-scale data analysis provide valuable information on multiple aspects of the rice transcriptome, which may be useful in understanding the responses of rice plants to changes in the N supply status of soil.

Development and comparative genomic mapping of Dasypyrum villosum 6V#4S-specific PCR markers using transcriptome data

Twenty-five Dasypyrum villosum 6V#4S-specific PCR markers were developed using transcriptome data and further assigned to comparative genomic maps of wheat chromosome 6A, 6B, and 6D and barley chromosome 6H contrasting their homologous genes in these genomes. Two Dasypyrum villosum accessions, D.v#2 and No. 1026 from England and Russia, respectively, contain Pm21 on chromosome 6V#2S and PmV on chromosome 6V#4S. Both genes confer high resistance to powdery mildew (PM) in wheat. Even though several molecular markers have been developed to detect Pm21 and PmV, only the MBH1 marker can simultaneously detect both Pm21 and PmV. In this study, we first used a high-throughput sequencing technique to obtain the transcriptome sequences of a wheat-D. villosum translocation line, Pm97033-which contains chromosome 6V#4S carrying the PmV locus, under wheat PM pathogen induction. Twenty-five 6V#4S chromosome-specific markers were developed. Three of them were able to clearly distinguish chromosomes 6V#4S and 6V#2S by product size, four amplified the product specific for chromosome 6V#4S only, and the remaining 18 markers identified chromosome 6VS in wheat backgrounds. Two different D. villosum accessions, their derived translocation lines and wheat varieties carrying different chromosome 6VS were identified using these specific markers. The 25 newly developed markers together with the known PM resistance gene Stpk-V were used to construct comparative genomic maps with the homoeologous chromosome arms of wheat and barley. The colinearity of the identified gene sequences amplified by the 25 markers among wheat chromosomes 6A, 6B, and 6D and barley chromosome 6H was not very conserved and interrupted frequently by inversion and insertion. Our markers have potential in marker assisted selection for PM resistance breeding, and for locating other potential important genes and cloning the PmV gene on chromosome 6V#4S.

Survey of High Throughput RNA-Seq Data Reveals Potential Roles for lncRNAs during Development and Stress Response in Bread Wheat

Long non-coding RNAs (lncRNAs) are a family of regulatory RNAs that play essential role in the various developmental processes and stress responses. Recent advances in sequencing technology and computational methods enabled identification and characterization of lncRNAs in certain plant species, but they are less known in Triticum aestivum (bread wheat). Herein, we analyzed 52 RNA seq data (>30 billion reads) and identified 44,698 lncRNAs in T. aestivum genome, which were characterized in comparison to the coding sequences (mRNAs). Similar to the mRNAs, lncRNAs were also derived from each sub-genome and chromosome, and showed tissue developmental stage specific and differential expression, as well. The modulated expression of lncRNAs during abiotic stresses like heat, drought, and salt indicated their putative role in stress response. The co-expression of lncRNAs with vital mRNAs including various transcription factors and enzymes involved in Abscisic acid (ABA) biosynthesis, and gene ontology mapping inferred their regulatory roles in numerous biological processes. A few lncRNAs were predicted as precursor (19 lncRNAs), while some as target mimics (1,047 lncRNAs) of known miRNAs involved in various regulatory functions. The results suggested numerous functions of lncRNAs in T. aestivum, and unfolded the opportunities for functional characterization of individual lncRNA in future studies.

Long noncoding miRNA gene represses wheat β-diketone waxes

The cuticle of terrestrial plants functions as a protective barrier against many biotic and abiotic stresses. In wheat and other Triticeae, β-diketone waxes are major components of the epicuticular layer leading to the bluish-white glaucous trait in reproductive-age plants. Glaucousness in durum wheat is controlled by a metabolic gene cluster at the WAX1 (W1) locus and a dominant suppressor INHIBITOR of WAX1 (Iw1) on chromosome 2B. The wheat D subgenome from progenitor Aegilops tauschii contains W2 and Iw2 paralogs on chromosome 2D. Here we identify the Iw1 gene from durum wheat and demonstrate the unique regulatory mechanism by which Iw1 acts to suppress a carboxylesterase-like protein gene, W1-COE, within the W1 multigene locus. Iw1 is a long noncoding RNA (lncRNA) containing an inverted repeat (IR) with >80% identity to W1-COE The Iw1 transcript forms a miRNA precursor-like long hairpin producing a 21-nt predominant miRNA, miRW1, and smaller numbers of related sRNAs associated with the nonglaucous phenotype. When Iw1 was introduced into glaucous bread wheat, miRW1 accumulated, W1-COE and its paralog W2-COE were down-regulated, and the phenotype was nonglaucous and β-diketone-depleted. The IR region of Iw1 has >94% identity to an IR region on chromosome 2 in Ae. tauschii that also produces miRW1 and lies within the marker-based location of Iw2 We propose the Iw loci arose from an inverted duplication of W1-COE and/or W2-COE in ancestral wheat to form evolutionarily young miRNA genes that act to repress the glaucous trait.