RNA-RNA interactions between respiratory syncytial virus and miR-26 and miR-27 are associated with regulation of cell cycle and antiviral immunity.

microRNAs (miRNAs) regulate nearly all physiological processes but our understanding of exactly how they function remains incomplete, particularly in the context of viral infections. Here, we adapt a biochemical method (CLEAR-CLIP) and analysis pipeline to identify targets of miRNAs in lung cells infected with Respiratory syncytial virus (RSV). We show that RSV binds directly to miR-26 and miR-27 through seed pairing and demonstrate that these miRNAs target distinct gene networks associated with cell cycle and metabolism (miR-27) and antiviral immunity (miR-26). Many of the targets are de-repressed upon infection and we show that the miR-27 targets most sensitive to miRNA inhibition are those associated with cell cycle. Finally, we demonstrate that high confidence chimeras map to long noncoding RNAs (lncRNAs) and pseudogenes in transcriptional regulatory regions. We validate that a proportion of miR-27 and Argonaute 2 (AGO2) is nuclear and identify a long non-coding RNA (lncRNA) as a miR-27 target that is linked to transcriptional regulation of nearby genes. This work expands the target networks of miR-26 and miR-27 to include direct interactions with RSV and lncRNAs and implicate these miRNAs in regulation of key genes that impact the viral life cycle associated with cell cycle, metabolism, and antiviral immunity.

Hybridization alters maternal and paternal genome contributions to early plant embryogenesis.

After fertilization, zygotic genome activation results in a transcriptionally competent embryo. Hybrid transcriptome experiments in Arabidopsis have concluded that the maternal and paternal genomes make equal contributions to zygotes and embryos, yet embryo defective (emb) mutants in the Columbia (Col) ecotype display early maternal effects. Here, we show that hybridization of Col with Landsberg erecta (Ler) or Cape Verde Islands (Cvi) ecotypes decreases the maternal effects of emb mutants. Reanalysis of Col/Ler and Col/Cvi transcriptomes confirmed equal parental contributions in Col/Cvi early embryos. By contrast, thousands of genes in Col/Ler zygotes and one-cell embryos were biallelic in one cross and monoallelic in the reciprocal cross, with analysis of intron reads pointing to active transcription as responsible for this parent-of-origin bias. Our analysis shows that, contrary to previous conclusions, the maternal and paternal genomes in Col/Ler zygotes are activated in an asymmetric manner. The decrease in maternal effects in hybrid embryos compared with those in isogenic Col along with differences in genome activation between Col/Cvi and Col/Ler suggest that neither of these hybrids accurately reflects the general trends of parent-of-origin regulation in Arabidopsis embryogenesis.

miR-7 is recruited to the high molecular weight RNA-induced silencing complex in CD8+ T cells upon activation and suppresses IL-2 signaling.

Increasing evidence suggests mammalian Argonaute (Ago) proteins partition into distinct complexes within cells, but there is still little biochemical or functional understanding of the miRNAs differentially associated with these complexes. In naïve T cells, Ago2 is found almost exclusively in low molecular weight (LMW) complexes which are associated with miRNAs but not their target mRNAs. Upon T-cell activation, a proportion of these Ago2 complexes move into a newly formed high molecular weight (HMW) RNA-induced silencing complex (RISC), which is characterized by the presence of the GW182 protein that mediates translational repression. Here, we demonstrate distinct partitioning of miRNAs and isomiRs in LMW versus HMW RISCs upon antigen-mediated activation of CD8+ T cells. We identify miR-7 as highly enriched in HMW RISC and demonstrate that miR-7 inhibition leads to increased production of IL-2 and up-regulation of the IL-2 receptor, the transferrin receptor, CD71 and the amino acid transporter, CD98. Our data support a model where recruitment of miR-7 to HMW RISC restrains IL-2 signaling and the metabolic processes regulated by IL-2.

Axonal precursor miRNAs hitchhike on endosomes and locally regulate the development of neural circuits.

