Genome-Scale Identification of Essential Metabolic Processes for Targeting the Plasmodium Liver Stage.

Plasmodium gene functions in mosquito and liver stages remain poorly characterized due to limitations in the throughput of phenotyping at these stages. To fill this gap, we followed more than 1,300 barcoded P. berghei mutants through the life cycle. We discover 461 genes required for efficient parasite transmission to mosquitoes through the liver stage and back into the bloodstream of mice. We analyze the screen in the context of genomic, transcriptomic, and metabolomic data by building a thermodynamic model of P. berghei liver-stage metabolism, which shows a major reprogramming of parasite metabolism to achieve rapid growth in the liver. We identify seven metabolic subsystems that become essential at the liver stages compared with asexual blood stages: type II fatty acid synthesis and elongation (FAE), tricarboxylic acid, amino sugar, heme, lipoate, and shikimate metabolism. Selected predictions from the model are individually validated in single mutants to provide future targets for drug development.

Functional Profiling of a Plasmodium Genome Reveals an Abundance of Essential Genes.

The genomes of malaria parasites contain many genes of unknown function. To assist drug development through the identification of essential genes and pathways, we have measured competitive growth rates in mice of 2,578 barcoded Plasmodium berghei knockout mutants, representing >50% of the genome, and created a phenotype database. At a single stage of its complex life cycle, P. berghei requires two-thirds of genes for optimal growth, the highest proportion reported from any organism and a probable consequence of functional optimization necessitated by genomic reductions during the evolution of parasitism. In contrast, extreme functional redundancy has evolved among expanded gene families operating at the parasite-host interface. The level of genetic redundancy in a single-celled organism may thus reflect the degree of environmental variation it experiences. In the case of Plasmodium parasites, this helps rationalize both the relative successes of drugs and the greater difficulty of making an effective vaccine.

A Genome-wide CRISPR Screen in Toxoplasma Identifies Essential Apicomplexan Genes.

Apicomplexan parasites are leading causes of human and livestock diseases such as malaria and toxoplasmosis, yet most of their genes remain uncharacterized. Here, we present the first genome-wide genetic screen of an apicomplexan. We adapted CRISPR/Cas9 to assess the contribution of each gene from the parasite Toxoplasma gondii during infection of human fibroblasts. Our analysis defines ∼200 previously uncharacterized, fitness-conferring genes unique to the phylum, from which 16 were investigated, revealing essential functions during infection of human cells. Secondary screens identify as an invasion factor the claudin-like apicomplexan microneme protein (CLAMP), which resembles mammalian tight-junction proteins and localizes to secretory organelles, making it critical to the initiation of infection. CLAMP is present throughout sequenced apicomplexan genomes and is essential during the asexual stages of the malaria parasite Plasmodium falciparum. These results provide broad-based functional information on T. gondii genes and will facilitate future approaches to expand the horizon of antiparasitic interventions.

Ancient Plasmodium genomes shed light on the history of human malaria.

Malaria-causing protozoa of the genus Plasmodium have exerted one of the strongest selective pressures on the human genome, and resistance alleles provide biomolecular footprints that outline the historical reach of these species1. Nevertheless, debate persists over when and how malaria parasites emerged as human pathogens and spread around the globe1,2. To address these questions, we generated high-coverage ancient mitochondrial and nuclear genome-wide data from P. falciparum, P. vivax and P. malariae from 16 countries spanning around 5,500 years of human history. We identified P. vivax and P. falciparum across geographically disparate regions of Eurasia from as early as the fourth and first millennia BCE, respectively; for P. vivax, this evidence pre-dates textual references by several millennia3. Genomic analysis supports distinct disease histories for P. falciparum and P. vivax in the Americas: similarities between now-eliminated European and peri-contact South American strains indicate that European colonizers were the source of American P. vivax, whereas the trans-Atlantic slave trade probably introduced P. falciparum into the Americas. Our data underscore the role of cross-cultural contacts in the dissemination of malaria, laying the biomolecular foundation for future palaeo-epidemiological research into the impact of Plasmodium parasites on human history. Finally, our unexpected discovery of P. falciparum in the high-altitude Himalayas provides a rare case study in which individual mobility can be inferred from infection status, adding to our knowledge of cross-cultural connectivity in the region nearly three millennia ago.

