Transcriptional heterogeneity and tightly regulated changes in gene expression during Plasmodium berghei sporozoite development.

Despite the critical role of Plasmodium sporozoites in malaria transmission, we still know little about the mechanisms underlying their development in mosquitoes. Here, we use single-cell RNA sequencing to characterize the gene expression profiles of 16,038 Plasmodium berghei sporozoites isolated throughout their development from midgut oocysts to salivary glands, and from forced salivation experiments. Our results reveal a succession of tightly regulated changes in gene expression occurring during the maturation of sporozoites and highlight candidate genes that could play important roles in oocyst egress, sporozoite motility, and the mechanisms underlying the invasion of mosquito salivary glands and mammalian hepatocytes. In addition, the single-cell data reveal extensive transcriptional heterogeneity among parasites isolated from the same anatomical site, suggesting that Plasmodium development in mosquitoes is asynchronous and regulated by intrinsic as well as environmental factors. Finally, our analyses show a decrease in transcriptional activity preceding the translational repression observed in mature sporozoites and associated with their quiescent state in salivary glands, followed by a rapid reactivation of the transcriptional machinery immediately upon salivation.

Whole-genome analysis of Malawian Plasmodium falciparum isolates identifies possible targets of allele-specific immunity to clinical malaria.

Individuals acquire immunity to clinical malaria after repeated Plasmodium falciparum infections. Immunity to disease is thought to reflect the acquisition of a repertoire of responses to multiple alleles in diverse parasite antigens. In previous studies, we identified polymorphic sites within individual antigens that are associated with parasite immune evasion by examining antigen allele dynamics in individuals followed longitudinally. Here we expand this approach by analyzing genome-wide polymorphisms using whole genome sequence data from 140 parasite isolates representing malaria cases from a longitudinal study in Malawi and identify 25 genes that encode possible targets of naturally acquired immunity that should be validated immunologically and further characterized for their potential as vaccine candidates.

Single-cell transcription analysis of Plasmodium vivax blood-stage parasites identifies stage- and species-specific profiles of expression.

Plasmodium vivax and P. falciparum, the parasites responsible for most human malaria worldwide, exhibit striking biological differences, which have important clinical consequences. Unfortunately, P. vivax, unlike P. falciparum, cannot be cultivated continuously in vitro, which limits our understanding of its biology and, consequently, our ability to effectively control vivax malaria. Here, we describe single-cell gene expression profiles of 9,215 P. vivax parasites from bloodstream infections of Aotus and Saimiri monkeys. Our results show that transcription of most P. vivax genes occurs during short periods of the intraerythrocytic cycle and that this pattern of gene expression is conserved in other Plasmodium species. However, we also identify a strikingly high proportion of species-specific transcripts in late schizonts, possibly associated with the specificity of erythrocyte invasion. Our findings provide new and robust markers of blood-stage parasites, including some that are specific to the elusive P. vivax male gametocytes, and will be useful for analyzing gene expression data from laboratory and field samples.

STRIDE: a command-line HMM-based identifier and sub-classifier of Plasmodium falciparum RIFIN and STEVOR variant surface antigen families.

BACKGROUND: RIFINs and STEVORs are variant surface antigens expressed by P. falciparum that play roles in severe malaria pathogenesis and immune evasion. These two highly diverse multigene families feature multiple paralogs, making their classification challenging using traditional bioinformatic methods. RESULTS: STRIDE (STevor and RIfin iDEntifier) is an HMM-based, command-line program that automates the identification and classification of RIFIN and STEVOR protein sequences in the malaria parasite Plasmodium falciparum. STRIDE is more sensitive in detecting RIFINs and STEVORs than available PFAM and TIGRFAM tools and reports RIFIN subtypes and the number of sequences with a FHEYDER amino acid motif, which has been associated with severe malaria pathogenesis. CONCLUSIONS: STRIDE will be beneficial to malaria research groups analyzing genome sequences and transcripts of clinical field isolates, providing insight into parasite biology and virulence.

