Looking for a needle in a haystack: de novo phenotypic target identification reveals Hippo pathway-mediated miR-202 regulation of egg production.

Understanding microRNA (miRNA) functions has been hampered by major difficulties in identifying their biological target(s). Currently, the main limitation is the lack of a suitable strategy to identify biologically relevant targets among a high number of putative targets. Here we provide a proof of concept of successful de novo (i.e. without prior knowledge of its identity) miRNA phenotypic target (i.e. target whose de-repression contributes to the phenotypic outcomes) identification from RNA-seq data. Using the medaka mir-202 knock-out (KO) model in which inactivation leads to a major organism-level reproductive phenotype, including reduced egg production, we introduced novel criteria including limited fold-change in KO and low interindividual variability in gene expression to reduce the list of 2853 putative targets to a short list of 5. We selected tead3b, a member of the evolutionarily-conserved Hippo pathway, known to regulate ovarian functions, due to its remarkably strong and evolutionarily conserved binding affinity for miR-202-5p. Deleting the miR-202-5p binding site in the 3' UTR of tead3b, but not of other Hippo pathway members sav1 and vgll4b, triggered a reduced egg production phenotype. This is one of the few successful examples of de novo functional assignment of a miRNA phenotypic target in vivo in vertebrates.

The Drosophila HP1 homolog Rhino is required for transposon silencing and piRNA production by dual-strand clusters.

Piwi-interacting RNAs (piRNAs) silence transposons and maintain genome integrity during germline development. In Drosophila, transposon-rich heterochromatic clusters encode piRNAs either on both genomic strands (dual-strand clusters) or predominantly one genomic strand (uni-strand clusters). Primary piRNAs derived from these clusters are proposed to drive a ping-pong amplification cycle catalyzed by proteins that localize to the perinuclear nuage. We show that the HP1 homolog Rhino is required for nuage organization, transposon silencing, and ping-pong amplification of piRNAs. rhi mutations virtually eliminate piRNAs from the dual-strand clusters and block production of putative precursor RNAs from both strands of the major 42AB dual-strand cluster, but not of transcripts or piRNAs from the uni-strand clusters. Furthermore, Rhino protein associates with the 42AB dual-strand cluster,but does not bind to uni-strand cluster 2 or flamenco. Rhino thus appears to promote transcription of dual-strand clusters, leading to production of piRNAs that drive the ping-pong amplification cycle.

Collapse of germline piRNAs in the absence of Argonaute3 reveals somatic piRNAs in flies.

Piwi-interacting RNAs (piRNAs) silence transposons in animal germ cells. piRNAs are thought to derive from long transcripts spanning transposon-rich genomic loci and to direct an autoamplification loop in which an antisense piRNA, bound to Aubergine or Piwi protein, triggers production of a sense piRNA bound to the PIWI protein Argonaute3 (Ago3). In turn, the new piRNA is envisioned to produce a second antisense piRNA. Here, we describe strong loss-of-function mutations in ago3, allowing a direct genetic test of this model. We find that Ago3 acts to amplify piRNA pools and to enforce on them an antisense bias, increasing the number of piRNAs that can act to silence transposons. We also detect a second, Ago3-independent piRNA pathway centered on Piwi. Transposons targeted by this second pathway often reside in the flamenco locus, which is expressed in somatic ovarian follicle cells, suggesting a role for piRNAs beyond the germline.

Amphioxus functional genomics and the origins of vertebrate gene regulation.

Vertebrates have greatly elaborated the basic chordate body plan and evolved highly distinctive genomes that have been sculpted by two whole-genome duplications. Here we sequence the genome of the Mediterranean amphioxus (Branchiostoma lanceolatum) and characterize DNA methylation, chromatin accessibility, histone modifications and transcriptomes across multiple developmental stages and adult tissues to investigate the evolution of the regulation of the chordate genome. Comparisons with vertebrates identify an intermediate stage in the evolution of differentially methylated enhancers, and a high conservation of gene expression and its cis-regulatory logic between amphioxus and vertebrates that occurs maximally at an earlier mid-embryonic phylotypic period. We analyse regulatory evolution after whole-genome duplications, and find that-in vertebrates-over 80% of broadly expressed gene families with multiple paralogues derived from whole-genome duplications have members that restricted their ancestral expression, and underwent specialization rather than subfunctionalization. Counter-intuitively, paralogues that restricted their expression increased the complexity of their regulatory landscapes. These data pave the way for a better understanding of the regulatory principles that underlie key vertebrate innovations.

