ChromaPipe: a pipeline for analysis, quality control and management for a DNA sequencing facility.

Optimizing and monitoring the data flow in high-throughput sequencing facilities is important for data input and output, for tracking the status of results for the users of the facility, and to guarantee a good, high-quality service. In a multi-user system environment with different throughputs, each user wants to access his/her data easily, track his/her sequencing history, analyze sequences and their quality, and apply some basic post-sequencing analysis, without the necessity of installing further software. Recently, Fiocruz established such a core facility as a "technological platform". Infrastructure includes a 48-capillary 3730 DNA Sequence Analyzer (Applied Biosystems) and supporting equipment. The service includes running samples for large-scale users, performing DNA sequencing reactions and runs for medium and small users, and participation in partial or full genome projects. We implemented a workflow that fulfills these requirements for small and high throughput users. Our implementation also includes the monitoring of data for continuous quality improvement (reports by plate, month and user) by the sequencing staff. For the user, different analyses of the chromatograms, such as visualization of good quality regions, as well as processing, such as comparisons or assemblies, are available. So far, 180 users have made use of the service, generating 155,000 sequences, 35% of which were produced for the BCG Moreau-RJ genome project. The pipeline (named ChromaPipe for Chromatogram Pipeline) is available for download by the scientific community at the url http://bioinfo.pdtis.fiocruz.br/ChromaPipe/. The support for assembly is also configured as a web service: http://bioinfo.pdtis.fiocruz.br/Assembly/.

In silico reconstruction of the amino acid metabolic pathways of Trypanosoma cruzi.

Trypanosoma cruzi is the epidemiological agent of Chagas' disease, affecting most of Central and South America, constituting a significant health and socio-economic problem. The parasite has a metabolism largely based on the consumption of amino acids, which participate in a diversity of metabolic pathways, leading to many crucial compounds for the survival of this parasite. Study of its enzymes has the potential to disclose new therapeutic targets and foster the development of new drugs. In this study, we employed computational approaches to reconstruct in silico the amino acid metabolic pathways of T. cruzi, aiming to link genomic information with functional information. For that, protein sequences from 570 EC classes belonging to 25 different pathways in general amino acid metabolism were downloaded from KEGG. A subset of 471 EC classes had at least one sequence deposited. Clustering of the proteins belonging to each EC class was performed using a similarity-based approach implemented in the tool AnEnPi. Reconstruction of the metabolic pathways comprising the amino acid metabolism of T. cruzi was performed by analyzing the output of BLASTP, using as query the dataset of predicted proteins of T. cruzi against all sequences of each individual cluster. This approach allowed us to identify 764 T. cruzi proteins probably involved in the metabolism of amino acids as well as the identification of several putative cases of analogy. Furthermore, we were able to identify several enzymatic activities of T. cruzi that were not previously included in KEGG.

Comparative genomics of two Leptospira interrogans serovars reveals novel insights into physiology and pathogenesis.

Leptospira species colonize a significant proportion of rodent populations worldwide and produce life-threatening infections in accidental hosts, including humans. Complete genome sequencing of Leptospira interrogans serovar Copenhageni and comparative analysis with the available Leptospira interrogans serovar Lai genome reveal that despite overall genetic similarity there are significant structural differences, including a large chromosomal inversion and extensive variation in the number and distribution of insertion sequence elements. Genome sequence analysis elucidates many of the novel aspects of leptospiral physiology relating to energy metabolism, oxygen tolerance, two-component signal transduction systems, and mechanisms of pathogenesis. A broad array of transcriptional regulation proteins and two new families of afimbrial adhesins which contribute to host tissue colonization in the early steps of infection were identified. Differences in genes involved in the biosynthesis of lipopolysaccharide O side chains between the Copenhageni and Lai serovars were identified, offering an important starting point for the elucidation of the organism's complex polysaccharide surface antigens. Differences in adhesins and in lipopolysaccharide might be associated with the adaptation of serovars Copenhageni and Lai to different animal hosts. Hundreds of genes encoding surface-exposed lipoproteins and transmembrane outer membrane proteins were identified as candidates for development of vaccines for the prevention of leptospirosis.

