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.

miR-9a regulates levels of both rhomboid mRNA and protein in the early Drosophila melanogaster embryo.

MicroRNAs can have subtle and combinatorial effects on the levels of the targets and pathways they act on. Studying the consequences of a single microRNA knockout often proves difficult as many such knockouts exhibit phenotypes only under stress conditions. This has often led to the hypothesis that microRNAs buffer the effects of intrinsic and environmental stochasticity on gene expression. Observing and understanding this buffering effect entails quantitative analysis of microRNA and target expression in single cells. To this end, we have employed single-molecule fluorescence in situ hybridization, immunofluorescence, and high-resolution confocal microscopy to investigate the effects of miR-9a loss on the expression of the serine-protease Rhomboid in Drosophila melanogaster early embryos. Our single-cell quantitative approach shows that spatially, the rhomboid mRNA pattern is identical in WT and miR-9a knockout embryos. However, we find that the number of mRNA molecules per cell is higher when miR-9a is absent, and the level and temporal accumulation of rhomboid protein shows a more dramatic increase in the miR-9a knockout. Specifically, we see accumulation of rhomboid protein in miR-9a mutants by stage 5, much earlier than in WT. The data, therefore, show that miR-9a functions in the regulation of rhomboid mRNA and protein levels. While further work is required to establish whether rhomboid is a direct target of miR-9 in Drosophila, our results further establish the miR-9 family microRNAs as conserved regulators of timing in neurogenic processes. This study shows the power of single-cell quantification as an experimental tool to study phenotypic consequences of microRNA mis-regulation.

The embryonic transcriptome of Parhyale hawaiensis reveals different dynamics of microRNAs and mRNAs during the maternal-zygotic transition.

Parhyale hawaiensis has emerged as the crustacean model of choice due to its tractability, ease of imaging, sequenced genome, and development of CRISPR/Cas9 genome editing tools. However, transcriptomic datasets spanning embryonic development are lacking, and there is almost no annotation of non-protein-coding RNAs, including microRNAs. We have sequenced microRNAs, together with mRNAs and long non-coding RNAs, in Parhyale using paired size-selected RNA-seq libraries at seven time-points covering important transitions in embryonic development. Focussing on microRNAs, we annotate 175 loci in Parhyale, 88 of which have no known homologs. We use these data to annotate the microRNAome of 37 crustacean genomes, and suggest a core crustacean microRNA set of around 61 sequence families. We examine the dynamic expression of microRNAs and mRNAs during the maternal-zygotic transition. Our data suggest that zygotic genome activation occurs in two waves in Parhyale with microRNAs transcribed almost exclusively in the second wave. Contrary to findings in other arthropods, we do not predict a general role for microRNAs in clearing maternal transcripts. These data significantly expand the available transcriptomics resources for Parhyale, and facilitate its use as a model organism for the study of small RNAs in processes ranging from embryonic development to regeneration.

Single-cell visualization of mir-9a and Senseless co-expression during Drosophila melanogaster embryonic and larval peripheral nervous system development.

The Drosophila melanogaster peripheral nervous system (PNS) comprises the sensory organs that allow the fly to detect environmental factors such as temperature and pressure. PNS development is a highly specified process where each sensilla originates from a single sensory organ precursor (SOP) cell. One of the major genetic orchestrators of PNS development is Senseless, which encodes a zinc finger transcription factor (Sens). Sens is both necessary and sufficient for SOP differentiation. Senseless expression and SOP number are regulated by the microRNA miR-9a. However, the reciprocal dynamics of Senseless and miR-9a are still obscure. By coupling single-molecule FISH with immunofluorescence, we are able to visualize transcription of the mir-9a locus and expression of Sens simultaneously. During embryogenesis, we show that the expression of mir-9a in SOP cells is rapidly lost as Senseless expression increases. However, this mutually exclusive expression pattern is not observed in the third instar imaginal wing disk, where some Senseless-expressing cells show active sites of mir-9a transcription. These data challenge and extend previous models of Senseless regulation and show complex co-expression dynamics between mir-9a and Senseless. The differences in this dynamic relationship between embryonic and larval PNS development suggest a possible switch in miR-9a function. Our work brings single-cell resolution to the understanding of dynamic regulation of PNS development by Senseless and miR-9a.

Rfam 14: expanded coverage of metagenomic, viral and microRNA families.