Various species of non-coding RNAs (ncRNAs) are enriched in specific subcellular compartments, but the mechanisms orchestrating their localization and their local functions remain largely unknown. We investigated both aspects using the elongating retinal ganglion cell axon and its tip, the growth cone, as models. We reveal that specific endogenous precursor microRNAs (pre-miRNAs) are actively trafficked to distal axons by hitchhiking primarily on late endosomes/lysosomes. Upon exposure to the axon guidance cue semaphorin 3A (Sema3A), pre-miRNAs are processed specifically within axons into newly generated miRNAs, one of which, in turn, silences the basal translation of tubulin beta 3 class III (TUBB3), but not amyloid beta precursor protein (APP). At the organismal level, these mature miRNAs are required for growth cone steering and a fully functional visual system. Overall, our results uncover a novel mode of ncRNA transport from one cytosolic compartment to another within polarized cells. They also reveal that newly generated miRNAs are critical components of a ncRNA-based signaling pathway that transduces environmental signals into the structural remodeling of subcellular compartments.

Gene expression differentiation in the reproductive tissues of Drosophila willistoni subspecies and their hybrids.

Early lineage diversification is central to understand what mutational events drive species divergence. Particularly, gene misregulation in interspecific hybrids can inform about what genes and pathways underlie hybrid dysfunction. In Drosophila hybrids, how regulatory evolution impacts different reproductive tissues remains understudied. Here, we generate a new genome assembly and annotation in Drosophila willistoni and analyse the patterns of transcriptome divergence between two allopatrically evolved D. willistoni subspecies, their male sterile and female fertile hybrid progeny across testis, male accessory gland, and ovary. Patterns of transcriptome divergence and modes of regulatory evolution were tissue-specific. Despite no indication for cell-type differences in hybrid testis, this tissue exhibited the largest magnitude of expression differentiation between subspecies and between parentals and hybrids. No evidence for anomalous dosage compensation in hybrid male tissues was detected nor was a differential role for the neo- and the ancestral arms of the D. willistoni X chromosome. Compared to the autosomes, the X chromosome appeared enriched for transgressively expressed genes in testis despite being the least differentiated in expression between subspecies. Evidence for fine genome clustering of transgressively expressed genes suggests a role of chromatin structure on hybrid gene misregulation. Lastly, transgressively expressed genes in the testis of the sterile male progeny were enriched for GO terms not typically associated with sperm function, instead hinting at anomalous development of the reproductive tissue. Our thorough tissue-level portrait of transcriptome differentiation between recently diverged D. willistoni subspecies and their hybrids provides a more nuanced view of early regulatory changes during speciation.

Translation Comes First: Ancient and Convergent Selection of Codon Usage Bias Across Prokaryotic Genomes.

Codon usage is the outcome of different evolutionary processes and can inform us about the conditions in which organisms live and evolve. Here, we present R_ENC', which is an improvement to the original S index developed by dos Reis et al. (2004). Our index is less sensitive to G+C content, which greatly affects synonymous codon usage in prokaryotes, making it better suited to detect selection acting on codon usage. We used R_ENC' to estimate the extent of selected codon usage bias in 1800 genomes representing 26 prokaryotic phyla. We found that Gammaproteobacteria, Betaproteobacteria, Actinobacteria, and Firmicutes are the phyla/subphyla showing more genomes with selected codon usage bias. In particular, we found that several lineages within Gammaproteobacteria and Firmicutes show a similar set of functional terms enriched in genes under selected codon usage bias, indicating convergent evolution. We also show that selected codon usage bias tends to evolve in genes coding for the translation machinery before other functional GO terms. Finally, we discuss the possibility to use R_ENC' to predict whether lineages evolved in copiotrophic or oligotrophic environments.

Multiscale analysis of the randomization limits of the chromosomal gene organization between Lepidoptera and Diptera.

How chromosome gene organization and gene content evolve among distantly related and structurally malleable genomes remains unresolved. This is particularly the case when considering different insect orders. We have compared the highly contiguous genome assemblies of the lepidopteran Danaus plexippus and the dipteran Drosophila melanogaster, which shared a common ancestor around 290 Ma. The gene content of 23 out of 30 D. plexippus chromosomes was significantly associated with one or two of the six chromosomal elements of the Drosophila genome, denoting common ancestry. Despite the phylogenetic distance, 9.6% of the 1-to-1 orthologues still reside within the same ancestral genome neighbourhood. Furthermore, the comparison D. plexippus-Bombyx mori indicated that the rates of chromosome repatterning are lower in Lepidoptera than in Diptera, although still within the same order of magnitude. Concordantly, 14 developmental gene clusters showed a higher tendency to retain full or partial clustering in D. plexippus, further supporting that the physical association between the SuperHox and NK clusters existed in the ancestral bilaterian. Our results illuminate the scope and limits of the evolution of the gene organization and content of the ancestral chromosomes to the Lepidoptera and Diptera while helping reconstruct portions of the genome in their most recent common ancestor.

