Computational discovery of hidden breaks in 28S ribosomal RNAs across eukaryotes and consequences for RNA Integrity Numbers

Trophic Position of the White Worm (Enchytraeus albidus) in the Context of Digestive Enzyme Genes Revealed by Transcriptomics Analysis

To assess the impact of Enchytraeidae (potworms) on the functioning of the decomposer system, knowledge of the feeding preferences of enchytraeid species is required. Different food preferences can be explained by variations in enzymatic activities among different enchytraeid species, as there are no significant differences in the morphology or anatomy of their alimentary tracts. However, it is crucial to distinguish between the contribution of microbial enzymes and the animal's digestive capacity. Here, we computationally analyzed the endogenous digestive enzyme genes in Enchytraeus albidus. The analysis was based on RNA-Seq of COI-monohaplotype culture (PL-A strain) specimens, utilizing transcriptome profiling to determine the trophic position of the species. We also corroborated the results obtained using transcriptomics data from genetically heterogeneous freeze-tolerant strains. Our results revealed that E. albidus expresses a wide range of glycosidases, including GH9 cellulases and a specific digestive SH3b-domain-containing i-type lysozyme, previously described in the earthworm Eisenia andrei. Therefore, E. albidus combines traits of both primary decomposers (primary saprophytophages) and secondary decomposers (sapro-microphytophages/microbivores) and can be defined as an intermediate decomposer. Based on assemblies of publicly available RNA-Seq reads, we found close homologs for these cellulases and i-type lysozymes in various clitellate taxa, including Crassiclitellata and Enchytraeidae.

Comprehensive map of ribosomal 2'-O-methylation and C/D box snoRNAs in Drosophila melanogaster

During their maturation, ribosomal RNAs (rRNAs) are decorated by hundreds of chemical modifications that participate in proper folding of rRNA secondary structures and therefore in ribosomal function. Along with pseudouridine, methylation of the 2'-hydroxyl ribose moiety (Nm) is the most abundant modification of rRNAs. The majority of Nm modifications in eukaryotes are placed by Fibrillarin, a conserved methyltransferase belonging to a ribonucleoprotein complex guided by C/D box small nucleolar RNAs (C/D box snoRNAs). These modifications impact interactions between rRNAs, tRNAs and mRNAs, and some are known to fine tune translation rates and efficiency. In this study, we built the first comprehensive map of Nm sites in Drosophila melanogaster rRNAs using two complementary approaches (RiboMethSeq and Nanopore direct RNA sequencing) and identified their corresponding C/D box snoRNAs by whole-transcriptome sequencing. We de novo identified 61 Nm sites, from which 55 are supported by both sequencing methods, we validated the expression of 106 C/D box snoRNAs and we predicted new or alternative rRNA Nm targets for 31 of them. Comparison of methylation level upon different stresses show only slight but specific variations, indicating that this modification is relatively stable in D. melanogaster. This study paves the way to investigate the impact of snoRNA-mediated 2'-O-methylation on translation and proteostasis in a whole organism.

Detection of Developmental Asexual Stage-Specific RNA Editing Events in Plasmodium falciparum 3D7 Malaria Parasite

Transcriptional variation has been studied but post-transcriptional modification due to RNA editing has not been investigated in Plasmodium. We investigated developmental stage-specific RNA editing in selected genes in Plasmodium falciparum 3D7. We detected extensive amination- and deamination-type RNA editing at 8, 16, 24, 32, 40, and 46 h in tightly synchronized Plasmodium. Most of the editing events were observed in 8 and 16 h ring-stage parasites. Extensive A-to-G deamination-type editing was detected more during the 16 h ring stage (25%) than the 8 h ring stage (20%). Extensive U-to-C amination-type editing was detected more during the 16 h ring stage (31%) than the 8 h ring stage (22%). In 28S, rRNA editing converted the loop structure to the stem structure. The hemoglobin binding activity of PF3D7_0216900 was also altered due to RNA editing. Among the expressed 28S rRNA genes, PF3D7_0532000 and PF3D7_0726000 expression was higher. Increased amounts of the transcripts of these two genes were found, particularly PF3D7_0726000 in the ring stage and PF3D7_0532000 in the trophozoite and schizont stages. Adenosine deaminase (ADA) expression did not correlate with the editing level. This first experimental report of RNA editing will help to identify the editing machinery that might be useful for antimalarial drug discovery and malaria control.

