Liisa Heikkinen

Global microRNA Expression Profiling of Caenorhabditis elegans Parkinson's Disease Models

MicroRNAs (miRNAs) play an important role in human brain development and maintenance. To search for miRNAs that may be involved in the pathogenesis of Parkinsons disease (PD), we utilized miRNA microarrays to identify potential gene expression changes in 115 annotated miRNAs in PD-associated Caenorhabditis elegans models that either overexpress human A53T α-synuclein or have mutations within the vesicular catecholamine transporter (cat-1) or parkin (pdr-1) ortholog. Here, we show that 12 specific miRNAs are differentially regulated in the animals overexpressing α-synuclein, five in cat-1, and three in the pdr-1 mutants. The family of miR-64 and miR-65 are co-underexpressed in the α-synuclein transgenic and cat-1 strains, and members of let-7 family co-underexpressed in the α-synuclein and pdr-1 strains; mdl-1 and ptc-1 genes are target candidates for miR-64 and miR-65 and are overexpressed in α-synuclein transgenic as well as miR-64/65 (tm3711) knockout animals. These results indicate that miRNAs are differentially expressed in C. elegans PD models and suggest a role for these molecules in disease pathogenesis.

The Draft Genome and Transcriptome of Panagrellus redivivus Are Shaped by the Harsh Demands of a Free-Living Lifestyle

Nematodes compose an abundant and diverse invertebrate phylum with members inhabiting nearly every ecological niche. Panagrellus redivivus (the “microworm”) is a free-living nematode frequently used to understand the evolution of developmental and behavioral processes given its phylogenetic distance to Caenorhabditis elegans. Here we report the de novo sequencing of the genome, transcriptome, and small RNAs of P. redivivus. Using a combination of automated gene finders and RNA-seq data, we predict 24,249 genes and 32,676 transcripts. Small RNA analysis revealed 248 microRNA (miRNA) hairpins, of which 63 had orthologs in other species. Fourteen miRNA clusters containing 42 miRNA precursors were found. The RNA interference, dauer development, and programmed cell death pathways are largely conserved. Analysis of protein family domain abundance revealed that P. redivivus has experienced a striking expansion of BTB domain-containing proteins and an unprecedented expansion of the cullin scaffold family of proteins involved in multi-subunit ubiquitin ligases, suggesting proteolytic plasticity and/or tighter regulation of protein turnover. The eukaryotic release factor protein family has also been dramatically expanded and suggests an ongoing evolutionary arms race with viruses and transposons. The P. redivivus genome provides a resource to advance our understanding of nematode evolution and biology and to further elucidate the genomic architecture leading to free-living lineages, taking advantage of the many fascinating features of this worm revealed by comparative studies.

Chronic Ethanol Exposure Increases Cytochrome P‐450 and Decreases Activated in Blocked Unfolded Protein Response Gene Family Transcripts in Caenorhabditis elegans

Ethanol is a widely consumed and rapidly absorbed toxin. While the physiological effects of ethanol consumption are well known, the underlying biochemical and molecular changes at the gene expression level in whole animals remain obscure. We exposed the model organism Caenorhabditis elegans to 0.2 M ethanol from the embryo to L4 larva stage and assayed gene expression changes in whole animals using RNA‐Seq and quantitative real‐time PCR. We observed gene expression changes in 1122 genes (411 up, 711 down). Cytochrome P‐450 (CYP) gene family members (12 of 78) were upregulated, whereas activated in blocked unfolded protein response (ABU) (7 of 15) were downregulated. Other detoxification gene family members were also regulated including four glutathione‐S‐transferases and three flavin monooxygenases. The results presented show specific gene expression changes following chronic ethanol exposure in C. elegans that indicate both persistent upregulation of detoxification response genes and downregulation of endoplasmic reticulum stress pathway genes.

Selective MicroRNA-Offset RNA Expression in Human Embryonic Stem Cells

