Rice lncRNAs Research Articles

Comprehensive transcriptomic analysis of three varieties with different brown planthopper-resistance identifies leaf sheath lncRNAs in rice

BackgroundLong non-coding RNAs (lncRNAs) have been brought great attention for their crucial roles in diverse biological processes. However, systematic identification of lncRNAs associated with specialized rice pest, brown planthopper (BPH), defense in rice remains unexplored.ResultsIn this study, a genome-wide high throughput sequencing analysis was performed using leaf sheaths of susceptible rice Taichung Native 1 (TN1) and resistant rice IR36 and R476 with and without BPH feeding. A total of 2283 lncRNAs were identified, of which 649 lncRNAs were differentially expressed. During BPH infestation, 84 (120 in total), 52 (70 in total) and 63 (94 in total) of differentially expressed lncRNAs were found only in TN1, IR36 and R476, respectively. Through analyzing their cis-, trans-, and target mimic-activities, not only the lncRNAs targeting resistance genes (NBS-LRR and RLKs) and transcription factors, but also the lncRNAs acting as the targets of the well-studied stress-related miRNAs (miR2118, miR528, and miR1320) in each variety were identified. Before the BPH feeding, 238 and 312 lncRNAs were found to be differentially expressed in TN1 vs. IR36 and TN1 vs. R476, respectively. Among their putative targets, the plant-pathogen interaction pathway was significantly enriched. It is speculated that the resistant rice was in a priming state by the regulation of lncRNAs. Furthermore, the lncRNAs extensively involved in response to BPH feeding were identified by Weighted Gene Co-expression Network Analysis (WGCNA), and the possible regulation networks of the key lncRNAs were constructed. These lncRNAs regulate different pathways that contribute to the basal defense and specific resistance of rice to the BPH.ConclusionIn summary, we identified the specific lncRNAs targeting the well-studied stress-related miRNAs, resistance genes, and transcription factors in each variety during BPH infestation. Additionally, the possible regulating network of the lncRNAs extensively responding to BPH feeding revealed by WGCNA were constructed. These findings will provide further understanding of the regulatory roles of lncRNAs in BPH defense, and lay a foundation for functional research on the candidate lncRNAs.

Comparative transcriptome-wide identification and differential expression of genes and lncRNAs in rice near-isogenic line (KW-Bph36-NIL) in response to BPH feeding.

Brown planthopper (BPH) is the most devastating pest of rice in Asia, causing substantial yield losses and has become a challenging task to be controlled under field conditions. Although extensive measures have been taken over the past decades, which resulted in the evolution of new resistant BPH strains. Therefore, besides other possible approaches, equipping host plants with resistant genes is the most effective and environment-friendly technique for BPH control. Here, we systematically analyzed transcriptome changes in the susceptible rice variety Kangwenqingzhan (KW) and the resistant near-isogenic line (NIL) KW-Bph36-NIL, through RNA-seq, depicting the differential expression profiles of mRNAs and long non-coding RNAs (lncRNAs) in rice before and after BPH feeding. We observed a proportion of genes (1.48%) and (2.74%) were altered in KW and NIL, respectively, indicating different responses of rice strains against BPH feeding. Nevertheless, we characterized 384 differentially expressed long non-coding RNAs (DELs) that can be impacted by the two strains by alternatively changing the expression patterns of the respective coding genes, suggesting their certain involvement in response to BPH feeding. In BPH invasion, KW and NIL responded differently by modifying the synthesis, storage, and transformation of intracellular substances, adjusting the nutrient accumulation and utilization inside and outside the cells. In addition, NIL expressed stronger resistance by acutely up-regulating genes and other transcription factors related to stress resistance and plant immunity. Altogether, our study elaborates valuable insights into the genome-wide DEGs and DELs expression profiles of rice under BPH invasion by high throughput sequencing and further suggests that NILs can be utilized in BPH resistance breeding programs in developing high-resistance rice lines.

Genome-wide identification and association analysis for virus-responsive lncRNAs in rice (Oryza sativa L.)

Genome-wide identification and association analysis for virus-responsive lncRNAs in rice (Oryza sativa L.)

LncRNA TCONS_00021861 is functionally associated with drought tolerance in rice (Oryza sativa L.) via competing endogenous RNA regulation