Cephalopod-omics: Emerging Fields and Technologies in Cephalopod Biology

Few animal groups can claim the level of wonder that cephalopods instill in the minds of researchers and the general public. Much of cephalopod biology, however, remains unexplored: the largest invertebrate brain, difficult husbandry conditions, and complex (meta-)genomes, among many other things, have hindered progress in addressing key questions. However, recent technological advancements in sequencing, imaging, and genetic manipulation have opened new avenues for exploring the biology of these extraordinary animals. The cephalopod molecular biology community is thus experiencing a large influx of researchers, emerging from different fields, accelerating the pace of research in this clade. In the first post-pandemic event at the Cephalopod International Advisory Council (CIAC) conference in April 2022, over 40 participants from all over the world met and discussed key challenges and perspectives for current cephalopod molecular biology and evolution. Our particular focus was on the fields of comparative and regulatory genomics, gene manipulation, single-cell transcriptomics, metagenomics, and microbial interactions. This article is a result of this joint effort, summarizing the latest insights from these emerging fields, their bottlenecks, and potential solutions. The article highlights the interdisciplinary nature of the cephalopod-omics community and provides an emphasis on continuous consolidation of efforts and collaboration in this rapidly evolving field.

Unveiling the Protein Components of the Secretory-Venom Gland and Venom of the Scorpion Centruroides possanii (Buthidae) through Omic Technologies

Centruroides possanii is a recently discovered species of "striped scorpion" found in Mexico. Certain species of Centruroides are known to be toxic to mammals, leading to numerous cases of human intoxications in the country. Venom components are thought to possess therapeutic potential and/or biotechnological applications. Hence, obtaining and analyzing the secretory gland transcriptome and venom proteome of C. possanii is relevant, and that is what is described in this communication. Since this is a newly described species, first, its LD50 to mice was determined and estimated to be 659 ng/g mouse weight. Using RNA extracted from this species and preparing their corresponding cDNA fragments, a transcriptome analysis was obtained on a Genome Analyzer (Illumina) using the 76-base pair-end sequencing protocol. Via high-throughput sequencing, 19,158,736 reads were obtained and ensembled in 835,204 sequences. Of them, 28,399 transcripts were annotated with Pfam. A total of 244 complete transcripts were identified in the transcriptome of C. possanii. Of these, 109 sequences showed identity to toxins that act on ion channels, 47 enzymes, 17 protease inhibitors (PINs), 11 defense peptides (HDPs), and 60 in other components. In addition, a sample of the soluble venom obtained from this scorpion was analyzed using an Orbitrap Velos apparatus, which allowed for identification by liquid chromatography followed by mass spectrometry (LC-MS/MS) of 70 peptides and proteins: 23 toxins, 27 enzymes, 6 PINs, 3 HDPs, and 11 other components. Until now, this work has the highest number of scorpion venom components identified through omics technologies. The main novel findings described here were analyzed in comparison with the known data from the literature, and this process permitted some new insights in this field.

Transcriptional and metabolomic investigation of the stress response in snow crab during simulated transport condition (Chionoecetes opilio)

The molecular mechanisms underlying the stress response are poorly described in crustaceans. This includes the snow crab (Chionoecetes opilio), a commercially important stenotherm species distributed throughout the northern hemisphere. A better understanding of the stress response in C. opilio is desperately needed for commercial and conservation purposes. The purpose of this study was to investigate the transcriptional and metabolomic response of C. opilio exposed to stressors. Crabs were randomly assigned to 24 or 72 h treatment groups where they were exposed to conditions simulating live transport (handling and air exposure). A control group was kept in cold (2 °C) and well‑oxygenated saltwater. The hepatopancreas of the crabs was sampled to perform RNA-sequencing and high-performance chemical isotope labeling metabolomics. Differential gene expression analyses showed that classic crustaceans' stress markers, such as crustacean hyperglycemic hormones and heat shock proteins, were overexpressed in response to stressors. Tyrosine decarboxylase was also up-regulated in stressed crabs, suggesting an implication of the catecholamines tyramine and octopamine in the stress response. Deregulated metabolites revealed that low oxygen was an important trigger in the stress response as intermediate metabolites of the tricarboxylic acid cycle (TCA) accumulated. Lactate, which accumulated unevenly between crabs could potentially be used to predict mortality. This study provides new information on how stressors affect crustaceans and provides a basis for the development of stress markers in C. opilio.

