The DEAD-box RNA helicase PfDOZI imposes opposing actions on RNA metabolism in Plasmodium falciparum

In malaria parasites, the regulation of mRNA translation, storage and degradation during development and life-stage transitions remains largely unknown. Here, we functionally characterized the DEAD-box RNA helicase PfDOZI in P. falciparum. Disruption of pfdozi enhanced asexual proliferation but reduced sexual commitment and impaired gametocyte development. By quantitative transcriptomics, we show that PfDOZI is involved in the regulation of invasion-related genes and sexual stage-specific genes during different developmental stages. PfDOZI predominantly participates in processing body-like mRNPs in schizonts but germ cell granule-like mRNPs in gametocytes to impose opposing actions of degradation and protection on different mRNA targets. We further show the formation of stress granule-like mRNPs during nutritional deprivation, highlighting an essential role of PfDOZI-associated mRNPs in stress response. We demonstrate that PfDOZI participates in distinct mRNPs to maintain mRNA homeostasis in response to life-stage transition and environmental changes by differentially executing post-transcriptional regulation on the target mRNAs.

MalariaSED: a deep learning framework to decipher the regulatory contributions of noncoding variants in malaria parasites

Malaria remains one of the deadliest infectious diseases. Transcriptional regulation effects of noncoding variants in this unusual genome of malaria parasites remain elusive. We developed a sequence-based, ab initio deep learning framework, MalariaSED, for predicting chromatin profiles in malaria parasites. The MalariaSED performance was validated by published ChIP-qPCR and TF motifs results. Applying MalariaSED to ~ 1.3 million variants shows that geographically differentiated noncoding variants are associated with parasite invasion and drug resistance. Further analysis reveals chromatin accessibility changes at Plasmodium falciparum rings are partly associated with artemisinin resistance. MalariaSED illuminates the potential functional roles of noncoding variants in malaria parasites.

A type II protein arginine methyltransferase regulates merozoite invasion in Plasmodium falciparum

Protein arginine methyltransferases (PRMTs) regulate many important cellular processes, such as transcription and RNA processing in model organisms but their functions in human malaria parasites are not elucidated. Here, we characterize PfPRMT5 in Plasmodium falciparum, which catalyzes symmetric dimethylation of histone H3 at R2 (H3R2me2s) and R8, and histone H4 at R3 in vitro. PfPRMT5 disruption results in asexual stage growth defects primarily due to lower invasion efficiency of the merozoites. Transcriptomic analysis reveals down-regulation of many transcripts related to invasion upon PfPRMT5 disruption, in agreement with H3R2me2s being an active chromatin mark. Genome-wide chromatin profiling detects extensive H3R2me2s marking of genes of different cellular processes, including invasion-related genes in wildtype parasites and PfPRMT5 disruption leads to the depletion of H3R2me2s. Interactome studies identify the association of PfPRMT5 with invasion-related transcriptional regulators such as AP2-I, BDP1, and GCN5. Furthermore, PfPRMT5 is associated with the RNA splicing machinery, and PfPRMT5 disruption caused substantial anomalies in RNA splicing events, including those for invasion-related genes. In summary, PfPRMT5 is critical for regulating parasite invasion and RNA splicing in this early-branching eukaryote.

Chemogenomic Profiling of a Plasmodium falciparum Transposon Mutant Library Reveals Shared Effects of Dihydroartemisinin and Bortezomib on Lipid Metabolism and Exported Proteins

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.

Large Mammals Have More Powerful Antibacterial Defenses Than Expected from Their Metabolic Rates

AbstractTerrestrial mammals span seven orders of magnitude in body size, ranging from the <2-g Etruscan pygmy shrew (Suncus etruscus) to the >3,900-kg African elephant (Loxodonta africana). Although body size profoundly affects the behavior, physiology, ecology, and evolution of species, how investment in functional immune defenses changes with body size across species is unknown. Here, we (1) developed a novel 12-point dilution curve approach to describe and compare antibacterial capacity against three bacterial species among >160 terrestrial species of mammals and (2) tested published predictions about the scaling of immune defenses. Our study focused on the safety factor hypothesis, which predicts that broad, early-acting immune defenses should scale hypermetrically with body mass. However, our three statistical approaches demonstrated that antibacterial activity in sera across mammals exhibits isometry; killing capacity did not change with body size across species. Intriguingly, this result indicates that the serum of a large mammal is less hospitable to bacteria than would be predicted by its metabolic rates. In other words, if metabolic rates underlie the rates of physiological reactions as postulated by the metabolic theory of ecology, large species should have disproportionately lower antibacterial capacity than small species, but they do not. These results have direct implications for effectively modeling the evolution of immune defenses and identifying potential reservoir hosts of pathogens.

A unique class of Zn2+-binding serine-based PBPs underlies cephalosporin resistance and sporogenesis in Clostridioides difficile

Treatment with β-lactam antibiotics, particularly cephalosporins, is a major risk factor for Clostridioides difficile infection. These broad-spectrum antibiotics irreversibly inhibit penicillin-binding proteins (PBPs), which are serine-based enzymes that assemble the bacterial cell wall. However, C. difficile has four different PBPs (PBP1-3 and SpoVD) with various roles in growth and spore formation, and their specific links to β-lactam resistance in this pathogen are underexplored. Here, we show that PBP2 (known to be essential for vegetative growth) is the primary bactericidal target for β-lactams in C. difficile. PBP2 is insensitive to cephalosporin inhibition, and this appears to be the main basis for cephalosporin resistance in this organism. We determine crystal structures of C. difficile PBP2, alone and in complex with β-lactams, revealing unique features including ligand-induced conformational changes and an active site Zn²⁺-binding motif that influences β-lactam binding and protein stability. The Zn²⁺-binding motif is also present in C. difficile PBP3 and SpoVD (which are known to be essential for sporulation), as well as in other bacterial taxa including species living in extreme environments and the human gut. We speculate that this thiol-containing motif and its cognate Zn²⁺ might function as a redox sensor to regulate cell wall synthesis for survival in adverse or anaerobic environments.

Probing the distinct chemosensitivity of Plasmodium vivax liver stage parasites and demonstration of 8-aminoquinoline radical cure activity in vitro

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.

