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.

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.

The piggyBac transposon-derived genes TPB1 and TPB6 mediate essential transposon-like excision during the developmental rearrangement of key genes in Tetrahymena thermophila.

Ciliated protozoans perform extreme forms of programmed somatic DNA rearrangement during development. The model ciliate Tetrahymena thermophila removes 34% of its germline micronuclear genome from somatic macronuclei by excising thousands of internal eliminated sequences (IESs), a process that shares features with transposon excision. Indeed, piggyBac transposon-derived genes are necessary for genome-wide IES excision in both Tetrahymena (TPB2 [Tetrahymena piggyBac-like 2] and LIA5) and Paramecium tetraurelia (PiggyMac). T. thermophila has at least three other piggyBac-derived genes: TPB1, TPB6, and TPB7 Here, we show that TPB1 and TPB6 excise a small, distinct set of 12 unusual IESs that disrupt exons. TPB1-deficient cells complete mating, but their progeny exhibit slow growth, giant vacuoles, and osmotic shock sensitivity due to retention of an IES in the vacuolar gene DOP1 (Dopey domain-containing protein). Unlike most IESs, TPB1-dependent IESs have piggyBac-like terminal inverted motifs that are necessary for excision. Transposon-like excision mediated by TPB1 and TPB6 provides direct evidence for a transposon origin of not only IES excision machinery but also IESs themselves. Our study highlights a division of labor among ciliate piggyBac-derived genes, which carry out mutually exclusive categories of excision events mediated by either transposon-like features or RNA-directed heterochromatin.

Detection of simple and complex de novo mutations with multiple reference sequences.

The characterization of de novo mutations in regions of high sequence and structural diversity from whole-genome sequencing data remains highly challenging. Complex structural variants tend to arise in regions of high repetitiveness and low complexity, challenging both de novo assembly, in which short reads do not capture the long-range context required for resolution, and mapping approaches, in which improper alignment of reads to a reference genome that is highly diverged from that of the sample can lead to false or partial calls. Long-read technologies can potentially solve such problems but are currently unfeasible to use at scale. Here we present Corticall, a graph-based method that combines the advantages of multiple technologies and prior data sources to detect arbitrary classes of genetic variant. We construct multisample, colored de Bruijn graphs from short-read data for all samples, align long-read-derived haplotypes and multiple reference data sources to restore graph connectivity information, and call variants using graph path-finding algorithms and a model for simultaneous alignment and recombination. We validate and evaluate the approach using extensive simulations and use it to characterize the rate and spectrum of de novo mutation events in 119 progeny from four Plasmodium falciparum experimental crosses, using long-read data on the parents to inform reconstructions of the progeny and to detect several known and novel nonallelic homologous recombination events.

Strain-specific genome evolution in Trypanosoma cruzi, the agent of Chagas disease.

The protozoan Trypanosoma cruzi almost invariably establishes life-long infections in humans and other mammals, despite the development of potent host immune responses that constrain parasite numbers. The consistent, decades-long persistence of T. cruzi in human hosts arises at least in part from the remarkable level of genetic diversity in multiple families of genes encoding the primary target antigens of anti-parasite immune responses. However, the highly repetitive nature of the genome-largely a result of these same extensive families of genes-have prevented a full understanding of the extent of gene diversity and its maintenance in T. cruzi. In this study, we have combined long-read sequencing and proximity ligation mapping to generate very high-quality assemblies of two T. cruzi strains representing the apparent ancestral lineages of the species. These assemblies reveal not only the full repertoire of the members of large gene families in the two strains, demonstrating extreme diversity within and between isolates, but also provide evidence of the processes that generate and maintain that diversity, including extensive gene amplification, dispersion of copies throughout the genome and diversification via recombination and in situ mutations. Gene amplification events also yield significant copy number variations in a substantial number of genes presumably not required for or involved in immune evasion, thus forming a second level of strain-dependent variation in this species. The extreme genome flexibility evident in T. cruzi also appears to create unique challenges with respect to preserving core genome functions and gene expression that sets this species apart from related kinetoplastids.

Comparative analyses of parasites with a comprehensive database of genome-scale metabolic models.