BackgroundWater deficit is an abiotic stress that retards plant growth and destabilizes crop production. Long non coding RNAs (lncRNAs) are a class of non-coding endogenous RNAs that participate in diverse cellular processes and stress responses in plants. lncRNAs could function as competing endogenous RNAs (ceRNA) and represent a novel layer of gene regulation. However, the regulatory mechanism of lncRNAs as ceRNA in drought stress response is yet unclear.ResultsIn this study, we performed transcriptome-wide identification of drought-responsive lncRNAs in rice. Thereafter, we constructed a lncRNA-mediated ceRNA network by analyzing competing relationships between mRNAs and lncRNAs based on ceRNA hypothesis. A drought responsive ceRNA network with 40 lncRNAs, 23 miRNAs and 103 mRNAs was obtained. Network analysis revealed TCONS_00021861/miR528-3p/YUCCA7 regulatory axis as a hub involved in drought response. The miRNA-target expression and interaction were validated by RT-qPCR and RLM-5’RACE. TCONS_00021861 showed significant positive correlation (r = 0.7102) with YUCCA7 and negative correlation with miR528-3p (r = -0.7483). Overexpression of TCONS_00021861 attenuated the repression of miR528-3p on YUCCA7, leading to increased IAA (Indole-3-acetic acid) content and auxin overproduction phenotypes.ConclusionsTCONS_00021861 could regulate YUCCA7 by sponging miR528-3p, which in turn activates IAA biosynthetic pathway and confer resistance to drought stress. Our findings provide a new perspective of the regulatory roles of lncRNAs as ceRNAs in drought resistance of rice.

Genome-Wide Analysis Identified a Set of Conserved lncRNAs Associated with Domestication-Related Traits in Rice

Crop domestication, which gives rise to a number of desirable agronomic traits, represents a typical model system of plant evolution. Numerous genomic evidence has proven that noncoding RNAs such as microRNAs and phasiRNAs, as well as protein-coding genes, are selected during crop domestication. However, limited data shows plant long noncoding RNAs (lncRNAs) are also involved in this biological process. In this study, we performed strand-specific RNA sequencing of cultivated rice Oryza sativa ssp. japonica and O. sativa ssp. indica, and their wild progenitor O. rufipogon. We identified a total of 8528 lncRNAs, including 4072 lncRNAs in O. rufipogon, 2091 lncRNAs in japonica rice, and 2365 lncRNAs in indica rice. The lncRNAs expressed in wild rice were revealed to be shorter in length and had fewer exon numbers when compared with lncRNAs from cultivated rice. We also identified a number of conserved lncRNAs in the wild and cultivated rice. The functional study demonstrated that several of these conserved lncRNAs are associated with domestication-related traits in rice. Our findings revealed the feature and conservation of lncRNAs during rice domestication and will further promote functional studies of lncRNAs in rice.

English

Rice seed development involves an intricate regulatory network that directly determines seed size and weight. Long noncoding RNAs (lncRNAs) have been defined as key regulators of gene expression involved in diverse biological processes. However, the function of lncRNAs in rice seed development is still poorly understood. We performed paired-end RNA sequencing of Nipponbare rice at 5, 10 and 15 DPA (days post anthesis) in two different environments (early and middle-season rice). A total of 382 lncRNAs were detected as differentially expressed among these stages, including 344 and 307 lncRNAs in early and middle-season rice, respectively, and 70.42% (269 of 382) of the lncRNAs were found in both environments. The results showed that environment had little effect on the expression of lncRNAs. Furthermore, there were 127, 172, and 31 DElncs (differentially expressed lncRNAs) and 154, 140, and 59 DElncs in early and middle-season rice, respectively, in comparisons of 10_DPA vs 5_DPA, 15_DPA vs 5_DPA and 15_DPA vs 10_DPA. This result indicated that the number and expression level of lncRNAs at 5 DAP were significantly different from those at 10 DAP and 15 DAP. Furthermore, GO pathway analysis of cis target genes of DElncs in 10_DPA vs 5_DPA and 15_DPA vs 5_DPA revealed that the significant GO pathways were extracellular region, nutrient reservoir activity and cell wall macromolecule catabolic process. Our study revealed dynamic expression of lncRNAs in three stages and systematically explored the differences in lncRNAs between early and middle-season rice, which could provide a valuable resource for future high-yield breeding. © 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers © 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers© 2021 Friends Science Publishers©

Identification and functional prediction of long non-coding RNAs of rice (Oryza sativa L.) at reproductive stage under salinity stress.

