Phylogenomic reconstruction of Cryptosporidium spp. captured directly from clinical samples reveals extensive genetic diversity.

Cryptosporidium is a leading cause of severe diarrhea and mortality in young children and infants in Africa and southern Asia. More than twenty Cryptosporidium species infect humans, of which C. parvum and C. hominis are the major agents causing moderate to severe diarrhea. Relatively few genetic markers are typically applied to genotype and/or diagnose Cryptosporidium. Most infections produce limited oocysts making it difficult to perform whole genome sequencing (WGS) directly from stool samples. Hence, there is an immediate need to apply WGS strategies to 1) develop high-resolution genetic markers to genotype these parasites more precisely, 2) to investigate endemic regions and detect the prevalence of different genotypes, and the role of mixed infections in generating genetic diversity, and 3) to investigate zoonotic transmission and evolution. To understand Cryptosporidium global population genetic structure, we applied Capture Enrichment Sequencing (CES-Seq) using 74,973 RNA-based 120 nucleotide baits that cover ~92% of the genome of C. parvum. CES-Seq is sensitive and successfully sequenced Cryptosporidium genomic DNA diluted up to 0.005% in human stool DNA. It also resolved mixed strain infections and captured new species of Cryptosporidium directly from clinical/field samples to promote genome-wide phylogenomic analyses and prospective GWAS studies.

The short interspersed repetitive element of Trypanosoma cruzi, SIRE, is part of VIPER, an unusual retroelement related to long terminal repeat retrotransposons.

The short interspersed repetitive element (SIRE) of Trypanosoma cruzi was first detected when comparing the sequences of loci that encode the TcP2beta genes. It is present in about 1,500-3,000 copies per genome, depending on the strain, and it is distributed in all chromosomes. An initial analysis of SIRE sequences from 21 genomic fragments allowed us to derive a consensus nucleotide sequence and structure for the element, consisting of three regions (I, II, and III) each harboring distinctive features. Analysis of 158 transcribed SIREs demonstrates that the consensus is highly conserved. The sequences of 51 cDNAs show that SIRE is included in the 3' end of several mRNAs, always transcribed from the sense strand, contributing the polyadenylation site in 63% of the cases. This study led to the characterization of VIPER (vestigial interposed retroelement), a 2,326-bp-long unusual retroelement. VIPER's 5' end is formed by the first 182 bp of SIRE, whereas its 3' end is formed by the last 220 bp of the element. Both SIRE moieties are connected by a 1,924-bp-long fragment that carries a unique ORF encoding a complete reverse transcriptase-RNase H gene whose 15 C-terminal amino acids derive from codons specified by SIRE's region II. The amino acid sequence of VIPER's reverse transcriptase-RNase H shares significant homology to that of long terminal repeat retrotransposons. The fact that SIRE and VIPER sequences are found only in the T. cruzi genome may be of relevance for studies concerning the evolution and the genome flexibility of this protozoan parasite.

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.

Analysis of the distribution of SIRE in the nuclear genome of Trypanosoma cruzi.

The short interspersed repetitive element (SIRE) of the nuclear genome of Trypanosoma cruzi was first detected when comparing the sequences of loci that encode the TcP2beta genes. The present study was designed to assess its distribution and organization in the nuclear genome of the parasite. Southern blots of genomic DNA from different strains demonstrated that each one possesses a defined and characteristic pattern of SIRE distribution. The conservation of the SIRE sequence in T. cruzi strains allowed the development of a rapid inter-SIRE PCR reaction that yields strain-specific amplicon profiles. In the T. cruzi CL Brener clone, we found 1500 copies of the element distributed in all chromosomes. 16 genomic fragments containing SIRE (SZs) were isolated and characterized. In fragments SZ10, SZ12 and SZ31, SIRE was linked to TcRel, a novel repeated sequence that constitutes the 3' end of vp85 genes. SIRE was also linked to an unknown open reading frame in fragments SZ14 and SZ23 which might be related to the subtelomeric regions of T. cruzi chromosomes. Further sequencing of SZ fragments revealed that SIRE was also linked to protein coding genes that have not yet been described in kinetoplastids such as the one coding for PRP22 helicase and a thimet oligopeptidase. To allow the rapid-generation genetic markers associated with SIRE, we developed a SIRE-bubble PCR reaction that provided several such markers for the construction of the physical map of chromosome XVI. The results herein demonstrate that SIRE-associated sites (SAS) may be of great help in physical mapping and interpretation of T. cruzi genomic sequence data.