De novo genome assembly and genome skims reveal LTRs dominate the genome of a limestone endemic Mountainsnail (Oreohelix idahoensis)

BACKGROUND

Calcareous outcrops, rocky areas composed of calcium carbonate (CaCO₃), often host a diverse, specialized, and threatened biomineralizing fauna. Despite the repeated evolution of physiological and morphological adaptations to colonize these mineral rich substrates, there is a lack of genomic resources for calcareous rock endemic species. This has hampered our ability to understand the genomic mechanisms underlying calcareous rock specialization and manage these threatened species.

RESULTS

Here, we present a new draft genome assembly of the threatened limestone endemic land snail Oreohelix idahoensis and genome skim data for two other Oreohelix species. The O. idahoensis genome assembly (scaffold N50: 404.19 kb; 86.6% BUSCO genes) is the largest (~ 5.4 Gb) and most repetitive mollusc genome assembled to date (85.74% assembly size). The repetitive landscape was unusually dominated by an expansion of long terminal repeat (LTR) transposable elements (57.73% assembly size) which have shaped the evolution genome size, gene composition through retrotransposition of host genes, and ectopic recombination. Genome skims revealed repeat content is more than 2-3 fold higher in limestone endemic O. idahoensis compared to non-calcareous Oreohelix species. Gene family size analysis revealed stress and biomineralization genes have expanded significantly in the O. idahoensis genome.

CONCLUSIONS

Hundreds of threatened land snail species are endemic to calcareous rock regions but there are very few genomic resources available to guide their conservation or determine the genomic architecture underlying CaCO₃ resource specialization. Our study provides one of the first high quality draft genomes of a calcareous rock endemic land snail which will serve as a foundation for the conservation genomics of this threatened species and for other groups. The high proportion and activity of LTRs in the O. idahoensis genome is unprecedented in molluscan genomics and sheds new light how transposable element content can vary across molluscs. The genomic resources reported here will enable further studies of the genomic mechanisms underlying calcareous rock specialization and the evolution of transposable element content across molluscs.

MicroRNAs are deeply linked to the emergence of the complex octopus brain

Soft-bodied cephalopods such as octopuses are exceptionally intelligent invertebrates with a highly complex nervous system that evolved independently from vertebrates. Because of elevated RNA editing in their nervous tissues, we hypothesized that RNA regulation may play a major role in the cognitive success of this group. We thus profiled messenger RNAs and small RNAs in three cephalopod species including 18 tissues of the Octopus vulgaris. We show that the major RNA innovation of soft-bodied cephalopods is an expansion of the microRNA (miRNA) gene repertoire. These evolutionarily novel miRNAs were primarily expressed in adult neuronal tissues and during the development and had conserved and thus likely functional target sites. The only comparable miRNA expansions happened, notably, in vertebrates. Thus, we propose that miRNAs are intimately linked to the evolution of complex animal brains.

Exploring LSU and ITS rDNA Sequences for Acanthamoeba Identification and Phylogeny

The identification and classification of strains of Acanthamoeba, a potentially pathogenic ubiquitous free-living amoeba, are largely based on the analysis of 18S rDNA sequences, currently delineating 23 genotypes, T1 to T23. In this study, the sequences of the ITS region, i.e., the 5.8S rDNA and the two internal transcribed spacers (ITS-1 and ITS-2), and those of the large subunit (LSU) rDNA of Acanthamoeba were recovered from amoeba genomes; the sequences are available in GenBank. The complete ITS–LSU sequences could be obtained for 15 strains belonging to 7 distinct lineages (T4A, T4D, T4F, T4G, T2, T5, and T18), and the site of the hidden break producing the 26Sα and 26Sβ was identified. For the other lines, either the LSU is partial (T2/T6, T7) or the ITS is fragmentary (T7, T10, T22). It is noteworthy that a number of sequences assigned to fungi turned out to actually be Acanthamoeba, only some of which could be affiliated with known genotypes. Analysis of the obtained sequences indicates that both ITS and LSU are promising for diagnostic and phylogenetic purposes.

Non-Coding RNAs and Oral Cancer: Small Molecules With Big Functions

