FISH on chromosomes derived from the snail model organism Biomphalaria glabrata

PIWI silencing mechanism involving the retrotransposon nimbus orchestrates resistance to infection with Schistosoma mansoni in the snail vector, Biomphalaria glabrata

Background

Schistosomiasis remains widespread in many regions despite efforts at its elimination. By examining changes in the transcriptome at the host-pathogen interface in the snail Biomphalaria glabrata and the blood fluke Schistosoma mansoni, we previously demonstrated that an early stress response in juvenile snails, manifested by induction of heat shock protein 70 (Hsp 70) and Hsp 90 and of the reverse transcriptase (RT) domain of the B. glabrata non-LTR- retrotransposon, nimbus, were critical for B. glabrata susceptibility to S. mansoni. Subsequently, juvenile B. glabrata BS-90 snails, resistant to S. mansoni at 25°C become susceptible by the F2 generation when maintained at 32°C, indicating an epigenetic response.

Methodology/Principal findings

To better understand this plasticity in susceptibility of the BS-90 snail, mRNA sequences were examined from S. mansoni exposed juvenile BS-90 snails cultured either at 25°C (non-permissive temperature) or 32°C (permissive). Comparative analysis of transcriptomes from snails cultured at the non-permissive and permissive temperatures revealed that whereas stress related transcripts dominated the transcriptome of susceptible BS-90 juvenile snails at 32°C, transcripts encoding proteins with a role in epigenetics, such as PIWI (BgPiwi), chromobox protein homolog 1 (BgCBx1), histone acetyltransferase (BgHAT), histone deacetylase (BgHDAC) and metallotransferase (BgMT) were highly expressed in those cultured at 25°C. To identify robust candidate transcripts that will underscore the anti-schistosome phenotype in B. glabrata, further validation of the differential expression of the above transcripts was performed by using the resistant BS-90 (25°C) and the BBO2 susceptible snail stock whose genome has now been sequenced and represents an invaluable resource for molecular studies in B. glabrata. A role for BgPiwi in B. glabrata susceptibility to S. mansoni, was further examined by using siRNA corresponding to the BgPiwi encoding transcript to suppress expression of BgPiwi, rendering the resistant BS-90 juvenile snail susceptible to infection at 25°C. Given transposon silencing activity of PIWI as a facet of its role as guardian of the integrity of the genome, we examined the expression of the nimbus RT encoding transcript at 120 min after infection of resistant BS90 piwi-siRNA treated snails. We observed that nimbus RT was upregulated, indicating that modulation of the transcription of the nimbus RT was associated with susceptibility to S. mansoni in BgPiwi-siRNA treated BS-90 snails. Furthermore, treatment of susceptible BBO2 snails with the RT inhibitor lamivudine, before exposure to S. mansoni, blocked S. mansoni infection concurrent with downregulation of the nimbus RT transcript and upregulation of the BgPiwi encoding transcript in the lamivudine-treated, schistosome-exposed susceptible snails.

Conclusions and significance

These findings support a role for the interplay of BgPiwi and nimbus in the epigenetic modulation of plasticity of resistance/susceptibility in the snail-schistosome relationship.

Progress is being made to eliminate schistosomiasis, a tropical disease that remains endemic in the tropics and neotropics. In 2020, WHO proposed controlling the snail population as part of a strategy toward reducing schistosomiasis, a vector borne disease, by 2025. The life cycle of the causative parasite is, however, complex and in the absence of vaccines, new drugs, and access to clean water and sanitation, reduction of schistosomiasis will remain elusive. To break the parasite’s life cycle during the snail stage of its development, a better understanding of the molecular basis of how schistosomes survive, or not, in the snail is required. By examining changes in the transcriptome at the host-pathogen interface in the snail Biomphalaria glabrata and Schistosoma mansoni, we showed that early stress response, manifested by the induction of Heat Shock Proteins (Hsps) and the RT domain of the non-LTR retrotransposon, nimbus, were critical for snail susceptibility. Subsequently, juvenile B. glabrata BS-90 snails, resistant to S. mansoni at 25°C were observed to become susceptible by the F2 generation when maintained at 32°C, indicating an epigenetic response. This study confirms these earlier results and shows an interplay between PIWI and nimbus in the anti-schistosome response in the snail host.

Differential spatial repositioning of activated genes in Biomphalaria glabrata snails infected with Schistosoma mansoni

Schistosomiasis is an infectious disease infecting mammals as the definitive host and fresh water snails as the intermediate host. Understanding the molecular and biochemical relationship between the causative schistosome parasite and its hosts will be key to understanding and ultimately treating and/or eradicating the disease. There is increasing evidence that pathogens that have co-evolved with their hosts can manipulate their hosts' behaviour at various levels to augment an infection. Bacteria, for example, can induce beneficial chromatin remodelling of the host genome. We have previously shown in vitro that Biomphalaria glabrata embryonic cells co-cultured with schistosome miracidia display genes changing their nuclear location and becoming up-regulated. This also happens in vivo in live intact snails, where early exposure to miracidia also elicits non-random repositioning of genes. We reveal differences in the nuclear repositioning between the response of parasite susceptible snails as compared to resistant snails and with normal or live, attenuated parasites. Interestingly, the stress response gene heat shock protein (Hsp) 70 is only repositioned and then up-regulated in susceptible snails with the normal parasite. This movement and change in gene expression seems to be controlled by the parasite. Other differences in the behaviour of genes support the view that some genes are responding to tissue damage, for example the ferritin genes move and are up-regulated whether the snails are either susceptible or resistant and upon exposure to either normal or attenuated parasite. This is the first time host genome reorganisation has been seen in a parasitic host and only the second time for any pathogen. We believe that the parasite elicits a spatio-epigenetic reorganisation of the host genome to induce favourable gene expression for itself and this might represent a fundamental mechanism present in the human host infected with schistosome cercariae as well as in other host-pathogen relationships.

