Identification of 18 new transcribed retrotransposons in Schistosoma mansoni

A beginner's guide to manual curation of transposable elements

Background

In the study of transposable elements (TEs), the generation of a high confidence set of consensus sequences that represent the diversity of TEs found in a given genome is a key step in the path to investigate these fascinating genomic elements. Many algorithms and pipelines are available to automatically identify putative TE families present in a genome. Despite the availability of these valuable resources, producing a library of high-quality full-length TE consensus sequences largely remains a process of manual curation. This know-how is often passed on from mentor-to-mentee within research groups, making it difficult for those outside the field to access this highly specialised skill.

Results

Our manuscript attempts to fill this gap by providing a set of detailed computer protocols, software recommendations and video tutorials for those aiming to manually curate TEs. Detailed step-by-step protocols, aimed at the complete beginner, are presented in the Supplementary Methods.

Conclusions

The proposed set of programs and tools presented here will make the process of manual curation achievable and amenable to all researchers and in special to those new to the field of TEs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-021-00259-7.

Species-Specific Proteins in the Oviducts of Snail Sibling Species: Proteotranscriptomic Study of Littorina fabalis and L. obtusata

Genus Littorina subgenus Neritrema (Mollusca, Caenogastropoda) includes the "obtusata" group of closely related species (Littorina obtusata and L. fabalis). The anatomy of the adult reproductive system (pallial oviduct) is the only reliable feature used for species identification in females of these species. Reproductive system anatomy and reproduction-associated proteins often diverge between sibling species. Despite being of high evolutionary interest, the molecular basis of this divergence remains poorly understood. We performed proteotranscriptomic comparison of oviducts of L. obtusata and L. fabalis by RNA-seq on Illumina HiSeq 2500 and two-dimensional protein electrophoresis (2D DIGE) with MS/MS identification of the species-specific proteins. The interspecies differences in the oviduct were associated with (1) metabolic proteins reflecting overall physiological differences between L. obtusata and L. fabalis, (2) receptor proteins, and (3) transcripts related to transposable elements (TEs). Various receptors identified may recognize a wide variety of ligands from pathogen-associated molecular patterns to specific carbohydrates on the sperm surface. Therefore, these may participate in immune defense as well as in sperm storage and regulation. Species-specificity of multiple TE sequences (coding for reverse transcriptase and ribonuclease H) may indicate the important role of these genomic elements in the Littorina species divergence, which has not been reported previously.

The function of twister ribozyme variants in non-LTR retrotransposition in Schistosoma mansoni

The twister ribozyme is widely distributed over numerous organisms and is especially abundant in Schistosoma mansoni, but has no confirmed biological function. Of the 17 non-LTR retrotransposons known in S. mansoni, none have thus far been associated with ribozymes. Here we report the identification of novel twister variant (T-variant) ribozymes and their function in S. mansoni non-LTR retrotransposition. We show that T-variant ribozymes are located at the 5′ end of Perere-3 non-LTR retrotransposons in the S. mansoni genome. T-variant ribozymes were demonstrated to be catalytically active in vitro. In reporter constructs, T-variants were shown to cleave in vivo, and cleavage of T-variants was sufficient for the translation of downstream reporter genes. Our analysis shows that the T-variants and Perere-3 are transcribed together. Target site duplications (TSDs); markers of target-primed reverse transcription (TPRT) and footmarks of retrotransposition, are located adjacent to the T-variant cleavage site and suggest that T-variant cleavage has taken place inS. mansoni. Sequence heterogeneity in the TSDs indicates that Perere-3 retrotransposition is not site-specific. The TSD sequences contribute to the 5′ end of the terminal ribozyme helix (P1 stem). Based on these results we conclude that T-variants have a functional role in Perere-3 retrotransposition.

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.

Transposons-Based Clonal Diversity in Trematode Involves Parts of CR1 (LINE) in Eu- and Heterochromatin

Trematode parthenitae have long been believed to form clonal populations, but clonal diversity has been discovered in this asexual stage of the lifecycle. Clonal polymorphism in the model species Himasthla elongata has been previously described, but the source of this phenomenon remains unknown. In this work, we traced cercarial clonal diversity using a simplified amplified fragment length polymorphism (SAFLP) method and characterised the nature of fragments in diverse electrophoretic bands. The repetitive elements were identified in both the primary sequence of the H. elongata genome and in the transcriptome data. Long-interspersed nuclear elements (LINEs) and long terminal repeat retrotransposons (LTRs) were found to represent an overwhelming majority of the genome and the transposon transcripts. Most sequenced fragments from SAFLP pattern contained the reverse transcriptase (RT, ORF2) domains of LINEs, and only a few sequences belonged to ORFs of LTRs and ORF1 of LINEs. A fragment corresponding to a CR1-like (LINE) spacer region was discovered and named CR1-renegade (CR1-rng). In addition to RT-containing CR1 transcripts, we found short CR1-rng transcripts in the redia transcriptome and short contigs in the mobilome. Probes against CR1-RT and CR1-rng presented strikingly different pictures in FISH mapping, despite both being fragments of CR1. In silico data and Southern blotting indicated that CR1-rng is not tandemly organised. CR1 involvement in clonal diversity is discussed.

Transposable Elements in the Genome of Human Parasite Schistosoma mansoni: A Review

