Identification of ascorbate- and salicylate-responsive miRNAs and verification of the spectral control of miR395 in Arabidopsis.

We assumed that miRNAs might regulate the physiological and biochemical processes in plants through their effects on the redox system and phytohormones. To check this hypothesis, the transcriptome profile of wild-type Arabidopsis and lines with decreased ascorbate (Asc), glutathione (GSH), or salicylate (Sal) levels were compared. GSH deficiency did not influence the miRNA expression, whereas lower levels of Asc and Sal reduced the accumulation of 9 and 44 miRNAs, respectively, but only four miRNAs were upregulated. Bioinformatics analysis revealed that their over-represented target genes are associated with the synthesis of nitrogen-containing and aromatic compounds, nucleic acids, and sulphate assimilation. Among them, the sulphate reduction-related miR395 - ATP-sulfurylase couple was selected to check the assumed modulating role of the light spectrum. A greater induction of the Asc- and Sal-responsive miR395 was observed under sulphur starvation in far-red light compared to white and blue light in wild-type and GSH-deficient Arabidopsis lines. Sal deficiency inhibited the induction of miR395 by sulphur starvation in blue light, whereas Asc deficiency greatly reduced it independently of the spectrum. Interestingly, sulphur starvation decreased only the level of ATP sulfurylase 4 among the miR395 target genes in far-red light. The expression level of ATP sulfurylase 3 was higher in far-red light than in blue light in wild-type and Asc-deficient lines. The results indicate the coordinated control of miRNAs by the redox and hormonal system since 11 miRNAs were affected by both Asc and Sal deficiency. This process can be modulated by light spectrum, as shown for miR395.

microRNA-seq of cartilage reveals an overabundance of miR-140-3p which contains functional isomiRs.

