lncRNAs in Plant and Animal Sexual Reproduction

Analysis of gonadal transcriptome reveals core long non-coding RNA-mRNA regulatory network in sea cucumber Apostichopus japonicus

Apostichopus japonicus is a representative temperate sea cucumber species, that mainly inhabits in coastal zone of the continental shelf. With high nutritional value and important medical value, A. japonicus become an important commercial aquaculture species and produce significant economic value in recent years. A. japonicus has no sexual dimorphism that can be used to distinguish female and male individuals by external appearance and morphology. The phenotype sex can be only detected by dissecting and observing gonad tissue, thus the breeding efficiency could be greatly reduced. This limitation has hindered the advancement of selective breeding programs and sea cucumber industry. To investigate the genetic basis of reproductive biology in A. japonicus, advanced sequencing techniques, such as next- and third-generation sequencing, have been employed to explore the roles of non-coding RNAs and other genetic factors, offering new insights into sex determination mechanisms. To further gain a deeper understanding of the knowledge underlying lncRNAs in gonadal differentiation, we conducted a comparative transcriptome sequencing analysis of gonadal tissues from both sexes. In our research, a total of 3990 novel lncRNAs and 1441 differentially expressed lncRNAs were identified between female and male gonads. Additionally, a molecular regulatory network indicating lncRNA-mRNA interactions was constructed based on transcriptional profiles, which provide insights into the potential cis- and trans- target genes of lncRNAs. The gonadal transcriptome analysis identified a number of novel long non-coding RNAs involved in female and male reproduction process. Both cis- and trans-acting regulatory networks indicating lncRNA-mRNA interaction were constructed based on transcriptional profiles. These findings provide new insights into the lncRNA-mediated regulation of reproductive biology in marine invertebrates, indicating the crucial roles of long non-coding sequences in regulating expression profiles. Further, the GO and KEGG enrichment analyses of cis- and trans- targeted mRNA for differentially expressed lncRNA indicated that sexual reproduction (GO:0019953), germ cell development (GO:0007281), and negative regulation of hormone secretion (GO:0046888) are potentially involved in gonadal differentiation through the regulation of long non-coding sequences. Notably, besides the classical reproduction related signaling pathway like Gonadotropin-releasing hormone (GnRH) secretion (ko04929), several regulatory pathways, such as Epidermal growth factor receptor (ErbB) signaling pathway (ko04012), TGF-beta signaling pathway (ko04350), and neurotrophin signaling pathway (ko04722) were also enriched and potentially involved in sex differentiation and gonadal development.

Deciphering the Vulnerability of Pollen to Heat Stress for Securing Crop Yields in a Warming Climate

Climate change is leading to more frequent and severe extreme temperature events, negatively impacting agricultural productivity and threatening global food security. Plant reproduction, the process fundamental to crop yield, is highly susceptible to heatwaves, which disrupt pollen development and ultimately affect seed-set and crop yields. Recent research has increasingly focused on understanding how pollen grains from various crops react to heat stress at the molecular and cellular levels. This surge in interest over the last decade has been driven by advances in genomic technologies, such as single-cell RNA sequencing, which holds significant potential for revealing the underlying regulatory reprogramming triggered by heat stress throughout the various stages of pollen development. This review focuses on how heat stress affects gene regulatory networks, including the heat stress response, the unfolded protein response, and autophagy, and discusses the impact of these changes on various stages of pollen development. It highlights the potential of pollen selection as a key strategy for improving heat tolerance in crops by leveraging the genetic variability among pollen grains. Additionally, genome-wide association studies and population screenings have shed light on the genetic underpinnings of traits in major crops that respond to high temperatures during male reproductive stages. Gene-editing tools like CRISPR/Cas systems could facilitate precise genetic modifications to boost pollen heat resilience. The information covered in this review is valuable for selecting traits and employing molecular genetic approaches to develop heat-tolerant genotypes.

Evolutionary comparison of lncRNAs in four cotton species and functional identification of LncR4682-PAS2-KCS19 module in fiber elongation

Long non-coding RNAs (lncRNAs) play an important role in various biological processes in plants. However, there have been few reports on the evolutionary signatures of lncRNAs in closely related cotton species. The lncRNA transcription patterns in two tetraploid cotton species and their putative diploid ancestors were compared in this paper. By performing deep RNA sequencing, we identified 280 429 lncRNAs from 21 tissues in four cotton species. lncRNA transcription evolves more rapidly than mRNAs, and exhibits more severe turnover phenomenon in diploid species compared to that in tetraploid species. Evolutionarily conserved lncRNAs exhibit higher expression levels, and lower tissue specificity compared with species-specific lncRNAs. Remarkably, tissue expression of homologous lncRNAs in Gossypium hirsutum and G. barbadense exhibited similar patterns, suggesting that these lncRNAs may be functionally conserved and selectively maintained during domestication. An orthologous lncRNA, lncR4682, was identified and validated in fibers of G. hirsutum and G. barbadense with the highest conservatism and expression abundance. Through virus-induced gene silencing in upland cotton, we found that lncR4682 and its target genes GHPAS2 and GHKCS19 positively regulated fiber elongation. In summary, the present study provides a systematic analysis of lncRNAs in four closely related cotton species, extending the understanding of transcriptional conservation of lncRNAs across cotton species. In addition, LncR4682-PAS2-KCS19 contributes to cotton fiber elongation by participating in the biosynthesis of very long-chain fatty acids.

An insight into the roles of ubiquitin-specific proteases in plants: development and growth, morphogenesis, and stress response

Ubiquitination is a highly conserved and dynamic post-translational modification in which protein substrates are modified by ubiquitin to influence their activity, localization, or stability. Deubiquitination enzymes (DUBs) counter ubiquitin signaling by removing ubiquitin from the substrates. Ubiquitin-specific proteases (UBPs), the largest subfamily of DUBs, are conserved in plants, serving diverse functions across various cellular processes, although members within the same group often exhibit functional redundancy. Here, we briefly review recent advances in understanding the biological roles of UBPs, particularly the molecular mechanism by which UBPs regulate plant development and growth, morphogenesis, and stress response, which sheds light on the mechanistic roles of deubiquitination in plants.

Transcriptome analysis revealed the characteristics and functions of long non-coding RNAs in the hypothalamus during sexual maturation in goats

The hypothalamus is an essential neuroendocrine area in animals that regulates sexual development. Long non-coding RNAs (lncRNAs) are hypothesized to regulate physiological processes related to animal reproduction. However, the regulatory mechanism by which lncRNAs participate in sexual maturity in goats is poorly known, particularly from birth to sexual maturation. In this study, RNAseq analysis was conducted on the hypothalamus of four developmental stages (1day (D1, n = 5), 2 months (M2, n = 5), 4 months (M4, n = 5), and 6 months (M6, n = 5)) of Jining grey goats. The results showed that a total of 237 differentially expressed lncRNAs (DELs) were identified in the hypothalamus. Among these, 221 DELs exhibited cis-regulatory effects on 693 target genes, while 24 DELs demonstrated trans-regulatory effects on 63 target genes. The target genes of these DELs are mainly involved in biological processes related to energy metabolism, signal transduction and hormone secretion, such as sphingolipid signaling pathway, adipocytokine signaling pathway, neurotrophic signaling pathway, glutamatergic synapse, P53 signaling pathway and GnRH signaling pathway. In addition, XR_001918477.1, TCONS_00077463, XR_001918760.1, and TCONS_00029048 and their potential target genes may play a crucial role in the process of goat sexual maturation. This study advances our understanding of lncRNA in hypothalamic tissue during sexual maturation in goats and will give a theoretical foundation for improving goat reproductive features.

Transcription of the Antisense Long Non-Coding RNA, SUPPRESSOR OF FEMINIZATION, Represses Expression of the Female-Promoting Gene FEMALE GAMETOPHYTE MYB in the Liverwort Marchantia polymorpha