Novel allelic variants of Mycobacteria isolated in Brazil as determined by PCR-restriction enzyme analysis of hsp65.

Human isolates of Mycobacterium collected in 16 different states of Brazil were submitted to PCR-restriction analysis (PRA) of a 439-bp fragment of the hsp65 gene with HaeIII and BstEII. Fourteen allelic variants not described in clinical isolates so far were observed among 36 (10%) of 356 Brazilian strains, including a new pattern for Mycobacterium scrofulaceum, M. intracellulare, and M. flavescens, two new patterns for M. fortuitum, three new patterns each for M. gordonae and M. terrae, and one new pattern for M. avium complex-like strains. Two unidentified strains each also presented a new pattern, strongly suggesting that Mycobacterium genotypes are distributed biogeographically. The PRA procedure was also performed with 43 reference isolates belonging to 34 species, adding a further six new patterns to the identification algorithm. A database containing the normalized restriction patterns of both enzymes was constructed. Patterns available on the Internet can be introduced into this database, which will make possible the comparison of genotypes from isolates from different parts of the world.

Parasite genome initiatives.

During 1993-1994, scientists from developing and developed countries planned and initiated a number of parasite genome projects and several consortiums for the mapping and sequencing of these medium-sized genomes were established, often based on already ongoing scientific collaborations. Financial and other support came from WHO/TDR, Wellcome Trust and other funding agencies. Thus, the genomes of Plasmodium falciparum, Schistosoma mansoni, Trypanosoma cruzi, Leishmania major, Trypanosoma brucei, Brugia malayi and other pathogenic nematodes are now under study. From an initial phase of network formation, mapping efforts and resource building (EST, GSS, phage, cosmid, BAC and YAC library constructions), sequencing was initiated in gene discovery projects but soon also on a small chromosome, and now on a fully fledged genome scale. Proteomics, functional analysis, genetic manipulation and microarray analysis are ongoing to different degrees in the respective genome initiatives, and as the funding for the whole genome sequencing becomes secured, most of the participating laboratories, apart from larger sequencing centres, become oriented to post-genomics. Bioinformatics networks are being expanded, including in developing countries, for data mining, annotation and in-depth analysis.

Detection of Mycobacterium leprae DNA by polymerase chain reaction in the blood of individuals, eight years after completion of anti-leprosy therapy.

Thirty eight patients with indeterminate leprosy (HI), at least 4 to 6 years after discharge from multibacillary (MB) or paucibacillary (PB) schemes of anti leprosy multidrug therapy (MDT), were submitted to traditional diagnostic procedures for leprosy and to polymerase chain reaction (PCR) analysis of different clinical samples for detection of Mycobacterium leprae DNA. No significant difference was observed for any of the parameters analyzed between PB or MB schemes of treatment and no indications were found for more efficient outcome of HI using the MB scheme. Remarkably, 18 (54.5%) of the individuals were PCR positive in at least one of the samples: positivity of PCR was highest in blood samples and four individuals were PCR positive in blood and some other sample. Upon comparison of PCR results with clinical and histopathological parameters, no correlation was found between PCR-positivity and eventual relapse. This is the first report on detection of M. leprae DNA in PB patients, more than half a decade after completion of MDT, suggesting that live bacilli are present and circulating much longer than expected, although reinfection of the individuals can not be excluded. Overall, we feel that because of the high sensitivity of the assay, extreme care should be taken about association of PCR results, efficacy of treatment and disease status.

The heterogeneity of choanomastigote-shaped trypanosomatids as analyzed by their kDNA minicircle size: taxonomic implicaitons.

