Whole-genome sequencing and microarray analysis of ex vivo Plasmodium vivax reveal selective pressure on putative drug resistance genes.

Plasmodium vivax causes 25-40% of malaria cases worldwide, yet research on this human malaria parasite has been neglected. Nevertheless, the recent publication of the P. vivax reference genome now allows genomics and systems biology approaches to be applied to this pathogen. We show here that whole-genome analysis of the parasite can be achieved directly from ex vivo-isolated parasites, without the need for in vitro propagation. A single isolate of P. vivax obtained from a febrile patient with clinical malaria from Peru was subjected to whole-genome sequencing (30× coverage). This analysis revealed over 18,261 single-nucleotide polymorphisms (SNPs), 6,257 of which were further validated using a tiling microarray. Within core chromosomal genes we find that one SNP per every 985 bases of coding sequence distinguishes this recent Peruvian isolate, designated IQ07, from the reference Salvador I strain obtained in 1972. This full-genome sequence of an uncultured P. vivax isolate shows that the same regions with low numbers of aligned sequencing reads are also highly variable by genomic microarray analysis. Finally, we show that the genes containing the largest ratio of nonsynonymous-to-synonymous SNPs include two AP2 transcription factors and the P. vivax multidrug resistance-associated protein (PvMRP1), an ABC transporter shown to be associated with quinoline and antifolate tolerance in Plasmodium falciparum. This analysis provides a data set for comparative analysis with important potential for identifying markers for global parasite diversity and drug resistance mapping studies.

Genome-wide nucleosome mapping of Plasmodium falciparum reveals histone-rich coding and histone-poor intergenic regions and chromatin remodeling of core and subtelomeric genes.

BACKGROUND: Epigenetic modifications of histones and regulation of chromatin structure have been implicated in regulation of virulence gene families in P. falciparum. To better understand chromatin-mediated gene regulation, we used a high-density oligonucleotide microarray to map the position and enrichment of nucleosomes across the entire genome of P. falciparum at three time points of the intra-erythrocytic developmental cycle (IDC) in vitro. We used an unmodified histone H4 antibody for chromatin immunoprecipitation of nucleosome-bound DNA. RESULTS: We observed generally low nucleosomal occupancy of intergenic regions and higher occupancy of protein coding regions. In contract to the overall small fluctuation of nucleosomal occupancy in most coding regions throughout the IDC, subtelomeric genes encoding surface proteins such as var and rif, as well as some core chromosomal genes such as transcription factors, showed large changes in chromatin structure. Telomeres harbored a region with the highest nucleosomal occupancy of the genome and also exhibited large changes with higher nucleosomal occupancy at schizont stages. While many of these subtelomeric genes were previously shown to be modified by H3K9 trimethylation, we also identified some housekeeping genes in core chromosome regions that showed extensive changes in chromatin structure but do not contain this modification. tRNA and basal transcription factor genes showed low nucleosomal occupancy at all times, suggesting of an open chromatin structure that might be permissive for constitutively high levels of expression. Generally, nucleosomal occupancy was not correlated with the steady-state mRNA levels. Several var genes were exceptions: the var gene with the highest expression level showed the lowest nucleosomal occupancy, and selection of parasites for var2CSA expression resulted in lower nucleosomal occupancy at the var2CSA locus. We identified nucleosome-free regions in intergenic regions that may serve as transcription start sites or transcription factor binding sites. Using the nucleosomal occupancy data as the baseline, we further mapped the genome-wide enrichment of H3K9 acetylation and detected general enrichment of this mark in intergenic regions. CONCLUSIONS: These data on nucleosome enrichment changes add to our understanding of the influence of chromatin structure on the regulation of gene expression. Histones are generally enriched in coding regions, and relatively poor in intergenic regions. Histone enrichment patterns allow for identification of new putative gene-coding regions. Most genes do not show correlation between chromatin structure and steady-state mRNA levels, indicating the dominant roles of other regulatory mechanisms. We present a genome-wide nucleosomal occupancy map, which can be used as a reference for future experiments of histone modification mapping.

A population study of the minicircles in Trypanosoma cruzi: predicting guide RNAs in the absence of empirical RNA editing.