A rationalized definition of general tumor suppressor microRNAs excludes miR-34a.

While several microRNAs (miRNAs) have been proposed to act as tumor suppressors, a consensual definition of tumor suppressing miRNAs is still missing. Similarly to coding genes, we propose that tumor suppressor miRNAs must show evidence of genetic or epigenetic inactivation in cancers, and exhibit an anti-tumorigenic (e.g., anti-proliferative) activity under endogenous expression levels. Here we observe that this definition excludes the most extensively studied tumor suppressor candidate miRNA, miR-34a. In analyzable cancer types, miR-34a does not appear to be down-regulated in primary tumors relatively to normal adjacent tissues. Deletion of miR-34a is occasionally found in human cancers, but it does not seem to be driven by an anti-tumorigenic activity of the miRNA, since it is not observed upon smaller, miR-34a-specific alterations. Its anti-proliferative action was observed upon large, supra-physiological transfection of synthetic miR-34a in cultured cells, and our data indicates that endogenous miR-34a levels do not have such an effect. Our results therefore argue against a general tumor suppressive function for miR-34a, providing an explanation to the lack of efficiency of synthetic miR-34a administration against solid tumors.

Synthetic miR-34a against solid tumours: a predictable failure.

A prolific scientific literature attributes pro- or anti-oncogenic properties to many human microRNAs ("miRNAs"). While many of these studies are based on unpersuasive analyses, one candidate suppressor tumour miRNA, miR-34a, appeared convincing enough to be administered to human patients in a clinical trial-with disappointing outcomes. Here, we review possible reasons for that failure, and their implications for other miRNAs.

Mechanistic analysis of the enhanced RNAi activity by 6-mCEPh-purine at the 5' end of the siRNA guide strand.

A key approach for improving siRNA efficacy is chemical modifications. Through an in silico screening of modifications at the 5'-end nucleobase of the guide strand, an adenine-derived compound called 6-(3-(2-carboxyethyl)phenyl)-purine (6-mCEPh-purine) was identified to improve the RNAi activity in cultured human cells and in vivo mouse models. Nevertheless, it remains unclear how this chemical modification enhances the siRNA potency. Here, we used a series of biochemical approaches to quantitatively evaluate the effect of the 6-mCEPh-purine modification at each step in the assembly of the RNAi effector complex called RISC. We found that the modification improves the formation of mature RISC at least in two different ways, by fixing the loading orientation of siRNA duplexes and increasing the stability of mature RISC after passenger strand ejection. Our data will provide a molecular platform for further development of chemically modified siRNA drugs.

Biochemistry-informed design selects potent siRNAs against SARS-CoV-2.

RNA interference (RNAi) offers an efficient way to repress genes of interest, and it is widely used in research settings. Clinical applications emerged more recently, with 5 approved siRNAs (the RNA guides of the RNAi effector complex) against human diseases. The development of siRNAs against the SARS-CoV-2 virus could therefore provide the basis of novel COVID-19 treatments, while being easily adaptable to future variants or to other, unrelated viruses. Because the biochemistry of RNAi is very precisely described, it is now possible to design siRNAs with high predicted activity and specificity using only computational tools. While previous siRNA design algorithms tended to rely on simplistic strategies (raising fully complementary siRNAs against targets of interest), our approach uses the most up-to-date mechanistic description of RNAi to allow mismatches at tolerable positions and to force them at beneficial positions, while optimizing siRNA duplex asymmetry. Our pipeline proposes 8 siRNAs against SARS-CoV-2, and ex vivo assessment confirms the high antiviral activity of 6 out of 8 siRNAs, also achieving excellent variant coverage (with several 3-siRNA combinations recognizing each correctly-sequenced variant as of September2022). Our approach is easily generalizable to other viruses as long as avariant genome database is available. With siRNA delivery procedures being currently improved, RNAi could therefore become an efficient and versatile antiviral therapeutic strategy.

Functional lability of RNA-dependent RNA polymerases in animals.