Oral cancer remains a major public concern with considerable socioeconomic impact in the world. Despite substantial advancements have been made in treating oral cancer, the five-year survival rate for oral cancer remained undesirable, and the molecular mechanisms underlying OSCC carcinogenesis have not been fully understood. Noncoding RNAs (ncRNAs) include transfer RNAs (tRNAs), as well as small RNAs such as microRNAs, and the long ncRNAs such as HOTAIR are a large segment of the transcriptome that do not have apparent protein-coding roles, but they have been verified to play important roles in diverse biological processes, including cancer cell development. Cell death, such as apoptosis, necrosis, and autophagy, plays a vital role in the progression of cancer. A better understanding of the regulatory relationships between ncRNAs and these various types of cancer cell death is therefore urgently required. The occurrence and development of oral cancer can be controlled by increasing or decreasing the expression of ncRNAs, a method which confers broad prospects for oral cancer treatment. Therefore, it is urgent for us to understand the influence of ncRNAs on the development of different modes of oral tumor death, and to evaluate whether ncRNAs have the potential to be used as biological targets for inducing cell death and recurrence of chemotherapy. The purpose of this review is to describe the impact of ncRNAs on cell apoptosis and autophagy in oral cancer in order to explore potential targets for oral cancer therapy.

Using total RNA quality metrics for time since deposition estimates in degrading bloodstains

The physicochemical changes occurring in biomolecules in degrading bloodstains can be used to approximate the time since deposition (TSD) of bloodstains. This would provide forensic scientists with critical information regarding the timeline of the events involving bloodshed. Our study aims to quantify the timewise degradation trends and temperature dependence found in total RNA from bloodstains without the use of amplification, expanding the scope of the RNA TSD research which has traditionally targeted mRNA and miRNA. Bovine blood with ACD-A anticoagulant was deposited and stored in plastic microcentrifuge tubes at 21 or 4°C and tested over different timepoints spanning 1 week. Total RNA was extracted from each sample and analyzed using automated high sensitivity gel electrophoresis. Nine RNA metrics were visually assessed and quantified using linear and mixed models. The RNA Integrity Number equivalent (RINe) and DV200 were not influenced by the addition of anticoagulant and demonstrated strong negative trends over time. The RINe model fit was high (R²  = 0.60), and while including the biological replicate as a random effect increased the fit for all RNA metrics, no significant differences were found between biological replicates stored at the same temperature for the RINe and DV200. This suggests that these standardized metrics can be directly compared between scenarios and individuals, with DV200 having an inflection point at approximately 28 h. This study provides a novel approach for blood TSD research, revealing metrics that are not affected by inter-individual variation, and improving our understanding of the rapid RNA degradation occurring in bloodstains.

Typical structure of rRNA coding genes in diplonemids points to two independent origins of the bizarre rDNA structures of euglenozoans

Background

Members of Euglenozoa (Discoba) are known for unorthodox rDNA organization. In Euglenida rDNA is located on extrachromosomal circular DNA. In Kinetoplastea and Euglenida the core of the large ribosomal subunit, typically formed by the 28S rRNA, consists of several smaller rRNAs. They are the result of the presence of additional internal transcribed spacers (ITSs) in the rDNA. Diplonemea is the third of the main groups of Euglenozoa and its members are known to be among the most abundant and diverse protists in the oceans. Despite that, the rRNA of only one diplonemid species, Diplonema papillatum, has been examined so far and found to exhibit continuous 28S rRNA. Currently, the rDNA organization has not been researched for any diplonemid. Herein we investigate the structure of rRNA genes in classical (Diplonemidae) and deep-sea diplonemids (Eupelagonemidae), representing the majority of known diplonemid diversity. The results fill the gap in knowledge about diplonemid rDNA and allow better understanding of the evolution of the fragmented structure of the rDNA in Euglenozoa.

Results

We used available genomic (culture and single-cell) sequencing data to assemble complete or almost complete rRNA operons for three classical and six deep-sea diplonemids. The rDNA sequences acquired for several euglenids and kinetoplastids were used to provide the background for the analysis. In all nine diplonemids, 28S rRNA seems to be contiguous, with no additional ITSs detected. Similarly, no additional ITSs were detected in basal prokinetoplastids. However, we identified five additional ITSs in the 28S rRNA of all analysed metakinetoplastids, and up to twelve in euglenids. Only three of these share positions, and they cannot be traced back to their common ancestor.

