High throughput single-cell genome sequencing gives insights into the generation and evolution of mosaic aneuploidy in Leishmania donovani.

Leishmania, a unicellular eukaryotic parasite, is a unique model for aneuploidy and cellular heterogeneity, along with their potential role in adaptation to environmental stresses. Somy variation within clonal populations was previously explored in a small subset of chromosomes using fluorescence hybridization methods. This phenomenon, termed mosaic aneuploidy (MA), might have important evolutionary and functional implications but remains under-explored due to technological limitations. Here, we applied and validated a high throughput single-cell genome sequencing method to study for the first time the extent and dynamics of whole karyotype heterogeneity in two clonal populations of Leishmania promastigotes representing different stages of MA evolution in vitro. We found that drastic changes in karyotypes quickly emerge in a population stemming from an almost euploid founder cell. This possibly involves polyploidization/hybridization at an early stage of population expansion, followed by assorted ploidy reduction. During further stages of expansion, MA increases by moderate and gradual karyotypic alterations, affecting a defined subset of chromosomes. Our data provide the first complete characterization of MA in Leishmania and pave the way for further functional studies.

The kinesin of the flagellum attachment zone in Leishmania is required for cell morphogenesis, cell division and virulence in the mammalian host.

Leishmania parasites possess a unique and complex cytoskeletal structure termed flagellum attachment zone (FAZ) connecting the base of the flagellum to one side of the flagellar pocket (FP), an invagination of the cell body membrane and the sole site for endocytosis and exocytosis. This structure is involved in FP architecture and cell morphogenesis, but its precise role and molecular composition remain enigmatic. Here, we characterized Leishmania FAZ7, the only known FAZ protein containing a kinesin motor domain, and part of a clade of trypanosomatid-specific kinesins with unknown functions. The two paralogs of FAZ7, FAZ7A and FAZ7B, display different localizations and functions. FAZ7A localizes at the basal body, while FAZ7B localizes at the distal part of the FP, where the FAZ structure is present in Leishmania. While null mutants of FAZ7A displayed normal growth rates, the deletion of FAZ7B impaired cell growth in both promastigotes and amastigotes of Leishmania. The kinesin activity is crucial for its function. Deletion of FAZ7B resulted in altered cell division, cell morphogenesis (including flagellum length), and FP structure and function. Furthermore, knocking out FAZ7B induced a mis-localization of two of the FAZ proteins, and disrupted the molecular organization of the FP collar, affecting the localization of its components. Loss of the kinesin FAZ7B has important consequences in the insect vector and mammalian host by reducing proliferation in the sand fly and pathogenicity in mice. Our findings reveal the pivotal role of the only FAZ kinesin as part of the factors important for a successful life cycle of Leishmania.

Universal highly efficient conditional knockout system in Leishmania, with a focus on untranscribed region preservation.

Trypanosomatids are divergent eukaryotes of high medical and economical relevance. Their biology exhibits original features that remain poorly understood; particularly, Leishmania is known for its high degree of genomic plasticity that makes genomic manipulation challenging. CRISPR-Cas9 has been applied successfully to these parasites providing a robust tool to study non-essential gene functions. Here, we have developed a versatile inducible system combining Di-Cre recombinase and CRISPR-Cas9 advantages. Cas9 is used to integrate the LoxP sequences, and the Cre-recombinase catalyses the recombination between LoxP sites, thereby excising the target gene. We used a Leishmania mexicana cell line expressing Di-Cre, Cas9, and T7 polymerase and then transfected donor DNAs and single guide RNAs as polymerase chain reaction (PCR) products. Because the location of LoxP sequences in the genomic DNA can interfere with the function and localisation of certain proteins of interest, we proposed to target the least transcribed regions upstream and/or downstream the gene of interest. To do so, we developed "universal" template plasmids for donor DNA cassettes with or without a tag, where LoxP sequences may be located either immediately upstream the ATG and downstream the stop codon of the gene of interest, or in the least transcribed areas of intergenic regions. Our methodology is fast, PCR-based (molecular cloning-free), highly efficient, versatile, and able to overcome the problems posed by genomic plasticity in Leishmania.

