Microbial impact to environmental toxicants Ni(II) and Co(II): Joint toxicity and cellular response in Paramecium.

Cobalt (Co) and Nickel (Ni) are increasingly found in our environment. We analysed their combined toxicity and uptake mechanisms in the early food chain by studying bacteria and the bacterivorous ciliate Paramecium as a primary consumer. We exposed both species to these metals to measure the toxicity, uptake and transfer of metals from bacteria to Paramecium. We found that Ni is more toxic than Co, and that toxicity increases for both metals when (i) food bacteria are absent and (ii) both metals are applied in combination. The cellular content in bacteria after exposure shows a concentration dependent bias for either Ni or Co. Comparing single treatment and joint exposure, bacteria show increased levels of both metals when these are both exposed. To imitate the basic level of the food chain, we fed these bacteria to paramecia. The cellular content shows a similar ratio of Nickel and Cobalt as in food bacteria. This is different to the direct application of both metals to paramecia, where Cobalt is enriched over Nickel. This indicates that bacteria can selectively pre-accumulate metals for introduction into the food chain. We also analysed the transcriptomic response of Paramecium to sublethal doses of Nickel and Cobalt to gain insight into their toxicity mechanisms. Gene ontology (GO) analysis indicates common deregulated pathways, such as ammonium transmembrane transport and ubiquitine-associated protein degradation. Many redox-related genes also show deregulation of gene expression, indicating cellular adaptation to increased RONS stress. This suggests that both metals may also target the same cellular pathways and this is consistent with the increased toxicity of both metals when used together. Our data reveal complex ecotoxicological pathways for these metals and highlights the different parameters for their fate in the ecosystem, in the food chain and their ecotoxicological risk after environmental contamination.

Species-Specific Duplication of Surface Antigen Genes in Paramecium.

Paramecium is a free-living ciliate that undergoes antigenic variation and still the functions of these variable surface antigen coats in this non-pathogenic ciliate remain elusive. Only a few surface antigen genes have been described, mainly in the two model species P. tetraurelia strain 51 and P. primaurelia strain 156. Given the lack of suitable sequence data to allow for phylogenetics and deeper sequence comparisons, we screened the genomes of six different Paramecium species for serotype genes and isolated 548 candidates. Our approach identified the subfamilies of the isogenes of individual serotypes that were mostly represented by intrachromosomal gene duplicates. These showed different duplication levels, and chromosome synteny suggested rather young duplication events after the emergence of the P. aurelia species complex, indicating a rapid evolution of surface antigen genes. We were able to identify the different subfamilies of the surface antigen genes with internal tandem repeats, which showed consensus motifs across species. The individual isogene families showed additional consensus motifs, indicating that the selection pressure holds individual amino acids constant in these repeats. This may be a hint of the receptor function of these antigens rather than a presentation of random epitopes, generating the variability of these surface molecules.

Manganese-Induced Toxicity in C. elegans: What Can We Learn from the Transcriptome?

Manganese (Mn) is an essential ubiquitous transition metal and, when occupationally or environmentally overexposed, a well-known risk factor for several neurological pathologies. However, the molecular mechanisms underlying Mn-induced neurotoxicity are largely unknown. In this study, addressing RNA-Seq analysis, bioavailability and survival assays, key pathways of transcriptional responses to Mn overexposure were investigated in the model organism Caenorhabditis elegans (C. elegans), providing insights into the Mn-induced cellular stress and damage response. Comparative transcriptome analyses identified a large number of differentially expressed genes (DEGs) in nematodes exposed to MnCl2, and functional annotation suggested oxidative nucleotide damage, unfolded protein response and innate immunity as major damage response pathways. Additionally, a time-dependent increase in the transcriptional response after MnCl2 exposure was identified by means of increased numbers of DEGs, indicating a time-dependent response and activation of the stress responses in Mn neurotoxicity. The data provided here represent a powerful transcriptomic resource in the field of Mn toxicity, and therefore, this study provides a useful basis for further planning of targeted mechanistic studies of Mn-induced neurotoxicity that are urgently needed in the face of increasing industrially caused environmental pollution with Mn.

Two Piwis with Ago-like functions silence somatic genes at the chromatin level.