Transposable elements (TEs) are DNA sequences able to transpose within the host genome and, consequently, influence the dynamics of evolution in the species. Among the possible effects, TEs insertions may alter the expression and coding patterns of genes, leading to genomic innovations. Gene-duplication events, resulting from DNA segmental duplication induced by TEs transposition, constitute another important mechanism that contributes to the plasticity of genomes. This review aims to cover the current knowledge regarding TEs in the genome of the parasite Schistosoma mansoni, an agent of schistosomiasis-a neglected tropical disease affecting at least 250 million people worldwide. In this context, the literature concerning TEs description and TEs impact on the genomic architecture for S. mansoni was revisited, displaying evidence of TEs influence on schistosome speciation-mediated by bursts of transposition-and in gene-duplication events related to schistosome-host coevolution processes, as well several instances of TEs contribution into the coding sequences of genes. These findings indicate the relevant role of TEs in the evolution of the S. mansoni genome.

Repeated Evolution Versus Common Ancestry: Sex Chromosome Evolution in the Haplochromine Cichlid Pseudocrenilabrus philander

Why sex chromosomes turn over and remain undifferentiated in some taxa, whereas they degenerate in others, is still an area of ongoing research. The recurrent occurrence of homologous and homomorphic sex chromosomes in distantly related taxa suggests their independent evolution or continued recombination since their first emergence. Fishes display a great diversity of sex-determining systems. Here, we focus on sex chromosome evolution in haplochromines, the most species-rich lineage of cichlid fishes. We investigate sex-specific signatures in the Pseudocrenilabrus philander species complex, which belongs to a haplochromine genus found in many river systems and ichthyogeographic regions in northern, eastern, central, and southern Africa. Using whole-genome sequencing and population genetic, phylogenetic, and read-coverage analyses, we show that one population of P. philander has an XX-XY sex-determining system on LG7 with a large region of suppressed recombination. However, in a second bottlenecked population, we did not find any sign of a sex chromosome. Interestingly, LG7 also carries an XX-XY system in the phylogenetically more derived Lake Malawi haplochromine cichlids. Although the genomic regions determining sex are the same in Lake Malawi cichlids and P. philander, we did not find evidence for shared ancestry, suggesting that LG7 evolved as sex chromosome at least twice in haplochromine cichlids. Hence, our work provides further evidence for the labile nature of sex determination in fishes and supports the hypothesis that the same genomic regions can repeatedly and rapidly be recruited as sex chromosomes in more distantly related lineages.

Metagenomic analysis of orange colored protrusions from the muscle of Queen Conch Lobatus gigas (Linnaeus, 1758)

The endangered marine gastropod, Lobatus gigas, is an important fishery resource in the Caribbean region. Microbiological and parasitological research of this species have been poorly addressed despite its role in ecological fitness, conservation status and prevention of potential pathogenic infections. This study identified taxonomic groups associated with orange colored protrusions in the muscle of queen conchs using histological analysis, 454 pyrosequencing, and a combination of PCR amplification and automated Sanger sequencing. The molecular approaches indicate that the etiological agent of the muscle protrusions is a parasite belonging to the subclass Digenea. Additionally, the scope of the molecular technique allowed the detection of bacterial and fungi clades in the assignment analysis. This is the first evidence of a digenean infection in the muscle of this valuable Caribbean resource.

Molecular cloning of Kuruma shrimp Marsupenaeus japonicus endonuclease-reverse transcriptase and its positive role in white spot syndrome virus and Vibrio alginolyticus infection

This study investigated the function of endonuclease-reverse transcriptase (mjERT) in Marsupenaeus japonicus. The 1129 bp cDNA sequence of mjERT was cloned from M. japonicus using rapid amplification of cDNA ends (RACE) PCR, and RT-qPCR analysis indicated that mjERT was highly expressed in the gills and hepatopancreas of M. japonicus. We also found that white spot syndrome virus (WSSV) or Vibrio alginolyticus challenge could enhance the expression of mjERT. When mjERT was inhibited, immune genes such as toll, p53, hemocyanin and tumor necrosis factor-α (TNF-α) were significantly down-regulated (P < .01) in the hemocytes of shrimp, while myosin was significantly up-regulated (P < .01). We demonstrated that mjERT is very important for the progression of WSSV infection and that the cumulative mortality of WSSV-infected and V. alginolyticus-infected shrimps was significantly increased following mjERT RNA interfere (RNAi). Apoptosis data provided information to suggest that mjERT-dsRNA challenge caused less apoptosis in hemocytes in both the disease-free and viral group. We also revealed that mjERT-dsRNA treatment resulted in a lower phagocytosis rate in the hemocytes of V. alginolyticus-challenged shrimp. Finally, we found that the absence of mjERT had an significantly negative impact upon shrimp phenoloxidase (PO) activity, superoxide dismutase (SOD) activity and total hemocyte count (THC) following WSSV or V. alginolyticus infection, indicating a regulative role for mjERT in the innate immunity of shrimp in response to pathogenic infection. In summary, we concluded that mjERT might promote the anti-WSSV immune response of shrimp by regulating apoptosis, PO activity, THC and SOD activity, and also exert a positive role in the immune response against V. alginolyticus by regulating phagocytosis, SOD activity, PO activity and THC.

Schistosoma mansoni Egg, Adult Male and Female Comparative Gene Expression Analysis and Identification of Novel Genes by RNA-Seq

Background

Schistosomiasis is one of the most prevalent parasitic diseases worldwide and is a public health problem. Schistosoma mansoni is the most widespread species responsible for schistosomiasis in the Americas, Middle East and Africa. Adult female worms (mated to males) release eggs in the hepatic portal vasculature and are the principal cause of morbidity. Comparative separate transcriptomes of female and male adult worms were previously assessed with using microarrays and Serial Analysis of Gene Expression (SAGE), thus limiting the possibility of finding novel genes. Moreover, the egg transcriptome was analyzed only once with limited bacterially cloned cDNA libraries.

Methodology/Principal findings