The architecture of a scrambled genome reveals massive levels of genomic rearrangement during development.

Programmed DNA rearrangements in the single-celled eukaryote Oxytricha trifallax completely rewire its germline into a somatic nucleus during development. This elaborate, RNA-mediated pathway eliminates noncoding DNA sequences that interrupt gene loci and reorganizes the remaining fragments by inversions and permutations to produce functional genes. Here, we report the Oxytricha germline genome and compare it to the somatic genome to present a global view of its massive scale of genome rearrangements. The remarkably encrypted genome architecture contains >3,500 scrambled genes, as well as >800 predicted germline-limited genes expressed, and some posttranslationally modified, during genome rearrangements. Gene segments for different somatic loci often interweave with each other. Single gene segments can contribute to multiple, distinct somatic loci. Terminal precursor segments from neighboring somatic loci map extremely close to each other, often overlapping. This genome assembly provides a draft of a scrambled genome and a powerful model for studies of genome rearrangement.

Genomes on the edge: programmed genome instability in ciliates.

Ciliates are an ancient and diverse group of microbial eukaryotes that have emerged as powerful models for RNA-mediated epigenetic inheritance. They possess extensive sets of both tiny and long noncoding RNAs that, together with a suite of proteins that includes transposases, orchestrate a broad cascade of genome rearrangements during somatic nuclear development. This Review emphasizes three important themes: the remarkable role of RNA in shaping genome structure, recent discoveries that unify many deeply diverged ciliate genetic systems, and a surprising evolutionary "sign change" in the role of small RNAs between major species groups.

Ancestral aneuploidy and stable chromosomal duplication resulting in differential genome structure and gene expression control in trypanosomatid parasites.

Aneuploidy is widely observed in both unicellular and multicellular eukaryotes, usually associated with adaptation to stress conditions. Chromosomal duplication stability is a tradeoff between the fitness cost of having unbalanced gene copies and the potential fitness gained from increased dosage of specific advantageous genes. Trypanosomatids, a family of protozoans that include species that cause neglected tropical diseases, are a relevant group to study aneuploidies. Their life cycle has several stressors that could select for different patterns of chromosomal duplications and/or losses, and their nearly universal use of polycistronic transcription increases their reliance on gene expansion/contraction, as well as post-transcriptional control as mechanisms for gene expression regulation. By evaluating the data from 866 isolates covering seven trypanosomatid genera, we have revealed that aneuploidy tolerance is an ancestral characteristic of trypanosomatids but has a reduced occurrence in a specific monophyletic clade that has undergone large genomic reorganization and chromosomal fusions. We have also identified an ancient chromosomal duplication that was maintained across these parasite's speciation, named collectively as the trypanosomatid ancestral supernumerary chromosome (TASC). TASC has most genes in the same coding strand, is expressed as a disomic chromosome (even having four copies), and has increased potential for functional variation, but it purges highly deleterious mutations more efficiently than other chromosomes. The evidence of stringent control over gene expression in this chromosome suggests that these parasites have adapted to mitigate the fitness cost associated with this ancient chromosomal duplication.

Comparative genomics of Cryptosporidium parvum reveals the emergence of an outbreak-associated population in Europe and its spread to the United States.