Small RNA molecules, including microRNAs (miRNAs), play critical roles in regulating pluripotency, proliferation and differentiation of embryonic stem cells. miRNA-offset RNAs (moRNAs) are similar in length to miRNAs, align to miRNA precursor (pre-miRNA) loci and are therefore believed to derive from processing of the pre-miRNA hairpin sequence. Recent next generation sequencing (NGS) studies have reported the presence of moRNAs in human neurons and cancer cells and in several tissues in mouse, including pluripotent stem cells. In order to gain additional knowledge about human moRNAs and their putative development-related expression, we applied NGS of small RNAs in human embryonic stem cells (hESCs) and fibroblasts. We found that certain moRNA isoforms are notably expressed in hESCs from loci coding for stem cell-selective or cancer-related miRNA clusters. In contrast, we observed only sparse moRNAs in fibroblasts. Consistent with earlier findings, most of the observed moRNAs derived from conserved loci and their expression did not appear to correlate with the expression of the adjacent miRNAs. We provide here the first report of moRNAs in hESCs, and their expression profile in comparison to fibroblasts. Moreover, we expand the repertoire of hESC miRNAs. These findings provide an expansion on the known repertoire of small non-coding RNA contents in hESCs.

Identification of phylogenetically conserved sequence motifs in microRNA 5' flanking sites from C. elegans and C. briggsae

BackgroundMicroRNAs (miRNAs) are small, noncoding RNA molecules that act as post-transcriptional regulators of gene expression. Studies concerning transcriptional regulation of miRNAs have so far concentrated on those located within the intergenic region of the genome and the search for putative promoters, thus leaving open the question of the existence of possible regulatory elements common to all miRNAs including those located in introns of protein coding genes.ResultsIn this study, we initially searched for motifs occurring in the area 1000 bp upstream from all miRNAs independent of their genomic location. We discovered a previously unknown sequence motif GANNNNGA that displayed a conserved distribution in the nematode worms Caenorhabditis elegans and Caenorhabditis briggsae. This motif had a peak occurrence at 500 bp upstream, with a sharp drop-off toward the miRNA start site. Further analysis indicated that this motif was locally restricted and not enriched 1000–5000 bp upstream or 0–2000 bp downstream of the miRNA start site. In addition, this motif was observed to be most abundant in the upstream sequences of two important miRNAs, mir-1 and mir-124. This abundance was also conserved in phylogenetically distant species including human and mouse.ConclusionThe results show that the motif GANNNNGA is conserved close to miRNA precursor start sites, suggesting that it may be involved in miRNA sequence recognition or regulation. This data provides important knowledge for the identification and computational prediction of miRNA sequences.

Methylmercury exposure increases lipocalin related (lpr) and decreases activated in blocked unfolded protein response (abu) genes and specific miRNAs in Caenorhabditis elegans.

Methylmercury (MeHg) is a persistent environmental and dietary contaminant that causes serious adverse developmental and physiologic effects at multiple cellular levels. In order to understand more fully the consequences of MeHg exposure at the molecular level, we profiled gene and miRNA transcripts from the model organism Caenorhabditis elegans. Animals were exposed to MeHg (10 μM) from embryo to larval 4 (L4) stage and RNAs were isolated. RNA-seq analysis on the Illumina platform revealed 541 genes up- and 261 genes down-regulated at a cutoff of 2-fold change and false discovery rate-corrected significance q < 0.05. Among the up-regulated genes were those previously shown to increase under oxidative stress conditions including hsp-16.11 (2.5-fold), gst-35 (10.1-fold), and fmo-2 (58.5-fold). In addition, we observed up-regulation of 6 out of 7 lipocalin related (lpr) family genes and down regulation of 7 out of 15 activated in blocked unfolded protein response (abu) genes. Gene Ontology enrichment analysis highlighted the effect of genes related to development and organism growth. miRNA-seq analysis revealed 6-8 fold down regulation of mir-37-3p, mir-41-5p, mir-70-3p, and mir-75-3p. Our results demonstrate the effects of MeHg on specific transcripts encoding proteins in oxidative stress responses and in ER stress pathways. Pending confirmation of these transcript changes at protein levels, their association and dissociation characteristics with interaction partners, and integration of these signals, these findings indicate broad and dynamic mechanisms by which MeHg exerts its harmful effects.

Evolutionary conservation and function of the human embryonic stem cell specific miR-302/367 cluster.

