Clustering of Schistosoma mansoni mRNA sequences and analysis of the most transcribed genes: implications in metabolism and biology of different developmental stages

The study of the Schistosoma mansoni genome, one of the etiologic agents of human schistosomiasis, is essential for a better understanding of the biology and development of this parasite. In order to get an overview of all S. mansoni catalogued gene sequences, we performed a clustering analysis of the parasite mRNA sequences available in public databases. This was made using softwares PHRAP and CAP3. The consensus sequences, generated after the alignment of cluster constituent sequences, allowed the identification by database homology searches of the most expressed genes in the worm. We analyzed these genes and looked for a correlation between their high expression and parasite metabolism and biology. We observed that the majority of these genes is related to the maintenance of basic cell functions, encoding genes whose products are related to the cytoskeleton, intracellular transport and energy metabolism. Evidences are presented here that genes for aerobic energy metabolism are expressed in all the developmental stages analyzed. Some of the most expressed genes could not be identified by homology searches and may have some specific functions in the parasite.

Characterization of new Schistosoma mansoni microsatellite loci in sequences obtained from public DNA databases and microsatellite enriched genomic libraries

In the last decade microsatellites have become one of the most useful genetic markers used in a large number of organisms due to their abundance and high level of polymorphism. Microsatellites have been used for individual identification, paternity tests, forensic studies and population genetics. Data on microsatellite abundance comes preferentially from microsatellite enriched libraries and DNA sequence databases. We have conducted a search in GenBank of more than 16,000 Schistosoma mansoni ESTs and 42,000 BAC sequences. In addition, we obtained 300 sequences from CA and AT microsatellite enriched genomic libraries. The sequences were searched for simple repeats using the RepeatMasker software. Of 16,022 ESTs, we detected 481 (3%) sequences that contained 622 microsatellites (434 perfect, 164 imperfect and 24 compounds). Of the 481 ESTs, 194 were grouped in 63 clusters containing 2 to 15 ESTs per cluster. Polymorphisms were observed in 16 clusters. The 287 remaining ESTs were orphan sequences. Of the 42,017 BAC end sequences, 1,598 (3.8%) contained microsatellites (2,335 perfect, 287 imperfect and 79 compounds). The 1,598 BAC end sequences 80 were grouped into 17 clusters containing 3 to 17 BAC end sequences per cluster. Microsatellites were present in 67 out of 300 sequences from microsatellite enriched libraries (55 perfect, 38 imperfect and 15 compounds). From all of the observed loci 55 were selected for having the longest perfect repeats and flanking regions that allowed the design of primers for PCR amplification. Additionally we describe two new polymorphic microsatellite loci.

The E.coli RuvAB proteins branch migrate Holliday junctions through heterologous DNA sequences in a reaction facilitated by SSB.

During genetic recombination a heteroduplex joint is formed between two homologous DNA molecules. The heteroduplex joint plays an important role in recombination since it accommodates sequence heterogeneities (mismatches, insertions or deletions) that lead to genetic variation. Two Escherichia coli proteins, RuvA and RuvB, promote the formation of heteroduplex DNA by catalysing the branch migration of crossovers, or Holliday junctions, which link recombining chromosomes. We show that RuvA and RuvB can promote branch migration through 1800 bp of heterologous DNA, in a reaction facilitated by the presence of E.coli single-stranded DNA binding (SSB) protein. Reaction intermediates, containing unpaired heteroduplex regions bound by SSB, were directly visualized by electron microscopy. In the absence of SSB, or when SSB was replaced by a single-strand binding protein from bacteriophage T4 (gene 32 protein), only limited heterologous branch migration was observed. These results show that the RuvAB proteins, which are induced as part of the SOS response to DNA damage, allow genetic recombination and the recombinational repair of DNA to occur in the presence of extensive lengths of heterology.

Biased gene conversion and GC-content evolution in the coding sequences of reptiles and vertebrates.

Mammalian and avian genomes are characterized by a substantial spatial heterogeneity of GC-content, which is often interpreted as reflecting the effect of local GC-biased gene conversion (gBGC), a meiotic repair bias that favors G and C over A and T alleles in high-recombining genomic regions. Surprisingly, the first fully sequenced nonavian sauropsid (i.e., reptile), the green anole Anolis carolinensis, revealed a highly homogeneous genomic GC-content landscape, suggesting the possibility that gBGC might not be at work in this lineage. Here, we analyze GC-content evolution at third-codon positions (GC3) in 44 vertebrates species, including eight newly sequenced transcriptomes, with a specific focus on nonavian sauropsids. We report that reptiles, including the green anole, have a genome-wide distribution of GC3 similar to that of mammals and birds, and we infer a strong GC3-heterogeneity to be already present in the tetrapod ancestor. We further show that the dynamic of coding sequence GC-content is largely governed by karyotypic features in vertebrates, notably in the green anole, in agreement with the gBGC hypothesis. The discrepancy between third-codon positions and noncoding DNA regarding GC-content dynamics in the green anole could not be explained by the activity of transposable elements or selection on codon usage. This analysis highlights the unique value of third-codon positions as an insertion/deletion-free marker of nucleotide substitution biases that ultimately affect the evolution of proteins.