The zoonotic parasite Cryptosporidium parvum is a global cause of gastrointestinal disease in humans and ruminants. Sequence analysis of the highly polymorphic gp60 gene enabled the classification of C. parvum isolates into multiple groups (e.g., IIa, IIc, Id) and a large number of subtypes. In Europe, subtype IIaA15G2R1 is largely predominant and has been associated with many water- and food-borne outbreaks. In this study, we generated new whole-genome sequence (WGS) data from 123 human- and ruminant-derived isolates collected in 13 European countries and included other available WGS data from Europe, Egypt, China, and the United States (n = 72) in the largest comparative genomics study to date. We applied rigorous filters to exclude mixed infections and analyzed a data set from 141 isolates from the zoonotic groups IIa (n = 119) and IId (n = 22). Based on 28,047 high-quality, biallelic genomic SNPs, we identified three distinct and strongly supported populations: Isolates from China (IId) and Egypt (IIa and IId) formed population 1; a minority of European isolates (IIa and IId) formed population 2; and the majority of European (IIa, including all IIaA15G2R1 isolates) and all isolates from the United States (IIa) clustered in population 3. Based on analyses of the population structure, population genetics, and recombination, we show that population 3 has recently emerged and expanded throughout Europe to then, possibly from the United Kingdom, reach the United States, where it also expanded. The reason(s) for the successful spread of population 3 remain elusive, although genes under selective pressure uniquely in this population were identified.

Piwi-interacting RNAs protect DNA against loss during Oxytricha genome rearrangement.

Genome duality in ciliated protozoa offers a unique system to showcase their epigenome as a model of inheritance. In Oxytricha, the somatic genome is responsible for vegetative growth, whereas the germline contributes DNA to the next sexual generation. Somatic nuclear development removes all transposons and other so-called "junk" DNA, which comprise ~95% of the germline. We demonstrate that Piwi-interacting small RNAs (piRNAs) from the maternal nucleus can specify genomic regions for retention in this process. Oxytricha piRNAs map primarily to the somatic genome, representing the ~5% of the germline that is retained. Furthermore, injection of synthetic piRNAs corresponding to normally deleted regions leads to their retention in later generations. Our findings highlight small RNAs as powerful transgenerational carriers of epigenetic information for genome programming.

PbARID-associated chromatin remodeling events are essential for gametocyte development in Plasmodium.

Gametocyte development of the Plasmodium parasite is a key step for transmission of the parasite. Male and female gametocytes are produced from a subpopulation of asexual blood-stage parasites, but the mechanisms that regulate the differentiation of sexual stages are still under investigation. In this study, we investigated the role of PbARID, a putative subunit of a SWI/SNF chromatin remodeling complex, in transcriptional regulation during the gametocyte development of P. berghei. PbARID expression starts in early gametocytes before the manifestation of male and female-specific features, and disruption of its gene results in the complete loss of gametocytes with detectable male features and the production of abnormal female gametocytes. ChIP-seq analysis of PbARID showed that it forms a complex with gSNF2, an ATPase subunit of the SWI/SNF chromatin remodeling complex, associating with the male cis-regulatory element, TGTCT. Further ChIP-seq of PbARID in gsnf2-knockout parasites revealed an association of PbARID with another cis-regulatory element, TGCACA. RIME and DNA-binding assays suggested that HDP1 is the transcription factor that recruits PbARID to the TGCACA motif. Our results indicated that PbARID could function in two chromatin remodeling events and paly essential roles in both male and female gametocyte development.

Selective whole-genome amplification reveals population genetics of Leishmania braziliensis directly from patient skin biopsies.

In Brazil, Leishmania braziliensis is the main causative agent of the neglected tropical disease, cutaneous leishmaniasis (CL). CL presents on a spectrum of disease severity with a high rate of treatment failure. Yet the parasite factors that contribute to disease presentation and treatment outcome are not well understood, in part because successfully isolating and culturing parasites from patient lesions remains a major technical challenge. Here we describe the development of selective whole genome amplification (SWGA) for Leishmania and show that this method enables culture-independent analysis of parasite genomes obtained directly from primary patient skin samples, allowing us to circumvent artifacts associated with adaptation to culture. We show that SWGA can be applied to multiple Leishmania species residing in different host species, suggesting that this method is broadly useful in both experimental infection models and clinical studies. SWGA carried out directly on skin biopsies collected from patients in Corte de Pedra, Bahia, Brazil, showed extensive genomic diversity. Finally, as a proof-of-concept, we demonstrated that SWGA data can be integrated with published whole genome data from cultured parasite isolates to identify variants unique to specific geographic regions in Brazil where treatment failure rates are known to be high. SWGA provides a relatively simple method to generate Leishmania genomes directly from patient samples, unlocking the potential to link parasite genetics with host clinical phenotypes.