Global analyses of the effect of different cellular contexts on microRNA targeting.

MicroRNA (miRNA) regulation clearly impacts animal development, but the extent to which development-with its resulting diversity of cellular contexts-impacts miRNA regulation is unclear. Here, we compared cohorts of genes repressed by the same miRNAs in different cell lines and tissues and found that target repertoires were largely unaffected, with secondary effects explaining most of the differential responses detected. Outliers resulting from differential direct targeting were often attributable to alternative 3' UTR isoform usage that modulated the presence of miRNA sites. More inclusive examination of alternative 3' UTR isoforms revealed that they influence ∼10% of predicted targets when comparing any two cell types. Indeed, considering alternative 3' UTR isoform usage improved prediction of targeting efficacy significantly beyond the improvements observed when considering constitutive isoform usage. Thus, although miRNA targeting is remarkably consistent in different cell types, considering the 3' UTR landscape helps predict targeting efficacy and explain differential regulation that is observed.

Disentangling sRNA-Seq data to study RNA communication between species.

Many organisms exchange small RNAs (sRNAs) during their interactions, that can target or bolster defense strategies in host-pathogen systems. Current sRNA-Seq technology can determine the sRNAs present in any symbiotic system, but there are very few bioinformatic tools available to interpret the results. We show that one of the biggest challenges comes from sequences that map equally well to the genomes of both interacting organisms. This arises due to the small size of the sRNAs compared to large genomes, and because a large portion of sequenced sRNAs come from genomic regions that encode highly conserved miRNAs, rRNAs or tRNAs. Here, we present strategies to disentangle sRNA-Seq data from samples of communicating organisms, developed using diverse plant and animal species that are known to receive or exchange RNA with their symbionts. We show that sequence assembly, both de novo and genome-guided, can be used for these sRNA-Seq data, greatly reducing the ambiguity of mapping reads. Even confidently mapped sequences can be misleading, so we further demonstrate the use of differential expression strategies to determine true parasite-derived sRNAs within host cells. We validate our methods on new experiments designed to probe the nature of the extracellular vesicle sRNAs from the parasitic nematode Heligmosomoides bakeri that get into mouse intestinal epithelial cells.

Secretion of an Argonaute protein by a parasitic nematode and the evolution of its siRNA guides.

Extracellular RNA has been proposed to mediate communication between cells and organisms however relatively little is understood regarding how specific sequences are selected for export. Here, we describe a specific Argonaute protein (exWAGO) that is secreted in extracellular vesicles (EVs) released by the gastrointestinal nematode Heligmosomoides bakeri, at multiple copies per EV. Phylogenetic and gene expression analyses demonstrate exWAGO orthologues are highly conserved and abundantly expressed in related parasites but highly diverged in free-living genus Caenorhabditis. We show that the most abundant small RNAs released from the nematode parasite are not microRNAs as previously thought, but rather secondary small interfering RNAs (siRNAs) that are produced by RNA-dependent RNA Polymerases. The siRNAs that are released in EVs have distinct evolutionary properties compared to those resident in free-living or parasitic nematodes. Immunoprecipitation of exWAGO demonstrates that it specifically associates with siRNAs from transposons and newly evolved repetitive elements that are packaged in EVs and released into the host environment. Together this work demonstrates molecular and evolutionary selectivity in the small RNA sequences that are released in EVs into the host environment and identifies a novel Argonaute protein as the mediator of this.

Danger signals activate a putative innate immune system during regeneration in a filamentous fungus.