Human Migration and the Spread of the Nematode Parasite Wuchereria bancrofti.

The human disease lymphatic filariasis causes the debilitating effects of elephantiasis and hydrocele. Lymphatic filariasis currently affects the lives of 90 million people in 52 countries. There are three nematodes that cause lymphatic filariasis, Brugia malayi, Brugia timori, and Wuchereria bancrofti, but 90% of all cases of lymphatic filariasis are caused solely by W. bancrofti (Wb). Here we use population genomics to reconstruct the probable route and timing of migration of Wb strains that currently infect Africa, Haiti, and Papua New Guinea (PNG). We used selective whole genome amplification to sequence 42 whole genomes of single Wb worms from populations in Haiti, Mali, Kenya, and PNG. Our results are consistent with a hypothesis of an Island Southeast Asia or East Asian origin of Wb. Our demographic models support divergence times that correlate with the migration of human populations. We hypothesize that PNG was infected at two separate times, first by the Melanesians and later by the migrating Austronesians. The migrating Austronesians also likely introduced Wb to Madagascar where later migrations spread it to continental Africa. From Africa, Wb spread to the New World during the transatlantic slave trade. Genome scans identified 17 genes that were highly differentiated among Wb populations. Among these are genes associated with human immune suppression, insecticide sensitivity, and proposed drug targets. Identifying the distribution of genetic diversity in Wb populations and selection forces acting on the genome will build a foundation to test future hypotheses and help predict response to current eradication efforts.

Optimization of parasite DNA enrichment approaches to generate whole genome sequencing data for Plasmodium falciparum from low parasitaemia samples.

BACKGROUND: Owing to the large amount of host DNA in clinical samples, generation of high-quality Plasmodium falciparum whole genome sequencing (WGS) data requires enrichment for parasite DNA. Enrichment is often achieved by leukocyte depletion of infected blood prior to storage. However, leukocyte depletion is difficult in low-resource settings and limits analysis to prospectively-collected samples. As a result, approaches such as selective whole genome amplification (sWGA) are being used to enrich for parasite DNA. However, sWGA has had limited success in generating reliable sequencing data from low parasitaemia samples. In this study, enzymatic digestion with MspJI prior to sWGA and whole genome sequencing was evaluated to determine whether this approach improved genome coverage compared to sWGA alone. The potential of sWGA to cause amplification bias in polyclonal infections was also examined. METHODS: DNA extracted from laboratory-created dried blood spots was treated with a modification-dependent restriction endonuclease, MspJI, and filtered via vacuum filtration. Samples were then selectively amplified using a previously reported sWGA protocol and subjected to WGS. Genome coverage statistics were compared between the optimized sWGA approach and the previously reported sWGA approach performed in parallel. Differential amplification by sWGA was assessed by comparing WGS data generated from lab-created mixtures of parasite isolates, from the same geographical region, generated with or without sWGA. RESULTS: MspJI digestion did not enrich for parasite DNA. Samples that underwent vacuum filtration (without MspJI digestion) prior to sWGA had the highest parasite DNA concentration and displayed greater genome coverage compared to MspJI + sWGA and sWGA alone, particularly for low parasitaemia samples. The optimized sWGA (filtration + sWGA) approach was successfully used to generate WGS data from 218 non-leukocyte depleted field samples from Malawi. Sequences from lab-created mixtures of parasites did not show evidence of differential amplification of parasite strains compared to directly sequenced samples. CONCLUSION: This optimized sWGA approach is a reliable method to obtain WGS data from non-leukocyte depleted, low parasitaemia samples. The absence of amplification bias in data generated from mixtures of isolates from the same geographic region suggests that this approach can be appropriately used for molecular epidemiological studies.

Recrudescence, Reinfection, or Relapse? A More Rigorous Framework to Assess Chloroquine Efficacy for Plasmodium vivax Malaria.

