Databases and computational methods for the identification of piRNA-related molecules: A survey

Piwi-interacting RNAs (piRNAs) are a class of small non-coding RNAs (ncRNAs) that plays important roles in many biological processes and major cancer diagnosis and treatment, thus becoming a hot research topic. This study aims to provide an in-depth review of computational piRNA-related research, including databases and computational models. Herein, we perform literature analysis and use comparative evaluation methods to summarize and analyze three aspects of computational piRNA-related research: (i) computational models for piRNA-related molecular identification tasks, (ii) computational models for piRNA-disease association prediction tasks, and (iii) computational resources and evaluation metrics for these tasks. This study shows that computational piRNA-related research has significantly progressed, exhibiting promising performance in recent years, whereas they also suffer from the emerging challenges of inconsistent naming systems and the lack of data. Different from other reviews on piRNA-related identification tasks that focus on the organization of datasets and computational methods, we pay more attention to the analysis of computational models, algorithms, and performances that aim to provide valuable references for computational piRNA-related identification tasks. This study will benefit the theoretical development and practical application of piRNAs by better understanding computational models and resources to investigate the biological functions and clinical implications of piRNA.

LncRNA CARMN facilitates odontogenic differentiation of dental pulp cells by impairing EZH2

OBJECTIVE

This study aimed to reveal the potential role of CARMN in odontogenic differentiation of dental pulp cells (DPCs).

METHODS

Laser capture microdissection was used to detect Carmn in DPCs and odontoblasts in P0 mice. After manipulating CARMN expression in odontogenic differentiation induced hDPCs, the state of odontogenic differentiation was evaluated by ALP staining, ARS, and related marker expression in qRT-PCR and western blotting. The subcutaneous transplantation of HA/β-TCP loaded with hDPCs was performed to verify CARMN's role in promoting odontogenic differentiation in vivo. RNAplex and RIP were employed to reveal potential mechanism of CARMN in hDPCs.

RESULTS

CARMN expressed more abundantly in odontoblasts than DPCs in P0 mice. CARMN expression boosted during in vitro odontogenic differentiation of hDPCs. CARMN overexpression enhanced odontogenic differentiation of hDPCs in vitro, while inhibition impaired the process. CARMN overexpression in HA/β-TCP composites promoted more mineralized nodule formation in vivo. CARMN knockdown led to soared EZH2, while CARMN overexpression brought about EZH2 inhibition. CARMN functioned via direct interaction with EZH2.

CONCLUSIONS

The results uncovered CARMN as a modulator during the odontogenic differentiation of DPCs. CARMN promoted odontogenic differentiation of DPCs by impairing EZH2.

DinoKnot: Duplex Interaction of Nucleic Acids With PseudoKnots

Interaction of nucleic acid molecules is essential for their functional roles in the cell and their applications in biotechnology. While simple duplex interactions have been studied before, the problem of efficiently predicting the minimum free energy structure of more complex interactions with possibly pseudoknotted structures remains a challenge. In this work, we introduce a novel and efficient algorithm for prediction of Duplex Interaction of Nucleic acids with pseudoKnots, DinoKnot follows the hierarchical folding hypothesis to predict the secondary structure of two interacting nucleic acid strands (both homo- and hetero-dimers). DinoKnot utilizes the structure of molecules before interaction as a guide to find their duplex structure allowing for possible base pair competitions. To showcase DinoKnots's capabilities we evaluated its predicted structures against (1) experimental results for SARS-CoV-2 genome and nine primer-probe sets, (2) a clinically verified example of a mutation affecting detection, and (3) a known nucleic acid interaction involving a pseudoknot. In addition, we compared our results against our closest competition, RNAcofold, further highlighting DinoKnot's strengths. We believe DinoKnot can be utilized for various applications including screening new variants for potential detection issues and supporting existing applications involving DNA/RNA interactions, adding structural considerations to the interaction to elicit functional information.

Signaling pathways and regulatory networks in quail skeletal muscle development: insights from whole transcriptome sequencing

Quail, as an advantageous avian model organism due to its compact size and short reproductive cycle, holds substantial potential for enhancing our understanding of skeletal muscle development. The quantity of skeletal muscle represents a vital economic trait in poultry production. Unraveling the molecular mechanisms governing quail skeletal muscle development is of paramount importance for optimizing meat and egg yield through selective breeding programs. However, a comprehensive characterization of the regulatory dynamics and molecular control underpinning quail skeletal muscle development remains elusive. In this study, through the application of HE staining on quail leg muscle sections, coupled with preceding fluorescence quantification PCR of markers indicative of skeletal muscle differentiation, we have delineated embryonic day 9 (E9) and embryonic day 14 (E14) as the start and ending points, respectively, of quail skeletal muscle differentiation. Then, we employed whole transcriptome sequencing to investigate the temporal expression profiles of leg muscles in quail embryos at the initiation of differentiation (E9) and upon completion of differentiation (E14). Our analysis revealed the expression patterns of 12,012 genes, 625 lncRNAs, 14,457 circRNAs, and 969 miRNAs in quail skeletal muscle samples. Differential expression analysis between the E14 and E9 groups uncovered 3,479 differentially expressed mRNAs, 124 lncRNAs, 292 circRNAs, and 154 miRNAs. Furthermore, enrichment analysis highlighted the heightened activity of signaling pathways related to skeletal muscle metabolism and intermuscular fat formation, such as the ECM-receptor interaction, focal adhesion, and PPAR signaling pathway during E14 skeletal muscle development. Conversely, the E9 stage exhibited a prevalence of pathways associated with myoblast proliferation, exemplified by cell cycle processes. Additionally, we constructed regulatory networks encompassing lncRNA‒mRNA, miRNA‒mRNA, lncRNA‒miRNA-mRNA, and circRNA-miRNA‒mRNA interactions, thus shedding light on their putative roles within quail skeletal muscle. Collectively, our findings illuminate the gene and non-coding RNA expression characteristics during quail skeletal muscle development, serving as a foundation for future investigations into the regulatory mechanisms governing non-coding RNA and quail skeletal muscle development in poultry production.

Genome-wide view and characterization of natural antisense transcripts in Cannabis Sativa L

Natural Antisense Transcripts (NATs) are a kind of complex regulatory RNAs that play crucial roles in gene expression and regulation. However, the NATs in Cannabis Sativa L., a widely economic and medicinal plant rich in cannabinoids remain unknown. In this study, we comprehensively predicted C. sativa NATs genome-wide using strand-specific RNA sequencing (ssRNA-Seq) data, and validated the expression profiles by strand-specific quantitative reverse transcription PCR (ssRT-qPCR). Consequently, a total of 307 NATs were predicted in C. sativa, including 104 cis- and 203 trans- NATs. Functional enrichment analysis demonstrated the potential involvement of the C. sativa NATs in DNA polymerase activity, RNA-DNA hybrid ribonuclease activity, and nucleic acid binding. Finally, 18 cis- and 376 trans- NAT-ST pairs were predicted to produce 621 cis- and 5,679 trans- small interfering RNA (nat-siRNAs), respectively. These nat-siRNAs were potentially involved in the biosynthesis of cannabinoids and cellulose. All these results will shed light on the regulation of NATs and nat-siRNAs in C. sativa.

Transcriptome analysis reveals the potential lncRNA-mRNA modules involved in genetic male sterility and fertility of Chinese cabbage (brassica rapa L. ssp. pekinensis)

BACKGROUND

Long non-coding RNAs (lncRNAs) play a crucial role in regulating gene expression vital for the growth and development of plants. Despite this, the role of lncRNAs in Chinese cabbage (Brassica rapa L. ssp. pekinensis) pollen development and male fertility remains poorly understood.

RESULTS

In this study, we characterized a recessive genic male sterile mutant (366-2 S), where the delayed degradation of tapetum and the failure of tetrad separation primarily led to the inability to form single microspores, resulting in male sterility. To analyze the role of lncRNAs in pollen development, we conducted a comparative lncRNA sequencing using anthers from the male sterile mutant line (366-2 S) and the wild-type male fertile line (366-2 F). We identified 385 differentially expressed lncRNAs between the 366-2 F and 366-2 S lines, with 172 of them potentially associated with target genes. To further understand the alterations in mRNA expression and explore potential lncRNA-target genes (mRNAs), we performed comparative mRNA transcriptome analysis in the anthers of 366-2 S and 366-2 F at two stages. We identified 1,176 differentially expressed mRNAs. Remarkably, GO analysis revealed significant enrichment in five GO terms, most notably involving mRNAs annotated as pectinesterase and polygalacturonase, which play roles in cell wall degradation. The considerable downregulation of these genes might contribute to the delayed degradation of tapetum in 366-2 S. Furthermore, we identified 15 lncRNA-mRNA modules through Venn diagram analysis. Among them, MSTRG.9997-BraA04g004630.3 C (β-1,3-glucanase) is associated with callose degradation and tetrad separation. Additionally, MSTRG.5212-BraA02g040020.3 C (pectinesterase) and MSTRG.13,532-BraA05g030320.3 C (pectinesterase) are associated with cell wall degradation of the tapetum, indicating that these three candidate lncRNA-mRNA modules potentially regulate pollen development.