Light pollution affects West Nile virus exposure risk across Florida

Emerging infectious diseases (EIDs) present global health threats, and their emergences are often linked to anthropogenic change. Artificial light at night (ALAN) is one form of anthropogenic change that spans beyond urban boundaries and may be relevant to EIDs through its influence on the behaviour and physiology of hosts and/or vectors. Although West Nile virus (WNV) emergence has been described as peri-urban, we hypothesized that exposure risk could also be influenced by ALAN in particular, which is testable by comparing the effects of ALAN on prevalence while controlling for other aspects of urbanization. By modelling WNV exposure among sentinel chickens in Florida, we found strong support for a nonlinear relationship between ALAN and WNV exposure risk in chickens with peak WNV risk occurring at low ALAN levels. Although our goal was not to discern how ALAN affected WNV relative to other factors, effects of ALAN on WNV exposure were stronger than other known drivers of risk (i.e. impervious surface, human population density). Ambient temperature in the month prior to sampling, but no other considered variables, strongly influenced WNV risk. These results indicate that ALAN may contribute to spatio-temporal changes in WNV risk, justifying future investigations of ALAN on other vector-borne parasites.

False Alarms in Consumer Genomics Add to Public Fear and Potential Health Care Burden

We have entered an era of direct-to-consumer (DTC) genomics. Patients have relayed many success stories of DTC genomics about finding causal mutations of genetic diseases before showing any symptoms and taking precautions. However, consumers may also take unnecessary medical actions based on false alarms of “pathogenic alleles”. The severity of this problem is not well known. Using publicly available data, we compared DTC microarray genotyping data with deep-sequencing data of 5 individuals and manually checked each inconsistently reported single nucleotide variants (SNVs). We estimated that, on average, a person would have ~5 “pathogenic” alleles reported due to wrongly reported genotypes if using a 23andMe genotyping microarray. We also found that the number of wrongly classified “pathogenic” alleles per person is at least as significant as those due to wrongly reported genotypes. We show that the scale of the false alarm problem could be large enough that the medical costs will become a burden to public health.

The human malaria parasite genome is configured into thousands of coexpressed linear regulatory units

The human malaria parasite Plasmodium falciparum thrives in radically different host environments in mosquitoes and humans, with only a limited set of transcription factors. The nature of regulatory elements or their target genes in the P. falciparum genome remains elusive. Here, we found that this eukaryotic parasite uses an efficient way to maximally use genetic and epigenetic regulation to form regulatory units (RUs) during blood infections. Genes located in the same RU tend to have the same pattern of expression over time and are associated with open chromatin along regulatory elements. To precisely define and quantify these RUs, a novel hidden Markov model was developed to capture the regulatory structure in a genome-wide fashion by integrating expression and epigenetic evidence. We successfully identified thousands of RUs and cross-validated with previous findings. We found more genes involved in red blood cell (RBC) invasion located in the same RU as the PfAP2-I (AP2-I) transcription factor, demonstrating that AP2-I is responsible for regulating RBC invasion. Our study has provided a regulatory mechanism for a compact eukaryotic genome and offers new insights into the in vivo transcriptional regulation of the P. falciparum intraerythrocytic stage.

An adaptable soft-mold embossing process for fabricating optically-accessible, microfeature-based culture systems and application toward liver stage antimalarial compound testing

Advanced cell culture methods for modeling organ-level structure have been demonstrated to replicate in vivo conditions more accurately than traditional in vitro cell culture. Given that the liver is particularly important to human health, several advanced culture methods have been developed to experiment with liver disease states, including infection with Plasmodium parasites, the causative agent of malaria. These models have demonstrated that intrahepatic parasites require functionally stable hepatocytes to thrive and robust characterization of the parasite populations' response to investigational therapies is dependent on high-content and high-resolution imaging (HC/RI). We previously reported abiotic confinement extends the functional longevity of primary hepatocytes in a microfluidic platform and set out to instill confinement in a microtiter plate platform while maintaining optical accessibility for HC/RI; with an end-goal of producing an improved P. vivax liver stage culture model. We developed a novel fabrication process in which a PDMS soft mold embosses hepatocyte-confining microfeatures into polystyrene, resulting in microfeature-based hepatocyte confinement (μHEP) slides and plates. Our process was optimized to form both microfeatures and culture wells in a single embossing step, resulting in a 100 μm-thick bottom ideal for HC/RI, and was found inexpensively amendable to microfeature design changes. Microfeatures improved intrahepatic parasite infection rates and μHEP systems were used to reconfirm the activity of reference antimalarials in phenotypic dose-response assays. RNAseq of hepatocytes in μHEP systems demonstrated microfeatures sustain hepatic differentiation and function, suggesting broader utility for preclinical hepatic assays; while our tailorable embossing process could be repurposed for developing additional organ models.

5-Aminolevulinate synthase catalysis: The catcher in heme biosynthesis

5-Aminolevulinate (ALA) synthase (ALAS), a homodimeric pyridoxal-5'-phosphate (PLP)-dependent enzyme, catalyzes the first step of heme biosynthesis in metazoa, fungi and α-proteobacteria. In this review, we focus on the advances made in unraveling the mechanism of the ALAS-catalyzed reaction during the past decade. The interplay between the PLP cofactor and the protein moiety determines and modulates the multi-intermediate reaction cycle of ALAS, which involves the decarboxylative condensation of two substrates, glycine and succinyl-CoA. Substrate binding and catalysis are rapid, and product (ALA) release dominates the overall ALAS kinetic mechanism. Interconversion between a catalytically incompetent, open conformation and a catalytically competent, closed conformation is linked to ALAS catalysis. Reversion to the open conformation, coincident with ALA dissociation, defines the slowest step of the reaction cycle. These findings were further substantiated by introducing seven mutations in the16-amino acid loop that gates the active site, yielding an ALAS variant with a greatly increased rate of catalytic turnover and heightened specificity constants for both substrates. Recently, molecular dynamics (MD) simulation analysis of various dimeric ALAS forms revealed that the seven active site loop mutations caused the proteins to adopt different conformations. In particular, the emergence of a β-strand in the mutated loop, which interacted with two preexisting β-strands to form an anti-parallel three-stranded β-sheet, conferred the murine heptavariant with a more stable open conformation and prompted faster product release than wild-type mALAS2. Moreover, the dynamics of the mALAS2 active site loop anti-correlated with that of the 35 amino acid C-terminal sequence. This led us to propose that this C-terminal extension, which is absent in prokaryotic ALASs, finely tunes mammalian ALAS activity. Based on the above results, we extend our previous proposal to include that discovery of a ligand inducing the mammalian C-terminal extension to fold offers a good prospect for the development of a new drug for X-linked protoporphyria and/or other porphyrias associated with enhanced ALAS activity and/or porphyrin accumulation.