Salinity adversely affects the yield and growth of rice (Oryza sativa L.) plants severely, particularly at reproductive stage. Long non-coding RNAs (lncRNAs) are key regulators of diverse molecular and cellular processes in plants. Till now, no systematic study has been reported for regulatory roles of lncRNAs in rice under salinity at reproductive stage. In this study, total80 RNA-seq data of Horkuch (salt-tolerant) and IR-29 (salt-sensitive) genotypes of rice were used and found 1626 and 2208 transcripts as putative high confidence lncRNAs, among which 1529 and 2103 were found to be novel putative lncRNAs in root and leaf tissue respectively. In Horkuch and IR-29, 14 and 16 lncRNAs were differentially expressed in root tissue while 18 and 63 lncRNAs were differentially expressed in leaf tissue. Interaction analysis among the lncRNAs, miRNAs and corresponding mRNAs indicated that these modules are involved in different biochemical pathways e.g. phenyl propanoid pathway during salinity stress in rice. Interestingly, two differentially expressed lncRNAs such as TCONS_00008914 and TCONS_00008749 were found as putative target mimics of known rice miRNAs. This study indicates that lncRNAs are involved in salinity adaptation of rice at reproductive stage through certain biochemical pathways.

Genome-wide identification and integrated analysis of lncRNAs in rice backcross introgression lines (BC2F12)

BackgroundDistant hybridization is an important way to create interspecific genetic variation and breed new varieties in rice. A lot of backcross introgression lines (BILs) had been constructed for the scientific issues in rice. However, studies on the critical regulatory factor lncRNA in cultivated rice, wild rice and their BIL progenies were poorly reported.ResultsHere, high-throughput RNA sequencing technology was used to explore the functional characteristics and differences of lncRNAs in O. sativa, O. longistaminata and their three BC2F12 progenies. A total of 1254 lncRNAs were screened out, and the number of differentially expressed lncRNAs between progenies and O. sativa were significantly less than that between progenies and O. longistaminata. Some lncRNAs regulated more than one mRNA, and 89.5% of lncRNAs regulated the expression of target genes through cis-acting. A total of 78 lncRNAs and 271 mRNAs were targeted by 280 miRNAs, and 22 lncRNAs were predicted to be the precursor of 20 microRNAs. Some miRNAs were found to target their own potential precursor lncRNAs. Over 50% of lncRNAs showed parental expression level dominance (ELD) in all three progenies, and most lncRNAs showed ELD-O. sativa rather than ELD-O. longistaminata. Further analysis showed that lncRNAs might regulate the expression of plant hormone-related genes and the adaptability of O. sativa, O. longistaminata and their progenies.ConclusionsTaken together, the above results provided valuable clues for elucidating the functional features and expression differences of lncRNAs between O. sativa, O. longistaminata and their BIL progenies, and expanded our understanding about the biological functions of lncRNAs in rice.

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

Global identification and analysis revealed differentially expressed lncRNAs associated with meiosis and low fertility in autotetraploid rice

BackgroundAutotetraploid rice is a useful germplasm for polyploid rice breeding. Our previous research showed that non-coding RNAs might be associated with low fertility in autotetraploid rice. However, little information is available on long non-coding RNAs (lncRNAs) involved in the low fertility of autotetraploid rice. In the present study, RNA-seq was employed to detect the differentially expressed meiosis-related lncRNAs in autotetraploid rice, and gene overexpression and knock out experiments were used to validate the potential function of candidate lncRNA.ResultsA total of 444 differentially expressed lncRNAs (DEL) were detected during anther and ovary meiosis in autotetraploid rice. Of these, 328 DEL were associated with the transposable elements, which displayed low expression levels during meiosis in autotetraploid rice. We used rapid amplification of cDNA ends (RACE) assay to validate 10 DEL and found that the lncRNAs were not assembly artifacts, and six of them were conserved in tetraploid rice. Moreover, 237 and 20 lncRNAs were associated with pollen mother cell (PMC) and embryo sac mother cell (EMC) meiosis in autotetraploid rice, respectively. The differential expressions of some meiosis-related targets and its DEL regulator, including MEL1 regulated by TCONS_00068868, LOC_Os12g41350 (meiotic asynaptic mutant 1) by TCONS_00057811 in PMC, and LOC_Os12g39420 by TCONS_00144592 in EMC, were confirmed by qRT-PCR. TCONS_00057811, TCONS_00055980 and TCONS_00130461 showed anther specific expression patterns and were found to be highly expressed during meiosis. CRISPR/Cas9 editing of lncRNA57811 displayed similar morphology compared to wild type. The overexpression of lncRNA57811 resulted in low pollen fertility (29.70%) and seed setting (33%) in rice.ConclusionThe differential expression levels of lncRNAs, associated with transposable elements and meiosis-regulated targets, might be endogenous noncoding regulators of pollen/embryo sac development that cause low fertility in autotetraploid rice. The results enhance our understanding about rice lncRNAs, and facilitate functional research in autotetraploid rice.