To compare the gene expression of S. mansoni eggs, females, and males, we performed RNA-Seq on these three parasite forms using 454/Roche technology and reconstructed the transcriptome using Trinity de novo assembly. The resulting contigs were mapped to the genome and were cross-referenced with predicted Smp genes and H3K4me3 ChIP-Seq public data. For the first time, we obtained separate, unbiased gene expression profiles for S. mansoni eggs and female and male adult worms, identifying enriched biological processes and specific enriched functions for each of the three parasite forms. Transcripts with no match to predicted genes were analyzed for their protein-coding potential and the presence of an encoded conserved protein domain. A set of 232 novel protein-coding genes with putative functions related to reproduction, metabolism, and cell biogenesis was detected, which contributes to the understanding of parasite biology.

Conclusions/Significance

Large-scale RNA-Seq analysis using de novo assembly associated with genome-wide information for histone marks in the vicinity of gene models constitutes a new approach to transcriptome analysis that has not yet been explored in schistosomes. Importantly, all data have been consolidated into a UCSC Genome Browser search- and download-tool (http://schistosoma.usp.br/). This database provides new ways to explore the schistosome genome and transcriptome and will facilitate molecular research on this important parasite.

Schistosomiasis is a public health problem caused by parasites of the genus Schistosoma, of which S. mansoni is the primary causative agent. The parasite has a complex life cycle; their sexual reproductive stage is dependent on female and male adult worms mating inside the mesenteric circulation of the human host, with the female releasing hundreds of eggs daily. This phase of the life cycle is responsible for the development of pathology, which is proportional to the number of eggs accumulating in the liver and intestine of the human host. Genome and transcriptome sequencing of this parasite represent important advances in schistosome research, but there is still a need for integrated analyses to better understand the biology of this parasite. In this study, we describe the first large-scale transcriptomes of eggs, and female and male adult worms, the parasite forms that are mainly responsible for the pathology of schistosomiasis. We were able to cross-reference the gene transcription regions with promoter regions, thus improving the gene annotations. Moreover, we identified the expression of novel protein-coding genes not yet described in the current genome annotation, advancing the biological knowledge regarding this parasite.

Accelerated evolution of schistosome genes coding for proteins located at the host-parasite interface

Study of proteins located at the host–parasite interface in schistosomes might provide clues about the mechanisms utilized by the parasite to escape the host immune system attack. Micro-exon gene (MEG) protein products and venom allergen-like (VAL) proteins have been shown to be present in schistosome secretions or associated with glands, which led to the hypothesis that they are important components in the molecular interaction of the parasite with the host. Phylogenetic and structural analysis of genes and their transcripts in these two classes shows that recent species-specific expansion of gene number for these families occurred separately in three different species of schistosomes. Enrichment of transposable elements in MEG and VAL genes in Schistosoma mansoni provides a credible mechanism for preferential expansion of gene numbers for these families. Analysis of the ratio between synonymous and nonsynonymous substitution rates (dN/dS) in the comparison between schistosome orthologs for the two classes of genes reveals significantly higher values when compared with a set of a control genes coding for secreted proteins, and for proteins previously localized in the tegument. Additional analyses of paralog genes indicate that exposure of the protein to the definitive host immune system is a determining factor leading to the higher than usual dN/dS values in those genes. The observation that two genes encoding S. mansoni vaccine candidate proteins, known to be exposed at the parasite surface, also display similar evolutionary dynamics suggests a broad response of the parasite to evolutionary pressure imposed by the definitive host immune system.

Endonuclease-containing Penelope retrotransposons in the bdelloid rotifer Adineta vaga exhibit unusual structural features and play a role in expansion of host gene families

Background

Penelope-like elements (PLEs) are an enigmatic group of retroelements sharing a common ancestor with telomerase reverse transcriptases. In our previous studies, we identified endonuclease-deficient PLEs that are associated with telomeres in bdelloid rotifers, small freshwater invertebrates best known for their long-term asexuality and high foreign DNA content. Completion of the high-quality draft genome sequence of the bdelloid rotifer Adineta vaga provides us with the opportunity to examine its genomic transposable element (TE) content, as well as TE impact on genome function and evolution.

Results

We performed an exhaustive search of the A. vaga genome assembly, aimed at identification of canonical PLEs combining both the reverse transcriptase (RT) and the GIY-YIG endonuclease (EN) domains. We find that the RT/EN-containing Penelope families co-exist in the A. vaga genome with the EN-deficient RT-containing Athena retroelements. Canonical PLEs are present at very low copy numbers, often as a single-copy, and there is no evidence that they might preferentially co-mobilize EN-deficient PLEs. We also find that Penelope elements can participate in expansion of A. vaga multigene families via trans-action of their enzymatic machinery, as evidenced by identification of intron-containing host genes framed by the Penelope terminal repeats and characteristic target-site duplications generated upon insertion. In addition, we find that Penelope open reading frames (ORFs) in several families have incorporated long stretches of coding sequence several hundred amino acids (aa) in length that are highly enriched in asparagine residues, a phenomenon not observed in other retrotransposons.

Conclusions

Our results show that, despite their low abundance and low transcriptional activity in the A. vaga genome, endonuclease-containing Penelope elements can participate in expansion of host multigene families. We conclude that the terminal repeats represent the cis-acting sequences required for mobilization of the intervening region in trans by the Penelope-encoded enzymatic activities. We also hypothesize that the unusual capture of long N-rich segments by the Penelope ORF occurs as a consequence of peculiarities of its replication mechanism. These findings emphasize the unconventional nature of Penelope retrotransposons, which, in contrast to all other retrotransposon types, are capable of dispersing intron-containing genes, thereby questioning the validity of traditional estimates of gene retrocopies in PLE-containing eukaryotic genomes.

Endonuclease domain of the Drosophila melanogaster R2 non-LTR retrotransposon and related retroelements: a new model for transposition

The molecular mechanisms of the transposition of non-long terminal repeat (non-LTR) retrotransposons are not well understood; the key questions of how the 3′-ends of cDNA copies integrate and how site-specific integration occurs remain unresolved. Integration depends on properties of the endonuclease (EN) domain of retrotransposons. Using the EN domain of the Drosophila R2 retrotransposon as a model for other, closely related non-LTR retrotransposons, we investigated the EN domain and found that it resembles archaeal Holliday-junction resolvases. We suggest that these non-LTR retrotransposons are co-transcribed with the host transcript. Combined with the proposed resolvase activity of the EN domain, this model yields a novel mechanism for site-specific retrotransposition within this class of retrotransposons, with resolution proceeding via a Holliday junction intermediate.

Identification and characterization of microRNAs and endogenous siRNAs in Schistosoma japonicum

BACKGROUND

Small endogenous non-coding RNAs (sncRNAs) such as small interfering RNA (siRNA), microRNA and other small RNA transcripts are derived from distinct loci in the genome and play critical roles in RNA-mediated gene silencing mechanisms in plants and metazoa. They are approximately 22 nucleotides long; regulate mRNA stability through perfect or imperfect match to the targets. The biological activities of sncRNAs have been related to many biological events, from resistance to microbe infections to cellular differentiation. The development of the zoonotic parasite Schistosoma japonicum parasite includes multiple steps of morphological alterations and biological differentiations, which provide a unique model for studies on the functions of small RNAs. Characterization of the genome-wide transcription of the sncRNAs will be a major step in understanding of the parasite biology. The objective of this study is to investigate the transcriptional profile and potential function of the small non-coding RNAs in the development of S. japanicum.

RESULTS

The endogenous siRNAs were found mainly derived from transposable elements (TE) or transposons and the natural antisense transcripts (NAT). In contrast to other organisms, the TE-derived siRNAs in S. japonicum were more predominant than other sncRNAs including microRNAs (miRNAs). Further, there were distinct length and 3'end variations in the sncRNAs, which were associated with the developmental differentiation of the parasite. Among the identified miRNA transcripts, there were 38 unique to S. japonicum and 16 that belonged to 13 miRNA families are common to other metazoan lineages. These miRNAs were either ubiquitously expressed, or they exhibited specific expression patterns related to the developmental stages or sex. Genes that encoded miRNAs are mainly located in clusters within the genome of S. japonicum. However, genes within one cluster could be differentially transcribed, which suggested that individual genes might be regulated by distinct mechanisms during parasite development.

CONCLUSIONS

Many miRNA and endogenous siRNA transcripts were identified in S. japonicum and the amount of siRNA was at least 4.4 and 1.6 times more than that of miRNA in both schistosomulum and adult worm stages respectively. SiRNAs are mainly derived from transposable elements (or transposons); while natural antisense transcripts (NAT)-derived siRNAs were much less. A majority of miRNA transcripts identified in the parasite were species-specific and the expression of certain miRNAs was found developmentally regulated. Both miRNA and siRNAs are potentially important regulators in the development of schistosomal parasites.

Bursts of transposition from non-long terminal repeat retrotransposon families of the RTE clade in Schistosoma mansoni

The genus Schistosoma is composed of blood flukes that infect vertebrates, from which three species are major causative agents of human schistosomiasis, a tropical disease that affects more than 200 million people. Current models of the recent evolution of Schistosoma indicate multiple events of migration and speciation from an Asian ancestral species. Transposable elements are important drivers of genome evolution and have been hypothesised to have an important role in speciation. In this work, we describe a comprehensive inventory of Schistosoma mansoni and Schistosoma japonicum retrotransposons, based on their recently published genomic data. We find a considerable difference in retrotransposon representation between the two species (22% and 13%, respectively). A large part of this difference can be attributed to higher representation of two previously described families of S. mansoni retrotransposons (SR2 and Perere-3/SR3), compared with the representation of their closest relative families in S. japonicum. A more detailed analysis suggests that these two S. mansoni families were the subject of recent bursts of transposition that were not paralleled by their S. japonicum counterparts. We hypothesise that these bursts could be a consequence of the evolutionary pressure resulting from migration of Schistosoma from Asia to Africa and their establishment in this new environment, helping both speciation and adaptation.

The genome of the blood fluke Schistosoma mansoni

Schistosoma mansoni is responsible for the neglected tropical disease schistosomiasis that affects 210 million people in 76 countries. We report here analysis of the 363 megabase nuclear genome of the blood fluke. It encodes at least 11,809 genes, with an unusual intron size distribution, and novel families of micro-exon genes that undergo frequent alternate splicing. As the first sequenced flatworm, and a representative of the lophotrochozoa, it offers insights into early events in the evolution of the animals, including the development of a body pattern with bilateral symmetry, and the development of tissues into organs. Our analysis has been informed by the need to find new drug targets. The deficits in lipid metabolism that make schistosomes dependent on the host are revealed, while the identification of membrane receptors, ion channels and more than 300 proteases, provide new insights into the biology of the life cycle and novel targets. Bioinformatics approaches have identified metabolic chokepoints while a chemogenomic screen has pinpointed schistosome proteins for which existing drugs may be active. The information generated provides an invaluable resource for the research community to develop much needed new control tools for the treatment and eradication of this important and neglected disease.

Identification and isolation of a retrotransposon from the freshwater sponge Lubomirskia baicalensis: implication in rapid evolution of endemic sponges

Transposons are mobile genetic elements that are found in all major branches of life. Similarities to retroviruses concerning genome structure and transposition mechanism suggest a familial relationship. Transposons are important evolutionary drivers that trigger genetic changes such as genomic rearrangement, alteration of gene expression, and gene duplication. And, indeed, now more than ever the effect of transposons on genome evolution represents a dynamic field of research. Since sponges (phylum Porifera) are the phylogenetically oldest still extant metazoan taxon, the study of poriferan mobile elements contributes to the understanding of the generation of phenotypic diversity and speciation at the base of the metazoan tree of life. This work describes the analyses of the first poriferan mobile genetic element so far identified, the long terminal repeats- retrotransposon Baikalum-1 of Lubomirskia baicalensis (Demospongiae; Ceractinomorpha). Baikalum-1 embraces a continuous open reading frame, putatively coding for a polyprotein that consists of nucleo capsid, protease, reverse transcriptase, RNase H, and integrase, all proteins/ enzymes characteristic of retrotransposons. Baikalum-1 was discovered in all freshwater sponge species endemic to Lake Baikal, as well as in cosmopolitan sponge species that inhabit a Lake Baikal-feeding rivulet. However, the same cosmopolitan species sampled from lakes and rivers (Siberian and European) with no direct contact to Lake Baikal did not contain this particular mobile genetic element. Thus, Baikalum-1 is probably the result of an evolutionarily ancient retroviral infection that spread exclusively amongst Baikalian sponge species. In addition, the retro-transposon is found in the vicinity of the silicatein-A1 gene. Silicateins are cathepsin-like proteins that catalyze the synthesis of poriferan siliceous skeletal elements (spicules). In L. baicalensis, the silicatein-A1 gene is flanked by two palindroms, probably remnants of transposons that might be responsible for the emergence of four different silicatein genes, uniquely present in freshwater but not marine sponges. Adaptation of sponges to the freshwater habitat (with its significantly higher silica content compared to the marine milieu) required the ability to evolve rapidly, which could be conferred by high transpositional activity, accompanied by duplication and diversification of the ancestral silicatein gene of marine species.

The Brazilian contribution to the study of the Schistosoma mansoni transcriptome

Brazilian scientists have made definitive contributions towards the understanding of the transcriptome of Schistosoma mansoni. The transcriptome, as a collection of expressed genes, was studied with two basic approaches, the sequencing of cloned cDNA ends and the production and sequencing of random fragments of cDNA. The data that was generated covers nearly the entire set of transcribed genes. These data have been important for gene discovery and annotation of gene function, annotation of the genome, construction of microarrays, identification of peptides in the study of the proteome and can be mined for splice variants, polymorphisms such as single nucleotide polymorphisms and microsatellites, for example. This review will describe the advances achieved with the contribution of Brazilian researchers to the understanding of the transcribed regions of the S. mansoni genome.

Evidence of multiple retrotransposons in two litopenaeid species

Retrotransposons encompass a specific class of mobile genetic elements that are widespread across eukaryotic genomes. The impact of the varied types of retrotransposons on these genomes is just beginning to be deciphered. In a step towards understanding their role in litopenaeid shrimp, we have herein identified nine non-LTR retrotransposons, among which several appear to exist outside the standard defined clades. Two Litopenaeus stylirostris elements were discovered through degenerate PCR amplification using previously defined non-LTR degenerate primers, and through primers designed from a RAPD-derived sequence. A third genomic L. stylirostris element was identified using specific priming from an amplification protocol. These three PCR-derived sequences showed conserved domains of the non-LTR reverse transcriptase gene. In silico searching of genome databases and subsequent contig construction yielded six non-LTR retrotransposons (both genomic and expressed) in the Litopenaeus vannamei genome that also exhibited the highly conserved domains found in our PCR-derived sequences. Phylogenetic placement among representatives from all non-LTR clades showed a possibly novel monophyletic group that included five of our nine sequences. This group, which included elements from both L. stylirostris and L. vannamei, appeared most closely related to the highly active RTE clade. Our remaining four sequences placed in the CR1 and I clades of retrotransposons, with one showing strong similarity to ancient Penelope elements. This research describes three newly discovered retrotransposons in the L. stylirostris genome. Phylogenetic analysis clusters these in a monophyletic grouping with retrotransposons previously described from two closely related species, L. vannamei and Penaeus monodon.

Physical localization of the retrotransposons Boudicca and Perere 03 in Schistosoma mansoni

Schistosoma mansoni is 1 of the causative agents of schistosomiasis, an endemic disease in 76 countries of the world. The study of its genome, estimated to be 270 Mb, is very important to understanding schistosome biology, the mechanisms of drug resistance, and immune evasion. Repetitive elements constitute more than 40% of the S. mansoni genome and may play a role in the parasite evolution. The retrotransposons Boudicca, a long terminal repeat (LTR), and Perere 03, a non-LTR, are present in a high number in the S. mansoni genome and were localized with the use of fluorescence in situ hybridization (FISH) and primed in situ labeling (PRINS). Bacterial artificial chromosomes (BAC) clones containing the retrotransposons Boudicca and Perere 03 were selected by bioinformatic analysis and used as probes in FISH. Using metaphase chromosomes from sporocysts and the FISH and PRINS techniques, we were able to map these retrotransposons. Perere 03 was localized in the euchromatic regions of the short arm of chromosome 2 and Boudicca in the euchromatic regions of the short arm of chromosomes 2 and Z.

Integration of reporter transgenes into Schistosoma mansoni chromosomes mediated by pseudotyped murine leukemia virus

The recent release of draft genome sequences of two of the major human schistosomes has underscored the pressing need to develop functional genomics approaches for these significant pathogens. The sequence information also makes feasible genome-scale investigation of transgene integration into schistosome chromosomes. Retrovirus-mediated transduction offers a means to establish transgenic lines of schistosomes, to elucidate schistosome gene function and expression, and to advance functional genomics approaches for these parasites. We investigated the utility of the Moloney murine leukemia retrovirus (MLV) pseudotyped with vesicular stomatitis virus glycoprotein (VSVG) for the transduction of Schistosoma mansoni and delivery of reporter transgenes into schistosome chromosomes. Schistosomula were exposed to virions of VSVG-pseudotyped MLV, after which genomic DNA was extracted from the transduced schistosomes. Southern hybridization analysis indicated the presence of proviral MLV retrovirus in the transduced schistosomes. Fragments of the MLV transgene and flanking schistosome sequences recovered using an anchored PCR-based approach demonstrated definitively that somatic transgenesis of schistosome chromosomes had taken place and, moreover, revealed widespread retrovirus integration into schistosome chromosomes. More specifically, MLV transgenes had inserted in the vicinity of genes encoding immunophilin, zinc finger protein Sma-Zic, and others, as well as near the endogenous schistosome retrotransposons, the fugitive and SR1. Proviral integration of the MLV transgene appeared to exhibit primary sequence site specificity, targeting a gGATcc-like motif. Reporter luciferase transgene activity driven by the schistosome actin gene promoter was expressed in the tissues of transduced schistosomula and adult schistosomes. Luciferase activity appeared to be developmentally expressed in schistosomula with increased activity observed after 1 to 2 wk in culture. These findings indicate the utility of VSVG-pseudotyped MLV for transgenesis of S. mansoni, herald a tractable pathway forward toward germline transgenesis and functional genomics of parasitic helminths, and provide the basis for comparative molecular pathogenesis studies of chromosomal lesions arising from retroviral integration into human compared with schistosome chromosomes.

Transgenesis of schistosomes: approaches employing mobile genetic elements

Draft genome sequences for Schistosoma mansoni and Schistosoma japonicum are now available. However, the identity and importance of most schistosome genes have yet to be determined. Recently, progress has been made towards the genetic manipulation and transgenesis of schistosomes. Both loss-of-function and gain-of-function approaches appear to be feasible in schistosomes based on findings described in the past 5 years. This review focuses on reports of schistosome transgenesis, specifically those dealing with the transformation of schistosomes with exogenous mobile genetic elements and/or their endogenous relatives for the genetic manipulation of schistosomes. Transgenesis mediated by mobile genetic elements offers a potentially tractable route to introduce foreign genes to schistosomes, a means to determine the importance of schistosome genes, including those that could be targeted in novel interventions and the potential to undertake large-scale forward genetics by insertional mutagenesis.

Use of a 44k oligoarray to explore the transcriptome of Schistosoma mansoni adult worms

Recent advances in the study of Schistosoma mansoni genome and transcriptome have led to a better description of the S. mansoni gene complement. In this work, we report the design and use of a new S. mansoni 60-mer oligonucleotide microarray platform with approximately 44,000 probes, based on all publicly available cDNA sequence data for S. mansoni and Schistosoma japonicum. The large number of probes combined with the extensive sequence annotation available allowed a comprehensive approach, where most of the S. mansoni transcriptome is represented. Hybridization with adult worm RNA pointed to a set of genes transcriptionally active in this stage of the parasite's life cycle. Interestingly, a large proportion (43%) of genes for which transcription was detected in adults is comprised of "no match" genes, i.e. S. mansoni genes with unknown function and no identifiable orthologs in GenBank. Moreover, detection of bi-directional transcription for 7% of the active "no match" genes in adults leads us to hypothesize a widespread production of antisense RNA in S. mansoni, with possible regulatory roles.

Schistosoma mansoni: the dicer gene and its expression

RNA interference (RNAi) is a gene silencing mechanism that plays an important role in regulating gene expression in many eukaryotes and has become a valuable molecular tool for analyzing gene function. Multi-domain nucleases called Dicer proteins play pivotal roles in RNAi. In this paper, we characterize the structure and expression of the Dicer gene from the platyhelminth parasite Schistosoma mansoni. The gene (SmDicer) is over 54kb long and comprises 30 exons that potentially encode a 2641 amino acid protein. This is the largest Dicer protein yet described. SmDicer contains all domains that are characteristic of metazoan dicers including an amino terminal helicase domain, DUF283, a PAZ domain, two RNAse III domains and an RNA binding domain. An examination of the available S. mansoni genome sequence suggests that the Dicer gene described here is the only Dicer gene in the parasite genome. SmDicer is expressed throughout schistosome development suggesting that RNAi technologies might be employed in deciphering gene function in all life stages of this parasite.

piggyBac transposon mediated transgenesis of the human blood fluke, Schistosoma mansoni

The transposon piggyBac from the genome of the cabbage looper moth Trichoplusia ni has been observed in the laboratory to jump into the genomes of key model and pathogenic eukaryote organisms including mosquitoes, planarians, human and other mammalian cells, and the malaria parasite Plasmodium falciparum. Introduction of exogenous transposons into schistosomes has not been reported but transposon-mediated transgenesis of schistosomes might supersede current methods for functional genomics of this important human pathogen. In the present study we examined whether the piggyBac transposon could deliver reporter transgenes into the genome of Schistosoma mansoni parasites. A piggyBac donor plasmid modified to encode firefly luciferase under control of schistosome gene promoters was introduced along with 7-methylguanosine capped RNAs encoding piggyBac transposase into cultured schistosomula by square wave electroporation. The activity of the helper transposase mRNA was confirmed by Southern hybridization analysis of genomic DNA from the transformed schistosomes, and hybridization signals indicated that the piggyBac transposon had integrated into numerous sites within the parasite chromosomes. piggyBac integrations were recovered by retrotransposon-anchored PCR, revealing characteristic piggyBac TTAA footprints in the vicinity of the endogenous schistosome retrotransposons Boudicca, SR1, and SR2. This is the first report of chromosomal integration of a transgene and somatic transgenesis of this important human pathogen, in this instance accomplished by mobilization of the piggyBac transposon.

Distribution and Phylogeny of Penelope-Like Elements in Eukaryotes

Penelope-like elements (PLEs) are a relatively little studied class of eukaryotic retroelements, distinguished by the presence of the GIY-YIG endonuclease domain, the ability of some representatives to retain introns, and the similarity of PLE-encoded reverse transcriptases to telomerases. Although these retrotransposons are abundant in many animal genomes, the reverse transcriptase moiety can also be found in several protists, fungi, and plants, indicating its ancient origin. A comprehensive phylogenetic analysis of PLEs was conducted, based on extended sequence alignments and a considerably expanded data set. PLEs exhibit the pattern of evolution similar to that of non-LTR retrotransposons, which form deep-branching clades dating back to the Precambrian era. However, PLEs seem to have experienced a much higher degree of lineage losses than non-LTR retrotransposons. It is suggested that PLEs and non-LTR retrotransposons are included into a larger eTPRT (eukaryotic target-primed) group of retroelements, characterized by 5' truncation, variable target-site duplication, and the potential of the 3' end to participate in formation of non-autonomous derivatives.

'Oming in on schistosomes: prospects and limitations for post-genomics

The recent release of version 3 of the Schistosoma mansoni genome assembly has made a wealth of information available to researchers. Here, progress made in schistosome genomics and post-genomics is considered. The current status of knowledge about the genome, transcriptome, proteome, glycome and immunome is summarized and recent publications briefly reviewed. The prospects for advances in understanding schistosome biology are highlighted. Most importantly, the limitations (which are mostly technical) that need to be addressed before the full potential of the genome database(s) can be realized are emphasized.

SmTRC1, a novel Schistosoma mansoni DNA transposon, discloses new families of animal and fungi transposons belonging to the CACTA superfamily

Background

The CACTA (also called En/Spm) superfamily of DNA-only transposons contain the core sequence CACTA in their Terminal Inverted Repeats (TIRs) and so far have only been described in plants. Large transcriptome and genome sequence data have recently become publicly available for Schistosoma mansoni, a digenetic blood fluke that is a major causative agent of schistosomiasis in humans, and have provided a comprehensive repository for the discovery of novel genes and repetitive elements. Despite the extensive description of retroelements in S. mansoni, just a single DNA-only transposon belonging to the Merlin family has so far been reported in this organism.

Results

We describe a novel S. mansoni transposon named SmTRC1, for S. mansoni Transposon Related to CACTA 1, an element that shares several characteristics with plant CACTA transposons. Southern blotting indicates approximately 30–300 copies of SmTRC1 in the S. mansoni genome. Using genomic PCR followed by cloning and sequencing, we amplified and characterized a full-length and a truncated copy of this element. RT-PCR using S. mansoni mRNA followed by cloning and sequencing revealed several alternatively spliced transcripts of this transposon, resulting in distinct ORFs coding for different proteins. Interestingly, a survey of complete genomes from animals and fungi revealed several other novel TRC elements, indicating new families of DNA transposons belonging to the CACTA superfamily that have not previously been reported in these kingdoms. The first three bases in the S. mansoni TIR are CCC and they are identical to those in the TIRs of the insects Aedes aegypti and Tribolium castaneum, suggesting that animal TRCs may display a CCC core sequence.

Conclusion

The DNA-only transposable element SmTRC1 from S. mansoni exhibits various characteristics, such as generation of multiple alternatively-spliced transcripts, the presence of terminal inverted repeats at the extremities of the elements flanked by direct repeats and the presence of a Transposase_21 domain, that suggest a distant relationship to CACTA transposons from Magnoliophyta. Several sequences from other Metazoa and Fungi code for proteins similar to those encoded by SmTRC1, suggesting that such elements have a common ancestry, and indicating inheritance through vertical transmission before separation of the Eumetazoa, Fungi and Plants.

Cloning the genes and DNA binding properties of High Mobility Group B1 (HMGB1) proteins from the human blood flukes Schistosoma mansoni and Schistosoma japonicum

The parasitic helminth Schistosoma mansoni contains three HMGB proteins, HMGB1, HMGB2 and HMGB3, of primary amino acid sequences highly similar to vertebrate proteins. In this report we describe the characterization of the HMGB1 proteins and their genes from S. mansoni and Schistosoma japonicum. The deduced amino acid sequences of HMGB1 proteins from both schistosome species are identical, and comprise 176 residues. The proteins contain the two evolutionarily highly conserved HMG-box domains, A and B, exhibiting 60% similarity to mammalian HMGB1. Unlike the human HMGB1 which contains an unbroken run of 30 glutamic or aspartic residues, the SmHMGB1 or SjHMGB1 proteins possess unusually short acidic C-terminal tails (5 acidic residues interrupted by 2 serines). Southern hybridization and DNA sequencing revealed a single copy HMGB1 gene, composed of 3 exons and two introns, in S. mansoni. The exon/intron boundaries are identical to those of the human HMGB1 gene, with the exception that the second exon of the SmHMGB1 gene which is not split into two exons as in the human HMGB1 gene. RNA blot analysis revealed that the SmHMGB1 gene is constitutively expressed in similar levels both in male and female worms. The single-sized mRNA for SmHMGB1 is consistent with the size derived from the cDNA. Although DNA binding properties of SmHMGB1 (or SjHMGB1) protein seem to be similar to those previously reported with human HMGB1, i.e., preferential binding to supercoiled DNA over linear DNA, specific recognition of DNA four-way junctions, DNA-induced supercoiling in the presence of topoisomerase I, and DNA bending, we have observed two important differences relative to those observed with the human HMGB1: (i) the inability of the isolated SmHMGB1 domain A to bend DNA (as revealed by T4 ligase-mediated circularization assay), and (ii) higher DNA supercoiling and bending potential of the SmHMGB1 protein as compared to its human counterpart. The latter finding may indicate that the long acidic C-tail of human HMGB1 has much stronger repressive role on DNA bending or DNA supercoiling by topoisomerase I at physiological ionic strength than the short C-tail of the SmHMGB1 protein. Considering the important role of HMGB1 in DNA replication, transcription, recombination, and in particularly, the mediation of inflammation responses in mammalian cells, further studies on schistosome HMGB proteins may provide valuable information related to schistosomiasis, where inflammation plays a critical role in this disease.

Identification of rDNA-specific non-LTR retrotransposons in Cnidaria

Ribosomal RNA genes are abundant repetitive sequences in most eukaryotes. Ribosomal DNA (rDNA) contains many insertions derived from mobile elements including non-long terminal repeat (non-LTR) retrotransposons. R2 is the well-characterized 28S rDNA-specific non-LTR retrotransposon family that is distributed over at least 4 bilaterian phyla. R2 is a large family sharing the same insertion specificity and classified into 4 clades (R2-A, -B, -C, and -D) based on the N-terminal domain structure and the phylogeny. There is no observation of horizontal transfer of R2; therefore, the origin of R2 dates back to before the split between protostomes and deuterostomes. Here, we in silico identified 1 R2 element from the sea anemone Nematostella vectensis and 2 R2-like retrotransposons from the hydrozoan Hydra magnipapillata. R2 from N. vectensis was inserted into the 28S rDNA like other R2, but the R2-like elements from H. magnipapillata were inserted into the specific sequence in the highly conserved region of the 18S rDNA. We designated the Hydra R2-like elements R8. R8 is inserted at 37 bp upstream from R7, another 18S rDNA-specific retrotransposon family. There is no obvious sequence similarity between targets of R2 and R8, probably because they recognize long DNA sequences. Domain structure and phylogeny indicate that R2 from N. vectensis is the member of the R2-D clade, and R8 from H. magnipapillata belongs to the R2-A clade despite its different sequence specificity. These results suggest that R2 had been generated before the split between cnidarians and bilaterians and that R8 is a retrotransposon family that changed its target from the 28S rDNA to the 18S rDNA.

Parasite histories and novel phylogenetic tools: Alternative approaches to inferring parasite evolution from molecular markers

Parasitological research is often contingent on the knowledge of the phylogeny/genealogy of the studied group. Although molecular phylogenetics has proved to be a powerful tool in such investigations, its application in the traditional fashion, based on a tree inference from the primary nucleotide sequences may, in many cases, be insufficient or even improper. These limitations are due to a number of factors, such as a scarcity/ambiguity of phylogenetic information in the sequences, an intricacy of gene relationships at low phylogenetic levels, or a lack of criteria when deciding among several competing coevolutionary scenarios. With respect to the importance of a precise and reliable phylogenetic background in many biological studies, attempts are being made to extend molecular phylogenetics with a variety of new data sources and methodologies. In this review, selected approaches potentially applicable to parasitological research are presented and their advantages as well as drawbacks are discussed. These issues include the usage of idiosyncratic markers (unique features with presumably low probability of homoplasy), such as insertion of mobile elements, gene rearrangements and secondary structure features; the problem of ancestral polymorphism and reticulate relationships at low phylogenetic levels; and the utility of a molecular clock to facilitate discrimination among alternative scenarios in host-parasite coevolution.

Characterization of SR3 reveals abundance of non-LTR retrotransposons of the RTE clade in the genome of the human blood fluke, Schistosoma mansoni

BACKGROUND

It is becoming apparent that perhaps as much as half of the genome of the human blood fluke Schistosoma mansoni is constituted of mobile genetic element-related sequences. Non-long terminal repeat (LTR) retrotransposons, related to the LINE elements of mammals, comprise much of this repetitive component of the schistosome genome. Of more than 12 recognized clades of non-LTR retrotransposons, only members of the CR1, RTE, and R2 clades have been reported from the schistosome genome.

RESULTS

Inspection of the nucleotide sequence of bacterial artificial chromosome number 49_J_14 from chromosome 1 of the genome of Schistosoma mansoni (GenBank AC093105) revealed the likely presence of several RTE-like retrotransposons. Among these, a new non-LTR retrotransposon designated SR3 was identified and is characterized here. Analysis of gene structure and phylogenetic analysis of both the reverse transcriptase and endonuclease domains of the mobile element indicated that SR3 represented a new family of RTE-like non-LTR retrotransposons. Remarkably, two full-length copies of SR3-like elements were present in BAC 49-J-14, and one of 3,211 bp in length appeared to be intact, indicating SR3 to be an active non-LTR retrotransposon. Both were flanked by target site duplications of 10-12 bp. Southern hybridization and bioinformatics analyses indicated the presence of numerous copies (probably >1,000) of SR3 interspersed throughout the genome of S. mansoni. Bioinformatics analyses also revealed SR3 to be transcribed in both larval and adult developmental stages of S. mansoni and to be also present in the genomes of the other major schistosome parasites of humans, Schistosoma haematobium and S. japonicum.

CONCLUSION

Numerous copies of SR3, a novel non-LTR retrotransposon of the RTE clade are present in the genome of S. mansoni. Non-LTR retrotransposons of the RTE clade including SR3 appear to have been remarkably successful in colonizing, and proliferation within the schistosome genome.