Rfam is a database of RNA families where each of the 3444 families is represented by a multiple sequence alignment of known RNA sequences and a covariance model that can be used to search for additional members of the family. Recent developments have involved expert collaborations to improve the quality and coverage of Rfam data, focusing on microRNAs, viral and bacterial RNAs. We have completed the first phase of synchronising microRNA families in Rfam and miRBase, creating 356 new Rfam families and updating 40. We established a procedure for comprehensive annotation of viral RNA families starting with Flavivirus and Coronaviridae RNAs. We have also increased the coverage of bacterial and metagenome-based RNA families from the ZWD database. These developments have enabled a significant growth of the database, with the addition of 759 new families in Rfam 14. To facilitate further community contribution to Rfam, expert users are now able to build and submit new families using the newly developed Rfam Cloud family curation system. New Rfam website features include a new sequence similarity search powered by RNAcentral, as well as search and visualisation of families with pseudoknots. Rfam is freely available at https://rfam.org.

Regulatory RNAs: A Universal Language for Inter-Domain Communication.

In eukaryotes, microRNAs (miRNAs) have roles in development, homeostasis, disease and the immune response. Recent work has shown that plant and mammalian miRNAs also mediate cross-kingdom and cross-domain communications. However, these studies remain controversial and are lacking critical mechanistic explanations. Bacteria do not produce miRNAs themselves, and therefore it is unclear how these eukaryotic RNA molecules could function in the bacterial recipient. In this review, we compare and contrast the biogenesis and functions of regulatory RNAs in eukaryotes and bacteria. As a result, we discovered several conserved features and homologous components in these distinct pathways. These findings enabled us to propose novel mechanisms to explain how eukaryotic miRNAs could function in bacteria. Further understanding in this area is necessary to validate the findings of existing studies and could facilitate the use of miRNAs as novel tools for the directed remodelling of the human microbiota.

Dynamical gene regulatory networks are tuned by transcriptional autoregulation with microRNA feedback.

Concepts from dynamical systems theory, including multi-stability, oscillations, robustness and stochasticity, are critical for understanding gene regulation during cell fate decisions, inflammation and stem cell heterogeneity. However, the prevalence of the structures within gene networks that drive these dynamical behaviours, such as autoregulation or feedback by microRNAs, is unknown. We integrate transcription factor binding site (TFBS) and microRNA target data to generate a gene interaction network across 28 human tissues. This network was analysed for motifs capable of driving dynamical gene expression, including oscillations. Identified autoregulatory motifs involve 56% of transcription factors (TFs) studied. TFs that autoregulate have more interactions with microRNAs than non-autoregulatory genes and 89% of autoregulatory TFs were found in dual feedback motifs with a microRNA. Both autoregulatory and dual feedback motifs were enriched in the network. TFs that autoregulate were highly conserved between tissues. Dual feedback motifs with microRNAs were also conserved between tissues, but less so, and TFs regulate different combinations of microRNAs in a tissue-dependent manner. The study of these motifs highlights ever more genes that have complex regulatory dynamics. These data provide a resource for the identification of TFs which regulate the dynamical properties of human gene expression.

Silencing miR-370-3p rescues funny current and sinus node function in heart failure.

Bradyarrhythmias are an important cause of mortality in heart failure and previous studies indicate a mechanistic role for electrical remodelling of the key pacemaking ion channel HCN4 in this process. Here we show that, in a mouse model of heart failure in which there is sinus bradycardia, there is upregulation of a microRNA (miR-370-3p), downregulation of the pacemaker ion channel, HCN4, and downregulation of the corresponding ionic current, If, in the sinus node. In vitro, exogenous miR-370-3p inhibits HCN4 mRNA and causes downregulation of HCN4 protein, downregulation of If, and bradycardia in the isolated sinus node. In vivo, intraperitoneal injection of an antimiR to miR-370-3p into heart failure mice silences miR-370-3p and restores HCN4 mRNA and protein and If in the sinus node and blunts the sinus bradycardia. In addition, it partially restores ventricular function and reduces mortality. This represents a novel approach to heart failure treatment.

Quo vadis microRNAs?

Since 2002, published miRNAs have been collected and named by the online repository miRBase. However, with 11 000 annual publications this has become challenging. Recently, four specialized miRNA databases were published, addressing particular needs for diverse scientific communities. This development provides major opportunities for the future of miRNA annotation and nomenclature.

The Transcription Factor-microRNA Regulatory Network during hESC-chondrogenesis.