BACKGROUND: The structurally complex network of minicircles and maxicircles comprising the mitochondrial DNA of kinetoplastids mirrors the complexity of the RNA editing process that is required for faithful expression of encrypted maxicircle genes. Although a few of the guide RNAs that direct this editing process have been discovered on maxicircles, guide RNAs are mostly found on the minicircles. The nuclear and maxicircle genomes have been sequenced and assembled for Trypanosoma cruzi, the causative agent of Chagas disease, however the complement of 1.4-kb minicircles, carrying four guide RNA genes per molecule in this parasite, has been less thoroughly characterised. RESULTS: Fifty-four CL Brener and 53 Esmeraldo strain minicircle sequence reads were extracted from T. cruzi whole genome shotgun sequencing data. With these sequences and all published T. cruzi minicircle sequences, 108 unique guide RNAs from all known T. cruzi minicircle sequences and two guide RNAs from the CL Brener maxicircle were predicted using a local alignment algorithm and mapped onto predicted or experimentally determined sequences of edited maxicircle open reading frames. For half of the sequences no statistically significant guide RNA could be assigned. Likely positions of these unidentified gRNAs in T. cruzi minicircle sequences are estimated using a simple Hidden Markov Model. With the local alignment predictions as a standard, the HMM had an ~85% chance of correctly identifying at least 20 nucleotides of guide RNA from a given minicircle sequence. Inter-minicircle recombination was documented. Variable regions contain species-specific areas of distinct nucleotide preference. Two maxicircle guide RNA genes were found. CONCLUSION: The identification of new minicircle sequences and the further characterization of all published minicircles are presented, including the first observation of recombination between minicircles. Extrapolation suggests a level of 4% recombinants in the population, supporting a relatively high recombination rate that may serve to minimize the persistence of gRNA pseudogenes. Characteristic nucleotide preferences observed within variable regions provide potential clues regarding the transcription and maturation of T. cruzi guide RNAs. Based on these preferences, a method of predicting T. cruzi guide RNAs using only primary minicircle sequence data was created.

Trypanosoma cruzi mitochondrial maxicircles display species- and strain-specific variation and a conserved element in the non-coding region.

BACKGROUND: The mitochondrial DNA of kinetoplastid flagellates is distinctive in the eukaryotic world due to its massive size, complex form and large sequence content. Comprised of catenated maxicircles that contain rRNA and protein-coding genes and thousands of heterogeneous minicircles encoding small guide RNAs, the kinetoplast network has evolved along with an extreme form of mRNA processing in the form of uridine insertion and deletion RNA editing. Many maxicircle-encoded mRNAs cannot be translated without this post-transcriptional sequence modification. RESULTS: We present the complete sequence and annotation of the Trypanosoma cruzi maxicircles for the CL Brener and Esmeraldo strains. Gene order is syntenic with Trypanosoma brucei and Leishmania tarentolae maxicircles. The non-coding components have strain-specific repetitive regions and a variable region that is unique for each strain with the exception of a conserved sequence element that may serve as an origin of replication, but shows no sequence identity with L. tarentolae or T. brucei. Alternative assemblies of the variable region demonstrate intra-strain heterogeneity of the maxicircle population. The extent of mRNA editing required for particular genes approximates that seen in T. brucei. Extensively edited genes were more divergent among the genera than non-edited and rRNA genes. Esmeraldo contains a unique 236-bp deletion that removes the 5'-ends of ND4 and CR4 and the intergenic region. Esmeraldo shows additional insertions and deletions outside of areas edited in other species in ND5, MURF1, and MURF2, while CL Brener has a distinct insertion in MURF2. CONCLUSION: The CL Brener and Esmeraldo maxicircles represent two of three previously defined maxicircle clades and promise utility as taxonomic markers. Restoration of the disrupted reading frames might be accomplished by strain-specific RNA editing. Elements in the non-coding region may be important for replication, transcription, and anchoring of the maxicircle within the kinetoplast network.

Differential transcription profiles in Trypanosoma cruzi associated with clinical forms of Chagas disease: Maxicircle NADH dehydrogenase subunit 7 gene truncation in asymptomatic patient isolates.