miR-140 is selectively expressed in cartilage. Deletion of the entire Mir140 locus in mice results in growth retardation and early-onset osteoarthritis-like pathology; however, the relative contribution of miR-140-5p or miR-140-3p to the phenotype remains to be determined. An unbiased small RNA sequencing approach identified miR-140-3p as significantly more abundant (>10-fold) than miR-140-5p in human cartilage. Analysis of these data identified multiple miR-140-3p isomiRs differing from the miRBase annotation at both the 5' and 3' end, with >99% having one of two seed sequences (5' bases 2-8). Canonical (miR-140-3p.2) and shifted (miR-140-3p.1) seed isomiRs were overexpressed in chondrocytes and transcriptomics performed to identify targets. miR-140-3p.1 and miR-140-3p.2 significantly down-regulated 694 and 238 genes, respectively, of which only 162 genes were commonly down-regulated. IsomiR targets were validated using 3'UTR luciferase assays. miR-140-3p.1 targets were enriched within up-regulated genes in rib chondrocytes of Mir140-null mice and within down-regulated genes during human chondrogenesis. Finally, through imputing the expression of miR-140 from the expression of the host gene WWP2 in 124 previously published data sets, an inverse correlation with miR-140-3p.1 predicted targets was identified. Together these data suggest the novel seed containing isomiR miR-140-3p.1 is more functional than original consensus miR-140-3p seed containing isomiR.

Mechanistic insights into non-coding Y RNA processing.

Y RNAs (84-112 nt) are non-coding RNAs transcribed by RNA polymerase III and are characterized by a distinctive secondary structure. Human Y RNAs interact with the autoimmune proteins SSB and RO60 that together form a ribonucleoprotein (RNP) complex termed RoRNP and Y RNAs also perform regulatory roles in DNA and RNA replication and stability, which has major implications for diseases including cancer. During cellular stress and apoptosis, Y RNAs are cleaved into 3' and 5' end fragments termed Y RNA-derived small RNAs (ysRNAs). Although some ysRNA functions in stress, apoptosis and cancer have been reported, their fundamental biogenesis has not been described. Here we report that 3' end RNY5 cleavage is structure dependent. In high throughput mutagenesis experiments, cleavage occurred between the 2nd and 3rd nt above a double stranded stem comprising high GC content. We demonstrate that an internal loop above stem S3 is critical for producing 3' end ysRNAs (31 nt) with mutants resulting in longer or no ysRNAs. We show a UGGGU sequence motif at position 22 of RNY5 is critical for producing 5' end ysRNAs (22-25 nt). We show that intact RO60 is critical for ysRNA biogenesis. We conclude that ribonuclease L (RNASEL) contributes to Y RNA cleavage in mouse embryonic fibroblasts but is not the only endoribonuclease important in human cells.

Molecular insights into an ancient form of Paget's disease of bone.

Paget's disease of bone (PDB) is a chronic skeletal disorder that can affect one or several bones in individuals older than 55 y of age. PDB-like changes have been reported in archaeological remains as old as Roman, although accurate diagnosis and natural history of the disease is lacking. Six skeletons from a collection of 130 excavated at Norton Priory in the North West of England, which dates to medieval times, show atypical and extensive pathological changes resembling contemporary PDB affecting as many as 75% of individual skeletons. Disease prevalence in the remaining collection is high, at least 16% of adults, with age at death estimations as low as 35 y. Despite these atypical features, paleoproteomic analysis identified sequestosome 1 (SQSTM1) or p62, a protein central to the pathological milieu of PDB, as one of the few noncollagenous human sequences preserved in skeletal samples. Targeted proteomic analysis detected >60% of the ancient p62 primary sequence, with Western blotting indicating p62 abnormalities, including in dentition. Direct sequencing of ancient DNA excluded contemporary PDB-associated SQSTM1 mutations. Our observations indicate that the ancient p62 protein is likely modified within its C-terminal ubiquitin-associated domain. Ancient miRNAs were remarkably preserved in an osteosarcoma from a skeleton with extensive disease, with miR-16 expression consistent with that reported in contemporary PDB-associated bone tumors. Our work displays the use of proteomics to inform diagnosis of ancient diseases such as atypical PDB, which has unusual features presumably potentiated by yet-unidentified environmental or genetic factors.

Gene expression during larval caste determination and differentiation in intermediately eusocial bumblebees, and a comparative analysis with advanced eusocial honeybees.

The queen-worker caste system of eusocial insects represents a prime example of developmental polyphenism (environmentally-induced phenotypic polymorphism) and is intrinsic to the evolution of advanced eusociality. However, the comparative molecular basis of larval caste determination and subsequent differentiation in the eusocial Hymenoptera remains poorly known. To address this issue within bees, we profiled caste-associated gene expression in female larvae of the intermediately eusocial bumblebee Bombus terrestris. In B. terrestris, female larvae experience a queen-dependent period during which their caste fate as adults is determined followed by a nutrition-sensitive period also potentially affecting caste fate but for which the evidence is weaker. We used mRNA-seq and qRT-PCR validation to isolate genes differentially expressed between each caste pathway in larvae at developmental stages before and after each of these periods. We show that differences in gene expression between caste pathways are small in totipotent larvae, then peak after the queen-dependent period. Relatively few novel (i.e., taxonomically-restricted) genes were differentially expressed between castes, though novel genes were significantly enriched in late-instar larvae in the worker pathway. We compared sets of caste-associated genes in B. terrestris with those reported from the advanced eusocial honeybee, Apis mellifera, and found significant but relatively low levels of overlap of gene lists between the two species. These results suggest both the existence of low numbers of shared toolkit genes and substantial divergence in caste-associated genes between Bombus and the advanced eusocial Apis since their last common eusocial ancestor.

Tobacco RNA-dependent RNA polymerase 1 affects the expression of defence-related genes in Nicotiana benthamiana upon Tomato leaf curl Gujarat virus infection.

MAIN CONCLUSION: RNA-dependent RNA polymerase 1 of Nicotiana tabacum modulates ToLCGV pathogenesis by influencing a number of defence-related genes in N. benthamiana plants. Key means of plants protecting themselves from the invading viruses is through RNA silencing. RNA-dependent RNA polymerase-1 (RDR1) is one of the crucial proteins of the RNA silencing pathway, which is induced after infection by viruses. RDR1 functions in the generation of small interfering RNAs (siRNAs) against the viral genome, thus it is antiviral in nature. Here, we used the transgenic Nicotiana benthamiana plant expressing N. tabacum NtRDR1 and observed reduced susceptibility towards Tomato leaf curl Gujarat virus (ToLCGV) infection compared to the wild-type N. benthamiana plants. To understand the reason for such reduced susceptibility, we prepared high-definition small RNA (sRNA) cDNA libraries from ToLCGV-infected wild-type N. benthamiana and NtRDR1 expressing N. benthamiana lines and carried out next-generation sequencing (NGS). We found that upon ToLCGV infection the majority of siRNAs generated from the host genome were of the 24 nucleotide (nt) class, while viral siRNAs (vsiRNAs) were of the 21-22-nt class, indicating that transcriptional gene silencing (TGS) is the major pathway for silencing of host genes while viral genes are silenced, predominantly, by post transcriptional gene silencing (PTGS) pathways. We estimated the changes in the expression of various defence-related genes, such as Constitutively Photomorphogenic-9 (COP9) signalosome (CSN) complex subunit-7, Pentatricopeptide repeat containing protein (PPRP), Laccase-3, Glutathione peroxidase-1 (GPX-1), Universal stress protein (USP) A-like protein, Heat shock transcription factor B4 (HSTF-B4), Auxin response factor-18 (ARF18), WRKY-6 and Short chain dehydrogenase reductase-3a. The differential expression of these genes might be linked with the enhanced tolerance of NtRDR1 N. benthamiana transgenic plants to ToLCGV. Our study suggests that reduced expression of subunit-7 of CSN complex and WRKY6, and increased expression of USPA-like protein might be linked with the reduced susceptibility of NtRDR1-transgenic N. benthamiana plants to ToLCGV.

LIMK, Cofilin 1 and actin dynamics involvement in fear memory processing.

Long-term memory has been associated with morphological changes in the brain, which in turn tightly correlate with changes in synaptic efficacy. Such plasticity is proposed to rely on dendritic spines as a neuronal canvas on which these changes can occur. Given the key role of actin cytoskeleton dynamics in spine morphology, major regulating factors of this process such as Cofilin 1 (Cfl1) and LIM kinase (LIMK), an inhibitor of Cfl1 activity, are prime molecular targets that may regulate dendritic plasticity. Using a contextual fear conditioning paradigm in mice, we found that pharmacological induction of depolymerization of actin filaments through the inhibition of LIMK causes an impairment in memory reconsolidation, as well as in memory consolidation. On top of that, Cfl1 activity is inhibited and its mRNA is downregulated in CA1 neuropil after re-exposure to the training context. Moreover, by pharmacological disruption of actin cytoskeleton dynamics, the process of memory extinction can either be facilitated or impaired. Our results lead to a better understanding of the role of LIMK, Cfl1 and actin cytoskeleton dynamics in the morphological and functional changes underlying the synaptic plasticity of the memory trace.

PAREsnip2: a tool for high-throughput prediction of small RNA targets from degradome sequencing data using configurable targeting rules.

Small RNAs (sRNAs) are short, non-coding RNAs that play critical roles in many important biological pathways. They suppress the translation of messenger RNAs (mRNAs) by directing the RNA-induced silencing complex to their sequence-specific mRNA target(s). In plants, this typically results in mRNA cleavage and subsequent degradation of the mRNA. The resulting mRNA fragments, or degradome, provide evidence for these interactions, and thus degradome analysis has become an important tool for sRNA target prediction. Even so, with the continuing advances in sequencing technologies, not only are larger and more complex genomes being sequenced, but also degradome and associated datasets are growing both in number and read count. As a result, existing degradome analysis tools are unable to process the volume of data being produced without imposing huge resource and time requirements. Moreover, these tools use stringent, non-configurable targeting rules, which reduces their flexibility. Here, we present a new and user configurable software tool for degradome analysis, which employs a novel search algorithm and sequence encoding technique to reduce the search space during analysis. The tool significantly reduces the time and resources required to perform degradome analysis, in some cases providing more than two orders of magnitude speed-up over current methods.

Acid sensing ion channel 2: A new potential player in the pathophysiology of multiple sclerosis.

Acid-sensing ion channels (ASICs) are proton-gated channels involved in multiple biological functions such as: pain modulation, mechanosensation, neurotransmission, and neurodegeneration. Earlier, we described the genetic association, within the Nuoro population, between Multiple Sclerosis (MS) and rs28936, located in ASIC2 3'UTR. Here we investigated the potential involvement of ASIC2 in MS inflammatory process. We induced experimental autoimmune encephalomyelitis (EAE) in wild-type (WT), knockout Asic1-/- and Asic2-/- mice and observed a significant reduction of clinical score in Asic1-/- mice and a significant reduction in the clinical score in Asic2-/- mice in a limited time window (i.e., at days 20-23 after immunization). Immunohistochemistry confirmed the reduction in adaptive immune cell infiltrates in the spinal cord of EAE Asic1-/- mice. Analysis of mechanical allodynia, showed a significant higher pain threshold in Asic2-/- mice under physiological conditions, before immunization, as compared to WT mice and Asic1-/- . A significant reduction in pain threshold was observed in all three strains of mice after immunization. More importantly, analysis of human autoptic brain tissue in MS and control samples showed an increase of ASIC2 mRNA in MS samples. Subsequently, in vitro luciferase reporter gene assays, showed that ASIC2 expression is under possible miRNA regulation, in a rs28936 allele-specific manner. Taken together, these findings suggest a potential role of ASIC2 in the pathophysiology of MS.

Comprehensive processing of high-throughput small RNA sequencing data including quality checking, normalization, and differential expression analysis using the UEA sRNA Workbench.

Recently, high-throughput sequencing (HTS) has revealed compelling details about the small RNA (sRNA) population in eukaryotes. These 20 to 25 nt noncoding RNAs can influence gene expression by acting as guides for the sequence-specific regulatory mechanism known as RNA silencing. The increase in sequencing depth and number of samples per project enables a better understanding of the role sRNAs play by facilitating the study of expression patterns. However, the intricacy of the biological hypotheses coupled with a lack of appropriate tools often leads to inadequate mining of the available data and thus, an incomplete description of the biological mechanisms involved. To enable a comprehensive study of differential expression in sRNA data sets, we present a new interactive pipeline that guides researchers through the various stages of data preprocessing and analysis. This includes various tools, some of which we specifically developed for sRNA analysis, for quality checking and normalization of sRNA samples as well as tools for the detection of differentially expressed sRNAs and identification of the resulting expression patterns. The pipeline is available within the UEA sRNA Workbench, a user-friendly software package for the processing of sRNA data sets. We demonstrate the use of the pipeline on a H. sapiens data set; additional examples on a B. terrestris data set and on an A. thaliana data set are described in the Supplemental Information A comparison with existing approaches is also included, which exemplifies some of the issues that need to be addressed for sRNA analysis and how the new pipeline may be used to do this.

Genomic responses to the socio-sexual environment in male Drosophila melanogaster exposed to conspecific rivals.

Socio-sexual environments have profound effects on fitness. Local sex ratios can alter the threat of sexual competition, to which males respond via plasticity in reproductive behaviors and ejaculate composition. In Drosophila melanogaster, males detect the presence of conspecific, same-sex mating rivals prior to mating using multiple, redundant sensory cues. Males that respond to rivals gain significant fitness benefits by altering mating duration and ejaculate composition. Here we investigated the underlying genome-wide changes involved. We used RNA-seq to analyze male transcriptomic responses 2, 26, and 50 h after exposure to rivals, a time period that was previously identified as encompassing the major facets of male responses to rivals. The results showed a strong early activation of multiple sensory genes in the head-thorax (HT), prior to the expression of any phenotypic differences. This gene expression response was reduced by 26 h, at the time of maximum phenotypic change, and shut off by 50 h. In the abdomen (A), fewer genes changed in expression and gene expression responses appeared to increase over time. The results also suggested that different sets of functionally equivalent genes might be activated in different replicates. This could represent a mechanism by which robustness is conferred upon highly plastic traits. Overall, our study reveals that mRNA-seq can identify subtle genomic signatures characteristic of flexible behavioral phenotypes.

The UEA sRNA Workbench (version 4.4): a comprehensive suite of tools for analyzing miRNAs and sRNAs.

Motivation: RNA interference, a highly conserved regulatory mechanism, is mediated via small RNAs (sRNA). Recent technical advances enabled the analysis of larger, complex datasets and the investigation of microRNAs and the less known small interfering RNAs. However, the size and intricacy of current data requires a comprehensive set of tools, able to discriminate the patterns from the low-level, noise-like, variation; numerous and varied suggestions from the community represent an invaluable source of ideas for future tools, the ability of the community to contribute to this software is essential. Results: We present a new version of the UEA sRNA Workbench, reconfigured to allow an easy insertion of new tools/workflows. In its released form, it comprises of a suite of tools in a user-friendly environment, with enhanced capabilities for a comprehensive processing of sRNA-seq data e.g. tools for an accurate prediction of sRNA loci (CoLIde) and miRNA loci (miRCat2), as well as workflows to guide the users through common steps such as quality checking of the input data, normalization of abundances or detection of differential expression represent the first step in sRNA-seq analyses. Availability and implementation: The UEA sRNA Workbench is available at: http://srna-workbench.cmp.uea.ac.uk. The source code is available at: https://github.com/sRNAworkbenchuea/UEA_sRNA_Workbench. Supplementary information: Supplementary data are available at Bioinformatics online.

Chromosomal-Level Assembly of the Asian Seabass Genome Using Long Sequence Reads and Multi-layered Scaffolding.

We report here the ~670 Mb genome assembly of the Asian seabass (Lates calcarifer), a tropical marine teleost. We used long-read sequencing augmented by transcriptomics, optical and genetic mapping along with shared synteny from closely related fish species to derive a chromosome-level assembly with a contig N50 size over 1 Mb and scaffold N50 size over 25 Mb that span ~90% of the genome. The population structure of L. calcarifer species complex was analyzed by re-sequencing 61 individuals representing various regions across the species' native range. SNP analyses identified high levels of genetic diversity and confirmed earlier indications of a population stratification comprising three clades with signs of admixture apparent in the South-East Asian population. The quality of the Asian seabass genome assembly far exceeds that of any other fish species, and will serve as a new standard for fish genomics.

Correction: Chromosomal-Level Assembly of the Asian Seabass Genome Using Long Sequence Reads and Multi-layered Scaffolding.

[This corrects the article DOI: 10.1371/journal.pgen.1005954.].

microRNAs associated with early neural crest development in Xenopus laevis.

BACKGROUND: The neural crest (NC) is a class of transitory stem cell-like cells unique to vertebrate embryos. NC cells arise within the dorsal neural tube where they undergo an epithelial to mesenchymal transition in order to migrate and differentiate throughout the developing embryo. The derivative cell types give rise to multiple tissues, including the craniofacial skeleton, peripheral nervous system and skin pigment cells. Several well-studied gene regulatory networks underpin NC development, which when disrupted can lead to various neurocristopathies such as craniofrontonasal dysplasia, DiGeorge syndrome and some forms of cancer. Small RNAs, such as microRNAs (miRNAs) are non-coding RNA molecules important in post-transcriptional gene silencing and critical for cellular regulation of gene expression. RESULTS: To uncover novel small RNAs in NC development we used high definition adapters and next generation sequencing of libraries derived from ectodermal explants of Xenopus laevis embryos induced to form neural and NC tissue. Ectodermal and blastula animal pole (blastula) stage tissues were also sequenced. We show that miR-427 is highly abundant in all four tissue types though in an isoform specific manner and we define a set of 11 miRNAs that are enriched in the NC. In addition, we show miR-301a and miR-338 are highly expressed in both the NC and blastula suggesting a role for these miRNAs in maintaining the stem cell-like phenotype of NC cells. CONCLUSION: We have characterised the miRNAs expressed in Xenopus embryonic explants treated to form ectoderm, neural or NC tissue. This has identified novel tissue specific miRNAs and highlighted differential expression of miR-427 isoforms.

Control of seminal fluid protein expression via regulatory hubs in Drosophila melanogaster.

Highly precise, yet flexible and responsive coordination of expression across groups of genes underpins the integrity of many vital functions. However, our understanding of gene regulatory networks (GRNs) is often hampered by the lack of experimentally tractable systems, by significant computational challenges derived from the large number of genes involved or from difficulties in the accurate identification and characterization of gene interactions. Here we used a tractable experimental system in which to study GRNs: the genes encoding the seminal fluid proteins that are transferred along with sperm (the 'transferome') in Drosophila melanogaster fruit flies. The products of transferome genes are core determinants of reproductive success and, to date, only transcription factors have been implicated in the modulation of their expression. Hence, as yet, we know nothing about the post-transcriptional mechanisms underlying the tight, responsive and precise regulation of this important gene set. We investigated this omission in the current study. We first used bioinformatics to identify potential regulatory motifs that linked the transferome genes in a putative interaction network. This predicted the presence of putative microRNA (miRNA) 'hubs'. We then tested this prediction, that post-transcriptional regulation is important for the control of transferome genes, by knocking down miRNA expression in adult males. This abolished the ability of males to respond adaptively to the threat of sexual competition, indicating a regulatory role for miRNAs in the regulation of transferome function. Further bioinformatics analysis then identified candidate miRNAs as putative regulatory hubs and evidence for variation in the strength of miRNA regulation across the transferome gene set. The results revealed regulatory mechanisms that can underpin robust, precise and flexible regulation of multiple fitness-related genes. They also help to explain how males can adaptively modulate ejaculate composition.

miRCat2: accurate prediction of plant and animal microRNAs from next-generation sequencing datasets.

MOTIVATION: MicroRNAs are a class of ∼21-22 nt small RNAs which are excised from a stable hairpin-like secondary structure. They have important gene regulatory functions and are involved in many pathways including developmental timing, organogenesis and development in eukaryotes. There are several computational tools for miRNA detection from next-generation sequencing datasets. However, many of these tools suffer from high false positive and false negative rates. Here we present a novel miRNA prediction algorithm, miRCat2. miRCat2 incorporates a new entropy-based approach to detect miRNA loci, which is designed to cope with the high sequencing depth of current next-generation sequencing datasets. It has a user-friendly interface and produces graphical representations of the hairpin structure and plots depicting the alignment of sequences on the secondary structure. RESULTS: We test miRCat2 on a number of animal and plant datasets and present a comparative analysis with miRCat, miRDeep2, miRPlant and miReap. We also use mutants in the miRNA biogenesis pathway to evaluate the predictions of these tools. Results indicate that miRCat2 has an improved accuracy compared with other methods tested. Moreover, miRCat2 predicts several new miRNAs that are differentially expressed in wild-type versus mutants in the miRNA biogenesis pathway. AVAILABILITY AND IMPLEMENTATION: miRCat2 is part of the UEA small RNA Workbench and is freely available from http://srna-workbench.cmp.uea.ac.uk/. CONTACT: v.moulton@uea.ac.uk or s.moxon@uea.ac.uk. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Ambient temperature regulates the expression of a small set of sRNAs influencing plant development through NF-YA2 and YUC2.

Plants substantially alter their developmental programme upon changes in the ambient temperature. The 21-24 nt small RNAs (sRNAs) are important gene expression regulators, which play a major role in development and adaptation. However, little is known about how the different sRNA classes respond to changes in the ambient temperature. We profiled the sRNA populations in four different tissues of Arabidopsis thaliana plants grown at 15°C, 21°C, and 27°C. We found that only a small fraction (0.6%) of the sRNA loci are ambient temperature-controlled. We identified thermoresponsive microRNAs and identified their target genes using degradome libraries. We verified that the target of the thermoregulated miR169, NF-YA2, is also ambient temperature-regulated. NF-YA2, as the component of the conserved transcriptional regulator NF-Y complex, binds the promoter of the flowering time regulator FT and the auxin biosynthesis gene YUC2. Other differentially expressed loci include thermoresponsive phased siRNA loci that target various auxin pathway genes and tRNA fragments. Furthermore, a temperature-dependent 24-nt heterochromatic siRNA locus in the promoter of YUC2 may contribute to the epigenetic regulation of auxin homeostasis. This holistic approach facilitated a better understanding of the role of different sRNA classes in ambient temperature adaptation of plants.

A non-canonical RNA silencing pathway promotes mRNA degradation in basal Fungi.

The increasing knowledge on the functional relevance of endogenous small RNAs (esRNAs) as riboregulators has stimulated the identification and characterization of these molecules in numerous eukaryotes. In the basal fungus Mucor circinelloides, an emerging opportunistic human pathogen, esRNAs that regulate the expression of many protein coding genes have been described. These esRNAs share common machinery for their biogenesis consisting of an RNase III endonuclease Dicer, a single Argonaute protein and two RNA-dependent RNA polymerases. We show in this study that, besides participating in this canonical dicer-dependent RNA interference (RNAi) pathway, the rdrp genes are involved in a novel dicer-independent degradation process of endogenous mRNAs. The analysis of esRNAs accumulated in wild type and silencing mutants demonstrates that this new rdrp-dependent dicer-independent regulatory pathway, which does not produce sRNA molecules of discrete sizes, controls the expression of target genes promoting the specific degradation of mRNAs by a previously unknown RNase. This pathway mainly regulates conserved genes involved in metabolism and cellular processes and signaling, such as those required for heme biosynthesis, and controls responses to specific environmental signals. Searching the Mucor genome for candidate RNases to participate in this pathway, and functional analysis of the corresponding knockout mutants, identified a new protein, R3B2. This RNase III-like protein presents unique domain architecture, it is specifically found in basal fungi and, besides its relevant role in the rdrp-dependent dicer-independent pathway, it is also involved in the canonical dicer-dependent RNAi pathway, highlighting its crucial role in the biogenesis and function of regulatory esRNAs. The involvement of RdRPs in RNA degradation could represent the first evolutionary step towards the development of an RNAi mechanism and constitutes a genetic link between mRNA degradation and post-transcriptional gene silencing.

Profile and functional analysis of small RNAs derived from Aspergillus fumigatus infected with double-stranded RNA mycoviruses.

BACKGROUND: Mycoviruses are viruses that naturally infect and replicate in fungi. Aspergillus fumigatus, an opportunistic pathogen causing fungal lung diseases in humans and animals, was recently shown to harbour several different types of mycoviruses. A well-characterised defence against virus infection is RNA silencing. The A. fumigatus genome encodes essential components of the RNA silencing machinery, including Dicer, Argonaute and RNA-dependent RNA polymerase (RdRP) homologues. Active silencing of double-stranded (ds)RNA and the generation of small RNAs (sRNAs) has been shown for several mycoviruses and it is anticipated that a similar mechanism will be activated in A. fumigatus isolates infected with mycoviruses. RESULTS: To investigate the existence and nature of A. fumigatus sRNAs, sRNA-seq libraries of virus-free and virus-infected isolates were created using Scriptminer adapters and compared. Three dsRNA viruses were investigated: Aspergillus fumigatus partitivirus-1 (AfuPV-1, PV), Aspergillus fumigatus chrysovirus (AfuCV, CV) and Aspergillus fumigatus tetramycovirus-1 (AfuTmV-1, NK) which were selected because they induce phenotypic changes such as coloration and sectoring. The dsRNAs of all three viruses, which included two conventionally encapsidated ones PV and CV and one unencapsidated example NK, were silenced and yielded characteristic vsiRNAs together with co-incidental silencing of host fungal genes which shared sequence homology with the viral genomes. CONCLUSIONS: Virus-derived sRNAs were detected and characterised in the presence of virus infection. Differentially expressed A. fumigatus microRNA-like (miRNA-like) sRNAs and small interfering RNAs (siRNAs) were detected and validated. Host sRNA loci which were differentially expressed as a result of virus infection were also identified. To our knowledge, this is the first study reporting the sRNA profiles of A. fumigatus isolates.