Integrative Analysis of the lncRNA and mRNA Transcriptome Revealed Genes and Pathways Potentially Involved in the Anther Abortion of Cotton (Gossypium hirsutum L.)

The roles of non-coding RNAs in male reproductive development and abiotic stress responses during this unique process in flowering plants

Successful male reproductive development is the guarantee for sexual reproduction of flowering plants. Male reproductive development is a complicated and multi-stage process that integrates physiological processes and adaptation and tolerance to a myriad of environmental stresses. This well-coordinated process is governed by genetic and epigenetic machineries. Non-coding RNAs (ncRNAs) play pleiotropic roles in the plant growth and development. The identification, characterization and functional analysis of ncRNAs and their target genes have opened a new avenue for comprehensively revealing the regulatory network of male reproductive development and its response to environmental stresses in plants. This review briefly addresses the types, origin, biogenesis and mechanisms of ncRNAs in plants, highlights important updates on the roles of ncRNAs in regulating male reproductive development and emphasizes the contribution of ncRNAs, especially miRNAs and lncRNAs, in responses to abiotic stresses during this unique process in flowering plants.

An analysis of lncRNAs related to fiber quality and the discovery of their target genes in a Gossypium hirsutum line with Gossypium mustelinum introgression

KEY MESSAGE

Analysis of fiber quality lncRNAs and their target genes from a pair of Gossypium mustelinum near-isogenic lines provide new prospects for improving the fiber quality of Upland cotton. Long noncoding RNAs (lncRNAs) are an important part of genome transcription and play roles in a wide range of biological processes in plants. In this research, a pair of near-isogenic cotton lines, namely, a Gossypium mustelinum introgression line (IL9) with outstanding fiber quality and its recurrent Upland cotton parent (PD94042), were used as the experimental materials. Cotton fibers were selected for lncRNA sequencing at 17 and 21 days post-anthesis. A total of 2693 differentially expressed genes were identified. In total, 5841 lncRNAs were ultimately screened, from which 163 differentially expressed lncRNAs were identified. Target genes of the lncRNAs were predicted by two different methods: cis and trans. Some of the target genes were related to cell components, membrane components, plant hormone signal transduction and catalytic metabolism, and the results indicated that there might also be important effects on the development of fiber. Four differentially expressed target genes related to fiber quality (Gomus.D05G015100, Gomus.A05G281300, Gomus.A12G023400 and Gomus.A10G226800) were screened through gene function annotation, and the functions of these four genes were verified through virus-induced gene silencing (VIGS). Compared to the negative controls, plants in which any of these four genes were silenced showed significant reductions in fiber strength. In addition, the plants in which the Gomus.A12G023400 gene was silenced showed a significant reduction in fiber uniformity, whereas the plants in which Gomus.A05G281300 was silenced showed a significant increase in fiber fineness as measured via micronaire. Our results showed that these genes play different roles during fiber development, impacting fiber quality.

Genomic survey of high-throughput RNA-Seq data implicates involvement of long intergenic non-coding RNAs (lincRNAs) in cytoplasmic male-sterility and fertility restoration in pigeon pea

BACKGROUND

Long-intergenic non-coding RNAs (lincRNAs) originate from intergenic regions and have no coding potential. LincRNAs have emerged as key players in the regulation of various biological processes in plant development. Cytoplasmic male-sterility (CMS) in association with restorer-of-fertility (Rf) systems makes it a highly reliable tool for exploring heterosis for producing commercial hybrid seeds. To date, there have been no reports of lincRNAs during pollen development in CMS and fertility restorer lines in pigeon pea.

OBJECTIVE

Identification of lincRNAs in the floral buds of cytoplasmic male-sterile (AKCMS11) and fertility restorer (AKPR303) pigeon pea lines.

METHODS

We employed a computational approach to identify lincRNAs in the floral buds of cytoplasmic male-sterile (AKCMS11) and fertility restorer (AKPR303) pigeon pea lines using RNA-Seq data.

RESULTS

We predicted a total of 2145 potential lincRNAs of which 966 were observed to be differentially expressed between the sterile and fertile pollen. We identified, 927 cis-regulated and 383 trans-regulated target genes of the lincRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the target genes revealed that these genes were specifically enriched in pathways like pollen and pollen tube development, oxidative phosphorylation, etc. We detected 23 lincRNAs that were co-expressed with 17 pollen-related genes with known functions. Fifty-nine lincRNAs were predicted to be endogenous target mimics (eTMs) for 25 miRNAs, and found to be associated with pollen development. The, lincRNA regulatory networks revealed that different lincRNA-miRNA-mRNA networks might be associated with CMS and fertility restoration.

CONCLUSION