Light pollution increases West Nile virus competence of a ubiquitous passerine reservoir species

Among the many anthropogenic changes that impact humans and wildlife, one of the most pervasive but least understood is light pollution. Although detrimental physiological and behavioural effects resulting from exposure to light at night are widely appreciated, the impacts of light pollution on infectious disease risk have not been studied. Here, we demonstrate that artificial light at night (ALAN) extends the infectious-to-vector period of the house sparrow (Passer domesticus), an urban-dwelling avian reservoir host of West Nile virus (WNV). Sparrows exposed to ALAN maintained transmissible viral titres for 2 days longer than controls but did not experience greater WNV-induced mortality during this window. Transcriptionally, ALAN altered the expression of gene regulatory networks including key hubs (OASL, PLBD1 and TRAP1) and effector genes known to affect WNV dissemination (SOCS). Despite mounting anti-viral immune responses earlier, transcriptomic signatures indicated that ALAN-exposed individuals probably experienced pathogen-induced damage and immunopathology, potentially due to evasion of immune effectors. A simple mathematical modelling exercise indicated that ALAN-induced increases of host infectious-to-vector period could increase WNV outbreak potential by approximately 41%. ALAN probably affects other host and vector traits relevant to transmission, and additional research is needed to advise the management of zoonotic diseases in light-polluted areas.

Iron Hack - A symposium/hackathon focused on porphyrias, Friedreich's ataxia, and other rare iron-related diseases

Background: Basic and clinical scientific research at the University of South Florida (USF) have intersected to support a multi-faceted approach around a common focus on rare iron-related diseases. We proposed a modified version of the National Center for Biotechnology Information’s (NCBI) Hackathon-model to take full advantage of local expertise in building “Iron Hack”, a rare disease-focused hackathon. As the collaborative, problem-solving nature of hackathons tends to attract participants of highly-diverse backgrounds, organizers facilitated a symposium on rare iron-related diseases, specifically porphyrias and Friedreich’s ataxia, pitched at general audiences.

Methods: The hackathon was structured to begin each day with presentations by expert clinicians, genetic counselors, researchers focused on molecular and cellular biology, public health/global health, genetics/genomics, computational biology, bioinformatics, biomolecular science, bioengineering, and computer science, as well as guest speakers from the American Porphyria Foundation (APF) and Friedreich’s Ataxia Research Alliance (FARA) to inform participants as to the human impact of these diseases.

Results: As a result of this hackathon, we developed resources that are relevant not only to these specific disease-models, but also to other rare diseases and general bioinformatics problems. Within two and a half days, “Iron Hack” participants successfully built collaborative projects to visualize data, build databases, improve rare disease diagnosis, and study rare-disease inheritance.

Conclusions: The purpose of this manuscript is to demonstrate the utility of a hackathon model to generate prototypes of generalizable tools for a given disease and train clinicians and data scientists to interact more effectively.

A SLC6 transporter cloned from the lion's mane jellyfish (Cnidaria, Scyphozoa) is expressed in neurons

In the course of recent comparative genomic studies conducted on nervous systems across the phylogeny, current thinking is leaning in favor of more heterogeneity among nervous systems than what was initially expected. The isolation and characterization of molecular components that constitute the cnidarian neuron is not only of interest to the physiologist but also, on a larger scale, to those who study the evolution of nervous systems. Understanding the function of those ancient neurons involves the identification of neurotransmitters and their precursors, the description of nutrients used by neurons for metabolic purposes and the identification of integral membrane proteins that bind to those compounds. Using a molecular cloning strategy targeting membrane proteins that are known to be present in all forms of life, we isolated a member of the solute carrier family 6 from the scyphozoan jellyfish Cyanea capillata. The phylogenetic analysis suggested that the new transporter sequence belongs to an ancestral group of the nutrient amino acid transporter subfamily and is part of a cluster of cnidarian sequences which may translocate the same substrate. We found that the jellyfish transporter is expressed in neurons of the motor nerve net of the animal. To this end, we established an in situ hybridization protocol for the tissues of C. capillata and developed a specific antibody to the jellyfish transporter. Finally, we showed that the gene that codes for the jellyfish transporter also expresses a long non-coding RNA. We hope that this research will contribute to studies that seek to understand what constitutes a neuron in species that belong to an ancient phylum.

Plasmodium vivax readiness to transmit: implication for malaria eradication

Background

The lack of a continuous long-term in vitro culture system for Plasmodium vivax severely limits our knowledge of pathophysiology of the most widespread malaria parasite. To gain direct understanding of P. vivax human infections, we used Next Generation Sequencing data mining to unravel parasite in vivo expression profiles for P. vivax, and P. falciparum as comparison.

Results

We performed cloud and local computing to extract parasite transcriptomes from publicly available raw data of human blood samples. We developed a Poisson Modelling (PM) method to confidently identify parasite derived transcripts in mixed RNAseq signals of infected host tissues. We successfully retrieved and reconstructed parasite transcriptomes from infected patient blood as early as the first blood stage cycle; and the same methodology did not recover any significant signal from controls. Surprisingly, these first generation blood parasites already show strong signature of transmission, which indicates the commitment from asexual-to-sexual stages. Further, we place the results within the context of P. vivax’s complex life cycle, by developing mathematical models for P. vivax and P. falciparum and using sensitivity analysis assess the relative epidemiological impact of possible early stage transmission.

Conclusion

The study uncovers the earliest onset of P. vivax blood pathogenesis and highlights the challenges of P. vivax eradication programs.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0669-4) contains supplementary material, which is available to authorized users.

Altered expression of K13 disrupts DNA replication and repair in Plasmodium falciparum

BACKGROUND

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.

RESULTS