Conclusions

Presented results indicate that independent origin of additional ITSs in euglenids and kinetoplastids seems to be the most likely. The reason for such unmatched fragmentation remains unknown, but for some reason euglenozoan ribosomes appear to be prone to 28S rRNA fragmentation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02014-9.

Evaluation of transcriptional biomarkers using a high-resolution regression approach: Concentration-dependence of selected transcripts in copper-exposed freshwater mussels (Anodonta anatina)

We tested concentration-dependence of selected gene transcripts (cat, gst, hsp70, hsp90, mt and sod) for evaluation as biomarkers of chemical stress. Contrary to the common approach of factorial designs and few exposure concentrations, we used regression across a high-resolution concentration series. Specifically, freshwater mussels (Anodonta anatina) were acutely (96 h) exposed to Cu (13 nominal concentrations, measuring 0.13-1 600 µg/L), and transcripts were measured by RT-qPCR. In digestive glands, cat, hsp90 and mt decreased with water Cu (p < 0.05), but response magnitudes saturated at < 2-fold decreases. In gills, gst, hsp70, hsp90 and mt increased with water Cu (p < 0.05). While hsp70, hsp90 and mt exceeded 2-fold increases within the exposure range, high Cu concentrations were required (38-160 µg/L). Although gill responses were generally more robust compared to digestive glands, overall small response magnitudes and moderate sensitivity may set limit for potential application as general biomarkers of chemical stress.

Bioinformatic Tools for the Analysis and Prediction of ncRNA Interactions

Noncoding RNAs (ncRNAs) play prominent roles in the regulation of gene expression via their interactions with other biological molecules such as proteins and nucleic acids. Although much of our knowledge about how these ncRNAs operate in different biological processes has been obtained from experimental findings, computational biology can also clearly substantially boost this knowledge by suggesting possible novel interactions of these ncRNAs with other molecules. Computational predictions are thus used as an alternative source of new insights through a process of mutual enrichment because the information obtained through experiments continuously feeds through into computational methods. The results of these predictions in turn shed light on possible interactions that are subsequently validated experimentally. This review describes the latest advances in databases, bioinformatic tools, and new in silico strategies that allow the establishment or prediction of biological interactions of ncRNAs, particularly miRNAs and lncRNAs. The ncRNA species described in this work have a special emphasis on those found in humans, but information on ncRNA of other species is also included.

Insights into the Molecular Basis of Genome Stability and Pristine Proteostasis in Naked Mole-Rats

The naked mole-rat (Heterocephalus glaber) is the longest-lived rodent, with a maximal reported lifespan of 37 years. In addition to its long lifespan - which is much greater than predicted based on its small body size (longevity quotient of ~4.2) - naked mole-rats are also remarkably healthy well into old age. This is reflected in a striking resistance to tumorigenesis and minimal declines in cardiovascular, neurological and reproductive function in older animals. Over the past two decades, researchers have been investigating the molecular mechanisms regulating the extended life- and health- span of this animal, and since the sequencing and assembly of the naked mole-rat genome in 2011, progress has been rapid. Here, we summarize findings from published studies exploring the unique molecular biology of the naked mole-rat, with a focus on mechanisms and pathways contributing to genome stability and maintenance of proteostasis during aging. We also present new data from our laboratory relevant to the topic and discuss our findings in the context of the published literature.

ACME dissociation: a versatile cell fixation-dissociation method for single-cell transcriptomics

Single-cell sequencing technologies are revolutionizing biology, but they are limited by the need to dissociate live samples. Here, we present ACME (ACetic-MEthanol), a dissociation approach for single-cell transcriptomics that simultaneously fixes cells. ACME-dissociated cells have high RNA integrity, can be cryopreserved multiple times, and are sortable and permeable. As a proof of principle, we provide single-cell transcriptomic data of different species, using both droplet-based and combinatorial barcoding single-cell methods. ACME uses affordable reagents, can be done in most laboratories and even in the field, and thus will accelerate our knowledge of cell types across the tree of life.

Sticky problems: extraction of nucleic acids from molluscs