The majority of individuals in the chronic phase of Chagas disease are asymptomatic (indeterminate form). Every year 2-3% of these individuals develop severe clinical manifestations (cardiac and digestive forms). In this study a Trypanosoma cruzi DNA microarray was used to compare the transcript profiles of six human isolates: three from asymptomatic and three from cardiac patients. Seven signals were expressed differentially between the two classes of isolates, including tryparedoxin, surface protease GP63, cyclophilin, some hypothetical proteins and the pre-edited maxicircle gene NADH dehydrogenase subunit 7 (ND7). The approximately 30-fold greater signal in cardiac strains for ND7 was the most pronounced of the group, and differential levels of pre-edited ND7 transcript confirmed the microarray analysis. The ND7 gene from asymptomatic isolates showed a deletion of 455bp from nt 222 to nt 677 relative to ND7 of the CL Brener reference strain. The ND7 gene structure correlated with disease manifestation for 20 isolates from clinically characterised, chronic phase patients. The ND7 lesion produces a truncated product that could impair the function of mitochondrial complex I. Possible links between the integrity of the electron transport chain and symptom presentation are discussed. We propose that ND7 and other genes of the pathway constitute valuable targets for PCR assays in the differential diagnosis of the infective T. cruzi strain. While this hypothesis requires validation by the examination of additional recent parasite isolates from patients with defined pathologies, the identification of specific molecular markers represents a promising advance in the association between parasite genetics and disease pathology.

Two hybridization events define the population structure of Trypanosoma cruzi.

Genetic variation in Trypanosoma cruzi is likely a key determinant in transmission and pathogenesis of Chagas disease. We have examined nine loci as markers for the extant T. cruzi strains. Four distinct alleles were found for each locus, corresponding to the sequence classes present in the homozygous discrete typing units (DTUs) I, IIa, IIb, and IIc. The alleles in DTUs IIa and IIc showed a spectrum of polymorphism ranging from DTU I-like to DTU IIb-like, in addition to DTU-specific sequence variation. DTUs IId and IIe were indistinguishable, showing DTU homozygosity at one locus and heterozygosity with DTU IIb and IIc allelic sequences at eight loci. Recombination between the DTU IIb and IIc alleles is evidenced from mosaic polymorphisms. These data imply that two discrete hybridization events resulted in the formation of the current DTUs. We propose a model in which a fusion between ancestral DTU I and IIb strains gave rise to a heterozygous hybrid that homogenized its genome to become the homozygous progenitor of DTUs IIa and IIc. The second hybridization between DTU IIb and IIc strains that generated DTUs IId and IIe resulted in extensive heterozygosity with subsequent recombination of parental genotypes.

Trypanosoma cruzi 5S rRNA arrays define five groups and indicate the geographic origins of an ancestor of the heterozygous hybrids.

Isolates of the etiological agent of Chagas disease, Trypanosoma cruzi, have been subdivided into six subgroups referred to as discrete typing units. The subgroups are related through two distinct hybridisation events: representatives of homozygous discrete typing units I and IIb fused to form discrete typing units IIa and IIc, whose homozygous genotypes have features of both ancestral types; a second fusion between strains of homozygous discrete typing units IIb and IIc created the heterozygous hybrid strains discrete typing units IId and IIe. The intergenic region of the tandemly repeated 5S rRNA array displays four variant sequence classes, allowing the discrimination of five discrete typing units. The genome project reference strain, CL Brener, is a hybrid discrete typing unit IIe strain that contains both discrete typing unit IIb and IIc classes of 5S rRNA repeats in distinct arrays present on different chromosomes. The CL Brener discrete typing unit IIb-type array contains approximately 193 repeated units, of which about one-third contain a 129 bp sequence that replaces a majority of the 5S rRNA sequence. The 129 bp 'invader' sequence was detected within the arrays of all hybrid discrete typing unit IId and IIe strains and in a subset of discrete typing unit IIb strains. This array invader replaces the internal promoter elements conserved in 5S rRNA. The discrete typing unit IIb Esmeraldo strain contains approximately 135 repeats and shows a region of homology to the array invader in the 5' flank of the array, but no evidence of the invading sequence element within the array. A survey of additional discrete typing unit IIb strains revealed a split within the subgroup, in which some strains contained invaded arrays and others were homogeneous for the 5S rRNA. The putative discrete typing unit IIb ancestor of the hybrid discrete typing units IId and IIe more closely resembles the extant Bolivian/Chilean IIb isolates than the Brazilian IIb isolates based on the correlation with the array invader.

