piRBase: integrating piRNA annotation in all aspects

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.

Mimicking bacterial infection in male mice changes sperm small RNA profiles and multigenerationally alters offspring behavior and physiology

Paternal pre-conceptual exposures, including stress, diet, substance abuse, parasite infection, and viral immune activation via Poly I:C, have been reported to influence the brains and behavior of offspring through sperm epigenetic changes. However, the effects of paternal (F0) pre-conceptual exposure to bacterial-induced immune activation on the behavior and physiology of F1 and F2 generations remain unexplored. We examined this using C57BL/6J mice. Eight-week-old males (F0) received a single intraperitoneal injection of the bacterial mimetic lipopolysaccharide (LPS: 5 mg/kg) or 0.9 % saline (vehicle control) before mating with naïve females at four weeks post-injection. Comprehensive behavioral assessments were conducted to investigate anxiety, social behaviors, depressive-like behaviors and cognition in both the F1 and F2 generations within the age range of 8 to 14 weeks. Results demonstrated that only female offspring of LPS-exposed fathers exhibited reduced anxiety levels in the light/dark box, large open field, and novelty-suppressed feeding test. These F1 female offspring also exhibited heightened sociability in the 3-chambered social interaction test and a reduced preference for saccharin in the saccharin preference test. Additionally, the F1 male offspring of LPS-challenged males demonstrated an increased total distance traveled in the light/dark box and a longer distance covered in the light zone. They also exhibited diminished preference for social novelty in the 3-chambered social interaction test and an elevated novel arm preference index in the Y-maze. In the F2 generation, male descendants of LPS-treated fathers showed reduced latency to feed in the novelty-suppressed feeding test. Additionally, the F2 generation of LPS-challenged fathers, but not the F1 generation, displayed enhanced immune response in both sexes after an acute LPS immune challenge (5 mg/kg). Analysis of sperm small noncoding RNA profiles from LPS-treated F0 mice revealed significant changes at 4 weeks after administration of LPS. These changes included three microRNAs, eight PIWI-interacting RNAs, and two transfer RNAs, exhibiting significant upregulation (mmu-miR-146a-5p, mmu-piR-27082 and mmu-piR-29102) or downregulation (mmu-miR-5110, mmu-miR-467e-3p, mmu-piR-22583, mmu-piR-23548, mmu-piR-36341, mmu-piR-50293, mmu-piR-16583, mmu-piR-36507, Mus_musculus_tRNA-Ile-AAT-2-1 and Mus_musculus_tRNA-Tyr-GTA-1-1). Additionally, we detected 52 upregulated small noncoding RNAs (including 9 miRNAs, 41 piRNAs, and 2 tRNAs) and 7 downregulated small noncoding RNAs (3 miRNAs, 3 piRNAs, and 1 tRNA) in the sperm of F1 offspring from LPS-treated males. These findings provide compelling evidence for the involvement of epigenetic mechanisms in the modulation of brain function and immunity, and associated behavioral and immunological traits, across generations, in response to bacterial infection.

piRNA associates with immune diseases

PIWI-interacting RNA (piRNA) is the most abundant small non-coding RNA in animal cells, typically 26-31 nucleotides in length and it binds with PIWI proteins, a subfamily of Argonaute proteins. Initially discovered in germ cells, piRNA is well known for its role in silencing transposons and maintaining genome integrity. However, piRNA is also present in somatic cells as well as in extracellular vesicles and exosomes. While piRNA has been extensively studied in various diseases, particular cancer, its function in immune diseases remains unclear. In this review, we summarize current research on piRNA in immune diseases. We first introduce the basic characteristics, biogenesis and functions of piRNA. Then, we review the association of piRNA with different types of immune diseases, including autoimmune diseases, immunodeficiency diseases, infectious diseases, and other immune-related diseases. piRNA is considered a promising biomarker for diseases, highlighting the need for further research into its potential mechanisms in disease pathogenesis.

Selection and application of small non-coding RNAs for normalizing RT-qPCR data of bovine preimplantation embryo conditioned medium

Small non-coding RNAs (sncRNAs) present in the conditioned medium (CM) of bovine preimplantation embryos are potential noninvasive biomarkers for assessing embryo quality. Accurate quantification of sncRNA levels in the spent CM is of utmost importance in this regard. RT-qPCR is considered as the gold standard for quantifying RNA. In order to standardize RT-qPCR data in the sample type under investigation, the use of suitable stable sncRNAs is essential. Here, we selected 10 sncRNAs from small RNA sequencing of CM samples derived from both bovine blastocysts and degenerate embryos, and evaluated their expression stability together with that of cel-miR-39 as a spike and the often-used U6 small nuclear RNA at different embryo developmental stages. In CM of 2-cell embryos, rsRNA-1044 showed the most stable expression, while tDR-1:32-Gly-CCC-1 was the most stable expressed sncRNA in CM of the stages beyond the 2-cell stage. Next, tDR-1:32-Gly-CCC-1 was used for normalizing the RT-qPCR data from the CM of blastocysts and degenerate embryos. Bta-miR-155 and tDR-39:75-Arg-CCG-2 were found to be significantly up-regulated in the CM of blastocysts compared to that of the degenerated embryos (P = 0.028 and P = 0.017, respectively), suggesting their expression levels are related to embryo development stage. In conclusion, tDR-1:32-Gly-CCC-1 can serve as a suitable reference sncRNA for normalization of RT-qPCR data of the CM from bovine blastocysts.

