Plasmodium falciparum, the causative agent of severe human malaria, employs antigenic variation to avoid host immunity. Antigenic variation is achieved by transcriptional switching amongst polymorphic var genes, enforced by epigenetic modification of chromatin. The histone-modifying 'sirtuin' enzymes PfSir2a and PfSir2b have been implicated in this process. Disparate patterns of var expression have been reported in patient isolates as well as in cultured strains. We examined var expression in three commonly used laboratory strains (3D7, NF54 and FCR-3) in parallel. NF54 parasites express significantly lower levels of var genes compared to 3D7, despite the fact that 3D7 was originally a clone of the NF54 strain. To investigate whether this was linked to the expression of sirtuins, genetic disruption of both sirtuins was attempted in all three strains. No dramatic changes in var gene expression occurred in NF54 or FCR-3 following PfSir2b disruption, contrasting with previous observations in 3D7. In 3D7, complementation of the PfSir2a genetic disruption resulted in a significant decrease in previously-elevated var gene expression levels, but with the continued expression of multiple var genes. Finally, rearranged chromosomes were observed in the 3D7 PfSir2a knockout line. Our results focus on the potential for parasite genetic background to contribute to sirtuin function in regulating virulence gene expression and suggest a potential role for sirtuins in maintaining genome integrity.

We have shown previously that enteral administration of 2, 4, 6-trinitrobenzene sulfonic acid in 10-d-old C57BL/6 pups produces an acute necrotizing enterocolitis with histopathological and inflammatory changes similar to human necrotizing enterocolitis (NEC). To determine whether murine neonatal 2, 4, 6-trinitrobenzene sulfonic acid (TNBS)-mediated intestinal injury could be used as a NEC model, we compared gene expression profiles of TNBS-mediated intestinal injury and NEC.

We sequenced and annotated the genomes of four P. vivax strains collected from disparate geographic locations, tripling the number of genome sequences available for this understudied parasite and providing the first genome-wide perspective of global variability in this species. We observe approximately twice as much SNP diversity among these isolates as we do among a comparable collection of isolates of P. falciparum, a malaria-causing parasite that results in higher mortality. This indicates a distinct history of global colonization and/or a more stable demographic history for P. vivax relative to P. falciparum, which is thought to have undergone a recent population bottleneck. The SNP diversity, as well as additional microsatellite and gene family variability, suggests a capacity for greater functional variation in the global population of P. vivax. These findings warrant a deeper survey of variation in P. vivax to equip disease interventions targeting the distinctive biology of this neglected but major pathogen.

Drug resistant strains of the malaria parasite, Plasmodium falciparum, have rendered chloroquine ineffective throughout much of the world. In parts of Africa and Asia, the coordinated shift from chloroquine to other drugs has resulted in the near disappearance of chloroquine-resistant (CQR) parasites from the population. Currently, there is no molecular explanation for this phenomenon. Herein, we employ metabolic quantitative trait locus mapping (mQTL) to analyze progeny from a genetic cross between chloroquine-susceptible (CQS) and CQR parasites. We identify a family of hemoglobin-derived peptides that are elevated in CQR parasites and show that peptide accumulation, drug resistance, and reduced parasite fitness are all linked in vitro to CQR alleles of the P. falciparum chloroquine resistance transporter (pfcrt). These findings suggest that CQR parasites are less fit because mutations in pfcrt interfere with hemoglobin digestion by the parasite. Moreover, our findings may provide a molecular explanation for the reemergence of CQS parasites in wild populations.

Downy mildews are the most speciose group of oomycetes and affect crops of great economic importance. So far, there is only a single deeply-sequenced downy mildew genome available, from Hyaloperonospora arabidopsidis. Further genomic resources for downy mildews are required to study their evolution, including pathogenicity effector proteins, such as RxLR effectors. Plasmopara halstedii is a devastating pathogen of sunflower and a potential pathosystem model to study downy mildews, as several Avr-genes and R-genes have been predicted and unlike Arabidopsis downy mildew, large quantities of almost contamination-free material can be obtained easily.

Transcriptional responses to small molecules can provide insights into drug mode of action (MOA). The capacity of the human malaria parasite, Plasmodium falciparum, to respond specifically to transcriptional perturbations has been unclear based on past approaches. Here, we present the most extensive profiling to date of the parasite's transcriptional responsiveness to thirty-one chemically and functionally diverse small molecules.