Analysis of the first and second internal transcribed spacer sequences of the ribosomal DNA in Biomphalaria tenagophila complex (Mollusca: Planorbidae)

The first and second internal transcribed spacer regions (ITS1 and ITS2) of the ribosomal DNA of Biomphalaria tenagophila complex (B. tenagophila, B. occidentalis, and B. t. guaibensis) were sequenced and compared. The alignment lengths of these regions were about 655 bp and 481 bp, respectively. Phylogenetic relationships among the Biomphalaria species were inferred by Maximum Parsimony and Neighbor-joining methods. The phylogenetic trees produced, in most of the cases, were in accordance with morphological systematics and other molecular data previously obtained by polymerase chain reaction and restriction fragment length polymorphism analysis. The present results provide support for the proposal that B. tenagophila represents a complex comprising B. tenagophila, B. occidentalis and B. t. guaibensis.

Evolution of vitellogenin genes: comparative analysis of the nucleotide sequences downstream of the transcription initiation site of four Xenopus laevis and one chicken gene.

Electron microscopic analysis of heteroduplexes between the most distantly related Xenopus vitellogenin genes (A genes X B genes) has revealed the distribution of homologous regions that have been preferentially conserved after the duplication events that gave rise to the multigene family in Xenopus laevis. DNA sequence analysis was limited to the region downstream of the transcription initiation site of the Xenopus genes A1, B1 and B2 and a comparison with the Xenopus A2 and the major chicken vitellogenin gene is presented. Within the coding regions of the first three exons, nucleotide substitutions resulting in amino acid changes accumulate at a rate similar to that observed in globin genes. This suggests that the duplication event which led to the formation of the A and B ancestral genes in Xenopus laevis occurred about 150 million years ago. Homologous exons of the A1-A2 and B1-B2 gene pairs, which formed about 30 million years ago, show a quite similar sequence divergence. In contrast, A1-A2 homologous introns seem to have evolved much faster than their B1-B2 counterparts.

Identification of estrogen-responsive DNA sequences by transient expression experiments in a human breast cancer cell line.

The expression of a hybrid gene formed by the promoter region of the Xenopus laevis vitellogenin gene B1 and the CAT coding region is regulated by estrogen when the gene is transfected into hormone-responsive MCF-7 cells. Furthermore, the 5' flanking region of the gene B1 alone can confer inducibility to heterologous promoters, although to a varying extent depending on the promoter used. Deletion mapping of he vitellogenin hormone-responsive sequences revealed that a 13 bp element 5'-AGTCACTGTGACC-3' at position -334 is essential for estrogen inducibility. We have shown previously that this 13 bp element is present upstream of several liver-specific estrogen-inducible genes.

Temporal patterns of fire sequences observed in Canton of Ticino (southern Switzerland)

Temporal dynamical analysis in fire sequences recorded from 1969 to 2008 in Canton Ticino (Switzerland) was carried out by using the Allan Factor statistics. The obtained results show the presence of daily periodicities, superimposed to two time-scaling regimes. The daily cycle vanishes for sequences of higher altitude fires, for which a single scaling behaviour is observed.

Athletic injuries of the extensor carpi ulnaris subsheath: MRI findings and utility of gadolinium-enhanced fat-saturated T1-weighted sequences with wrist pronation and supination

Purpose: To report the magnetic resonance imaging (MRI) findings in athletic injuries of the extensor carpi ulnaris (ECU) subsheath, assessing the utility of gadolinium-enhanced (Gd) fat-saturated (FS) T1-weighted sequences with wrist pronation and supination. Methods and Materials: Sixteen patients (13 males, 3 females; mean age 30.3 years) with athletic injuries of the ECU subsheath sustained between January 2003 and June 2009 were included in this retrospective study. Initial and follow‑up 1.5-T wrist MRIs were performed with transverse T1-weighted and STIR sequences in pronation, and Gd FS T1-weighted sequences with wrist pronation and supination. Two radiologists assessed the type of injury (A to C), ECU tendon stability, associated lesions and rated pulse sequences using a three-point scale: 1 = poor, 2 = good and 3 = excellent. Results: Gd-enhanced FS T1-weighted transverse sequences in supination (2.63) and pronation (2.56) were most valuable, compared with STIR (2.19) and T1 weighted (1.94). Nine type A, one type B and six type C injuries were found. There were trends towards diminution in size, signal intensity and enhancement of associated pouches on follow‑up MRI and tendon stabilisation within the ulnar groove. Conclusion: Gd-enhanced FS T1-weighted sequences with wrist pronation and supination are most valuable in assessing and follow‑up athletic injuries of the ECU subsheath on 1.5-T MRI.