Background: Plasmodium vivax resistance to chloroquine (CQ) has been reported worldwide, although the World Health Organization clinical drug efficacy studies protocol does not permit classification of patient outcomes. Methods: We enrolled 40 patients with P. vivax malaria in northeastern Cambodia, where >17% treatment failures were previously reported. Patients were treated with CQ (30 mg/kg) and followed for 2 months, with frequent clinical examination and capillary blood sample collection for microscopy, molecular parasite detection and genotyping, and drug concentration measurements. Reinfections were prevented by relocating patients to a transmission-free area. Results: P. vivax parasites were eliminated in all patients by day 3. Genomic analyses revealed that all clones in polyclonal infections were cleared at the same rate, indicating their equal susceptibility to CQ. CQ blood concentrations were below the therapeutic level in all recurrent infections (24 of 40 patients), which were efficiently cleared by a second course of CQ treatment. Genotyping (128 SNPs barcode) and sequences of entire parasite genome (Whole-Genome Sequencing, Illumina) indicated that two thirds (6 of 8) of the recurrent parasites resulted from heterologous relapses whose 50% are from by sibling/recombinant clones. Conclusions: No evidence of CQ resistance was observed. Our data suggest that P. vivax antimalarial drug resistance is likely overestimated and that the current guidelines for clinical drug studies of P. vivax malaria need to be revised.

Therapeutic and Transmission-Blocking Efficacy of Dihydroartemisinin/Piperaquine and Chloroquine against Plasmodium vivax Malaria, Cambodia.

We assessed the efficacy of standard 3-day courses of chloroquine and dihydroartemisinin/piperaquine against Plasmodium vivax malaria. Compared with chloroquine, dihydroartemisinin/piperaquine was faster in clearing asexual P. vivax parasites and blocking human-to-mosquito transmission. This drug combination was also more effective in preventing potential recurrences for >2 months.

Single-cell genomics for dissection of complex malaria infections.

Most malaria infections contain complex mixtures of distinct parasite lineages. These multiple-genotype infections (MGIs) impact virulence evolution, drug resistance, intra-host dynamics, and recombination, but are poorly understood. To address this we have developed a single-cell genomics approach to dissect MGIs. By combining cell sorting and whole-genome amplification (WGA), we are able to generate high-quality material from parasite-infected red blood cells (RBCs) for genotyping and next-generation sequencing. We optimized our approach through analysis of >260 single-cell assays. To quantify accuracy, we decomposed mixtures of known parasite genotypes and obtained highly accurate (>99%) single-cell genotypes. We applied this validated approach directly to infections of two major malaria species, Plasmodium falciparum, for which long term culture is possible, and Plasmodium vivax, for which no long-term culture is feasible. We demonstrate that our single-cell genomics approach can be used to generate parasite genome sequences directly from patient blood in order to unravel the complexity of P. vivax and P. falciparum infections. These methods open the door for large-scale analysis of within-host variation of malaria infections, and reveal information on relatedness and drug resistance haplotypes that is inaccessible through conventional sequencing of infections.

Replication fork integrity and intra-S phase checkpoint suppress gene amplification.

Gene amplification is a phenotype-causing form of chromosome instability and is initiated by DNA double-strand breaks (DSBs). Cells with mutant p53 lose G1/S checkpoint and are permissive to gene amplification. In this study we show that mammalian cells become proficient for spontaneous gene amplification when the function of the DSB repair protein complex MRN (Mre11/Rad50/Nbs1) is impaired. Cells with impaired MRN complex experienced severe replication stress and gained substrates for gene amplification during replication, as evidenced by the increase of replication-associated single-stranded breaks that were converted to DSBs most likely through replication fork reversal. Impaired MRN complex directly compromised ATM/ATR-mediated checkpoints and allowed cells to progress through cell cycle in the presence of DSBs. Such compromised intra-S phase checkpoints promoted gene amplification independently from mutant p53. Finally, cells adapted to endogenous replication stress by globally suppressing genes for DNA replication and cell cycle progression. Our results indicate that the MRN complex suppresses gene amplification by stabilizing replication forks and by securing DNA damage response to replication-associated DSBs.