Sexual differentiation is a fundamental process in the life cycles of land plants, ensuring successful sexual reproduction and thereby contributing to species diversity and survival. In the dioicous liverwort Marchantia polymorpha, this process is governed by an autosomal sex-differentiation locus comprising FEMALE GAMETOPHYTE MYB (FGMYB), a female-promoting gene, and SUPPRESSOR OF FEMINIZATION (SUF), an antisense strand-encoded long non-coding RNA (lncRNA). SUF is specifically transcribed in male plants and suppresses the expression of FGMYB, leading to male differentiation. However, the molecular mechanisms underlying this process remain elusive. Here, we show that SUF acts through its transcription to suppress FGMYB expression. Transgene complementation analysis using CRISPR/Cas9D10A-based large-deletion mutants identified a genomic region sufficient for the sex differentiation switch function in the FGMYB-SUF locus. Inserting a transcriptional terminator sequence into the SUF-transcribed region resulted in the loss of SUF function and allowed expression of FGMYB in genetically male plants, leading to conversion of the sex phenotype from male to female. Partial deletions of SUF had no obvious impact on its function. Replacement of the FGMYB sequence with that of an unrelated gene did not affect the ability of SUF transcription to suppress sense-strand expression. Taken together, our findings suggest that the process of SUF transcription, rather than the resulting transcripts, is required for controlling sex differentiation in M. polymorpha.

Emerging roles and mechanisms of lncRNAs in fruit and vegetables

With the development of genome sequencing technologies, many long non-coding RNAs (lncRNAs) have been identified in fruit and vegetables. lncRNAs are primarily transcribed and spliced by RNA polymerase II (Pol II) or plant-specific Pol IV/V, and exhibit limited evolutionary conservation. lncRNAs intricately regulate various aspects of fruit and vegetables, including pigment accumulation, reproductive tissue development, fruit ripening, and responses to biotic and abiotic stresses, through diverse mechanisms such as gene expression modulation, interaction with hormones and transcription factors, microRNA regulation, and involvement in alternative splicing. This review presents a comprehensive overview of lncRNA classification, basic characteristics, and, most importantly, recent advances in understanding their functions and regulatory mechanisms.

Identifying long non-coding RNAs involved in heat stress response during wheat pollen development

Introduction

Wheat is a staple food crop for over one-third of the global population. However, the stability of wheat productivity is threatened by heat waves associated with climate change. Heat stress at the reproductive stage can result in pollen sterility and failure of grain development.

Methods

This study used transcriptome data analysis to explore the specific expression of long non-coding RNAs (lncRNAs) in response to heat stress during pollen development in four wheat cultivars.

Results and discussion

We identified 11,054 lncRNA-producing loci, of which 5,482 lncRNAs showed differential expression in response to heat stress. Heat-responsive lncRNAs could target protein-coding genes in cis and trans and in lncRNA-miRNA-mRNA regulatory networks. Gene ontology analysis predicted that target protein-coding genes of lncRNAs regulate various biological processes such as hormonal responses, protein modification and folding, response to stress, and biosynthetic and metabolic processes. We also noted some paired lncRNA/protein-coding gene modules and some lncRNA-miRNA-mRNA regulatory modules shared in two or more wheat cultivars. These modules were related to regulating plant responses to heat stress, such as heat-shock proteins and transcription factors, and protein domains, such as MADS-box, Myc-type, and Alpha crystallin/Hsp20 domain.

Conclusion

Our results provide the basic knowledge and molecular resources for future functional studies investigating wheat reproductive development under heat stress.

Long noncoding RNAs in human reproductive processes and diseases

Infertility has become a global disease burden. Although assisted reproductive technologies are widely used, the assisted reproduction birth rate is no more than 30% worldwide. Therefore, understanding the mechanisms of reproduction can provide new strategies to improve live birth rates and clinical outcomes of enhanced implantation. Long noncoding RNAs (lncRNAs) have been reported to exert regulatory roles in various biological processes and diseases in many species. In this review, we especially focus on the role of lncRNAs in human reproduction. We summarize the function and mechanisms of lncRNAs in processes vital to reproduction, such as spermatogenesis and maturation, sperm motility and morphology, follicle development and maturation, embryo development and implantation. Then, we highlight the importance and diverse potential of lncRNAs as good diagnostic molecular biomarkers and therapeutic targets for infertility treatment.

Transcriptome analysis reveals core lncRNA-mRNA networks regulating melanization and biomineralization in Patinopecten yessoensis shell-infested by Polydora

BACKGROUND

Patinopecten yessoensis, a large and old molluscan group, has been one of the most important aquaculture shellfish in Asian countries because of its high economic value. However, the aquaculture of the species has recently been seriously affected by the frequent outbreaks of Polydora disease, causing great economic losses. Long non-coding RNAs (lncRNAs) exhibit exhibit crucial effects on diverse biological processes, but still remain poorly studied in scallops, limiting our understanding of the molecular regulatory mechanism of P. yessoensis in response to Polydora infestation.

RESULTS

In this study, a high-throughput transcriptome analysis was conducted in the mantles of healthy and Polydora-infected P. yessoensis by RNA sequencing. A total of 19,133 lncRNAs with 2,203 known and 16,930 novel were identified. The genomic characterizations of lncRNAs showed shorter sequence and open reading frame (ORF) length, fewer number of exons and lower expression levels in comparison with mRNAs. There were separately 2280 and 1636 differentially expressed mRNAs and lncRNAs (DEGs and DELs) detected in diseased individuals. The target genes of DELs were determined by both co-location and co-expression analyses. Functional enrichment analysis revealed that DEGs involved in melanization and biomineralization were significantly upregulated; further, obviously increased melanin granules were observed in epithelial cells of the edge mantle in diseased scallops by histological and TEM study, indicating the crucial role of melanizaiton and biomineralization in P. yessoensis to resist against Polydora infestation. Moreover, many key genes, such as Tyrs, Frizzled, Wnts, calmodulins, Pifs, perlucin, laccase, shell matrix protein, mucins and chitins, were targeted by DELs. Finally, a core lncRNA-mRNA interactive network involved in melanization and biomineralization was constructed and validated by qRT-PCR.

CONCLUSIONS

This work provides valuable resources for studies of lncRNAs in scallops, and adds a new insight into the molecular regulatory mechanisms of P. yessoensis defending against Polydora infestation, which will contribute to Polydora disease control and breeding of disease-resistant varieties in molluscs.

Conserved and non-conserved RNA-target modules in plants: lessons for a better understanding of Marchantia development

A wide variety of functional regulatory non-coding RNAs (ncRNAs) have been identified as essential regulators of plant growth and development. Depending on their category, ncRNAs are not only involved in modulating target gene expression at the transcriptional and post-transcriptional levels but also are involved in processes like RNA splicing and RNA-directed DNA methylation. To fulfill their molecular roles properly, ncRNAs must be precisely processed by multiprotein complexes. In the case of small RNAs, DICER-LIKE (DCL) proteins play critical roles in the production of mature molecules. Land plant genomes contain at least four distinct classes of DCL family proteins (DCL1-DCL4), of which DCL1, DCL3 and DCL4 are also present in the genomes of bryophytes, indicating the early divergence of these genes. The liverwort Marchantia polymorpha has become an attractive model species for investigating the evolutionary history of regulatory ncRNAs and proteins that are responsible for ncRNA biogenesis. Recent studies on Marchantia have started to uncover the similarities and differences in ncRNA production and function between the basal lineage of bryophytes and other land plants. In this review, we summarize findings on the essential role of regulatory ncRNAs in Marchantia development. We provide a comprehensive overview of conserved ncRNA-target modules among M. polymorpha, the moss Physcomitrium patens and the dicot Arabidopsis thaliana, as well as Marchantia-specific modules. Based on functional studies and data from the literature, we propose new connections between regulatory pathways involved in Marchantia's vegetative and reproductive development and emphasize the need for further functional studies to understand the molecular mechanisms that control ncRNA-directed developmental processes.