CONCLUSION

This study lays the foundation for understanding the roles of lncRNAs in pollen development and for elucidating their molecular mechanisms in regulating male sterility in Chinese cabbage.

Screening and identification of lncRNAs in preadipocyte differentiation in sheep

Studies of preadipocyte differentiation and fat deposition in sheep have mainly focused on functional genes, and with no emphasis placed on the role that long non-coding RNAs (lncRNAs) may have on the activity of those genes. Here, the expression profile of lncRNAs in ovine preadipocyte differentiation was investigated and the differentially expressed lncRNAs were screened on day 0 (D0), day 2(D2) and day 8(D8) of ovine preadipocyte differentiation, with their target genes being predicted. The competing endogenous RNA (ceRNA) regulatory network was constructed by GO and KEGG enrichment analysis for functional annotation, and some differentially expressed lncRNAs were randomly selected to verify the RNA-Seq results by RT-qPCR. In the study, a total of 2517 novel lncRNAs and 3943 known lncRNAs were identified from ovine preadipocytes at the three stages of differentiation, with the highest proportion being intergenic lncRNAs. A total of 3455 lncRNAs were expressed at all three stages of preadipocyte differentiation, while 214, 226 and 228 lncRNAs were uniquely expressed at day 0, day 2 and day 8, respectively. By comparing the expression of the lncRNAs between the three stages of differentiation stages, a total of 405, 272 and 359 differentially expressed lncRNAs were found in D0-vs-D2, D0-vs-D8, and D2-vs-D8, respectively. Functional analysis revealed that the differentially expressed lncRNAs were enriched in signaling pathways related to ovine preadipocyte differentiation, such as mitogen-activated protein kinase (MAPK) pathway, the phosphoinositide 3-kinase protein kinase B (PI3K-Akt) pathway, and the transforming growth factor beta (TGF-β) pathway. In summary, lncRNAs from preadipocytes at different stages of differentiation in sheep were identified and screened using RNA-Seq technology, and the regulatory mechanisms of lncRNAs in preadipocyte differentiation and lipid deposition were explored. This study provides a theoretical reference for revealing the roles of lncRNAs in ovine preadipocyte differentiation and also offers a theoretical basis for further understanding the regulatory mechanisms of ovine preadipocyte differentiation.

Transcriptional landscape of pathogen-responsive lncRNAs in tomato unveils the role of hydrolase encoding genes in response to Botrytis cinerea invasion

LncRNAs have gained increasing attention owing to their important regulatory roles on growth and stress responses of plants. However, the mechanisms underlying the functions of lncRNAs in fruit-pathogen interaction are still largely unknown. In this study, a total of 273 lncRNAs responding to Botrytis cinerea infection were identified in tomato fruit, among which a higher percentage of antisense lncRNAs were targeted to the genes enriched in hydrolase activity. To ascertain the roles of these lncRNAs, seven hydrolase-related transcripts were transiently knocked-down by virus-induced gene silencing. Silencing of lncRNACXE20 reduced the expression level of a carboxylesterase gene, further enhancing the resistance of tomato to B. cinerea. In contrast, silencing of lncRNACHI, lncRNAMMP, lncRNASBT1.9 and lncRNAPME1.9 impaired the resistance to B. cinerea, respectively. Further RT-qPCR assay and enzymatic activity detection displayed that the attenuated resistance of lncRNAMMP and lncRNASBT1.9-silenced plants was associated with the inhibition on the expression of JA-related genes, while the decreased resistance of lncRNACHI-silenced plants resulted in reduced chitinase activity. Collectively, these results may provide references for deciphering the mechanisms underlying specific lncRNAs to interfere with B. cinerea infection by regulating the expression of defence-related genes or affecting hydrolase activity.

Expression profiles of lncRNAs, miRNAs, and mRNAs and interaction analysis indicate their potential involvement during testicular fusion in Spodoptera litura

Testicular fusion of Spodoptera litura occures during metamorphosis, which benefits sperms development. Previous research identified involvement of ECM-integrin interaction pathways, MMPs in testicular fusion, but the regulatory mechanism remains unclear. RNA-seq was performed to analyze long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) in testes, aiming to uncover potential regulatory mechanisms of testicular fusion. 2150 lncRNAs, 2742 targeted mRNAs, and 347 miRNAs were identified in testes at three different developmental stages. Up-regulated DElncRNAs and DEmRNAs, as well as down-regulated DEmiRNAs, were observed during testicular fusion, while the opposite expression pattern was observed after fusion. Enrichment analysis of DEmRNAs revealed that cAMP signal pathway, ECM remodeling enzymes, ECM-integrin interaction pathways, and cell adhesion molecules were potentially associated with testicular fusion. The identified DElncRNA-DEmiRNA-DEmRNA regulatory network related to cAMP signal pathway, ECM remodeling enzymes suggests their roles during testicular fusion. Our research will provide new targets for studying the mechanism of testicular fusion.

How to do RNA-RNA Interaction Prediction? A Use-Case Driven Handbook Using IntaRNA

Computational prediction of RNA-RNA interactions (RRI) is a central methodology for the specific investigation of inter-molecular RNA interactions and regulatory effects of non-coding RNAs like eukaryotic microRNAs or prokaryotic small RNAs. Available methods can be classified according to their underlying prediction strategies, each implicating specific capabilities and restrictions often not transparent to the non-expert user. Within this work, we review seven classes of RRI prediction strategies and discuss the advantages and limitations of respective tools, since such knowledge is essential for selecting the right tool in the first place.Among the RRI prediction strategies, accessibility-based approaches have been shown to provide the most reliable predictions. Here, we describe how IntaRNA, as one of the state-of-the-art accessibility-based tools, can be applied in various use cases for the task of computational RRI prediction. Detailed hands-on examples for individual RRI predictions as well as large-scale target prediction scenarios are provided. We illustrate the flexibility and capabilities of IntaRNA through the examples. Each example is designed using real-life data from the literature and is accompanied by instructions on interpreting the respective results from IntaRNA output. Our use-case driven instructions enable non-expert users to comprehensively understand and utilize IntaRNA's features for effective RRI predictions.

Comparative RNA Genomics

Over the last quarter of a century it has become clear that RNA is much more than just a boring intermediate in protein expression. Ancient RNAs still appear in the core information metabolism and comprise a surprisingly large component in bacterial gene regulation. A common theme with these types of mostly small RNAs is their reliance of conserved secondary structures. Large-scale sequencing projects, on the other hand, have profoundly changed our understanding of eukaryotic genomes. Pervasively transcribed, they give rise to a plethora of large and evolutionarily extremely flexible non-coding RNAs that exert a vastly diverse array of molecule functions. In this chapter we provide a-necessarily incomplete-overview of the current state of comparative analysis of non-coding RNAs, emphasizing computational approaches as a means to gain a global picture of the modern RNA world.

Discovering Novel Bacterial Small RNA by RNA-seq Analysis Toolkit ANNOgesic

ANNOgesic is an RNA-seq analysis pipeline that can detect sRNAs and many other genomic features in bacteria and archaea. In addition to listing sRNA candidates, ANNOgesic also generates various formats of data files for visual examination and downstream experimental design. Based on validations from previous studies, the sRNA predictions are accurate and reliable. In this chapter, we outline the sRNA detection algorithm, important parameters used, step-by-step execution, and data interpretation with a B. pertussis study as an example. Following those procedures, novel sRNA can be revealed by ANNOgesic.

Rocket-miR, a translational launchpad for miRNA-based antimicrobial drug development

IMPORTANCE

Antimicrobial-resistant infections contribute to millions of deaths worldwide every year. In particular, the group of bacteria collectively known as ESKAPE (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter sp.) pathogens are of considerable medical concern due to their virulence and exceptional ability to develop antibiotic resistance. New kinds of antimicrobial therapies are urgently needed to treat patients for whom existing antibiotics are ineffective. The Rocket-miR application predicts targets of human miRNAs in bacterial and fungal pathogens, rapidly identifying candidate miRNA-based antimicrobials. The application's target audience are microbiologists that have the laboratory resources to test the application's predictions. The Rocket-miR application currently supports 24 recognized human pathogens that are relevant to numerous diseases including cystic fibrosis, chronic obstructive pulmonary disease (COPD), urinary tract infections, and pneumonia. Furthermore, the application code was designed to be easily extendible to other human pathogens that commonly cause hospital-acquired infections.

The response of LncRNAs associated with photosynthesis-and pigment synthesis-related genes to green light in Chlamydomonas reinhardtii