Oviductal extracellular vesicles miRNA cargo varies in response to embryos and their quality

BACKGROUND

Increasing evidence points to an active role of oviductal extracellular vesicles (oEVs) in the early embryo-maternal dialogue. However, it remains unclear whether oEVs contribute to the recognition of the presence of embryos and their quality in the oviduct. Hence, we examined whether the molecular cargo of oEVs secreted by bovine oviduct epithelial cells (BOEC) differs depending on the presence of good (≥ 8 cells, G) or poor (< 8 cells, P) quality embryos. In addition, differences in RNA profiles between G and P embryos were analyzed in attempt to distinguish oEVs and embryonic EVs cargos.

METHODS

For this purpose, primary BOEC were co-cultured with in vitro produced embryos (IVP) 53 h post fertilization as follows: BOEC with G embryos (BGE); BOEC with P embryos (BPE); G embryos alone (GE); P embryos alone (PE); BOEC alone (B) and medium control (M). After 24 h of co-culture, conditioned media were collected from all groups and EVs were isolated and characterized. MicroRNA profiling of EVs and embryos was performed by small RNA-sequencing.

RESULTS

In EVs, 84 miRNAs were identified, with 8 differentially abundant (DA) miRNAs for BGE vs. B and 4 for BPE vs. B (P-value < 0.01). In embryos, 187 miRNAs were identified, with 12 DA miRNAs for BGE vs. BPE, 3 for G vs. P, 8 for BGE vs. GE, and 11 for BPE vs. PE (P-value < 0.01).

CONCLUSIONS

These results indicated that oEVs are involved in the oviductal-embryo recognition and pointed to specific miRNAs with signaling and supporting roles during early embryo development.

PIWI-Interacting RNAs: A Pivotal Regulator in Neurological Development and Disease

PIWI-interacting RNAs (piRNAs), a class of small non-coding RNAs (sncRNAs) with 24-32 nucleotides (nt), were initially identified in the reproductive system. Unlike microRNAs (miRNAs) or small interfering RNAs (siRNAs), piRNAs normally guide P-element-induced wimpy testis protein (PIWI) families to slice extensively complementary transposon transcripts without the seed pairing. Numerous studies have shown that piRNAs are abundantly expressed in the brain, and many of them are aberrantly regulated in central neural system (CNS) disorders. However, the role of piRNAs in the related developmental and pathological processes is unclear. The elucidation of piRNAs/PIWI would greatly improve the understanding of CNS development and ultimately lead to novel strategies to treat neural diseases. In this review, we summarized the relevant structure, properties, and databases of piRNAs and their functional roles in neural development and degenerative disorders. We hope that future studies of these piRNAs will facilitate the development of RNA-based therapeutics for CNS disorders.

tRNAGlu-derived fragments from embryonic extracellular vesicles modulate bovine embryo hatching

Transfer RNA-derived small RNAs (tsRNAs) have been shown to be involved in early embryo development and repression of endogenous retroelements in embryos and stem cells. However, it is unknown whether tsRNAs also regulate embryo hatching. In this study, we mined the sequencing data of a previous experiment in which we demonstrated that the microRNA (miRNA) cargo of preimplantation embryonic extracellular vesicles (EVs) influences embryo development. We thus profiled the tsRNA cargo of EVs secreted by blastocysts and non-blastocysts. The majority of tsRNAs was identified as tRNA halves originating from the 5´ ends of tRNAs. Among the 148 differentially expressed tsRNAs, the 19 nt tRNA fragment (tRF) tDR-14:32-Glu-CTC-1 was found to be significantly up-regulated in EVs derived from non-blastocysts. RT-qPCR assays confirmed its significant up-regulation in non-blastocyst embryos and their conditioned medium compared to the blastocyst group (P < 0.05). Inhibition of tDR-14:32-Glu-CTC-1 by supplementing antagomirs to the conditioned medium improved embryo hatching (P < 0.05). Transcriptomic analysis of embryos treated with tDR-14:32-Glu-CTC-1 antagomirs further showed differential expression of genes that are associated with embryo hatching and implantation. In summary, tDR-14:32-Glu-CTC-1 is up-regulated in non-blastocyst embryos and their secretions, and inhibition of tDR-14:32-Glu-CTC-1 promotes embryo hatching, while influencing embryo implantation-related genes and pathways. These results indicate that embryonic EVs containing specific tRFs may regulate preimplantation embryo development.

NCAD v1.0: a database for non-coding variant annotation and interpretation

