Genome-wide identification, putative functionality and interactions between lncRNAs and miRNAs in Brassica species

Transgressive and parental dominant gene expression and cytosine methylation during seed development in Brassica napus hybrids

KEY MESSAGE

Transcriptomic and epigenomic profiling of gene expression and small RNAs during seed and seedling development reveals expression and methylation dominance levels with implications on early stage heterosis in oilseed rape. The enhanced performance of hybrids through heterosis remains a key aspect in plant breeding; however, the underlying mechanisms are still not fully elucidated. To investigate the potential role of transcriptomic and epigenomic patterns in early expression of hybrid vigor, we investigated gene expression, small RNA abundance and genome-wide methylation in hybrids from two distant Brassica napus ecotypes during seed and seedling developmental stages using next-generation sequencing. A total of 31117, 344, 36229 and 7399 differentially expressed genes, microRNAs, small interfering RNAs and differentially methylated regions were identified, respectively. Approximately 70% of the differentially expressed or methylated features displayed parental dominance levels where the hybrid followed the same patterns as the parents. Via gene ontology enrichment and microRNA-target association analyses during seed development, we found copies of reproductive, developmental and meiotic genes with transgressive and paternal dominance patterns. Interestingly, maternal dominance was more prominent in hypermethylated and downregulated features during seed formation, contrasting to the general maternal gamete demethylation reported during gametogenesis in angiosperms. Associations between methylation and gene expression allowed identification of putative epialleles with diverse pivotal biological functions during seed formation. Furthermore, most differentially methylated regions, differentially expressed siRNAs and transposable elements were in regions that flanked genes without differential expression. This suggests that differential expression and methylation of epigenomic features may help maintain expression of pivotal genes in a hybrid context. Differential expression and methylation patterns during seed formation in an F₁ hybrid provide novel insights into genes and mechanisms with potential roles in early heterosis.

Gene expression profiling reveals transcription factor networks and subgenome bias during Brassica napus seed development

We profiled the global gene expression landscape across the reproductive lifecycle of Brassica napus. Comparative analysis of this nascent amphidiploid revealed the contribution of each subgenome to plant reproduction. Whole-genome transcription factor networks identified BZIP11 as a transcriptional regulator of early B. napus seed development. Knockdown of BZIP11 using RNA interference resulted in a similar reduction in gene activity of predicted gene targets, and a reproductive-lethal phenotype. Global mRNA profiling revealed lower accumulation of Cⁿ subgenome transcripts relative to the Aⁿ subgenome. Subgenome-specific transcription factor networks identified distinct transcription factor families enriched in each of the Aⁿ and Cⁿ subgenomes early in seed development. Analysis of laser-microdissected seed subregions further reveal subgenome expression dynamics in the embryo, endosperm and seed coat of early stage seeds. Transcription factors predicted to be regulators encoded by the Aⁿ subgenome are expressed primarily in the seed coat, whereas regulators encoded by the Cⁿ subgenome were expressed primarily in the embryo. Data suggest subgenome bias are characteristic features of the B. napus seed throughout development, and that such bias might not be universal across the embryo, endosperm and seed coat of the developing seed. Transcriptional networks spanning both the Aⁿ and Cⁿ genomes of the whole B. napus seed can identify valuable targets for seed development research and that -omics level approaches to studying gene regulation in B. napus can benefit from both broad and high-resolution analyses.

Genome-wide prediction of cauliflower miRNAs and lncRNAs and their roles in post-transcriptional gene regulation