The quality of light is an important abiotic factor that affects the growth and development of green plants. Ultraviolet, red, blue, and far-red light all have demonstrated roles in regulating green plant growth and development, as well as light morphogenesis. However, the mechanism underlying photosynthetic organism responses to green light throughout the life of them are not clear. In this study, we exposed the unicellular green alga Chlamydomonas reinhardtii to green light and analyzed the dynamics of transcriptome changes. Based on the whole transcriptome data from C. reinhardtii, a total of 9974 differentially expressed genes (DEGs) were identified under green light. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that these DEGs were mainly related to "carboxylic acid metabolic process," "enzyme activity," "carbon metabolism," and "photosynthesis and other processes." At the same time, 253 differentially expressed long non-coding RNAs (DELs) were characterized as green light responsive. We also made a detailed analysis of the responses of photosynthesis- and pigment synthesis-related genes in C. reinhardtii to green light and found that these genes exhibited obvious dynamic expression. Lastly, we constructed a co-expression regulatory network, comprising 49 long non-coding RNAs (lncRNAs) and 20 photosynthesis and pigment related genes, of which 9 mRNAs were also the predicted trans/cis-targets of 8 lncRNAs, these results suggested that lncRNAs may affect the expression of mRNAs related to photosynthesis and pigment synthesis. Our findings give a preliminary explanation of the response mechanism of C. reinhardtii to green light at the transcriptional level.

TargetRNA3: predicting prokaryotic RNA regulatory targets with machine learning

Small regulatory RNAs pervade prokaryotes, with the best-studied family of these non-coding genes corresponding to trans-acting regulators that bind via base pairing to their message targets. Given the increasing frequency with which these genes are being identified, it is important that methods for illuminating their regulatory targets keep pace. Using a machine learning approach, we investigate thousands of interactions between small RNAs and their targets, and we interrogate more than a hundred features indicative of these interactions. We present a new method, TargetRNA3, for predicting targets of small RNA regulators and show that it outperforms existing approaches. TargetRNA3 is available at https://cs.wellesley.edu/~btjaden/TargetRNA3 .

Transcriptome and Metabolome Analyses Reveal Sugar and Acid Accumulation during Apricot Fruit Development

The apricot (Prunus armeniaca L.) is a fruit that belongs to the Rosaceae family; it has a unique flavor and is of important economic and nutritional value. The composition and content of soluble sugars and organic acids in fruit are key factors in determining the flavor quality. However, the molecular mechanism of sugar and acid accumulation in apricots remains unclear. We measured sucrose, fructose, glucose, sorbitol, starch, malate, citric acid, titratable acid, and pH, and investigated the transcriptome profiles of three apricots (the high-sugar cultivar 'Shushanggan', common-sugar cultivar 'Sungold', and low-sugar cultivar 'F43') at three distinct developmental phases. The findings indicated that 'Shushanggan' accumulates a greater amount of sucrose, glucose, fructose, and sorbitol, and less citric acid and titratable acid, resulting in a better flavor; 'Sungold' mainly accumulates more sucrose and less citric acid and starch for the second flavor; and 'F43' mainly accumulates more titratable acid, citric acid, and starch for a lesser degree of sweetness. We investigated the DEGs associated with the starch and sucrose metabolism pathways, citrate cycle pathway, glycolysis pathway, and a handful of sugar transporter proteins, which were considered to be important regulators of sugar and acid accumulation. Additionally, an analysis of the co-expression network of weighted genes unveiled a robust correlation between the brown module and sucrose, glucose, and fructose, with VIP being identified as a hub gene that interacted with four sugar transporter proteins (SLC35B3, SLC32A, SLC2A8, and SLC2A13), as well as three structural genes for sugar and acid metabolism (MUR3, E3.2.1.67, and CSLD). Furthermore, we found some lncRNAs and miRNAs that regulate these genes. Our findings provide clues to the functional genes related to sugar metabolism, and lay the foundation for the selection and cultivation of high-sugar apricots in the future.

A Hitchhiker's guide to RNA-RNA structure and interaction prediction tools

RNA biology has risen to prominence after a remarkable discovery of diverse functions of noncoding RNA (ncRNA). Most untranslated transcripts often exert their regulatory functions into RNA-RNA complexes via base pairing with complementary sequences in other RNAs. An interplay between RNAs is essential, as it possesses various functional roles in human cells, including genetic translation, RNA splicing, editing, ribosomal RNA maturation, RNA degradation and the regulation of metabolic pathways/riboswitches. Moreover, the pervasive transcription of the human genome allows for the discovery of novel genomic functions via RNA interactome investigation. The advancement of experimental procedures has resulted in an explosion of documented data, necessitating the development of efficient and precise computational tools and algorithms. This review provides an extensive update on RNA-RNA interaction (RRI) analysis via thermodynamic- and comparative-based RNA secondary structure prediction (RSP) and RNA-RNA interaction prediction (RIP) tools and their general functions. We also highlighted the current knowledge of RRIs and the limitations of RNA interactome mapping via experimental data. Then, the gap between RSP and RIP, the importance of RNA homologues, the relationship between pseudoknots, and RNA folding thermodynamics are discussed. It is hoped that these emerging prediction tools will deepen the understanding of RNA-associated interactions in human diseases and hasten treatment processes.

The Pleiotropic Phenotypes Caused by an hfq Null Mutation in Vibrio harveyi

Hfq is a global regulator and can be involved in multiple cellular processes by assisting small regulatory RNAs (sRNAs) to target mRNAs. To gain insight into the virulence regulation of Hfq in Vibrio harveyi, the hfq null mutant, ∆hfq, was constructed in V. harveyi strain 345. Compared with the wild-type strain, the mortality of pearl gentian sharply declined from 80% to 0% in ∆hfq when infected with a dose that was 7.5-fold the median lethal dose (LD50). Additionally, ∆hfq led to impairments of bacterial growth, motility, and biofilm formation and resistance to reactive oxygen species, chloramphenicol, and florfenicol. A transcriptome analysis indicated that the expression of 16.39% genes on V. harveyi 345 were significantly changed after the deletion of hfq. Without Hfq, the virulence-related pathways, including flagellar assembly and bacterial chemotaxis, were repressed. Moreover, eleven sRNAs, including sRNA0405, sRNA0078, sRNA0419, sRNA0145, and sRNA0097, which, respectively, are involved in chloramphenicol/florfenicol resistance, outer membrane protein synthesis, electron transport, amino acid metabolism, and biofilm formation, were significantly down-regulated. In general, Hfq contributes to the virulence of V. harveyi 345 probably via positively regulating bacterial motility and biofilm formation. It is involved in flagellar assembly and bacterial chemotaxis by binding sRNAs and regulating the target mRNAs.

Diverse Regulatory Manners and Potential Roles of lncRNAs in the Developmental Process of Asian Honey Bee (Apis cerana) Larval Guts

Long non-coding RNAs (lncRNAs) are crucial modulators in a variety of biological processes, such as gene expression, development, and immune defense. However, little is known about the function of lncRNAs in the development of Asian honey bee (Apis cerana) larval guts. Here, on the basis of our previously obtained deep-sequencing data from the 4-, 5-, and 6-day-old larval guts of A. cerana workers (Ac4, Ac5, and Ac6 groups), an in-depth transcriptome-wide investigation was conducted to decipher the expression pattern, regulatory manners, and potential roles of lncRNAs during the developmental process of A. cerana worker larval guts, followed by the verification of the relative expression of differentially expressed lncRNAs (DElncRNAs) and the targeting relationships within a competing endogenous RNA (ceRNA) axis. In the Ac4 vs. Ac5 and Ac5 vs. Ac6 comparison groups, 527 and 498 DElncRNAs were identified, respectively, which is suggestive of the dynamic expression of lncRNAs during the developmental process of larval guts. A cis-acting analysis showed that 330 and 393 neighboring genes of the aforementioned DElncRNAs were respectively involved in 29 and 32 functional terms, such as cellular processes and metabolic processes; these neighboring genes were also respectively engaged in 246 and 246 pathways such as the Hedgehog signaling pathway and the Wnt signaling pathway. Additionally, it was found that 79 and 76 DElncRNAs as potential antisense lncRNAs may, respectively, interact with 72 and 60 sense-strand mRNAs. An investigation of competing endogenous RNA (ceRNA) networks suggested that 75 (155) DElncRNAs in the Ac4 vs. Ac5 (Ac5 vs. Ac6) comparison group could target 7 (5) DEmiRNAs and further bind to 334 (248) DEmRNAs, which can be annotated to 33 (29) functional terms and 186 (210) pathways, including 12 (16) cellular- and humoral-immune pathways (lysosome pathway, necroptosis, MAPK signaling pathway, etc.) and 11 (10) development-associated signaling pathways (Wnt, Hippo, AMPK, etc.). The RT-qPCR detection of five randomly selected DElncRNAs confirmed the reliability of the used sequencing data. Moreover, the results of a dual-luciferase reporter assay were indicative of the binding relationship between MSTRG.11294.1 and miR-6001-y and between miR-6001-y and ncbi_107992440. These results demonstrate that DElncRNAs are likely to modulate the developmental process of larval guts via the regulation of the source genes' transcription, interaction with mRNAs, and ceRNA networks. Our findings not only yield new insights into the developmental mechanism underlying A. cerana larval guts, but also provide a candidate ceRNA axis for further functional dissection.