RNA interference (RNAi) requires RNA-dependent RNA polymerases (RdRPs) in many eukaryotes, and RNAi amplification constitutes the only known function for eukaryotic RdRPs. Yet in animals, classical model organisms can elicit RNAi without possessing RdRPs, and only nematode RNAi was shown to require RdRPs. Here we show that RdRP genes are much more common in animals than previously thought, even in insects, where they had been assumed not to exist. RdRP genes were present in the ancestors of numerous clades, and they were subsequently lost at a high frequency. In order to probe the function of RdRPs in a deuterostome (the cephalochordate Branchiostoma lanceolatum), we performed high-throughput analyses of small RNAs from various Branchiostoma developmental stages. Our results show that Branchiostoma RdRPs do not appear to participate in RNAi: we did not detect any candidate small RNA population exhibiting classical siRNA length or sequence features. Our results show that RdRPs have been independently lost in dozens of animal clades, and even in a clade where they have been conserved (cephalochordates) their function in RNAi amplification is not preserved. Such a dramatic functional variability reveals an unexpected plasticity in RNA silencing pathways.

microRNA target prediction programs predict many false positives.

According to the current view, each microRNA regulates hundreds of genes. Computational tools aim at identifying microRNA targets, usually selecting evolutionarily conserved microRNA binding sites. While the false positive rates have been evaluated for some prediction programs, that information is rarely put forward in studies making use of their predictions. Here, we provide evidence that such predictions are often biologically irrelevant. Focusing on miR-223-guided repression, we observed that it is often smaller than inter-individual variability in gene expression among wild-type mice, suggesting that most predicted targets are functionally insensitive to that microRNA. Furthermore, we found that human haplo-insufficient genes tend to bear the most highly conserved microRNA binding sites. It thus appears that biological functionality of microRNA binding sites depends on the dose-sensitivity of their host gene and that, conversely, it is unlikely that every predicted microRNA target is dose-sensitive enough to be functionally regulated by microRNAs. We also observed that some mRNAs can efficiently titrate microRNAs, providing a reason for microRNA binding site conservation for inefficiently repressed targets. Finally, many conserved microRNA binding sites are conserved in a microRNA-independent fashion: Sequence elements may be conserved for other reasons, while being fortuitously complementary to microRNAs. Collectively, our data suggest that the role of microRNAs in normal and pathological conditions has been overestimated due to the frequent overlooking of false positive rates.

Prospects and challenges of multi-omics data integration in toxicology.

Exposure of cells or organisms to chemicals can trigger a series of effects at the regulatory pathway level, which involve changes of levels, interactions, and feedback loops of biomolecules of different types. A single-omics technique, e.g., transcriptomics, will detect biomolecules of one type and thus can only capture changes in a small subset of the biological cascade. Therefore, although applying single-omics analyses can lead to the identification of biomarkers for certain exposures, they cannot provide a systemic understanding of toxicity pathways or adverse outcome pathways. Integration of multiple omics data sets promises a substantial improvement in detecting this pathway response to a toxicant, by an increase of information as such and especially by a systemic understanding. Here, we report the findings of a thorough evaluation of the prospects and challenges of multi-omics data integration in toxicological research. We review the availability of such data, discuss options for experimental design, evaluate methods for integration and analysis of multi-omics data, discuss best practices, and identify knowledge gaps. Re-analyzing published data, we demonstrate that multi-omics data integration can considerably improve the confidence in detecting a pathway response. Finally, we argue that more data need to be generated from studies with a multi-omics-focused design, to define which omics layers contribute most to the identification of a pathway response to a toxicant.

Argonaute proteins regulate HIV-1 multiply spliced RNA and viral production in a Dicer independent manner.

Argonaute (Ago) proteins associate with microRNAs (miRNAs) to form the core of the RNA-induced silencing complex (RISC) that mediates post-transcriptional gene silencing of target mRNAs. As key players in anti-viral defense, Ago proteins are thought to have the ability to interact with human immunodeficiency virus type 1 (HIV-1) RNA. However, the role of this interaction in regulating HIV-1 replication has been debated. Here, we used high throughput sequencing of RNA isolated by cross-linking immunoprecipitation (HITS-CLIP) to explore the interaction between Ago2 and HIV-1 RNA in infected cells. By only considering reads of 50 nucleotides length in our analysis, we identified more than 30 distinct binding sites for Ago2 along the viral RNA genome. Using reporter assays, we found four binding sites, located near splice donor sites, capable of repressing Luciferase gene expression in an Ago-dependent manner. Furthermore, inhibition of Ago1 and Ago2 levels in cells expressing HIV-1 led to an increase of viral multiply spliced transcripts and to a strong reduction in the extracellular CAp24 level. Depletion of Dicer did not affect these activities. Our results highlight a new role of Ago proteins in the control of multiply spliced HIV-1 transcript levels and viral production, independently of the miRNA pathway.

Silencing of X-Linked MicroRNAs by Meiotic Sex Chromosome Inactivation.

