Publications

Mass drug administration with praziquantel (PZQ) monotherapy is considered the mainstay for control and elimination of the parasites causing schistosomiasis in humans. This drug shows imperfect cure rates in the field, and parasites showing reduced PZQ response can be selected in the laboratory, but the extent of resistance in Schistosoma mansoni populations is unknown. We examined the genetic basis of the variation in response in a PZQ-selected S. mansoni population (SmLE-PZQ-R) in which 35% of the parasitic worms survive high-dose PZQ (73 micrograms per milliliter) treatment. We used genome-wide association to map loci underlying PZQ response and identified a transient receptor potential (Sm.TRPMPZQ) channel (Smp_246790) within the major chromosome 3 peak that is activated by nanomolar concentrations of PZQ. The PZQ response showed recessive inheritance and marker-assisted selection of parasites at a single Sm.TRPMPZQ SNP that produced populations of PZQ-enriched resistant (PZQ-ER) and PZQ-enriched sensitive (PZQ-ES) parasites, exhibiting >377-fold difference in PZQ response. The PZQ-ER parasites survived treatment in rodents at higher frequencies compared with PZQ-ES, and resistant parasites exhibited 2.25-fold lower expression of Sm.TRPMPZQ relative to sensitive parasites. Specific chemical blockers of Sm.TRPMPZQ enhanced PZQ resistance, whereas Sm.TRPMPZQ activators increased sensitivity. We surveyed Sm.TRPMPZQ sequence variations in 259 parasites from different global sites and identified one nonsense mutation that resulted in a truncated protein with no PZQ binding site. Our results demonstrate that Sm.TRPMPZQ underlies variation in PZQ responses in S. mansoni and provides an approach for monitoring emerging PZQ-resistant alleles in schistosome elimination programs.

The emergence of mutant K13-mediated artemisinin (ART) resistance in Plasmodium falciparum malaria parasites has led to widespread treatment failures across Southeast Asia. In Africa, K13-propeller genotyping confirms the emergence of the R561H mutation in Rwanda and highlights the continuing dominance of wild-type K13 elsewhere. Using gene editing, we show that R561H, along with C580Y and M579I, confer elevated in vitro ART resistance in some African strains, contrasting with minimal changes in ART susceptibility in others. C580Y and M579I cause substantial fitness costs, which may slow their dissemination in high-transmission settings, in contrast with R561H that in African 3D7 parasites is fitness neutral. In Cambodia, K13 genotyping highlights the increasing spatio-temporal dominance of C580Y. Editing multiple K13 mutations into a panel of Southeast Asian strains reveals that only the R561H variant yields ART resistance comparable to C580Y. In Asian Dd2 parasites C580Y shows no fitness cost, in contrast with most other K13 mutations tested, including R561H. Editing of point mutations in ferredoxin or mdr2, earlier associated with resistance, has no impact on ART susceptibility or parasite fitness. These data underline the complex interplay between K13 mutations, parasite survival, growth and genetic background in contributing to the spread of ART resistance.