Traditional molecular methods and omics-techniques across molluscan taxonomy increasingly inform biology of Mollusca. Recovery of DNA and RNA for such studies is challenged by common biological properties of the highly diverse molluscs. Molluscan biomineralization, adhesive structures and mucus involve polyphenolic proteins and mucopolysaccharides that hinder DNA extraction or copurify to inhibit enzyme-catalysed molecular procedures. DNA extraction methods that employ the detergent hexadecyltrimethylammoniumbromide (CTAB) to remove these contaminants importantly facilitate molecular-level study of molluscs. Molluscan pigments may stain DNA samples and interfere with spectrophotometry, necessitating gel electrophoresis or fluorometry for accurate quantification. RNA can reliably be extracted but the 'hidden break' in 28S rRNA of molluscs (like most protostomes) causes 18S and 28S rRNA fragments to co-migrate electrophoretically. This challenges the standard quality control based on the ratio of 18S and 28S rRNA, developed for deuterostome animals. High-AT content in molluscan rRNA prevents the effective purification of polyadenylated mRNA. Awareness of these matters aids the continuous expansion of molecular malacology, enabling work also with museum specimens and next-generation sequencing, with the latter imposing unprecedented demands on DNA quality. Alternative methods to extract nucleic acids from molluscs are available from literature and, importantly, from communications with others who study the molecular biology of molluscs. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Lack of support for Deuterostomia prompts reinterpretation of the first Bilateria

The bilaterally symmetric animals (Bilateria) are considered to comprise two monophyletic groups, Protostomia (Ecdysozoa and the Lophotrochozoa) and Deuterostomia (Chordata and the Xenambulacraria). Recent molecular phylogenetic studies have not consistently supported deuterostome monophyly. Here, we compare support for Protostomia and Deuterostomia using multiple, independent phylogenomic datasets. As expected, Protostomia is always strongly supported, especially by longer and higher-quality genes. Support for Deuterostomia, however, is always equivocal and barely higher than support for paraphyletic alternatives. Conditions that cause tree reconstruction errors-inadequate models, short internal branches, faster evolving genes, and unequal branch lengths-coincide with support for monophyletic deuterostomes. Simulation experiments show that support for Deuterostomia could be explained by systematic error. The branch between bilaterian and deuterostome common ancestors is, at best, very short, supporting the idea that the bilaterian ancestor may have been deuterostome-like. Our findings have important implications for the understanding of early animal evolution.

Ion channel profiling of the Lymnaea stagnalis ganglia via transcriptome analysis

BACKGROUND

The pond snail Lymnaea stagnalis (L. stagnalis) has been widely used as a model organism in neurobiology, ecotoxicology, and parasitology due to the relative simplicity of its central nervous system (CNS). However, its usefulness is restricted by a limited availability of transcriptome data. While sequence information for the L. stagnalis CNS transcripts has been obtained from EST libraries and a de novo RNA-seq assembly, the quality of these assemblies is limited by a combination of low coverage of EST libraries, the fragmented nature of de novo assemblies, and lack of reference genome.

RESULTS

In this study, taking advantage of the recent availability of a preliminary L. stagnalis genome, we generated an RNA-seq library from the adult L. stagnalis CNS, using a combination of genome-guided and de novo assembly programs to identify 17,832 protein-coding L. stagnalis transcripts. We combined our library with existing resources to produce a transcript set with greater sequence length, completeness, and diversity than previously available ones. Using our assembly and functional domain analysis, we profiled L. stagnalis CNS transcripts encoding ion channels and ionotropic receptors, which are key proteins for CNS function, and compared their sequences to other vertebrate and invertebrate model organisms. Interestingly, L. stagnalis transcripts encoding numerous putative Ca2+ channels showed the most sequence similarity to those of Mus musculus, Danio rerio, Xenopus tropicalis, Drosophila melanogaster, and Caenorhabditis elegans, suggesting that many calcium channel-related signaling pathways may be evolutionarily conserved.

CONCLUSIONS

Our study provides the most thorough characterization to date of the L. stagnalis transcriptome and provides insights into differences between vertebrates and invertebrates in CNS transcript diversity, according to function and protein class. Furthermore, this study provides a complete characterization of the ion channels of Lymnaea stagnalis, opening new avenues for future research on fundamental neurobiological processes in this model system.

Bulinus truncatus transcriptome – a resource to enable molecular studies of snail and schistosome biology

