Combining a transcriptomic approach and a targeted metabolomics approach for deciphering the molecular bases of compatibility phenotype in the snail Biomphalaria glabrata toward Schistosoma mansoni

Biomphalaria glabrata is a freshwater snail and the obligatory intermediate host of Schistosoma mansoni parasite, the etiologic agent of intestinal Schistosomiasis, in South America and Caribbean. Interestingly in such host-parasite interactions, compatibility varies between populations, strains or individuals. This observed compatibility polymorphism is based on a complex molecular-matching-phenotype, the molecular bases of which have been investigated in numerous studies, notably by comparing between different strains or geographical isolates or clonal selected snail lines. Herein we propose to decipher the constitutive molecular support of this interaction in selected non-clonal resistant and susceptible snail strain originating from the same natural population from Brazil and thus having the same genetic background. Thanks to a global RNAseq transcriptomic approach on whole snail, we identified a total of 328 differentially expressed genes between resistant and susceptible phenotypes among which 129 were up-regulated and 199 down-regulated. Metabolomic studies were used to corroborate the RNAseq results. The activation of immune genes and specific metabolic pathways in resistant snails might provide them with the capacity to better respond to parasite infection.

Subtelomeric plasticity contributes to gene family expansion in the human parasitic flatworm Schistosoma mansoni

BACKGROUND

The genomic region that lies between the telomere and chromosome body, termed the subtelomere, is heterochromatic, repeat-rich, and frequently undergoes rearrangement. Within this region, large-scale structural changes enable gene diversification, and, as such, large multicopy gene families are often found at the subtelomere. In some parasites, genes associated with proliferation, invasion, and survival are often found in these regions, where they benefit from the subtelomere's highly plastic, rapidly changing nature. The increasing availability of complete (or near complete) parasite genomes provides an opportunity to investigate these typically poorly defined and overlooked genomic regions and potentially reveal relevant gene families necessary for the parasite's lifestyle.

RESULTS

Using the latest chromosome-scale genome assembly and hallmark repeat richness observed at chromosome termini, we have identified and characterised the subtelomeres of Schistosoma mansoni, a metazoan parasitic flatworm that infects over 250 million people worldwide. Approximately 12% of the S. mansoni genome is classified as subtelomeric, and, in line with other organisms, we find these regions to be gene-poor but rich in transposable elements. We find that S. mansoni subtelomeres have undergone extensive interchromosomal recombination and that these sites disproportionately contribute to the 2.3% of the genome derived from segmental duplications. This recombination has led to the expansion of subtelomeric gene clusters containing 103 genes, including the immunomodulatory annexins and other gene families with unknown roles. The largest of these is a 49-copy plexin domain-containing protein cluster, exclusively expressed in the tegument-the tissue located at the host-parasite physical interface-of intramolluscan life stages.

CONCLUSIONS

We propose that subtelomeric regions act as a genomic playground for trial-and-error of gene duplication and subsequent divergence. Owing to the importance of subtelomeric genes in other parasites, gene families implicated in this subtelomeric expansion within S. mansoni warrant further characterisation for a potential role in parasitism.

Automated ChIPmentation procedure on limited biological material of the human blood fluke Schistosoma mansoni

In living cells, the genetic information stored in the DNA sequence is always associated with chromosomal and extra-chromosomal epigenetic information. Chromatin is formed by the DNA and associated proteins, in particular histones. Covalent histone modifications are important bearers of epigenetic information and as such have been increasingly studied since about the year 2000. One of the principal techniques to gather information about the association between DNA and modified histones is chromatin immunoprecipitation (ChIP), also combined with massive sequencing (ChIP-Seq). Automated ChIPmentation procedure is a convenient alternative to native chromatin immunoprecipitation (N-ChIP). It is now routinely used for ChIP-Seq in many model species, using in general roughly 10 6 cells per experiment. Such high cell numbers are sometimes difficult to produce. Using the human parasite Schistosoma mansoni, whose production requires sacrificing animals and should therefore be kept to a minimum, we show here that automated ChIPmentation is suitable for limited biological material. We define the operational limit as ≥20,000 Schistosoma cells with 30,000-300,000 cells as optimum. We also present a streamlined protocol for the preparation of ChIP input libraries.

Cooperation and cheating orchestrate Vibrio assemblages and polymicrobial synergy in oysters infected with OsHV-1 virus

Polymicrobial infections threaten the health of humans and animals but remain understudied in natural systems. We recently described the Pacific Oyster Mortality Syndrome (POMS), a polymicrobial disease affecting oyster production worldwide. In the French Atlantic coast, the disease involves coinfection with ostreid herpesvirus 1 (OsHV-1) and virulent Vibrio. However, it is unknown whether consistent Vibrio populations are associated with POMS in different regions, how Vibrio contribute to POMS, and how they interact with OsHV-1 during pathogenesis. By connecting field-based approaches in a Mediterranean ecosystem, laboratory infection assays and functional genomics, we uncovered a web of interdependencies that shape the structure and function of the POMS pathobiota. We show that Vibrio harveyi and Vibrio rotiferianus are predominant in OsHV-1-diseased oysters and that OsHV-1 drives the partition of the Vibrio community observed in the field. However only V. harveyi synergizes with OsHV-1 by promoting mutual growth and accelerating oyster death. V. harveyi shows high-virulence potential and dampens oyster cellular defenses through a type 3 secretion system, making oysters a more favorable niche for microbe colonization. In addition, V. harveyi produces a key siderophore called vibrioferrin. This important resource promotes the growth of V. rotiferianus, which cooccurs with V. harveyi in diseased oysters, and behaves as a cheater by benefiting from V. harveyi metabolite sharing. Our data show that cooperative behaviors contribute to synergy between bacterial and viral coinfecting partners. Additional cheating behaviors further shape the polymicrobial consortium. Controlling cooperative behaviors or countering their effects opens avenues for mitigating polymicrobial diseases.