We have predicted miRNAs, their targets and lncRNAs from the genome of Brassica oleracea along with their functional annotation. Selected miRNAs and their targets are experimentally validated. Roles of these non-coding RNAs in post-transcriptional gene regulation are also deciphered. Cauliflower (Brassica oleracea var. Botrytis) is an important vegetable crop for its dietary and medicinal values with rich source of vitamins, dietary fibers, flavonoids and antioxidants. MicroRNAs (miRNAs) are small non-coding RNAs (ncRNAs), which regulate gene expression by inhibiting translation or by degrading messenger RNAs (mRNAs). On the other hand, long non-coding RNAs (lncRNAs) are responsible for the up regulation and the down regulation of transcription. Although the genome of cauliflower is reported, yet the roles of these ncRNAs in post-transcriptional gene regulation (PTGR) remain elusive. In this study, we have computationally predicted 355 miRNAs, of which 280 miRNAs are novel compared to miRBase 22.1. All the predicted miRNAs belong to 121 different families. We have also identified 934 targets of 125 miRNAs along with their functional annotation. These targets are further classified into biological processes, molecular functions and cellular components. Moreover, we have predicted 634 lncRNAs, of which 61 are targeted by 30 novel miRNAs. Randomly chosen 10 miRNAs and 10 lncRNAs are experimentally validated. Five miRNA targets including squamosa promoter-binding-like protein 9, homeobox-leucine zipper protein HDG12-like, NAC domain-containing protein 100, CUP-SHAPED COTYLEDON 1 and kinesin-like protein NACK2 of four miRNAs including bol-miR156a, bol-miR162a, bol-miR164d and bol-miR2673 are also experimentally validated. We have built network models of interactions between miRNAs and their target mRNAs, as well as between miRNAs and lncRNAs. Our findings enhance the knowledge of non-coding genome of cauliflower and their roles in PTGR, and might play important roles in improving agronomic traits of this economically important crop.

Novel lncRNA lncRNA001074 participates in the low salinity-induced response in the sea cucumber Apostichopus japonicus by targeting the let-7/NKAα axis

Salinity fluctuations have severe impacts on sea cucumbers and therefore important consequences in sea cucumber farming. The responses of sea cucumbers to salinity changes are reflected in the expression profiles of multiple genes and non-coding RNAs (ncRNAs). The microRNA (let-7) which is a developmental regulator, the ion transporter gene sodium potassium ATPase gene (NKAα), and the long ncRNA lncRNA001074 were previously shown to be involved in responses to salinity changes in various marine species. To better understand the relationship between ncRNAs and target genes, the let-7/NKAα/lncRNA001074 predicted interaction was investigated in this study using luciferase reporter assays and gene knockdowns in the sea cucumber Apostichopus japonicus. The results showed that NKAα was the target gene of let-7 and NKAα expression levels were inversely correlated with let-7 expression based on the luciferase reporter assays and western blots. The let-7 abundance was negatively regulated by lncRNA001074 and NKAα both in vitro and in vivo. Knockdown of lncRNA001074 led to let-7 overexpression. These results demonstrated that lncRNA001074 binds to the 3'-UTR binding site of let-7 in a regulatory manner. Furthermore, the expression profiles of let-7, NKAα, and lncRNA001074 were analyzed in sea cucumbers after the knockdown of each of these genes. The results found that lncRNA001074 competitively bound let-7 to suppress NKAα expression under low salinity conditions. The downregulation of let-7, in conjunction with the upregulation of lncRNA001074 and NKAα, may be essential for the response to low salinity change in sea cucumbers. Therefore, the dynamic balance of the lncRNA001074, NKAα, and let-7 network might be a potential response mechanism to salinity change in sea cucumbers.

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

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

Comparative analysis of long noncoding RNAs in angiosperms and characterization of long noncoding RNAs in response to heat stress in Chinese cabbage

Long noncoding RNAs (lncRNAs) are widely present in different species and play critical roles in response to abiotic stresses. However, the functions of lncRNAs in Chinese cabbage under heat stress remain unknown. Here, we first conducted a global comparative analysis of 247,242 lncRNAs among 37 species. The results indicated that lncRNAs were poorly conserved among different species, and only 960 lncRNAs were homologous to 524 miRNA precursors. We then carried out lncRNA sequencing for a genome-wide analysis of lncRNAs and their target genes in Chinese cabbage at different stages of heat treatment. In total, 18,253 lncRNAs were identified, of which 1229 differentially expressed (DE) lncRNAs were characterized as being heat-responsive. The ceRNA network revealed that 38 lncRNAs, 16 miRNAs, and 167 mRNAs were involved in the heat response in Chinese cabbage. Combined analysis of the cis- and trans-regulated genes indicated that the targets of DE lncRNAs were significantly enriched in the "protein processing in endoplasmic reticulum" and "plant hormone signal transduction" pathways. Furthermore, the majority of HSP and PYL genes involved in these two pathways exhibited similar expression patterns and responded to heat stress rapidly. Based on the networks of DE lncRNA-mRNAs, 29 and 22 lncRNAs were found to interact with HSP and PYL genes, respectively. Finally, the expression of several critical lncRNAs and their targets was verified by qRT-PCR. Overall, we conducted a comparative analysis of lncRNAs among 37 species and performed a comprehensive analysis of lncRNAs in Chinese cabbage. Our findings expand the knowledge of lncRNAs involved in the heat stress response in Chinese cabbage, and the identified lncRNAs provide an abundance of resources for future comparative and functional studies.