Plasmodium parasites, the causative agents of malaria, have evolved a unique cell division cycle in the clinically relevant asexual blood stage of infection1. DNA replication commences approximately halfway through the intracellular development following invasion and parasite growth. The schizont stage is associated with multiple rounds of DNA replication and nuclear division without cytokinesis, resulting in a multinucleated cell. Nuclei divide asynchronously through schizogony, with only the final round of DNA replication and segregation being synchronous and coordinated with daughter cell assembly2,3. However, the control mechanisms for this divergent mode of replication are unknown. Here, we show that the Plasmodium-specific kinase PfCRK4 is a key cell-cycle regulator that orchestrates multiple rounds of DNA replication throughout schizogony in Plasmodium falciparum. PfCRK4 depletion led to a complete block in nuclear division and profoundly inhibited DNA replication. Quantitative phosphoproteomic profiling identified a set of PfCRK4-regulated phosphoproteins with greatest functional similarity to CDK2 substrates, particularly proteins involved in the origin of replication firing. PfCRK4 was required for initial and subsequent rounds of DNA replication during schizogony and, in addition, was essential for development in the mosquito vector. Our results identified an essential S-phase promoting factor of the unconventional P. falciparum cell cycle. PfCRK4 is required for both a prolonged period of the intraerythrocytic stage of Plasmodium infection, as well as for transmission, revealing a broad window for PfCRK4-targeted chemotherapeutics.

Competitive outcomes between co-infecting malaria parasite lines can reveal fitness disparities in blood stage growth. Blood stage fitness costs often accompany the evolution of drug resistance, with the expectation that relatively fitter parasites will be more likely to spread in populations. With the recent emergence of artemisinin resistance, it is important to understand the relative competitive fitness of the metabolically active asexual blood stage parasites. Genetically distinct drug resistant parasite clones with independently evolved sets of mutations are likely to vary in asexual proliferation rate, contributing to their chance of transmission to the mosquito vector.

Determining the genetic basis of fitness is central to understanding evolution and transmission of microbial pathogens. In human malaria parasites (Plasmodium falciparum), most experimental work on fitness has focused on asexual blood stage parasites, because this stage can be easily cultured, although the transmission of malaria requires both female Anopheles mosquitoes and vertebrate hosts. We explore a powerful approach to identify the genetic determinants of parasite fitness across both invertebrate and vertebrate life-cycle stages of P. falciparum. This combines experimental genetic crosses using humanized mice, with selective whole genome amplification and pooled sequencing to determine genome-wide allele frequencies and identify genomic regions under selection across multiple lifecycle stages. We applied this approach to genetic crosses between artemisinin resistant (ART-R, kelch13-C580Y) and ART-sensitive (ART-S, kelch13-WT) parasites, recently isolated from Southeast Asian patients. Two striking results emerge: we observed (i) a strong genome-wide skew (>80%) towards alleles from the ART-R parent in the mosquito stage, that dropped to ~50% in the blood stage as selfed ART-R parasites were selected against; and (ii) repeatable allele specific skews in blood stage parasites with particularly strong selection (selection coefficient (s) ≤ 0.18/asexual cycle) against alleles from the ART-R parent at loci on chromosome 12 containing MRP2 and chromosome 14 containing ARPS10. This approach robustly identifies selected loci and has strong potential for identifying parasite genes that interact with the mosquito vector or compensatory loci involved in drug resistance.

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.

Malaria liver stages represent an ideal therapeutic target with a bottleneck in parasite load and reduced clinical symptoms; however, current in vitro pre-erythrocytic (PE) models for Plasmodium vivax and P. falciparum lack the efficiency necessary for rapid identification and effective evaluation of new vaccines and drugs, especially targeting late liver-stage development and hypnozoites. Herein we report the development of a 384-well plate culture system using commercially available materials, including cryopreserved primary human hepatocytes. Hepatocyte physiology is maintained for at least 30 days and supports development of P. vivax hypnozoites and complete maturation of P. vivax and P. falciparum schizonts. Our multimodal analysis in antimalarial therapeutic research identifies important PE inhibition mechanisms: immune antibodies against sporozoite surface proteins functionally inhibit liver stage development and ion homeostasis is essential for schizont and hypnozoite viability. This model can be implemented in laboratories in disease-endemic areas to accelerate vaccine and drug discovery research.

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.

Improved control of Plasmodium vivax malaria can be achieved with the discovery of new antimalarials with radical cure efficacy, including prevention of relapse caused by hypnozoites residing in the liver of patients. We screened several compound libraries against P. vivax liver stages, including 1565 compounds against mature hypnozoites, resulting in one drug-like and several probe-like hits useful for investigating hypnozoite biology. Primaquine and tafenoquine, administered in combination with chloroquine, are currently the only FDA-approved antimalarials for radical cure, yet their activity against mature P. vivax hypnozoites has not yet been demonstrated in vitro. By developing an extended assay, we show both drugs are individually hypnozonticidal and made more potent when partnered with chloroquine, similar to clinically relevant combinations. Post-hoc analyses of screening data revealed excellent performance of ionophore controls and the high quality of single point assays, demonstrating a platform able to support screening of greater compound numbers. A comparison of P. vivax liver stage activity data with that of the P. cynomolgi blood, P. falciparum blood, and P. berghei liver stages reveals overlap in schizonticidal but not hypnozonticidal activity, indicating that the delivery of new radical curative agents killing P. vivax hypnozoites requires an independent and focused drug development test cascade.