Targeted insertion and reporter transgene activity at a gene safe harbor of the human blood fluke, Schistosoma mansoni

The identification and characterization of genomic safe harbor sites (GSHs) can facilitate consistent transgene activity with minimal disruption to the host cell genome. We combined computational genome annotation and chromatin structure analysis to predict the location of four GSHs in the human blood fluke, Schistosoma mansoni, a major infectious pathogen of the tropics. A transgene was introduced via CRISPR-Cas-assisted homology-directed repair into one of the GSHs in the egg of the parasite. Gene editing efficiencies of 24% and transgene-encoded fluorescence of 75% of gene-edited schistosome eggs were observed. The approach advances functional genomics for schistosomes by providing a tractable path for generating transgenics using homology-directed, repair-catalyzed transgene insertion. We also suggest that this work will serve as a roadmap for the development of similar approaches in helminths more broadly.

Single cell RNA sequencing reveals hemocyte heterogeneity in Biomphalaria glabrata: Plasticity over diversity

The freshwater snail Biomphalaria glabrata is an intermediate host of Schistosoma mansoni, the agent of human intestinal schistosomiasis. However, much is to be discovered about its innate immune system that appears as a complex black box, in which the immune cells (called hemocytes) play a major role in both cellular and humoral response towards pathogens. Until now, hemocyte classification has been based exclusively on cell morphology and ultrastructural description and depending on the authors considered from 2 to 5 hemocyte populations have been described. In this study, we proposed to evaluate the hemocyte heterogeneity at the transcriptomic level. To accomplish this objective, we used single cell RNA sequencing (scRNAseq) technology coupled to a droplet-based system to separate hemocytes and analyze their transcriptome at a unique cell level in naive Biomphalaria glabrata snails. We were able to demonstrate the presence of 7 hemocyte transcriptomic populations defined by the expression of specific marker genes. As a result, scRNAseq approach showed a high heterogeneity within hemocytes, but provides a detailed description of the different hemocyte transcriptomic populations in B. glabrata supported by distinct cellular functions and lineage trajectory. As a main result, scRNAseq revealed the 3 main population as a super-group of hemocyte diversity but, on the contrary, a great hemocytes plasticity with a probable capacity of hemocytes to engage to different activation pathways. This work opens a new field of research to understand the role of hemocytes particularly in response to pathogens, and towards S. mansoni parasites.

Early life microbial exposures shape the Crassostrea gigas immune system for lifelong and intergenerational disease protection

BACKGROUND

The interaction of organisms with their surrounding microbial communities influences many biological processes, a notable example of which is the shaping of the immune system in early life. In the Pacific oyster, Crassostrea gigas, the role of the environmental microbial community on immune system maturation - and, importantly, protection from infectious disease - is still an open question.

RESULTS

Here, we demonstrate that early life microbial exposure durably improves oyster survival when challenged with the pathogen causing Pacific oyster mortality syndrome (POMS), both in the exposed generation and in the subsequent one. Combining microbiota, transcriptomic, genetic, and epigenetic analyses, we show that the microbial exposure induced changes in epigenetic marks and a reprogramming of immune gene expression leading to long-term and intergenerational immune protection against POMS.

CONCLUSIONS

We anticipate that this protection likely extends to additional pathogens and may prove to be an important new strategy for safeguarding oyster aquaculture efforts from infectious disease. tag the videobyte/videoabstract in this section Video Abstract.

Hit-and-Run Epigenetic Editing for Vectors of Snail-Borne Parasitic Diseases

Snail-borne parasitic diseases represent an important challenge to human and animal health. Control strategies that target the intermediate snail host has proved very effective. Epigenetic mechanisms are involved in developmental processes and therefore play a fundamental role in developmental variation. DNA methylation is an important epigenetic information carrier in eukaryotes that plays a major role in the control of chromatin structure. Epigenome editing tools have been instrumental to demonstrate functional importance of this mark for gene expression in vertebrates. In invertebrates, such tools are missing, and the role of DNA methylation remains unknown. Here we demonstrate that methylome engineering can be used to modify in vivo the CpG methylation level of a target gene in the freshwater snail Biomphalaria glabrata, intermediate host of the human parasite Schistosoma mansoni. We used a dCas9-SunTag-DNMT3A complex and synthetic sgRNA to transfect B. glabrata embryos and observed an increase of CpG methylation at the target site in 50% of the hatching snails.

From Methylome to Integrative Analysis of Tissue Specificity