Genome-wide identification of Aedes albopictus long noncoding RNAs and their association with dengue and Zika virus infection

The Asian tiger mosquito, Aedes albopictus (Ae. albopictus), is an important vector that transmits arboviruses such as dengue (DENV), Zika (ZIKV) and Chikungunya virus (CHIKV). Long noncoding RNAs (lncRNAs) are known to regulate various biological processes. Knowledge on Ae. albopictus lncRNAs and their functional role in virus-host interactions are still limited. Here, we identified and characterized the lncRNAs in the genome of an arbovirus vector, Ae. albopictus, and evaluated their potential involvement in DENV and ZIKV infection. We used 148 public datasets, and identified a total of 10, 867 novel lncRNA transcripts, of which 5,809, 4,139, and 919 were intergenic, intronic and antisense respectively. The Ae. albopictus lncRNAs shared many characteristics with other species such as short length, low GC content, and low sequence conservation. RNA-sequencing of Ae. albopictus cells infected with DENV and ZIKV showed that the expression of lncRNAs was altered upon virus infection. Target prediction analysis revealed that Ae. albopictus lncRNAs may regulate the expression of genes involved in immunity and other metabolic and cellular processes. To verify the role of lncRNAs in virus infection, we generated mutations in lncRNA loci using CRISPR-Cas9, and discovered that two lncRNA loci mutations, namely XLOC_029733 (novel lncRNA transcript id: lncRNA_27639.2) and LOC115270134 (known lncRNA transcript id: XR_003899061.1) resulted in enhancement of DENV and ZIKV replication. The results presented here provide an important foundation for future studies of lncRNAs and their relationship with virus infection in Ae. albopictus.

Ae. albopictus is an important vector of arboviruses such as dengue and Zika viruses. Studies on virus-host interaction at gene expression and molecular level are crucial especially in devising methods to inhibit virus replication in Aedes mosquitoes. Previous reports have shown that, besides protein-coding genes, noncoding RNAs such as lncRNAs are also involved in virus-host interaction. In this study, we report a comprehensive catalog of novel lncRNA transcripts in the genome of Ae. albopictus. We also show that the expression of lncRNAs was altered upon infection with dengue and Zika. Additionally, depletion of certain lncRNAs resulted in increased replication of dengue and Zika; hence, suggesting potential association of lncRNAs in virus infection. Results of this study provide a new avenue to the investigation of mosquito-virus interactions, especially in the aspect of noncoding genes.

Long Non-Coding RNAs, the Dark Matter: An Emerging Regulatory Component in Plants

Long non-coding RNAs (lncRNAs) are pervasive transcripts of longer than 200 nucleotides and indiscernible coding potential. lncRNAs are implicated as key regulatory molecules in various fundamental biological processes at transcriptional, post-transcriptional, and epigenetic levels. Advances in computational and experimental approaches have identified numerous lncRNAs in plants. lncRNAs have been found to act as prime mediators in plant growth, development, and tolerance to stresses. This review summarizes the current research status of lncRNAs in planta, their classification based on genomic context, their mechanism of action, and specific bioinformatics tools and resources for their identification and characterization. Our overarching goal is to summarize recent progress on understanding the regulatory role of lncRNAs in plant developmental processes such as flowering time, reproductive growth, and abiotic stresses. We also review the role of lncRNA in nutrient stress and the ability to improve biotic stress tolerance in plants. Given the pivotal role of lncRNAs in various biological processes, their functional characterization in agriculturally essential crop plants is crucial for bridging the gap between phenotype and genotype.

Triple-Negative Breast Cancer Comparison With Canine Mammary Tumors From Light Microscopy to Molecular Pathology

