Population genomics reveals the expansion of highly inbred Plasmodium vivax lineages in the main malaria hotspot of Brazil.

BACKGROUND: Plasmodium vivax is a neglected human malaria parasite that causes significant morbidity in the Americas, the Middle East, Asia, and the Western Pacific. Population genomic approaches remain little explored to map local and regional transmission pathways of P. vivax across the main endemic sites in the Americas, where great progress has been made towards malaria elimination over the past decades. METHODOLOGY/PRINCIPAL FINDINGS: We analyze 38 patient-derived P. vivax genome sequences from Mâncio Lima (ML)-the Amazonian malaria hotspot next to the Brazil-Peru border-and 24 sequences from two other sites in Acre State, Brazil, a country that contributes 23% of malaria cases in the Americas. We show that the P. vivax population of ML is genetically diverse (π = 4.7 × 10-4), with a high polymorphism particularly in genes encoding proteins putatively involved in red blood cell invasion. Paradoxically, however, parasites display strong genome-wide linkage disequilibrium, being fragmented into discrete lineages that are remarkably stable across time and space, with only occasional recombination between them. Using identity-by-descent approaches, we identified a large cluster of closely related sequences that comprises 16 of 38 genomes sampled in ML over 26 months. Importantly, we found significant ancestry sharing between parasites at a large geographic distance, consistent with substantial gene flow between regional P. vivax populations. CONCLUSIONS/SIGNIFICANCE: We have characterized the sustained expansion of highly inbred P. vivax lineages in a malaria hotspot that can seed regional transmission. Potential source populations in hotspots represent a priority target for malaria elimination in the Amazon.

Genome-wide identification of evolutionarily conserved Small Heat-Shock and eight other proteins bearing α-crystallin domain-like in kinetoplastid protists.

Small Heat-Shock Proteins (sHSPs) and other proteins bearing alpha-crystallin domains (ACD) participate in defense against heat and oxidative stress and play important roles in cell cycle, cytoskeleton dynamics, and immunological and pathological mechanisms in eukaryotes. However, little is known about these proteins in early-diverging lineages of protists such as the kinetoplastids. Here, ACD-like proteins (ACDp) were investigated in genomes of 61 species of 12 kinetoplastid genera, including Trypanosoma spp. (23 species of mammals, reptiles and frogs), Leishmania spp. (mammals and lizards), trypanosomatids of insects, Phytomonas spp. of plants, and bodonids. Comparison of ACDps based on domain architecture, predicted tertiary structure, phylogeny and genome organization reveals a kinetoplastid evolutionarily conserved repertoire, which diversified prior to trypanosomatid adaptation to parasitic life. We identified 9 ACDp orthologs classified in 8 families of TryACD: four previously recognized (HSP20, Tryp23A, Tryp23B and ATOM69), and four characterized for the first time in kinetoplastids (TryACDP, TrySGT1, TryDYX1C1 and TryNudC). A single copy of each ortholog was identified in each genome alongside TryNudC1/TrypNudC2 homologs and, overall, ACDPs were under strong selection pressures at main phylogenetic lineages. Transcripts of all ACDPs were identified across the life stages of T. cruzi, T. brucei and Leishmania spp., but proteomic profiles suggested that most ACDPs may be species- and stage-regulated. Our findings establish the basis for functional studies, and provided evolutionary and structural support for an underestimated repertoire of ACDps in the kinetoplastids.

Human migration and the spread of malaria parasites to the New World.

We examined the mitogenomes of a large global collection of human malaria parasites to explore how and when Plasmodium falciparum and P. vivax entered the Americas. We found evidence of a significant contribution of African and South Asian lineages to present-day New World malaria parasites with additional P. vivax lineages appearing to originate from Melanesia that were putatively carried by the Australasian peoples who contributed genes to Native Americans. Importantly, mitochondrial lineages of the P. vivax-like species P. simium are shared by platyrrhine monkeys and humans in the Atlantic Forest ecosystem, but not across the Amazon, which most likely resulted from one or a few recent human-to-monkey transfers. While enslaved Africans were likely the main carriers of P. falciparum mitochondrial lineages into the Americas after the conquest, additional parasites carried by Australasian peoples in pre-Columbian times may have contributed to the extensive diversity of extant local populations of P. vivax.

Expanded repertoire of kinetoplast associated proteins and unique mitochondrial DNA arrangement of symbiont-bearing trypanosomatids.

In trypanosomatids, the kinetoplast is the portion of the single mitochondrion that is connected to the basal body and contains the kDNA, a network composed by circular and interlocked DNA. The kDNA packing is conducted by Kinetoplast Associated Proteins (KAPs), which are similar to eukaryotic histone H1. In symbiont-harboring trypanosomatids (SHTs) such as Angomonas deanei and Strigomonas culicis, a ß-proteobacterium co-evolves with the host in a mutualistic relationship. The prokaryote confers nutritional benefits to the host and affects its cell structure. Atomic force microscopy showed that the topology of isolated kDNA networks is quite similar in the two SHT species. Ultrastructural analysis using high-resolution microscopy techniques revealed that the DNA fibrils are more compact in the kinetoplast region that faces the basal body and that the presence of the symbiotic bacterium does not interfere with kDNA topology. However, RT-PCR data revealed differences in the expression of KAPs in wild-type protozoa as compared to aposymbiotic cells. Immunolocalization showed that different KAPs present distinct distributions that are coincident in symbiont-bearing and in symbiont-free cells. Although KAP4 and KAP7 are shared by all trypanosomatid species, the expanded repertoire of KAPs in SHTs can be used as phylogenetic markers to distinguish different genera.