A set of RNA-seq and microarray experiments confirmed that the expression profile of k13 is specifically altered during the early ring and early trophozoite stages of the mutant intraerythrocytic development cycle. The down-regulation of k13 transcripts in this mutant during the early ring stage is associated with a transcriptome advance towards a more trophozoite-like state. To discover the specific downstream effect of k13 dysregulation, we developed a new computational method to search for differential gene expression while accounting for the temporal sequence of transcription. We found that the strongest biological signature of the transcriptome shift is an up-regulation of DNA replication and repair genes during the early ring developmental stage and a down-regulation of DNA replication and repair genes during the early trophozoite stage; by contrast, the expressions of housekeeping genes are unchanged. This effect, due to k13 dysregulation, is antagonistic, such that k13 levels are negatively correlated with DNA replication and repair gene expression.

CONCLUSION

Our results support a role for k13 as a stress response regulator consistent with the hypothesis that artemisinins mode of action is oxidative stress and k13 as a functional homolog of Keap1 which in humans regulates DNA replication and repair genes in response to oxidative stress.

Unraveling the Plasmodium vivax sporozoite transcriptional journey from mosquito vector to human host

Malaria parasites transmitted by mosquito bite are remarkably efficient in establishing human infections. The infection process requires roughly 30 minutes and is highly complex as quiescent sporozoites injected with mosquito saliva must be rapidly activated in the skin, migrate through the body, and infect the liver. This process is poorly understood for Plasmodium vivax due to low infectivity in the in vitro models. To study this skin-to-liver-stage of malaria, we used quantitative bioassays coupled with transcriptomics to evaluate parasite changes linked with mammalian microenvironmental factors. Our in vitro phenotyping and RNA-seq analyses revealed key microenvironmental relationships with distinct biological functions. Most notable, preservation of sporozoite quiescence by exposure to insect-like factors coupled with strategic activation limits untimely activation of invasion-associated genes to dramatically increase hepatocyte invasion rates. We also report the first transcriptomic analysis of the P. vivax sporozoite interaction in salivary glands identifying 118 infection-related differentially-regulated Anopheles dirus genes. These results provide important new insights in malaria parasite biology and identify priority targets for antimalarial therapeutic interventions to block P. vivax infection.

Uncovering the essential genes of the human malaria parasite Plasmodium falciparum by saturation mutagenesis

Severe malaria is caused by the apicomplexan parasite Plasmodium falciparum. Despite decades of research, the distinct biology of these parasites has made it challenging to establish high-throughput genetic approaches to identify and prioritize therapeutic targets. Using transposon mutagenesis of P. falciparum in an approach that exploited its AT-rich genome, we generated more than 38,000 mutants, saturating the genome and defining mutability and fitness costs for over 87% of genes. Of 5399 genes, our study defined 2680 genes as essential for optimal growth of asexual blood stages in vitro. These essential genes are associated with drug resistance, represent leading vaccine candidates, and include approximately 1000 Plasmodium-conserved genes of unknown function. We validated this approach by testing proteasome pathways for individual mutants associated with artemisinin sensitivity.

In-depth phenotypic characterization of reticulocyte maturation using mass cytometry

Progress towards an in-depth understanding of the final steps of the erythroid lineage development is paramount for many hematological diseases. We have characterized the final stages of reticulocyte maturation from bone marrow to peripheral blood using for the first time single-cell Mass Cytometry (CyTOF). We were able to measure the expression of 31 surface markers within a single red blood cell (RBC). We demonstrate the validity of CyTOF for RBC phenotyping by confirming the progressive reduction of transferrin receptor 1 (CD71) during reticulocyte maturation to mature RBC. We highlight the high-dimensional nature of mass cytometry data by correlating the expression of multiple proteins on individual RBCs. We further describe a more drastic reduction pattern for a component of the alpha4/beta1 integrin CD49d at the very early steps of reticulocyte maturation in bone marrow and directly linked with the mitochondria remnants clearance pattern. The enhanced and accurate RBC phenotyping potential of CyTOF described herein could be beneficial to decipher RBC preferences, as well as still not well understood receptor-ligand interaction of some hemotropic parasites such as the malaria causing agent Plasmodium vivax.

A comprehensive model for assessment of liver stage therapies targeting Plasmodium vivax and Plasmodium falciparum

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.

Malaria infected red blood cells release small regulatory RNAs through extracellular vesicles

The parasite Plasmodium falciparum causes the most severe form of malaria. Cell communication between parasites is an important mechanism to control population density and differentiation. The infected red blood cells (iRBCs) release small extracellular vesicles (EVs) that transfer cargoes between cells. The EVs synchronize the differentiation of the asexual parasites into gametocytes to initiate the transmission to the mosquito. Beside their role in parasite communication, EVs regulate vascular function. So far, the exact cargoes responsible for cellular communication remain unknown. We isolated EVs from cultured iRBCs to determine their small RNA content. We identified several types of human and plasmodial regulatory RNAs. While the miRNAs and tRNA-derived fragments were the most abundant human RNAs, we also found Y-RNAs, vault RNAs, snoRNAs and piRNAs. Interestingly, we found about 120 plasmodial RNAs, including mRNAs coding for exported proteins and proteins involved in drug resistance, as well as non-coding RNAs, such as rRNAs, small nuclear (snRNAs) and tRNAs. These data show, that iRBC-EVs carry small regulatory RNAs. A role in cellular communication is possible since the RNAs were transferred to endothelial cells. Furthermore, the presence of Plasmodium RNAs, in EVs suggests that they may be used as biomarker to track and detect disease.

Lysophosphatidylcholine Regulates Sexual Stage Differentiation in the Human Malaria Parasite Plasmodium falciparum

Transmission represents a population bottleneck in the Plasmodium life cycle and a key intervention target of ongoing efforts to eradicate malaria. Sexual differentiation is essential for this process, as only sexual parasites, called gametocytes, are infective to the mosquito vector. Gametocyte production rates vary depending on environmental conditions, but external stimuli remain obscure. Here, we show that the host-derived lipid lysophosphatidylcholine (LysoPC) controls P. falciparum cell fate by repressing parasite sexual differentiation. We demonstrate that exogenous LysoPC drives biosynthesis of the essential membrane component phosphatidylcholine. LysoPC restriction induces a compensatory response, linking parasite metabolism to the activation of sexual-stage-specific transcription and gametocyte formation. Our results reveal that malaria parasites can sense and process host-derived physiological signals to regulate differentiation. These data close a critical knowledge gap in parasite biology and introduce a major component of the sexual differentiation pathway in Plasmodium that may provide new approaches for blocking malaria transmission.