Genomic Analysis Reveals a Common Breakpoint in Amplifications of the Plasmodium vivax Multidrug Resistance 1 Locus in Thailand.

In regions of coendemicity for Plasmodium falciparum and Plasmodium vivax where mefloquine is used to treat P. falciparum infection, drug pressure mediated by increased copy numbers of the multidrug resistance 1 gene (pvmdr1) may select for mefloquine-resistant P. vivax Surveillance is not undertaken routinely owing in part to methodological challenges in detection of gene amplification. Using genomic data on 88 P. vivax samples from western Thailand, we identified pvmdr1 amplification in 17 isolates, all exhibiting tandem copies of a 37.6-kilobase pair region with identical breakpoints. A novel breakpoint-specific polymerase chain reaction assay was designed to detect the amplification. The assay demonstrated high sensitivity, identifying amplifications in 13 additional, polyclonal infections. Application to 132 further samples identified the common breakpoint in all years tested (2003-2015), with a decline in prevalence after 2012 corresponding to local discontinuation of mefloquine regimens. Assessment of the structure of pvmdr1 amplification in other geographic regions will yield information about the population-specificity of the breakpoints and underlying amplification mechanisms.

Population genomics of the filarial nematode parasite Wuchereria bancrofti from mosquitoes.

Wuchereria bancrofti is a parasitic nematode and the primary cause of lymphatic filariasis--a disease specific to humans. W. bancrofti currently infects over 90 million people throughout the tropics and has been acknowledged by the world health organization as a vulnerable parasite. Current research has focused primarily on the clinical manifestations of disease and little is known about the evolutionary history of W. bancrofti. To improve upon knowledge of the evolutionary history of W. bancrofti, we whole genome sequenced 13 W. bancrofti larvae. We circumvent many of the difficulties of multiple infections by sampling larvae directly from mosquitoes that were experimentally inoculated with infected blood. To begin, we used whole genome data to reconstruct the historical population size. Our results support a history of fluctuating population sizes that can be correlated with human migration and fluctuating mosquito abundances. Next, we reconstructed the putative pedigree of W. bancrofti worms within an infection using the kinship coefficient. We deduced that there are full-sib and half-sib relationships residing within the same larval cohort. Through combined analysis of the mitochondrial and nuclear genomes we concluded that this is likely a results of polyandrous mating, the first time reported for W. bancrofti. Lastly, we scanned the genomes for signatures of natural selection. Annotation of putative selected regions identified proteins that may have aided in a parasitic life style or may have evolved to protect against current drug treatments. We discuss our results in the greater context of understanding the biology of an animal with a unique life history and ecology.

Whole-genome sequencing reveals absence of recent gene flow and separate demographic histories for Anopheles punctulatus mosquitoes in Papua New Guinea.

Anopheles mosquitoes are the vectors of several human diseases including malaria. In many malaria endemic areas, several species of Anopheles coexist, sometimes in the form of related sibling species that are morphologically indistinguishable. Determining the size and organization of Anopheles populations, and possible ongoing gene flow among them is important for malaria control and, in particular, for monitoring the spread of insecticide resistance alleles. However, these parameters have been difficult to evaluate in most Anopheles species due to the paucity of genetic data available. Here, we assess the extent of contemporary gene flow and historical variations in population size by sequencing and de novo assembling the genomes of wild-caught mosquitoes from four species of the Anopheles punctulatus group of Papua New Guinea. Our analysis of more than 50 Mb of orthologous DNA sequences revealed no evidence of contemporary gene flow among these mosquitoes. In addition, investigation of the demography of two of the An. punctulatus species revealed distinct population histories. Overall, our analyses suggest that, despite their similarities in morphology, behaviour and ecology, contemporary sympatric populations of An. punctulatus are evolving independently.