DNA methylation is the most studied epigenetic mark in both plants and animals. The gold standard for assaying genome-wide DNA methylation at single-base resolution is whole-genome bisulfite sequencing (WGBS). Here, we describe an improved procedure for WGBS and original bioinformatic workflows applied to unravel tissue-specific variations of the methylome in relation to gene expression and accumulation of secondary metabolites in the medicinal plant Catharanthus roseus.

The methylome of Biomphalaria glabrata and other mollusks: enduring modification of epigenetic landscape and phenotypic traits by a new DNA methylation inhibitor

BACKGROUND

5-Methylcytosine (5mC) is an important epigenetic mark in eukaryotes. Little information about its role exists for invertebrates. To investigate the contribution of 5mC to phenotypic variation in invertebrates, alteration of methylation patterns needs to be produced. Here, we apply new non-nucleoside DNA methyltransferase inhibitors (DNMTi) to introduce aleatory changes into the methylome of mollusk species.

RESULTS

Flavanone inhibitor Flv1 was efficient in reducing 5mC in the freshwater snails Biomphalaria glabrata and Physa acuta, and to a lesser degree, probably due to lower stability in sea water, in the oyster Crassostrea gigas. Flv1 has no toxic effects and significantly decreased the 5mC level in the treated B. glabrata and in its offspring. Drug treatment triggers significant variation in the shell height in both generations. A reduced representation bisulfite-sequencing method called epiGBS corroborates hypomethylation effect of Flv1 in both B. glabrata generations and identifies seven Differential Methylated Regions (DMR) out of 32 found both in Flv1-exposed snails and its progeny, from which 5 were hypomethylated, demonstrating a multigenerational effect. By targeted bisulfite sequencing, we confirmed hypomethylation in a locus and show that it is associated with reduced gene expression.

CONCLUSIONS

Flv1 is a new and efficient DNMTi that can be used to induce transient and heritable modifications of the epigenetic landscape and phenotypic traits in mollusks, a phylum of the invertebrates in which epigenetics is understudied.

RNAi suppression of DNA methylation affects the drought stress response and genome integrity in transgenic poplar

Trees are long-lived organisms that continuously adapt to their environments, a process in which epigenetic mechanisms are likely to play a key role. Via downregulation of the chromatin remodeler DECREASED IN DNA METHYLATION 1 (DDM1) in poplar (Populus tremula × Populus alba) RNAi lines, we examined how DNA methylation coordinates genomic and physiological responses to moderate water deficit. We compared the growth and drought response of two RNAi-ddm1 lines to wild-type (WT) trees under well-watered and water deficit/rewatering conditions, and analyzed their methylomes, transcriptomes, mobilomes and phytohormone contents in the shoot apical meristem. The RNAi-ddm1 lines were more tolerant to drought-induced cavitation but did not differ in height or stem diameter growth. About 5000 differentially methylated regions were consistently detected in both RNAi-ddm1 lines, colocalizing with 910 genes and 89 active transposable elements. Under water deficit conditions, 136 differentially expressed genes were found, including many involved in phytohormone pathways; changes in phytohormone concentrations were also detected. Finally, the combination of hypomethylation and drought led to the mobility of two transposable elements. Our findings suggest major roles for DNA methylation in regulation of genes involved in hormone-related stress responses, and the maintenance of genome integrity through repression of transposable elements.

Satellite-Like W-Elements: Repetitive, Transcribed, and Putative Mobile Genetic Factors with Potential Roles for Biology and Evolution of Schistosoma mansoni

A large portion of animal and plant genomes consists of noncoding DNA. This part includes tandemly repeated sequences and gained attention because it offers exciting insights into genome biology. We investigated satellite-DNA elements of the platyhelminth Schistosoma mansoni, a parasite with remarkable biological features. Schistosoma mansoni lives in the vasculature of humans causing schistosomiasis, a disease of worldwide importance. Schistosomes are the only trematodes that have evolved separate sexes, and the sexual maturation of the female depends on constant pairing with the male. The schistosome karyotype comprises eight chromosome pairs, males are homogametic (ZZ) and females are heterogametic (ZW). Part of the repetitive DNA of S. mansoni are W-elements (WEs), originally discovered as female-specific satellite DNAs in the heterochromatic block of the W-chromosome. Based on new genome and transcriptome data, we performed a reanalysis of the W-element families (WEFs). Besides a new classification of 19 WEFs, we provide first evidence for stage-, sex-, pairing-, gonad-, and strain-specific/preferential transcription of WEs as well as their mobile nature, deduced from autosomal copies of full-length and partial WEs. Structural analyses suggested roles as sources of noncoding RNA-like hammerhead ribozymes, for which we obtained functional evidence. Finally, the variable WEF occurrence in different schistosome species revealed remarkable divergence. From these results, we propose that WEs potentially exert enduring influence on the biology of S. mansoni. Their variable occurrence in different strains, isolates, and species suggests that schistosome WEs may represent genetic factors taking effect on variability and evolution of the family Schistosomatidae.

No pre-zygotic isolation mechanisms between Schistosoma haematobium and Schistosoma bovis parasites: From mating interactions to differential gene expression