Plasmodium falciparum CRK4 directs continuous rounds of DNA replication during schizogony

Plasmodium parasites, the causative agents of malaria, have evolved a unique cell division cycle in the clinically relevant asexual blood-stage of infection¹. 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 assembly^2,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 the multiple rounds of DNA replication throughout schizogony in P. 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 origin of replication firing. PfCRK4 was required for the 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 blood-stage of malaria infection, as well as for transmission, revealing a broad window for PfCRK4-targeted chemotherapeutics.

Punctuated chromatin states regulate Plasmodium falciparum antigenic variation at the intron and 2 kb upstream regions

Background

Understanding the regulation mechanism of var gene expression is crucial for explaining antigenic variation in Plasmodium falciparum. Recent work observed that while all var genes produce transcripts, only a few var genes exhibit high expression levels. However, the global regulation of var expression and the relationship between epigenetic and genetic control remains to be established.

Result

We have systematically reanalyzed the existing genomic data including chromatin configurations and gene expressions; and for the first time used robust statistical methods to show that the intron and 2 kb upstream regions of each endogenous var gene always maintain high chromatin accessibility, with high potential to bind transcription factors (TFs). The levels of transcripts for different var gene family members are associated with this chromatin accessibility. Any given var gene thus shows punctuated chromatin states throughout the asexual life cycle. This is demonstrated by three independent transcript datasets. Chromatin accessibility in the var intron and 2 kb upstream regions are also positively correlated with their GC content, suggesting the level of var genes silencing might be encoded in their intron sequences. Interestingly, both var intron and 2 kb upstream regions exhibit higher chromatin accessibility when the genes have relatively lower transcription levels, suggesting a punctuated repressive function for these regulatory elements.

Conclusion

By integrating and analyzing epigenomic, genomic and transcriptomic data, our work reveals a novel distal element in var control. We found dynamic modulations of specific epigenetic marks around the var intron and distal upstream regions are involved in the general var gene expression patterns in malarial antigenic variation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3005-7) contains supplementary material, which is available to authorized users.

Quantitative insertion-site sequencing (QIseq) for high throughput phenotyping of transposon mutants

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.

Ancient human sialic acid variant restricts an emerging zoonotic malaria parasite

Plasmodium knowlesi is a zoonotic parasite transmitted from macaques causing malaria in humans in Southeast Asia. Plasmodium parasites bind to red blood cell (RBC) surface receptors, many of which are sialylated. While macaques synthesize the sialic acid variant N-glycolylneuraminic acid (Neu5Gc), humans cannot because of a mutation in the enzyme CMAH that converts N-acetylneuraminic acid (Neu5Ac) to Neu5Gc. Here we reconstitute CMAH in human RBCs for the reintroduction of Neu5Gc, which results in enhancement of P. knowlesi invasion. We show that two P. knowlesi invasion ligands, PkDBPβ and PkDBPγ, bind specifically to Neu5Gc-containing receptors. A human-adapted P. knowlesi line invades human RBCs independently of Neu5Gc, with duplication of the sialic acid-independent invasion ligand, PkDBPα and loss of PkDBPγ. Our results suggest that absence of Neu5Gc on human RBCs limits P. knowlesi invasion, but that parasites may evolve to invade human RBCs through the use of sialic acid-independent pathways.

Phenotypic and molecular analysis of the effect of 20-hydroxyecdysone on the human filarial parasite Brugia malayi

A homologue of the ecdysone receptor has been identified and shown to be responsive to 20-hydroxyecdysone in Brugia malayi. However, the role of this master regulator of insect development has not been delineated in filarial nematodes. Gravid adult female B. malayi cultured in the presence of 20-hydroxyecdysone produced significantly more microfilariae and abortive immature progeny than control worms, implicating the ecdysone receptor in regulation of embryogenesis and microfilarial development. Transcriptome analyses identified 30 genes whose expression was significantly up-regulated in 20-hydroxyecdysone-treated parasites compared with untreated controls. Of these, 18% were identified to be regulating transcription. A comparative proteomic analysis revealed 932 proteins to be present in greater amounts in extracts of 20-hydroxyecdysone-treated adult females than in extracts prepared from worms cultured in the absence of the hormone. Of the proteins exhibiting a greater than two-fold difference in the 20-hydroxyecdysone-treated versus untreated parasite extracts, 16% were involved in transcriptional regulation. RNA interference (RNAi) phenotype analysis of Caenorhabditis elegans orthologs revealed that phenotypes involved in developmental processes associated with embryogenesis were significantly over-represented in the transcripts and proteins that were up-regulated by exposure to 20-hydroxyecdysone. Taken together, the transcriptomic, proteomic and phenotypic data suggest that the filarial ecdysone receptor may play a role analogous to that in insects, where it serves as a regulator of egg development.

Shared elements of host-targeting pathways among apicomplexan parasites of differing lifestyles

Apicomplexans are a diverse group of obligate parasites occupying different intracellular niches that require modification to meet the needs of the parasite. To efficiently manipulate their environment, apicomplexans translocate numerous parasite proteins into the host cell. Whereas some parasites remain contained within a parasitophorous vacuole membrane (PVM) throughout their developmental cycle, others do not, a difference that affects the machinery needed for protein export. A signal-mediated pathway for protein export into the host cell has been characterized in Plasmodium parasites, which maintain the PVM. Here, we functionally demonstrate an analogous host-targeting pathway involving organellar staging prior to secretion in the related bovine parasite, Babesia bovis, a parasite that destroys the PVM shortly after invasion. Taking into account recent identification of a similar signal-mediated pathway in the coccidian parasite Toxoplasma gondii, we suggest a model in which this conserved pathway has evolved in multiple steps from signal-mediated trafficking to specific secretory organelles for controlled secretion to a complex protein translocation process across the PVM.

Parasite Calcineurin Regulates Host Cell Recognition and Attachment by Apicomplexans