Conspicuous chloroplast with light harvesting-photosystem I/II megacomplex in marine Prorocentrum cordatum.

Marine photosynthetic (micro)organisms drive multiple biogeochemical cycles and display a large diversity. Among them, the bloom-forming, free-living dinoflagellate Prorocentrum cordatum CCMP 1329 (formerly P. minimum) stands out with its distinct cell biological features. Here, we obtained insights into the structural properties of the chloroplast and the photosynthetic machinery of P. cordatum using microscopic and proteogenomic approaches. High-resolution FIB/SEM analysis revealed a single large chloroplast (∼40% of total cell volume) with a continuous barrel-like structure, completely lining the inner face of the cell envelope and enclosing a single reticular mitochondrium, the Golgi apparatus, as well as diverse storage inclusions. Enriched thylakoid membrane fractions of P. cordatum were comparatively analyzed with those of the well-studied model-species Arabidopsis (Arabidopsis thaliana) using 2D BN DIGE. Strikingly, P. cordatum possessed a large photosystem-light harvesting megacomplex (>1.5 MDa), which is dominated by photosystems I and II (PSI, PSII), chloroplast complex I, and chlorophyll a-b binding light harvesting complex proteins. This finding parallels the absence of grana in its chloroplast and distinguishes from the predominant separation of PSI and PSII complexes in A. thaliana, indicating a different mode of flux balancing. Except for the core elements of the ATP synthase and the cytb6f-complex, the composition of the other complexes (PSI, PSII, and pigment-binding proteins, PBPs) of P. cordatum differed markedly from those of A. thaliana. Furthermore, a high number of PBPs was detected, accounting for a large share of the total proteomic data (∼65%) and potentially providing P. cordatum with flexible adaptation to changing light regimes.

A widespread coral-infecting apicomplexan with chlorophyll biosynthesis genes.

Apicomplexa is a group of obligate intracellular parasites that includes the causative agents of human diseases such as malaria and toxoplasmosis. Apicomplexans evolved from free-living phototrophic ancestors, but how this transition to parasitism occurred remains unknown. One potential clue lies in coral reefs, of which environmental DNA surveys have uncovered several lineages of uncharacterized basally branching apicomplexans1,2. Reef-building corals have a well-studied symbiotic relationship with photosynthetic Symbiodiniaceae dinoflagellates (for example, Symbiodinium3), but the identification of other key microbial symbionts of corals has proven to be challenging4,5. Here we use community surveys, genomics and microscopy analyses to identify an apicomplexan lineage-which we informally name 'corallicolids'-that was found at a high prevalence (over 80% of samples, 70% of genera) across all major groups of corals. Corallicolids were the second most abundant coral-associated microeukaryotes after the Symbiodiniaceae, and are therefore core members of the coral microbiome. In situ fluorescence and electron microscopy confirmed that corallicolids live intracellularly within the tissues of the coral gastric cavity, and that they possess apicomplexan ultrastructural features. We sequenced the genome of the corallicolid plastid, which lacked all genes for photosystem proteins; this indicates that corallicolids probably contain a non-photosynthetic plastid (an apicoplast6). However, the corallicolid plastid differs from all other known apicoplasts because it retains the four ancestral genes that are involved in chlorophyll biosynthesis. Corallicolids thus share characteristics with both their parasitic and their free-living relatives, which suggests that they are evolutionary intermediates and implies the existence of a unique biochemistry during the transition from phototrophy to parasitism.

Genomic insights into the cellular specialization of predation in raptorial protists.