Protozoan parasites cause diverse diseases with large global impacts. Research on the pathogenesis and biology of these organisms is limited by economic and experimental constraints. Accordingly, studies of one parasite are frequently extrapolated to infer knowledge about another parasite, across and within genera. Model in vitro or in vivo systems are frequently used to enhance experimental manipulability, but these systems generally use species related to, yet distinct from, the clinically relevant causal pathogen. Characterization of functional differences among parasite species is confined to post hoc or single target studies, limiting the utility of this extrapolation approach. To address this challenge and to accelerate parasitology research broadly, we present a functional comparative analysis of 192 genomes, representing every high-quality, publicly-available protozoan parasite genome including Plasmodium, Toxoplasma, Cryptosporidium, Entamoeba, Trypanosoma, Leishmania, Giardia, and other species. We generated an automated metabolic network reconstruction pipeline optimized for eukaryotic organisms. These metabolic network reconstructions serve as biochemical knowledgebases for each parasite, enabling qualitative and quantitative comparisons of metabolic behavior across parasites. We identified putative differences in gene essentiality and pathway utilization to facilitate the comparison of experimental findings and discovered that phylogeny is not the sole predictor of metabolic similarity. This knowledgebase represents the largest collection of genome-scale metabolic models for both pathogens and eukaryotes; with this resource, we can predict species-specific functions, contextualize experimental results, and optimize selection of experimental systems for fastidious species.

Identity-by-descent with uncertainty characterises connectivity of Plasmodium falciparum populations on the Colombian-Pacific coast.

Characterising connectivity between geographically separated biological populations is a common goal in many fields. Recent approaches to understanding connectivity between malaria parasite populations, with implications for disease control efforts, have used estimates of relatedness based on identity-by-descent (IBD). However, uncertainty around estimated relatedness has not been accounted for. IBD-based relatedness estimates with uncertainty were computed for pairs of monoclonal Plasmodium falciparum samples collected from five cities on the Colombian-Pacific coast where long-term clonal propagation of P. falciparum is frequent. The cities include two official ports, Buenaventura and Tumaco, that are separated geographically but connected by frequent marine traffic. Fractions of highly-related sample pairs (whose classification using a threshold accounts for uncertainty) were greater within cities versus between. However, based on both highly-related fractions and on a threshold-free approach (Wasserstein distances between parasite populations) connectivity between Buenaventura and Tumaco was disproportionally high. Buenaventura-Tumaco connectivity was consistent with transmission events involving parasites from five clonal components (groups of statistically indistinguishable parasites identified under a graph theoretic framework). To conclude, P. falciparum population connectivity on the Colombian-Pacific coast abides by accessibility not isolation-by-distance, potentially implicating marine traffic in malaria transmission with opportunities for targeted intervention. Further investigations are required to test this hypothesis. For the first time in malaria epidemiology (and to our knowledge in ecological and epidemiological studies more generally), we account for uncertainty around estimated relatedness (an important consideration for studies that plan to use genotype versus whole genome sequence data to estimate IBD-based relatedness); we also use threshold-free methods to compare parasite populations and identify clonal components. Threshold-free methods are especially important in analyses of malaria parasites and other recombining organisms with mixed mating systems where thresholds do not have clear interpretation (e.g. due to clonal propagation) and thus undermine the cross-comparison of studies.

The Paramecium histone chaperone Spt16-1 is required for Pgm endonuclease function in programmed genome rearrangements.