The kDNA minicircle size was analyzed in 15 species of choanomastigote-shaped trypanosomatids and four main groups of species were identified: (1) "Crithidia" deanei, "C." desouzai and "Herpetomonas" roitmani, which presented the largest molecules (> or = 3,800 bp), (2) "C." oncopelti with minicircles of at least four different sizes within 1,300-2,650 bp, (3) C. fasciculata, C. guilhermei and C. luciliae, having at least two sizes of minicircle (2,650 bp and 3,000 bp) and (4) a heterogeneous group of species presenting minicircles of a single size, including several Crithidia spp. (having 1,600 bp or 1,700 bp minicircles) and two Proteomonas spp. presenting the smallest minicircles (about 1,500 bp). These results were compared with other observations and discussed from a taxonomic point of view.

Usefulness of IS6110-restriction fragment length polymorphism typing of Brazilian strains of Mycobacterium tuberculosis and comparison with an international fingerprint database.

Strains of Mycobacterium tuberculosis isolated from 219 different tuberculosis patients, 115 from patients residing in Rio de Janeiro, 79 from Rio Grande do Sul and the remaining from other regions of the country, were analyzed by IS6110-restriction fragment length polymorphism fingerprinting. The IS6110-DNA patterns from these strains were highly polymorphic: 174 different patterns were observed and 25 patterns were shared by 70 isolates (32%). Most strains (93.4%) had multicopy patterns and only 17% of clustered strains had less than six IS6110 copies. Strain clustering was significantly higher for isolates from Rio Grande do Sul (36.7%) in comparison with strains from Rio de Janeiro (22.6%), but only when using high stringency during cluster analysis. Upon screening of an international database containing 3,970 fingerprints of M. tuberculosis strains, 15% of the patterns of Brazilian strains (21% of the strains) were identical to a fingerprint of an isolate from another country and one particular eight-band pattern forming the largest Brazilian cluster was detected in seven additional countries, suggesting that international transmission of tuberculosis from and to Brazil could be occurring frequently. Alternatively,preferential use of certain IS6110 integration sites could also be important in high-copy number strains, having important consequences for the use of databases for epidemiological studies on a large scale.

Gene expression, amino acid conservation, and hydrophobicity are the main factors shaping codon preferences in Mycobacterium tuberculosis and Mycobacterium leprae.

Mycobacterium tuberculosis and Mycobacterium leprae are the ethiological agents of tuberculosis and leprosy, respectively. After performing extensive comparisons between genes from these two GC-rich bacterial species, we were able to construct a set of 275 homologous genes. Since these two bacterial species also have a very low growth rate, translational selection could not be so determinant in their codon preferences as it is in other fast-growing bacteria. Indeed, principal-components analysis of codon usage from this set of homologous genes revealed that the codon choices in M. tuberculosis and M. leprae are correlated not only with compositional constraints and translational selection, but also with the degree of amino acid conservation and the hydrophobicity of the encoded proteins. Finally, significant correlations were found between GC3 and synonymous distances as well as between synonymous and nonsynonymous distances.

Physical mapping of a 670-kb region of chromosomes XVI and XVII from the human protozoan parasite Trypanosoma cruzi encompassing the genes for two immunodominant antigens.

As part of the Trypanosoma cruzi Genome Initiative, we have mapped a large portion of the chromosomal bands XVI (2.3 Mb) and XVII (2.6 Mb) containing the highly repetitive and immunodominant antigenic gene families h49 and jl8. Restriction mapping of the isolated chromosomal bands and hybridization with chromosome specific gene probes showed that genes h49 and jl8 are located in a pair of size-polymorphic homologous chromosomes. To construct the integrated map of the chromosomes harboring the h49 and jl8 loci, we used YAC, cosmid, and lambda phage overlapping clones, and long range restriction analysis using a variety of probes (i.e., known gene sequences, ESTs, polymorphic repetitive sequences, anonymous sequences, STSs generated from the YAC ends). The total length covered by the YAC contig was approximately 670 kb, and its map agreed and was complementary to the one obtained by long-range restriction fragment analysis. Average genetic marker spacing in a 105 kb region around h49 and jl8 genes was estimated to be 6.2 kb/marker. We have detected some polymorphism in the H49/JL8 antigens-encoding chromosomes, affecting also the coding regions. The physical map of this region, together with the isolation of specific chromosome markers, will contribute in the global effort to sequence the nuclear genome of this parasite.