Despite advances in high-throughput sequencing and bioinformatics, molecular investigations of snail intermediate hosts that transmit parasitic trematodes are scant. Here, we report the first transcriptome for Bulinus truncatus – a key intermediate host of Schistosoma haematobium – a blood fluke that causes urogenital schistosomiasis in humans. We assembled this transcriptome from short- and long-read RNA-sequence data. From this transcriptome, we predicted 12,998 proteins, 58% of which had orthologs in Biomphalaria glabrata – an intermediate host of Schistosoma mansoni – a blood fluke that causes hepato-intestinal schistosomiasis. We predicted that select protein groups are involved in signal transduction, cell growth and death, the immune system, environmental adaptation and/or the excretory/secretory system, suggesting roles in immune responses, pathogen defence and/or parasite-host interactions. The transcriptome of Bu. truncatus provides a useful resource to underpin future molecular investigations of this and related snail species, and its interactions with pathogens including S. haematobium. The present resource should enable comparative investigations of other molluscan hosts of socioeconomically important parasites in the future.

•

First transcriptome to represent Bulinus truncatus – a snail intermediate host of Schistosoma haematobium.

•

Select protein groups of Bu. truncatus are inferred to associate with innate immune responses against pathogens.

•

Transcriptome provides a resource for future studies of parasite-host interactions and snail-host resistance to pathogens.

Detection of target-site and metabolic resistance to pyrethroids in the bed bug Cimex lectularius in Berlin, Germany

Knockdown-resistance (kdr) against pyrethroids in bed bugs (Cimex lectularis) is associated with the presence of several point mutations in the voltage-sensitive sodium channel α-subunit gene and/or an increased metabolic detoxification by cytochrome P450 monooxygenases (CYPs). In the present study, pyrosequencing assays were developed to quantify the presence of the kdr substitutions (V419L or L925I substitution) in bed bugs in Berlin, Germany. In 14 of 17 bed bug field strains, pyrosequencing revealed the presence of the substitution L925I with allele frequencies between 30% and 100%. One field strain additionally carried the substitution V419L with allele frequencies of 40% in males and 96% in females. In seven of the 17 field strains, mRNA levels of four CYP genes were examined using RT-qPCR. Relative to a susceptible laboratory reference strain, five field strains showed significantly higher mRNA levels of cyp397a1 with 7.1 to 56-fold increases. One of these strains additionally showed a 4.9-fold higher mRNA level of cyp398a1 compared to the reference strain, while cyp4cm1 and cyp6dn1 showed no significant differences. Our findings indicate that multiple resistance mechanisms are present in German C. lectularius populations simultaneously.

Shedding light: a phylotranscriptomic perspective illuminates the origin of photosymbiosis in marine bivalves

BACKGROUND

Photosymbiotic associations between metazoan hosts and photosynthetic dinoflagellates are crucial to the trophic and structural integrity of many marine ecosystems, including coral reefs. Although extensive efforts have been devoted to study the short-term ecological interactions between coral hosts and their symbionts, long-term evolutionary dynamics of photosymbiosis in many marine animals are not well understood. Within Bivalvia, the second largest class of mollusks, obligate photosymbiosis is found in two marine lineages: the giant clams (subfamily Tridacninae) and the heart cockles (subfamily Fraginae), both in the family Cardiidae. Morphologically, giant clams show relatively conservative shell forms whereas photosymbiotic fragines exhibit a diverse suite of anatomical adaptations including flattened shells, leafy mantle extensions, and lens-like microstructural structures. To date, the phylogenetic relationships between these two subfamilies remain poorly resolved, and it is unclear whether photosymbiosis in cardiids originated once or twice.

RESULTS

In this study, we establish a backbone phylogeny for Cardiidae utilizing RNASeq-based transcriptomic data from Tridacninae, Fraginae and other cardiids. A variety of phylogenomic approaches were used to infer the relationship between the two groups. Our analyses found conflicting gene signals and potential rapid divergence among the lineages. Overall, results support a sister group relationship between Tridacninae and Fraginae, which diverged during the Cretaceous. Although a sister group relationship is recovered, ancestral state reconstruction using maximum likelihood-based methods reveals two independent origins of photosymbiosis, one at the base of Tridacninae and the other within a symbiotic Fraginae clade.

CONCLUSIONS

The newly revealed common ancestry between Tridacninae and Fraginae brings a possibility that certain genetic, metabolic, and/or anatomical exaptations existed in their last common ancestor, which promoted both lineages to independently establish photosymbiosis, possibly in response to the modern expansion of reef habitats.