Species usually develop reproductive isolation mechanisms allowing them to avoid interbreeding. These preventive barriers can act before reproduction, "pre-zygotic barriers", or after reproduction, "post-zygotic barriers". Pre-zygotic barriers prevent unfavourable mating, while post-zygotic barriers determine the viability and selective success of the hybrid offspring. Hybridization in parasites and the underlying reproductive isolation mechanisms maintaining their genetic integrity have been overlooked. Using an integrated approach this work aims to quantify the relative importance of pre-zygotic barriers in Schistosoma haematobium x S. bovis crosses. These two co-endemic species cause schistosomiasis, one of the major debilitating parasitic diseases worldwide, and can hybridize naturally. Using mate choice experiments we first tested if a specific mate recognition system exists between both species. Second, using RNA-sequencing we analysed differential gene expression between homo- and hetero-specific pairing in male and female adult parasites. We show that homo- and hetero-specific pairing occurs randomly between these two species, and few genes in both sexes are affected by hetero-specific pairing. This suggests that i) mate choice is not a reproductive isolating factor, and that ii) no pre-zygotic barrier except spatial isolation "by the final vertebrate host" seems to limit interbreeding between these two species. Interestingly, among the few genes affected by the pairing status of the worms, some can be related to pathways affected during male and female interactions and may also present interesting candidates for species isolation mechanisms and hybridization in schistosome parasites.

New Insights Into Biomphalysin Gene Family Diversification in the Vector Snail Biomphalaria glabrata

Aerolysins initially characterized as virulence factors in bacteria are increasingly found in massive genome and transcriptome sequencing data from metazoans. Horizontal gene transfer has been demonstrated as the main way of aerolysin-related toxins acquisition in metazoans. However, only few studies have focused on their potential biological functions in such organisms. Herein, we present an extensive characterization of a multigene family encoding aerolysins - named biomphalysin - in Biomphalaria glabrata snail, the intermediate host of the trematode Schistosoma mansoni. Our results highlight that duplication and domestication of an acquired bacterial toxin gene in the snail genome result in the acquisition of a novel and diversified toxin family. Twenty-three biomphalysin genes were identified. All are expressed and exhibited a tissue-specific expression pattern. An in silico structural analysis was performed to highlight the central role played by two distinct domains i) a large lobe involved in the lytic function of these snail toxins which constrained their evolution and ii) a small lobe which is structurally variable between biomphalysin toxins and that matched to various functional domains involved in moiety recognition of targets cells. A functional approach suggests that the repertoire of biomphalysins that bind to pathogens, depends on the type of pathogen encountered. These results underline a neo-and sub-functionalization of the biomphalysin toxins, which have the potential to increase the range of effectors in the snail’s immune arsenal.

Diverging patterns of introgression from Schistosoma bovis across S. haematobium African lineages

Hybridization is a fascinating evolutionary phenomenon that raises the question of how species maintain their integrity. Inter-species hybridization occurs between certain Schistosoma species that can cause important public health and veterinary issues. In particular hybrids between Schistosoma haematobium and S. bovis associated with humans and animals respectively are frequently identified in Africa. Recent genomic evidence indicates that some S. haematobium populations show signatures of genomic introgression from S. bovis. Here, we conducted a genomic comparative study and investigated the genomic relationships between S. haematobium, S. bovis and their hybrids using 19 isolates originating from a wide geographical range over Africa, including samples initially classified as S. haematobium (n = 11), S. bovis (n = 6) and S. haematobium x S. bovis hybrids (n = 2). Based on a whole genomic sequencing approach, we developed 56,181 SNPs that allowed a clear differentiation of S. bovis isolates from a genomic cluster including all S. haematobium isolates and a natural S. haematobium-bovis hybrid. All the isolates from the S. haematobium cluster except the isolate from Madagascar harbored signatures of genomic introgression from S. bovis. Isolates from Corsica, Mali and Egypt harbored the S. bovis-like Invadolysin gene, an introgressed tract that has been previously detected in some introgressed S. haematobium populations from Niger. Together our results highlight the fact that introgression from S. bovis is widespread across S. haematobium and that the observed introgression is unidirectional.

Hybridization is a fascinating evolutionary phenomenon that raises the question of how species maintain their integrity. Inter-species hybridization occurs between certain Schistosoma species that can cause important public health and veterinary issues. In particular hybrids between Schistosoma haematobium and S. bovis associated with humans and animals respectively are frequently identified in Africa. Recent genomic evidence indicates that some S. haematobium populations show signatures of genomic introgression from S. bovis. Here we conducted a comparative genomic study to assess the genomic diversity within S. haematobium and S. bovis species and genetic differentation at the genome scale between these two sister species over the African continent. We also investigated traces of possible ancient introgression from one species to another. We found that S. haematobium display low genetic diversity compared to S. bovis. We also found that most S. haematobium samples harbor signature of past introgression with S. bovis at some genomic positions. Our results strongly suggest that introgression occurred long time ago and that such introgression is unidirectional from S. bovis within S. haematobium. Such introgresssion event(s) result in diverging patterns of genomic introgression across S. haematobium lineages.

A simple ATAC-seq protocol for population epigenetics