During the pachytene stage of meiosis in male mammals, the X and Y chromosomes are transcriptionally silenced by Meiotic Sex Chromosome Inactivation (MSCI). MSCI is conserved in therian mammals and is essential for normal male fertility. Transcriptomics approaches have demonstrated that in mice, most or all protein-coding genes on the X chromosome are subject to MSCI. However, it is unclear whether X-linked non-coding RNAs behave in a similar manner. The X chromosome is enriched in microRNA (miRNA) genes, with many exhibiting testis-biased expression. Importantly, high expression levels of X-linked miRNAs (X-miRNAs) have been reported in pachytene spermatocytes, indicating that these genes may escape MSCI, and perhaps play a role in the XY-silencing process. Here we use RNA FISH to examine X-miRNA expression in the male germ line. We find that, like protein-coding X-genes, X-miRNAs are expressed prior to prophase I and are thereafter silenced during pachynema. X-miRNA silencing does not occur in mouse models with defective MSCI. Furthermore, X-miRNAs are expressed at pachynema when present as autosomally integrated transgenes. Thus, we conclude that silencing of X-miRNAs during pachynema in wild type males is MSCI-dependent. Importantly, misexpression of X-miRNAs during pachynema causes spermatogenic defects. We propose that MSCI represents a chromosomal mechanism by which X-miRNAs, and other potential X-encoded repressors, can be silenced, thereby regulating genes with critical late spermatogenic functions.

Cnidarian microRNAs frequently regulate targets by cleavage.

In bilaterians, which comprise most of extant animals, microRNAs (miRNAs) regulate the majority of messenger RNAs (mRNAs) via base-pairing of a short sequence (the miRNA "seed") to the target, subsequently promoting translational inhibition and transcript instability. In plants, many miRNAs guide endonucleolytic cleavage of highly complementary targets. Because little is known about miRNA function in nonbilaterian animals, we investigated the repertoire and biological activity of miRNAs in the sea anemone Nematostella vectensis, a representative of Cnidaria, the sister phylum of Bilateria. Our work uncovers scores of novel miRNAs in Nematostella, increasing the total miRNA gene count to 87. Yet only a handful are conserved in corals and hydras, suggesting that microRNA gene turnover in Cnidaria greatly exceeds that of other metazoan groups. We further show that Nematostella miRNAs frequently direct the cleavage of their mRNA targets via nearly perfect complementarity. This mode of action resembles that of small interfering RNAs (siRNAs) and plant miRNAs. It appears to be common in Cnidaria, as several of the miRNA target sites are conserved among distantly related anemone species, and we also detected miRNA-directed cleavage in Hydra. Unlike in bilaterians, Nematostella miRNAs are commonly coexpressed with their target transcripts. In light of these findings, we propose that post-transcriptional regulation by miRNAs functions differently in Cnidaria and Bilateria. The similar, siRNA-like mode of action of miRNAs in Cnidaria and plants suggests that this may be an ancestral state.

Issues in current microRNA target identification methods.

microRNAs are currently believed to control a large diversity of physiologic processes, through the collective repression of thousands of target genes. Both experimental and computational analyses indeed suggest that each microRNA regulates tens or hundreds of genes. But some observations suggest that the phenotypic consequences of many published miRNA/mRNA interactions are dubious. For example, the reported amplitude of miRNA-guided repression is very small, while biologic processes tend to be robust to small changes in gene expression. We recently showed, on one particular miRNA, that for most predicted targets, miRNA-guided repression is even smaller than inter-individual variability among wild-type specimens. We also put forward several sources of computational false positives. These issues are generally neglected by the scientific community, probably resulting in the frequent publication of irreproducible or misinterpreted results regarding microRNA function. We propose novel types of analyses, easily accessible to the community, that could help improve microRNA target identification.

microRNAs and the evolution of complex multicellularity: identification of a large, diverse complement of microRNAs in the brown alga Ectocarpus.

There is currently convincing evidence that microRNAs have evolved independently in at least six different eukaryotic lineages: animals, land plants, chlorophyte green algae, demosponges, slime molds and brown algae. MicroRNAs from different lineages are not homologous but some structural features are strongly conserved across the eukaryotic tree allowing the application of stringent criteria to identify novel microRNA loci. A large set of 63 microRNA families was identified in the brown alga Ectocarpus based on mapping of RNA-seq data and nine microRNAs were confirmed by northern blotting. The Ectocarpus microRNAs are highly diverse at the sequence level with few multi-gene families, and do not tend to occur in clusters but exhibit some highly conserved structural features such as the presence of a uracil at the first residue. No homologues of Ectocarpus microRNAs were found in other stramenopile genomes indicating that they emerged late in stramenopile evolution and are perhaps specific to the brown algae. The large number of microRNA loci in Ectocarpus is consistent with the developmental complexity of many brown algal species and supports a proposed link between the emergence and expansion of microRNA regulatory systems and the evolution of complex multicellularity.