Many similar characteristics in human and dog cancers including, spontaneous development, clinical presentation, tumor heterogeneity, disease progression, and response to standard therapies have promoted the approval of this comparative model as an alternative to mice. Breast cancer represents the second most frequent neoplasm in humans after lung cancer. Triple-negative breast cancers (TNBC) constitute around 15% of all cases of breast cancer and do not express estrogen receptor (ER), progesterone receptor (PR), and do not overexpress human epidermal growth factor receptor 2 (HER2). As a result, they do not benefit from hormonal or trastuzumab-based therapy. Patients with TNBC have worse overall survival than patients with non-TNBC. Lehmann and collaborators described six different molecular subtypes of TNBC which further demonstrated its transcriptional heterogeneity. This six TNBC subtype classification has therapeutic implications. Breast cancer is the second most frequent neoplasm in sexually intact female dogs after skin cancer. Canine mammary tumors are a naturally occurring heterogeneous group of cancers that have several features in common with human breast cancer (HBC). These similarities include etiology, signaling pathway activation, and histological classification. Molecularly CMTs are more like TNBCs, and therefore dogs are powerful spontaneous models of cancer to test new therapeutic approaches, particularly for human TNBCs. More malignant tumors of the breast are more often ER and PR negative in both humans and dogs. Promising breast cancer biomarkers in both humans and canines are cancer-associated stroma (CAS), circulating tumor cells and tumor DNA (ctDNA), exosomes and miRNAs, and metabolites.

Identification of lncRNAs in response to infection by Plasmodiophora brassicae in Brassica napus and development of lncRNA-based SSR markers

Clubroot resistance in spring canola has been introgressed from different Brassica sources; however, molecular mechanism underlying this resistance, especially the involvement of long non-coding RNAs (lncRNAs), is yet to be understood. We identified 464 differentially expressed (DE) lncRNAs from the roots of clubroot-resistant canola, carrying resistance on chromosome BnaA03, and susceptible canola lines challenged with Plasmodiophora brassicae pathotype 3. Pathway enrichment analysis showed that most of the target genes regulated by these DE lncRNAs belonged to plant-pathogen interaction and hormone signaling, as well as primary and secondary metabolic pathways. Comparative analysis of these lncRNAs with 530 previously reported DE lncRNAs, identified using resistance located on BnaA08, detected 12 lncRNAs that showed a similar trend of upregulation in both types of resistant lines; these lncRNAs probably play a fundamental role in clubroot resistance. We identified SSR markers within 196 DE lncRNAs. Genotyping of two DH populations carrying resistance on BnaA03 identified a marker capable of detecting the resistance in 98% of the DH lines. To our knowledge, this is the first report of the identification of SSRs within lncRNAs responsive to P. brassicae infection, demonstrating the potential use of lncRNAs in the breeding of Brassica crops.

Characterizing the Role of the miR156-SPL Network in Plant Development and Stress Response

MicroRNA (miRNA) is a short, single-stranded, non-coding RNA found in eukaryotic cells that can regulate the expression of many genes at the post-transcriptional level. Among various plant miRNAs with diverse functions, miR156 plays a key role in biological processes, including developmental regulation, immune response, metabolic regulation, and abiotic stress. MiRNAs have become the regulatory center for plant growth and development. MicroRNA156 (miR156) is a highly conserved and emerging tool for the improvement of plant traits, including crop productivity and stress tolerance. Fine-tuning of squamosa promoter biding-like (SPL) gene expression might be a useful strategy for crop improvement. Here, we studied the regulation of the miR156 module and its interaction with SPL factors to understand the developmental transition of various plant species. Furthermore, this review provides a strong background for plant biotechnology and is an important source of information for further molecular breeding to optimize farming productivity.

Splicing conservation signals in plant long noncoding RNAs