In Paramecium tetraurelia, a large proportion of the germline genome is reproducibly removed from the somatic genome after sexual events via a process involving small (s)RNA-directed heterochromatin formation and DNA excision and repair. How germline limited DNA sequences are specifically recognized in the context of chromatin remains elusive. Here, we use a reverse genetics approach to identify factors involved in programmed genome rearrangements. We have identified a P. tetraurelia homolog of the highly conserved histone chaperone Spt16 subunit of the FACT complex, Spt16-1, and show its expression is developmentally regulated. A functional GFP-Spt16-1 fusion protein localized exclusively in the nuclei where genome rearrangements take place. Gene silencing of Spt16-1 showed it is required for the elimination of all germline-limited sequences, for the survival of sexual progeny, and for the accumulation of internal eliminated sequence (ies)RNAs, an sRNA population produced when elimination occurs. Normal accumulation of 25 nt scanRNAs and deposition of silent histone marks H3K9me3 and H3K27me3 indicated that Spt16-1 does not regulate the scanRNA-directed heterochromatin pathway involved in the early steps of DNA elimination. We further show that Spt16-1 is required for the correct nuclear localization of the PiggyMac (Pgm) endonuclease, which generates the DNA double-strand breaks required for DNA elimination. Thus, Spt16-1 is essential for Pgm function during programmed genome rearrangements. We propose a model in which Spt16-1 mediates interactions between the excision machinery and chromatin, facilitating endonuclease access to DNA cleavage sites during genome rearrangements.

Structure of the Plasmodium-interspersed repeat proteins of the malaria parasite.

The deadly symptoms of malaria occur as Plasmodium parasites replicate within blood cells. Members of several variant surface protein families are expressed on infected blood cell surfaces. Of these, the largest and most ubiquitous are the Plasmodium-interspersed repeat (PIR) proteins, with more than 1,000 variants in some genomes. Their functions are mysterious, but differential pir gene expression associates with acute or chronic infection in a mouse malaria model. The membership of the PIR superfamily, and whether the family includes Plasmodium falciparum variant surface proteins, such as RIFINs and STEVORs, is controversial. Here we reveal the structure of the extracellular domain of a PIR from Plasmodium chabaudi We use structure-guided sequence analysis and molecular modeling to show that this fold is found across PIR proteins from mouse- and human-infective malaria parasites. Moreover, we show that RIFINs and STEVORs are not PIRs. This study provides a structure-guided definition of the PIRs and a molecular framework to understand their evolution.

A molecular barcode to inform the geographical origin and transmission dynamics of Plasmodium vivax malaria.

Although Plasmodium vivax parasites are the predominant cause of malaria outside of sub-Saharan Africa, they not always prioritised by elimination programmes. P. vivax is resilient and poses challenges through its ability to re-emerge from dormancy in the human liver. With observed growing drug-resistance and the increasing reports of life-threatening infections, new tools to inform elimination efforts are needed. In order to halt transmission, we need to better understand the dynamics of transmission, the movement of parasites, and the reservoirs of infection in order to design targeted interventions. The use of molecular genetics and epidemiology for tracking and studying malaria parasite populations has been applied successfully in P. falciparum species and here we sought to develop a molecular genetic tool for P. vivax. By assembling the largest set of P. vivax whole genome sequences (n = 433) spanning 17 countries, and applying a machine learning approach, we created a 71 SNP barcode with high predictive ability to identify geographic origin (91.4%). Further, due to the inclusion of markers for within population variability, the barcode may also distinguish local transmission networks. By using P. vivax data from a low-transmission setting in Malaysia, we demonstrate the potential ability to infer outbreak events. By characterising the barcoding SNP genotypes in P. vivax DNA sourced from UK travellers (n = 132) to ten malaria endemic countries predominantly not used in the barcode construction, we correctly predicted the geographic region of infection origin. Overall, the 71 SNP barcode outperforms previously published genotyping methods and when rolled-out within new portable platforms, is likely to be an invaluable tool for informing targeted interventions towards elimination of this resilient human malaria.

Unraveling the antitrypanosomal mechanism of benznidazole and related 2-nitroimidazoles: From prodrug activation to DNA damage.