The antimalarial activity of the frontline drug artemisinin involves generation of reactive oxygen species (ROS) leading to oxidative damage of parasite proteins. To achieve homeostasis and maintain protein quality control in the overwhelmed parasite, the ubiquitin-proteasome system kicks in. Even though molecular markers for artemisinin resistance like pfkelch13 have been identified, the intricate network of mechanisms driving resistance remains to be elucidated. Here, we report a forward genetic screening strategy that enables a broader identification of genetic factors responsible for altering sensitivity to dihydroartemisinin (DHA) and a proteasome inhibitor, bortezomib (BTZ). Using a library of isogenic piggyBac mutants in P. falciparum, we defined phenotype-genotype associations influencing drug responses and highlighted shared mechanisms between the two processes, which mainly included proteasome-mediated degradation and the lipid metabolism genes. Additional transcriptomic analysis of a DHA/BTZ-sensitive piggyBac mutant showed it is possible to find differences between the two response mechanisms on the specific components for regulation of the exportome. Our results provide further insight into the molecular mechanisms of antimalarial drug resistance. IMPORTANCE Malaria control is seriously threatened by the emergence and spread of Plasmodium falciparum resistance to the leading antimalarial, artemisinin. The potent killing activity of artemisinin results from oxidative damage unleashed by free heme activation released by hemoglobin digestion. Although the ubiquitin-proteasome system is considered critical for parasite survival of this toxicity, the diverse genetic changes linked to artemisinin resistance are complex and, so far, have not included the ubiquitin-proteasome system. In this study, we use a systematic forward genetic approach by screening a library of P. falciparum random piggyBac mutants to decipher the genetic factors driving malaria parasite responses to the oxidative stress caused by antimalarial drugs. This study compares phenotype-genotype associations influencing dihydroartemisinin responses with the proteasome inhibitor bortezomib to delineate the role of ubiquitin-proteasome system. Our study highlights shared and unique pathways from the complex array of molecular processes critical for P. falciparum survival resulting from the oxidative damage of artemisinin.

Fusarium oxysporum (Fo) belongs to a group of soil-borne hyphomycetes that are taxonomically collated in the Fusarium oxysporum Species Complex (FOSC). Hitherto, those infecting bananas were placed in the forma specialis cubense (Foc). Recently, however, these genetically different Foc lineages were recognized as new Fusarium spp. placed in the Fusarium of Banana Complex (FOBC). A member of this complex F. odoratissimum II-5 that uniquely comprises the so-called Tropical Race 4 (TR4), is a major problem sweeping through production zones of Cavendish banana in several regions of the world. Because of this, there is an urgent need for a phenotyping method that allows the screening for resistance to TR4 of large numbers of banana genotypes. Most Fusarium species produce three types of spores: macroconidia, microconidia and the persistent chlamydospores that can contaminate soils for many years. Inoculum production has been an important bottleneck for efficient phenotyping due to the low or variable number of conidia and the elaborate laboratory procedures requiring specific infrastructure. Here, we report a rapid, simple and high-yielding spore production method for nine F. oxysporum formae speciales as well as the biocontrol species Fo47 and Fo618-12. For Fusarium spp. causing Fusarium wilt or Panama disease of banana, we used the protocol for four species comprising the recognized physiological races, including Tropical Race 4 (TR4). We subsequently tested the produced inoculum in comparative inoculation trials on banana plants to evaluate their efficiency. All assays resulted in typical symptoms within 10 weeks; significant differences in final disease ratings were observed, depending on inoculum concentration. Pouring inoculum directly onto banana plants showed the most consistent and reproducible results, as expressed in external wilting, internal discoloration and determined by real-time PCR assays on entire rhizomes. Moreover, this method allows the inoculation of 250 plants per hour by one individual thereby facilitating the phenotyping of large mutant and breeding populations.