The ability to respond to injury is a biological process shared by organisms of different kingdoms that can even result in complete regeneration of a part or structure that was lost. Due to their immobility, multicellular fungi are prey to various predators and are therefore constantly exposed to mechanical damage. Nevertheless, our current knowledge of how fungi respond to injury is scarce. Here we show that activation of injury responses and hyphal regeneration in the filamentous fungus Trichoderma atroviride relies on the detection of two danger or alarm signals. As an early response to injury, we detected a transient increase in cytosolic free calcium ([Ca2+]c) that was promoted by extracellular ATP, and which is likely regulated by a mechanism of calcium-induced calcium-release. In addition, we demonstrate that the mitogen activated protein kinase Tmk1 plays a key role in hyphal regeneration. Calcium- and Tmk1-mediated signaling cascades activated major transcriptional changes early following injury, including induction of a set of regeneration associated genes related to cell signaling, stress responses, transcription regulation, ribosome biogenesis/translation, replication and DNA repair. Interestingly, we uncovered the activation of a putative fungal innate immune response, including the involvement of HET domain genes, known to participate in programmed cell death. Our work shows that fungi and animals share danger-signals, signaling cascades, and the activation of the expression of genes related to immunity after injury, which are likely the result of convergent evolution.

Arabidopsis thaliana miRNAs promote embryo pattern formation beginning in the zygote.

miRNAs are essential regulators of cell identity, yet their role in early embryo development in plants remains largely unexplored. To determine the earliest stage at which miRNAs act to promote pattern formation in embryogenesis, we examined a series of mutant alleles in the Arabidopsis thaliana miRNA biogenesis enzymes DICER-LIKE 1 (DCL1), SERRATE (SE), and HYPONASTIC LEAVES 1 (HYL1). Cellular and patterning defects were observed in dcl1, se and hyl1 embryos from the zygote through the globular stage of embryogenesis. To identify miRNAs that are expressed in early embryogenesis, we sequenced mRNAs from globular stage Columbia wild type (wt) and se-1 embryos, and identified transcripts potentially corresponding to 100 miRNA precursors. Considering genome location and transcript increase between wt and se-1, 39 of these MIRNAs are predicted to be bona fide early embryo miRNAs. Among these are conserved miRNAs such as miR156, miR159, miR160, miR161, miR164, miR165, miR166, miR167, miR168, miR171, miR319, miR390 and miR394, as well as miRNAs whose function has never been characterized. Our analysis demonstrates that miRNAs promote pattern formation beginning in the zygote, and provides a comprehensive dataset for functional studies of individual miRNAs in Arabidopsis embryogenesis.

Transcriptome diversity among rice root types during asymbiosis and interaction with arbuscular mycorrhizal fungi.

Root systems consist of different root types (RTs) with distinct developmental and functional characteristics. RTs may be individually reprogrammed in response to their microenvironment to maximize adaptive plasticity. Molecular understanding of such specific remodeling--although crucial for crop improvement--is limited. Here, RT-specific transcriptomes of adult rice crown, large and fine lateral roots were assessed, revealing molecular evidence for functional diversity among individual RTs. Of the three rice RTs, crown roots displayed a significant enrichment of transcripts associated with phytohormones and secondary cell wall (SCW) metabolism, whereas lateral RTs showed a greater accumulation of transcripts related to mineral transport. In nature, arbuscular mycorrhizal (AM) symbiosis represents the default state of most root systems and is known to modify root system architecture. Rice RTs become heterogeneously colonized by AM fungi, with large laterals preferentially entering into the association. However, RT-specific transcriptional responses to AM symbiosis were quantitatively most pronounced for crown roots despite their modest physical engagement in the interaction. Furthermore, colonized crown roots adopted an expression profile more related to mycorrhizal large lateral than to noncolonized crown roots, suggesting a fundamental reprogramming of crown root character. Among these changes, a significant reduction in SCW transcripts was observed that was correlated with an alteration of SCW composition as determined by mass spectrometry. The combined change in SCW, hormone- and transport-related transcript profiles across the RTs indicates a previously overlooked switch of functional relationships among RTs during AM symbiosis, with a potential impact on root system architecture and functioning.

RNA-mediated communication between helminths and their hosts: The missing links.