LinearCoFold and LinearCoPartition: linear-time algorithms for secondary structure prediction of interacting RNA molecules

Many RNAs function through RNA-RNA interactions. Fast and reliable RNA structure prediction with consideration of RNA-RNA interaction is useful, however, existing tools are either too simplistic or too slow. To address this issue, we present LinearCoFold, which approximates the complete minimum free energy structure of two strands in linear time, and LinearCoPartition, which approximates the cofolding partition function and base pairing probabilities in linear time. LinearCoFold and LinearCoPartition are orders of magnitude faster than RNAcofold. For example, on a sequence pair with combined length of 26,190 nt, LinearCoFold is 86.8× faster than RNAcofold MFE mode, and LinearCoPartition is 642.3× faster than RNAcofold partition function mode. Surprisingly, LinearCoFold and LinearCoPartition's predictions have higher PPV and sensitivity of intermolecular base pairs. Furthermore, we apply LinearCoFold to predict the RNA-RNA interaction between SARS-CoV-2 genomic RNA (gRNA) and human U4 small nuclear RNA (snRNA), which has been experimentally studied, and observe that LinearCoFold's prediction correlates better with the wet lab results than RNAcofold's.

Integrated microRNA and whole-transcriptome sequencing reveals the involvement of small and long non-coding RNAs in the fiber growth of ramie plant

BACKGROUND

MicroRNAs (miRNAs) and long non-coding RNAs (lncRNAs) are the two main types of non-coding RNAs that play crucial roles in plant growth and development. However, their specific roles in the fiber growth of ramie plant (Boehmeria nivea L. Gaud) remain largely unknown.

METHODS

In this study, we performed miRNA and whole-transcriptome sequencing of two stem bark sections exhibiting different fiber growth stages to determine the expression profiles of miRNAs, lncRNAs, and protein-encoding genes.

RESULTS

Among the identified 378 miRNAs and 6,839 lncRNAs, 88 miRNAs and 1,288 lncRNAs exhibited differential expression. Bioinformatics analysis revealed that 29 and 228 differentially expressed protein-encoding genes were targeted by differentially expressed miRNAs and lncRNAs, respectively, constituting eight putative competing endogenous RNA networks. lncR00022274 exhibited downregulated expression in barks with growing fibers. It also had an antisense overlap with the MYB gene, BntWG10016451, whose overexpression drastically increased the xylem fiber number and secondary wall thickness of fibers in the stems of transgenic Arabidopsis, suggesting the potential association of lncR00022274-BntWG10016451 expression with fiber growth.

CONCLUSIONS

These findings provide insights into the roles of ncRNAs in the regulation of fiber growth in ramie, which can be used for the biotechnological improvement of its fiber yield and quality in the future.

Regulatory Roles of Long Non-Coding RNAs Relevant to Antioxidant Enzymes and Immune Responses of Apis cerana Larvae Following Ascosphaera apis Invasion

Long non-coding RNAs (lncRNAs) play an essential part in controlling gene expression and a variety of biological processes such as immune defense and stress-response. However, whether and how lncRNAs regulate responses of Apis cerana larvae to Ascosphaera apis invasion has remained unclear until now. Here, the identification and structural analysis of lncRNAs in the guts of A. cerana worker larvae were conducted, and the expression profile of larval lncRNAs during the A. apis infection process was then analyzed, followed by an investigation of the regulatory roles of differentially expressed lncRNAs (DElncRNAs) in the host response. In total, 76 sense lncRNAs, 836 antisense lncRNAs, 184 intron lncRNAs, 362 bidirectional lncRNAs, and 2181 intron lncRNAs were discovered in the larval guts. Additionally, 30 known and 9 novel lncRNAs were potential precursors for 36 and 11 miRNAs, respectively. In the three comparison groups, 386, 351, and 272 DElncRNAs were respectively identified, indicating the change in the overall expression pattern of host lncRNAs following the A. apis invasion. Analysis of cis-acting effect showed that DElncRNAs in the 4-, 5-, and 6-day-old comparison groups putatively regulated 55, 30, and 20 up- and down-stream genes, respectively, which were involved in a series of crucial functional terms and pathways, such as MAPK signaling pathway, and cell process. Analysis showed that 31, 8, and 11 DElncRNAs as potential antisense lncRNAs may interact with 26, 8, and 9 sense-strand mRNAs. Moreover, investigation of the competing endogenous RNA (ceRNA) network indicated that 148, 283, and 257 DElncRNAs were putatively regulated. The expression of target genes by targeting corresponding DEmiRNAs included those associated with antioxidant enzymes and immune responses. These results suggested that DElncRNAs played a potential part in the larval guts responding to the A. apis infection through a cis-acting manner and ceRNA mechanisms. Our findings deepen our understanding of interactions between A. cerana larvae and A. apis and offer a basis for clarifying the DElncRNA-mediated mechanisms underlying the host response to fungal invasion.

Identification of a novel immune-related long noncoding RNA in carp primary macrophages associated with bisphenol A' s immunoregulatory effects

Increasing evidence suggests that long non-coding RNAs (lncRNAs) play pivotal roles in various biological processes. However, current studies on lncRNAs mostly focus on mammalian species, with little research on the functional roles of lncRNAs in teleost fish. Here, we identified a novel intergenic lncRNA (linc-93.2) in the head kidney primary macrophages of common carp (Cyprinus carpio) after exposure to a typical environmental endocrine disrupting chemical, bisphenol A (BPA). As a result, linc-93.2 was more than 3,619 bp in length and predominantly localized to the nucleus of primary macrophages other than cytoplasm, with the highest expression level in spleen followed by head kidney among different organs. Bioinformatic analysis predicted a cis-target gene, dennd1b, and 20 trans-target genes including hsp70, gna13 and rasgap, were potentially regulated by linc-93.2; NFκB and estrogen receptor (ERα) binding sites were located in the promoter region upstream of its transcription start site, which together suggested the involvement of linc-93.2 in immune and neurological functions in fish. Based on that, the expression level of linc-93.2 was determined in macrophages following acute lipopolysaccharide (LPS) and BPA treatments, both of which significantly induced linc-93.2 and IL-1β expression in cells. Moreover, a NF-κB inhibitor PDTC significantly reduced linc-93.2 expression in macrophages, but co-exposure of macrophages to PDTC with BPA or LPS could significantly rescue linc-93.2 expression, consistent with the observation on that LPS or BPA alone significantly induced both linc-93.2 and its target gene expression. Interestingly, linc-93.2 and its target gene expression was significantly suppressed by an ER antagonist ICI 182,780, however, the co-exposure of macrophages to ICI 182,780 with BPA failed to attenuate their declined expression. Overall, the current study demonstrated that linc-93.2, a novel immune-related lncRNA, may participate in the immune processes of common carp macrophages via the NF-κB and ER pathway. The results presented in this study enhance our understanding of the immunotoxin mechanisms of BPA in teleost fish.

Comprehensive transcriptional analysis of pig facial skin development

Background

Skin development is a complex process that is influenced by many factors. Pig skin is used as an ideal material for xenografts because it is more anatomically and physiologically similar to human skin. It has been shown that the skin development of different pig breeds is different, and some Chinese pig breeds have the characteristics of skin thickness and facial skin folds, but the specific regulatory mechanism of this skin development is not yet clear.

Methods

In this study, the facial skin of Chenghua sows in the four developmental stages of postnatal Day 3 (D3) , Day 90 (D90) , Day 180 (D180), and Year 3 (Y3) were used as experimental materials, and RNA sequencing (RNA-seq) analysis was used to explore the changes in RNA expression in skin development at the four developmental stages, determine the differentially expressed messenger RNAs (mRNAs), long noncoding RNAs (lncRNAs), microRNAs (miRNAs), and circular RNAs (circRNAs), and perform functional analysis of related genes by Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses.

Results

A pairwise comparison of the four developmental stages identified several differentially expressed genes (DEGs) and found that the number of differentially expressed RNAs (DE RNAs) increased with increasing developmental time intervals. Elastin (ELN) is an important component of the skin. Its content affects the relaxation of the epidermis and dermal connection, and its expression is continuously downregulated during the four developmental stages. The functions of DEGs at different developmental stages were examined by performing GO and KEGG analyses, and the GO terms and enrichment pathways of mRNAs, lncRNAs, miRNAs, and circRNAs highly overlapped, among which the PPAR signaling pathway, a classical pathway for skin development, was enriched by DEGs of D3 vs. D180, D90 vs. D180 and D180 vs. Y3. In addition, we constructed lncRNA-miRNA-mRNA and circRNA-miRNA interaction networks and found genes that may be associated with skin development, but their interactions need further study.

Conclusions

We identified a number of genes associated with skin development, performed functional analyses on some important DEGs and constructed interaction networks that facilitate further studies of skin development.