miRNA clusters define a group of related miRNAs closely localized in the genome with an evolution that remains poorly understood. The miR-302/367 cluster represents a single polycistronic transcript that produces five precursor miRNAs. The cluster is highly expressed and essential for maintenance of human embryonic stem cells. We found the cluster to be highly conserved and present in most mammals. In primates, seed sequence and miRNA structure are conserved, but inter-precursor sequences are evolving. Insertions of new miRNAs, deletions of individual miRNAs, and a cluster duplication observed in different species suggest an actively evolving cluster. Core transcriptional machinery consisting of NANOG and OCT-4 transcription factors that define stem cells are present upstream of the miR-302/367 cluster. Interestingly, we found the miR-302/367 cluster flanking region to be enriched as a target site of other miRNAs suggesting a mechanism for feedback control. Analysis of miR-302 and miR-367 targets demonstrated concordance of gene set enrichment groups at high gene ontology levels. This cluster also expresses isomiRs providing another means of establishing sequence diversity. Finally, using three different kidney tumor datasets, we observed consistent expression of miR-302 family members in normal tissue while adjacent tumor tissue showed a significant lack of expression. Clustering expression levels of miR-302 validated target genes showed a significant correlation between miR-302/367 cluster miRNAs and a subset of validated gene targets in healthy and adjacent tumor tissues. Taken together, our data show a highly conserved and still evolving miRNA cluster that may have additional unrecognized functions.

De novo transcriptome assembly and developmental mode specific gene expression of Pygospio elegans

Species with multiple different larval developmental modes are interesting models for the study of mechanisms underlying developmental mode transitions and life history evolution. Pygospio elegans, a small, tube‐dwelling polychaete worm commonly found in estuarine and marine habitats around the northern hemisphere, is one species with variable developmental modes. To provide new genomic resources for studying P. elegans and to address the differences in gene expression between individuals producing offspring with different larval developmental modes, we performed whole transcriptome Illumina RNA sequencing of adult worms from two populations and prepared a de novo assembly of the P. elegans transcriptome. The transcriptome comprises 66,233 unigenes, of which 33,807 contain predicted coding sequences, 26,448 have at least one functional annotation, and 3,076 are classified as putative long non‐coding RNAs. We found more than 8,000 unigenes significantly differentially expressed between adult worms from populations producing either planktonic or benthic larvae. This comprehensive transcriptome resource for P. elegans adds to the available genomic data for annelids and can be used to uncover mechanisms allowing developmental variation in this and potentially other marine invertebrate species.

DNA Methylation and Potential for Epigenetic Regulation in Pygospio elegans.

Transitions in developmental mode are common evolutionarily, but how and why they occur is not understood. Developmental mode describes larval phenotypes, including morphology, ecology and behavior of larvae, which typically are generalized across different species. The polychaete worm Pygospio elegans is one of few species polymorphic in developmental mode, with multiple larval phenotypes, providing a possibility to examine the potential mechanisms allowing transitions in developmental mode. We investigated the presence of DNA methylation in P. elegans, and, since maternal provisioning is a key factor determining eventual larval phenotype, we compared patterns of DNA methylation in females during oogenesis in this species. We demonstrate that intragenic CpG site DNA methylation and many relevant genes necessary for DNA methylation occur in P. elegans. Methylation-sensitive AFLP analysis showed that gravid females with offspring differing in larval developmental mode have significantly different methylation profiles and that the females with benthic larvae and non-reproductive females from the same location also differ in their epigenetic profiles. Analysis of CpG sites in transcriptome data supported our findings of DNA methylation in this species and showed that CpG observed/expected ratios differ among females gravid with embryos destined to different developmental modes. The differences in CpG site DNA methylation patterns seen among the samples suggest a potential for epigenetic regulation of gene expression (through DNA methylation) in this species.

miRToolsGallery: a tag-based and rankable microRNA bioinformatics resources database portal

Abstract Hundreds of bioinformatics tools have been developed for MicroRNA (miRNA) investigations including those used for identification, target prediction, structure and expression profile analysis. However, finding the correct tool for a specific application requires the tedious and laborious process of locating, downloading, testing and validating the appropriate tool from a group of nearly a thousand. In order to facilitate this process, we developed a novel database portal named miRToolsGallery. We constructed the portal by manually curating > 950 miRNA analysis tools and resources. In the portal, a query to locate the appropriate tool is expedited by being searchable, filterable and rankable. The ranking feature is vital to quickly identify and prioritize the more useful from the obscure tools. Tools are ranked via different criteria including the PageRank algorithm, date of publication, number of citations, average of votes and number of publications. miRToolsGallery provides links and data for the comprehensive collection of currently available miRNA tools with a ranking function which can be adjusted using different criteria according to specific requirements. Database URL: http://www.mirtoolsgallery.org