Bilharzia is a parasitic disease endemic in many parts of the world. The schistosoma parasite that causes Bilharzia infects humans but uses a fresh water snail as a secondary host. These two organisms have co-evolved together, and as such the parasite will have mechanisms to overcome the host defences. Understanding this delicately balanced relationship is fundamental to controlling or eradicating the disease. We have studied how this parasite can influence how the DNA within the snail behaves. We have shown snail genes have specific locations within the cell nuclei. Further, we have revealed that specific snail genes related to a schistosome infection change to a new non-random nuclear location as they are turned on or up-regulated. We have snail strains that are susceptible or resistant to the infection of parasites and we can also take live parasites and make them unable to complete an infection by irradiating them. In this unique study, we have shown a gene that is involved in stress pathways moves to a new nuclear location and becomes turned on, but only in susceptible snails, infected with fully functional parasite. Our data suggest that this gene is regulated by the parasite, which has control over the host's DNA, so that the gene is moved to an area where it can be actively expressed. We have uncovered a novel mechanism whereby the spatial organization of a host organism is interfered with by a pathogen. This type of control is probably found in other host-pathogen relationships.

Fine mapping of the pond snail left-right asymmetry (chirality) locus using RAD-Seq and fibre-FISH

The left-right asymmetry of snails, including the direction of shell coiling, is determined by the delayed effect of a maternal gene on the chiral twist that takes place during early embryonic cell divisions. Yet, despite being a well-established classical problem, the identity of the gene and the means by which left-right asymmetry is established in snails remain unknown. We here demonstrate the power of new genomic approaches for identification of the chirality gene, "D". First, heterozygous (Dd) pond snails Lymnaea stagnalis were self-fertilised or backcrossed, and the genotype of more than six thousand offspring inferred, either dextral (DD/Dd) or sinistral (dd). Then, twenty of the offspring were used for Restriction-site-Associated DNA Sequencing (RAD-Seq) to identify anonymous molecular markers that are linked to the chirality locus. A local genetic map was constructed by genotyping three flanking markers in over three thousand snails. The three markers lie either side of the chirality locus, with one very tightly linked (<0.1 cM). Finally, bacterial artificial chromosomes (BACs) were isolated that contained the three loci. Fluorescent in situ hybridization (FISH) of pachytene cells showed that the three BACs tightly cluster on the same bivalent chromosome. Fibre-FISH identified a region of greater that ∼0.4 Mb between two BAC clone markers that must contain D. This work therefore establishes the resources for molecular identification of the chirality gene and the variation that underpins sinistral and dextral coiling. More generally, the results also show that combining genomic technologies, such as RAD-Seq and high resolution FISH, is a robust approach for mapping key loci in non-model systems.

Identification and characterisation of functional expressed sequence tags-derived simple sequence repeat (eSSR) markers for genetic linkage mapping of Schistosoma mansoni juvenile resistance and susceptibility loci in Biomphalaria glabrata

Biomphalaria glabrata susceptibility to Schistosoma mansoni has a strong genetic component, offering the possibility for investigating host-parasite interactions at the molecular level, perhaps leading to novel control approaches. The identification, mapping and molecular characterisation of genes that influence the outcome of parasitic infection in the intermediate snail host is, therefore, seen as fundamental to the control of schistosomiasis. To better understand the evolutionary processes driving disease resistance/susceptibility phenotypes, we previously identified polymorphic random amplification of polymorphic DNA and genomic simple sequence repeats from B. glabrata. In the present study we identified and characterised polymorphic expressed simple sequence repeats markers (Bg-eSSR) from existing B. glabrata expressed sequence tags. Using these markers, and with previously identified genomic simple sequence repeats, genetic linkage mapping for parasite refractory and susceptibility phenotypes, the first known for B. glabrata, was initiated. Data mining of 54,309 expressed sequence tag, produced 660 expressed simple sequence repeats of which dinucleotide motifs (TA)n were the most common (37.88%), followed by trinucleotide (29.55%), mononucleotide (18.64%) and tetranucleotide (10.15%). Penta- and hexanucleotide motifs represented <3% of the Bg-eSSRs identified. While the majority (71%) of Bg-eSSRs were monomorphic between resistant and susceptible snails, several were, however, useful for the construction of a genetic linkage map based on their inheritance in segregating F2 progeny snails derived from crossing juvenile BS-90 and NMRI snails. Polymorphic Bg-eSSRs assorted into six linkage groups at a logarithm of odds score of 3. Interestingly, the heritability of four markers (Prim1_910, Prim1_771, Prim6_1024 and Prim7_823) with juvenile snail resistance were, by t-test, significant (P<0.05) while an allelic marker, Prim24_524, showed linkage with the juvenile snail susceptibility phenotype. On the basis of our results it is possible that the gene(s) controlling juvenile resistance and susceptibility to S. mansoni infection in B. glabrata are not only on the same linkage group but lie within a short distance (42cM) of each other.