The application of whole genome sequencing is expanding in clinical diagnostics across various genetic disorders, and the significance of non-coding variants in penetrant diseases is increasingly being demonstrated. Therefore, it is urgent to improve the diagnostic yield by exploring the pathogenic mechanisms of variants in non-coding regions. However, the interpretation of non-coding variants remains a significant challenge, due to the complex functional regulatory mechanisms of non-coding regions and the current limitations of available databases and tools. Hence, we develop the non-coding variant annotation database (NCAD, http://www.ncawdb.net/), encompassing comprehensive insights into 665,679,194 variants, regulatory elements, and element interaction details. Integrating data from 96 sources, spanning both GRCh37 and GRCh38 versions, NCAD v1.0 provides vital information to support the genetic diagnosis of non-coding variants, including allele frequencies of 12 diverse populations, with a particular focus on the population frequency information for 230,235,698 variants in 20,964 Chinese individuals. Moreover, it offers prediction scores for variant functionality, five categories of regulatory elements, and four types of non-coding RNAs. With its rich data and comprehensive coverage, NCAD serves as a valuable platform, empowering researchers and clinicians with profound insights into non-coding regulatory mechanisms while facilitating the interpretation of non-coding variants.

sRNAflow: A Tool for the Analysis of Small RNA-Seq Data

The analysis of small RNA sequencing data across a range of biofluids is a significant research area, given the diversity of RNA types that hold potential diagnostic, prognostic, and predictive value. The intricate task of segregating the complex mixture of small RNAs from both human and other species, including bacteria, fungi, and viruses, poses one of the most formidable challenges in the analysis of small RNA sequencing data, currently lacking satisfactory solutions. This study introduces sRNAflow, a user-friendly bioinformatic tool with a web interface designed for the analysis of small RNAs obtained from biological fluids. Tailored to the unique requirements of such samples, the proposed pipeline addresses various challenges, including filtering potential RNAs from reagents and environment, classifying small RNA types, managing small RNA annotation overlap, conducting differential expression assays, analysing isomiRs, and presenting an approach to identify the sources of small RNAs within samples. sRNAflow also encompasses an alternative alignment-free analysis of RNA-seq data, featuring clustering and initial RNA source identification using BLAST. This comprehensive approach facilitates meaningful comparisons of results between different analytical methods.

piOxi database: a web resource of germline and somatic tissue piRNAs identified by chemical oxidation

PIWI-interacting RNAs (piRNAs) are a class of small non-coding RNAs that are highly expressed and extensively studied from the germline. piRNAs associate with PIWI proteins to maintain DNA methylation for transposon silencing and transcriptional gene regulation for genomic stability. Mature germline piRNAs have distinct characteristics including a 24- to 32-nucleotide length and a 2'-O-methylation signature at the 3' end. Although recent studies have identified piRNAs in somatic tissues, they remain poorly characterized. For example, we recently demonstrated notable expression of piRNA in the murine soma, and while overall expression was lower than that of the germline, unique characteristics suggested tissue-specific functions of this class. While currently available databases commonly use length and association with PIWI proteins to identify piRNA, few have included a chemical oxidation method that detects piRNA based on its 3' modification. This method leads to reproducible and rigorous data processing when coupled with next-generation sequencing and bioinformatics analysis. Here, we introduce piOxi DB, a user-friendly web resource that provides a comprehensive analysis of piRNA, generated exclusively through sodium periodate treatment of small RNA. The current version of piOxi DB includes 435 749 germline and 9828 somatic piRNA sequences robustly identified from M. musculus, M. fascicularis and H. sapiens. The database provides species- and tissue-specific data that are further analyzed according to chromosome location and correspondence to gene and repetitive elements. piOxi DB is an informative tool to assist broad research applications in the fields of RNA biology, cancer biology, environmental toxicology and beyond. Database URL:  https://pioxidb.dcmb.med.umich.edu/.

Paternal immune activation by Poly I:C modulates sperm noncoding RNA profiles and causes transgenerational changes in offspring behavior

Paternal pre-conceptual environmental experiences, such as stress and diet, can affect offspring brain and behavioral phenotypes via epigenetic modifications in sperm. Furthermore, maternal immune activation due to infection during gestation can reprogram offspring behavior and brain functioning in adulthood. However, the effects of paternal pre-conceptual exposure to immune activation on the behavior and physiology of offspring (F1) and grand-offspring (F2) are not currently known. We explored effects of paternal pre-conceptual exposure to viral-like immune activation on F1 and F2 behavioral and physiological phenotypes using a C57BL/6J mouse model. Males were treated with a single injection (intraperitoneal) of the viral mimetic polyinosinic:polycytidylic acid (Poly I:C: 12 mg/kg) then bred with naïve female mice four weeks after the Poly I:C (or 0.9% saline control) injection. The F1 offspring of Poly I:C treated fathers displayed increased depression-like behavior in the Porsolt swim test, an altered stress response in the novelty-suppressed feeding test, and significant transcriptomic changes in their hippocampus. Additionally, the F1 male offspring of Poly I:C treated F0 males showed significantly increased immune responsivity after a Poly I:C immune challenge (12 mg/kg). Furthermore, the F2 male grand-offspring took longer to enter and travelled significantly shorter distances in the light zone of the light/dark box. An analysis of the small noncoding RNA profiles in sperm from Poly I:C treated males and their male offspring revealed significant effects of Poly I:C on the sperm microRNA content at the time of conception and on the sperm PIWI-interacting RNA content of the male offspring. Notably, eight miRNAs with an FDR < 0.05 (miR-141-3p, miR-126b-5p, miR-669o-5p, miR-10b-3p, miR-471-5p, miR-463-5p, miR-148b-3p, and miR-181c-5p) were found to be significantly downregulated in the sperm of Poly I:C treated males. Collectively, we demonstrate that paternal pre-conceptual exposure to a viral immune challenge results in both intergenerational and transgenerational effects on brain and behavior that may be mediated by alterations in the sperm small noncoding RNA content.

Novel roles of PIWI proteins and PIWI-interacting RNAs in human health and diseases

Non-coding RNA has aroused great research interest recently, they play a wide range of biological functions, such as regulating cell cycle, cell proliferation, and intracellular substance metabolism. Piwi-interacting RNAs (piRNAs) are emerging small non-coding RNAs that are 24-31 nucleotides in length. Previous studies on piRNAs were mainly limited to evaluating the binding to the PIWI protein family to play the biological role. However, recent studies have shed more lights on piRNA functions; aberrant piRNAs play unique roles in many human diseases, including diverse lethal cancers. Therefore, understanding the mechanism of piRNAs expression and the specific functional roles of piRNAs in human diseases is crucial for developing its clinical applications. Presently, research on piRNAs mainly focuses on their cancer-specific functions but lacks investigation of their expressions and epigenetic modifications. This review discusses piRNA's biogenesis and functional roles and the recent progress of functions of piRNA/PIWI protein complexes in human diseases. Video Abstract.

LSTM4piRNA: Efficient piRNA Detection in Large-Scale Genome Databases Using a Deep Learning-Based LSTM Network

Piwi-interacting RNAs (piRNAs) are a new class of small, non-coding RNAs, crucial in the regulation of gene expression. Recent research has revealed links between piRNAs, viral defense mechanisms, and certain human cancers. Due to their clinical potential, there is a great interest in identifying piRNAs from large genome databases through efficient computational methods. However, piRNAs lack conserved structure and sequence homology across species, which makes piRNA detection challenging. Current detection algorithms heavily rely on manually crafted features, which may overlook or improperly use certain features. Furthermore, there is a lack of suitable computational tools for analyzing large-scale databases and accurately identifying piRNAs. To address these issues, we propose LSTM4piRNA, a highly efficient deep learning-based method for predicting piRNAs in large-scale genome databases. LSTM4piRNA utilizes a compact LSTM network that can effectively analyze RNA sequences from extensive datasets to detect piRNAs. It can automatically learn the dependencies among RNA sequences, and regularization is further integrated to reduce the generalization error. Comprehensive performance evaluations based on piRNAs from the piRBase database demonstrate that LSTM4piRNA outperforms current advanced methods and is well-suited for analysis with large-scale databases.

Detection of a 7SL RNA-derived small non-coding RNA using Molecular Beacons in vitro and in cells

Small non-coding RNAs (sncRNA) are involved in many steps of the gene expression cascade and regulate processing and expression of mRNAs by the formation of ribonucleoprotein complexes (RNP) such as the RNA-induced silencing complex (RISC). By analyzing small RNA Seq data sets, we identified a sncRNA annotated as piR-hsa-1254, which is likely derived from the 3'-end of 7SL RNA2 (RN7SL2), herein referred to as snc7SL RNA. The 7SL RNA is an abundant long non-coding RNA polymerase III transcript and serves as structural component of the cytoplasmic signal recognition particle (SRP). To evaluate a potential functional role of snc7SL RNA, we aimed to define its cellular localization by live cell imaging. Therefore, a Molecular Beacon (MB)-based method was established to compare the subcellular localization of snc7SL RNA with its precursor 7SL RNA. We designed and characterized several MBs in vitro and tested those by live cell fluorescence microscopy. Using a multiplex approach, we show that 7SL RNA localizes mainly to the endoplasmic reticulum (ER), as expected for the SRP, whereas snc7SL RNA predominately localizes to the nucleus. This finding suggests a fundamentally different function of 7SL RNA and its derivate snc7SL RNA.

Pan-cancer landscape of immunology PIWI-interacting RNAs

PIWI-interacting RNAs (piRNAs), an emergent type of non-coding RNAs during oncogenesis, play critical roles in regulating tumor microenvironment. Systematic analysis of piRNAs' roles in modulating immune pathways is important for tumor immunotherapy. In this study, in-depth analysis of piRNAs was performed to develop an integrated computational algorithm, the immunology piRNA (ImmPI) pipeline, for uncovering the global expression landscape of piRNAs and identifying their regulatory roles in immune pathways. The immunology piRNAs show a tendency towards overexpression patterns in immune cells, causing perturbations in tumors, being significantly associated with infiltration of immune cells, and having prognostic value. The ImmPI score can contribute to prioritizing tumor-related piRNAs and distinguish two subtypes of SKCM (immune-cold and hot phenotypes), as characterized by different prognoses, immunogenicity and antitumor immunity. Finally, we developed an interactive web resource (ImmPI portal: http://www.hbpding.com/ImmPi) for the biomedical research community, with several useful modules to browse, visualize, and download the results of immunology piRNAs analysis. Overall, our work provides a comprehensive landscape of piRNAs across multiple cancer types and sheds light on their regulatory and functional roles in tumor immunity. These findings pave the way for future research and development of piRNA-based immunotherapies for cancer treatment.

Interactions of piRNAs with the mRNA of Candidate Genes in Esophageal Squamous Cell Carcinoma

Recently, a database of human piRNAs (piwi-interacting RNAs) was created, which allows the study of the binding of many piRNAs to the mRNAs of genes involved in many diseases, including cancer. In the present work, we identified the piRNAs that can interact with candidate esophageal squamous cell carcinoma (ESCC) genes. The binding of 480 thousand piRNAs with the mRNAs of 66 candidate ESCC genes was studied. Bioinformatic studies found that piRNAs bind only to the mRNAs of nine candidate genes: AURKA, BMP7, GCOM1, ERCC1, MTHFR, SASH1, SIX4, SULT1A1, and TP53. It has been shown that piRNAs can bind to mRNA by overlapping nucleotide sequences in limited 3'UTR and 5'UTR regions called clusters of binding sites (BSs). The existence of clusters of piRNA BSs significantly reduces the proportion of the nucleotide sequences of these sites in the mRNA of target genes. Competition between piRNAs occurs for binding to the mRNA of target genes. Individual piRNAs and groups of piRNAs that have separate BSs and clusters of BSs in the mRNAs of two or more candidate genes have been identified in the mRNAs of these genes. This organization of piRNAs BSs indicates the interdependence of the expression of candidate genes through piRNAs. Significant differences in the ability of genes to interact with piRNAs prevent the side effects of piRNAs on genes with a lack of the ability to bind such piRNAs. Individual piRNAs and sets of piRNAs are proposed and recommended for the diagnosis and therapy of ESCC.

iPiDA-SWGCN: Identification of piRNA-disease associations based on Supplementarily Weighted Graph Convolutional Network

Accurately identifying potential piRNA-disease associations is of great importance in uncovering the pathogenesis of diseases. Recently, several machine-learning-based methods have been proposed for piRNA-disease association detection. However, they are suffering from the high sparsity of piRNA-disease association network and the Boolean representation of piRNA-disease associations ignoring the confidence coefficients. In this study, we propose a supplementarily weighted strategy to solve these disadvantages. Combined with Graph Convolutional Networks (GCNs), a novel predictor called iPiDA-SWGCN is proposed for piRNA-disease association prediction. There are three main contributions of iPiDA-SWGCN: (i) Potential piRNA-disease associations are preliminarily supplemented in the sparse piRNA-disease network by integrating various basic predictors to enrich network structure information. (ii) The original Boolean piRNA-disease associations are assigned with different relevance confidence to learn node representations from neighbour nodes in varying degrees. (iii) The experimental results show that iPiDA-SWGCN achieves the best performance compared with the other state-of-the-art methods, and can predict new piRNA-disease associations.

A systematic review of non-coding RNA genes with differential expression profiles associated with autism spectrum disorders

AIMS

To identify differential expression of shorter non-coding RNA (ncRNA) genes associated with autism spectrum disorders (ASD).

BACKGROUND

ncRNA are functional molecules that derive from non-translated DNA sequence. The HUGO Gene Nomenclature Committee (HGNC) have approved ncRNA gene classes with alignment to the reference human genome. One subset is microRNA (miRNA), which are highly conserved, short RNA molecules that regulate gene expression by direct post-transcriptional repression of messenger RNA. Several miRNA genes are implicated in the development and regulation of the nervous system. Expression of miRNA genes in ASD cohorts have been examined by multiple research groups. Other shorter classes of ncRNA have been examined less. A comprehensive systematic review examining expression of shorter ncRNA gene classes in ASD is timely to inform the direction of research.

METHODS

We extracted data from studies examining ncRNA gene expression in ASD compared with non-ASD controls. We included studies on miRNA, piwi-interacting RNA (piRNA), small NF90 (ILF3) associated RNA (snaR), small nuclear RNA (snRNA), small nucleolar RNA (snoRNA), transfer RNA (tRNA), vault RNA (vtRNA) and Y RNA. The following electronic databases were searched: Cochrane Library, EMBASE, PubMed, Web of Science, PsycINFO, ERIC, AMED and CINAHL for papers published from January 2000 to May 2022. Studies were screened by two independent investigators with a third resolving discrepancies. Data was extracted from eligible papers.

RESULTS

Forty-eight eligible studies were included in our systematic review with the majority examining miRNA gene expression alone. Sixty-four miRNA genes had differential expression in ASD compared to controls as reported in two or more studies, but often in opposing directions. Four miRNA genes had differential expression in the same direction in the same tissue type in at least 3 separate studies. Increased expression was reported in miR-106b-5p, miR-155-5p and miR-146a-5p in blood, post-mortem brain, and across several tissue types, respectively. Decreased expression was reported in miR-328-3p in bloods samples. Seven studies examined differential expression from other classes of ncRNA, including piRNA, snRNA, snoRNA and Y RNA. No individual ncRNA genes were reported in more than one study. Six studies reported differentially expressed snoRNA genes in ASD. A meta-analysis was not possible because of inconsistent methodologies, disparate tissue types examined, and varying forms of data presented.

CONCLUSION

There is limited but promising evidence associating the expression of certain miRNA genes and ASD, although the studies are of variable methodological quality and the results are largely inconsistent. There is emerging evidence associating differential expression of snoRNA genes in ASD. It is not currently possible to say whether the reports of differential expression in ncRNA may relate to ASD aetiology, a response to shared environmental factors linked to ASD such as sleep and nutrition, other molecular functions, human diversity, or chance findings. To improve our understanding of any potential association, we recommend improved and standardised methodologies and reporting of raw data. Further high-quality research is required to shine a light on possible associations, which may yet yield important information.

A ligation-independent sequencing method reveals tRNA-derived RNAs with blocked 3' termini

Despite the numerous sequencing methods available, the vast diversity in size and chemical modifications of RNA molecules makes the capture of the full spectrum of cellular RNAs a difficult task. By combining quasi-random hexamer priming with a custom template switching strategy, we developed a method to construct sequencing libraries from RNA molecules of any length and with any type of 3' terminal modification, allowing the sequencing and analysis of virtually all RNA species. Ligation-independent detection of all types of RNA (LIDAR) is a simple, effective tool to comprehensively characterize changes in small non-coding RNAs and mRNAs simultaneously, with performance comparable to separate dedicated methods. With LIDAR, we comprehensively characterized the coding and non-coding transcriptome of mouse embryonic stem cells, neural progenitor cells, and sperm. LIDAR detected a much larger variety of tRNA-derived RNAs (tDRs) compared to traditional ligation-dependent sequencing methods, and uncovered the presence of tDRs with blocked 3' ends that had previously escaped detection. Our findings highlight the potential of LIDAR to systematically detect all RNAs in a sample and uncover new RNA species with potential regulatory functions.

IntiCom-DB: A Manually Curated Database of Inter-Tissue Communication Molecules and Their Communication Routes

Inter-tissue communication (ITC) is critical for maintaining the physiological functions of multiple tissues and is closely related to the onset and development of various complex diseases. Nevertheless, there is no well-organized data resource for known ITC molecules with explicit ITC routes from source tissues to target tissues. To address this issue, in this work, we manually reviewed nearly 190,000 publications and identified 1408 experimentally supported ITC entries in which the ITC molecules, their communication routes, and their functional annotations were included. To facilitate our work, these curated ITC entries were incorporated into a user-friendly database named IntiCom-DB. This database also enables visualization of the expression abundances of ITC proteins and their interaction partners. Finally, bioinformatics analyses on these data revealed common biological characteristics of the ITC molecules. For example, tissue specificity scores of ITC molecules at the protein level are often higher than those at the mRNA level in the target tissues. Moreover, the ITC molecules and their interaction partners are more abundant in both the source tissues and the target tissues. IntiCom-DB is freely available as an online database. As the first comprehensive database of ITC molecules with explicit ITC routes to the best of our knowledge, we hope that IntiCom-DB will benefit future ITC-related studies.

Endogenous piRNAs Can Interact with the Omicron Variant of the SARS-CoV-2 Genome

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which caused the COVID-19 pandemic, can still infect populations in many countries around the globe. The Omicron strain is the most mutated variant of SARS-CoV-2. The high transmissibility of the strain and its ability to evade immunity necessitate a priority study of its properties in order to quickly create effective means of preventing its spread. The current research aimed to examine the in silico interaction between PIWI-interacting RNAs (piRNAs) and the SARS-CoV-2 genome (gRNA) to identify endogenous piRNAs and propose synthetic piRNAs with strong antiviral activity for drug development. This study used validated bioinformatic approaches regarding the interaction of more than eight million piRNAs with the SARS-CoV-2 genome. The piRNAs’ binding sites (BSs) in the 5′UTR were located with overlapping nucleotide sequences termed clusters of BSs. Several BSs clusters have been found in the nsp3, nsp7, RNA-dependent RNA polymerase, endoRNAse, S surface glycoprotein, ORF7a, and nucleocapsid. Sixteen synthetic piRNAs that interact with gRNA have been proposed with free binding energy ranging from −170 kJ/mol to −175 kJ/mol, which can be used to create drugs that suppress the reproduction of SARS-CoV-2.

Role of piRNA biogenesis and its neuronal function in the development of neurodegenerative diseases

Neurodegenerative diseases, such as Alzheimer's disease (AD), Parkinson's disease (PD), and amyotrophic lateral sclerosis (ALS), are caused by neuronal loss and dysfunction. Despite remarkable improvements in our understanding of these pathogeneses, serious worldwide problems with significant public health burdens are remained. Therefore, new efficient diagnostic and therapeutic strategies are urgently required. PIWI-interacting RNAs (piRNAs) are a major class of small non-coding RNAs that silence gene expression through transcriptional and post-transcriptional processes. Recent studies have demonstrated that piRNAs, originally found in the germ line, are also produced in non-gonadal somatic cells, including neurons, and further revealed the emerging roles of piRNAs, including their roles in neurodevelopment, aging, and neurodegenerative diseases. In this review, we aimed to summarize the current knowledge regarding the piRNA roles in the pathophysiology of neurodegenerative diseases. In this context, we first reviewed on recent updates on neuronal piRNA functions, including biogenesis, axon regeneration, behavior, and memory formation, in humans and mice. We also discuss the aberrant expression and dysregulation of neuronal piRNAs in neurodegenerative diseases, such as AD, PD, and ALS. Moreover, we review pioneering preclinical studies on piRNAs as biomarkers and therapeutic targets. Elucidation of the mechanisms underlying piRNA biogenesis and their functions in the brain would provide new perspectives for the clinical diagnosis and treatment of AD and various neurodegenerative diseases.

Characteristics of piRNAs and their comparative profiling in testes of sheep with different fertility

As a novel class of small RNAs, piRNAs are highly expressed in the animal gonads and their main known role is to inhibit transposon activity for ensuring the correctness and integrity of genome. In order to explore the characteristics of piRNAs in sheep testis and their possible regulatory roles on male reproduction, deep sequencing technology was used to sequence small RNAs and identify piRNAs in testes of sheep. The length of piRNAs in sheep testes showed a unimodal distribution between 26 and 31 nt, with a peak at 29 nt. These piRNAs exhibited obvious ping-pong signature and strand specificity. In the genome, they were mainly aligned to CDS, intron, repetitive sequence regions and unannotated regions. Furthermore, in transposon analysis, piRNAs were aligned predominantly to LINE, SINE, and LTR types of retrotransposon in sheep testes, and the piRNAs derived from each type showed obvious ping-pong signature. The piRNA clusters identified in sheep testes were mainly distributed on chromosomes 3, 7, 15, 17, 18 and 20. The results combining semen determination with pathway enrichment analysis implied that differentially expressed piRNAs between the testes of rams with different fertility might participate in spermatogenesis by regulating multiple pathways closely related to stabilization of blood-testis barrier and renewal and differentiation of spermatogonial stem cell. Taken together, the study provided new insights into the characteristics, origin and expression patterns of piRNAs in sheep testes tissue, which would help us better understand the role of piRNAs in sheep reproduction.

piRNAs may regulate expression of candidate genes of esophageal adenocarcinoma

Elucidation of ways to regulate the expression of candidate cancer genes will contribute to the development of methods for cancer diagnosis and therapy. The aim of the present study was to show the role of piRNAs as efficient regulators of mRNA translation of esophageal adenocarcinoma (EAC) candidate genes. We used bioinformatic methods to study the interaction characteristics of up to 200 thousand piRNAs with mRNAs of 38 candidate EAC genes. The piRNAs capable of binding to mRNA of AR, BTG3, CD55, ERBB3, FKBP5, FOXP1, LEP, SEPP1, SMAD4, and TP53 genes with high free energy by the formation of hydrogen bonds between canonical (G-C, A-U) and noncanonical (G-U, A-C) piRNA and mRNA nucleotide pairs were revealed. The organization of piRNA binding sites (BSs) in the mRNA of candidate genes was found to overlap nucleotide sequences to form clusters. Clusters of piRNA BSs were detected in the 5'-untranslated region, coding domain sequence, and 3'-untranslated region of mRNA. Due to the formation of piRNA binding site clusters, compaction of BSs occurs and competition between piRNAs for binding to mRNA of candidate EAC genes occurs. Associations of piRNA and candidate genes were selected for use as markers for the diagnosis of EAC.

NPInter v5.0: ncRNA interaction database in a new era

Noncoding RNAs (ncRNAs) play key regulatory roles in biological processes by interacting with other biomolecules. With the development of high-throughput sequencing and experimental technologies, extensive ncRNA interactions have been accumulated. Therefore, we updated the NPInter database to a fifth version to document these interactions. ncRNA interaction entries were doubled from 1 100 618 to 2 596 695 by manual literature mining and high-throughput data processing. We integrated global RNA-DNA interactions from iMARGI, ChAR-seq and GRID-seq, greatly expanding the number of RNA-DNA interactions (from 888 915 to 8 329 382). In addition, we collected different types of RNA interaction between SARS-CoV-2 virus and its host from recently published studies. Long noncoding RNA (lncRNA) expression specificity in different cell types from tumor single cell RNA-seq (scRNA-seq) data were also integrated to provide a cell-type level view of interactions. A new module named RBP was built to display the interactions of RNA-binding proteins with annotations of localization, binding domains and functions. In conclusion, NPInter v5.0 (http://bigdata.ibp.ac.cn/npinter5/) provides informative and valuable ncRNA interactions for biological researchers.

iPiDA-GCN: Identification of piRNA-disease associations based on Graph Convolutional Network

Motivation

Piwi-interacting RNAs (piRNAs) play a critical role in the progression of various diseases. Accurately identifying the associations between piRNAs and diseases is important for diagnosing and prognosticating diseases. Although some computational methods have been proposed to detect piRNA-disease associations, it is challenging for these methods to effectively capture nonlinear and complex relationships between piRNAs and diseases because of the limited training data and insufficient association representation.

Results

With the growth of piRNA-disease association data, it is possible to design a more complex machine learning method to solve this problem. In this study, we propose a computational method called iPiDA-GCN for piRNA-disease association identification based on graph convolutional networks (GCNs). The iPiDA-GCN predictor constructs the graphs based on piRNA sequence information, disease semantic information and known piRNA-disease associations. Two GCNs (Asso-GCN and Sim-GCN) are used to extract the features of both piRNAs and diseases by capturing the association patterns from piRNA-disease interaction network and two similarity networks. GCNs can capture complex network structure information from these networks, and learn discriminative features. Finally, the full connection networks and inner production are utilized as the output module to predict piRNA-disease association scores. Experimental results demonstrate that iPiDA-GCN achieves better performance than the other state-of-the-art methods, benefitted from the discriminative features extracted by Asso-GCN and Sim-GCN. The iPiDA-GCN predictor is able to detect new piRNA-disease associations to reveal the potential pathogenesis at the RNA level. The data and source code are available at http://bliulab.net/iPiDA-GCN/.

First characterization of PIWI-interacting RNA clusters in a cichlid fish with a B chromosome

BACKGROUND

B chromosomes are extra elements found in several eukaryote species. Usually, they do not express a phenotype in the host. However, advances in bioinformatics over the last decades have allowed us to describe several genes and molecular functions related to B chromosomes. These advances enable investigations of the relationship between the B chromosome and the host to understand how this element has been preserved in genomes. However, considering that transposable elements (TEs) are highly abundant in this supernumerary chromosome, there is a lack of knowledge concerning the dynamics of TE control in B-carrying cells. Thus, the present study characterized PIWI-interacting RNA (piRNA) clusters and pathways responsible for silencing the mobilization of TEs in gonads of the cichlid fish Astatotilapia latifasciata carrying the B chromosome.

RESULTS

Through small RNA-seq and genome assembly, we predicted and annotated piRNA clusters in the A. latifasciata genome for the first time. We observed that these clusters had biased expression related to sex and the presence of the B chromosome. Furthermore, three piRNA clusters, named curupira, were identified in the B chromosome. Two of them were expressed exclusively in gonads of samples with the B chromosome. The composition of these curupira sequences was derived from LTR, LINE, and DNA elements, representing old and recent transposition events in the A. latifasciata genome and the B chromosome. The presence of the B chromosome also affected the expression of piRNA pathway genes. The mitochondrial cardiolipin hydrolase-like (pld6) gene is present in the B chromosome, as previously reported, and an increase in its expression was detected in gonads with the B chromosome.

CONCLUSIONS

Due to the high abundance of TEs in the B chromosome, it was possible to investigate the origin of piRNA from these jumping genes. We hypothesize that the B chromosome has evolved its own genomic guardians to prevent uncontrolled TE mobilization. Furthermore, we also detected an expression bias in the presence of the B chromosome over A. latifasciata piRNA clusters and pathway genes.

In Silico Study of piRNA Interactions with the SARS-CoV-2 Genome

A prolonged pandemic with numerous human casualties requires a rapid search for means to control the various strains of SARS-CoV-2. Since only part of the human population is affected by coronaviruses, there are probably endogenous compounds preventing the spread of these viral pathogens. It has been shown that piRNA (PIWI-interacting RNAs) interact with the mRNA of human genes and can block protein synthesis at the stage of translation. Estimated the effects of piRNA on SARS-CoV-2 genomic RNA (gRNA) in silico. A cluster of 13 piRNA binding sites (BS) in the SARS-CoV-2 gRNA region encoding the oligopeptide was identified. The second cluster of BSs 39 piRNAs also encodes the oligopeptide. The third cluster of 24 piRNA BS encodes the oligopeptide. Twelve piRNAs were identified that strongly interact with the gRNA. Based on the identified functionally important endogenous piRNAs, synthetic piRNAs (spiRNAs) are proposed that will suppress the multiplication of the coronavirus even more strongly. These spiRNAs and selected endogenous piRNAs have little effect on human 17494 protein-coding genes, indicating a low probability of side effects. The piRNA and spiRNA selection methodology created for the control of SARS-CoV-2 (NC_045512.2) can be used to control all strains of SARS-CoV-2.

PIWI-interacting RNAs in human diseases: databases and computational models

PIWI-interacting RNAs (piRNAs) are short 21-35 nucleotide molecules that comprise the largest class of non-coding RNAs and found in a large diversity of species including yeast, worms, flies, plants and mammals including humans. The most well-understood function of piRNAs is to monitor and protect the genome from transposons particularly in germline cells. Recent data suggest that piRNAs may have additional functions in somatic cells although they are expressed there in far lower abundance. Compared with microRNAs (miRNAs), piRNAs have more limited bioinformatics resources available. This review collates 39 piRNA specific and non-specific databases and bioinformatics resources, describes and compares their utility and attributes and provides an overview of their place in the field. In addition, we review 33 computational models based upon function: piRNA prediction, transposon element and mRNA-related piRNA prediction, cluster prediction, signature detection, target prediction and disease association. Based on the collection of databases and computational models, we identify trends and potential gaps in tool development. We further analyze the breadth and depth of piRNA data available in public sources, their contribution to specific human diseases, particularly in cancer and neurodegenerative conditions, and highlight a few specific piRNAs that appear to be associated with these diseases. This briefing presents the most recent and comprehensive mapping of piRNA bioinformatics resources including databases, models and tools for disease associations to date. Such a mapping should facilitate and stimulate further research on piRNAs.

HIV-1 Tat and cocaine coexposure impacts piRNAs to affect astrocyte energy metabolism

Aim: To understand the effect of HIV infection and cocaine exposure on piRNA expression in human primary astrocytes. Materials & methods: We used small RNA sequencing analysis to investigate the impacts of HIV-1 Tat and cocaine coexposure on the expression of piRNAs in human primary astrocytes. Results: We identified 27,700 piRNAs and analyzed them by small RNA next-generation sequencing. A total of 239 piRNAs were significantly altered by HIV-1 Tat and cocaine coexposure. We also identified PIWIL1, PIWIL2, PIWIL3 and PIWIL4 as interacting partners of piRNAs that were affected by cocaine and HIV-1 Tat coexposure. Epigenetic changes in the expression levels of these piRNA targets were associated with Kyoto Encyclopedia of Genes and Genomes pathways of energy metabolism and neurodegeneration. Conclusion: These findings provide evidence that cocaine exposure and HIV infection affect the expression levels of piRNA, PIWIL1, PIWIL2, PIWIL3 and PIWIL4.

PIWI-interacting RNAs and human testicular function

Small noncoding RNAs (sncRNAs) are pieces of RNA with a length below 200 bp and represent a diverse group of RNAs having many different biological functions. The best described subtype is the microRNAs which primarily function in posttranscriptional gene regulation and appear essential for most physiological processes. Of particular interest for the germline is the PIWI-interacting RNAs (piRNAs) which are a class of sncRNA of 21-35 bp in length that are almost exclusively found in germ cells. Recently, it has become clear that piRNAs are essential for testicular function, and in this perspective, we outline the current knowledge of piRNAs in humans. Although piRNAs appear unique to germ cells, they have also been described in various somatic cancers and biofluids. Here, we discuss the potential function of piRNAs in somatic tissues and whether detection in biofluids may be used as a biomarker for testicular function. This article is categorized under: Reproductive System Diseases > Genetics/Genomics/Epigenetics Reproductive System Diseases > Molecular and Cellular Physiology.