Thus, this study provides valuable information by highlighting the functions of lincRNAs as regulators during pollen development in pigeon pea and utilization in hybrid seed production.

Research progress on the roles of lncRNAs in plant development and stress responses

Long non-coding RNAs (lncRNAs) are RNAs of more than 200 nucleotides in length that are not (or very rarely) translated into proteins. In eukaryotes, lncRNAs regulate gene expression at the transcriptional, post-transcriptional, and epigenetic levels. lncRNAs are categorized according to their genomic position and molecular mechanism. This review summarized the characteristics and mechanisms of plant lncRNAs involved in vegetative growth, reproduction, and stress responses. Our discussion and model provide a theoretical basis for further studies of lncRNAs in plant breeding.

Whole-transcriptome analyses identify key differentially expressed mRNAs, lncRNAs, and miRNAs associated with male sterility in watermelon

Male sterility is a valuable trait for watermelon breeding, as watermelon hybrids exhibit obvious heterosis. However, the underlying regulatory mechanism is still largely unknown, especially regarding the related non-coding genes. In the present study, approximately 1035 differentially expressed genes (DEGs), as well as 80 DE-lncRNAs and 10 DE-miRNAs, were identified, with the overwhelming majority down-regulated in male-sterile floral buds. Enrichment analyses revealed that the general phenylpropanoid pathway as well as its related metabolisms was predicted to be altered in a mutant compared to its fertile progenitor. Meanwhile, the conserved genetic pathway DYT1-TDF1-AMS-MS188-MS1, as well as the causal gene ClAMT1 for the male-sterile mutant Se18, was substantially disrupted during male reproductive development. In addition, some targets of the key regulators AMS and MS188 in tapetum development were also down-regulated at a transcriptional level, such as ABCG26 (Cla004479), ACOS5 (Cla022956), CYP703A2 (Cla021151), PKSA (Cla021099), and TKPR1 (Cla002563). Considering lncRNAs may act as functional endogenous target mimics of miRNAs, competitive endogenous RNA networks were subsequently constructed, with the most complex one containing three DE-miRNAs, two DE-lncRNAs, and 21 DEGs. Collectively, these findings not only contribute to a better understanding of genetic regulatory networks underlying male sterility in watermelon, but also provide valuable candidates for future research.

Investigating nicotine pathway-related long non-coding RNAs in tobacco

Long non-coding RNAs (lncRNAs) are transcripts longer than 200 bp with low or no protein-coding ability, which play essential roles in various biological processes in plants. Tobacco is an ideal model plant for studying nicotine biosynthesis and metabolism, and there is little research on lncRNAs in this field. Therefore, how to take advantage of the mature tobacco system to profoundly investigate the lncRNAs involved in the nicotine pathway is intriguing. By exploiting 549 public RNA-Seq datasets of tobacco, 30,212 lncRNA candidates were identified, including 24,084 large intervening non-coding RNAs (lincRNAs), 5,778 natural antisense transcripts (NATs) and 350 intronic non-coding RNAs (incRNAs). Compared with protein-coding genes, lncRNAs have distinct properties in terms of exon number, sequence length, A/U content, and tissue-specific expression pattern. lincRNAs showed an asymmetric evolutionary pattern, with a higher proportion (68.71%) expressed from the Nicotiana sylvestris (S) subgenome. We predicted the potential cis/trans-regulatory effects on protein-coding genes. One hundred four lncRNAs were detected as precursors of 30 known microRNA (miRNA) family members, and 110 lncRNAs were expected to be the potential endogenous target mimics for 39 miRNAs. By combining the results of weighted gene co-expression network analysis with the differentially expressed gene analysis of topping RNA-seq data, we constructed a sub-network containing eight lncRNAs and 25 nicotine-related coding genes. We confirmed that the expression of seven lncRNAs could be affected by MeJA treatment and may be controlled by the transcription factor NtMYC2 using a quantitative PCR assay and gene editing. The results suggested that lncRNAs are involved in the nicotine pathway. Our findings further deepened the understanding of the features and functions of lncRNAs and provided new candidates for regulating nicotine biosynthesis in tobacco.

Comparative study on abortion characteristics of Nsa CMS and Pol CMS and analysis of long non-coding RNAs related to pollen abortion in Brassica napus

