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.

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.

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.

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.

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.

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.

Sensing external signals and transducing these into intracellular responses requires a molecular signaling system that is crucial for every living organism. Two important eukaryotic signal transduction pathways that are often interlinked are G-protein signaling and phospholipid signaling. Heterotrimeric G-protein subunits activated by G-protein-coupled receptors (GPCRs) are typical stimulators of phospholipid signaling enzymes such as phosphatidylinositol phosphate kinases (PIPKs) or phospholipase C (PLC). However, a direct connection between the two pathways likely exists in oomycetes and slime molds, as they possess a unique class of GPCRs that have a PIPK as an accessory domain. In principle, these so-called GPCR-PIPKs have the capacity of perceiving an external signal (via the GPCR domain) that, via PIPK, directly activates downstream phospholipid signaling. Here we reveal the sporadic occurrence of GPCR-PIPKs in all eukaryotic supergroups, except for plants. Notably, all species having GPCR-PIPKs are unicellular microorganisms that favor aquatic environments. Phylogenetic analysis revealed that GPCR-PIPKs are likely ancestral to eukaryotes and significantly expanded in the last common ancestor of oomycetes. In addition to GPCR-PIPKs, we identified five hitherto-unknown classes of GPCRs with accessory domains, four of which are universal players in signal transduction. Similarly to GPCR-PIPKs, this enables a direct coupling between extracellular sensing and downstream signaling. Overall, our findings point to an ancestral signaling system in eukaryotes where GPCR-mediated sensing is directly linked to downstream responses.IMPORTANCE G-protein-coupled receptors (GPCRs) are central sensors that activate eukaryotic signaling and are the primary targets of human drugs. In this report, we provide evidence for the widespread though limited presence of a novel class of GPCRs in a variety of unicellular eukaryotes. These include free-living organisms and organisms that are pathogenic for plants, animals, and humans. The novel GPCRs have a C-terminal phospholipid kinase domain, pointing to a direct link between sensing external signals via GPCRs and downstream intracellular phospholipid signaling. Genes encoding these receptors were likely present in the last common eukaryotic ancestor and were lost during the evolution of higher eukaryotes. We further describe five other types of GPCRs with a catalytic accessory domain, the so-called GPCR-bigrams, four of which may potentially have a role in signaling. These findings shed new light onto signal transduction in microorganisms and provide evidence for alternative eukaryotic signaling pathways.

Banana is the most popular and most exported fruit and also a major food crop for millions of people around the world. Despite its importance and the presence of serious disease threats, research into this crop is limited. One of those is Panama disease or Fusarium wilt. In the previous century Fusarium wilt wiped out the "Gros Michel" based banana industry in Central America. The epidemic was eventually quenched by planting "Cavendish" bananas. However, 50 years ago the disease recurred, but now on "Cavendish" bananas. Since then the disease has spread across South-East Asia, to the Middle-East and the Indian subcontinent and leaped into Africa. Here, we report the presence of Fusarium oxysporum f.sp. cubense Tropical Race 4 (Foc TR4) in "Cavendish" plantations in Laos, Myanmar, and Vietnam. A combination of classical morphology, DNA sequencing, and phenotyping assays revealed a very close relationship between the Foc TR4 strains in the entire Greater Mekong Subregion (GMS), which is increasingly prone to intensive banana production. Analyses of single-nucleotide polymorphisms enabled us to initiate a phylogeography of Foc TR4 across three geographical areas-GMS, Indian subcontinent, and the Middle East revealing three distinct Foc TR4 sub-lineages. Collectively, our data place these new incursions in a broader agroecological context and underscore the need for awareness campaigns and the implementation of validated quarantine measures to prevent further international dissemination of Foc TR4.

Phytophthora species are oomycete plant pathogens with such major social and economic impact that genome sequences have been determined for Phytophthora infestans, P. sojae and P. ramorum. Pepsin-like aspartic proteinases (APs) are produced in a wide variety of species (from bacteria to humans) and contain conserved motifs and landmark residues. APs fulfil critical roles in infectious organisms and their host cells. Annotation of Phytophthora APs would provide invaluable information for studies into their roles in the physiology of Phytophthora species and interactions with their hosts.

