Jannotti Passos, 2010
Jannotti Passos, 2010
Jannotti Passos, 2010
Mem Inst Oswaldo Cruz, Rio de Janeiro, Vol. 105(4): 504-511, July 2010
Mitochondrial DNA of Biomphalaria tenagophila, a mollusc intermediate host of Schistosoma mansoni in Brazil,
was sequenced and characterised. The genome size found for B. tenagophila was 13,722 bp and contained 13 mes
senger RNAs, 22 transfer RNAs (tRNA) and two ribosomal RNAs (rRNA). In addition to sequencing, the
mitochon
drial DNA (mtDNA) genome organization of B. tenagophila was analysed based on its content and localization of
both coding and non-coding regions, regions of gene overlap and tRNA nucleotide sequences. Sequences of protein,
rRNA 12S and rRNA 16S nucleotides as well as gene organization were compared between B. tenagophila and Biomphalaria glabrata, as the latter is the most important S. mansoni intermediate host in Brazil. Differences between
such species were observed regarding rRNA composition. The complete sequence of the B. tenagophila mitochon
drial genome was deposited in GenBank (accession EF433576). Furthermore, phylogenetic relationships were esti
mated among 28 mollusc species, which had their complete mitochondrial genome deposited in GenBank, using the
neighbour-joining method, maximum parsimony and maximum likelihood bootstrap. B. tenagophila was positioned
at a branch close to B. glabrata and Pulmonata molluscs, collectively comprising a paraphyletic group, contrary to
Opistobranchia, which was positioned at a single branch and constituted a monophyletic group.
Key words: Biomphalaria tenagophila - mitochondrial genome - phylogenetic
phylum has a long geological history and approximately 35,000 species are known. Biom
phalaria (Gastropoda) molluscs are known to have originated during the Jurassic era and endured a wide range of
environmental changes, resulting in the development of
survival mechanisms such as self-fecundation, diapause,
aestivation, capacity to survive into deep trenches and
high prolificacy (Paraense 1955), making them currently
adapted to different environmental conditions.
Biomphalaria tenagophila is an important species
concerning schistosomiasis transmission in Southern
Brazil. This species has been the focus of investigations,
as it comprises susceptible (Paraense & Corra 1978)
populations and one population that is resistant to Schis
tosoma mansoni, which is from the Biological Reserve
of Taim in Rio Grande do Sul, Brazil.
These molluscs have already been challenged by different S. mansoni strains and variable number of miracidia, yet remained resistant to infection (Santos et al.
505
Molecular cloning of PCR products (COI-16S, 16SCyb, Cyb-12S, 12S-COI) - Molecular fragments generated by PCR were cloned in different plasmids according to size: the 2,200 bp fragment (COI-16S) was cloned
in TOPO TA cloning (Invitrogen - Life Technologies),
the 2,800 bp fragment (Cyb-12S) was cloned in pGEM-T
EASY (Promega) and the fragments of 4,500 bp (16SCyb e 12S-COI) in TOPO XL PCR cloning (Invitrogen
- Life Technologies). The protocol followed the manufacturers instructions. The TOPO 10 F` strain from
Escherichia coli bacteria was transformed with recombinant plasmids.
Genome assembly and sequencing analysis - The
shotgun sequencing reads obtained from 521 DNA templates, including those obtained using the 54 primers
specifically designed for primer walking were analysed
using PHRED/PHRAP/CONSED (Ewing et al. 1998, Ewing & Green 1998, Gordon et al. 1998) in order to read
sequencer trace data, write the base calls and quality values, assemble shotgun DNA sequence data and to visually
inspect, edit and finish the alignments. A PHRED quality
value of 20, which corresponds to an error probability of
1/100 in the base call at each position of the read and a
minimal overlap of 50 bp with 92% identity were used as
initial constraints in the clustering for genome assembly.
Tandem Repeat Finder software was used to find and
analyse repeats in the B. tenagophila assembled genome.
Repeat Masker (AFA Smit & P Green, unpublished observations) together with Repbase databank (http://www.
girinst.org) were used to identify known repeats and
mask low complexity regions in the mollusc genome.
Analyses of the B. tenagophila assembled genome
using different algorithms such as BLAST and FASTA
were carried out using different databases including:
the protein database from National Center for Biotechnology Information (NCBI), the database of Clusters of
Orthologous Groups of proteins (http://www.ncbi.nlm.
nih.gov/COG/), a simple modular architecture research
tool (SMART) (http://smart.emblheidelberg.de/), protein families database of alignments and HMMs (http://
www.sanger.ac.uk/Software/Pfam/) and Gene Ontology (http://www.geneontology.org/). Locally compiled
mollusca and ribosomal RNA DBs were built by parsing public databases.
The annotation and graphical output of the B. tena
gophila genomic sequence was performed with an ARTEMIS (http://www.sanger.ac.uk/Software/Artemis/) and
in-house developed PERL scripts to analyse and format
the results. All in silico searches and analyses described
above, together with specific pipelines and PERL scripts
developed during this work, were run in a local server using in-house copies of databases and software.
The search for tRNA genes in genomic sequences
was performed using a combinatorial approach using
ab initio predictions made by tRNAscan-SE software
together with similarity searches against a locally compiled tRNA database.
The visualisation of comparative genome analysis between B. tenagophila and B. glabrata was performed using ACT software (http://www.sanger.ac.uk/Software/).
506
Mem Inst Oswaldo Cruz, Rio de Janeiro, Vol. 105(4), July 2010
TABLE I
Specific primers
Primers
16SL1
16SL2
COIL2
COIR4
CybLF
CybLR
12SLR
12SLF
Sequences
5GCACAGTCAAGGTACTGCGGCCATTTAAAAAAAATATCACTGG
5AATCTCAAATGTATCCTGCCCAGTGATATTTTTTAAATGG
5GGCCTGGTTGGAACTGGGTTATCTTTACTTATTCGTTTAGAGC
5CCAAGCTCTAAACGAATAAGTAAAGATAACCCAGTTCCAACCAG
5AAAGGAAGAACTAAGCCTTTAGGAGATTTAAACCATTTGG
5TTAAATGGGATTTTTTCCAAATGGTTTAAATCTCCTAAAGG
5ATAATTAATGAATTTCTTTGATTTTGTGCACACATTGCCC
5AAAGAAATTCATTAATTATATATTTTTTAGTTACTTTTAAGTCC
RESULTS
mtDNA content - B. tenagophila mtDNA was completely sequenced and characterised. The genome size
encompasses 13,722 bp. This genome codes for two rRNAs, 13 protein genes and a set of 22 tRNAs. The proportion of AT was 75.78%. Nine non-coding regions ranging
from 1-54 bp consisted of a total of 110 bp. The largest
non-coding regions were located between the NAD1 and
NAD4L genes and between COIII and tRNAIle, consisting of 54 and 42 bp, respectively. Nineteen overlapping
regions were found in B. tenagophila mtDNA and the
largest overlap was found between NAD4L and Cytb
genes comprising 29 bp.
The tRNAscan software identified six out of the 22
tRNA (Asp, Asn, Thr, Val, Ala, Pro). The other tRNAs
were identified by nucleotide alignment with its correspondent in B. glabrata and B. tenagophila tRNA had an
average size of 63 bp.
The genes found in the mitochondrial genome coding for proteins were analysed based on the mitochondrial genetic code of invertebrates. Most protein-coding
genes were located on the plus strand (NAD5, NAD1,
NAD4L, Cytb, COII, NAD4, NAD2, COI, NAD6).
Genes coding the ATP6, ATP8, NAD3 and COIII proteins were located on the minus strand.
Regarding protein size, COII, NAD2, NAD4 and
NAD4L were longer in B. tenagophila than in B. glabra
ta. The protein with the largest size difference between
the two species proteins was NAD1, which was 11 additional amino acids shorter in B. glabrata. The sizes of
NAD1 and NAD6 were smaller in B. tenagophila than in
B. glabrata. For ATP8, ATP6, NAD2, NAD4 and COIII,
the termination codon is likely completed by adding an
A during mRNA transcription. The least conserved protein was ATP8, with 60% identity in the two species and
the most conserved protein was COI, with 98% identity
between these two species.
Four mtDNA initiation codons were found in B.
tenagophila, namely: ATT (in 5 proteins), ATA (in 5
proteins), ATG (in 2 protein) and TTG (in 1 protein).
The stop codons were TAA (in 6 proteins) and TAG (in
2 proteins), which have been described for mitochondrial genes (Wolstenholme 1992). Incomplete stop codons
were found in five proteins (ATP6, ATP8, COI, COII,
NAD3, NAD4).
507
TABLE II
Molluscs used in the phylogenetic study
Mollusc
Accession (GenBank)
Abbreviation
Lampsilis ornata
Haliotis rubra
Lophiotoma cerithiformis
Nautilus macromphalus
Lottia digitalis
Ilyanassa obsoleta
Aplysia californica
Octopus ocellatus
Crassostrea virginica
Mytilus trossulus
Watasenia scintillans
Sepioteuthis lessoniana
Sepia officinalis
Roboastra europaea
Placopecten magellanicus
Mytilus galloprovincialis
Albinaria coerulea
Todarodes pacificus
Octopus vulgaris
Biomphalaria glabrata
Graptacme eborea
Mytilus edulis
Siphonodentalium lobatum
Crassostrea gigas
Cepaea nemoralis
Loligo bleekeri
Pupa strigosa
Katharina tunicata
AY365193
AY588938
DQ284754
DQ472026
DQ238599
DQ238598
AY569552
AB240156
AY905542
AY823625
AB240152
AB240154
AB240155
AY083457
DQ088274
AY497292
X83390
AB158364
AB158363
AY380531
AY484748
AY484747
AY342055
AF177226
U23045
AB029616
AB028237
U09810
Lam orn
Hal rub
Lop cer
Nau mac
Lot dig
Ily obs
Apl cal
Oct oce
Cra vir
Myt tro
Wat sci
Sep les
Sep off
Rob eur
Pla mag
Myt gal
Alb coe
Tod pac
Oct vul
Bio gla
Gra ebo
Myt edu
Sip lob
Cra gig
Cep nem
Lol ble
Pup str
Kat tun
Phylogenomic analysis - mtDNA amino acid sequences of 29 organisms from the Mollusca phylum were
aligned (except for ATP8) and used for phylogenomic
positioning analysis of B. tenagophila and B. glabrata in
comparison with other molluscs.
Regardless the phylogenetic methodology, the trees
showed similar topology, comprising five groups. Group
I: MP bootstrap, 94% (Fig. 2), ML bootstrap, 99% (Fig.
3) encompassed six molluscs of the Bivalvia class (Protobranchia sub class, Pteriomorphia super order), one
species of the Gastropoda class (Prosobranchia) and
one species of the Scaphopoda class (Gadelida). In the
NJ tree (Fig. 4), this group clustered the same molluscs
when compared with MP and ML trees, except for the
508
Mem Inst Oswaldo Cruz, Rio de Janeiro, Vol. 105(4), July 2010
Fig. 1: gene rearrangements in the Biomphalaria tenagophila and Biomphalaria glabrata mitochondrial genomes. No transfer RNAs (tRNA)
was found among ND6, ND5, ND1, ND4L and Cytb genes, but the following tRNAs were identified: between Cytb and COII genes for Asp, Cys
and Phe; between COII and ATP8, six tRNAs (Tyr, Trp, Gly, His, Gln, Leu); between ATP8 and ATP6, tRNA Asn; between ATP6 and rRNA 12S,
two tRNAs (Arg and Glu); between rRNA 12S and NAD3, tRNA Met; between NAD3 and NAD4, two tRNAs for Ser; between NAD4 and COIII,
tRNAThr; between COIII and NAD2, tRNA Ile; between NAD2 and COI, tRNA Lys; between COI and rRNA 16S, tRNAVal and between rRNA 16S
and NAD6, three tRNAs (Leu, Ala, Pro).
DISCUSSION
509
Fig. 2: the phylogeny of the mollusc. The tree was inferred from the
well aligned positions of the mitochondrial protein sequences (excluding ATP8), using maximum parsimony method. The bootstrap
values are shown in the tree for 100 replicates.
Fig. 3: the phylogeny of the mollusc. The tree was inferred from the
well aligned positions of the mitochondrial protein sequences (excluding ATP8), using maximum likelihood method. The bootstrap values
are shown in the tree for 100 replicates.
510
Mem Inst Oswaldo Cruz, Rio de Janeiro, Vol. 105(4), July 2010
de Freitas JR, Boschi MB, dos Santos MB 1985. Suscetibilidade
de hbridos de Biomphalaria tenagophila cepa LE (BH) do
Schistosoma mansoni. Rev Inst Med Trop Sao Paulo 27: 6-12.
DeJong RJ, Emery AM, Adema CM 2004. The mitochondrial genome
of Biomphalaria glabrata (Gastropoda, Basommatophora), intermediate host of Schistosoma mansoni. J Parasitol 90: 991-997.
Ewing B, Green P 1998. Base-calling of automated sequencer traces
using phred. II. Error probabilities. Genome Res 8: 186-194.
Ewing B, Hillier L, Wendl MC, Green P 1998. Base-calling of automated sequencer traces using phred. I. Accuracy assessment.
Genome Res 8: 175-185.
Felsenstein J 1985. Confidence limits on phylogenies: an approach
using the bootstrap. Evolution 39: 783-791.
Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R 1994. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Bio
technol 3: 294-299.
Foster PG, Jermiin LS, Hickey DA 1997. Nucleotide composition bias
affects amino acid content in proteins coded by animal mitochondria. J Mol Evol 44: 282-288.
Gordon D, Abajian C, Green P 1998. Consed: a graphical tool for sequence finishing. Genome Res 8: 195-202.
Grande C, Templado J, Cervera JL, Zardoya R 2004. Molecular phylogeny
of euthyneura (Mollusca: Gastropoda). Mol Biol Evol 21: 303-313.
Fig. 4: the phylogeny of the mollusc. The tree was inferred from the
well aligned positions of the mitochondrial protein sequences (excluding ATP8), using neighbor-joining method. The bootstrap values
are shown in the tree for 100 replicates.
REFERENCES
Avise JC 1986. Mitochondrial DNA and the evolutionary genetics
of higher animals. Philos Trans R Soc Lond B Biol Sci 312:
325-342.
Barbosa L 2001. Transferncia de resistncia ao Schistosoma mansoni em Biomphalaria tenagophila por transplante do rgo he
matopoitico de linhagem resistente (Taim, RS) para linhagem
suscetvel (Cabo Frio, RJ), MSc Thesis, Universidade Federal
de Minas Gerais, Belo Horizonte, 40 pp.
Barbosa L, Caldeira RL, Carvalho OS, Vidigal TH, Jannotti-Passos
LK, Coelho PM 2006. Resistance to Schistosoma mansoni by
transplantation of APO Biomphalaria tenagophila. Parasite Im
munol 28: 209-212.
Bezerra FSM, Nogueira-Machado JA, Martins-Souza RL, Chaves
MM, Correa RF, Coelho PMZ 2003. Effect of gamma radiation
on the activity of hemocytes and on the course of Schistosoma
mansoni infection in resistant Biomphalaria tenagophila snails.
Mem Inst Oswaldo Cruz 98: 73-75.
Boore JL, Brown WM 1994. Mitochondrial genomes and phylogeny
of mollusks. Nautilus 108 (Suppl. 2): 61-78.
Brusca RC, Brusca GJ 2007. Invertebrados, 2nd ed., Guanabara
Koogan, Rio de Janeiro, 1098 pp.
Cheng S, Chang SY, Gravitt P, Respess R 1994. Long PCR. Nature
369: 684-685.
Coelho JR, Bezerra FSM 2006. Compatibility of Biomphalaria
tenagophila with Schistosoma mansoni: a study of homologous
plasma transference. Mem Inst Oswaldo Cruz 101: 111-112.
511
Simon C, Paabo S, Kocher TD, Wilson AC 1990. Evolution of mitochondrial ribosomal RNA in insects as shown by the polymerase chain reaction. In M Clegg, S OBrien (eds.), Molecular
evolution, UCLA Symposia on Molecular and Cellular Biology, Liss, p. 235-244.
Rosa FM, Caldeira RL, Carvalho OS, Godard ALB, Coelho PMZ 2004.
Dominant character of the molecular marker of Biomphalaria
tenagophila strain (Mollusca: Planorbidae) resistant to Schisto
soma mansoni. Mem Inst Oswaldo Cruz 99: 85-87.