Cytoplasmic male sterile system (CMS) is one of the important methods for the utilization of heterosisin Brassica napus. The involvement of long non-coding RNAs (lncRNAs) in anther and pollen development in B.napus has been recognized, but there is little data on the involvement of lncRNAs in pollen abortion in different types of rapeseed CMS. The present study compared the cytological, physiological and biochemical characteristics of Nsa CMS (1258A) and Pol CMS (P5A) during pollen abortion, and high-throughput sequencing of flower buds of different sizes before and after pollen abortion. The results showed that insufficient energy supply was an important physiological basis for 1258A and P5A pollen abortion, and 1258A had excessive ROS (reactive oxygen species) accumulation in the stage of pollen abortion. Functional analysis showed that Starch and sucrose metabolism and Sulfur metabolism were significantly enriched before and after pollen abortion in 1258A and P5A, and a large number of genes were down-regulated. In 1258A, 227 lncRNAs had cis-targeting regulation, and 240 cis-target genes of the lncRNAs were identified. In P5A, 116 lncRNAs had cis-targeting regulation, and 101 cis-target genes of the lncRNAs were identified. There were five lncRNAs cis-target genes in 1258A and P5A during pollen abortion, and LOC106445716 encodes β-D-glucopyranosyl abscisate β-glucosidase and could regulate pollen abortion. Taken together, this study, provides a new perspective for lncRNAs to participate in the regulation of Nsa CMS and Pol CMS pollen abortion.

Role of long non-coding RNAs in rice reproductive development

Rice is a staple crop, feeding over half of the global population. The future demand of population growth and climate change requires substantial rice improvement. Recent advances in rice genomics have highlighted the vital role of the non-coding part of the genome. The protein-coding regions account for only a tiny portion of the eukaryotic genome, and most of the genomic regions transcribe copious amounts of non-coding RNAs. Of these, the long non-coding RNAs, including linear non-coding RNAs (lncRNAs) and circular non-coding RNAs (circRNAs), have been shown to play critical roles in various developmental processes by regulating the expression of genes and functions of proteins at transcriptional, post-transcriptional and post-translational levels. With the advances in next-generation sequencing technologies, a substantial number of long non-coding RNAs have been found to be expressed in plant reproductive organs in a cell- and tissue-specific manner suggesting their reproductive development-related functions. Accumulating evidence points towards the critical role of these non-coding RNAs in flowering, anther, and pollen development, ovule and seed development and photoperiod and temperature regulation of male fertility. In this mini review, we provide a brief overview of the role of the linear and circular long non-coding RNAs in rice reproductive development and control of fertility and crop yield.

Identification of long non-coding RNAs and microRNAs involved in anther development in the tropical Camellia oleifera

BACKGROUND

Explored the molecular science of anther development is important for improving productivity and overall yield of crops. Although the role of regulatory RNAs, including long non-coding RNAs (lncRNAs) and microRNAs (miRNAs), in regulating anther development has been established, their identities and functions in Camellia oleifera, an important industrial crop, have yet not been clearly explored. Here, we report the identification and characterization of genes, lncRNAs and miRNAs during three stages of the tropical C. oleifera anther development by single-molecule real-time sequencing, RNA sequencing and small RNA sequencing, respectively.

RESULTS

These stages, viz. the pollen mother cells stage, tetrad stage and uninucleate pollen stage, were identified by analyzing paraffin sections of floral buds during rapid expansion periods. A total of 18,393 transcripts, 414 putative lncRNAs and 372 miRNAs were identified, of which 5,324 genes, 115 lncRNAs, and 44 miRNAs were differentially accumulated across three developmental stages. Of these, 44 and 92 genes were predicted be regulated by 37 and 30 differentially accumulated lncRNAs and miRNAs, respectively. Additionally, 42 differentially accumulated lncRNAs were predicted as targets of 27 miRNAs. Gene ontology enrichment indicated that potential target genes of lncRNAs were enriched in photosystem II, regulation of autophagy and carbohydrate phosphatase activity, which are essential for anther development. Functional annotation of genes targeted by miRNAs indicated that they are relevant to transcription and metabolic processes that play important roles in microspore development. An interaction network was built with 2 lncRNAs, 6 miRNAs and 10 mRNAs. Among these, miR396 and miR156 family were up-regulated, while their targets, genes (GROWTH REGULATING FACTORS and SQUAMOSA PROMOTER BINDING PROTEIN-LIKE genes) and lncRNAs, were down-regulated. Further, the trans-regulated targets of these lncRNAs, like wall-associated kinase2 and phosphomannose isomerase1, are involved in pollen wall formation during anther development.

CONCLUSIONS

This study unravels lncRNAs, miRNAs and miRNA-lncRNA-mRNA networks involved in development of anthers of the tropical C. oleifera lays a theoretical foundation for further elucidation of regulatory roles of lncRNAs and miRNAs in anther development.

Transcriptome of miRNA during inhibition of Bombyx mori nuclear polyhedrosis virus by geldanamycin in BmN cells

