Publications

Characterising connectivity between geographically separated biological populations is a common goal in many fields. Recent approaches to understanding connectivity between malaria parasite populations, with implications for disease control efforts, have used estimates of relatedness based on identity-by-descent (IBD). However, uncertainty around estimated relatedness has not been accounted for. IBD-based relatedness estimates with uncertainty were computed for pairs of monoclonal Plasmodium falciparum samples collected from five cities on the Colombian-Pacific coast where long-term clonal propagation of P. falciparum is frequent. The cities include two official ports, Buenaventura and Tumaco, that are separated geographically but connected by frequent marine traffic. Fractions of highly-related sample pairs (whose classification using a threshold accounts for uncertainty) were greater within cities versus between. However, based on both highly-related fractions and on a threshold-free approach (Wasserstein distances between parasite populations) connectivity between Buenaventura and Tumaco was disproportionally high. Buenaventura-Tumaco connectivity was consistent with transmission events involving parasites from five clonal components (groups of statistically indistinguishable parasites identified under a graph theoretic framework). To conclude, P. falciparum population connectivity on the Colombian-Pacific coast abides by accessibility not isolation-by-distance, potentially implicating marine traffic in malaria transmission with opportunities for targeted intervention. Further investigations are required to test this hypothesis. For the first time in malaria epidemiology (and to our knowledge in ecological and epidemiological studies more generally), we account for uncertainty around estimated relatedness (an important consideration for studies that plan to use genotype versus whole genome sequence data to estimate IBD-based relatedness); we also use threshold-free methods to compare parasite populations and identify clonal components. Threshold-free methods are especially important in analyses of malaria parasites and other recombining organisms with mixed mating systems where thresholds do not have clear interpretation (e.g. due to clonal propagation) and thus undermine the cross-comparison of studies.

Vaccine and antiviral development against SARS-CoV-2 infection or COVID-19 disease would benefit from validated small animal models. Here, we show that transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) by the human cytokeratin 18 promoter (K18 hACE2) represent a susceptible rodent model. K18 hACE2 transgenic mice succumbed to SARS-CoV-2 infection by day 6, with virus detected in lung airway epithelium and brain. K18 ACE2 transgenic mice produced a modest TH1/2/17 cytokine storm in the lung and spleen that peaked by day 2, and an extended chemokine storm that was detected in both lungs and brain. This chemokine storm was also detected in the brain at day 6. K18 hACE2 transgenic mice are, therefore, highly susceptible to SARS-CoV-2 infection and represent a suitable animal model for the study of viral pathogenesis, and for identification and characterization of vaccines (prophylactic) and antivirals (therapeutics) for SARS-CoV-2 infection and associated severe COVID-19 disease.

The first experimental crosses carried out with the human malaria parasite Plasmodium falciparum played a key role in determining the genetic loci responsible for drug resistance, virulence, invasion, growth rate, and transmission. These crosses relied on splenectomized chimpanzees to complete the liver stage of the parasite's life cycle and the subsequent transition to asexual blood stage culture followed by cloning of recombinant progeny in vitro. Crosses can now be routinely carried out using human-liver-chimeric mice infused with human erythrocytes to generate hundreds of unique recombinant progeny for genetic linkage mapping, bulk segregant analysis, and high-throughput 'omics readouts. The high number of recombinant progeny should allow for unprecedented power and efficiency in the execution of a systems genetics approach to study P. falciparum biology.

Infectious coronavirus (CoV) disease 2019 (COVID-19) emerged in the city of Wuhan (China) in December 2019, causing a pandemic that has dramatically impacted public health and socioeconomic activities worldwide. A previously unknown coronavirus, severe acute respiratory syndrome CoV-2 (SARS-CoV-2), has been identified as the causative agent of COVID-19. To date, there are no U.S. Food and Drug Administration (FDA)-approved vaccines or therapeutics available for the prevention or treatment of SARS-CoV-2 infection and/or associated COVID-19 disease, which has triggered a large influx of scientific efforts to develop countermeasures to control SARS-CoV-2 spread. To contribute to these efforts, we have developed an infectious cDNA clone of the SARS-CoV-2 USA-WA1/2020 strain based on the use of a bacterial artificial chromosome (BAC). Recombinant SARS-CoV-2 (rSARS-CoV-2) was readily rescued by transfection of the BAC into Vero E6 cells. Importantly, BAC-derived rSARS-CoV-2 exhibited growth properties and plaque sizes in cultured cells comparable to those of the natural SARS-CoV-2 isolate. Likewise, rSARS-CoV-2 showed levels of replication similar to those of the natural isolate in nasal turbinates and lungs of infected golden Syrian hamsters. This is, to our knowledge, the first BAC-based reverse genetics system for the generation of infectious rSARS-CoV-2 that displays features in vivo similar to those of a natural viral isolate. This SARS-CoV-2 BAC-based reverse genetics will facilitate studies addressing several important questions in the biology of SARS-CoV-2, as well as the identification of antivirals and development of vaccines for the treatment of SARS-CoV-2 infection and associated COVID-19 disease.IMPORTANCE The pandemic coronavirus (CoV) disease 2019 (COVID-19) caused by severe acute respiratory syndrome CoV-2 (SARS-CoV-2) is a major threat to global human health. To date, there are no approved prophylactics or therapeutics available for COVID-19. Reverse genetics is a powerful approach to understand factors involved in viral pathogenesis, antiviral screening, and vaccine development. In this study, we describe the feasibility of generating recombinant SARS-CoV-2 (rSARS-CoV-2) by transfection of a single bacterial artificial chromosome (BAC). Importantly, rSARS-CoV-2 possesses the same phenotype as the natural isolate in vitro and in vivo This is the first description of a BAC-based reverse genetics system for SARS-CoV-2 and the first time that an rSARS-CoV-2 isolate has been shown to be phenotypically identical to a natural isolate in a validated animal model of SARS-CoV-2 infection. The BAC-based reverse genetics approach will facilitate the study of SARS-CoV-2 and the development of prophylactics and therapeutics for the treatment of COVID-19.