Transcriptional landscape of pathogen-responsive lncRNAs in rice unveils the role of ALEX1 in jasmonate pathway and disease resistance.

SummaryPlant defence is multilayered and is essential for surviving in a changing environment. The discovery of long noncoding RNAs (lncRNAs) has dramatically extended our understanding of post‐transcriptional gene regulation in diverse biological processes. However, the expression profile and function of lncRNAs in disease resistance are still largely unknown, especially in monocots. Here, we performed strand‐specific RNA sequencing of rice leaves infected by Xanthomonas oryzae pv. Oryzae (Xoo) in different time courses and systematically identified 567 disease‐responsive rice lncRNAs. Target analyses of these lncRNAs showed that jasmonate (JA) pathway was significantly enriched. To reveal the interaction between lncRNAs and JA‐related genes, we studied the coexpression of them and found 39 JA‐related protein‐coding genes to be interplayed with 73 lncRNAs, highlighting the potential modulation of lncRNAs in JA pathway. We subsequently identified an lncRNA, ALEX1, whose expression is highly induced by Xoo infection. A T‐DNA insertion line constructed using enhancer trap system showed a higher expression of ALEX1 and exerted a significant resistance to rice bacterial blight. Functional study revealed that JA signalling is activated and the endogenous content of JA and JA‐Ile is increased. Overexpressing ALEX1 in rice further confirmed the activation of JA pathway and resistance to bacterial blight. Our findings reveal the expression of pathogen‐responsive lncRNAs in rice and provide novel insights into the connection between lncRNAs and JA pathway in the regulation of plant disease resistance.

Stress-responsive regulation of long non-coding RNA polyadenylation in Oryza sativa.

Recently, long non-coding RNAs (lncRNAs) have been demonstrated to be involved in many biological processes of plants; however, a systematic study on transcriptional and, in particular, post-transcriptional regulation of stress-responsive lncRNAs in Oryza sativa (rice) is lacking. We sequenced three types of RNA libraries (poly(A)+, poly(A)- and nuclear RNAs) under four abiotic stresses (cold, heat, drought and salt). Based on an integrative bioinformatics approach and ~200 high-throughput data sets, ~170 of which have been published, we revealed over 7000 lncRNAs, nearly half of which were identified for the first time. Notably, we found that the majority of the ~500 poly(A) lncRNAs that were differentially expressed under stress were significantly downregulated, but approximately 25% were found to have upregulated non-poly(A) forms. Moreover, hundreds of lncRNAs with downregulated polyadenylation (DPA) tend to be highly conserved, show significant nuclear retention and are co-expressed with protein-coding genes that function under stress. Remarkably, these DPA lncRNAs are significantly enriched in quantitative trait loci (QTLs) for stress tolerance or development, suggesting their potential important roles in rice growth under various stresses. In particular, we observed substantially accumulated DPA lncRNAs in plants exposed to drought and salt, which is consistent with the severe reduction of RNA 3'-end processing factors under these conditions. Taken together, the results of this study reveal that polyadenylation and subcellular localization of many rice lncRNAs are likely to be regulated at the post-transcriptional level. Our findings strongly suggest that many upregulated/downregulated lncRNAs previously identified by traditional RNA-seq analyses need to be carefully reviewed to assess the influence of post-transcriptional modification.

PLncPRO for prediction of long non-coding RNAs (lncRNAs) in plants and its application for discovery of abiotic stress-responsive lncRNAs in rice and chickpea.

Long non-coding RNAs (lncRNAs) make up a significant portion of non-coding RNAs and are involved in a variety of biological processes. Accurate identification/annotation of lncRNAs is the primary step for gaining deeper insights into their functions. In this study, we report a novel tool, PLncPRO, for prediction of lncRNAs in plants using transcriptome data. PLncPRO is based on machine learning and uses random forest algorithm to classify coding and long non-coding transcripts. PLncPRO has better prediction accuracy as compared to other existing tools and is particularly well-suited for plants. We developed consensus models for dicots and monocots to facilitate prediction of lncRNAs in non-model/orphan plants. The performance of PLncPRO was quite better with vertebrate transcriptome data as well. Using PLncPRO, we discovered 3714 and 3457 high-confidence lncRNAs in rice and chickpea, respectively, under drought or salinity stress conditions. We investigated different characteristics and differential expression under drought/salinity stress conditions, and validated lncRNAs via RT-qPCR. Overall, we developed a new tool for the prediction of lncRNAs in plants and showed its utility via identification of lncRNAs in rice and chickpea.