BACKGROUND: Predation is a fundamental mechanism for organisms to acquire energy, and various species have evolved diverse tools to enhance their hunting abilities. Among protozoan predators, raptorial Haptorian ciliates are particularly fascinating as they possess offensive extrusomes known as toxicysts, which are rapidly discharged upon prey contact. However, our understanding of the genetic processes and specific toxins involved in toxicyst formation and discharge is still limited. RESULTS: In this study, we investigated the predation strategies and subcellular structures of seven Haptoria ciliate species and obtained their genome sequences using single-cell sequencing technology. Comparative genomic analysis revealed distinct gene duplications related to membrane transport proteins and hydrolytic enzymes in Haptoria, which play a crucial role in the production and discharge of toxicysts. Transcriptomic analysis further confirmed the abundant expression of genes related to membrane transporters and cellular toxins in Haptoria compared to Trichostomatia. Notably, polyketide synthases (PKS) and L-amino acid oxidases (LAAO) were identified as potentially toxin genes that underwent extensive duplication events in Haptoria. CONCLUSIONS: Our results shed light on the evolutionary and genomic adaptations of Haptorian ciliates for their predation strategies in evolution and provide insights into their toxic mechanisms.

Comparative genome-wide identification and characterization of SET domain-containing and JmjC domain-containing proteins in piroplasms.

BACKGROUND: SET domain-containing histone lysine methyltransferases (HKMTs) and JmjC domain-containing histone demethylases (JHDMs) are essential for maintaining dynamic changes in histone methylation across parasite development and infection. However, information on the HKMTs and JHDMs in human pathogenic piroplasms, such as Babesia duncani and Babesia microti, and in veterinary important pathogens, including Babesia bigemina, Babesia bovis, Theileria annulata and Theileria parva, is limited. RESULTS: A total of 38 putative KMTs and eight JHDMs were identified using a comparative genomics approach. Phylogenetic analysis revealed that the putative KMTs can be divided into eight subgroups, while the JHDMs belong to the JARID subfamily, except for BdJmjC1 (BdWA1_000016) and TpJmjC1 (Tp Muguga_02g00471) which cluster with JmjC domain only subfamily members. The motifs of SET and JmjC domains are highly conserved among piroplasm species. Interspecies collinearity analysis provided insight into the evolutionary duplication events of some SET domain and JmjC domain gene families. Moreover, relative gene expression analysis by RT‒qPCR demonstrated that the putative KMT and JHDM gene families were differentially expressed in different intraerythrocytic developmental stages of B. duncani, suggesting their role in Apicomplexa parasite development. CONCLUSIONS: Our study provides a theoretical foundation and guidance for understanding the basic characteristics of several important piroplasm KMT and JHDM families and their biological roles in parasite differentiation.

Autochthonous Leishmaniasis Caused by Leishmania tropica, Identified by Using Whole-Genome Sequencing, Sri Lanka.

Cutaneous leishmaniasis is atypical in Sri Lanka because Leishmania donovani, which typically causes visceral disease, is the causative agent. The origins of recently described hybrids between L. donovani and other Leishmania spp. usually responsible for cutaneous leishmaniasis remain unknown. Other endemic dermotropic Leishmania spp. have not been reported in Sri Lanka. Genome analysis of 27 clinical isolates from Sri Lanka and 32 Old World Leishmania spp. strains found 8 patient isolates clustered with L. tropica and 19 with L. donovani. The L. tropica isolates from Sri Lanka shared markers with strain LtK26 reported decades ago in India, indicating they were not products of recent interspecies hybridization. Because L. tropica was isolated from patients with leishmaniasis in Sri Lanka, our findings indicate L. donovani is not the only cause of cutaneous leishmaniasis in Sri Lanka and potentially explains a haplotype that led to interspecies dermotropic L. donovani hybrids.

Read, Write, Adapt: Challenges and Opportunities during Kinetoplastid Genome Replication.

The genomes of all organisms are read throughout their growth and development, generating new copies during cell division and encoding the cellular activities dictated by the genome's content. However, genomes are not invariant information stores but are purposefully altered in minor and major ways, adapting cellular behaviour and driving evolution. Kinetoplastids are eukaryotic microbes that display a wide range of such read-write genome activities, in many cases affecting critical aspects of their biology, such as host adaptation. Here we discuss the range of read-write genome changes found in two well-studied kinetoplastid parasites, Trypanosoma brucei and Leishmania, focusing on recent work that suggests such adaptive genome variation is linked to novel strategies the parasites use to replicate their unconventional genomes.

Chromosome-level genome of Schistosoma haematobium underpins genome-wide explorations of molecular variation.