Long noncoding RNAs (lncRNAs) have recently emerged as prominent regulators of gene expression in eukaryotes. LncRNAs often drive the modification and maintenance of gene activation or gene silencing states via chromatin conformation rearrangements. In plants, lncRNAs have been shown to participate in gene regulation, and are essential to processes such as vernalization and photomorphogenesis. Despite their prominent functions, only over a dozen lncRNAs have been experimentally and functionally characterized. Similar to its animal counterparts, the rates of sequence divergence are much higher in plant lncRNAs than in protein coding mRNAs, making it difficult to identify lncRNA conservation using traditional sequence comparison methods. Beyond this, little is known about the evolutionary patterns of lncRNAs in plants. Here, we characterized the splicing conservation of lncRNAs in Brassicaceae. We generated a whole-genome alignment of 16 Brassica species and used it to identify synthenic lncRNA orthologs. Using a scoring system trained on transcriptomes from A. thaliana and B. oleracea, we identified splice sites across the whole alignment and measured their conservation. Our analysis revealed that 17.9% (112/627) of all intergenic lncRNAs display splicing conservation in at least one exon, an estimate that is substantially higher than previous estimates of lncRNA conservation in this group. Our findings agree with similar studies in vertebrates, demonstrating that splicing conservation can be evidence of stabilizing selection. We provide conclusive evidence for the existence of evolutionary deeply conserved lncRNAs in plants and describe a generally applicable computational workflow to identify functional lncRNAs in plants.

Comparative Transcriptome Analysis Reveals New lncRNAs Responding to Salt Stress in Sweet Sorghum

Long non-coding RNAs (lncRNAs) can enhance plant stress resistance by regulating the expression of functional genes. Sweet sorghum is a salt-tolerant energy crop. However, little is known about how lncRNAs in sweet sorghum respond to salt stress. In this study, we identified 126 and 133 differentially expressed lncRNAs in the salt-tolerant M-81E and the salt-sensitive Roma strains, respectively. Salt stress induced three new lncRNAs in M-81E and inhibited two new lncRNAs in Roma. These lncRNAs included lncRNA13472, lncRNA11310, lncRNA2846, lncRNA26929, and lncRNA14798, which potentially function as competitive endogenous RNAs (ceRNAs) that influence plant responses to salt stress by regulating the expression of target genes related to ion transport, protein modification, transcriptional regulation, and material synthesis and transport. Additionally, M-81E had a more complex ceRNA network than Roma. This study provides new information regarding lncRNAs and the complex regulatory network underlying salt-stress responses in sweet sorghum.

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.

Genome-wide analysis of long non-coding RNAs (lncRNAs) in two contrasting rapeseed (Brassica napus L.) genotypes subjected to drought stress and re-watering

BACKGROUND

Drought stress is a major abiotic factor that affects rapeseed (Brassica napus L.) productivity. Though previous studies indicated that long non-coding RNAs (lncRNAs) play a key role in response to drought stress, a scheme for genome-wide identification and characterization of lncRNAs' response to drought stress is still lacking, especially in the case of B. napus. In order to further understand the molecular mechanism of the response of B. napus to drought stress, we compared changes in the transcriptome between Q2 (a drought-tolerant genotype) and Qinyou8 (a drought-sensitive genotype) responding drought stress and rehydration treatment at the seedling stage.

RESULTS

A total of 5546 down-regulated and 6997 up-regulated mRNAs were detected in Q2 compared with 7824 and 10,251 in Qinyou8, respectively; 369 down-regulated and 108 up- regulated lncRNAs were detected in Q2 compared with 449 and 257 in Qinyou8, respectively. LncRNA-mRNA interaction network analysis indicated that the co-expression network of Q2 was composed of 145 network nodes and 5175 connections, while the co-expression network of Qinyou8 was composed of 305 network nodes and 22,327 connections. We further identified 34 transcription factors (TFs) corresponding to 126 differentially expressed lncRNAs in Q2, and 45 TFs corresponding to 359 differentially expressed lncRNAs in Qinyou8. Differential expression analysis of lncRNAs indicated that up- and down-regulated mRNAs co-expressed with lncRNAs participated in different metabolic pathways and were involved in different regulatory mechanisms in the two genotypes. Notably, some lncRNAs were co-expressed with BnaC07g44670D, which are associated with plant hormone signal transduction. Additionally, some mRNAs co-located with XLOC_052298, XLOC_094954 and XLOC_012868 were mainly categorized as signal transport and defense/stress response.

CONCLUSIONS

The results of this study increased our understanding of expression characterization of rapeseed lncRNAs in response to drought stress and re-watering, which would be useful to provide a reference for the further study of the function and action mechanisms of lncRNAs under drought stress and re-watering.

The green peach aphid gut contains host plant microRNAs identified by comprehensive annotation of Brassica oleracea small RNA data