Use of PCR-restriction fragment length polymorphism analysis of the hsp65 gene for rapid identification of mycobacteria in Brazil.

Polymerase chain reaction amplification of part of the gene coding for the heat shock protein hsp65 followed by restriction enzyme analysis (PRA) is a recently described tool for rapid identification of mycobacteria. In this study, the speed and simplicity of PRA for identification of isolates of mycobacteria from patients with clinical symptoms of tuberculosis was evaluated and compared with identification results obtained by commercially available methods. Established PRA patterns were observed for nineteen isolates of Mycobacterium tuberculosis, eleven belonging to the complex M. avium-intracellulare, four of M. kansasii, one of M. fortuitum, one of M. abscessus, three of M. gordonae and one of the recently described species M. lentiflavum, as identified by commercially available methods. Two isolates of M. fortuitum and one of M. gordonae had unique and so far undescribed PRA patterns, suggesting geographically-related intra-species variation within the hsp65 sequence. We propose the inclusion of these new patterns in the PRA identification algorithm and have defined more accurately the molecular weight values of the restriction fragments. This is the first report on the isolation of M. lentiflavum in Brazil suggesting that identification by means of PRA could be useful for detection of mycobacterial species that are usually unnoticed. Where the use of several commercial techniques in combination was necessary for correct identification, PRA demonstrated to be a simple technique with good cost-benefit for characterization of all mycobacterial isolates in this study.

Cystic fibrosis: low frequency of DF508 mutation in 2 population samples from Rio de Janeiro, Brazil.

Blood samples from 44 unrelated cystic fibrosis (CF) patients from Rio de Janeiro, Brazil, were analyzed for the 8 European CF mutations. Six homozygous and 15 heterozygous carriers of the DF508 mutation were found, corresponding to 47.7% of CF patients (allele frequency 0.3068). The G542X and G551D mutations were also observed with allele frequencies of 0.0227 and 0.0114, respectively. An analysis of the DF508 mutation in 291 randomly chosen, healthy individuals was performed, and only 3 heterozygous carriers were identified. These results show that the frequency of the DF508 allele in Rio de Janeiro is much lower than the world average; this may be due to the extremely heterogeneous ethnic admixture of the study population. By combining the results of these 2 different samples (CF patients and random population) and admixture data from Rio de Janeiro, we can estimate the CF incidence in this population to be 1:3542 individuals. However, taking into account the Rio de Janeiro ethnic admixture, we can find an estimate of 1:6902 individuals.

Genomics and the biology of parasites.

Despite the advances of modern medicine, the threat of chronic illness, disfigurement, or death that can result from parasitic infection still affects the majority of the world population, retarding economic development. For most parasitic diseases, current therapeutics often leave much to be desired in terms of administration regime, toxicity, or effectiveness and potential vaccines are a long way from market. Our best prospects for identifying new targets for drug, vaccine, and diagnostics development and for dissecting the biological basis of drug resistance, antigenic diversity, infectivity and pathology lie in parasite genome analysis, and international mapping and gene discovery initiatives are under way for a variety of protozoan and helminth parasites. These are far from ideal experimental organisms, and the influence of biological and genomic characteristics on experimental approaches is discussed, progress is reviewed and future prospects are examined.

DNA fingerprinting of Mycobacterium tuberculosis from patients with and without AIDS in Rio de Janeiro.