RecQ helicases in the malaria parasite Plasmodium falciparum affect genome stability, gene expression patterns and DNA replication dynamics.

The malaria parasite Plasmodium falciparum has evolved an unusual genome structure. The majority of the genome is relatively stable, with mutation rates similar to most eukaryotic species. However, some regions are very unstable with high recombination rates, driving the generation of new immune evasion-associated var genes. The molecular factors controlling the inconsistent stability of this genome are not known. Here we studied the roles of the two putative RecQ helicases in P. falciparum, PfBLM and PfWRN. When PfWRN was knocked down, recombination rates increased four-fold, generating chromosomal abnormalities, a high rate of chimeric var genes and many microindels, particularly in known 'fragile sites'. This is the first identification of a gene involved in suppressing recombination and maintaining genome stability in Plasmodium. By contrast, no change in mutation rate appeared when the second RecQ helicase, PfBLM, was mutated. At the transcriptional level, however, both helicases evidently modulate the transcription of large cohorts of genes, with several hundred genes-including a large proportion of vars-showing deregulated expression in each RecQ mutant. Aberrant processing of stalled replication forks is a possible mechanism underlying elevated mutation rates and this was assessed by measuring DNA replication dynamics in the RecQ mutant lines. Replication forks moved slowly and stalled at elevated rates in both mutants, confirming that RecQ helicases are required for efficient DNA replication. Overall, this work identifies the Plasmodium RecQ helicases as major players in DNA replication, antigenic diversification and genome stability in the most lethal human malaria parasite, with important implications for genome evolution in this pathogen.

Identification of the centromeres of Leishmania major: revealing the hidden pieces.

Leishmania affects millions of people worldwide. Its genome undergoes constitutive mosaic aneuploidy, a type of genomic plasticity that may serve as an adaptive strategy to survive distinct host environments. We previously found high rates of asymmetric chromosome allotments during mitosis that lead to the generation of such ploidy. However, the underlying molecular events remain elusive. Centromeres and kinetochores most likely play a key role in this process, yet their identification has failed using classical methods. Our analysis of the unconventional kinetochore complex recently discovered in Trypanosoma brucei (KKTs) leads to the identification of a Leishmania KKT gene candidate (LmKKT1). The GFP-tagged LmKKT1 displays "kinetochore-like" dynamics of intranuclear localization throughout the cell cycle. By ChIP-Seq assay, one major peak per chromosome is revealed, covering a region of 4 ±2 kb. We find two largely conserved motifs mapping to 14 of 36 chromosomes while a higher density of retroposons are observed in 27 of 36 centromeres. The identification of centromeres and of a kinetochore component of Leishmania chromosomes opens avenues to explore their role in mosaic aneuploidy.

The nucleoporin Mlp2 is involved in chromosomal distribution during mitosis in trypanosomatids.

Nucleoporins are evolutionary conserved proteins mainly involved in the constitution of the nuclear pores and trafficking between the nucleus and cytoplasm, but are also increasingly viewed as main actors in chromatin dynamics and intra-nuclear mitotic events. Here, we determined the cellular localization of the nucleoporin Mlp2 in the 'divergent' eukaryotes Leishmania major and Trypanosoma brucei. In both protozoa, Mlp2 displayed an atypical localization for a nucleoporin, essentially intranuclear, and preferentially in the periphery of the nucleolus during interphase; moreover, it relocated at the mitotic spindle poles during mitosis. In T. brucei, where most centromeres have been identified, TbMlp2 was found adjacent to the centromeric sequences, as well as to a recently described unconventional kinetochore protein, in the periphery of the nucleolus, during interphase and from the end of anaphase onwards. TbMlp2 and the centromeres/kinetochores exhibited a differential migration towards the poles during mitosis. RNAi knockdown of TbMlp2 disrupted the mitotic distribution of chromosomes, leading to a surprisingly well-tolerated aneuploidy. In addition, diploidy was restored in a complementation assay where LmMlp2, the orthologue of TbMlp2 in Leishmania, was expressed in TbMlp2-RNAi-knockdown parasites. Taken together, our results demonstrate that Mlp2 is involved in the distribution of chromosomes during mitosis in trypanosomatids.

First efficient CRISPR-Cas9-mediated genome editing in Leishmania parasites.