Apicomplexans invade a variety of metazoan host cells through mechanisms involving host cell receptor engagement and secretion of parasite factors to facilitate cellular attachment. We find that the parasite homolog of calcineurin, a calcium-regulated phosphatase complex central to signal transduction in eukaryotes, also contributes to host cell invasion by the malaria parasite Plasmodium falciparum and related Toxoplasma gondii. Using reverse-genetic and chemical-genetic approaches, we determine that calcineurin critically regulates and stabilizes attachment of extracellular P. falciparum to host erythrocytes before intracellular entry and has similar functions in host cell engagement by T. gondii. Calcineurin-mediated Plasmodium invasion is strongly associated with host receptors required for host cell recognition, and calcineurin function distinguishes this form of receptor-mediated attachment from a second mode of host-parasite adhesion independent of host receptors. This specific role of calcineurin in coordinating physical interactions with host cells highlights an ancestral mechanism for parasitism used by apicomplexans.

Malaria. A forward genetic screen identifies erythrocyte CD55 as essential for Plasmodium falciparum invasion

Efforts to identify host determinants for malaria have been hindered by the absence of a nucleus in erythrocytes, which precludes genetic manipulation in the cell in which the parasite replicates. We used cultured red blood cells derived from hematopoietic stem cells to carry out a forward genetic screen for Plasmodium falciparum host determinants. We found that CD55 is an essential host factor for P. falciparum invasion. CD55-null erythrocytes were refractory to invasion by all isolates of P. falciparum because parasites failed to attach properly to the erythrocyte surface. Thus, CD55 is an attractive target for the development of malaria therapeutics. Hematopoietic stem cell-based forward genetic screens may be valuable for the identification of additional host determinants of malaria pathogenesis.

A Plasmodium falciparum histone deacetylase regulates antigenic variation and gametocyte conversion

The asexual forms of the malaria parasite Plasmodium falciparum are adapted for chronic persistence in human red blood cells, continuously evading host immunity using epigenetically regulated antigenic variation of virulence-associated genes. Parasite survival on a population level also requires differentiation into sexual forms, an obligatory step for further human transmission. We reveal that the essential nuclear gene, P. falciparum histone deacetylase 2 (PfHda2), is a global silencer of virulence gene expression and controls the frequency of switching from the asexual cycle to sexual development. PfHda2 depletion leads to dysregulated expression of both virulence-associated var genes and PfAP2-g, a transcription factor controlling sexual conversion, and is accompanied by increases in gametocytogenesis. Mathematical modeling further indicates that PfHda2 has likely evolved to optimize the parasite's infectious period by achieving low frequencies of virulence gene expression switching and sexual conversion. This common regulation of cellular transcriptional programs mechanistically links parasite transmissibility and virulence.

Functional analysis of sirtuin genes in multiple Plasmodium falciparum strains

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.

Distinctive expansion of potential virulence genes in the genome of the oomycete fish pathogen Saprolegnia parasitica

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.

Fish are an increasingly important source of animal protein globally, with aquaculture production rising dramatically over the past decade. Saprolegnia is a fungal-like oomycete and one of the most destructive fish pathogens, causing millions of dollars in losses to the aquaculture industry annually. Saprolegnia has also been linked to a worldwide decline in wild fish and amphibian populations. Here we describe the genome sequence of the first animal pathogenic oomycete and compare the genome content with the available plant pathogenic oomycetes. We found that Saprolegnia lacks the large effector families that are hallmarks of plant pathogenic oomycetes, showing evolutionary adaptation to the host. Moreover, Saprolegnia harbors pathogenesis-related genes that were derived by lateral gene transfer from the host and other animal pathogens. The retrotransposon LINE family also appears to be acquired from animal lineages. By transcriptome analysis we show a high rate of allelic variation, which reveals rapidly evolving genes and potentially adaptive evolutionary mechanisms coupled to selective pressures exerted by the animal host. The genome and transcriptome data, as well as subsequent biochemical analyses, provided us with insight in the disease process of Saprolegnia at a molecular and cellular level, providing us with targets for sustainable control of Saprolegnia.

Eukaryotic virulence determinants utilize phosphoinositides at the ER and host cell surface

Similar to bacteria, eukaryotic pathogens may utilize common strategies of pathogenic secretion, because effector proteins from the oomycete Phytophthora infestans and virulence determinants from the human malaria parasite Plasmodium falciparum share a functionally equivalent host-cell-targeting motif (RxLR-dEER in P. infestans and RxLxE/D/Q in P. falciparum). Here we summarize recent studies that reveal that the malarial motif may function differently than previously envisioned. Binding of the lipid phosphatidylinositol 3-phosphate [PI(3)P] is a critical step in accessing the host for both pathogens, but occurs in different locations. Nanomolar affinity for PI(3)P by these short amino acid motifs suggests that a newly identified mechanism of phosphoinositide binding that unexpectedly occurs in secretory locations has been exploited for virulence by diverse eukaryotic pathogens.

The malaria parasite Plasmodium vivax exhibits greater genetic diversity than Plasmodium falciparum

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.

A PIP Gets the plasmodium protein export pathway going

Survival of blood stage malaria parasites requires extensive host cell remodeling, which is facilitated by secretion of parasite proteins via a dedicated protein export pathway. In a recent Cell paper, Bhattacharjee et al., (2012) describe PI(3)P binding as one of the first steps in targeting parasite proteins to the host cell.

Massively parallel sequencing technology in pathogenic microbes

Next-Generation Sequencing (NGS) methods have revolutionized various aspects of genomics including transcriptome analysis. Digital expression analysis is all set to replace analog expression analysis that uses microarray chips through their cost-effectiveness, reproducibility, accuracy, and speed. The last 2 years have seen a surge in the development of statistical methods and software tools for analysis and visualization of NGS data. Large amounts of NGS data are available for pathogenic fungi and oomycetes. As the analysis results start pouring in, it brings about a paradigm shift in the understanding of host pathogen interactions with discovery of new transcripts, splice variants, mutations, regulatory elements, and epigenetic controls. Here we describe the core technology of the new sequencing platforms, the methodology of data analysis, and different aspects of applications.

Mechanisms and evolution of virulence in oomycetes