Applying 'omics technologies in chemicals risk assessment: Report of an ECETOC workshop.

Prevailing knowledge gaps in linking specific molecular changes to apical outcomes and methodological uncertainties in the generation, storage, processing, and interpretation of 'omics data limit the application of 'omics technologies in regulatory toxicology. Against this background, the European Centre for Ecotoxicology and Toxicology of Chemicals (ECETOC) convened a workshop Applying 'omics technologies in chemicals risk assessment that is reported herein. Ahead of the workshop, multi-expert teams drafted frameworks on best practices for (i) a Good-Laboratory Practice-like context for collecting, storing and curating 'omics data; (ii) the processing of 'omics data; and (iii) weight-of-evidence approaches for integrating 'omics data. The workshop participants confirmed the relevance of these Frameworks to facilitate the regulatory applicability and use of 'omics data, and the workshop discussions provided input for their further elaboration. Additionally, the key objective (iv) to establish approaches to connect 'omics perturbations to phenotypic alterations was addressed. Generally, it was considered promising to strive to link gene expression changes and pathway perturbations to the phenotype by mapping them to specific adverse outcome pathways. While further work is necessary before gene expression changes can be used to establish safe levels of substance exposure, the ECETOC workshop provided important incentives towards achieving this goal.

Coding and non-coding variants in the SHOX2 gene in patients with early-onset atrial fibrillation.

Atrial fibrillation (AF) is the most prevalent cardiac arrhythmia with a strong genetic component. Molecular pathways involving the homeodomain transcription factor Shox2 control the development and function of the cardiac conduction system in mouse and zebrafish. Here we report the analysis of human SHOX2 as a potential susceptibility gene for early-onset AF. To identify causal variants and define the underlying mechanisms, results from 378 patients with early-onset AF before the age of 60 years were analyzed and compared to 1870 controls or reference datasets. We identified two missense mutations (p.G81E, p.H283Q), that were predicted as damaging. Transactivation studies using SHOX2 targets and phenotypic rescue experiments in zebrafish demonstrated that the p.H283Q mutation severely affects SHOX2 pacemaker function. We also demonstrate an association between a 3'UTR variant c.*28T>C of SHOX2 and AF (p = 0.00515). Patients carrying this variant present significantly longer PR intervals. Mechanistically, this variant creates a functional binding site for hsa-miR-92b-5p. Circulating hsa-miR-92b-5p plasma levels were significantly altered in AF patients carrying the 3'UTR variant (p = 0.0095). Finally, we demonstrate significantly reduced SHOX2 expression levels in right atrial appendages of AF patients compared to patients with sinus rhythm. Together, these results suggest a genetic contribution of SHOX2 in early-onset AF.

Advancing the use of noncoding RNA in regulatory toxicology: Report of an ECETOC workshop.

The European Centre for the Ecotoxicology and Toxicology of Chemicals (ECETOC) organised a workshop to discuss the state-of-the-art research on noncoding RNAs (ncRNAs) as biomarkers in regulatory toxicology and as analytical and therapeutic agents. There was agreement that ncRNA expression profiling data requires careful evaluation to determine the utility of specific ncRNAs as biomarkers. To advance the use of ncRNA in regulatory toxicology, the following research priorities were identified: (1) Conduct comprehensive literature reviews to identify possibly suitable ncRNAs and areas of toxicology where ncRNA expression profiling could address prevailing scientific deficiencies. (2) Develop consensus on how to conduct ncRNA expression profiling in a toxicological context. (3) Conduct experimental projects, including, e.g., rat (90-day) oral toxicity studies, to evaluate the toxicological relevance of the expression profiles of selected ncRNAs. Thereby, physiological ncRNA expression profiles should be established, including the biological variability of healthy individuals. To substantiate the relevance of key ncRNAs for cell homeostasis or pathogenesis, molecular events should be dose-dependently linked with substance-induced apical effects. Applying a holistic approach, knowledge on ncRNAs, 'omics and epigenetics technologies should be integrated into adverse outcome pathways to improve the understanding of the functional roles of ncRNAs within a regulatory context.