Small RNAs have been discovered in a wide variety of extracellular environments and are now thought to participate in communication between cells and even between different organisms and species. Helminths are parasitic worms that generally reside in extracellular niches in their hosts and can establish chronic infection through the release of immunomodulatory factors. Recent work has demonstrated that Extracellular RNA (exRNA) may be another class of immunomodulator secreted by helminths. Here we will detail what is known about small RNA pathways in helminth pathogens (focusing on nematodes) and mammalian hosts. We will then explore the computational challenges with identifying RNA-RNA interactions between 2 different species and the paradigm of RNA-RNA co-evolution that accompanies this. Finally we explore the lingering questions that require further investigation to understand the properties of exRNA that would enable it to function as an immunomodulator.

The endonuclease activity of Mili fuels piRNA amplification that silences LINE1 elements.

Piwi proteins and Piwi-interacting RNAs (piRNAs) have conserved functions in transposon silencing. The murine Piwi proteins Mili and Miwi2 (also called Piwil2 and Piwil4, respectively) direct epigenetic LINE1 and intracisternal A particle transposon silencing during genome reprogramming in the embryonic male germ line. Piwi proteins are proposed to be piRNA-guided endonucleases that initiate secondary piRNA biogenesis; however, the actual contribution of their endonuclease activities to piRNA biogenesis and transposon silencing remain unknown. To investigate the role of Piwi-catalysed endonucleolytic activity, we engineered point mutations in mice that substitute the second aspartic acid to an alanine in the DDH catalytic triad of Mili and Miwi2, generating the Mili(DAH) and Miwi2(DAH) alleles, respectively. Analysis of Mili-bound piRNAs from homozygous Mili(DAH) fetal gonadocytes revealed a failure of transposon piRNA amplification, resulting in the marked reduction of piRNA bound within Miwi2 ribonuclear particles. We find that Mili-mediated piRNA amplification is selectively required for LINE1, but not intracisternal A particle, silencing. The defective piRNA pathway in Mili(DAH) mice results in spermatogenic failure and sterility. Surprisingly, homozygous Miwi2(DAH) mice are fertile, transposon silencing is established normally and no defects in secondary piRNA biogenesis are observed. In addition, the hallmarks of piRNA amplification are observed in Miwi2-deficient gonadocytes. We conclude that cycles of intra-Mili secondary piRNA biogenesis fuel piRNA amplification that is absolutely required for LINE1 silencing.

Zebrafish miR-1 and miR-133 shape muscle gene expression and regulate sarcomeric actin organization.

microRNAs (miRNAs) represent approximately 4% of the genes in vertebrates, where they regulate deadenylation, translation, and decay of the target messenger RNAs (mRNAs). The integrated role of miRNAs to regulate gene expression and cell function remains largely unknown. Therefore, to identify the targets coordinately regulated by muscle miRNAs in vivo, we performed gene expression arrays on muscle cells sorted from wild type, dicer mutants, and single miRNA knockdown embryos. Our analysis reveals that two particular miRNAs, miR-1 and miR-133, influence gene expression patterns in the zebrafish embryo where they account for >54% of the miRNA-mediated regulation in the muscle. We also found that muscle miRNA targets (1) tend to be expressed at low levels in wild-type muscle but are more highly expressed in dicer mutant muscle, and (2) are enriched for actin-related and actin-binding proteins. Loss of dicer function or down-regulation of miR-1 and miR-133 alters muscle gene expression and disrupts actin organization during sarcomere assembly. These results suggest that miR-1 and miR-133 actively shape gene expression patterns in muscle tissue, where they regulate sarcomeric actin organization.

Improving microRNA target prediction with gene expression profiles.

BACKGROUND: Mammalian genomes encode for thousands of microRNAs, which can potentially regulate the majority of protein-coding genes. They have been implicated in development and disease, leading to great interest in understanding their function, with computational methods being widely used to predict their targets. Most computational methods rely on sequence features, thermodynamics, and conservation filters; essentially scanning the whole transcriptome to predict one set of targets for each microRNA. This has the limitation of not considering that the same microRNA could have different sets of targets, and thus different functions, when expressed in different types of cells. RESULTS: To address this problem, we combine popular target prediction methods with expression profiles, via machine learning, to produce a new predictor: TargetExpress. Using independent data from microarrays and high-throughput sequencing, we show that TargetExpress outperforms existing methods, and that our predictions are enriched in functions that are coherent with the added expression profile and literature reports. CONCLUSIONS: Our method should be particularly useful for anyone studying the functions and targets of miRNAs in specific tissues or cells. TargetExpress is available at: http://targetexpress.ceiabreulab.org/ .