Most sRNA biogenesis mechanisms involve either RNAse III cleavage or ping-pong amplification by different Piwi proteins harbouring slicer activity. Here, we follow the question why the mechanism of transgene-induced silencing in the ciliate Paramecium needs both Dicer activity and two Ptiwi proteins. This pathway involves primary siRNAs produced from non-translatable transgenes and secondary siRNAs from targeted endogenous loci. Our data does not indicate any signatures from ping-pong amplification but Dicer cleavage of long dsRNA. Ptiwi13 and 14 prefer different sub-cellular localizations and different preferences for primary and secondary siRNAs but do not load them mutually exclusive. Both Piwis enrich for antisense RNAs and show a general preference for uridine-rich sRNAs along the entire sRNA length. In addition, Ptiwi14-loaded siRNAs show a 5´-U signature. Our data indicates both Ptiwis and 2´-O-methylation contributing to strand selection of Dicer cleaved siRNAs. This unexpected function of the two distinct vegetative Piwis extends the increasing knowledge of the diversity of Piwi functions in diverse silencing pathways. We describe an unusual mode of action of Piwi proteins extending not only the great variety of Piwi-associated RNAi pathways but moreover raising the question whether this could have been the primordial one.

The complete mitochondrial genome of the photosymbiotic sea slug Berghia stephanieae (Valdés, 2005) (Gastropoda, Nudibranchia).

Berghia stephanieae (Nudibranchia, Cladobranchia) is a photosymbiotic sea slug that feeds exclusively on sea anemones from the genus Exaiptasia. It then specifically incorporates dinoflagellates belonging to the Symbiodiniaceae obtained from their prey. Here, we present the complete mitochondrial genome sequence of B. stephanieae combining Oxford Nanopore long read and Illumina short-read sequencing data. The mitochondrial genome has a total length of 14,786 bp, it contains the 13 protein-encoding genes, 23 tRNAs, and two rRNAs and is similar to other nudibranchs except for the presence of a duplicated tRNA-Ser 1.

Dual-Seq reveals genome and transcriptome of Caedibacter taeniospiralis, obligate endosymbiont of Paramecium.

Interest in host-symbiont interactions is continuously increasing, not only due to the growing recognition of the importance of microbiomes. Starting with the detection and description of novel symbionts, attention moves to the molecular consequences and innovations of symbioses. However, molecular analysis requires genomic data which is difficult to obtain from obligate intracellular and uncultivated bacteria. We report the identification of the Caedibacter genome, an obligate symbiont of the ciliate Paramecium. The infection does not only confer the host with the ability to kill other cells but also renders them immune against this effect. We obtained the C. taeniospiralis genome and transcriptome by dual-Seq of DNA and RNA from infected paramecia. Comparison of codon usage and expression level indicates that genes necessary for a specific trait of this symbiosis, i.e. the delivery of an unknown toxin, result from horizontal gene transfer hinting to the relevance of DNA transfer for acquiring new characters. Prediction of secreted proteins of Caedibacter as major agents of contact with the host implies, next to several toxin candidates, a rather uncharacterized secretome which appears to be highly adapted to this symbiosis. Our data provides new insights into the molecular establishment and evolution of this obligate symbiosis and for the pathway characterization of toxicity and immunity.

Comparative Analysis of Biochemical Biases by Ligation- and Template-Switch-Based Small RNA Library Preparation Protocols.

BACKGROUND: Small RNAs are key players in the regulation of gene expression and differentiation. However, many different classes of small RNAs (sRNAs) have been described with distinct biogenesis pathways and, as a result, with different biochemical properties. To analyze sRNAs by deep sequencing, complementary DNA synthesis requires manipulation of the RNA molecule itself. Thus, enzymatic activities during library preparation bias the library content owing to biochemical criteria. METHODS: We compared 4 different manipulations of RNA for library preparation: (a) a ligation-based procedure allowing only 5'-mono-phosphorylated RNA to enter the library, (b) a ligation-based procedure allowing additional 5'-triphosphates and Cap structures, (c) a ligation-independent, template-switch-based library preparation, and (d) a template-switch-based library preparation allowing 3'-phosphorylated RNAs to enter the library. RESULTS: Our data show large differences between ligation-dependent and ligation-independent libraries in terms of their preference for individual sRNA classes such as microRNAs (miRNAs), Piwi-interacting RNAs (piRNAs), and transfer RNA fragments. Moreover, the miRNA composition is different between both procedures, and more microRNA isoforms (isomiRs) can be identified after pyrophosphatase treatment. piRNAs are enriched in template-switch libraries, and this procedure apparently includes more different RNA species. CONCLUSIONS: Our data indicate that miRNAomics from both methods will hardly be comparable. Ligation-based libraries enrich for canonical miRNAs, which thus may be suitable methods for miRNAomics. Template-switch libraries contain increased numbers and different compositions of fragments and long RNAs. Following different interests for other small RNA species, ligation-independent libraries appear to show a more realistic sRNA landscape with lower bias against biochemical modifications.

Exogenous RNAi mechanisms contribute to transcriptome adaptation by phased siRNA clusters in Paramecium.