Protozoan pathogens that cause leishmaniasis in humans are relatively refractory to genetic manipulation. In this work, we implemented the CRISPR-Cas9 system in Leishmania parasites and demonstrated its efficient use for genome editing. The Cas9 endonuclease was expressed under the control of the Dihydrofolate Reductase-Thymidylate Synthase (DHFR-TS) promoter and the single guide RNA was produced under the control of the U6snRNA promoter and terminator. As a proof of concept, we chose to knockout a tandemly repeated gene family, the paraflagellar rod-2 locus. We were able to obtain null mutants in a single round of transfection. In addition, we confirmed the absence of off-target editions by whole genome sequencing of two independent clones. Our work demonstrates that CRISPR-Cas9-mediated gene knockout represents a major improvement in comparison with existing methods. Beyond gene knockout, this genome editing tool opens avenues for a multitude of functional studies to speed up research on leishmaniasis.

Single-molecule analysis reveals that DNA replication dynamics vary across the course of schizogony in the malaria parasite Plasmodium falciparum.

The mechanics of DNA replication and cell cycling are well-characterized in model organisms, but less is known about these basic aspects of cell biology in early-diverging Apicomplexan parasites, which do not divide by canonical binary fission but undergo unconventional cycles. Schizogony in the malaria parasite, Plasmodium, generates ~16-24 new nuclei via independent, asynchronous rounds of genome replication prior to cytokinesis and little is known about the control of DNA replication that facilitates this. We have characterised replication dynamics in P. falciparum throughout schizogony, using DNA fibre labelling and combing to visualise replication forks at a single-molecule level. We show that origins are very closely spaced in Plasmodium compared to most model systems, and that replication dynamics vary across the course of schizogony, from faster synthesis rates and more widely-spaced origins through to slower synthesis rates and closer-spaced origins. This is the opposite of the pattern usually seen across S-phase in human cells, when a single genome is replicated. Replication forks also appear to stall at an unusually high rate throughout schizogony. Our work explores Plasmodium DNA replication in unprecedented detail and opens up tremendous scope for analysing cell cycle dynamics and developing interventions targetting this unique aspect of malaria biology.

Single-molecule analysis of DNA replication reveals novel features in the divergent eukaryotes Leishmania and Trypanosoma brucei versus mammalian cells.

Leishmania and Trypanosoma are unicellular parasites that possess markedly original biological features as compared to other eukaryotes. The Leishmania genome displays a constitutive 'mosaic aneuploidy', whereas in Trypanosoma brucei, the megabase-sized chromosomes are diploid. We accurately analysed DNA replication parameters in three Leishmania species and Trypanosoma brucei as well as mouse embryonic fibroblasts (MEF). Active replication origins were visualized at the single molecule level using DNA molecular combing. More than one active origin was found on most DNA fibres, showing that the chromosomes are replicated from multiple origins. Inter-origin distances (IODs) were measured and found very large in trypanosomatids: the mean IOD was 160 kb in T. brucei and 226 kb in L. mexicana. Moreover, the progression of replication forks was faster than in any other eukaryote analyzed so far (mean velocity 1.9 kb/min in T. brucei and 2.4-2.6 kb/min in Leishmania). The estimated total number of active DNA replication origins in trypanosomatids is ~170. Finally, 14.4% of unidirectional replication forks were observed in T. brucei, in contrast to 1.5-1.7% in Leishmania and 4% in MEF cells. The biological significance of these original features is discussed.

Constitutive mosaic aneuploidy is a unique genetic feature widespread in the Leishmania genus.

Using fluorescence in situ hybridization, we determined the ploidy of four species of Leishmania: Leishmania infantum, Leishmania donovani, Leishmania tropica and Leishmania amazonensis. We found that each cell in a strain possesses a combination of mono-, di- and trisomies for all chromosomes; ploidy patterns were different among all strains/species. These results extend those we previously described in Leishmania major, demonstrating that mosaic aneuploidy is a genetic feature widespread to the Leishmania genus. In addition to the genetic consequences induced by this mosaicism, the apparent absence of alternation between haploid/diploid stages questions the modality of genetic exchange in Leishmania sp.