Isolates of Mycobacterium tuberculosis derived from patients with AIDS from a single hospital in Rio de Janeiro were typed using a standardized RFLP technique detecting IS6110 polymorphism. Nineteen isolates were obtained from 15 different patients. Eleven distinct IS6110 patterns were found, with 4 banding patterns shared by 2 patients. The clustering value of 53% was much higher in comparison with clustering of M. tuberculosis strains from TB patients without clinical signs for HIV infection from randomly selected health centers. We present these results as preliminary data on M. tuberculosis strain polymorphism in Brazil and on the higher risk for recent transmission amongst patients with AIDS.

DNA fingerprinting of mycobacterium tuberculis from patients with and without AIDS in Rio de Janeiro

Isolates of Mycobacterium tuberculosis derived from patients with AIDS from a single hospital in Rio de Janeiro were typed using a standardized RFLP technique detecting IS6110 polymorphism. Nineteen isolates were obtained from 15 different patients. Eleven distinct IS6110 patterns were found, with 4 banding patterns shared by 2 patients. The clustering value of 53 percent was much higher in comparison with clustering of M. tuberculosis strains from TB patients without clinical signs for HIV infection from randomly selected health centers. We present these results as preliminary data on M. tuberculosis strain polymorphism in Brazil and on the higher risk for recent transmission amongst patients with AIDS.

The changing face of the epidemiology of tuberculosis due to molecular strain typing--a review.

About one third of the world population is infected with tubercle bacilli, causing eight million new cases of tuberculosis (TB) and three million deaths each year. After years of lack of interest in the disease, World Health Organization recently declared TB a global emergency and it is clear that there is need for more efficient national TB programs and newly defined research priorities. A more complete epidemiology of tuberculosis will lead to a better identification of index cases and to a more efficient treatment of the disease. Recently, new molecular tools became available for the identification of strains of Mycobacterium tuberculosis (M. tuberculosis), allowing a better recognition of transmission routes of defined strains. Both a standardized restriction-fragment-length-polymorphism-based methodology for epidemiological studies on a large scale and deoxyribonucleic acids (DNA) amplification-based methods that allow rapid detection of outbreaks with multidrug-resistant (MDR) strains, often characterized by high mortality rates, have been developed. This review comments on the existing methods of DNA-based recognition of M. tuberculosis strains and their peculiarities. It also summarizes literature data on the application of molecular fingerprinting for detection of outbreaks of M. tuberculosis, for identification of index cases, for study of interaction between TB and infection with the human immuno-deficiency virus, for analysis of the behavior of MDR strains, for a better understanding of risk factors for transmission of TB within communities and for population-based studies of TB transmission within and between countries.

Use of the polymerase chain reaction in the diagnosis of leprosy.

One of the main limitations for successful epidemiological control of leprosy is the lack of a method for its diagnosis in subclinical cases. Because of the long incubation period of the disease, liberation and spread of Mycobacterium leprae during subclinical stages-principally in cases of untreated multibacillary forms of leprosy-constitute the main source of infection. This report describes the use of the polymerase chain reaction (PCR) for the detection of M. leprae in different types of tissue samples (blood, lymph, nasal secretion and hair) from an individual who was suspected of having leprosy. Although no conclusive diagnosis could be made by traditional diagnostic methods, the individual was found to be infected with M. leprae after amplification of the bacterial DNA.

The Trypanosoma cruzi genome project: nuclear karyotype and gene mapping of clone CL Brener.

By using improved pulsed field gel electrophoresis conditions, the molecular karyotype of the reference clone CL Brener selected for Trypanosoma cruzi genome project was established. A total of 20 uniform chromosomal bands ranging in size from 0.45 to 3.5 Megabase pairs (Mbp) were resolved in a single run. The weighted sum of the chromosomal bands was approximately 87 Mbp. Chromoblots were hybridized with 39 different homologous probes, 13 of which identified single chromosomes. Several markers showed linkage and four different linkage groups were identified, each comprising two markers. Densitometric analysis suggests that most of the chromosomal bands contain two or more chromosomes representing either homologous chromosomes and/or heterologous chromosomes with similar sizes.