Single-cell analyses of polyclonal Plasmodium vivax infections and their consequences on parasite transmission.

Most Plasmodium vivax infections contain genetically distinct parasites, but the consequences of this polyclonality on the development of asexual parasites, their sexual differentiation, and their transmission remain unknown. We describe infections of Saimiri monkeys with two strains of P. vivax and the analyses of 117,350 parasites characterized by single cell RNA sequencing and individually genotyped. In our model, consecutive inoculations fail to establish polyclonal infections. By contrast, simultaneous inoculations of two strains lead to sustained polyclonal infections, although without detectable differences in parasite regulation or sexual commitment. Analyses of sporozoites dissected from mosquitoes fed on coinfected monkeys show that all genotypes are successfully transmitted to mosquitoes. However, after sporozoite inoculation, not all genotypes contribute to the subsequent blood infections, highlighting an important bottleneck during pre-erythrocytic development. Overall, these studies provide new insights on the mechanisms regulating the establishment of polyclonal P. vivax infections and their consequences for disease transmission.

Gene expression analyses reveal differences in children's response to malaria according to their age.

In Bandiagara, Mali, children experience on average two clinical malaria episodes per year. However, even in the same transmission area, the number of uncomplicated symptomatic infections, and their parasitemia, can vary dramatically among children. We simultaneously characterize host and parasite gene expression profiles from 136 Malian children with symptomatic falciparum malaria and examine differences in the relative proportion of immune cells and parasite stages, as well as in gene expression, associated with infection and or patient characteristics. Parasitemia explains much of the variation in host and parasite gene expression, and infections with higher parasitemia display proportionally more neutrophils and fewer T cells, suggesting parasitemia-dependent neutrophil recruitment and/or T cell extravasation to secondary lymphoid organs. The child's age also strongly correlates with variations in gene expression: Plasmodium falciparum genes associated with age suggest that older children carry more male gametocytes, while variations in host gene expression indicate a stronger innate response in younger children and stronger adaptive response in older children. These analyses highlight the variability in host responses and parasite regulation during P. falciparum symptomatic infections and emphasize the importance of considering the children's age when studying and treating malaria infections.

Tick hemocytes have a pleiotropic role in microbial infection and arthropod fitness.

Uncovering the complexity of systems in non-model organisms is critical for understanding arthropod immunology. Prior efforts have mostly focused on Dipteran insects, which only account for a subset of existing arthropod species in nature. Here we use and develop advanced techniques to describe immune cells (hemocytes) from the clinically relevant tick Ixodes scapularis at a single-cell resolution. We observe molecular alterations in hemocytes upon feeding and infection with either the Lyme disease spirochete Borrelia burgdorferi or the rickettsial agent Anaplasma phagocytophilum. We reveal hemocyte clusters exhibiting defined signatures related to immunity, metabolism, and proliferation. Depletion of phagocytic hemocytes affects hemocytin and astakine levels, two I. scapularis hemocyte markers, impacting blood-feeding, molting behavior, and bacterial acquisition. Mechanistically, astakine alters hemocyte proliferation, whereas hemocytin affects the c-Jun N-terminal kinase (JNK) signaling pathway in I. scapularis. Altogether, we discover a role for tick hemocytes in immunophysiology and provide a valuable resource for comparative biology in arthropods.

Cathepsin B is a novel gender-dependent determinant of cholesterol absorption from the intestine.

We used a mouse C57BL/6J×CASA/Rk intercross to map a locus on chromosome 14 that displayed a gender-dependent effect on cholesterol absorption from the intestine. Studies in congenic animals revealed a complex locus with multiple operating genetic determinants resulting in alternating gender-dependent phenotypic effects. Fine-mapping narrowed the locus to a critical 6.3 Mb interval. Female subcongenics, but not males, of the critical interval displayed a decrease of 33% in cholesterol absorption. RNA-Seq analysis of female subcongenic jejunum revealed that cysteine protease cathepsin B (Ctsb) is a candidate to explain the interval effect. Consistent with the phenotype in critical interval subcongenics, female Ctsb knockout mice, but not males, displayed a decrease of 31% in cholesterol absorption. Although studies in Ctsb knockouts revealed a gender-dependent effect on cholesterol absorption, further fine-mapping dismissed a role for Ctsb in determining the effect of the critical 6.3 Mb interval on cholesterol absorption.