We describe here a protocol for the generation of sequence-ready libraries for population epigenomics studies, and the analysis of alignment results. We show that the protocol can be used to monitor chromatin structure changes in populations when exposed to environmental cues. The protocol is a streamlined version of the Assay for transposase accessible chromatin with high-throughput sequencing (ATAC-seq) that provides a positive display of accessible, presumably euchromatic regions. The protocol is straightforward and can be used with small individuals such as daphnia and schistosome worms, and probably many other biological samples of comparable size (~10,000 cells), and it requires little molecular biology handling expertise.

Developmental Methylome of the Medicinal Plant Catharanthus roseus Unravels the Tissue-Specific Control of the Monoterpene Indole Alkaloid Pathway by DNA Methylation

Catharanthus roseus produces a wide spectrum of monoterpene indole alkaloids (MIAs). MIA biosynthesis requires a tightly coordinated pathway involving more than 30 enzymatic steps that are spatio-temporally and environmentally regulated so that some MIAs specifically accumulate in restricted plant parts. The first regulatory layer involves a complex network of transcription factors from the basic Helix Loop Helix (bHLH) or AP2 families. In the present manuscript, we investigated whether an additional epigenetic layer could control the organ-, developmental- and environmental-specificity of MIA accumulation. We used Whole-Genome Bisulfite Sequencing (WGBS) together with RNA-seq to identify differentially methylated and expressed genes among nine samples reflecting different plant organs and experimental conditions. Tissue specific gene expression was associated with specific methylation signatures depending on cytosine contexts and gene parts. Some genes encoding key enzymatic steps from the MIA pathway were found to be simultaneously differentially expressed and methylated in agreement with the corresponding MIA accumulation. In addition, we found that transcription factors were strikingly concerned by DNA methylation variations. Altogether, our integrative analysis supports an epigenetic regulation of specialized metabolisms in plants and more likely targeting transcription factors which in turn may control the expression of enzyme-encoding genes.

Contribution of Viral Genomic Diversity to Oyster Susceptibility in the Pacific Oyster Mortality Syndrome

Juvenile Pacific oysters (Crassostrea gigas) are subjected to recurrent episodes of mass mortalities that constitute a threat for the oyster industry. This mortality syndrome named “Pacific Oyster Mortality Syndrome” (POMS) is a polymicrobial disease whose pathogenesis is initiated by a primary infection by a variant of an Ostreid herpes virus named OsHV-1 μVar. The characterization of the OsHV-1 genome during different disease outbreaks occurring in different geographic areas has revealed the existence of a genomic diversity for OsHV-1 μVar. However, the biological significance of this diversity is still poorly understood. To go further in understanding the consequences of OsHV-1 diversity on POMS, we challenged five biparental families of oysters to two different infectious environments on the French coasts (Atlantic and Mediterranean). We observed that the susceptibility to POMS can be different among families within the same environment but also for the same family between the two environments. Viral diversity analysis revealed that Atlantic and Mediterranean POMS are caused by two distinct viral populations. Moreover, we observed that different oyster families are infected by distinct viral populations within a same infectious environment. Altogether these results suggest that the co-evolutionary processes at play between OsHV-1 μVar and oyster populations have selected a viral diversity that could facilitate the infection process and the transmission in oyster populations. These new data must be taken into account in the development of novel selective breeding programs better adapted to the oyster culture environment.

A simple ATAC-seq protocol for population epigenetics

We describe here a protocol for the generation of sequence-ready libraries for population epigenomics studies. The protocol is a streamlined version of the Assay for transposase accessible chromatin with high-throughput sequencing (ATAC-seq) that provides a positive display of accessible, presumably euchromatic regions. The protocol is straightforward and can be used with small individuals such as daphnia and schistosome worms, and probably many other biological samples of comparable size, and it requires little molecular biology handling expertise.

Correction to: Measuring Histone Modifications in the Human Parasite Schistosoma mansoni

Correction to: Chapter 9 in: David J. Timson (ed.), Schistosoma mansoni: Methods and Protocols, Methods in Molecular Biology, vol. 2151.

Natural resistance to Fasciola hepatica (Trematoda) in Pseudosuccinea columella snails: A review from literature and insights from comparative "omic" analyses

The snail Pseudosuccinea columella is one of the main vectors of the medically-important trematode Fasciola hepatica. In Cuba, the existence of natural P. columella populations that are either susceptible or resistant to F. hepatica infection offers a unique snail-parasite for study of parasite-host compatibility and immune function in gastropods. Here, we review all previous literature on this system and present new "omic" data that provide a molecular baseline of both P. columella phenotypes from naïve snails. Comparison of whole snail transcriptomes (RNAseq) and the proteomes of the albumen gland (2D-electrophoresis, MS) revealed that resistant and susceptible strains differed mainly in an enrichment of particular biological processes/functions and a greater abundance of proteins/transcripts associated with immune defense/stress response in resistant snails. These results indicate a differential allocation of molecular resources to self-maintenance and survival in resistant P. columella that may cause enhanced responsiveness to stressors (i.e. F. hepatica infection or tolerance to variations in environmental pH/total water hardness), possibly as trade-off against reproduction and the ecological cost of resistance previously suggested in resistant populations of P. columella.

On the Importance to Acknowledge Transposable Elements in Epigenomic Analyses