Non-coding RNAs-mediated environmental surveillance determines male fertility in plants

Plants are continuously exposed to environmental stresses leading to significant yield losses. With the changing climatic conditions, the intensity and duration of these stresses are expected to increase, posing a severe threat to crop productivity worldwide. Male gametogenesis is one of the most sensitive developmental stages. Exposure to environmental stresses during this stage leads to male sterility and yield loss. Elucidating the underlying molecular mechanism of environment-affected male sterility is essential to address this challenge. High-throughput RNA sequencing studies, loss-of-function phenotypes of sRNA biogenesis genes and functional genomics studies with non-coding RNAs have started to unveil the roles of small RNAs, long non-coding RNAs and the complex regulatory interactions between them in regulating male fertility under different growth regimes. Here, we discuss the current understanding of the non-coding RNA-mediated environmental stress surveillance and regulation of male fertility in plants. The candidate ncRNAs emerging from these studies can be leveraged to generate environment-sensitive male sterile lines for hybrid breeding or mitigate the impact of climate change on male fertility, as the situation demands.

Comprehensive analysis of lncRNA-mRNA regulatory network in Populus associated with xylem development

Long noncoding RNAs (lncRNAs) play essential roles in numerous biological processes in plants, such as regulating the gene expression. However, only a few studies have looked into their potential functions in xylem development. High-throughput sequencing of P. euramericana 'Zhonglin46' developing and mature xylem was performed in this study. Through sequencing analysis, 14,028 putative lncRNA transcripts were identified, including 4525 differentially expressed lncRNAs (DELs). Additional research revealed that in mature xylem, a total of 2320 DELs were upregulated and 2205 were downregulated compared to developing xylem. Meanwhile, there were a total of 8122 differentially expressed mRNAs (DEMs) that were upregulated and 16,424 that were downregulated in mature xylem compared with developing xylem. The cis- and trans-target genes of DELs were analyzed for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, which indicated that these DELs participate in controlling the phenylpropanoid and lignin biosynthesis pathway as well as the starch and sucrose metabolism pathway. Among the cis-regulated DELs, LNC_006291, LNC_006292, and LNC_006532 all participate in regulating multiple HCT gene family membranes. As targets, POPTR_001G045900v3 (CCR2) and POPTR_018G063500v3 (SUS) both have only one cis-regulatory lncRNA, referred to as LNC_000057 and LNC_006212, respectively. Moreover, LNC_004484 and two DELs named LNC_008014 and LNC_010781 were revealed to be important nodes in the co-expression network of trans-lncRNAs and mRNAs associated to the lignin biosynthesis pathway and cellulose and xylan biosynthetic pathways, respectively. Finally, quantitative real-time PCR (qRT-PCR) was used to confirme 34 pairs of lncRNA-mRNA. Taken together, these findings may help to clarify the regulatory role that lncRNAs play in xylem development and wood formation.

Transcriptome analysis reveals mRNAs and long non-coding RNAs associated with fecundity in the hypothalamus of high-and low-fecundity goat

As an important organ that coordinates the neuroendocrine system, the hypothalamus synthesizes and secretes reproductive hormones that act on the goat organism, thereby precisely regulating follicular development and reproductive processes in goats. However, it is still elusive to explore the mechanism of hypothalamic effects on goat fertility alone. Therefore, RNA-seq was used to analyze the gene expression in hypothalamic tissues of goats in high fertility group (HFG: litter size per litter ≥2) and low fertility group (LFG: litter size per litter = 1), and identified the differential lncRNAs and mRNAs and their associated pathways related to their fertility. The results showed that a total of 23 lncRNAs and 57 mRNAs were differentially expressed in the hypothalamic tissue of high and low fertility goats. GO terms and KEGG functional annotation suggest that DE lncRNAs and DE mRNAs were significantly enriched in hormone-related pathways regulating ovarian development, hormone synthesis and secretion, regulation of reproductive processes, Estrogen signaling pathway, Oxytocin signaling pathway and GnRH signaling pathway. And we constructed a co-expression network of lncRNAs and target genes, and identified reproduction-related genes such as NMUR2, FEZF1, and WT1. The sequencing results of the hypothalamic transcriptome have broadened our understanding of lncRNA and mRNA in goat hypothalamic tissue and provided some new insights into the molecular mechanisms of follicle development and regulation of its fertility in goats.

Genome-wide analysis reveals the crucial role of lncRNAs in regulating the expression of genes controlling pollen development

The genomic location and stage-specific expression pattern of many long non-coding RNAs reveal their critical role in regulating protein-coding genes crucial in pollen developmental progression and male germ line specification. Long non-coding RNAs (lncRNAs) are transcripts longer than 200 bp with no apparent protein-coding potential. Multiple investigations have revealed high expression of lncRNAs in plant reproductive organs in a cell and tissue-specific manner. However, their potential role as essential regulators of molecular processes involved in sexual reproduction remains largely unexplored. We have used developing field mustard (Brassica rapa) pollen as a model system for investigating the potential role of lncRNAs in reproductive development. Reference-based transcriptome assembly performed to update the existing genome annotation identified novel expressed protein-coding genes and long non-coding RNAs (lncRNAs), including 4347 long intergenic non-coding RNAs (lincRNAs, 1058 expressed) and 2,045 lncRNAs overlapping protein-coding genes on the opposite strand (lncNATs, 780 expressed). The analysis of expression profiles reveals that lncRNAs are significant and stage-specific contributors to the gene expression profile of developing pollen. Gene co-expression networks accompanied by genome location analysis identified 38 cis-acting lincRNA, 31 cis-acting lncNAT, 7 trans-acting lincRNA and 14 trans-acting lncNAT to be substantially co-expressed with target protein-coding genes involved in biological processes regulating pollen development and male lineage specification. These findings provide a foundation for future research aiming at developing strategies to employ lncRNAs as regulatory tools for gene expression control during reproductive development.

RNA-seq analysis of differentially expressed LncRNAs from leishmaniasis patients compared to uninfected humans

Leishmaniasis is a parasitic disease that seriously endangers human health. Furthermore, among the parasitic diseases, leishmaniasis is the third most common cause of death after malaria and schistosomiasis. However, the potential function of LncRNAs in leishmaniasis remain unclear. This study aimed to explore the differentially expressed LncRNAs in leishmaniasis. The sera of leishmaniasis patients and uninfected persons for controls were obtained and analyzed by high-throughput sequencing. Moreover, the expression of key LncRNAs was detected by qPCR. The results showed that 970 differentially expressed LncRNAs and 1692 differentially expressed mRNAs were screened compared to control groups. Then, 520 target genes were identified by using bioinformation analysis and the ENCORI database. The bioinformatics analysis revealed that the differentially expressed target genes were enriched in autophagy animal, FoxO signaling pathway, mTOR signaling pathway, and apoptosis, et al. Among those differentially expressed LncRNAs, nine key LncRNAs were selected (MALAT1, NUTM2A-AS1, and LINC00963 had high expression; LINC00622, MAPKAPK5-AS1, LINC02289, XPC-AS1, ZFAS1 and SNHG5 had low expression) by qPCR. This study suggests that different expressions of LncRNAs may involve in the potential function in leishmaniasis and provide a novel insight for diagnosis of this zoonotic disease.

Combined population transcriptomic and genomic analysis reveals cis-regulatory differentiation of non-coding RNAs in maize