Intragenomic spliced leader RNA array analysis of kinetoplastids reveals unexpected transcribed region diversity in Trypanosoma cruzi.

The spliced leader RNA gene (SL RNA) repeat is present in large multicopy arrays and has been used as a marker for the diversity of kinetoplastid protozoans. Intra-array variation could affect conclusions made using a randomly isolated repeat as a marker. We examined the Leishmania major (Friedlin) and Trypanosoma cruzi (CL Brener) genome projects for SL RNA repeat sequences in order to assess their homogeneity and the possible effects of sequence variation on taxonomic interpretation. Of the dozens of distinct sequence classes examined, no single copy would bias clustering analyses with regard to other closely related species or isolates. Six dimorphic sites within the T. cruzi transcribed region were found to be linked and are predicted to yield a heterogeneous SL RNA population. The variation that exists among the repeats paints a picture of the broad mechanisms of array maintenance and evolution where site-specific mutations in a single repeat may be spread throughout the array and recombined with existing repeats to create new sequence classes, all occurring under selective pressure to maintain or increase the fitness of the cell line in which these events occur.

First-generation linkage map of the gray, short-tailed opossum, Monodelphis domestica, reveals genome-wide reduction in female recombination rates.

The gray, short-tailed opossum, Monodelphis domestica, is the most extensively used, laboratory-bred marsupial resource for basic biologic and biomedical research worldwide. To enhance the research utility of this species, we are building a linkage map, using both anonymous markers and functional gene loci, that will enable the localization of quantitative trait loci (QTL) and provide comparative information regarding the evolution of mammalian and other vertebrate genomes. The current map is composed of 83 loci distributed among eight autosomal linkage groups and the X chromosome. The autosomal linkage groups appear to encompass a very large portion of the genome, yet span a sex-average distance of only 633.0 cM, making this the most compact linkage map known among vertebrates. Most surprising, the male map is much larger than the female map (884.6 cM vs. 443.1 cM), a pattern contrary to that in eutherian mammals and other vertebrates. The finding of genome-wide reduction in female recombination in M. domestica, coupled with recombination data from two other, distantly related marsupial species, suggests that reduced female recombination might be a widespread metatherian attribute. We discuss possible explanations for reduced female recombination in marsupials as a consequence of the metatherian characteristic of determinate paternal X chromosome inactivation.

The determinants of Chagas disease: connecting parasite and host genetics.

As a consequence of infection by Trypanosoma cruzi, 30% of victims may develop chronic Chagas disease, which presents a spectrum of pathology including cardiomyopathy, megacolon and megaesophagus. The outcome of infection in a particular individual is the result of a set of complex interactions among the host genetic background, environmental and social factors, and the genetic composition of the parasite, all of which can be complicated by mixed infections and re-infections. Initially we consider what is known about the genetic structure and biological properties of the protozoan. Currently, six distinct subgroups have been characterized by different combinations of four distinct genotypic classes. The recent demonstration of genetic exchange via non-meiotic cell fusion illustrates a mechanism by which maintained heterogeneous polyploidy may have been generated in these parasites. Subsequently, we consider factors in humans and in experimental mouse-infection and tissue culture models that have contributed to our understanding of the host's susceptibility or resistance to disease. Identification of the direct players in host-pathogen interactions at the establishment and chronic phases of the disease is perhaps the best hope of a clinical handle for treatment. At some point in the future, these disparate areas of study will have to come together. It is to be hoped that this scientific fusion will result in better prognosis and treatment of Chagas disease.

Evidence for multiple hybrid groups in Trypanosoma cruzi.