Production of full-length cDNA sequences by sequencing and analysis of expressed sequence tags from Schistosoma mansoni

The number of sequences generated by genome projects has increased exponentially, but gene characterization has not followed at the same rate. Sequencing and analysis of full-length cDNAs is an important step in gene characterization that has been used nowadays by several research groups. In this work, we have selected Schistosoma mansoni clones for full-length sequencing, using an algorithm that investigates the presence of the initial methionine in the parasite sequence based on the positions of alignment start between two sequences. BLAST searches to produce such alignments have been performed using parasite expressed sequence tags produced by Minas Gerais Genome Network against sequences from the database Eukaryotic Cluster of Orthologous Groups (KOG). This procedure has allowed the selection of clones representing 398 proteins which have not been deposited as S. mansoni complete CDS in any public database. Dedicated sequencing of 96 of such clones with reads from both 5' and 3' ends has been performed. These reads have been assembled using PHRAP, resulting in the production of 33 full-length sequences that represent novel S. mansoni proteins. These results shall contribute to construct a more complete view of the biology of this important parasite.

Challenges of phylogenetic analyses of aDNA sequences

One of the crucial steps of authentication of aDNA sequences is phylogenetic consistency. Amplified sequences should fit into the phylogenetic framework of their supposed origin. An inherent property of aDNA sequences however, is their short sequence length. Additionally, genes for aDNA studies are often chosen by their preservation potential rather than by phylogenetically informative content. This poses potential challenges regarding their analyses, and might result in an inaccurate reflection of the supposed phylogenetic history of the sequence or organism under study. In this paper some fundamental problems of phylogenetic analysis and interpretation of aDNA datasets are discussed. Suggestions for character sampling and treatment of missing data are made. The publication is the result of a talk from the 1st PAMINSA Meeting in Rio de Janeiro, July 2005.

Re-mapping the molecular features of the human immunodeficiency virus type 1 and human T-cell lymphotropic virus type 1 Brazilian sequences using a bioinformatics unit established in Salvador, Bahia, Brazil, to give support to the viral epidemiology studies

The analysis of genetic data for human immunodeficiency virus type 1 (HIV-1) and human T-cell lymphotropic virus type 1 (HTLV-1) is essential to improve treatment and public health strategies as well as to select strains for vaccine programs. However, the analysis of large quantities of genetic data requires collaborative efforts in bioinformatics, computer biology, molecular biology, evolution, and medical science. The objective of this study was to review and improve the molecular epidemiology of HIV-1 and HTLV-1 viruses isolated in Brazil using bioinformatic tools available in the Laboratório Avançado de Sáude Pública (Lasp) bioinformatics unit. The analysis of HIV-1 isolates confirmed a heterogeneous distribution of the viral genotypes circulating in the country. The Brazilian HIV-1 epidemic is characterized by the presence of multiple subtypes (B, F1, C) and B/F1 recombinant virus while, on the other hand, most of the HTLV-1 sequences were classified as Transcontinental subgroup of the Cosmopolitan subtype. Despite the high variation among HIV-1 subtypes, protein glycosylation and phosphorylation domains were conserved in the pol, gag, and env genes of the Brazilian HIV-1 strains suggesting constraints in the HIV-1 evolution process. As expected, the functional protein sites were highly conservative in the HTLV-1 env gene sequences. Furthermore, the presence of these functional sites in HIV-1 and HTLV-1 strains could help in the development of vaccines that pre-empt the viral escape process.

Mapping the molecular characteristics of Brazilian human T-cell lymphotropic virus type 1 Env (gp46) and Pol amino acid sequences for vaccine design

This study was carried out to evaluate the molecular pattern of all available Brazilian human T-cell lymphotropic virus type 1 Env (n = 15) and Pol (n = 43) nucleotide sequences via epitope prediction, physico-chemical analysis, and protein potential sites identification, giving support to the Brazilian AIDS vaccine program. In 12 previously described peptides of the Env sequences we found 12 epitopes, while in 4 peptides of the Pol sequences we found 4 epitopes. The total variation on the amino acid composition was 9 and 17% for human leukocyte antigen (HLA) class I and class II Env epitopes, respectively. After analyzing the Pol sequences, results revealed a total amino acid variation of 0.75% for HLA-I and HLA-II epitopes. In 5 of the 12 Env epitopes the physico-chemical analysis demonstrated that the mutations magnified the antigenicity profile. The potential protein domain analysis of Env sequences showed the loss of a CK-2 phosphorylation site caused by D197N mutation in one epitope, and a N-glycosylation site caused by S246Y and V247I mutations in another epitope. Besides, the analysis of selection pressure have found 8 positive selected sites (w = 9.59) using the codon-based substitution models and maximum-likelihood methods. These studies underscore the importance of this Env region for the virus fitness, for the host immune response and, therefore, for the development of vaccine candidates.