Nitroheterocycles represent an important class of compound used to treat trypanosomiasis. They often function as prodrugs and can undergo type I nitroreductase (NTR1)-mediated activation before promoting their antiparasitic activities although the nature of these downstream effects has yet to be determined. Here, we show that in an NTR1-dependent process, benznidazole promotes DNA damage in the nuclear genome of Trypanosoma brucei, providing the first direct link between activation of this prodrug and a downstream trypanocidal mechanism. Phenotypic and protein expression studies revealed that components of the trypanosome's homologous recombination (HR) repair pathway (TbMRE11, γH2A, TbRAD51) cooperate to resolve the benznidazole-induced damage, indicating that the prodrug-induced lesions are most likely double stand DNA breaks, while the sequence/recruitment kinetics of these factors parallels that in other eukaryotes HR systems. When extended to other NTR1-activated 2-nitroimidazoles, some were shown to promote DNA damage. Intriguingly, the lesions induced by these required TbMRE11 and TbCSB activities to fix leading us to postulate that TbCSB may operate in systems other than the transcription-coupled nucleotide excision repair pathway. Understanding how existing trypanosomal drugs work will aid future drug design and help unlock novel reactions/pathways that could be exploited as targets for therapeutic intervention.

The non-LTR retrotransposons of Entamoeba histolytica: genomic organization and biology.

Genome sequence analysis of Entamoeba species revealed various classes of transposable elements. While E. histolytica and E. dispar are rich in non-long terminal repeat (LTR) retrotransposons, E. invadens contains predominantly DNA transposons. Non-LTR retrotransposons of E. histolytica constitute three families of long interspersed nuclear elements (LINEs), and their short, nonautonomous partners, SINEs. They occupy ~ 11% of the genome. The EhLINE1/EhSINE1 family is the most abundant and best studied. EhLINE1 is 4.8 kb, with two ORFs that encode functions needed for retrotransposition. ORF1 codes for the nucleic acid-binding protein, and ORF2 has domains for reverse transcriptase (RT) and endonuclease (EN). Most copies of EhLINEs lack complete ORFs. ORF1p is expressed constitutively, but ORF2p is not detected. Retrotransposition could be demonstrated upon ectopic over expression of ORF2p, showing that retrotransposition machinery is functional. The newly retrotransposed sequences showed a high degree of recombination. In transcriptomic analysis, RNA-Seq reads were mapped to individual EhLINE1 copies. Although full-length copies were transcribed, no full-length 4.8 kb transcripts were seen. Rather, sense transcripts mapped to ORF1, RT and EN domains. Intriguingly, there was strong antisense transcription almost exclusively from the RT domain. These unique features of EhLINE1 could serve to attenuate retrotransposition in E. histolytica.

Diversification of small RNA amplification mechanisms for targeting transposon-related sequences in ciliates.

The silencing of repetitive transposable elements (TEs) is ensured by signal amplification of the initial small RNA trigger, which occurs at distinct steps of TE silencing in different eukaryotes. How such a variety of secondary small RNA biogenesis mechanisms has evolved has not been thoroughly elucidated. Ciliated protozoa perform small RNA-directed programmed DNA elimination of thousands of TE-related internal eliminated sequences (IESs) in the newly developed somatic nucleus. In the ciliate Paramecium, secondary small RNAs are produced after the excision of IESs. In this study, we show that in another ciliate, Tetrahymena, secondary small RNAs accumulate at least a few hours before their derived IESs are excised. We also demonstrate that DNA excision is dispensable for their biogenesis in this ciliate. Therefore, unlike in Paramecium, small RNA amplification occurs before IES excision in Tetrahymena This study reveals the remarkable diversity of secondary small RNA biogenesis mechanisms, even among ciliates with similar DNA elimination processes, and thus raises the possibility that the evolution of TE-targeting small RNA amplification can be traced by investigating the DNA elimination mechanisms of ciliates.

Comparative 3D genome organization in apicomplexan parasites.

The positioning of chromosomes in the nucleus of a eukaryotic cell is highly organized and has a complex and dynamic relationship with gene expression. In the human malaria parasite Plasmodium falciparum, the clustering of a family of virulence genes correlates with their coordinated silencing and has a strong influence on the overall organization of the genome. To identify conserved and species-specific principles of genome organization, we performed Hi-C experiments and generated 3D genome models for five Plasmodium species and two related apicomplexan parasites. Plasmodium species mainly showed clustering of centromeres, telomeres, and virulence genes. In P. falciparum, the heterochromatic virulence gene cluster had a strong repressive effect on the surrounding nuclear space, while this was less pronounced in Plasmodium vivax and Plasmodium berghei, and absent in Plasmodium yoelii In Plasmodium knowlesi, telomeres and virulence genes were more dispersed throughout the nucleus, but its 3D genome showed a strong correlation with gene expression. The Babesia microti genome showed a classical Rabl organization with colocalization of subtelomeric virulence genes, while the Toxoplasma gondii genome was dominated by clustering of the centromeres and lacked virulence gene clustering. Collectively, our results demonstrate that spatial genome organization in most Plasmodium species is constrained by the colocalization of virulence genes. P. falciparum and P. knowlesi, the only two Plasmodium species with gene families involved in antigenic variation, are unique in the effect of these genes on chromosome folding, indicating a potential link between genome organization and gene expression in more virulent pathogens.