Long intergenic non-coding RNA (lincRNA), cis-acting expression quantitative trait locus (cis-eQTL), maize, regulatory evolution. The law of genetic variation during domestication explains the evolutionary mechanism and provides a theoretical basis for improving existing varieties of maize. Previous studies focused on exploiting regulatory variations controlling the expression of protein-coding genes rather than of non-protein-coding genes. Here, we examined the genetic and evolutionary features of long non-coding RNAs from intergenic regions (long intergenic non-coding RNAs, lincRNAs) using population-scale transcriptome data and identified 1168 lincRNAs with cis-acting expression quantitative trait loci (cis-eQTLs). We found that lincRNAs are more likely to be regulated by cis-eQTLs, which exert stronger effects than the protein-coding genes. During maize domestication and improvement, upregulated alleles of lincRNAs, which originated from both standing variation and new mutation, accumulate more frequently and show larger effect sizes than the coding genes. A stronger signature of genetic differentiation was observed in their regulatory regions compared to those of randomly sampled lincRNAs. In addition, we found that cis-regulatory differentiation of lincRNAs is related to the sequence conservation of lincRNA transcripts. Non-conserved lincRNAs more tend to gain upregulated alleles and show a stronger relationship with selected traits than conserved lincRNAs between maize and its wild relatives. Our findings in maize improve the understanding of cis-regulatory variation in lincRNA genes during domestication and improvement and provide an effective approach for prioritizing candidates for further investigation.

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.

Single-cell RNA-seq analyses show that long non-coding RNAs are conspicuously expressed in Schistosoma mansoni gamete and tegument progenitor cell populations