Whole Transcriptome Analysis Reveals the Global Molecular Responses of mRNAs, lncRNAs, miRNAs, circRNAs, and Their ceRNA Networks to Salinity Stress in Hong Kong Oysters, Crassostrea hongkongensis

The Hong Kong oyster, Crassostrea hongkongensis, is an estuarine bivalve with remarkable commercial value in South China, and the increase of salinity in estuaries during the dry season has posed a major threat to the oyster farming. To explore the global transcriptional response to salinity stress, a whole-transcriptome analysis was performed with the gills of oysters in 6‰, 18‰, and 30‰ filtered seawater. Overall, 2243, 194, 371, and 167 differentially expressed mRNAs (DEmRNAs), differentially expressed long non-coding RNAs (DElncRNAs), differentially expressed circular RNAs (DEcircRNAs), and differentially expressed microRNAs (DEmiRNAs) were identified, respectively. Based on GO enrichment and KEGG pathway analysis, these important DEmRNAs, DElncRNAs, DEcircRNAs, and DEmiRNAs were predicted to be mainly involved in amino acids metabolism, microtubule movement, and immune defense. This demonstrated the complexity of dynamic transcriptomic profiles of C. hongkongensis in response to salinity fluctuation. The regulatory relationships of DEmiRNAs-DEmRNAs, DElncRNAs-DEmiRNAs, and DEcircRNAs-DEmiRNAs were also predicted, and finally, a circRNA-associated competing endogenous RNA (ceRNA) network was constructed, consisting of six DEcircRNAs, eight DEmiRNAs, and five DEmRNAs. The key roles of taurine and hypotaurine metabolism and phenylalanine metabolism were highlighted in this ceRNA network, which was consistent with the major contribution of free amino acids to intracellular osmolality and cell volume regulation. Collectively, this study provides comprehensive data, contributing to the exploration of coding and non-coding RNAs in C. hongkongensis salinity response. The results would benefit the understanding of the response mechanism of bivalves against salinity fluctuation, and provide clues for genetic improvement of C. hongkongensis with hyper-salinity tolerance.

Multi-omics of Circular RNAs and Their Responses to Hormones in Moso Bamboo (Phyllostachys edulis)

Circular RNAs (circRNAs) are endogenous non-coding RNAs with covalently closed structures, which have important functions in plants. However, their biogenesis, degradation, and function upon treatment with gibberellins (GAs) and auxins (1-naphthaleneacetic acid, NAA) remain unknown. Here, we systematically identified and characterized the expression patterns, evolutionary conservation, genomic features, and internal structures of circRNAs using RNase R-treated libraries from moso bamboo (Phyllostachys edulis) seedlings. Moreover, we investigated the biogenesis of circRNAs dependent on both cis- and trans-regulation. We explored the function of circRNAs, including their roles in regulating microRNA (miRNA)-related genes and modulating the alternative splicing of their linear counterparts. Importantly, we developed a customized degradome sequencing approach to detect miRNA-mediated cleavage of circRNAs. Finally, we presented a comprehensive view of the participation of circRNAs in the regulation of hormone metabolism upon treatment of bamboo seedlings with GA and NAA. Collectively, our study provides insights into the biogenesis, function, and miRNA-mediated degradation of circRNAs in moso bamboo.

LncRNAs regulate the cyclic growth and development of hair follicles in Dorper sheep

Introduction

Hair follicles in Dorper sheep are characterized by seasonal cyclic growth and development, consequently resulting in hair shedding during spring. The cyclic growth and development of hair follicles are regulated by several influencing factors such as photoperiods, hormones, age of the animal, genes, long non-coding RNAs (lncRNAs), and signaling pathways.

Methods

In the present study, skin samples of five shedding sheep (S), used as experimental animals, and three non-shedding sheep (N), used as controls, were collected at three time points (September 27, 2019; January 3, 2020; and March 17, 2020) for RNA sequencing (RNA-seq) technology. Nine different groups (S1-vs-S2, S1-vs-S3, S2-vs-S3, N1- vs-N2, N1-vs-N3, N2-vs-N3, S1-vs-N1, S2-vs-N2, and S3-vs-N3) were compared using FDR < 0.05 and log 21 FC >as thresholds to assess the differences in the expression of lncRNAs.

Results and discussion

In total, 395 differentially expressed (DE) lncRNAs were screened. Cluster heatmap analysis identified two types of expression patterns, namely, high expression during the anagen phase (A pattern) and high expression during the telogen phase (T pattern). Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that the target genes were largely enriched in the Estrogen signaling pathway, PI3K-Akt signaling pathway, Fc gamma R-mediated phagocytosis, and cell adhesion molecules (CAMs), which are associated with hair follicle cyclic growth and development-related pathways. In addition, 17 pairs of lncRNAs-target genes related to hair follicle cyclic growth and development were screened, and a regulatory network was constructed. Altogether, candidate lncRNAs and their regulated target genes were screened that contributed to sheep hair follicle cyclic growth and development. We believe these findings will provide useful insights into the underlying regulatory mechanisms.

Evaluation of lncRNA Expression Pattern and Potential Role in Heart Failure Pathology

During the last few decades, the morbidity and mortality of heart failure (HF) have remained on an upward trend. Despite the advances in therapeutic and diagnostic measures, there are still many aspects requiring further research. This study is aimed at finding potential long noncoding RNAs (lncRNAs) that could aid with the diagnosis and treatment of HF. We performed RNA sequencing on the peripheral blood of healthy controls as well as HF patients. The expression of lncRNAs was validated by RT-qPCR. Bioinformatic analysis was performed to investigate the possible mechanism of differentially expressed lncRNAs and mRNAs. The diagnostic value of lncRNAs was analysed by ROC analysis. Finally, a total of 207 mRNAs and 422 lncRNAs were identified. GO and KEGG pathway analyses revealed that biological pathways such as immune response, regulation of cell membrane, and transcriptional regulatory process were associated with the pathological progress of HF. The lncRNA-mRNA coexpression network was conducted, and several mRNAs were identified as key potential pathological targets, while lncRNA CHST11, MIR29B2CHG, CR381653.1, and FP236383.2 presented a potential diagnostic value for HF. These findings provide novel insights for the underlying mechanisms and possible therapeutic targets for HF.

C-RCPred: a multi-objective algorithm for interactive secondary structure prediction of RNA complexes integrating user knowledge and SHAPE data