Unraveling the 3D genome of human malaria parasites.

The chromosomes within the eukaryotic cell nucleus are highly dynamic and adopt complex hierarchical structures. Understanding how this three-dimensional (3D) nuclear architectureaffects gene regulation, cell cycle progression and disease pathogenesis are important biological questions in development and disease. Recently, many genome-wide technologies including chromosome conformation capture (3C) and 3C-based methodologies (4C, 5C, and Hi-C) have been developed to investigate 3D chromatin structure. In this review, we introduce 3D genome methodologies, with a focus on their application for understanding the nuclear architecture of the human malaria parasite, Plasmodium falciparum. An increasing amount of evidence now suggests that gene regulation in the parasite is largely regulated by epigenetic mechanisms and nuclear reorganization. Here, we explore the 3D genome architecture of P. falciparum, including local and global chromatin structure. In addition, molecular components important for maintaining 3D chromatin organization including architectural proteins and long non-coding RNAs are discussed. Collectively, these studies contribute to our understanding of how the plasticity of 3D genome architecture regulates gene expression and cell cycle progression in this deadly parasite.

Transfection of Capsaspora owczarzaki, a close unicellular relative of animals.

How animals emerged from their unicellular ancestor remains a major evolutionary question. New genome data from the closest unicellular relatives of animals have provided important insights into the evolution of animal multicellularity. We know that the unicellular ancestor of animals had an unexpectedly complex genetic repertoire, including many genes that are key to animal development and multicellularity. Thus, assessing the function of these genes among unicellular relatives of animals is key to understanding how they were co-opted at the onset of the Metazoa. However, such analyses have been hampered by the lack of genetic tools. Progress has been made in choanoflagellates and teretosporeans, two of the three lineages closely related to animals, whereas no tools are yet available for functional analysis in the third lineage: the filastereans. Importantly, filastereans have a striking repertoire of genes involved in transcriptional regulation and other developmental processes. Here, we describe a reliable transfection method for the filasterean Capsaspora owczarzaki We also provide a set of constructs for visualising subcellular structures in live cells. These tools convert Capsaspora into a unique experimentally tractable organism to use to investigate the origin and evolution of animal multicellularity.

Eye-like ocelloids are built from different endosymbiotically acquired components.

Multicellularity is often considered a prerequisite for morphological complexity, as seen in the camera-type eyes found in several groups of animals. A notable exception exists in single-celled eukaryotes called dinoflagellates, some of which have an eye-like 'ocelloid' consisting of subcellular analogues to a cornea, lens, iris, and retina. These planktonic cells are uncultivated and rarely encountered in environmental samples, obscuring the function and evolutionary origin of the ocelloid. Here we show, using a combination of electron microscopy, tomography, isolated-organelle genomics, and single-cell genomics, that ocelloids are built from pre-existing organelles, including a cornea-like layer made of mitochondria and a retinal body made of anastomosing plastids. We find that the retinal body forms the central core of a network of peridinin-type plastids, which in dinoflagellates and their relatives originated through an ancient endosymbiosis with a red alga. As such, the ocelloid is a chimaeric structure, incorporating organelles with different endosymbiotic histories. The anatomical complexity of single-celled organisms may be limited by the components available for differentiation, but the ocelloid shows that pre-existing organelles can be assembled into a structure so complex that it was initially mistaken for a multicellular eye. Although mitochondria and plastids are acknowledged chiefly for their metabolic roles, they can also be building blocks for greater structural complexity.