Non-human primate models will expedite therapeutics and vaccines for coronavirus disease 2019 (COVID-19) to clinical trials. Here, we compare acute severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in young and old rhesus macaques, baboons and old marmosets. Macaques had clinical signs of viral infection, mild to moderate pneumonitis and extra-pulmonary pathologies, and both age groups recovered in two weeks. Baboons had prolonged viral RNA shedding and substantially more lung inflammation compared with macaques. Inflammation in bronchoalveolar lavage was increased in old versus young baboons. Using techniques including computed tomography imaging, immunophenotyping, and alveolar/peripheral cytokine response and immunohistochemical analyses, we delineated cellular immune responses to SARS-CoV-2 infection in macaque and baboon lungs, including innate and adaptive immune cells and a prominent type-I interferon response. Macaques developed T-cell memory phenotypes/responses and bystander cytokine production. Old macaques had lower titres of SARS-CoV-2-specific IgG antibody levels compared with young macaques. Acute respiratory distress in macaques and baboons recapitulates the progression of COVID-19 in humans, making them suitable as models to test vaccines and therapies.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the viral pathogen responsible for the current coronavirus disease 2019 (COVID-19) pandemic. As of 19 May 2021, John Hopkins University's COVID-19 tracking platform reported 3.3 million deaths associated with SARS-CoV-2 infection. Currently, the World Health Organization has granted emergency use listing (EUL) to six COVID-19 vaccine candidates. However, much of the pathogenesis observed during SARS-CoV-2 infection remains elusive. To gain insight into the contribution of individual accessory open reading frame (ORF) proteins in SARS-CoV-2 pathogenesis, we used our recently described reverse-genetics system approach to successfully engineer recombinant SARS-CoV-2 (rSARS-CoV-2) constructs; we removed individual viral ORF3a, -6, -7a, -7b, and -8 proteins from them, and we characterized the resulting recombinant viruses in vitro and in vivo. Our results indicate differences in plaque morphology, with ORF-deficient (ΔORF) viruses producing smaller plaques than those of the wild type (rSARS-CoV-2/WT). However, growth kinetics of ΔORF viruses were like those of rSARS-CoV-2/WT. Interestingly, infection of K18 human angiotensin-converting enzyme 2 (hACE2) transgenic mice with the ΔORF rSARS-CoV-2s identified ORF3a and ORF6 as the major contributors of viral pathogenesis, while ΔORF7a, ΔORF7b, and ΔORF8 rSARS-CoV-2s induced pathology comparable to that of rSARS-CoV-2/WT. This study demonstrates the robustness of our reverse-genetics system to generate rSARS-CoV-2 constructs and the major role for ORF3a and ORF6 in viral pathogenesis, providing important information for the generation of attenuated forms of SARS-CoV-2 for their implementation as live attenuated vaccines for the treatment of SARS-CoV-2 infection and associated COVID-19. IMPORTANCE Despite great efforts put forward worldwide to combat the current coronavirus disease 2019 (COVID-19) pandemic, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to be a human health and socioeconomic threat. Insights into the pathogenesis of SARS-CoV-2 and the contribution of viral proteins to disease outcome remain elusive. Our study aims (i) to determine the contribution of SARS-CoV-2 accessory open reading frame (ORF) proteins to viral pathogenesis and disease outcome and (ii) to develop a synergistic platform combining our robust reverse-genetics system to generate recombinant SARS-CoV-2 constructs with a validated rodent model of infection and disease. We demonstrate that SARS-CoV-2 ORF3a and ORF6 contribute to lung pathology and ultimately disease outcome in K18 hACE2 transgenic mice, while ORF7a, ORF7b, and ORF8 have little impact on disease outcome. Moreover, our combinatory platform serves as a foundation for generating attenuated forms of the virus to develop live attenuated vaccines for the treatment of SARS-CoV-2.

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

Human schistosomiasis is a debilitating, life-threatening disease affecting more than 229 million people in as many as 78 countries. There is only one drug of choice effective against all three major species of Schistosoma, praziquantel (PZQ). However, as with many monotherapies, evidence for resistance is emerging in the field and can be selected for in the laboratory. Previously used therapies include oxamniquine (OXA), but shortcomings such as drug resistance and affordability resulted in discontinuation. Employing a genetic, biochemical and molecular approach, a sulfotransferase (SULT-OR) was identified as responsible for OXA drug resistance. By crystallizing SmSULT- OR with OXA, the mode of action of OXA was determined. This information allowed a rational approach to novel drug design. Our team approach with schistosome biologists, medicinal chemists, structural biologists and geneticists has enabled us to develop and test novel drug derivatives of OXA to treat this disease. Using an iterative process for drug development, we have successfully identified derivatives that are effective against all three species of the parasite. One derivative CIDD-0149830 kills 100% of all three human schistosome species within 5 days. The goal is to generate a second therapeutic with a different mode of action that can be used in conjunction with praziquantel to overcome the ever-growing threat of resistance and improve efficacy. The ability and need to design, screen, and develop future, affordable therapeutics to treat human schistosomiasis is critical for successful control program outcomes.

Schistosomiasis control efforts in Nigeria received a boost in 2016 when Merck Group made the largest single donation of praziquantel to an African country. We examined urine samples from 2,023 school age children from 15 locations in 10 states and an Internally Displaced Person's (IDP) camp in Nigeria. We recorded an overall Schistosoma haematobium prevalence of 10.4% in the 10 states that ranged between 6 - 37%, while prevalence in the IDP camp was 2.9%. The highest infection prevalence (37%) recorded was from the population in Wasai Dam area in Minjibir (Kano State), while five locations had no positive urine samples. We observed heavy intensity of infection (≥ 50 eggs/10 ml urine) in 87.9% of infected samples and co-occurrence of the eggs of S. haematobium and S. mansoni in urine for two participants. The overall prevalence we recorded is slightly above the national average (9.5%) reported in 2015. Our findings indicate that despite the ongoing administration of praziquantel in Nigeria, urogenital schistosomiasis is still prevalent with heavy intensity of infection. Large-scale epidemiological monitoring is required to monitor the efficacy of schistosomiasis control in Nigeria.