Genome-wide diversity and differentiation in New World populations of the human malaria parasite Plasmodium vivax.

BACKGROUND: The Americas were the last continent colonized by humans carrying malaria parasites. Plasmodium falciparum from the New World shows very little genetic diversity and greater linkage disequilibrium, compared with its African counterparts, and is clearly subdivided into local, highly divergent populations. However, limited available data have revealed extensive genetic diversity in American populations of another major human malaria parasite, P. vivax. METHODS: We used an improved sample preparation strategy and next-generation sequencing to characterize 9 high-quality P. vivax genome sequences from northwestern Brazil. These new data were compared with publicly available sequences from recently sampled clinical P. vivax isolates from Brazil (BRA, total n = 11 sequences), Peru (PER, n = 23), Colombia (COL, n = 31), and Mexico (MEX, n = 19). PRINCIPAL FINDINGS/CONCLUSIONS: We found that New World populations of P. vivax are as diverse (nucleotide diversity π between 5.2 × 10-4 and 6.2 × 10-4) as P. vivax populations from Southeast Asia, where malaria transmission is substantially more intense. They display several non-synonymous nucleotide substitutions (some of them previously undescribed) in genes known or suspected to be involved in antimalarial drug resistance, such as dhfr, dhps, mdr1, mrp1, and mrp-2, but not in the chloroquine resistance transporter ortholog (crt-o) gene. Moreover, P. vivax in the Americas is much less geographically substructured than local P. falciparum populations, with relatively little between-population genome-wide differentiation (pairwise FST values ranging between 0.025 and 0.092). Finally, P. vivax populations show a rapid decline in linkage disequilibrium with increasing distance between pairs of polymorphic sites, consistent with very frequent outcrossing. We hypothesize that the high diversity of present-day P. vivax lineages in the Americas originated from successive migratory waves and subsequent admixture between parasite lineages from geographically diverse sites. Further genome-wide analyses are required to test the demographic scenario suggested by our data.

Using mitochondrial genome sequences to track the origin of imported Plasmodium vivax infections diagnosed in the United States.

Although the geographic origin of malaria cases imported into the United States can often be inferred from travel histories, these histories may be lacking or incomplete. We hypothesized that mitochondrial haplotypes could provide region-specific molecular barcodes for tracing the origin of imported Plasmodium vivax infections. An analysis of 348 mitochondrial genomes from worldwide parasites and new sequences from 69 imported malaria cases diagnosed across the United States allowed for a geographic assignment of most infections originating from the Americas, southeast Asia, east Asia, and Melanesia. However, mitochondrial lineages from Africa, south Asia, central Asia, and the Middle East, which altogether contribute the vast majority of imported malaria cases in the United States, were closely related to each other and could not be reliably assigned to their geographic origins. More mitochondrial genomes are required to characterize molecular barcodes of P. vivax from these regions.

Repertoire, genealogy and genomic organization of cruzipain and homologous genes in Trypanosoma cruzi, T. cruzi-like and other trypanosome species.

Trypanosoma cruzi, the agent of Chagas disease, is a complex of genetically diverse isolates highly phylogenetically related to T. cruzi-like species, Trypanosoma cruzi marinkellei and Trypanosoma dionisii, all sharing morphology of blood and culture forms and development within cells. However, they differ in hosts, vectors and pathogenicity: T. cruzi is a human pathogen infective to virtually all mammals whilst the other two species are non-pathogenic and bat restricted. Previous studies suggest that variations in expression levels and genetic diversity of cruzipain, the major isoform of cathepsin L-like (CATL) enzymes of T. cruzi, correlate with levels of cellular invasion, differentiation, virulence and pathogenicity of distinct strains. In this study, we compared 80 sequences of genes encoding cruzipain from 25 T. cruzi isolates representative of all discrete typing units (DTUs TcI-TcVI) and the new genotype Tcbat and 10 sequences of homologous genes from other species. The catalytic domain repertoires diverged according to DTUs and trypanosome species. Relatively homogeneous sequences are found within and among isolates of the same DTU except TcV and TcVI, which displayed sequences unique or identical to those of TcII and TcIII, supporting their origin from the hybridization between these two DTUs. In network genealogies, sequences from T. cruzi clustered tightly together and closer to T. c. marinkellei than to T. dionisii and largely differed from homologues of T. rangeli and T. b. brucei. Here, analysis of isolates representative of the overall biological and genetic diversity of T. cruzi and closest T. cruzi-like species evidenced DTU- and species-specific polymorphisms corroborating phylogenetic relationships inferred with other genes. Comparison of both phylogenetically close and distant trypanosomes is valuable to understand host-parasite interactions, virulence and pathogenicity. Our findings corroborate cruzipain as valuable target for drugs, vaccine, diagnostic and genotyping approaches.