Extensive research has characterized distinct exogenous RNAi pathways interfering in gene expression during vegetative growth of the unicellular model ciliate Paramecium. However, role of RNAi in endogenous transcriptome regulation, and environmental adaptation is unknown. Here, we describe the first genome-wide profiling of endogenous sRNAs in context of different transcriptomic states (serotypes). We developed a pipeline to identify, and characterize 2602 siRNA producing clusters (SRCs). Our data show no evidence that SRCs produce miRNAs, and in contrast to other species, no preference for strand specificity of siRNAs. Interestingly, most SRCs overlap coding genes and a separate group show siRNA phasing along the entire open reading frame, suggesting that the mRNA transcript serves as a source for siRNAs. Integrative analysis of siRNA abundance and gene expression levels revealed surprisingly that mRNA and siRNA show negative as well as positive associations. Two RNA-dependent RNA Polymerase mutants, RDR1 and RDR2, show a drastic loss of siRNAs especially in phased SRCs accompanied with increased mRNA levels. Importantly, most SRCs depend on both RDRs, reminiscent to primary siRNAs in the RNAi against exogenous RNA, indicating mechanistic overlaps between exogenous and endogenous RNAi contributing to flexible transcriptome adaptation.

The sncRNA Zoo: a repository for circulating small noncoding RNAs in animals.

The repertoire of small noncoding RNAs (sncRNAs), particularly miRNAs, in animals is considered to be evolutionarily conserved. Studies on sncRNAs are often largely based on homology-based information, relying on genomic sequence similarity and excluding actual expression data. To obtain information on sncRNA expression (including miRNAs, snoRNAs, YRNAs and tRNAs), we performed low-input-volume next-generation sequencing of 500 pg of RNA from 21 animals at two German zoological gardens. Notably, none of the species under investigation were previously annotated in any miRNA reference database. Sequencing was performed on blood cells as they are amongst the most accessible, stable and abundant sources of the different sncRNA classes. We evaluated and compared the composition and nature of sncRNAs across the different species by computational approaches. While the distribution of sncRNAs in the different RNA classes varied significantly, general evolutionary patterns were maintained. In particular, miRNA sequences and expression were found to be even more conserved than previously assumed. To make the results available for other researchers, all data, including expression profiles at the species and family levels, and different tools for viewing, filtering and searching the data are freely available in the online resource ASRA (Animal sncRNA Atlas) at https://www.ccb.uni-saarland.de/asra/.

Next Generation Sequencing Analysis of Total Small Noncoding RNAs from Low Input RNA from Dried Blood Sampling.

Circulating miRNAs are favored for biomarker candidates as they can reflect tissue specific miRNA dysregulation in disease contexts. Moreover, they have the additional advantage that they can be monitored in a minimally invasive manner. Blood-borne miRNAs are therefore currently characterized to identify, describe, and validate their potential suitability as biomarkers; however, sampling and as well miRNA detection methods limit these studies in terms of sensitivity but also practicability in clinical, at-home, or low-resource sampling of high-quality circulating RNA samples. We describe here a novel and innovative method of circulating RNA microsampling from minimal volume dried blood samples with direct enrichment for small RNA fractions in combination with ligation free library preparation. We evaluated crucial parameters for efficient library preparation from low RNA inputs of 50 pg for efficient dissection not only of miRNAs but also isomiRs, piRNAs, and lincRNAs. We compared these data to classical microarrays and characterize the technical reproducibility and its sensitivity. We demonstrate and evaluate a method for easy low resource sampling and NGS analysis of circulating RNAs providing a powerful tool for massive cohort and remote patient monitoring.

Environmental Temperature Controls Accumulation of Transacting siRNAs Involved in Heterochromatin Formation.

Genes or alleles can interact by small RNAs in a homology dependent manner meaning that short interfering (siRNAs) can act in trans at the chromatin level producing stable and heritable silencing phenotypes. Because of the puzzling data on endogenous paramutations, their impact contributing to adaptive evolution in a Lamarckian manner remains unknown. An increasing number of studies characterizes the underlying siRNA accumulation pathways using transgene experiments. Also in the ciliate Paramecium tetraurelia, we induce trans silencing on the chromatin level by injection of truncated transgenes. Here, we characterize the efficiency of this mechanism at different temperatures showing that silencing of the endogenous genes is temperature dependent. Analyzing different transgene constructs at different copy numbers, we dissected whether silencing efficiency is due to varying precursor RNAs or siRNA accumulation. Our data shows that silencing efficiency correlates with more efficient accumulation of primary siRNAs at higher temperatures rather than higher expression of precursor RNAs. Due to higher primary levels, secondary siRNAs also show temperature dependency and interestingly increase their relative proportion to primary siRNAs. Our data shows that efficient trans silencing on the chromatin level in P. tetraurelia depends on environmental parameters, thus being an important epigenetic factor limiting regulatory effects of siRNAs.