Human embryonic stem cells (ESCs) offer a promising therapeutic approach for osteoarthritis (OA). The unlimited source of cells capable of differentiating to chondrocytes has potential for repairing damaged cartilage or to generate disease models via gene editing. However their use is limited by the efficiency of chondrogenic differentiation. An improved understanding of the transcriptional and post-transcriptional regulation of chondrogenesis will enable us to improve hESC chondrogenic differentiation protocols. Small RNA-seq and whole transcriptome sequencing was performed on distinct stages of hESC-directed chondrogenesis. This revealed significant changes in the expression of several microRNAs including upregulation of known cartilage associated microRNAs and those transcribed from the Hox complexes, and the downregulation of pluripotency associated microRNAs. Integration of miRomes and transcriptomes generated during hESC-directed chondrogenesis identified key functionally related clusters of co-expressed microRNAs and protein coding genes, associated with pluripotency, primitive streak, limb development and extracellular matrix. Analysis identified regulators of hESC-directed chondrogenesis such as miR-29c-3p with 10 of its established targets identified as co-regulated 'ECM organisation' genes and miR-22-3p which is highly co-expressed with ECM genes and may regulate these genes indirectly by targeting the chondrogenic regulators SP1 and HDAC4. We identified several upregulated transcription factors including HOXA9/A10/D13 involved in limb patterning and RELA, JUN and NFAT5, which have targets enriched with ECM associated genes. We have developed an unbiased approach for integrating transcriptome and miRome using protein-protein interactions, transcription factor regulation and miRNA target interactions and identified key regulatory networks prominent in hESC chondrogenesis.

A guide to naming human non-coding RNA genes.

Research on non-coding RNA (ncRNA) is a rapidly expanding field. Providing an official gene symbol and name to ncRNA genes brings order to otherwise potential chaos as it allows unambiguous communication about each gene. The HUGO Gene Nomenclature Committee (HGNC, www.genenames.org) is the only group with the authority to approve symbols for human genes. The HGNC works with specialist advisors for different classes of ncRNA to ensure that ncRNA nomenclature is accurate and informative, where possible. Here, we review each major class of ncRNA that is currently annotated in the human genome and describe how each class is assigned a standardised nomenclature.

Restoring the constitutional alignment with a restrictive kinematic protocol improves quantitative soft-tissue balance in total knee arthroplasty: a randomized controlled trial.

AIMS: It is unknown whether kinematic alignment (KA) objectively improves knee balance in total knee arthroplasty (TKA), despite this being the biomechanical rationale for its use. This study aimed to determine whether restoring the constitutional alignment using a restrictive KA protocol resulted in better quantitative knee balance than mechanical alignment (MA). METHODS: We conducted a randomized superiority trial comparing patients undergoing TKA assigned to KA within a restrictive safe zone or MA. Optimal knee balance was defined as an intercompartmental pressure difference (ICPD) of 15 psi or less using a pressure sensor. The primary endpoint was the mean intraoperative ICPD at 10° of flexion prior to knee balancing. Secondary outcomes included balance at 45° and 90°, requirements for balancing procedures, and presence of tibiofemoral lift-off. RESULTS: A total of 63 patients (70 knees) were randomized to KA and 62 patients (68 knees) to MA. Mean ICPD at 10° flexion in the KA group was 11.7 psi (SD 13.1) compared with 32.0 psi in the MA group (SD 28.9), with a mean difference in ICPD between KA and MA of 20.3 psi (p < 0.001). Mean ICPD in the KA group was significantly lower than in the MA group at 45° and 90°, respectively (25.2 psi MA vs 14.8 psi KA, p = 0.004; 19.1 psi MA vs 11.7 psi KA, p < 0.002, respectively). Overall, participants in the KA group were more likely to achieve optimal knee balance (80% vs 35%; p < 0.001). Bone recuts to achieve knee balance were more likely to be required in the MA group (49% vs 9%; p < 0.001). More participants in the MA group had tibiofemoral lift-off (43% vs 13%; p < 0.001). CONCLUSION: This study provides persuasive evidence that restoring the constitutional alignment with KA in TKA results in a statistically significant improvement in quantitative knee balance, and further supports this technique as a viable alternative to MA. Cite this article: Bone Joint J. 2020;102-B(1):117-124.

Enhanced genome assembly and a new official gene set for Tribolium castaneum.