We present an optimized probe design for copy number variation (CNV) and SNP genotyping in the Plasmodium falciparum genome. We demonstrate that variable length and isothermal probes are superior to static length probes. We show that sample preparation and hybridization conditions mitigate the effects of host DNA contamination in field samples. The microarray and workflow presented can be used to identify CNVs and SNPs with 95% accuracy in a single hybridization, in field samples containing up to 92% human DNA contamination.

Although maps of intracellular interactions are increasingly well characterized, little is known about large-scale maps of host-pathogen protein interactions. The investigation of host-pathogen interactions can reveal features of pathogenesis and provide a foundation for the development of drugs and disease prevention strategies. A compilation of experimentally verified interactions between HIV-1 and human proteins and a set of HIV-dependency factors (HDF) allowed insights into the topology and intricate interplay between viral and host proteins on a large scale. We found that targeted and HDF proteins appear predominantly in rich-clubs, groups of human proteins that are strongly intertwined among each other. These assemblies of proteins may serve as an infection gateway, allowing the virus to take control of the human host by reaching protein pathways and diversified cellular functions in a pronounced and focused way. Particular transcription factors and protein kinases facilitate indirect interactions between HDFs and viral proteins. Discerning the entanglement of directly targeted and indirectly interacting proteins may uncover molecular and functional sites that can provide novel perspectives on the progression of HIV infection and highlight new avenues to fight this virus.

Failure to establish an appropriate balance between pro- and anti-inflammatory immune responses is believed to contribute to pathogenesis of severe malaria. To determine whether this balance is maintained by classical regulatory T cells (CD4(+) FOXP3(+) CD127(-/low); Tregs) we compared cellular responses between Gambian children (n = 124) with severe Plasmodium falciparum malaria or uncomplicated malaria infections. Although no significant differences in Treg numbers or function were observed between the groups, Treg activity during acute disease was inversely correlated with malaria-specific memory responses detectable 28 days later. Thus, while Tregs may not regulate acute malarial inflammation, they may limit memory responses to levels that subsequently facilitate parasite clearance without causing immunopathology. Importantly, we identified a population of FOXP3(-), CD45RO(+) CD4(+) T cells which coproduce IL-10 and IFN-gamma. These cells are more prevalent in children with uncomplicated malaria than in those with severe disease, suggesting that they may be the regulators of acute malarial inflammation.

Multiple kelch13 alleles conferring artemisinin resistance (ART-R) are currently spreading through Southeast Asian malaria parasite populations, providing a unique opportunity to observe an ongoing soft selective sweep, investigate why resistance alleles have evolved multiple times and determine fundamental population genetic parameters for Plasmodium We sequenced kelch13 (n = 1,876), genotyped 75 flanking SNPs, and measured clearance rate (n = 3,552) in parasite infections from Western Thailand (2001-2014). We describe 32 independent coding mutations including common mutations outside the kelch13 propeller associated with significant reductions in clearance rate. Mutations were first observed in 2003 and rose to 90% by 2014, consistent with a selection coefficient of ∼0.079. ART-R allele diversity rose until 2012 and then dropped as one allele (C580Y) spread to high frequency. The frequency with which adaptive alleles arise is determined by the rate of mutation and the population size. Two factors drive this soft sweep: (1) multiple kelch13 amino-acid mutations confer resistance providing a large mutational target-we estimate the target is 87-163 bp. (2) The population mutation parameter (Θ = 2Neμ) can be estimated from the frequency distribution of ART-R alleles and is ∼5.69, suggesting that short term effective population size is 88 thousand to 1.2 million. This is 52-705 times greater than Ne estimated from fluctuation in allele frequencies, suggesting that we have previously underestimated the capacity for adaptive evolution in Plasmodium Our central conclusions are that retrospective studies may underestimate the complexity of selective events and the Ne relevant for adaptation for malaria is considerably higher than previously estimated.