Eukaryotic genomes comprise a large proportion of repeated sequences, an important fraction of which are transposable elements (TEs). TEs are mobile elements that have a significant impact on genome evolution and on gene functioning. Although some TE insertions could provide adaptive advantages to species, transposition is a highly mutagenic event that has to be tightly controlled to ensure its viability. Genomes have evolved sophisticated mechanisms to control TE activity, the most important being epigenetic silencing. However, the epigenetic control of TEs can also affect genes located nearby that can become epigenetically regulated. It has been proposed that the combination of TE mobilization and the induced changes in the epigenetic landscape could allow a rapid phenotypic adaptation to global environmental changes. In this review, we argue the crucial need to take into account the repeated part of genomes when studying the global impact of epigenetic modifications on an organism. We emphasize more particularly why it is important to carefully consider TEs and what bioinformatic tools can be used to do so.

Sympatric versus allopatric evolutionary contexts shape differential immune response in Biomphalaria / Schistosoma interaction

Selective pressures between hosts and their parasites can result in reciprocal evolution or adaptation of specific life history traits. Local adaptation of resident hosts and parasites should lead to increase parasite infectivity/virulence (higher compatibility) when infecting hosts from the same location (in sympatry) than from a foreign location (in allopatry). Analysis of geographic variations in compatibility phenotypes is the most common proxy used to infer local adaptation. However, in some cases, allopatric host-parasite systems demonstrate similar or greater compatibility than in sympatry. In such cases, the potential for local adaptation remains unclear. Here, we study the interaction between Schistosoma and its vector snail Biomphalaria in which such discrepancy in local versus foreign compatibility phenotype has been reported. Herein, we aim at bridging this gap of knowledge by comparing life history traits (immune cellular response, host mortality, and parasite growth) and molecular responses in highly compatible sympatric and allopatric Schistosoma/Biomphalaria interactions originating from different geographic localities (Brazil, Venezuela and Burundi). We found that despite displaying similar prevalence phenotypes, sympatric schistosomes triggered a rapid immune suppression (dual-RNAseq analyses) in the snails within 24h post infection, whereas infection by allopatric schistosomes (regardless of the species) was associated with immune cell proliferation and triggered a non-specific generalized immune response after 96h. We observed that, sympatric schistosomes grow more rapidly. Finally, we identify miRNAs differentially expressed by Schistosoma mansoni that target host immune genes and could be responsible for hijacking the host immune response during the sympatric interaction. We show that despite having similar prevalence phenotypes, sympatric and allopatric snail-Schistosoma interactions displayed strong differences in their immunobiological molecular dialogue. Understanding the mechanisms allowing parasites to adapt rapidly and efficiently to new hosts is critical to control disease emergence and risks of Schistosomiasis outbreaks.

Schistosomiasis, the second most widespread human parasitic disease after malaria, is caused by helminth parasites of the genus Schistosoma. More than 200 million people in 74 countries suffer from the pathological, and societal consequences of this disease. To complete its life cycle, the parasite requires an intermediate host, a freshwater snail of the genus Biomphalaria for its transmission. Given the limited options for treating Schistosoma mansoni infections in humans, much research has focused on developing methods to control transmission by its intermediate snail host. Biomphalaria glabrata. Comparative studies have shown that infection of the snail triggers complex cellular and humoral immune responses resulting in significant variations in parasite infectivity and snail susceptibility, known as the so-called polymorphism of compatibility. However, studies have mostly focused on characterizing the immunobiological mechanisms in sympatric interactions. Herein we used a combination of molecular and phenotypic approaches to compare the effect of infection in various sympatric and allopatric evolutionary contexts, allowing us to better understand the mechanisms of host-parasite local adaptation. Learning more about the immunobiological interactions between B. glabrata and S. mansoni could have important socioeconomic and public health impacts by changing the way we attempt to eradicate parasitic diseases and prevent or control schistosomiasis in the field.

Immune-suppression by OsHV-1 viral infection causes fatal bacteraemia in Pacific oysters

Infectious diseases are mostly explored using reductionist approaches despite repeated evidence showing them to be strongly influenced by numerous interacting host and environmental factors. Many diseases with a complex aetiology therefore remain misunderstood. By developing a holistic approach to tackle the complexity of interactions, we decipher the complex intra-host interactions underlying Pacific oyster mortality syndrome affecting juveniles of Crassostrea gigas, the main oyster species exploited worldwide. Using experimental infections reproducing the natural route of infection and combining thorough molecular analyses of oyster families with contrasted susceptibilities, we demonstrate that the disease is caused by multiple infection with an initial and necessary step of infection of oyster haemocytes by the Ostreid herpesvirus OsHV-1 µVar. Viral replication leads to the host entering an immune-compromised state, evolving towards subsequent bacteraemia by opportunistic bacteria. We propose the application of our integrative approach to decipher other multifactorial diseases that affect non-model species worldwide.

Pacific oyster mortality syndrome is a poorly understood cause of mortality in commercially important oyster species. Here, the authors use multiple infection experiments to show that the syndrome is caused by sequential infection by herpesvirus and opportunistic bacteria.

The anti-fecundity effect of 5-azacytidine (5-AzaC) on Schistosoma mansoni is linked to dis-regulated transcription, translation and stem cell activities