BACKGROUND: The red flour beetle Tribolium castaneum has emerged as an important model organism for the study of gene function in development and physiology, for ecological and evolutionary genomics, for pest control and a plethora of other topics. RNA interference (RNAi), transgenesis and genome editing are well established and the resources for genome-wide RNAi screening have become available in this model. All these techniques depend on a high quality genome assembly and precise gene models. However, the first version of the genome assembly was generated by Sanger sequencing, and with a small set of RNA sequence data limiting annotation quality. RESULTS: Here, we present an improved genome assembly (Tcas5.2) and an enhanced genome annotation resulting in a new official gene set (OGS3) for Tribolium castaneum, which significantly increase the quality of the genomic resources. By adding large-distance jumping library DNA sequencing to join scaffolds and fill small gaps, the gaps in the genome assembly were reduced and the N50 increased to 4753kbp. The precision of the gene models was enhanced by the use of a large body of RNA-Seq reads of different life history stages and tissue types, leading to the discovery of 1452 novel gene sequences. We also added new features such as alternative splicing, well defined UTRs and microRNA target predictions. For quality control, 399 gene models were evaluated by manual inspection. The current gene set was submitted to Genbank and accepted as a RefSeq genome by NCBI. CONCLUSIONS: The new genome assembly (Tcas5.2) and the official gene set (OGS3) provide enhanced genomic resources for genetic work in Tribolium castaneum. The much improved information on transcription start sites supports transgenic and gene editing approaches. Further, novel types of information such as splice variants and microRNA target genes open additional possibilities for analysis.

miRBase: from microRNA sequences to function.

miRBase catalogs, names and distributes microRNA gene sequences. The latest release of miRBase (v22) contains microRNA sequences from 271 organisms: 38 589 hairpin precursors and 48 860 mature microRNAs. We describe improvements to the database and website to provide more information about the quality of microRNA gene annotations, and the cellular functions of their products. We have collected 1493 small RNA deep sequencing datasets and mapped a total of 5.5 billion reads to microRNA sequences. The read mapping patterns provide strong support for the validity of between 20% and 65% of microRNA annotations in different well-studied animal genomes, and evidence for the removal of >200 sequences from the database. To improve the availability of microRNA functional information, we are disseminating Gene Ontology terms annotated against miRBase sequences. We have also used a text-mining approach to search for microRNA gene names in the full-text of open access articles. Over 500 000 sentences from 18 542 papers contain microRNA names. We score these sentences for functional information and link them with 12 519 microRNA entries. The sentences themselves, and word clouds built from them, provide effective summaries of the functional information about specific microRNAs. miRBase is publicly and freely available at http://mirbase.org/.

Large-scale profiling of noncoding RNA function in yeast.

Noncoding RNAs (ncRNAs) are emerging as key regulators of cellular function. We have exploited the recently developed barcoded ncRNA gene deletion strain collections in the yeast Saccharomyces cerevisiae to investigate the numerous ncRNAs in yeast with no known function. The ncRNA deletion collection contains deletions of tRNAs, snoRNAs, snRNAs, stable unannotated transcripts (SUTs), cryptic unstable transcripts (CUTs) and other annotated ncRNAs encompassing 532 different individual ncRNA deletions. We have profiled the fitness of the diploid heterozygous ncRNA deletion strain collection in six conditions using batch and continuous liquid culture, as well as the haploid ncRNA deletion strain collections arrayed individually onto solid rich media. These analyses revealed many novel environmental-specific haplo-insufficient and haplo-proficient phenotypes providing key information on the importance of each specific ncRNA in every condition. Co-fitness analysis using fitness data from the heterozygous ncRNA deletion strain collection identified two ncRNA groups required for growth during heat stress and nutrient deprivation. The extensive fitness data for each ncRNA deletion strain has been compiled into an easy to navigate database called Yeast ncRNA Analysis (YNCA). By expanding the original ncRNA deletion strain collection we identified four novel essential ncRNAs; SUT527, SUT075, SUT367 and SUT259/691. We defined the effects of each new essential ncRNA on adjacent gene expression in the heterozygote background identifying both repression and induction of nearby genes. Additionally, we discovered a function for SUT527 in the expression, 3' end formation and localization of SEC4, an essential protein coding mRNA. Finally, using plasmid complementation we rescued the SUT075 lethal phenotype revealing that this ncRNA acts in trans. Overall, our findings provide important new insights into the function of ncRNAs.

Small RNAs: Big Impact on Plant Development.