Extent, causes, and consequences of small RNA expression variation in human adipose tissue.

Small RNAs are functional molecules that modulate mRNA transcripts and have been implicated in the aetiology of several common diseases. However, little is known about the extent of their variability within the human population. Here, we characterise the extent, causes, and effects of naturally occurring variation in expression and sequence of small RNAs from adipose tissue in relation to genotype, gene expression, and metabolic traits in the MuTHER reference cohort. We profiled the expression of 15 to 30 base pair RNA molecules in subcutaneous adipose tissue from 131 individuals using high-throughput sequencing, and quantified levels of 591 microRNAs and small nucleolar RNAs. We identified three genetic variants and three RNA editing events. Highly expressed small RNAs are more conserved within mammals than average, as are those with highly variable expression. We identified 14 genetic loci significantly associated with nearby small RNA expression levels, seven of which also regulate an mRNA transcript level in the same region. In addition, these loci are enriched for variants significant in genome-wide association studies for body mass index. Contrary to expectation, we found no evidence for negative correlation between expression level of a microRNA and its target mRNAs. Trunk fat mass, body mass index, and fasting insulin were associated with more than twenty small RNA expression levels each, while fasting glucose had no significant associations. This study highlights the similar genetic complexity and shared genetic control of small RNA and mRNA transcripts, and gives a quantitative picture of small RNA expression variation in the human population.

The RNAi machinery regulates growth and development in the filamentous fungus Trichoderma atroviride.

The RNAi machinery is generally involved in genome protection in filamentous fungi; however, the physiological role of RNAi has been poorly studied in fungal models. Here, we report that in the filamentous fungus Trichoderma atroviride, the products of the dcr2 and rdr3 genes control reproductive development, because mutations in these genes affect conidiation. In addition, Dcr1 together with Dcr2 control vegetative growth since Δdcr1, Δdcr2 and Δdcr1Δdcr2 present morphological alterations. Whole-genome transcriptional analysis of WT, Δdcr1, Δdcr2 and Δdcr1Δdcr2 show that each Dicer controls different biological processes, such as development or metabolism, which could explain the lack of conidiation in the mutants. Finally, we observed sRNAs that are differentially expressed in the WT and Δdcr2. The expression of some of these sRNAs correlates with the expression of differential transcripts, suggesting that these mRNAs may contain the corresponding targets. Together these data show that in T. atroviride, the RNAi machinery plays a central role in endogenous processes such as development and fitness, beyond controlling genome protection against invasive nucleic acids as reported for other fungi.

Post-developmental microRNA expression is required for normal physiology, and regulates aging in parallel to insulin/IGF-1 signaling in C. elegans.

Regulation of gene expression by microRNAs (miRNAs) is essential for normal development, but the roles of miRNAs in the physiology of adult animals are poorly understood. We have isolated a conditional allele of DGCR8/pash-1, which allows reversible and rapid inactivation of miRNA synthesis in vivo in Caenorhabditis elegans. This is a powerful new tool that allows dissection of post-developmental miRNA functions. We demonstrate that continuous synthesis of miRNAs is dispensable for cellular viability but critical for the physiology of adult animals. Loss of miRNA synthesis in the adult reduces lifespan and results in rapid aging. The insulin/IGF-1 signaling pathway is a critical determinant of lifespan, and is modulated by miRNAs. We find that although miRNA expression is required for some mechanisms of lifespan extension, it is not essential for the longevity of animals lacking insulin/IGF-1 signaling. Further, misregulated insulin/IGF-1 signaling cannot account for the reduced lifespan caused by disruption of miRNA synthesis. We show that miRNAs act in parallel with insulin/IGF-1 signaling to regulate a shared set of downstream genes important for physiological processes that determine lifespan. We conclude that coordinated transcriptional and post-transcriptional regulation of gene expression promotes longevity.