Uncontrolled host immunological reactions directed against tissue-trapped eggs precipitate a potentially lethal, pathological cascade responsible for schistosomiasis. Blocking schistosome egg production, therefore, presents a strategy for simultaneously reducing immunopathology as well as limiting disease transmission in endemic or emerging areas. We recently demonstrated that the ribonucleoside analogue 5-azacytidine (5-AzaC) inhibited Schistosoma mansoni oviposition, egg maturation and ovarian development. While these anti-fecundity effects were associated with a loss of DNA methylation, other molecular processes affected by 5-AzaC were not examined at the time. By comparing the transcriptomes of 5-AzaC-treated females to controls, we provide evidence that this ribonucleoside analogue also modulates other crucial aspects of schistosome egg-laying biology. For example, S. mansoni gene products associated with amino acid-, carbohydrate-, fatty acid-, nucleotide- and tricarboxylic acid (TCA)- homeostasis are all dysregulated in 5-AzaC treated females. To validate the metabolic pathway most significantly affected by 5-AzaC, amino acid metabolism, nascent protein synthesis was subsequently quantified in adult schistosomes. Here, 5-AzaC inhibited this process by 68% ±16.7% (SEM) in male- and 81% ±4.8% (SEM) in female-schistosomes. Furthermore, the transcriptome data indicated that adult female stem cells were also affected by 5-AzaC. For instance, 40% of transcripts associated with proliferating schistosome cells were significantly down-regulated by 5-AzaC. This finding correlated with a considerable reduction (95%) in the number of 5-ethynyl-2'-deoxyuridine (EdU) positive cells found in 5-AzaC-treated females. In addition to protein coding genes, the effect that 5-AzaC had on repetitive element expression was also assessed. Here, 46 repeats were found differentially transcribed between 5-AzaC-treated and control females with long terminal repeat (LTR) and DNA transposon classes being amongst the most significant. This study demonstrates that the anti-fecundity activity of 5-AzaC affects more than just DNA methylation in schistosome parasites. Further characterisation of these processes may reveal novel targets for schistosomiasis control.

Notos - a galaxy tool to analyze CpN observed expected ratios for inferring DNA methylation types

BACKGROUND

DNA methylation patterns store epigenetic information in the vast majority of eukaryotic species. The relatively high costs and technical challenges associated with the detection of DNA methylation however have created a bias in the number of methylation studies towards model organisms. Consequently, it remains challenging to infer kingdom-wide general rules about the functions and evolutionary conservation of DNA methylation. Methylated cytosine is often found in specific CpN dinucleotides, and the frequency distributions of, for instance, CpG observed/expected (CpG o/e) ratios have been used to infer DNA methylation types based on higher mutability of methylated CpG.

RESULTS

Predominantly model-based approaches essentially founded on mixtures of Gaussian distributions are currently used to investigate questions related to the number and position of modes of CpG o/e ratios. These approaches require the selection of an appropriate criterion for determining the best model and will fail if empirical distributions are complex or even merely moderately skewed. We use a kernel density estimation (KDE) based technique for robust and precise characterization of complex CpN o/e distributions without a priori assumptions about the underlying distributions.

CONCLUSIONS

We show that KDE delivers robust descriptions of CpN o/e distributions. For straightforward processing, we have developed a Galaxy tool, called Notos and available at the ToolShed, that calculates these ratios of input FASTA files and fits a density to their empirical distribution. Based on the estimated density the number and shape of modes of the distribution is determined, providing a rational for the prediction of the number and the types of different methylation classes. Notos is written in R and Perl.

Successful Lung Transplantation From Hepatitis C Positive Donor to Seronegative Recipient

Lung transplantation using RNA+ hepatitis C (HCV+) donors to seronegative recipients is not currently performed due to the very high risk of transmission. Previous reports have shown poor survival when this practice was applied. The emergence of new direct-acting antiviral drugs (DAA) suggests a high chance of sustained virologic response in immunocompetent patients. We report here successful transplantation of lungs from HCV+ donor to HCV- recipient. The recipient was an HCV- patient with chronic lung allograft dysfunction. Viral transmission occurred early posttransplant but excellent clinical outcomes were observed including elimination of HCV after 12 weeks of treatment using DAAs.

A multistrain approach to studying the mechanisms underlying compatibility in the interaction between Biomphalaria glabrata and Schistosoma mansoni

In recent decades, numerous studies have sought to better understand the mechanisms underlying the compatibility between Biomphalaria glabrata and Schistosoma mansoni. The developments of comparative transcriptomics, comparative genomics, interactomics and more targeted approaches have enabled researchers to identify a series of candidate genes. However, no molecular comparative work has yet been performed on multiple populations displaying different levels of compatibility. Here, we seek to fill this gap in the literature. We focused on B. glabrata FREPs and S. mansoni SmPoMucs, which were previously demonstrated to be involved in snail/schistosome compatibility. We studied the expression and polymorphisms of these factors in combinations of snail and schistosome isolates that display different levels of compatibility. We found that the polymorphism and expression levels of FREPs and SmPoMucs could be linked to the compatibility level of S. mansoni. These data and our complementary results obtained by RNA-seq of samples from various snail strains indicate that the mechanism of compatibility is much more complex than previously thought, and that it is likely to be highly variable within and between populations. This complexity must be taken into account if we hope to identify the molecular pathways that are most likely to be good targets for strategies aimed at blocking transmission of the parasite through the snail intermediate host.