A role for parasite genetic variability in the spectrum of Chagas disease is emerging but not yet evident, in part due to an incomplete understanding of the population structure of Trypanosoma cruzi. To investigate further the observed genotypic variation at the sequence and chromosomal levels in strains of standard and field-isolated T. cruzi we have undertaken a comparative analysis of 10 regions of the genome from two isolates representing T. cruzi I (Dm28c and Silvio X10) and two from T. cruzi II (CL Brener and Esmeraldo). Amplified regions contained intergenic (non-coding) sequences from tandemly repeated genes. Multiple nucleotide polymorphisms correlated with the T. cruzi I/T. cruzi II classification. Two intergenic regions had useful polymorphisms for the design of classification probes to test on genomic DNA from other known isolates. Two adjacent nucleotide polymorphisms in HSP 60 correlated with the T. cruzi I and T. cruzi II distinction. 1F8 nucleotide polymorphisms revealed multiple subdivisions of T. cruzi II: subgroups IIa and IIc displayed the T. cruzi I pattern; subgroups IId and IIe possessed both the I and II patterns. Furthermore, isolates from subgroups IId and IIe contained the 1F8 polymorphic markers on different chromosome bands supporting a genetic exchange event that resulted in chromosomes V and IX of T. cruzi strain CL Brener. Based on these analyses, T. cruzi I and subgroup IIb appear to be pure lines, while subgroups IIa/IIc and IId/IIe are hybrid lines. These data demonstrate for the first time that IIa/IIc are hybrid, consistent with the hypothesis that genetic recombination has occurred more than once within the T. cruzi lines.

A systems-based analysis of Plasmodium vivax lifecycle transcription from human to mosquito.

BACKGROUND: Up to 40% of the world's population is at risk for Plasmodium vivax malaria, a disease that imposes a major public health and economic burden on endemic countries. Because P. vivax produces latent liver forms, eradication of P. vivax malaria is more challenging than it is for P. falciparum. Genetic analysis of P. vivax is exceptionally difficult due to limitations of in vitro culture. To overcome the barriers to traditional molecular biology in P. vivax, we examined parasite transcriptional changes in samples from infected patients and mosquitoes in order to characterize gene function, define regulatory sequences and reveal new potential vaccine candidate genes. PRINCIPAL FINDINGS: We observed dramatic changes in transcript levels for various genes at different lifecycle stages, indicating that development is partially regulated through modulation of mRNA levels. Our data show that genes involved in common biological processes or molecular machinery are co-expressed. We identified DNA sequence motifs upstream of co-expressed genes that are conserved across Plasmodium species that are likely binding sites of proteins that regulate stage-specific transcription. Despite their capacity to form hypnozoites we found that P. vivax sporozoites show stage-specific expression of the same genes needed for hepatocyte invasion and liver stage development in other Plasmodium species. We show that many of the predicted exported proteins and members of multigene families show highly coordinated transcription as well. CONCLUSIONS: We conclude that high-quality gene expression data can be readily obtained directly from patient samples and that many of the same uncharacterized genes that are upregulated in different P. vivax lifecycle stages are also upregulated in similar stages in other Plasmodium species. We also provide numerous examples of how systems biology is a powerful method for determining the likely function of genes in pathogens that are neglected due to experimental intractability.

Use of high-density tiling microarrays to identify mutations globally and elucidate mechanisms of drug resistance in Plasmodium falciparum.

BACKGROUND: The identification of genetic changes that confer drug resistance or other phenotypic changes in pathogens can help optimize treatment strategies, support the development of new therapeutic agents, and provide information about the likely function of genes. Elucidating mechanisms of phenotypic drug resistance can also assist in identifying the mode of action of uncharacterized but potent antimalarial compounds identified in high-throughput chemical screening campaigns against Plasmodium falciparum. RESULTS: Here we show that tiling microarrays can detect de novo a large proportion of the genetic changes that differentiate one genome from another. We show that we detect most single nucleotide polymorphisms or small insertion deletion events and all known copy number variations that distinguish three laboratory isolates using readily accessible methods. We used the approach to discover mutations that occur during the selection process after transfection. We also elucidated a mechanism by which parasites acquire resistance to the antimalarial fosmidomycin, which targets the parasite isoprenoid synthesis pathway. Our microarray-based approach allowed us to attribute in vitro derived fosmidomycin resistance to a copy number variation event in the pfdxr gene, which enables the parasite to overcome fosmidomycin-mediated inhibition of isoprenoid biosynthesis. CONCLUSIONS: We show that newly emerged single nucleotide polymorphisms can readily be detected and that malaria parasites can rapidly acquire gene amplifications in response to in vitro drug pressure. The ability to define comprehensively genetic variability in P. falciparum with a single overnight hybridization creates new opportunities to study parasite evolution and improve the treatment and control of malaria.

Discovery and barcoding by analysis of spliced leader RNA gene sequences of new isolates of Trypanosomatidae from Heteroptera in Costa Rica and Ecuador.

