Integrative analysis of genome‐wide lncRNA and mRNA expression in newly synthesized Brassica hexaploids

Polyploidization, as a significant evolution force, has been considered to facilitate plant diversity. The expression levels of lncRNAs and how they control the expression of protein‐coding genes in allopolyploids remain largely unknown. In this study, lncRNA expression profiles were compared between Brassica hexaploid and its parents using a high‐throughput sequencing approach. A total of 2,725, 1,672, and 2,810 lncRNAs were discovered in Brassica rapa, Brassica carinata, and Brassica hexaploid, respectively. It was also discovered that 725 lncRNAs were differentially expressed between Brassica hexaploid and its parents, and 379 lncRNAs were nonadditively expressed in this hexaploid. LncRNAs have multiple expression patterns between Brassica hexaploid and its parents and show paternal parent‐biased expression. These lncRNAs were found to implement regulatory functions directly in the long‐chain form, and acted as precursors or targets of miRNAs. According to the prediction of the targets of differentially expressed lncRNAs, 109 lncRNAs were annotated, and their target genes were involved in the metabolic process, pigmentation, reproduction, exposure to stimulus, biological regulation, and so on. Compared with the paternal parent, differentially expressed lncRNAs between Brassica hexaploid and its maternal parent participated in more regulation pathways. Additionally, 61 lncRNAs were identified as putative targets of known miRNAs, and 15 other lncRNAs worked as precursors of miRNAs. Some conservative motifs of lncRNAs from different groups were detected, which indicated that these motifs could be responsible for their regulatory roles. Our findings may provide a reference for the further study of the function and action mechanisms of lncRNAs during plant evolution.     

Sl‐lncRNA15492 interacts with Sl‐miR482a and affects Solanum lycopersicum immunity against Phytophthora infestans

Long non‐coding RNAs (lncRNAs) are involved in the resistance of plants to infection by pathogens via interactions with microRNAs (miRNAs). Long non‐coding RNAs are cleaved by miRNAs to produce phased small interfering RNAs (phasiRNAs), which, as competing endogenous RNAs (ceRNAs), function as decoys for mature miRNAs, thus inhibiting their expression, and contain pre‐miRNA sequences to produce mature miRNAs. However, whether lncRNAs and miRNAs mediate other molecular mechanisms during plant resistance to pathogens is unknown. In this study, as a positive regulator, Sl‐lncRNA15492 from tomato (Solanum lycopersicum Zaofen No. 2) plants affected tomato resistance to Phytophthora infestans. Gain‐ and loss‐of‐function experiments and RNA ligase‐mediated 5′‐amplification of cDNA ends (RLM‐5′ RACE) also revealed that Sl‐miR482a was negatively involved in tomato resistance by targeting Sl‐NBS‐LRR genes and that silencing of Sl‐NBS‐LRR1 decreased tomato resistance. Sl‐lncRNA15492 inhibited the expression of mature Sl‐miR482a, whose precursor was located within the antisense sequence of Sl‐lncRNA15492. Further degradome analysis and additional RLM‐5′ RACE experiments verified that mature Sl‐miR482a could also cleave Sl‐lncRNA15492. These results provide a mechanism by which lncRNAs might inhibit precursor miRNA expression through antisense strands of lncRNAs, and demonstrate that Sl‐lncRNA15492 and Sl‐miR482a mutually inhibit the maintenance of Sl‐NBS‐LRR1 homeostasis during tomato resistance to P. infestans.     

Thank you

The importance of long non‐coding RNAs (lncRNAs) in plant development has been established, but a systematic analysis of lncRNAs expressed during pollen development and fertilization has been elusive. We performed a time series of RNA‐seq experiments at five developmental stages during pollen development and three different time points after pollination in Brassica rapa and identified 12 051 putative lncRNAs. A comprehensive view of dynamic lncRNA expression networks underpinning pollen development and fertilization was provided. B. rapa lncRNAs share many common characteristics of lncRNAs: relatively short length, low expression but specific in narrow time windows, and low evolutionary conservation. Gene modules and key lncRNAs regulating reproductive development such as exine formation were uncovered. Forty‐seven cis‐acting lncRNAs and 451 trans‐acting lncRNAs were revealed to be highly coexpressed with their target protein‐coding genes. Of particular importance are the discoveries of 14 lncRNAs that were highly coexpressed with 10 function‐known pollen‐associated coding genes. Fifteen lncRNAs were predicted as endogenous target mimics for 13 miRNAs, and two lncRNAs were proved to be functional target mimics for miR160 after experimental verification and shown to function in pollen development. Our study provides the systematic identification of lncRNAs during pollen development and fertilization in B. rapa and forms the foundation for future genetic, genomic, and evolutionary studies.     

Genetic dissection of the gene coexpression network underlying photosynthesis in Populus

Photosynthesis is a key reaction that ultimately generates the carbohydrates needed to form woody tissues in trees. However, the genetic regulatory network of protein‐encoding genes (PEGs) and regulatory noncoding RNAs (ncRNAs), including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), underlying the photosynthetic pathway is unknown. Here, we integrated data from coexpression analysis, association studies (additive, dominance and epistasis), and expression quantitative trait nucleotide (eQTN) mapping to dissect the causal variants and genetic interaction network underlying photosynthesis in Populus. We initially used 30 PEGs, 6 miRNAs and 12 lncRNAs to construct a coexpression network based on the tissue‐specific gene expression profiles of 15 Populus samples. Then, we performed association studies using a natural population of 435 unrelated Populus tomentosa individuals, and identified 72 significant associations (P ≤ 0.001, q ≤ 0.05) with diverse additive and dominance patterns underlying photosynthesis‐related traits. Analysis of epistasis and eQTNs revealed that the complex genetic interactions in the coexpression network contribute to phenotypes at various levels. Finally, we demonstrated that heterologously expressing the most highly linked gene (PtoPsbX1) in this network significantly improved photosynthesis in Arabidopsis thaliana, pointing to the functional role of PtoPsbX1 in the photosynthetic pathway. This study provides an integrated strategy for dissecting a complex genetic interaction network, which should accelerate marker‐assisted breeding efforts to genetically improve woody plants.     

Identification of lncRNA‐associated differential subnetworks in oesophageal squamous cell carcinoma by differential co‐expression analysis

Differential expression analysis has led to the identification of important biomarkers in oesophageal squamous cell carcinoma (ESCC). Despite enormous contributions, it has not harnessed the full potential of gene expression data, such as interactions among genes. Differential co‐expression analysis has emerged as an effective tool that complements differential expression analysis to provide better insight of dysregulated mechanisms and indicate key driver genes. Here, we analysed the differential co‐expression of lncRNAs and protein‐coding genes (PCGs) between normal oesophageal tissue and ESCC tissues, and constructed a lncRNA‐PCG differential co‐expression network (DCN). DCN was characterized as a scale‐free, small‐world network with modular organization. Focusing on lncRNAs, a total of 107 differential lncRNA‐PCG subnetworks were identified from the DCN by integrating both differential expression and differential co‐expression. These differential subnetworks provide a valuable source for revealing lncRNA functions and the associated dysfunctional regulatory networks in ESCC. Their consistent discrimination suggests that they may have important roles in ESCC and could serve as robust subnetwork biomarkers. In addition, two tumour suppressor genes (AL121899.1 and ELMO2), identified in the core modules, were validated by functional experiments. The proposed method can be easily used to investigate differential subnetworks of other molecules in other cancers.