Oomycetes in the class Saprolegniomycetidae of the Eukaryotic kingdom Stramenopila have evolved as severe pathogens of amphibians, crustaceans, fish and insects, resulting in major losses in aquaculture and damage to aquatic ecosystems. We have sequenced the 63 Mb genome of the fresh water fish pathogen, Saprolegnia parasitica. Approximately 1/3 of the assembled genome exhibits loss of heterozygosity, indicating an efficient mechanism for revealing new variation. Comparison of S. parasitica with plant pathogenic oomycetes suggests that during evolution the host cellular environment has driven distinct patterns of gene expansion and loss in the genomes of plant and animal pathogens. S. parasitica possesses one of the largest repertoires of proteases (270) among eukaryotes that are deployed in waves at different points during infection as determined from RNA-Seq data. In contrast, despite being capable of living saprotrophically, parasitism has led to loss of inorganic nitrogen and sulfur assimilation pathways, strikingly similar to losses in obligate plant pathogenic oomycetes and fungi. The large gene families that are hallmarks of plant pathogenic oomycetes such as Phytophthora appear to be lacking in S. parasitica, including those encoding RXLR effectors, Crinkler's, and Necrosis Inducing-Like Proteins (NLP). S. parasitica also has a very large kinome of 543 kinases, 10% of which is induced upon infection. Moreover, S. parasitica encodes several genes typical of animals or animal-pathogens and lacking from other oomycetes, including disintegrins and galactose-binding lectins, whose expression and evolutionary origins implicate horizontal gene transfer in the evolution of animal pathogenesis in S. parasitica.

Albugo candida is a biotrophic oomycete that parasitizes various species of Brassicaceae, causing a disease (white blister rust) with remarkable convergence in behaviour to unrelated rusts of basidiomycete fungi.

Plasmodium falciparum exhibits resistance to the artemisinin component of the frontline antimalarial treatment Artemisinin-based Combination Therapy in South East Asia. Millions of lives will be at risk if artemisinin resistance (ART-R) spreads to Africa. Single non-synonymous mutations in the propeller region of PF3D7_1343700,"K13" are implicated in resistance. In this work, we use transcriptional profiling to characterize a laboratory-generated k13 insertional mutant previously demonstrated to have increased sensitivity to artemisinins to explore the functional role of k13.

The spread of Plasmodium falciparum multidrug resistance highlights the urgency to discover new targets and chemical scaffolds. Unfortunately, lack of experimentally validated functional information about most P. falciparum genes remains a strategic hurdle. Chemogenomic profiling is an established tool for classification of drugs with similar mechanisms of action by comparing drug fitness profiles in a collection of mutants. Inferences of drug mechanisms of action and targets can be obtained by associations between shifts in drug fitness and specific genetic changes in the mutants. In this screen, P. falciparum, piggyBac single insertion mutants were profiled for altered responses to antimalarial drugs and metabolic inhibitors to create chemogenomic profiles. Drugs targeting the same pathway shared similar response profiles and multiple pairwise correlations of the chemogenomic profiles revealed novel insights into drugs' mechanisms of action. A mutant of the artemisinin resistance candidate gene - "K13-propeller" gene (PF3D7_1343700) exhibited increased susceptibility to artemisinin drugs and identified a cluster of 7 mutants based on similar enhanced responses to the drugs tested. Our approach of chemogenomic profiling reveals artemisinin functional activity, linked by the unexpected drug-gene relationships of these mutants, to signal transduction and cell cycle regulation pathways.

Genetic screening using random transposon insertions has been a powerful tool for uncovering biology in prokaryotes, where whole-genome saturating screens have been performed in multiple organisms. In eukaryotes, such screens have proven more problematic, in part because of the lack of a sensitive and robust system for identifying transposon insertion sites. We here describe quantitative insertion-site sequencing, or QIseq, which uses custom library preparation and Illumina sequencing technology and is able to identify insertion sites from both the 5' and 3' ends of the transposon, providing an inbuilt level of validation. The approach was developed using piggyBac mutants in the human malaria parasite Plasmodium falciparum but should be applicable to many other eukaryotic genomes. QIseq proved accurate, confirming known sites in >100 mutants, and sensitive, identifying and monitoring sites over a >10,000-fold dynamic range of sequence counts. Applying QIseq to uncloned parasites shortly after transfections revealed multiple insertions in mixed populations and suggests that >4000 independent mutants could be generated from relatively modest scales of transfection, providing a clear pathway to genome-scale screens in P. falciparum QIseq was also used to monitor the growth of pools of previously cloned mutants and reproducibly differentiated between deleterious and neutral mutations in competitive growth. Among the mutants with fitness defects was a mutant with a piggyBac insertion immediately upstream of the kelch protein K13 gene associated with artemisinin resistance, implying mutants in this gene may have competitive fitness costs. QIseq has the potential to enable the scale-up of piggyBac-mediated genetics across multiple eukaryotic systems.