Mitochondrial genome sequences reveal deep divergences among Anopheles punctulatus sibling species in Papua New Guinea.

BACKGROUND: Members of the Anopheles punctulatus group (AP group) are the primary vectors of human malaria in Papua New Guinea. The AP group includes 13 sibling species, most of them morphologically indistinguishable. Understanding why only certain species are able to transmit malaria requires a better comprehension of their evolutionary history. In particular, understanding relationships and divergence times among Anopheles species may enable assessing how malaria-related traits (e.g. blood feeding behaviours, vector competence) have evolved. METHODS: DNA sequences of 14 mitochondrial (mt) genomes from five AP sibling species and two species of the Anopheles dirus complex of Southeast Asia were sequenced. DNA sequences from all concatenated protein coding genes (10,770 bp) were then analysed using a Bayesian approach to reconstruct phylogenetic relationships and date the divergence of the AP sibling species. RESULTS: Phylogenetic reconstruction using the concatenated DNA sequence of all mitochondrial protein coding genes indicates that the ancestors of the AP group arrived in Papua New Guinea 25 to 54 million years ago and rapidly diverged to form the current sibling species. CONCLUSION: Through evaluation of newly described mt genome sequences, this study has revealed a divergence among members of the AP group in Papua New Guinea that would significantly predate the arrival of humans in this region, 50 thousand years ago. The divergence observed among the mtDNA sequences studied here may have resulted from reproductive isolation during historical changes in sea-level through glacial minima and maxima. This leads to a hypothesis that the AP sibling species have evolved independently for potentially thousands of generations. This suggests that the evolution of many phenotypes, such as insecticide resistance will arise independently in each of the AP sibling species studied here.

Neanderthals in central Asia and Siberia.

Morphological traits typical of Neanderthals began to appear in European hominids at least 400,000 years ago and about 150,000 years ago in western Asia. After their initial appearance, such traits increased in frequency and the extent to which they are expressed until they disappeared shortly after 30,000 years ago. However, because most fossil hominid remains are fragmentary, it can be difficult or impossible to determine unambiguously whether a fossil is of Neanderthal origin. This limits the ability to determine when and where Neanderthals lived. To determine how far to the east Neanderthals ranged, we determined mitochondrial DNA (mtDNA) sequences from hominid remains found in Uzbekistan and in the Altai region of southern Siberia. Here we show that the DNA sequences from these fossils fall within the European Neanderthal mtDNA variation. Thus, the geographic range of Neanderthals is likely to have extended at least 2,000 km further to the east than commonly assumed.

A genome-wide association study identifies novel risk loci for type 2 diabetes.

Type 2 diabetes mellitus results from the interaction of environmental factors with a combination of genetic variants, most of which were hitherto unknown. A systematic search for these variants was recently made possible by the development of high-density arrays that permit the genotyping of hundreds of thousands of polymorphisms. We tested 392,935 single-nucleotide polymorphisms in a French case-control cohort. Markers with the most significant difference in genotype frequencies between cases of type 2 diabetes and controls were fast-tracked for testing in a second cohort. This identified four loci containing variants that confer type 2 diabetes risk, in addition to confirming the known association with the TCF7L2 gene. These loci include a non-synonymous polymorphism in the zinc transporter SLC30A8, which is expressed exclusively in insulin-producing beta-cells, and two linkage disequilibrium blocks that contain genes potentially involved in beta-cell development or function (IDE-KIF11-HHEX and EXT2-ALX4). These associations explain a substantial portion of disease risk and constitute proof of principle for the genome-wide approach to the elucidation of complex genetic traits.