Trypanosomatid diversity in Heteroptera was sampled using a culture-independent approach based on amplification and sequencing of Spliced Leader RNA gene repeats from environmental samples. By combining the data collected herein with that of previous work, the prevalence of parasites was found to be 22%-23%. Out of approximately 170 host species investigated nearly 60 were found to harbor trypanosomatids. The parasites found were grouped by cluster analysis into 48 typing units. Most of these were well separated from the known groups and, therefore, likely represent new trypanosomatid species. The sequences for each typing unit serve as barcodes to facilitate their recognition in the future. As the sampled host species represent a minor fraction of potential hosts, the entire trypanosomatid diversity is far greater than described thus far. Investigations of trypanosomatid diversity, host-specificity, and biogeography have become feasible using the approach described herein.

Development, inheritance, and linkage-group assignment of 60 novel microsatellite markers for the gray, short-tailed opossum Monodelphis domestica.

Short-tandem-repeat (SSR) or microsatellite polymorphisms are some of the most extensively employed genetic markers in contemporary linkage mapping studies. To date, only a limited number of microsatellites have been isolated in the gray, short-tailed opossum Monodelphis domestica, a South American marsupial widely used for comparative biological and biomedical research. To increase the number of potentially useful mapping markers, we screened 2 microsatellite-enriched genomic libraries containing alternatively (CA)n or (GA)n repeats. A total of 184 clones were sequenced, from which 60 polymorphic microsatellite markers were successfully optimized. The efficiency of this enrichment protocol for M. domestica microsatellite isolation is discussed, and suggestions to improve the outcome are made. All 60 loci showed high allelic diversity, with allele numbers ranging from 2 to 10 in a subset of 33 unrelated animals. Normal Mendelian inheritance was confirmed for all loci by analyzing allelic segregation in 5 two-generation families. One microsatellite appeared to be X linked, and null alleles were found in 5 others. Two-point linkage analyses were implemented using the data on the 5 families, leading to the assignment of 59 of these loci to the existing linkage groups. The 60 novel microsatellites developed in this study will contribute significantly to the M. domestica linkage map, and further QTL mapping studies.

The genome sequence of Trypanosoma cruzi, etiologic agent of Chagas disease.

Whole-genome sequencing of the protozoan pathogen Trypanosoma cruzi revealed that the diploid genome contains a predicted 22,570 proteins encoded by genes, of which 12,570 represent allelic pairs. Over 50% of the genome consists of repeated sequences, such as retrotransposons and genes for large families of surface molecules, which include trans-sialidases, mucins, gp63s, and a large novel family (>1300 copies) of mucin-associated surface protein (MASP) genes. Analyses of the T. cruzi, T. brucei, and Leishmania major (Tritryp) genomes imply differences from other eukaryotes in DNA repair and initiation of replication and reflect their unusual mitochondrial DNA. Although the Tritryp lack several classes of signaling molecules, their kinomes contain a large and diverse set of protein kinases and phosphatases; their size and diversity imply previously unknown interactions and regulatory processes, which may be targets for intervention.

Diversity of insect trypanosomatids assessed from the spliced leader RNA and 5S rRNA genes and intergenic regions.

We have determined the sequences of 5S rRNA and spliced leader (SL) RNA genes, and adjacent intergenic regions for representatives of all known trypanosomatid genera parasitizing insects. The genetic loci have been analyzed separately as well as by a combined approach. Several isolates, assigned by morphology to different genera (Leptomonas spp., Blastocrithidia spp.), seem to belong to a single species with an unexpectedly wide host and geographical range. An unnamed trypanosomatid isolated from rats in Egypt was found to belong to the genus Herpetomonas, so far associated with insect hosts only. It is closely related to Herpetomonas ztiplika, a parasite of a blood-sucking biting midge. Apparently several different trypanosomatid species can infect one insect species, as exemplified by Leptomonas sp. PL and Wallaceina sp. Wsd, which were isolated from different specimens of Salda littoralis on the same locality and day. However, since the same species of Leptomonas was obtained from insect hosts belonging to different genera, some insect trypanosomatids may have low host specificity. Our data revealed additional discrepancies between molecular phylogenetic data and cell morphology, rendering current trypanosomatid taxonomy unreliable.