Population genomics of bulk malaria infections is unable to examine intrahost evolution; therefore, most work has focused on the role of recombination in generating genetic variation. We used single-cell sequencing protocol for low-parasitaemia infections to generate 406 near-complete single Plasmodium vivax genomes from 11 patients sampled during sequential febrile episodes. Parasite genomes contain hundreds of de novo mutations, showing strong signatures of selection, which are enriched in the ApiAP2 family of transcription factors, known targets of adaptation. Comparing 315 P. falciparum single-cell genomes from 15 patients with our P. vivax data, we find broad complementary patterns of de novo mutation at the gene and pathway level, revealing the importance of within-host evolution during malaria infections.

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the pathogen responsible of coronavirus disease 2019 (COVID-19), has devastated public health services and economies worldwide. Despite global efforts to contain the COVID-19 pandemic, SARS-CoV-2 is now found in over 200 countries and has caused an upward death toll of over 1 million human lives as of November 2020. To date, only one Food and Drug Administration (FDA)-approved therapeutic drug (Remdesivir) and a monoclonal antibody, MAb (Bamlanivimab) are available for the treatment of SARS-CoV-2. As with other viruses, studying SARS-CoV-2 requires the use of secondary approaches to detect the presence of the virus in infected cells. To overcome this limitation, we have generated replication-competent recombinant (r)SARS-CoV-2 expressing fluorescent (Venus or mCherry) or bioluminescent (Nluc) reporter genes. Vero E6 cells infected with reporter-expressing rSARS-CoV-2 can be easily detected via fluorescence or luciferase expression and display a good correlation between reporter gene expression and viral replication. Moreover, rSARS-CoV-2 expressing reporter genes have comparable plaque sizes and growth kinetics to those of wild-type virus, rSARS-CoV-2/WT. We used these reporter-expressing rSARS-CoV-2 to demonstrate their feasibility to identify neutralizing antibodies (NAbs) or antiviral drugs. Our results demonstrate that reporter-expressing rSARS-CoV-2 represent an excellent option to identify therapeutics for the treatment of SARS-CoV-2, where reporter gene expression can be used as valid surrogates to track viral infection. Moreover, the ability to manipulate the viral genome opens the feasibility of generating viruses expressing foreign genes for their use as vaccines for the treatment of SARS-CoV-2 infection.IMPORTANCE Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the pathogen that causes coronavirus disease 2019 (COVID-19), has significantly impacted the human health and economic status worldwide. There is an urgent need to identify effective prophylactics and therapeutics for the treatment of SARS-CoV-2 infection and associated COVID-19 disease. The use of fluorescent- or luciferase-expressing reporter expressing viruses has significantly advanced viral research. Here, we generated recombinant (r)SARS-CoV-2 expressing fluorescent (Venus and mCherry) or luciferase (Nluc) reporter genes and demonstrate that they represent an excellent option to track viral infections in vitro. Importantly, reporter-expressing rSARS-CoV-2 display similar growth kinetics and plaque phenotype that their wild-type counterpart (rSARS-CoV-2/WT), demonstrating their feasibility to identify drugs and/or neutralizing antibodies (NAbs) for the therapeutic treatment of SARS-CoV-2. Henceforth, these reporter-expressing rSARS-CoV-2 can be used to interrogate large libraries of compounds and/or monoclonal antibodies (MAb), in high-throughput screening settings, to identify those with therapeutic potential against SARS-CoV-2.

Genetic crosses are most powerful for linkage analysis when progeny numbers are high, parental alleles segregate evenly and numbers of inbred progeny are minimized. We previously developed a novel genetic crossing platform for the human malaria parasite Plasmodium falciparum, an obligately sexual, hermaphroditic protozoan, using mice carrying human hepatocytes (the human liver-chimeric FRG NOD huHep mouse) as the vertebrate host. We report on two genetic crosses-(1) an allopatric cross between a laboratory-adapted parasite (NF54) of African origin and a recently patient-derived Asian parasite, and (2) a sympatric cross between two recently patient-derived Asian parasites. We generated 144 unique recombinant clones from the two crosses, doubling the number of unique recombinant progeny generated in the previous 30 years. The allopatric African/Asian cross has minimal levels of inbreeding and extreme segregation distortion, while in the sympatric Asian cross, inbred progeny predominate and parental alleles segregate evenly. Using simulations, we demonstrate that these progeny provide the power to map small-effect mutations and epistatic interactions. The segregation distortion in the allopatric cross slightly erodes power to detect linkage in several genome regions. We greatly increase the power and the precision to map biomedically important traits with these new large progeny panels.