Urogenital schistosomiasis is caused by the blood fluke Schistosoma haematobium and is one of the most neglected tropical diseases worldwide, afflicting > 100 million people. It is characterised by granulomata, fibrosis and calcification in urogenital tissues, and can lead to increased susceptibility to HIV/AIDS and squamous cell carcinoma of the bladder. To complement available treatment programs and break the transmission of disease, sound knowledge and understanding of the biology and ecology of S. haematobium is required. Hybridisation/introgression events and molecular variation among members of the S. haematobium-group might effect important biological and/or disease traits as well as the morbidity of disease and the effectiveness of control programs including mass drug administration. Here we report the first chromosome-contiguous genome for a well-defined laboratory line of this blood fluke. An exploration of this genome using transcriptomic data for all key developmental stages allowed us to refine gene models (including non-coding elements) and annotations, discover 'new' genes and transcription profiles for these stages, likely linked to development and/or pathogenesis. Molecular variation within S. haematobium among some geographical locations in Africa revealed unique genomic 'signatures' that matched species other than S. haematobium, indicating the occurrence of introgression events. The present reference genome (designated Shae.V3) and the findings from this study solidly underpin future functional genomic and molecular investigations of S. haematobium and accelerate systematic, large-scale population genomics investigations, with a focus on improved and sustained control of urogenital schistosomiasis.

Massive colonization of protein-coding exons by selfish genetic elements in Paramecium germline genomes.

Ciliates are unicellular eukaryotes with both a germline genome and a somatic genome in the same cytoplasm. The somatic macronucleus (MAC), responsible for gene expression, is not sexually transmitted but develops from a copy of the germline micronucleus (MIC) at each sexual generation. In the MIC genome of Paramecium tetraurelia, genes are interrupted by tens of thousands of unique intervening sequences called internal eliminated sequences (IESs), which have to be precisely excised during the development of the new MAC to restore functional genes. To understand the evolutionary origin of this peculiar genomic architecture, we sequenced the MIC genomes of 9 Paramecium species (from approximately 100 Mb in Paramecium aurelia species to >1.5 Gb in Paramecium caudatum). We detected several waves of IES gains, both in ancestral and in more recent lineages. While the vast majority of IESs are single copy in present-day genomes, we identified several families of mobile IESs, including nonautonomous elements acquired via horizontal transfer, which generated tens to thousands of new copies. These observations provide the first direct evidence that transposable elements can account for the massive proliferation of IESs in Paramecium. The comparison of IESs of different evolutionary ages indicates that, over time, IESs shorten and diverge rapidly in sequence while they acquire features that allow them to be more efficiently excised. We nevertheless identified rare cases of IESs that are under strong purifying selection across the aurelia clade. The cases examined contain or overlap cellular genes that are inactivated by excision during development, suggesting conserved regulatory mechanisms. Similar to the evolution of introns in eukaryotes, the evolution of Paramecium IESs highlights the major role played by selfish genetic elements in shaping the complexity of genome architecture and gene expression.

Expression and biochemical characterization of the putative insulinase enzyme PF11_0189 found in the Plasmodium falciparum genome.

PF11_0189 is a putative insulin degrading enzyme present in Plasmodium falciparum genome. The catalytic domain of PF11_0189 is about 27 kDa. Substrate specificity study shows PF11_0189 acts upon different types of proteins. The substrate specificity is found to be highest when insulin is used as a substrate. Metal dependency study shows highest dependency of PF11_0189 towards zinc metal for its proteolytic activity. Chelation of zinc metal with EDTA shows complete absence of PF11_0189 activity. Peptide inhibitors, P-70 and P-121 from combinatorial peptide library prepared against PF11_0189 show inhibition with an IC50 value of 4.8 µM and 7.5 µM respectively. A proven natural anti-malarial peptide cyclosporin A shows complete inhibition against PF11_0189 with an IC50 value of 0.75 µM suggesting PF11_0189 as a potential target for peptide inhibitors. The study implicates that PF11_0189 is a zinc metalloprotease involved in catalysis of insulin. The study gives a preliminary insight into the mechanism of complications arising from glucose abnormalities during severe malaria.