Bombyx mori nuclear polyhedrosis virus (BmNPV) is one of several viruses that cause great harm to the sericulture industry, and its pathogenic mechanism is still being explored. Geldanamycin (GA), a kind of HSP90 inhibitor, has been verified to suppress BmNPV proliferation. However, the molecular mechanism by which GA inhibits BmNPV is unclear. MicroRNAs (miRNAs) have been shown to play a key role in regulating virus proliferation and host-pathogen interactions. In this study, BmN cells infected with BmNPV were treated by GA and DMSO for 72 h, respectively, then transcriptome analysis of miRNA was performed from the GA group and the control group. As a result, a total of 29 miRNAs were differentially expressed (DE), with 13 upregulated and 16 downregulated. Using bioinformatics analysis, it was found that the target genes of DEmiRNAs were involved in ubiquitin-mediated proteolysis, phagosome, proteasome, endocytosis pathways, and so on. Six DEmiRNAs were verified by quantitative reverse-transcription polymerase chain reaction. DElong noncoding RNA (DElncRNA)-DEmiRNA-messenger RNA (mRNA) regulatory networks involved in apoptosis and immune pathways were constructed in GA-treated BmN cells, which included 12 DEmiRNA, 132 DElncRNA, and 69 mRNAs. This regulatory network enriched the functional role of miRNA in the BmNPV-silkworm interactions and improved our understanding of the molecular mechanism of HSP90 inhibitors on BmNPV proliferation.

Transcriptomic analysis of saffron at different flowering stages using RNA sequencing uncovers cytochrome P450 genes involved in crocin biosynthesis

Saffron is a well-known Chinese traditional herb, and crocin biosynthesis is related to the yield and quality of saffron. This study aimed to screen differentially expressed genes (DEGs) in saffron at different flowering stages and identify cytochrome P450 (CYP) genes involved in crocin biosynthesis. Saffron samples at different flowering stages were used for RNA sequencing, and DEGs between the samples at three days before the flowering stage (- 3da) and two days after the flowering stage (+ 2da) were screened. Thereafter, significantly differentially expressed CYP genes were identified, and CYP gene expression at different flowering stages and in various tissues of saffron was determined using real-time quantitative polymerase chain reaction (RT-qPCR). After sequencing and analysis, 1508 DEGs between the samples at - 3da and + 2da were identified, including 487 upregulated and 1021 downregulated genes, which were enriched in 16 biological processes, 5 cellular components, 3 molecular functions, and 11 KEGG pathways, including protein processing in endoplasmic reticulum, pentose and glucuronate interconversions, starch and sucrose metabolism, estrogen signaling pathway, and mitogen-activated protein kinase signaling pathway. In addition, 12 significantly differentially expressed CYP genes were identified. The RT-qPCR results showed that CYP76C4, CYP72A15, CYP72A219, CYP97B2, CYP714C2, CYP71A1, CYP94C1, and CYP86A8 were all expressed in the pistils, and CYP72A219, CYP72A15, CYP97B2, CYP71A1, and CYP86A8 were highly expressed in the pistils. Our study established a transcriptome library of saffron and found that CYP72A219, CYP72A15, CYP97B2, CYP71A1, and CYP86A8 may be candidates involved crocin biosynthesis in saffron.

Transcriptome analysis reveals the encystment-related lncRNA expression profile and coexpressed mRNAs in Pseudourostyla cristata

Ciliated protozoans form dormant cysts for survival under adverse conditions. The molecular mechanisms regulating this process are critical for understanding how single-celled eukaryotes adapt to the environment. Despite the accumulated data on morphology and gene coding sequences, the molecular mechanism by which lncRNAs regulate ciliate encystment remains unknown. Here, we first detected and analyzed the lncRNA expression profile and coexpressed mRNAs in dormant cysts versus vegetative cells in the hypotrich ciliate Pseudourostyla cristata by high-throughput sequencing and qRT-PCR. A total of 853 differentially expressed lncRNAs were identified. Compared to vegetative cells, 439 and 414 lncRNAs were upregulated and downregulated, respectively, while 47 lncRNAs were specifically expressed in dormant cysts. A lncRNA-mRNA coexpression network was constructed, and the possible roles of lncRNAs were screened. Three of the identified lncRNAs, DN12058, DN20924 and DN30855, were found to play roles in fostering encystment via their coexpressed mRNAs. These lncRNAs can regulate a variety of physiological activities that are essential for encystment, including autophagy, protein degradation, the intracellular calcium concentration, microtubule-associated dynein and microtubule interactions, and cell proliferation inhibition. These findings provide the first insight into the potentially functional lncRNAs and their coexpressed mRNAs involved in the dormancy of ciliated protozoa and contribute new evidence for understanding the molecular mechanisms regulating encystment.