Schistosomiasis is the second most widespread human tropical parasitic disease after malaria. It is caused by flatworms of the genus Schistosoma, and poses a considerable threat for human health in numerous Asian, African and South American countries. The World Health Organization has set the goal of eradicating schistosomiasis by 2025. However, no vaccine is available, and we currently have only one drug (praziquantel) that can effectively and efficiently treat the disease. As treatment by mass drug administration would enhance the risk of drug resistance in schistosome parasites, complementary strategies to fight this parasitic disease are urgently needed. Freshwater snails of the Biomphalaria genus act as intermediate hosts in the transmission of the schistosome species. Thus, learning more about the mechanisms of the interaction between these snails and the schistosomes could critically facilitate the identification of potential new candidate molecules that may be targeted to prevent schistosome transmission in the field.

Effects of a parental exposure to diuron on Pacific oyster spat methylome

Environmental epigenetic is an emerging field that studies the cause-effect relationship between environmental factors and heritable trait via an alteration in epigenetic marks. This field has received much attentions since the impact of environmental factors on different epigenetic marks have been shown to be associated with a broad range of phenotypic disorders in natural ecosystems. Chemical pollutants have been shown to affect immediate epigenetic information carriers of several aquatic species but the heritability of the chromatin marks and the consequences for long term adaptation remain open questions. In this work, we investigated the impact of the diuron herbicide on the DNA methylation pattern of spat from exposed Crassotrea gigas genitors. This oyster is one of the most important mollusk species produced worldwide and a key coastal economic resource in France. The whole genome bisulfite sequencing (WGBS, BS-Seq) was applied to obtain a methylome at single nucleotide resolution on DNA extracted from spat issued from diuron exposed genitors comparatively to control spat. We showed that the parental diuron exposure has an impact on the DNA methylation pattern of its progeny. Most of the differentially methylated regions occurred within coding sequences and we showed that this change in methylation level correlates with RNA level only in a very small group of genes. Although the DNA methylation profile is variable between individuals, we showed conserved DNA methylation patterns in response to parental diuron exposure. This relevant result opens perspectives for the setting of new markers based on epimutations as early indicators of marine pollutions.

Survival in sensitized lung transplant recipients with perioperative desensitization

Donor-specific HLA antibodies (DSA) have an adverse effect on short-term and long-term lung transplant outcomes. We implemented a perioperative strategy to treat DSA-positive recipients, leading to equivalent rejection and graft survival outcomes. Pretransplant DSA were identified to HLA-A, B, C, DR and DQ antigens. DSA-positive patients were transplanted if panel reactive antibody (PRA) ≥30% or medically urgent and desensitized with perioperative plasma exchange, intravenous immune globulin, antithymocyte globulin (ATG), and mycophenolic acid (MPA). PRA-positive/DSA-negative recipients received MPA. Unsensitized patients received routine cyclosporine, azathioprine and prednisone without ATG. From 2008-2011, 340 lung-only first transplants were performed: 53 DSA-positive, 93 PRA-positive/DSA-negative and 194 unsensitized. Thirty-day survival was 96 %/99%/96% in the three groups, respectively. One-year graft survival was 89%/88%/86% (p = 0.47). DSA-positive and PRA-positive/DSA-negative patients were less likely to experience any ≥ grade 2 acute rejection (9% and 9% vs. 18% unsensitized p = 0.04). Maximum predicted forced expiratory volume (1 s) (81%/74%/76%, p = NS) and predicted forced vital capacity (81%/77%/78%, respectively, p = NS) were equivalent between groups. With the application of this perioperative treatment protocol, lung transplantation can be safely performed in DSA/PRA-positive patients, with similar outcomes to unsensitized recipients.

Rate of cyp51A mutation in Aspergillus fumigatus among lung transplant recipients with targeted prophylaxis

OBJECTIVES

The most common mechanism of azole (itraconazole and voriconazole) resistance in Aspergillus fumigatus is a mutation at the cyp51A locus. The aim of our study was to determine the rate of cyp51A mutations in lung transplant recipients (LTR) undergoing targeted antifungal prophylaxis with 12 weeks of voriconazole.

METHODS

We conducted a prospective study that included 22 LTR with A. fumigatus between October 2008 and November 2011. Of those, 10 LTR were colonized with A. fumigatus and 12 had invasive pulmonary aspergillosis.

RESULTS

Four patients were found to have A. fumigatus isolates with a cyp51A mutation, two had colonization and two had invasive pulmonary aspergillosis. The remaining 18 LTR had WT cyp51A A. fumigatus isolates. All A. fumigatus isolates (except one due to mixed growth) were tested for antifungal susceptibility. A total of nine LTR were exposed to azoles prior to A. fumigatus isolation for a median duration of 249 (IQR 99-524) days. Azole exposure preceded the isolation of two mutant isolates and seven WT isolates. None of the cyp51A mutant isolates conferred phenotypic resistance to azoles.

CONCLUSIONS

Targeted antifungal prophylaxis in LTR did not lead to cyp51A resistance mutations in this cohort. Data on larger cohorts who receive universal antifungal prophylaxis are needed.