Human schistosomiasis is a disease which globally affects over 229 million people. Three major species affecting humans are Schistosoma mansoni, S. haematobium and S. japonicum. Previous treatment of S. mansoni includes the use of oxamniquine (OXA), a prodrug that is enzymatically activated in S. mansoni but is ineffective against S. haematobium and S. japonicum. The OXA activating enzyme was identified and crystallized, as being a S. mansoni sulfotransferase (SmSULT). S. haematobium and S. japonicum possess homologs of SmSULT (ShSULT and SjSULT) begging the question; why does oxamniquine fail to kill S. haematobium and S. japonicum adult worms? Investigation of the molecular structures of the sulfotransferases indicates that structural differences, specifically in OXA contact residues, do not abrogate OXA binding in the active sites as previously hypothesized. Data presented argue that the ability of SULTs to sulfate and thus activate OXA and its derivatives is linked to the ability of OXA to fit in the binding pocket to allow the transfer of a sulfur group.

Currently there is only one method of treatment for human schistosomiasis, the drug praziquantel. Strong selective pressure has caused a serious concern for a rise in resistance to praziquantel leading to the necessity for additional pharmaceuticals, with a distinctly different mechanism of action, to be used in combination therapy with praziquantel. Previous treatment of Schistosoma mansoni included the use of oxamniquine (OXA), a prodrug that is enzymatically activated in S. mansoni but is ineffective against S. haematobium and S. japonicum. The oxamniquine activating enzyme was identified as a S. mansoni sulfotransferase (SmSULT-OR). Structural data have allowed for directed drug development in reengineering oxamniquine to be effective against S. haematobium and S. japonicum. Guided by data from X-ray crystallographic studies and Schistosoma worm killing assays on oxamniquine, our structure-based drug design approach produced a robust SAR program that tested over 300 derivatives and identified several new lead compounds with effective worm killing in vitro. Previous studies resulted in the discovery of compound CIDD-0066790, which demonstrated broad-species activity in killing of schistosome species. As these compounds are racemic mixtures, we tested and demonstrate that the R enantiomer CIDD-007229 kills S. mansoni, S. haematobium and S. japonicum better than the parent drug (CIDD-0066790). The search for derivatives that kill better than CIDD-0066790 has resulted in a derivative (CIDD- 149830) that kills 100% of S. mansoni, S. haematobium and S. japonicum adult worms within 7 days. We hypothesize that the difference in activation and thus killing by the derivatives is due to the ability of the derivative to fit in the binding pocket of each sulfotransferase (SmSULT-OR, ShSULT-OR, SjSULT-OR) and to be efficiently sulfated. The purpose of this research is to develop a second drug to be used in conjunction with praziquantel to treat the major human species of Schistosoma. Collectively, our findings show that CIDD-00149830 and CIDD-0072229 are promising novel drugs for the treatment of human schistosomiasis and strongly support further development and in vivo testing.

BACKGROUND: Tracking and understanding artemisinin resistance is key for preventing global setbacks in malaria eradication efforts. The ring-stage survival assay (RSA) is the current gold standard for in vitro artemisinin resistance phenotyping. However, the RSA has several drawbacks: it is relatively low throughput, has high variance due to microscopy readout, and correlates poorly with the current benchmark for in vivo resistance, patient clearance half-life post-artemisinin treatment. Here a modified RSA is presented, the extended Recovery Ring-stage Survival Assay (eRRSA), using 15 cloned patient isolates from Southeast Asia with a range of patient clearance half-lives, including parasite isolates with and without kelch13 mutations.

METHODS: Plasmodium falciparum cultures were synchronized with single layer Percoll during the schizont stage of the intraerythrocytic development cycle. Cultures were left to reinvade to early ring-stage and parasitaemia was quantified using flow cytometry. Cultures were diluted to 2% haematocrit and 0.5% parasitaemia in a 96-well plate to start the assay, allowing for increased throughput and decreased variability between biological replicates. Parasites were treated with 700 nM of dihydroartemisinin or 0.02% dimethyl sulfoxide (DMSO) for 6 h, washed three times in drug-free media, and incubated for 66 or 114 h, when samples were collected and frozen for PCR amplification. A SYBR Green-based quantitative PCR method was used to quantify the fold-change between treated and untreated samples.