While the role of proteins in determining cell identity has been extensively studied, the contribution of small noncoding RNA molecules such as miRNAs and siRNAs has been also recognised. miRNAs bind to complementary sites in target mRNA molecules to trigger the degradation or translational inhibition of those targets. Recent studies have revealed that miRNAs play pivotal roles in key developmental processes such as patterning of the embryo, meristem, leaf, and flower. Furthermore, these miRNAs have been recruited throughout plant evolution into pathways that create diverse plant organ forms and shapes. This review focuses on the roles of miRNAs in establishing plant cell identity during key plant development processes and creating morphological diversity during plant evolution.

Abundant expression of somatic transposon-derived piRNAs throughout Tribolium castaneum embryogenesis.

BACKGROUND: Piwi-interacting RNAs (piRNAs) are a class of short (~26-31-nucleotide) non-protein-coding RNAs expressed in the metazoan germline. The piRNA pathway in arthropods is best understood in the ovary of Drosophila melanogaster, where it acts to silence active transposable elements (TEs). Maternal loading of piRNAs in oocytes is further required for the inheritance of piRNA-mediated transposon defence. However, our understanding of the diversity, evolution and function of the piRNA complement beyond drosophilids is limited. The red flour beetle, Tribolium castaneum, is an emerging model organism separated from Drosophila by ~ 350 million years of evolution that displays a number of features ancestral to arthropods, including short germ embryogenesis. Here, we characterize the maternally deposited and zygotically expressed small RNA and mRNA complements throughout T. castaneum embryogenesis. RESULTS: We find that beetle oocytes and embryos of all stages are abundant in heterogeneous ~ 28-nucleotide RNAs. These small RNAs originate from discrete genomic loci enriched in TE sequences and display the molecular signatures of transposon-derived piRNAs. In addition to the maternally loaded primary piRNAs, Tribolium embryos produce secondary piRNAs by the cleavage of zygotically activated TE transcripts via the ping-pong mechanism. The two Tribolium piRNA pathway effector proteins, Tc-Piwi/Aub and Tc-Ago3, are also expressed throughout the soma of early embryos. CONCLUSIONS: Our results show that the piRNA pathway in Tribolium is not restricted to the germline, but also operates in the embryo and may act to antagonize zygotically activated transposons. Taken together, these data highlight a functional divergence of the piRNA pathway between insects.

The house spider genome reveals an ancient whole-genome duplication during arachnid evolution.

BACKGROUND: The duplication of genes can occur through various mechanisms and is thought to make a major contribution to the evolutionary diversification of organisms. There is increasing evidence for a large-scale duplication of genes in some chelicerate lineages including two rounds of whole genome duplication (WGD) in horseshoe crabs. To investigate this further, we sequenced and analyzed the genome of the common house spider Parasteatoda tepidariorum. RESULTS: We found pervasive duplication of both coding and non-coding genes in this spider, including two clusters of Hox genes. Analysis of synteny conservation across the P. tepidariorum genome suggests that there has been an ancient WGD in spiders. Comparison with the genomes of other chelicerates, including that of the newly sequenced bark scorpion Centruroides sculpturatus, suggests that this event occurred in the common ancestor of spiders and scorpions, and is probably independent of the WGDs in horseshoe crabs. Furthermore, characterization of the sequence and expression of the Hox paralogs in P. tepidariorum suggests that many have been subject to neo-functionalization and/or sub-functionalization since their duplication. CONCLUSIONS: Our results reveal that spiders and scorpions are likely the descendants of a polyploid ancestor that lived more than 450 MYA. Given the extensive morphological diversity and ecological adaptations found among these animals, rivaling those of vertebrates, our study of the ancient WGD event in Arachnopulmonata provides a new comparative platform to explore common and divergent evolutionary outcomes of polyploidization events across eukaryotes.

RNAcentral: a comprehensive database of non-coding RNA sequences.

RNAcentral is a database of non-coding RNA (ncRNA) sequences that aggregates data from specialised ncRNA resources and provides a single entry point for accessing ncRNA sequences of all ncRNA types from all organisms. Since its launch in 2014, RNAcentral has integrated twelve new resources, taking the total number of collaborating database to 22, and began importing new types of data, such as modified nucleotides from MODOMICS and PDB. We created new species-specific identifiers that refer to unique RNA sequences within a context of single species. The website has been subject to continuous improvements focusing on text and sequence similarity searches as well as genome browsing functionality. All RNAcentral data is provided for free and is available for browsing, bulk downloads, and programmatic access at http://rnacentral.org/.