Like all organisms, aphids, plant sap-sucking insects that house a bacterial endosymbiont called Buchnera, are members of a species interaction network. Ecological interactions across such networks can result in phenotypic change in network members mediated by molecular signals, like microRNAs. Here, we interrogated small RNA data from the aphid, Myzus persicae, to determine the source of reads that did not map to the aphid or Buchnera genomes. Our analysis revealed that the pattern was largely explained by reads that mapped to the host plant, Brassica oleracea, and a facultative symbiont, Regiella. To start elucidating the function of plant small RNA in aphid gut, we annotated 213 unique B. oleracea miRNAs; 32/213 were present in aphid gut as mature and star miRNAs. Next, we predicted targets in the B. oleracea and M. persicae genomes for these 32 plant miRNAs. We found that plant targets were enriched for genes associated with transcription, while the distribution of targets in the aphid genome was similar to the functional distribution of all genes in the aphid genome. We discuss the potential of plant miRNAs to regulate aphid gene expression and the mechanisms involved in processing, export and uptake of plant miRNAs by aphids.

Identification of conserved and novel miRNAs responsive to heat stress in flowering Chinese cabbage using high-throughput sequencing

Plant microRNAs (miRNAs) are noncoding and endogenous key regulators that play significant functions in regulating plant responses to stress, and plant growth and development. Heat stress is a critical abiotic stress that reduces the yield and quality of flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee). However, limited information is available on whether miRNAs are involved in the regulation of heat stress in B. campestris. A high-throughput sequencing approach was used to identify novel and conserved heat-responsive miRNAs in four small RNA libraries of flowering Chinese cabbage using leaves collected at 0 h, 1 h, 6 h and 12 h after a 38 °C heat-stress treatment. The analysis identified 41 conserved miRNAs (belonging to 19 MIR families), of which MIR156, MIR159, MIR168, MIR171 and MIR1885 had the most abundant molecules. Prediction and evaluation of novel miRNAs using the unannotated reads resulted in 18 candidate miRNAs. Differential expression analysis showed that most of the identified miRNAs were downregulated in heat-treated groups. To better understand functional importance, bioinformatic analysis predicted 432 unique putative target miRNAs involved in cells, cell parts, catalytic activity, cellular processes and abiotic stress responses. Furthermore, the Kyoto Encyclopedia of Genes and Genomes maps of flowering Chinese cabbage identified the significant role of miRNAs in stress adaptation and stress tolerance, and in several mitogen-activated protein kinases signaling pathways including cell death. This work presents a comprehensive study of the miRNAs for understanding the regulatory mechanisms and their participation in the heat stress of flowering Chinese cabbage.

Non-coding RNAs: Functional roles in the regulation of stress response in Brassica crops

Brassica crops face a combination of different abiotic and biotic stresses in the field that can reduce plant growth and development by affecting biochemical and morpho-physiological processes. Emerging evidence suggests that non-coding RNAs (ncRNAs), especially microRNAs (miRNAs) and long ncRNAs (lncRNAs), play a significant role in the modulation of gene expression in response to plant stresses. Recent advances in computational and experimental approaches are of great interest for identifying and functionally characterizing ncRNAs. While progress in this field is limited, numerous ncRNAs involved in the regulation of gene expression in response to stress have been reported in Brassica. In this review, we summarize the modes of action and functions of stress-related miRNAs and lncRNAs in Brassica as well as the approaches used to identify ncRNAs.

Genome-Wide Analysis of Coding and Long Non-Coding RNAs Involved in Cuticular Wax Biosynthesis in Cabbage (Brassica oleracea L. var. capitata)

Cuticular wax is a mixture of very long chain fatty acids (VLCFAs) and their derivatives, which determines vital roles for plant growth. In cabbage, the cuticular wax content of leaf blades is an important trait influencing morphological features of the head. Understanding the molecular basis of cuticular wax biosynthesis can help breeders develop high quality cabbage varieties. Here, we characterize a cabbage non-wax glossy (nwgl) plant, which exhibits glossy green phenotype. Cryo-scanning electron microscope analysis showed abnormal wax crystals on the leaf surfaces of nwgl plants. Cuticular wax composition analyzed by GC-MS displayed severely decreased in total wax loads, and individual wax components in nwgl leaves. We delimited the NWGL locus into a 99-kb interval between the at004 marker and the end of chromosome C08 through fine mapping. By high-throughput RNA sequencing, we identified 1247 differentially expressed genes (DEGs) and 148 differentially expressed lncRNAs in nwgl leaves relative to the wild-type. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis revealed that the DEGs and cis-regulated target genes for differentially expressed lncRNAs were significantly enriched in wax and lipid biosynthetic or metabolic processes. Our results provide the novel foundation to explore the complex molecular basis of cuticular wax biosynthesis.