Schistosoma mansoni is a flatworm that causes schistosomiasis, a neglected tropical disease that affects over 200 million people worldwide. New therapeutic targets are needed with only one drug available for treatment and no vaccine. Long non-coding RNAs (lncRNAs) are transcripts longer than 200 nucleotides with low or no protein-coding potential. In other organisms, they have been shown as involved with reproduction, stem cell maintenance and drug resistance, and they tend to exhibit tissue-specific expression patterns. S. mansoni expresses thousands of lncRNA genes; however, the cell type expression patterns of lncRNAs in the parasite remain uncharacterized. Here, we have re-analyzed publicly available single-cell RNA-sequencing (scRNA-seq) data obtained from adult S. mansoni to identify the lncRNAs signature of adult schistosome cell types. A total of 8023 lncRNAs (79% of all lncRNAs) were detected. Analyses of the lncRNAs expression profiles in the cells using statistically stringent criteria were performed to identify 74 lncRNA gene markers of cell clusters. Male gamete and tegument progenitor lineages clusters contained most of the cluster-specific lncRNA markers. We also identified lncRNA markers of specific neural clusters. Whole-mount in situ hybridization (WISH) and double fluorescence in situ hybridization were used to validate the cluster-specific expression of 13 out of 16 selected lncRNA genes (81%) in the male and female adult parasite tissues; for one of these 16 gene loci, probes for two different lncRNA isoforms were used, which showed differential isoform expression in testis and ovary. An atlas of the expression profiles across the cell clusters of all lncRNAs detected in our analysis is available as a public website resource (http://verjolab.usp.br:8081). The results presented here give strong support to a tissue-specific expression and to a regulated expression program of lncRNAs in S. mansoni. This will be the basis for further exploration of lncRNA genes as potential therapeutic targets.

Transcriptome-guided annotation and functional classification of long non-coding RNAs in Arabidopsis thaliana

Long non-coding RNAs (lncRNAs) are a prominent class of eukaryotic regulatory genes. Despite the numerous available transcriptomic datasets, the annotation of plant lncRNAs remains based on dated annotations that have been historically carried over. We present a substantially improved annotation of Arabidopsis thaliana lncRNAs, generated by integrating 224 transcriptomes in multiple tissues, conditions, and developmental stages. We annotate 6764 lncRNA genes, including 3772 that are novel. We characterize their tissue expression patterns and find 1425 lncRNAs are co-expressed with coding genes, with enriched functional categories such as chloroplast organization, photosynthesis, RNA regulation, transcription, and root development. This improved transcription-guided annotation constitutes a valuable resource for studying lncRNAs and the biological processes they may regulate.

Genome-Wide Analysis of Long Non-coding RNAs Involved in Nodule Senescence in Medicago truncatula

Plant long non-coding RNAs (lncRNAs) are widely accepted to play crucial roles during diverse biological processes. In recent years, thousands of lncRNAs related to the establishment of symbiosis, root nodule organogenesis and nodule development have been identified in legumes. However, lncRNAs involved in nodule senescence have not been reported. In this study, senescence-related lncRNAs were investigated in Medicago truncatula nodules by high-throughput strand-specific RNA-seq. A total of 4576 lncRNAs and 126 differentially expressed lncRNAs (DElncRNAs) were identified. We found that more than 60% lncRNAs were associated with transposable elements, especially TIR/Mutator and Helitron DNA transposons families. In addition, 49 DElncRNAs were predicted to be the targets of micro RNAs. Functional analysis showed that the largest sub-set of differently expressed target genes of DElncRNAs were associated with the membrane component. Of these, nearly half genes were related to material transport, suggesting that an important function of DElncRNAs during nodule senescence is the regulation of substance transport across membranes. Our findings will be helpful for understanding the functions of lncRNAs in nodule senescence and provide candidate lncRNAs for further research.

Transcriptome and DNA methylome provide insights into the molecular regulation of drought stress in sea buckthorn

Sea buckthorn is a typical drought-resistant tree species. However, there is a general lack of understanding of the pattern of DNA methylation linked with sea buckthorn responses to drought, and its relationship with drought tolerance mechanisms. In this study, we performed whole-transcriptome RNA sequencing and methylome sequencing in response to drought stress to explore differentially expressed mRNAs, miRNAs, lncRNAs and circRNAs in sea buckthorn leaves. Based on predicted DE pairs, we constructed a competitive endogenous RNA network, which revealed potential transcriptional regulatory roles in response to drought stress. The results of methylome sequencing revealed that the DNA methylation level was increased in sea buckthorn leaves under drought stress. We identified 13,405 differentially methylated regions between CK and TR. We found one DMR-associated DEG (Vacuolar-sorting receptor 6) involved in the ABA accumulation pathway. In addition, two DNA methyltransferases (HrMET1 and HrDRM1) were closely associated with drought-induced hypermethylation in sea buckthorn. Together, we firstly conducted a comprehensive transcriptomic and epigenetic analysis of sea buckthorn under drought stress, providing a resource for further study of the potential functions of genes, miRNAs, lncRNAs, circRNAs and DNA methyltransferases.

Genome-wide identification of lncRNAs, miRNAs, mRNAs and their regulatory networks involved in tapping panel dryness in rubber tree (Hevea brasiliensis)

Noncoding RNAs (ncRNAs) play pivotal roles in various biological processes in plants. However, the role of ncRNAs in tapping panel dryness (TPD) of rubber tree (Hevea brasiliensis Muell. Arg.) is largely unknown. Here, the whole transcriptome analyses of bark tissues from healthy and TPD trees were performed to identify differentially expressed long ncRNAs (DELs), microRNAs/miRNAs (DEMs), genes (DEGs) and their regulatory networks involved in TPD. A total of 263 DELs, 174 DEMs and 1574 DEGs were identified in the bark of TPD tree compared with that of healthy tree. Kyoto Encyclopedia of Genes and Genomes analysis revealed that most of the DEGs and targets of DELs and DEMs were mainly enriched in metabolic pathways, biosynthesis of secondary metabolites and plant hormone signal transduction. Additionally, the majority of DEGs and DELs related to rubber biosynthesis were downregulated in TPD trees. Furthermore, 98 DEGs and 44 DELs were targeted by 54 DEMs, 190 DEGs were identified as putative targets of 56 DELs, and 2 and 44 DELs were predicted as precursors and endogenous target mimics of 2 and 6 DEMs, respectively. Based on these, the DEL-DEM-DEG regulatory network involved in TPD was constructed, and 13 hub DELs, 3 hub DEMs and 2 hub DEGs were identified. The results provide novel insights into the regulatory roles of ncRNAs underlying TPD and lay a foundation for future functional characterization of long ncRNAs, miRNAs and genes involved in TPD in rubber tree.

Biological Parts for Engineering Abiotic Stress Tolerance in Plants

It is vital to ramp up crop production dramatically by 2050 due to the increasing global population and demand for food. However, with the climate change projections showing that droughts and heatwaves becoming common in much of the globe, there is a severe threat of a sharp decline in crop yields. Thus, developing crop varieties with inbuilt genetic tolerance to environmental stresses is urgently needed. Selective breeding based on genetic diversity is not keeping up with the growing demand for food and feed. However, the emergence of contemporary plant genetic engineering, genome-editing, and synthetic biology offer precise tools for developing crops that can sustain productivity under stress conditions. Here, we summarize the systems biology-level understanding of regulatory pathways involved in perception, signalling, and protective processes activated in response to unfavourable environmental conditions. The potential role of noncoding RNAs in the regulation of abiotic stress responses has also been highlighted. Further, examples of imparting abiotic stress tolerance by genetic engineering are discussed. Additionally, we provide perspectives on the rational design of abiotic stress tolerance through synthetic biology and list various bioparts that can be used to design synthetic gene circuits whose stress-protective functions can be switched on/off in response to environmental cues.

Long Intergenic Noncoding RNA (lincRNA) Discovery from Non-Strand-Specific RNA-Seq Data

Long noncoding RNAs (lncRNAs) are transcripts over 200 base pairs in length without discernible protein coding potential. Long intergenic noncoding RNAs (lincRNAs) constitute a subset of lncRNAs, which do not overlap protein coding genes. Here we describe a detailed pipeline for lincRNA discovery from publicly available non-stranded RNA-Seq datasets. The pipeline presented can be applied to any plant species for which RNA-Seq data and a reference genome sequence are available.

High-Throughput Evolutionary Comparative Analysis of Long Intergenic Noncoding RNAs in Multiple Organisms

Comparative genomic and transcriptomic analyses can help prioritize and facilitate the functional analysis of long noncoding RNAs (lncRNAs). Evolinc-II is a bioinformatic pipeline that automates comparative analyses, searching for sequence and structural conservation for thousands of lncRNAs at once. In addition, Evolinc-II takes a phylogenetic approach to infer key evolutionary events that may have occurred during the emergence of each query lncRNA. Here, we describe how to use command line or GUI (CyVerse's Discovery Environment) versions of Evolinc-II to identify lncRNA homologs and prioritize them for functional analysis.

Two genes related to reproductive development in the juvenile prawn, Macrobrachium nipponense: Molecular characterization and transcriptional response to nanoplastic exposure

Plastic pollution in the form of nanoplastics poses a global threat to aquatic ecosystems and the organisms inhabiting them. However, few studies have been conducted on the effects of nanoplastic exposure on reproductive development in crustaceans. In order to address this issue, juvenile oriental river prawns (Macrobrachium nipponense) were exposed to different concentrations of 75-nm polystyrene nanoplastics (0, 5, 10, 20, 40 mg/L) for 28 days. In order to study the regulation of reproduction-related genes in the presence of nanoplastics, the Wee1 protein kinase gene (Wee1) and OTU domain ubiquitin aldehyde binding protein gene (OTUB) were selected. In this study, for the first time, the full-length cDNA of Mn-Wee1 and Mn-OTUB were cloned from M. nipponense. Homologous alignments revealed that Mn-Wee1 had a highly conserved function-critical sequence, and that Mn-OTUB was more closely related to OTUB1 than OTUB2. With increasing concentration of nanoplastics, the expression of both genes increased initially, then decreased. The inhibition of expression of Wee1 and OTUB occurred in 40 mg/L group, respectively. Analysis of the data also indicated that nanoplastic exposure might have differing effects on gene expression in M. nipponense male and female reproductive organs.

Circular RNAs Repertoire and Expression Profile during Brassica rapa Pollen Development

Circular RNAs (circRNAs) are covalently closed RNA molecules generated by the back-splicing of exons from linear precursor mRNAs. Though various linear RNAs have been shown to play important regulatory roles in many biological and developmental processes, little is known about the role of their circular counterparts. In this study, we performed high-throughput RNA sequencing to delineate the expression profile and potential function of circRNAs during the five stages of pollen development in Brassica rapa. A total of 1180 circRNAs were detected in pollen development, of which 367 showed stage-specific expression patterns. Functional enrichment and metabolic pathway analysis showed that the parent genes of circRNAs were mainly involved in pollen-related molecular and biological processes such as mitotic and meiotic cell division, DNA processes, protein synthesis, protein modification, and polysaccharide biosynthesis. Moreover, by predicting the circRNA–miRNA network from our differentially expressed circRNAs, we found 88 circRNAs with potential miRNA binding sites, suggesting their role in post-transcriptional regulation of the genes. Finally, we confirmed the back-splicing sites of nine selected circRNAs using divergent primers and Sanger sequencing. Our study presents the systematic analysis of circular RNAs during pollen development and forms the basis of future studies for unlocking complex gene regulatory networks underpinning reproduction in flowering plants.

License to Regulate: Noncoding RNA Special Agents in Plant Meiosis and Reproduction

The complexity of the subcellular processes that take place during meiosis requires a significant remodeling of cellular metabolism and dynamic changes in the organization of chromosomes and the cytoskeleton. Recently, investigations of meiotic transcriptomes have revealed additional noncoding RNA factors (ncRNAs) that directly or indirectly influence the course of meiosis. Plant meiosis is the point at which almost all known noncoding RNA-dependent regulatory pathways meet to influence diverse processes related to cell functioning and division. ncRNAs have been shown to prevent transposon reactivation, create germline-specific DNA methylation patterns, and affect the expression of meiosis-specific genes. They can also influence chromosome-level processes, including the stimulation of chromosome condensation, the definition of centromeric chromatin, and perhaps even the regulation of meiotic recombination. In many cases, our understanding of the mechanisms underlying these processes remains limited. In this review, we will examine how the different functions of each type of ncRNA have been adopted in plants, devoting attention to both well-studied examples and other possible functions about which we can only speculate for now. We will also briefly discuss the most important challenges in the investigation of ncRNAs in plant meiosis.

Identification and characterization of immune-related lncRNAs and lncRNA-miRNA-mRNA networks of Paralichthys olivaceus involved in Vibrio anguillarum infection

BACKGROUND

Long non-coding RNAs (lncRNAs) structurally resemble mRNAs and exert crucial effects on host immune defense against pathogen infection. Japanese flounder (Paralichthys olivaceus) is an economically important marine fish susceptible to Vibrio anguillarum infection. To date, study on lncRNAs in flounder is scarce.

RESULTS

Here, we reported the first systematic identification and characterization of flounder lncRNAs induced by V. anguillarum infection at different time points. A total of 2,368 lncRNAs were identified, 414 of which were differentially expressed lncRNAs (DElncRNAs) that responded significantly to V. anguillarum infection. For these DElncRNAs, 3,990 target genes (named DETGs) and 42 target miRNAs (named DETmiRs) were identified based on integrated analyses of lncRNA-mRNA and lncRNA-miRNA expressions, respectively. The DETGs were enriched in a cohort of functional pathways associated with immunity. In addition to modulating mRNAs, 36 DElncRNAs were also found to act as competitive endogenous RNAs (ceRNAs) that regulate 37 DETGs through 16 DETmiRs. The DETmiRs, DElncRNAs, and DETGs formed ceRNA regulatory networks consisting of 114 interacting DElncRNAs-DETmiRs-DETGs trinities spanning 10 immune pathways.

CONCLUSIONS

This study provides a comprehensive picture of lncRNAs involved in V. anguillarum infection. The identified lncRNAs and ceRNA networks add new insights into the anti-bacterial immunity of flounder.

From "Dark Matter" to "Star": Insight Into the Regulation Mechanisms of Plant Functional Long Non-Coding RNAs

Long non-coding RNAs (lncRNAs) play a vital role in a variety of biological functions in plant growth and development. In this study, we provided an overview of the molecular mechanisms of lncRNAs in interacting with other biomolecules with an emphasis on those lncRNAs validated only by low-throughput experiments. LncRNAs function through playing multiple roles, including sponger for sequestering RNA or DNA, guider or decoy for recruiting or hijacking transcription factors or peptides, and scaffold for binding with chromatin modification complexes, as well as precursor of microRNAs or small interfering RNAs. These regulatory roles have been validated in several plant species with a comprehensive list of 73 lncRNA-molecule interaction pairs in 16 plant species found so far, suggesting their commonality in the plant kingdom. Such initial findings of a small number of functional plant lncRNAs represent the beginning of what is to come as lncRNAs with unknown functions were found in orders of magnitude more than proteins.

Study of Sexual-Linked Genes (OGI and MeGI) on the Performance of Androecious Persimmons (Diospyros kaki Thunb.)

It is reported that the production of floral sexual phenotype in hexaploid monoecious persimmon (Diospyros kaki) is closely related to a pseudogene called OGI, and a short interspersed nuclear element (SINE)-like insertion (named Kali) in the OGI promoter leads to the gene silence. As a result, DNA methylation level of MeGI promoter determines the development of male or female flowers. However, the molecular mechanism in androecious D. kaki, which only bear male flowers, remains elusive. Here, real-time quantitative polymerase chain reaction (RT-qPCR), molecular cloning, and bisulfite PCR sequencing technique were carried out using 87 materials, including 56 androecious resources, 15 monoecious, and 16 gynoecious cultivars, to investigate the performance of OGI and MeGI on the specific androecious type of D. kaki in China. In conclusion, the Kali insertion was exactly located in the OGI promoter region, and the OGI gene and the Kali sequence were existing and conserved in androecious D. kaki. Meanwhile, we also demonstrated that the MeGI gene was widespread in our investigated samples. Ultimately, our result convincingly provided evidence that the low expression of OGI is probably ascribed to the presence of Kali displaying strong methylation in the OGI promoter, and low expression of MeGI, as well as high DNA methylation level, in the promoter was closely connected with the production of male flowers; this result was consistent with the monoecious persimmon model. Our findings provide predominant genetic aspects for investigation into androecious D. kaki, and future perfecting the sex-determining mechanisms in persimmon.

In Silico Analysis of Common Long Noncoding RNAs in Schistosoma mansoni and Schistosoma haematobium

Background

Schistosomiasis caused by Schistosoma parasites is one of the most common parasitic infections worldwide. Genetic regulation of the genus Schistosoma, which has different developmental stages throughout its life, is quite complex. In these parasites, thousands of long noncoding RNAs (lncRNAs) estimated to be functional were identified. Identifying the transcripts expressed in common and detecting their functions for better understanding of the role of these lncRNAs require a comparative study.

Methods

Assembled RNA-seq datasets belonging to S. mansoni and S. haematobium were obtained from the National Center for Biotechnology. A basic local alignment search tool (BLASTN) analysis was conducted against previously constructed lncRNA library to identify the common lncRNAs between two species. LncRNAs target genes and their gene ontology annotation was performed.

Results

In S. mansoni and S. haematobium, 5132 and 3589 lncRNA transcripts were detected, respectively. These two species had 694 lncRNAs in common. A significant number of lncRNAs was determined to be transcribed from sex chromosomes. The frequently expressed lncRNAs appear to be involved in metabolic and biological regulation processes.

Conclusions

These two species share similar lncRNAs; thus, this finding is a clue that they might have similar functions. In sexual development, they especially might play important roles. Our results will provide important clues to further studies about interactions between human hosts and parasites and the infection mechanisms of Schistosoma parasites.

Barley Anther and Meiocyte Transcriptome Dynamics in Meiotic Prophase I

In flowering plants, successful germinal cell development and meiotic recombination depend upon a combination of environmental and genetic factors. To gain insights into this specialized reproductive development program we used short- and long-read RNA-sequencing (RNA-seq) to study the temporal dynamics of transcript abundance in immuno-cytologically staged barley (Hordeum vulgare) anthers and meiocytes. We show that the most significant transcriptional changes in anthers occur at the transition from pre-meiosis to leptotene-zygotene, which is followed by increasingly stable transcript abundance throughout prophase I into metaphase I-tetrad. Our analysis reveals that the pre-meiotic anthers are enriched in long non-coding RNAs (lncRNAs) and that entry to meiosis is characterized by their robust and significant down regulation. Intriguingly, only 24% of a collection of putative meiotic gene orthologs showed differential transcript abundance in at least one stage or tissue comparison. Argonautes, E3 ubiquitin ligases, and lys48 specific de-ubiquitinating enzymes were enriched in prophase I meiocyte samples. These developmental, time-resolved transcriptomes demonstrate remarkable stability in transcript abundance in meiocytes throughout prophase I after the initial and substantial reprogramming at meiosis entry and the complexity of the regulatory networks involved in early meiotic processes.

Computational methods for annotation of plant regulatory non-coding RNAs using RNA-seq

Plant transcriptome encompasses numerous endogenous, regulatory non-coding RNAs (ncRNAs) that play a major biological role in regulating key physiological mechanisms. While studies have shown that ncRNAs are extremely diverse and ubiquitous, the functions of the vast majority of ncRNAs are still unknown. With ever-increasing ncRNAs under study, it is essential to identify, categorize and annotate these ncRNAs on a genome-wide scale. The use of high-throughput RNA sequencing (RNA-seq) technologies provides a broader picture of the non-coding component of transcriptome, enabling the comprehensive identification and annotation of all major ncRNAs across samples. However, the detection of known and emerging class of ncRNAs from RNA-seq data demands complex computational methods owing to their unique as well as similar characteristics. Here, we discuss major plant endogenous, regulatory ncRNAs in an RNA sample followed by computational strategies applied to discover each class of ncRNAs using RNA-seq. We also provide a collection of relevant software packages and databases to present a comprehensive bioinformatics toolbox for plant ncRNA researchers. We assume that the discussions in this review will provide a rationale for the discovery of all major categories of plant ncRNAs.

Evidences for a role of two Y-specific genes in sex determination in Populus deltoides

Almost all plants in the genus Populus are dioecious (i.e. trees are either male or female), but it is unknown whether dioecy evolved in a common ancestor or independently in different subgenera. Here, we sequence the small peritelomeric X- and Y-linked regions of P. deltoides chromosome XIX. Two genes are present only in the Y-linked region. One is a duplication of a non-Y-linked, female-specifically expressed response regulator, which produces siRNAs that block this gene's expression, repressing femaleness. The other is an LTR/Gypsy transposable element family member, which generates long non-coding RNAs. Overexpression of this gene in A. thaliana promotes androecium development. We also find both genes in the sex-determining region of P. simonii, a different poplar subgenus, which suggests that they are both stable components of poplar sex-determining systems. By contrast, only the duplicated response regulator gene is present in the sex-linked regions of P. davidiana and P. tremula. Therefore, findings in our study suggest dioecy may have evolved independently in different poplar subgenera.

Genome-Wide Network of lncRNA-mRNA During Ovine Oocyte Development From Germinal Vesicle to Metaphase II in vitro

Long non-coding RNA (lncRNA) is involved in many biological processes, and it has been closely investigated. However, research into the role of lncRNA in ovine ovarian development is scant and poorly understood, particularly in relation to the molecular mechanisms of ovine oocyte maturation. In the current study, RNA sequencing was performed with germinal vesicle (GV) and in vitro matured metaphase II (MII) stage oocytes, isolated from ewes. Through the use of bioinformatic analysis, abundant candidate lncRNAs in stage-specific ovine oocytes were identified, and their trans- and cis-regulatory effects were deeply dissected using computational prediction software. Functional enrichment analysis of these lncRNAs revealed that they were involved in the regulation of many key signaling pathways during ovine oocyte development, which was reflected by their targeted genes. From this study, multiple lncRNA-mRNA networks were presumed to be involved in key signaling pathways regarding ovine oocyte maturation and meiotic resumption. In particular, one novel lncRNA (MSTRG.17927) appeared to mediate the regulation of phosphatidylinositol 3-kinase signaling (PI3K) signaling during ovine oocyte maturation. Therefore, this research offers novel insights into the molecular mechanisms underlying ovine oocyte meiotic maturation regulated by lncRNA-mRNA networks from a genome-wide perspective.

Plant pan-genomes are the new reference

Recent years have seen a surge in plant genome sequencing projects and the comparison of multiple related individuals. The high degree of genomic variation observed led to the realization that single reference genomes do not represent the diversity within a species, and led to the expansion of the pan-genome concept. Pan-genomes represent the genomic diversity of a species and includes core genes, found in all individuals, as well as variable genes, which are absent in some individuals. Variable gene annotations often show similarities across plant species, with genes for biotic and abiotic stress commonly enriched within variable gene groups. Here we review the growth of pan-genomics in plants, explore the origins of gene presence and absence variation, and show how pan-genomes can support plant breeding and evolution studies.

Global Role of Crop Genomics in the Face of Climate Change

The development of climate change resilient crops is necessary if we are to meet the challenge of feeding the growing world's population. We must be able to increase food production despite the projected decrease in arable land and unpredictable environmental conditions. This review summarizes the technological and conceptual advances that have the potential to transform plant breeding, help overcome the challenges of climate change, and initiate the next plant breeding revolution. Recent developments in genomics in combination with high-throughput and precision phenotyping facilitate the identification of genes controlling critical agronomic traits. The discovery of these genes can now be paired with genome editing techniques to rapidly develop climate change resilient crops, including plants with better biotic and abiotic stress tolerance and enhanced nutritional value. Utilizing the genetic potential of crop wild relatives (CWRs) enables the domestication of new species and the generation of synthetic polyploids. The high-quality crop plant genome assemblies and annotations provide new, exciting research targets, including long non-coding RNAs (lncRNAs) and cis-regulatory regions. Metagenomic studies give insights into plant-microbiome interactions and guide selection of optimal soils for plant cultivation. Together, all these advances will allow breeders to produce improved, resilient crops in relatively short timeframes meeting the demands of the growing population and changing climate.

Genome-wide identification and comparison of differentially expressed profiles of miRNAs and lncRNAs with associated ceRNA networks in the gonads of Chinese soft-shelled turtle, Pelodiscus sinensis

Background

The gonad is the major factor affecting animal reproduction. The regulatory mechanism of the expression of protein-coding genes involved in reproduction still remains to be elucidated. Increasing evidence has shown that ncRNAs play key regulatory roles in gene expression in many life processes. The roles of microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) in reproduction have been investigated in some species. However, the regulatory patterns of miRNA and lncRNA in the sex biased expression of protein coding genes remains to be elucidated. In this study, we performed an integrated analysis of miRNA, messenger RNA (mRNA), and lncRNA expression profiles to explore their regulatory patterns in the female ovary and male testis of Pelodiscus sinensis.

Results

We identified 10,446 mature miRNAs, 20,414 mRNAs and 28,500 lncRNAs in the ovaries and testes, and 633 miRNAs, 11,319 mRNAs, and 10,495 lncRNAs showed differential expression. A total of 2814 target genes were identified for miRNAs. The predicted target genes of these differentially expressed (DE) miRNAs and lncRNAs included abundant genes related to reproductive regulation. Furthermore, we found that 189 DEmiRNAs and 5408 DElncRNAs showed sex-specific expression. Of these, 3 DEmiRNAs and 917 DElncRNAs were testis-specific, and 186 DEmiRNAs and 4491 DElncRNAs were ovary-specific. We further constructed complete endogenous lncRNA-miRNA-mRNA networks using bioinformatics, including 103 DEmiRNAs, 636 DEmRNAs, and 1622 DElncRNAs. The target genes for the differentially expressed miRNAs and lncRNAs included abundant genes involved in gonadal development, including Wt1, Creb3l2, Gata4, Wnt2, Nr5a1, Hsd17, Igf2r, H2afz, Lin52, Trim71, Zar1, and Jazf1.

Conclusions

In animals, miRNA and lncRNA as master regulators regulate reproductive processes by controlling the expression of mRNAs. Considering their importance, the identified miRNAs, lncRNAs, and their targets in P. sinensis might be useful for studying the molecular processes involved in sexual reproduction and genome editing to produce higher quality aquaculture animals. A thorough understanding of ncRNA-based cellular regulatory networks will aid in the improvement of P. sinensis reproductive traits for aquaculture.

Dynamic Expression of Long Non-Coding RNAs Throughout Parasite Sexual and Neural Maturation in Schistosoma japonicum

Schistosoma japonicum is a flatworm that causes schistosomiasis, a neglected tropical disease. S. japonicum RNA-Seq analyses has been previously reported in the literature on females and males obtained during sexual maturation from 14 to 28 days post-infection in mouse, resulting in the identification of protein-coding genes and pathways, whose expression levels were related to sexual development. However, this work did not include an analysis of long non-coding RNAs (lncRNAs). Here, we applied a pipeline to identify and annotate lncRNAs in 66 S. japonicum RNA-Seq publicly available libraries, from different life-cycle stages. We also performed co-expression analyses to find stage-specific lncRNAs possibly related to sexual maturation. We identified 12,291 S. japonicum expressed lncRNAs. Sequence similarity search and synteny conservation indicated that some 14% of S. japonicum intergenic lncRNAs have synteny conservation with S. mansoni intergenic lncRNAs. Co-expression analyses showed that lncRNAs and protein-coding genes in S. japonicum males and females have a dynamic co-expression throughout sexual maturation, showing differential expression between the sexes; the protein-coding genes were related to the nervous system development, lipid and drug metabolism, and overall parasite survival. Co-expression pattern suggests that lncRNAs possibly regulate these processes or are regulated by the same activation program as that of protein-coding genes.

Genome-wide identification of lncRNAs during hickory (Carya cathayensis) flowering

Non-coding RNAs with lengths greater than 200 bp are known as long non-coding RNAs (lncRNAs), and these RNAs play important role in gene regulation and plant development. However, to date, little is known regarding the role played by lncRNAs during flowering in hickory (Carya cathayensis). Here, we performed whole transcriptome RNA-sequencing of samples from hickory female and male floral buds, in which three samples (H0311PF, H0318PF, and H0402PF) represent pre-flowering, flowering, and post-flowering, respectively, while eight male samples collected from May 8th to June 13th as this time course are the key stage for male floral bud differentiation. We identified 2163 lncRNAs in hickory during flowering, containing 213 intronic, 1488 intergenic, and 462 antisense lncRNAs. We noticed that 510 and 648 lncRNAs were differentially expressed corresponding to female and male floral buds, respectively. And some of the lncRNAs were in a tightly tissue-specific or stage-specific manner. To further understand the roles of the lncRNAs, we predicted the function of the lncRNAs in cis- and trans-acting modes. The results showed that 924 lncRNAs were cis-correlated with 1536 protein-coding genes, while 1207 lncRNAs co-expressed (trans-acting) with 7432 protein-coding genes (R > 0.95 or R < - 0.95). These lncRNAs were all enriched in flower development-associated biological processes, i.e., circadian rhythm, vernalization response, response to gibberellin, inflorescence development, floral organ development, etc. To further understand the relationships between lncRNAs and floral-core genes, we build a co-expressing lncRNA-mRNA flowering network. We classified these floral genes into different pathway (photoperiod, vernalization, gibberellin, autonomous, and sucrose pathway) according to their particular functions. We found a set of lncRNAs that preferentially expressed in these pathways. The network showed that some lncRNAs (i.e., XLOC_038669, XLOC_017938) functioned in a particular flowering time pathway, while others (i.e., XLOC_011251, XLOC_04110) were involved in multiple pathway. Furthermore, some lncRNAs (i.e., XLOC_038669, XLOC_009597, and XLOC_049539) played roles in single or multiple pathways via interaction with each other. This study provides a genome-wide survey of hickory flower-related lncRNAs and will contribute to further understanding of the molecular mechanism underpinning flowering in hickory.

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.

Step-by-Step Bioinformatics Analysis of Schistosoma mansoni Long Non-coding RNA Sequences

In the last few years, long non-coding RNAs (lncRNAs) have been widely studied in humans, and their relevance for physiological and pathological conditions has been demonstrated. In parasites, there are only a few works, such as in Plasmodium falciparum, where it was shown that an lncRNA regulates the expression of a gene associated with immune system evasion, also indicating the relevance of understanding the role of this class of RNAs in parasites. In Schistosoma mansoni, in the last 2 years, there were four published articles related to the annotation of lncRNAs in different life cycle stages using RNA-Seq libraries. In order to make this process of lncRNA identification and annotation more accessible to biologists with no bioinformatics training, considering the growing number of S. mansoni RNA-Seq libraries publicly available from different sources, such as ovary tissues from bi-sex and single-sex infections, and the potential of lncRNAs as therapeutic targets, we provide this step-by-step protocol of lncRNA identification and quantification. This guide includes the download of RNA-Seq libraries from a public database and reads processing and mapping against the genome, transcript reconstruction, novel lncRNA identification, transcripts expression level determination, and the identification of differentially expressed lncRNAs.

Weighted Gene Co-Expression Analyses Point to Long Non-Coding RNA Hub Genes at Different Schistosoma mansoni Life-Cycle Stages

Long non-coding RNAs (lncRNAs) (>200 nt) are expressed at levels lower than those of the protein-coding mRNAs, and in all eukaryotic model species where they have been characterized, they are transcribed from thousands of different genomic loci. In humans, some four dozen lncRNAs have been studied in detail, and they have been shown to play important roles in transcriptional regulation, acting in conjunction with transcription factors and epigenetic marks to modulate the tissue-type specific programs of transcriptional gene activation and repression. In Schistosoma mansoni, around 10,000 lncRNAs have been identified in previous works. However, the limited number of RNA-sequencing (RNA-seq) libraries that had been previously assessed, together with the use of old and incomplete versions of the S. mansoni genome and protein-coding transcriptome annotations, have hampered the identification of all lncRNAs expressed in the parasite. Here we have used 633 publicly available S. mansoni RNA-seq libraries from whole worms at different stages (n = 121), from isolated tissues (n = 24), from cell-populations (n = 81), and from single-cells (n = 407). We have assembled a set of 16,583 lncRNA transcripts originated from 10,024 genes, of which 11,022 are novel S. mansoni lncRNA transcripts, whereas the remaining 5,561 transcripts comprise 120 lncRNAs that are identical to and 5,441 lncRNAs that have gene overlap with S. mansoni lncRNAs already reported in previous works. Most importantly, our more stringent assembly and filtering pipeline has identified and removed a set of 4,293 lncRNA transcripts from previous publications that were in fact derived from partially processed mRNAs with intron retention. We have used weighted gene co-expression network analyses and identified 15 different gene co-expression modules. Each parasite life-cycle stage has at least one highly correlated gene co-expression module, and each module is comprised of hundreds to thousands lncRNAs and mRNAs having correlated co-expression patterns at different stages. Inspection of the top most highly connected genes within the modules’ networks has shown that different lncRNAs are hub genes at different life-cycle stages, being among the most promising candidate lncRNAs to be further explored for functional characterization.

CNIT: a fast and accurate web tool for identifying protein-coding and long non-coding transcripts based on intrinsic sequence composition

As more and more high-throughput data has been produced by next-generation sequencing, it is still a challenge to classify RNA transcripts into protein-coding or non-coding, especially for poorly annotated species. We upgraded our original coding potential calculator, CNCI (Coding-Non-Coding Index), to CNIT (Coding-Non-Coding Identifying Tool), which provides faster and more accurate evaluation of the coding ability of RNA transcripts. CNIT runs ∼200 times faster than CNCI and exhibits more accuracy compared with CNCI (0.98 versus 0.94 for human, 0.95 versus 0.93 for mouse, 0.93 versus 0.92 for zebrafish, 0.93 versus 0.92 for fruit fly, 0.92 versus 0.88 for worm, and 0.98 versus 0.85 for Arabidopsis transcripts). Moreover, the AUC values of 11 animal species and 27 plant species showed that CNIT was capable of obtaining relatively accurate identification results for almost all eukaryotic transcripts. In addition, a mobile-friendly web server is now freely available at http://cnit.noncode.org/CNIT.

Revisiting the Female Germline and Its Expanding Toolbox

The Arabidopsis thaliana ovule arises as a female reproductive organ composed solely of somatic diploid cells. Among them, one cell will acquire a unique identity and initiate female germline development. In this review we explore the complex network that facilitates differentiation of this single cell, and consider how it becomes committed to a distinct developmental program. We highlight recent progress towards understanding the role of intercellular communication, cell competency, and cell-cycle regulation in the ovule primordium, and we discuss the possibility that distinct pathways restrict germline development at different stages. Importantly, these recent findings suggest a renaissance in plant ovule research, restoring the female germline as an attractive model to study cell communication and cell fate establishment in multicellular organs.

MCRiceRepGP: a framework for the identification of genes associated with sexual reproduction in rice

Rice is an important cereal crop, being a staple food for over half of the world's population, and sexual reproduction resulting in grain formation underpins global food security. However, despite considerable research efforts, many of the genes, especially long intergenic non-coding RNA (lincRNA) genes, involved in sexual reproduction in rice remain uncharacterized. With an increasing number of public resources becoming available, information from different sources can be combined to perform gene functional annotation. We report the development of MCRiceRepGP, a machine learning framework which integrates heterogeneous evidence and employs multicriteria decision analysis and machine learning to predict coding and lincRNA genes involved in sexual reproduction in rice. The rice genome was reannotated using deep-sequencing transcriptomic data from reproduction-associated tissue/cell types identifying previously unannotated putative protein-coding genes and lincRNAs. MCRiceRepGP was used for genome-wide discovery of sexual reproduction associated coding and lincRNA genes. The protein-coding and lincRNA genes identified have distinct expression profiles, with a large proportion of lincRNAs reaching maximum expression levels in the sperm cells. Some of the genes are potentially linked to male- and female-specific fertility and heat stress tolerance during the reproductive stage. MCRiceRepGP can be used in combination with other genome-wide studies, such as genome-wide association studies, giving greater confidence that the genes identified are associated with the biological process of interest. As more data, especially about mutant plant phenotypes, become available, the power of MCRiceRepGP will grow, providing researchers with a tool to identify candidate genes for future experiments. MCRiceRepGP is available as a web application (http://mcgplannotator.com/MCRiceRepGP/).