RNAs can interact with other molecules in their environment, such as ions, proteins or other RNAs, to form complexes with important biological roles. The prediction of the structure of these complexes is therefore an important issue and a difficult task. We are interested in RNA complexes composed of several (more than two) interacting RNAs. We show how available knowledge on the considered RNAs can help predict their secondary structure. We propose an interactive tool for the prediction of RNA complexes, called C-RCPRed, that considers user knowledge and probing data (which can be generated experimentally or artificially). C-RCPred is based on a multi-objective optimization algorithm. Through an extensive benchmarking procedure, which includes state-of-the-art methods, we show the efficiency of the multi-objective approach and the positive impact of considering user knowledge and probing data on the prediction results. C-RCPred is freely available as an open-source program and web server on the EvryRNA website (https://evryrna.ibisc.univ-evry.fr).

Combined analysis of mRNA-miRNA from testis tissue in Tibetan sheep with different FecB genotypes

Testis size is important for identifying breeding animals with adequate sperm production. The aim of this study was to survey the expression profile of mRNA and miRNA in testis tissue from rams carrying different FecB genotypes, including the wild-type and heterozygous genotypes in Tibetan sheep. Comparative transcriptome profiles for ovine testes were established for wild-type and heterozygote Tibetan sheep by next-generation sequencing. RNA-seq results identified 3,910 (2,034 up- and 1,876 downregulated) differentially expressed (DE) genes and 243 (158 up- and 85 downregulated) DE microRNAs (miRNAs) in wild-type vs heterozygote sheep, respectively. Combined analysis of mRNA-seq and miRNA-seq revealed that 20 miRNAs interacted with 48 true DE target genes in wild-type testes compared to heterozygous genotype testes. These results provide evidence for a functional series of genes operating in Tibetan sheep testis. In addition, quantitative real-time PCR analysis showed that the expression trends of randomly selected DE genes in testis tissues from different genotypes were consistent with high-throughput sequencing results.

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.

Whole-Transcriptome Sequencing Reveals a ceRNA Regulatory Network Associated with the Process of Periodic Albinism under Low Temperature in Baiye No. 1 (Camellia sinensis)

The young shoots of the tea plant Baiye No. 1 display an albino phenotype in the early spring under low environmental temperatures, and the leaves re-green like those of common tea cultivars during the warm season. Periodic albinism is precisely regulated by a complex gene network that leads to metabolic differences and enhances the nutritional value of tea leaves. Here, we identified messenger RNAs (mRNAs), long noncoding RNAs (lncRNAs), circular RNAs (circRNAs), and microRNAs (miRNAs) to construct competing endogenous RNA (ceRNA) regulatory networks. We performed whole-transcriptome sequencing of 12 samples from four periods (Bud, leaves not expanded; Alb, albino leaves; Med, re-greening leaves; and Gre, green leaves) and identified a total of 6325 differentially expressed mRNAs (DEmRNAs), 667 differentially expressed miRNAs (DEmiRNAs), 1702 differentially expressed lncRNAs (DElncRNAs), and 122 differentially expressed circRNAs (DEcircRNAs). Furthermore, we constructed ceRNA networks on the basis of co-differential expression analyses which comprised 112, 35, 38, and 15 DEmRNAs, DEmiRNAs, DElncRNAs, and DEcircRNAs, respectively. Based on the regulatory networks, we identified important genes and their interactions with lncRNAs, circRNAs, and miRNAs during periodic albinism, including the ceRNA regulatory network centered on miR5021x, the GAMYB-miR159-lncRNA regulatory network, and the NAC035-miR319x-circRNA regulatory network. These regulatory networks might be involved in the response to cold stress, photosynthesis, chlorophyll synthesis, amino acid synthesis, and flavonoid accumulation. Our findings provide novel insights into ceRNA regulatory mechanisms involved in Baiye No. 1 during periodic albinism and will aid future studies of the molecular mechanisms underlying albinism mutants.

LncRNA expression analysis by comparative transcriptomics among closely related poplars and their regulatory roles in response to salt stress

Long noncoding RNAs (lncRNAs) play crucial roles in regulating key biological processes; however, our knowledge of lncRNAs' roles in plant adaptive evolution is still limited. Here, we determined the divergence of conserved lncRNAs in closely related poplar species that were either tolerant or sensitive to salt stress by comparative transcriptome analysis. Among the 34,363 identified lncRNAs, approximately 3% were shared among poplar species with conserved sequences but diversified in their function, copy number, originating genomic region and expression patterns. Further cluster analysis revealed that the conserved lncRNAs showed more similar expression patterns within salt-tolerant poplars (P. euphratica and P. pruinosa) than between salt-tolerant and salt-sensitive poplars. Among these lncRNAs, the antisense lncRNA lncERF024 was induced by salt and differentiated expression between salt-sensitive and salt-tolerant poplars. Overexpression of lncERF024 in P. alba var. pyramidalis enhanced poplar tolerance to salt stress. Furthermore, RNA pull-down and RNA-seq analysis showed that numerous candidate genes or proteins associated with stress response and photosynthesis might be involved in salt resistance in PeulncERF024-OE poplars. Altogether, our study provided novel insight into how the diversification of lncRNA expression contributes to plant adaptation traits and showed that lncERF024 may be involved in the regulation both of gene expression and protein function conferring salt tolerance in Populus.

Transcriptome Analysis Reveals the Profile of Long Non-Coding RNAs during Myogenic Differentiation in Goats

The long non-coding RNAs (lncRNAs) are emerging as essential regulators of the growth and development of skeletal muscles. However, little is known about the expression profiles of lncRNAs during the proliferation and differentiation of skeletal muscle satellite cells (MuSCs) in goats. In this study, we investigate potential regulatory lncRNAs that govern muscle development by performing lncRNA expression profiling analysis during the proliferation (cultured in the growth medium, GM) and differentiation (cultured in the differentiation medium, DM1/DM5) of MuSCs. In total, 1001 lncRNAs were identified in MuSC samples, and 314 differentially expressed (DE) (FDR < 0.05, |log2FC| > 1) lncRNAs were screened by pairwise comparisons from three comparison groups (GM-vs-DM1, GM-vs-DM5, DM1-vs-DM5). Moreover, we identified the cis-, trans-, and antisense-regulatory target genes of DE lncRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed that these target genes were significantly enriched in muscle development-related GO terms and KEGG pathways. In addition, the network of interactions between DE lncRNAs and their target genes was identified, which included well-known myogenesis regulators such as Myogenic differentiation 1 (MyoD), Myogenin (MyoG), and Myosin heavy chain (MyHC). Meanwhile, competing endogenous RNA (ceRNA) network analysis showed that 237 DE lncRNAs could bind to 329 microRNAs (miRNAs), while miRNAs could target 564 mRNAs. Together, our results provide a genome-wide resource of lncRNAs that may contribute to myogenic differentiation in goats and lay the groundwork for future investigation into their functions during skeletal muscle development.

Conjoint Analysis of Genome-Wide lncRNA and mRNA Expression during the Salicylic Acid Response in Populus × euramericana

Long noncoding RNAs (lncRNAs) participate in a wide range of biological processes, but lncRNAs in plants remain largely unknown; in particular, we lack a systematic identification of plant lncRNAs involved in hormone responses. To explore the molecular mechanism of the response of poplar to salicylic acid (SA), the changes in protective enzymes, which are closely related to plant resistance induced by exogenous SA, were studied, and the expression of mRNA and lncRNA were determined by high-throughput RNA sequencing. The results showed that the activities of phenylalanine ammonia lyase (PAL) and polyphenol oxidase (PPO), in the leaves of Populus × euramericana, were significantly increased by exogenous SA application. High-throughput RNA sequencing showed that 26,366 genes and 5690 lncRNAs were detected under the different treatment conditions: SA and H₂O application. Among these, 606 genes and 49 lncRNAs were differentially expressed. According to target prediction, lncRNAs and target genes involved in light response, stress response, plant disease resistance, and growth and development, were differentially expressed in SA-treated leaves. Interaction analysis showed that lncRNA-mRNA interactions, following exogenous SA, were involved in the response of poplar leaves to the external environment. Our study provides a comprehensive view of Populus × euramericana lncRNAs and offers insights into the potential functions and regulatory interactions of SA-responsive lncRNAs, thus forming the foundation for future functional analysis of SA-responsive lncRNAs in Populus × euramericana.

Identification of long noncoding RNAs involved in plumule-vernalization of Chinese cabbage

Vernalization is a phenomenon in which plants must undergo a period of continuous low temperatures to change from the vegetative growth stage to the reproductive growth stage. Chinese cabbage is a heading vegetable, and flowering time is an essential developmental trait. Premature vernalization leads to premature bolting, which causes a loss of product value and yield. While research into vernalization has provided a wealth of information, a complete understanding of the molecular mechanism for controlling vernalization requirements has not yet been elucidated. In this study, using high-throughput RNA sequencing, we analyzed the plumule-vernalization response of mRNA and long noncoding RNA in the bolting-resistant Chinese cabbage double haploid (DH) line 'Ju Hongxin' (JHX). A total of 3382 lncRNAs were identified, of which 1553 differentially expressed (DE) lncRNAs were characterized as plumule-vernalization responses. The ceRNA network revealed that 280 ceRNA pairs participated in the plumule-vernalization reaction of Chinese cabbage. Through identifying DE lncRNAs in Chinese cabbage and analyzing anti-, cis-, and trans-functional analysis, some candidate lncRNAs related to vernalization promoting flowering of Chinese cabbage and their regulated mRNA genes were found. Moreover, the expression of several critical lncRNAs and their targets was verified using qRT-PCR. Furthermore, we identified the candidate plumule-vernalization-related long noncoding RNAs that regulate BrFLCs in Chinese cabbage, which was interesting and different from previous studies and was a new discovery. Our findings expand the knowledge of lncRNAs in the vernalization of Chinese cabbage, and the identified lncRNAs provide rich resources for future comparative and functional studies.

MicroRNA Target Identification: Revisiting Accessibility and Seed Anchoring

By pairing to messenger RNAs (mRNAs for short), microRNAs (miRNAs) regulate gene expression in animals and plants. Accurately identifying which mRNAs interact with a given miRNA and the precise location of the interaction sites is crucial to reaching a more complete view of the regulatory network of an organism. Only a few experimental approaches, however, allow the identification of both within a single experiment. Computational predictions of miRNA–mRNA interactions thus remain generally the first step used, despite their drawback of a high rate of false-positive predictions. The major computational approaches available rely on a diversity of features, among which anchoring the miRNA seed and measuring mRNA accessibility are the key ones, with the first being universally used, while the use of the second remains controversial. Revisiting the importance of each is the aim of this paper, which uses Cross-Linking, Ligation, And Sequencing of Hybrids (CLASH) datasets to achieve this goal. Contrary to what might be expected, the results are more ambiguous regarding the use of the seed match as a feature, while accessibility appears to be a feature worth considering, indicating that, at least under some conditions, it may favour anchoring by miRNAs.

Identification and analysis of long non-coding RNAs that are involved in response to GCRV infection in grass carp (Ctenopharyngodon idella)

Long noncoding RNAs (lncRNAs) play important roles in many biological processes including the immune response against virus infection. However, their roles in grass carp reovirus (GCRV) pathogenicity are largely unknown. In this study, the next-generation sequencing (NGS) technology was used to analyze the profiles of lncRNAs in GCRV-infected and mock-infected grass carp kidney (CIK) cells. Our results showed that 37 lncRNAs and 1039 mRNA transcripts exhibited differential expression in CIK cells after GCRV infection compared with the mock infection. Functional analysis through the gene ontology and Kyoto Encyclopedia of Genes and Genomes databases (KEGG) indicated that target genes of the differentially expressed lncRNAs were mainly enriched in the biological processes - biological regulation, cellular process, metabolic process and regulation of the biological process, such as MAPK signaling pathway and Notch signaling. Furthermore, we observed that the lncRNA3076 (ON693852) was markedly upregulated after the GCRV infection. In addition, silencing lncRNA3076 decreased the GCRV replication, which indicates that it might play an important role in the replication of GCRV.

Two Novel lncRNAs Regulate Primordial Germ Cell Development in Zebrafish

Long noncoding RNAs (lncRNAs) are regulatory transcripts in various biological processes. However, the role of lncRNAs in germline development remains poorly understood, especially for fish primordial germ cell (PGC) development. In this study, the lncRNA profile of zebrafish PGC was revealed by single cell RNA-sequencing and bioinformatic prediction. We established the regulation network of lncRNA-mRNA associated with PGC development, from which we identified three novel lncRNAs-lnc172, lnc196, and lnc304-highly expressing in PGCs and gonads. Fluorescent in situ hybridization indicated germline-specific localization of lnc196 and lnc304 in the cytoplasm and nucleus of spermatogonia, spermatocyte, and occyte, and they were co-localized with vasa in the cytoplasm of the spermatogonia. By contrast, lnc172 was localized in the cytoplasm of male germline, myoid cells and ovarian somatic cells. Loss- and gain-of-function experiments demonstrated that knockdown and PGC-specific overexpression of lnc304 as well as universal overexpression of lnc172 significantly disrupted PGC development. In summary, the present study revealed the lncRNA profile of zebrafish PGC and identified two novel lncRNAs associated with PGC development, providing new insights for understanding the regulatory mechanism of PGC development.

Genome-Wide Analysis of microRNAs and Their Target Genes in Dongxiang Wild Rice (Oryza rufipogon Griff.) Responding to Salt Stress

Rice (Oryza sativa) is a staple food for more than half of the world's population, and its production is critical for global food security. Moreover, rice yield decreases when exposed to abiotic stresses, such as salinity, which is one of the most detrimental factors for rice production. According to recent trends, as global temperatures continue to rise due to climate change, more rice fields may become saltier. Dongxiang wild rice (Oryza rufipogon Griff., DXWR) is a progenitor of cultivated rice and has a high tolerance to salt stress, making it useful for studying the regulatory mechanisms of salt stress tolerance. However, the regulatory mechanism of miRNA-mediated salt stress response in DXWR remains unclear. In this study, miRNA sequencing was performed to identify miRNAs and their putative target genes in response to salt stress in order to better understand the roles of miRNAs in DXWR salt stress tolerance. A total of 874 known and 476 novel miRNAs were identified, and the expression levels of 164 miRNAs were found to be significantly altered under salt stress. The stem-loop quantitative real-time PCR (qRT-PCR) expression levels of randomly selected miRNAs were largely consistent with the miRNA sequencing results, suggesting that the sequencing results were reliable. The gene ontology (GO) analysis indicated that the predicted target genes of salt-responsive miRNAs were involved in diverse biological pathways of stress tolerance. This study contributes to our understanding of DXWR salt tolerance mechanisms regulated by miRNAs and may ultimately improve salt tolerance in cultivated rice breeding using genetic methods in the future.

A molecular atlas reveals the tri-sectional spinning mechanism of spider dragline silk

The process of natural silk production in the spider major ampullate (Ma) gland endows dragline silk with extraordinary mechanical properties and the potential for biomimetic applications. However, the precise genetic roles of the Ma gland during this process remain unknown. Here, we performed a systematic molecular atlas of dragline silk production through a high-quality genome assembly for the golden orb-weaving spider Trichonephila clavata and a multiomics approach to defining the Ma gland tri-sectional architecture: Tail, Sac, and Duct. We uncovered a hierarchical biosynthesis of spidroins, organic acids, lipids, and chitin in the sectionalized Ma gland dedicated to fine silk constitution. The ordered secretion of spidroins was achieved by the synergetic regulation of epigenetic and ceRNA signatures for genomic group-distributed spidroin genes. Single-cellular and spatial RNA profiling identified ten cell types with partitioned functional division determining the tri-sectional organization of the Ma gland. Convergence analysis and genetic manipulation further validated that this tri-sectional architecture of the silk gland was analogous across Arthropoda and inextricably linked with silk formation. Collectively, our study provides multidimensional data that significantly expand the knowledge of spider dragline silk generation and ultimately benefit innovation in spider-inspired fibers.

Hyperbaric Oxygen Treatment in Spinal Cord Injury Recovery: Profiling Long Noncoding RNAs

STUDY DESIGN

A functional, transcriptome, and long noncoding RNAs (lncRNAs) expression analysis in the spinal cord of mice after hyperbaric oxygen (HBO) treatment.

OBJECTIVE

We aimed to explore the mechanism by which HBO treats spinal cord injury (SCI) at the level of lncRNAs.

SUMMARY OF BACKGROUND DATA

Immense amounts of research have established that HBO treatment promotes the recovery of neurological function after SCI. The mechanism of action remains to be clarified.

METHODS

High-throughput RNA sequencing, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes enrichment analysis were used to profile lncRNA expression and analyze biological function in the spinal cords of mice from sham-operated, SCI, and HBO-treated groups. The differential expression of lncRNA between the groups was assessed using real-time quantitative polymerase chain reaction.

RESULTS

Differential expression across 577 lncRNAs was identified among the three groups. GO analysis showed that free ubiquitin chain polymerization, ubiquitin homeostasis, DNA replication, synthesis of RNA primer, single-stranded telomeric DNA binding, and alpha-amylase activity were significantly enriched. Kyoto Encyclopedia of Genes and Genomes enrichment analysis displayed that vitamin B6 metabolism, one carbon pool by folate, DNA replication, lysine degradation, beta-alanine metabolism, fanconi anemia pathway, and Notch signal pathway were the main pathways with enrichment significance. LncRNAs NONMMUT 092674.1, NONMMUT042986.2, and NONMMUT018850.2 showed significantly different expression between the SCI and the other two groups (P<0.05, <0.01).

CONCLUSIONS

This study is the first to determine the expression profiles of lncRNAs in the injured spinal cord after HBO treatment. We identified several important dysregulated lncRNAs in this setting. These results help us better understand the mechanism by which HBO treats SCI and provide new potential therapeutic targets for SCI.

RNA sequencing revealed the multi-stage transcriptome transformations during the development of gallbladder cancer associated with chronic inflammation

Gallbladder cancer (GBC) is a highly malignant tumor with extremely poor prognosis. Previous studies have suggested that the carcinogenesis and progression of GBC is a multi-stage and multi-step process, but most of them focused on the genome changes. And a few studies just compared the transcriptome differences between tumor tissues and adjacent noncancerous tissues. The transcriptome changes, relating to every stage of GBC evolution, have rarely been studied. We selected three cases of normal gallbladder, four cases of gallbladder with chronic inflammation induced by gallstones, five cases of early GBC, and five cases of advanced GBC, using next-generation RNA sequencing to reveal the changes in mRNAs and lncRNAs expression during the evolution of GBC. In-depth analysis of the sequencing data indicated that transcriptome changes from normal gallbladder to gallbladder with chronic inflammation were distinctly related to inflammation, lipid metabolism, and sex hormone metabolism; transcriptome changes from gallbladder with chronic inflammation to early GBC were distinctly related to immune activities and connection between cells; and the transcriptome changes from early GBC to advanced GBC were distinctly related to transmembrane transport of substances and migration of cells. Expression profiles of mRNAs and lncRNAs change significantly during the evolution of GBC, in which lipid-based metabolic abnormalities play an important promotive role, inflammation and immune activities play a key role, and membrane proteins are very highlighted molecular changes.

Small RNA Targets: Advances in Prediction Tools and High-Throughput Profiling

MicroRNAs (miRNAs) are an abundant class of small non-coding RNAs that regulate gene expression at the post-transcriptional level. They are suggested to be involved in most biological processes of the cell primarily by targeting messenger RNAs (mRNAs) for cleavage or translational repression. Their binding to their target sites is mediated by the Argonaute (AGO) family of proteins. Thus, miRNA target prediction is pivotal for research and clinical applications. Moreover, transfer-RNA-derived fragments (tRFs) and other types of small RNAs have been found to be potent regulators of Ago-mediated gene expression. Their role in mRNA regulation is still to be fully elucidated, and advancements in the computational prediction of their targets are in their infancy. To shed light on these complex RNA-RNA interactions, the availability of good quality high-throughput data and reliable computational methods is of utmost importance. Even though the arsenal of computational approaches in the field has been enriched in the last decade, there is still a degree of discrepancy between the results they yield. This review offers an overview of the relevant advancements in the field of bioinformatics and machine learning and summarizes the key strategies utilized for small RNA target prediction. Furthermore, we report the recent development of high-throughput sequencing technologies, and explore the role of non-miRNA AGO driver sequences.

Transcriptome analysis reveals the spinal expression profiles of non-coding RNAs involved in anorectal malformations in rat fetuses

BACKGROUND

Despite improvements in anorectal malformation (ARM) therapy, patients might still experience post-operative problems such as fecal incontinence, constipation, and soiling. In particular, the dysplasia of the lumbosacral spinal cord in ARM patients is a major disorder that affects fecal function post-operation. However, the pathological mechanisms involved are still unclear.

METHODS

The non-coding RNAs (ncRNAs) in the lumbosacral spinal cord of fetal rats with ethylenethiourea-induced ARM were identified using RNA sequencing (RNA-seq) and examined to determine their potential function. The lumbosacral spinal cord was isolated on embryonic day 17 for subsequent RNA extraction and RNA-seq. The transcriptome data was analyzed using bioinformatics analysis, followed by validation using quantitative reverse transcription PCR.

RESULTS

Compared to the control group, 26 differentially expressed microRNAs (DEMs; 22 upregulated, 4 downregulated) and 112 differentially expressed long non-coding RNAs (63 upregulated, 49 downregulated) were identified in the ARM group. Several DEMs related to development, namely miR-200a-3p, miR-200b-3p, miR-200c-3p, miR-200a-5p, and miR-429, were selected for further analysis. Notably, compared to the control, the relative expression of miR-200 family members was highly upregulated in ARM fetal rats. Furthermore, GO and KEGG enrichment and miRNA-transcription factor-lncRNA/mRNA network analysis was explored to show molecular mechanism underlying DEMs.

CONCLUSIONS

Our findings suggest the involvement of ncRNAs, especially the miR-200 family members, in the pathogenesis of lumbosacral spinal cord dysplasia in ARM fetal rats.

Identification and characterization of lncRNAs and the interaction of lncRNA-mRNA in Epinephelus coioides induced with Singapore grouper iridovirus infection

Singapore grouper iridovirus (SGIV) is a highly pathogenic double-stranded DNA virus, and the fatality rate of SGIV-infected grouper is more than 90%. Up to now, there is no effective methods to control the disease. Long non-coding RNAs (lncRNAs) might play an important role in individual growth and development, immune regulation and other life processes. In this study, lncRNAs were identified in Epinephelus coioides, an important economic aquaculture marine fish in China and Southeast Asia, and the regulatory relationships of lncRNAs and mRNA response to SGIV infection were analyzed. A total of 11,678 lncRNAs were identified and classified from the spleen and GS (grouper spleen) cells. 105 differentially expressed lncRNAs (DElncRNAs) were detected during SGIV infection. The lncRNAs and the regulated mRNAs were analyzed using co-expression network, lncRNA target gene annotation and GO enrichment. At 24 and 48 h after SGIV infection, 118 and 339 lncRNA-mRNA pairs in GS cells were detected, and 728 and 688 differentially expressed lncRNA-mRNA pairs in spleen were obtained, respectively. GO and KEGG were used to predict the DE lncRNAs' target genes, and deduce the DE lncRNAs-affected signaling pathways. In GS cells, lncRNAs might participate in cell part, binding and catalytic activity; and lncRNAs might be involved in immune system process and transcription factor activity in spleen. These data demonstrated that lncRNAs could regulate the expression of immune-related genes response to viral infection, and providing a new insight into understanding the complexity of immune regulatory networks mediated by lncRNAs during viral infection in teleost fish.

Whole transcriptome analysis of HCT-8 cells infected by Cryptosporidium parvum

BACKGROUND

Cryptosporidium species are zoonotic protozoans that are important causes of diarrhoeal disease in both humans and animals. Non-coding RNAs (ncRNAs) play an important role in the innate immune defense against Cryptosporidium infection, but the underlying molecular mechanisms in the interaction between human ileocecal adenocarcinoma (HCT-8) cells and Cryptosporidium species have not been entirely revealed.

METHODS

The expression profiles of messenger RNAs (mRNAs), long non-coding RNAs (lncRNAs), microRNAs (miRNAs) and circular RNAs (circRNAs) in the early phase of infection of HCT-8 cells with Cryptosporidium parvum and at 3 and 12 h post infection were analyzed using the RNA-sequencing technique. The biological functions of differentially expressed RNAs (dif-RNAs) were discovered through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. The targeting relationships between three ncRNAs and mRNAs were analyzed using bioinformatics methods, followed by building a competing endogenous RNA (ceRNA) regulatory network centered on miRNAs.

RESULTS

After strictly filtering the raw data, our analysis revealed 393 dif-lncRNAs, 69 dif-miRNAs and 115 dif-mRNAs at 3 hpi, and 450 dif-lncRNAs, 129 dif-miRNAs, 117 dif-mRNAs and one dif-circRNA at 12 hpi. Of these, 94 dif-lncRNAs, 24 dif-miRNAs and 22 dif-mRNAs were detected at both post-infection time points. Eleven dif-lncRNAs, seven dif-miRNAs, eight dif-mRNAs and one circRNA were randomly selected and confirmed using the quantitative real-time PCR. Bioinformatics analyses showed that the dif-mRNAs were significantly enriched in nutritional absorption, metabolic processes and metabolism-related pathways, while the dif-lncRNAs were mainly involved in the pathways related to the infection and pathogenicity of C. parvum (e.g. tight junction protein) and immune-related pathways (e.g. cell adhesion molecules). In contrast, dif-miRNAs and dif-circRNA were significantly enriched in apoptosis and apoptosis-related pathways. Among the downregulated RNAs, the miRNAs has-miR-324-3p and hsa-miR-3127-5p appear to be crucial miRNAs which could negatively regulate circRNA, lncRNA and mRNA.

CONCLUSIONS

The whole transcriptome profiles of HCT-8 cells infected with C. parvum were obtained in this study. The results of the GO and KEGG pathway analyses suggest significant roles for these dif-RNAs during the course of C. parvum infection. A ceRNA regulation network containing miRNA at its center was constructed for the first time, with hsa-miR-324-3p and hsa-miR-3127-5p being the crucial miRNAs. These findings provide novel insights into the responses of human intestinal epithelial cells to C. parvum infection.

Maternal NO2 exposure disturbs the long noncoding RNA expression profile in the lungs of offspring in time-series patterns

Gestation is a sensitive window to nitrogen dioxide (NO₂) exposure, which may disturb fetal lung development and lung function later in life. Animal and epidemiological studies indicated that long noncoding RNAs (lncRNAs) participate in abnormal lung development induced by environmental pollutant exposure. In the present study, pregnant C57BL/6J mice were exposed to 2.5 ppm NO₂ (mimicking indoor occupational exposure) or clean air, and lncRNAs expression profiles in the lungs of offspring mice were determined by lncRNA-seq on embryonic day 13.5 (E13.5), E18.5, postnatal day 1 (P1), and P14. The lung histopathology examination of offspring was performed, followed by weighted gene coexpression network analysis (WGCNA), prediction of lncRNAs-target genes, and the biological processes enrichment analysis of lncRNAs. Our results indicated that maternal NO₂ exposure induced hypoalveolarization on P14 and differentially expressed lncRNAs showed a time-series pattern. Following WGCNA and enrichment analysis, 2 modules participated in development-related pathways. Importantly, the expressions of related genes were altered, some of which were confirmed to be related to abnormal vascular development and even lung diseases. The research points out that the maternal NO₂ exposure leads to abnormal lung development in offspring that might be related to altered lncRNAs expression profiles with time-series-pattern.

A multi-network comparative analysis of whole-transcriptome and translatome reveals the effect of high-fat diet on APP/PS1 mice and the intervention with Chinese medicine

Different studies on the effects of high-fat diet (HFD) on Alzheimer’s disease (AD) pathology have reported conflicting findings. Our previous studies showed HFD could moderate neuroinflammation and had no significant effect on amyloid-β levels or contextual memory on AD mice. To gain more insights into the involvement of HFD, we performed the whole-transcriptome sequencing and ribosome footprints profiling. Combined with competitive endogenous RNA analysis, the transcriptional regulation mechanism of HFD on AD mice was systematically revealed from RNA level. Mmu-miR-450b-3p and mmu-miR-6540-3p might be involved in regulating the expression of Th and Ddc expression. MiR-551b-5p regulated the expression of a variety of genes including Slc18a2 and Igfbp3. The upregulation of Pcsk9 expression in HFD intervention on AD mice might be closely related to the increase of cholesterol in brain tissues, while Huanglian Jiedu Decoction significantly downregulated the expression of Pcsk9. Our data showed the close connection between the alterations of transcriptome and translatome under the effect of HFD, which emphasized the roles of translational and transcriptional regulation were relatively independent. The profiled molecular responses in current study might be valuable resources for advanced understanding of the mechanisms underlying the effect of HFD on AD.