Genome-wide analysis of long non-coding RNAs unveils the regulatory roles in the heat tolerance of Chinese cabbage (Brassica rapa ssp.chinensis)

Long non-coding RNAs (lncRNAs) mediate important epigenetic regulation in various biological processes related to the stress response in plants. However, the systematic analysis of the lncRNAs expressed in Brassica rapa under heat stress has been elusive. In this study, we performed a genome-wide analysis of the lncRNA expression profiles in non-heading Chinese cabbage leaves using strand-specific RNA-sequencing. A total of 4594 putative lncRNAs were identified with a comprehensive landscape of dynamic lncRNA expression networks under heat stress. Co-expression networks of the interactions among the differentially expressed lncRNAs, mRNAs and microRNAs revealed that several phytohormones were associated with heat tolerance, including salicylic acid (SA) and brassinosteroid (BR) pathways. Of particular importance is the discovery of 25 lncRNAs that were highly co-expressed with 10 heat responsive genes. Thirty-nine lncRNAs were predicted as endogenous target mimics (eTMs) for 35 miRNAs, and five of them were validated to be involved in the heat tolerance of Chinese cabbage. Heat responsive lncRNA (TCONS_00048391) is an eTM for bra-miR164a, that could be a sponge for miRNA binding and may be a competing endogenous RNA (ceRNA) for the target gene NAC1 (Bra030820), affecting the expression of bra-miR164a in Chinese cabbage. Thus, these findings provide new insights into the functions of lncRNAs in heat tolerance and highlight a set of candidate lncRNAs for further studies in non-heading Chinese cabbage.

Potential regulatory network in the PSG10P/miR-19a-3p/IL1RAP pathway is possibly involved in preeclampsia pathogenesis

Preeclampsia (PE), a pregnancy-specific disorder, is a leading cause of perinatal maternal-fetal mortality and morbidity. Impaired cell migration and invasion of trophoblastic cells and an imbalanced systemic maternal inflammatory response have been proposed as potential mechanisms of PE pathogenesis. Comparative analysis between PE placentas and normal placentas profiled differentially expressed miRNAs, lncRNAs, and mRNAs, including miR-19a-3p (miRNA), PSG10P (lncRNA), and IL1RAP (mRNA). This study was conducted to investigate their potential roles in PE pathogenesis. The expression of miR-19a-3p, PSG10P, and IL1RAP was examined in PE and normal placentas using RT-qPCR. An in vitro experiment was performed in human trophoblast HET8/SVneo and TEV-1 cells cultured in normoxic and hypoxic conditions. MiR-19a-3p targets were identified using Targetscan, miRanda, and PicTar analysis as well as luciferase reporter assays. The mouse model of PE was conducted using sFlt-1 for in vivo tests. Lower levels of miR-19a-3p, but higher levels of PSG10P and IL1RAP were observed in PE placentas and the trophoblast cells in hypoxia. Luciferase reporter assays confirmed that PSG10P and IL1RAP were both direct targets of miR-19a-3p. Exposure to hypoxia inhibited cell viability, migration, and invasion of HET8/SVneo and TEV-1 cells. Knocking out PSG10P and IL1RAP or overexpressing miR-19a-3p rescued the inhibition caused by hypoxia. In vivo experiments showed that IL1RAP promoted the expression of caspase-3, a key apoptosis enzyme, but inhibited MMP9, which is responsible for degrading the extracellular matrix, suggesting a significant role of IL1RAP in cell proliferation, migration, and invasion. miR-19a-3p, PSG10P, and IL1RAP were all found to be involved in PE pathogenesis. With a common targeting region in their sequences, a regulatory network in the PSG10P/miR-19a-3p/IL1RAP pathway may contribute to PE pathogenesis during pregnancy.