Towards the physical map of the Trypanosoma cruzi nuclear genome: construction of YAC and BAC libraries of the reference clone T. cruzi CL-Brener

Strategies to construct the physical map of the Trypanosoma cruzi nuclear genome have to capitalize on three main advantage of the parasite genome, namely (a) its small size, (b) the fact that all chromosomes can be defined, and many of them can be isolated by pulse field gel electrophoresis, and (c) the fact that simple Southern blots of electrophoretic karyotypes can be used to map sequence tagged sites and expressed sequence tags to chromosomal bands. A major drawback to cope with is the complexity of T. cruzi genetics, that hinders the construction of a comprehensive genetic map. As a first step towards physical mapping, we report the construction and partial characterization of a T. cruzi CL-Brener genomic library in yeast artificial chromosomes (YACs) that consists of 2.770 individual YACs with a mean insert size of 365 kb encompassing around 10 genomic equivalents. Two libraries in bacterial artificial chromosomes (BACs) have been constructed, BACI and BACII. Both libraries represent about three genome equivalents. A third BAC library (BAC III) is being constructed. YACs and BACs are invaluable tools for physical mapping. More generally, they have to be considered as a common resource for research in Chagas disease.

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.

Towards the physical map of the Trypanosoma cruzi nuclear genome: construction of YAC and BAC libraries of the reference clone T. cruzi CL-Brener.

Strategies to construct the physical map of the Trypanosoma cruzi nuclear genome have to capitalize on three main advantages of the parasite genome, namely (a) its small size, (b) the fact that all chromosomes can be defined, and many of them can be isolated by pulse field gel electrophoresis, and (c) the fact that simple Southern blots of electrophoretic karyotypes can be used to map sequence tagged sites and expressed sequence tags to chromosomal bands. A major drawback to cope with is the complexity of T. cruzi genetics, that hinders the construction of a comprehensive genetic map. As a first step towards physical mapping, we report the construction and partial characterization of a T. cruzi CL-Brener genomic library in yeast artificial chromosomes (YACs) that consists of 2,770 individual YACs with a mean insert size of 365 kb encompassing around 10 genomic equivalents. Two libraries in bacterial artificial chromosomes (BACs) have been constructed, BACI and BACII. Both libraries represent about three genome equivalents. A third BAC library (BAC III) is being constructed. YACs and BACs are invaluable tools for physical mapping. More generally, they have to be considered as a common resource for research in Chagas disease.

Detection of polymorphism in the Trypanosoma cruzi TcP2 beta gene family by single strand conformational analysis (SSCA).

Single strand conformation analysis (SSCA) is a technique that has been used to detect point mutations. We explored its usefulness in the analysis of four different members of the Trypanosoma cruzi TcP2 beta gene family and its suitability for detection of polymorphism in different parasite strains. The availability of primers covering a 97-bp sequence at the 5' end of the genes allowed assessment of the effect of a single base substitution, while the analysis of a 321 bp long sequence permitted the evaluation of sequences differing in several bases. PCR products were analysed under four different electrophoretic conditions: with or without the addition of 10% glycerol in a 6% polyacrylamide gel run at room temperature or at 4 degrees C. Shifts in mobility were radically dependent on the migration condition. Both 97-bp and 321-bp amplicons were best resolved at 4 degrees C, without glycerol. Amplification products derived from total genomic DNA showed a pattern that resembled closely a combination of the products derived from the cloned genes. The results herein demonstrate the usefulness of SSCA to differentiate forms of a complex protozoan gene family, and to scan its polymorphic nature. Furthermore, due to the remarkable sensitivity of the technique it can generate genomic markers, such as Sequence Tagged Sites (STS), of great need in the T. cruzi genome project.

Cloning and sequence analysis of the TcP2 beta cDNA variants of Trypanosoma cruzi.

P2 beta acidic ribosomal protein is a relevant antigen in Chagas' disease. cDNA cloning demonstrated that Trypanosoma cruzi expresses at least six types of TcP2 beta transcripts that differed by point nucleotides substitutions. The distribution and type of mutations seemed to follow the typical structural organization of eukaryotic acidic ribosomal proteins. Most of the synonymous changes clustered in the amino-terminus, suggesting that conservation of this domain was crucial for an equivalent functional ability of each TcP2 beta variant to bind to the ribosome. Interestingly, most amino acid changes in the central-globular and hinge domains caused polymorphism in putative phosphorylatable sites.