Many destructive diseases of plants and animals are caused by oomycetes, a group of eukaryotic pathogens important to agricultural, ornamental, and natural ecosystems. Understanding the mechanisms underlying oomycete virulence and the genomic processes by which those mechanisms rapidly evolve is essential to developing effective long-term control measures for oomycete diseases. Several common mechanisms underlying oomycete virulence, including protein toxins and cell-entering effectors, have emerged from comparing oomycetes with different genome characteristics, parasitic lifestyles, and host ranges. Oomycete genomes display a strongly bipartite organization in which conserved housekeeping genes are concentrated in syntenic gene-rich blocks, whereas virulence genes are dispersed into highly dynamic, repeat-rich regions. There is also evidence that key virulence genes have been acquired by horizontal transfer from other eukaryotic and prokaryotic species.

De novo sequence assembly of Albugo candida reveals a small genome relative to other biotrophic oomycetes

Background

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.

Results

A recent genome analysis of the oomycete Hyaloperonospora arabidopsidis suggests that a reduction in the number of genes encoding secreted pathogenicity proteins, enzymes for assimilation of inorganic nitrogen and sulphur represent a genomic signature for the evolution of obligate biotrophy. Here, we report a draft reference genome of a major crop pathogen Albugo candida (another obligate biotrophic oomycete) with an estimated genome of 45.3 Mb. This is very similar to the genome size of a necrotrophic oomycete Pythium ultimum (43 Mb) but less than half that of H. arabidopsidis (99 Mb). Sequencing of A. candida transcripts from infected host tissue and zoosporangia combined with genome-wide annotation revealed 15,824 predicted genes. Most of the predicted genes lack significant similarity with sequences from other oomycetes. Most intriguingly, A. candida appears to have a much smaller repertoire of pathogenicity-related proteins than H. arabidopsidis including genes that encode RXLR effector proteins, CRINKLER-like genes, and elicitins. Necrosis and Ethylene inducing Peptides were not detected in the genome of A. candida. Putative orthologs of tat-C, a component of the twin arginine translocase system, were identified from multiple oomycete genera along with proteins containing putative tat-secretion signal peptides.

Conclusion

Albugo candida has a comparatively small genome amongst oomycetes, retains motility of sporangial inoculum, and harbours a much smaller repertoire of candidate effectors than was recently reported for H. arabidopsidis. This minimal gene repertoire could indicate a lack of expansion, rather than a reduction, in the number of genes that signify the evolution of biotrophy in oomycetes.

Genome evolution following host jumps in the Irish potato famine pathogen lineage

Many plant pathogens, including those in the lineage of the Irish potato famine organism Phytophthora infestans, evolve by host jumps followed by specialization. However, how host jumps affect genome evolution remains largely unknown. To determine the patterns of sequence variation in the P. infestans lineage, we resequenced six genomes of four sister species. This revealed uneven evolutionary rates across genomes with genes in repeat-rich regions showing higher rates of structural polymorphisms and positive selection. These loci are enriched in genes induced in planta, implicating host adaptation in genome evolution. Unexpectedly, genes involved in epigenetic processes formed another class of rapidly evolving residents of the gene-sparse regions. These results demonstrate that dynamic repeat-rich genome compartments underpin accelerated gene evolution following host jumps in this pathogen lineage.

Signatures of adaptation to obligate biotrophy in the Hyaloperonospora arabidopsidis genome

Many oomycete and fungal plant pathogens are obligate biotrophs, which extract nutrients only from living plant tissue and cannot grow apart from their hosts. Although these pathogens cause substantial crop losses, little is known about the molecular basis or evolution of obligate biotrophy. Here, we report the genome sequence of the oomycete Hyaloperonospora arabidopsidis (Hpa), an obligate biotroph and natural pathogen of Arabidopsis thaliana. In comparison with genomes of related, hemibiotrophic Phytophthora species, the Hpa genome exhibits dramatic reductions in genes encoding (i) RXLR effectors and other secreted pathogenicity proteins, (ii) enzymes for assimilation of inorganic nitrogen and sulfur, and (iii) proteins associated with zoospore formation and motility. These attributes comprise a genomic signature of evolution toward obligate biotrophy.

Recognition of Phytophthora infestans Avr4 by potato R4 is triggered by C-terminal domains comprising W motifs

SUMMARY Oomycete RXLR-dEER effector proteins are rapidly evolving proteins with the selective pressure targeted predominantly at their C-terminal ends. The majority of RXLR-dEER proteins have recognizable motifs of 21-30 amino acids in the C-terminal domain that are named after conserved amino acid residues at fixed positions within the respective motifs. In this article, it is reported that the Phytophthora infestans RXLR-dEER protein Avr4 contains three W motifs and one Y motif in its C-terminal domain. Agroinfection assays using constructs encoding modified forms of PiAvr4 have shown that the region containing the W2 motif, in combination with either the W1 or W3 motif, triggers a necrotic response in potato plants carrying the resistance gene R4. By mining the superfamily of avirulence homologues (Avh) deduced from three sequenced Phytophthora genomes, several Avh proteins were identified as homologues of PiAvr4: six in P. infestans, one in P. ramorum and seven in P. sojae. One very close homologue of PiAvr4 was cloned from the sibling species, P. mirabilis. This species is not pathogenic on potato but, similar to PiAvr4, PmirAvh4 triggered a necrotic response on potato clones carrying R4, but not on clones lacking R4. Genes encoding RXLR-dEER effectors are often located in regions showing genome rearrangements. Alignment of the genomic region harbouring PiAvr4 with syntenic regions in P. sojae and P. ramorum revealed that PiAvr4 is located on a 100-kb indel block and is surrounded by transposable elements.

RXLR-mediated entry of Phytophthora sojae effector Avr1b into soybean cells does not require pathogen-encoded machinery