RESULTS: 15 cloned patient isolates from Southeast Asia with a range of patient clearance half-lives were assayed using the eRRSA. Due to the large number of pyknotic and dying parasites at 66 h post-exposure (72 h sample), parasites were grown for an additional cell cycle (114 h post-exposure, 120 h sample), which drastically improved correlation with patient clearance half-life compared to the 66 h post-exposure sample. A Spearman correlation of - 0.8393 between fold change and patient clearance half-life was identified in these 15 isolates from Southeast Asia, which is the strongest correlation reported to date.

CONCLUSIONS: eRRSA drastically increases the efficiency and accuracy of in vitro artemisinin resistance phenotyping compared to the traditional RSA, which paves the way for extensive in vitro phenotyping of hundreds of artemisinin resistant parasites.

In high-transmission regions, we expect parasite lineages within complex malaria infections to be unrelated due to parasite inoculations from different mosquitoes. This project was designed to test this prediction. We generated 485 single-cell genome sequences from fifteen P. falciparum malaria patients from Chikhwawa, Malawi-an area of intense transmission. Patients harbored up to seventeen unique parasite lineages. Surprisingly, parasite lineages within infections tend to be closely related, suggesting that superinfection by repeated mosquito bites is rarer than co-transmission of parasites from a single mosquito. Both closely and distantly related parasites comprise an infection, suggesting sequential transmission of complex infections between multiple hosts. We identified tetrads and reconstructed parental haplotypes, which revealed the inbred ancestry of infections and non-Mendelian inheritance. Our analysis suggests strong barriers to secondary infection and outbreeding amongst malaria parasites from a high transmission setting, providing unexpected insights into the biology and transmission of malaria.

Do mutations required for adaptation occur de novo, or are they segregating within populations as standing genetic variation? This question is key to understanding adaptive change in nature, and has important practical consequences for the evolution of drug resistance. We provide evidence that alleles conferring resistance to oxamniquine (OXA), an antischistosomal drug, are widespread in natural parasite populations under minimal drug pressure and predate OXA deployment. OXA has been used since the 1970s to treat Schistosoma mansoni infections in the New World where S. mansoni established during the slave trade. Recessive loss-of-function mutations within a parasite sulfotransferase (SmSULT-OR) underlie resistance, and several verified resistance mutations, including a deletion (p.E142del), have been identified in the New World. Here we investigate sequence variation in SmSULT-OR in S. mansoni from the Old World, where OXA has seen minimal usage. We sequenced exomes of 204 S. mansoni parasites from West Africa, East Africa and the Middle East, and scored variants in SmSULT-OR and flanking regions. We identified 39 non-synonymous SNPs, 4 deletions, 1 duplication and 1 premature stop codon in the SmSULT-OR coding sequence, including one confirmed resistance deletion (p.E142del). We expressed recombinant proteins and used an in vitro OXA activation assay to functionally validate the OXA-resistance phenotype for four predicted OXA-resistance mutations. Three aspects of the data are of particular interest: (i) segregating OXA-resistance alleles are widespread in Old World populations (4.29-14.91% frequency), despite minimal OXA usage, (ii) two OXA-resistance mutations (p.W120R, p.N171IfsX28) are particularly common (>5%) in East African and Middle-Eastern populations, (iii) the p.E142del allele has identical flanking SNPs in both West Africa and Puerto Rico, suggesting that parasites bearing this allele colonized the New World during the slave trade and therefore predate OXA deployment. We conclude that standing variation for OXA resistance is widespread in S. mansoni.

Introgression among parasite species has the potential to transfer traits of biomedical importance across species boundaries. The parasitic blood fluke Schistosoma haematobium causes urogenital schistosomiasis in humans across sub-Saharan Africa. Hybridization with other schistosome species is assumed to occur commonly, because genetic crosses between S. haematobium and livestock schistosomes, including S. bovis, can be staged in the laboratory, and sequencing of mtDNA and rDNA amplified from microscopic miracidia larvae frequently reveals markers from different species. However, the frequency, direction, age, and genomic consequences of hybridization are unknown. We hatched miracidia from eggs and sequenced the exomes from 96 individual S. haematobium miracidia from infected patients from Niger and the Zanzibar archipelago. These data revealed no evidence for contemporary hybridization between S. bovis and S. haematobium in our samples. However, all Nigerien S. haematobium genomes sampled show hybrid ancestry, with 3.3-8.2% of their nuclear genomes derived from S. bovis, providing evidence of an ancient introgression event that occurred at least 108-613 generations ago. Some S. bovis-derived alleles have spread to high frequency or reached fixation and show strong signatures of directional selection; the strongest signal spans a single gene in the invadolysin gene family (Chr. 4). Our results suggest that S. bovis/S. haematobium hybridization occurs rarely but demonstrate profound consequences of ancient introgression from a livestock parasite into the genome of S. haematobium, the most prevalent schistosome species infecting humans.