Effector proteins secreted by oomycete and fungal pathogens have been inferred to enter host cells, where they interact with host resistance gene products. Using the effector protein Avr1b of Phytophthora sojae, an oomycete pathogen of soybean (Glycine max), we show that a pair of sequence motifs, RXLR and dEER, plus surrounding sequences, are both necessary and sufficient to deliver the protein into plant cells. Particle bombardment experiments demonstrate that these motifs function in the absence of the pathogen, indicating that no additional pathogen-encoded machinery is required for effector protein entry into host cells. Furthermore, fusion of the Avr1b RXLR-dEER domain to green fluorescent protein (GFP) allows GFP to enter soybean root cells autonomously. The conclusion that RXLR and dEER serve to transduce oomycete effectors into host cells indicates that the >370 RXLR-dEER-containing proteins encoded in the genome sequence of P. sojae are candidate effectors. We further show that the RXLR and dEER motifs can be replaced by the closely related erythrocyte targeting signals found in effector proteins of Plasmodium, the protozoan that causes malaria in humans. Mutational analysis of the RXLR motif shows that the required residues are very similar in the motifs of Plasmodium and Phytophthora. Thus, the machinery of the hosts (soybean and human) targeted by the effectors may be very ancient.

Conserved C-terminal motifs required for avirulence and suppression of cell death by Phytophthora sojae effector Avr1b

The sequenced genomes of oomycete plant pathogens contain large superfamilies of effector proteins containing the protein translocation motif RXLR-dEER. However, the contributions of these effectors to pathogenicity remain poorly understood. Here, we show that the Phytophthora sojae effector protein Avr1b can contribute positively to virulence and can suppress programmed cell death (PCD) triggered by the mouse BAX protein in yeast, soybean (Glycine max), and Nicotiana benthamiana cells. We identify three conserved motifs (K, W, and Y) in the C terminus of the Avr1b protein and show that mutations in the conserved residues of the W and Y motifs reduce or abolish the ability of Avr1b to suppress PCD and also abolish the avirulence interaction of Avr1b with the Rps1b resistance gene in soybean. W and Y motifs are present in at least half of the identified oomycete RXLR-dEER effector candidates, and we show that three of these candidates also suppress PCD in soybean. Together, these results indicate that the W and Y motifs are critical for the interaction of Avr1b with host plant target proteins and support the hypothesis that these motifs are critical for the functions of the very large number of predicted oomycete effectors that contain them.

An integrated BAC and genome sequence physical map of Phytophthora sojae

Phytophthora spp. are serious pathogens that threaten numerous cultivated crops, trees, and natural vegetation worldwide. The soybean pathogen P. sojae has been developed as a model oomycete. Here, we report a bacterial artificial chromosome (BAC)-based, integrated physical map of the P. sojae genome. We constructed two BAC libraries, digested 8,681 BACs with seven restriction enzymes, end labeled the digested fragments with four dyes, and analyzed them with capillary electrophoresis. Fifteen data sets were constructed from the fingerprints, using individual dyes and all possible combinations, and were evaluated for contig assembly. In all, 257 contigs were assembled from the XhoI data set, collectively spanning approximately 132 Mb in physical length. The BAC contigs were integrated with the draft genome sequence of P. sojae by end sequencing a total of 1,440 BACs that formed a minimal tiling path. This enabled the 257 contigs of the BAC map to be merged with 207 sequence scaffolds to form an integrated map consisting of 79 superscaffolds. The map represents the first genome-wide physical map of a Phytophthora sp. and provides a valuable resource for genomics and molecular biology research in P. sojae and other Phytophthora spp. In one illustration of this value, we have placed the 350 members of a superfamily of putative pathogenicity effector genes onto the map, revealing extensive clustering of these genes.

Comparative analysis of Phytophthora genes encoding secreted proteins reveals conserved synteny and lineage-specific gene duplications and deletions

Comparative analysis of two Phytophthora genomes revealed overall colinearity in four genomic regions consisting of a 1.5-Mb sequence of Phytophthora sojae and a 0.9-Mb sequence of P. ramorum. In these regions with conserved synteny, the gene order is largely similar; however, genome rearrangements also have occurred. Deletions and duplications often were found in association with genes encoding secreted proteins, including effectors that are important for interaction with host plants. Among secreted protein genes, different evolutionary patterns were found. Elicitin genes that code for a complex family of highly conserved Phytophthora-specific elicitors show conservation in gene number and order, and often are clustered. In contrast, the race-specific elicitor gene Avrlb-1 appeared to be missing from the region with conserved synteny, as were its five homologs that are scattered over the four genomic regions. Some gene families encoding secreted proteins were found to be expanded in one species compared with the other. This could be the result of either repeated gene duplications in one species or specific deletions in the other. These different evolutionary patterns may shed light on the functions of these secreted proteins in the biology and pathology of the two Phytophthora spp.

Phytophthora genome sequences uncover evolutionary origins and mechanisms of pathogenesis

Draft genome sequences have been determined for the soybean pathogen Phytophthora sojae and the sudden oak death pathogen Phytophthora ramorum. Oömycetes such as these Phytophthora species share the kingdom Stramenopila with photosynthetic algae such as diatoms, and the presence of many Phytophthora genes of probable phototroph origin supports a photosynthetic ancestry for the stramenopiles. Comparison of the two species' genomes reveals a rapid expansion and diversification of many protein families associated with plant infection such as hydrolases, ABC transporters, protein toxins, proteinase inhibitors, and, in particular, a superfamily of 700 proteins with similarity to known oömycete avirulence genes.

Nonneutral GC3 and retroelement codon mimicry in Phytophthora

Phytophthora is a genus entirely comprised of destructive plant pathogens. It belongs to the Stramenopila, a unique branch of eukaryotes, phylogenetically distinct from plants, animals, or fungi. Phytophthora genes show a strong preference for usage of codons ending with G or C (high GC3). The presence of high GC3 in genes can be utilized to differentiate coding regions from noncoding regions in the genome. We found that both selective pressure and mutation bias drive codon bias in Phytophthora. Indicative for selection pressure is the higher GC3 value of highly expressed genes in different Phytophthora species. Lineage specific GC increase of noncoding regions is reminiscent of whole-genome mutation bias, whereas the elevated Phytophthora GC3 is primarily a result of translation efficiency-driven selection. Heterogeneous retrotransposons exist in Phytophthora genomes and many of them vary in their GC content. Interestingly, the most widespread groups of retroelements in Phytophthora show high GC3 and a codon bias that is similar to host genes. Apparently, selection pressure has been exerted on the retroelement's codon usage, and such mimicry of host codon bias might be beneficial for the propagation of retrotransposons.