Genome instability drives epistatic adaptation in the human pathogen Leishmania

R-loops acted on by RNase H1 influence DNA replication timing and genome stability in Leishmania

Genomes in eukaryotes normally undergo DNA replication in a choreographed temporal order, resulting in early and late replicating chromosome compartments. Leishmania, a human protozoan parasite, displays an unconventional DNA replication program in which the timing of DNA replication completion is chromosome size-dependent: larger chromosomes complete replication later then smaller ones. Here we show that both R-loops and RNase H1, a ribonuclease that resolves RNA-DNA hybrids, accumulate in Leishmania major chromosomes in a pattern that reflects their replication timing. Furthermore, we demonstrate that such differential organisation of R-loops, RNase H1 and DNA replication timing across the parasite's chromosomes correlates with size-dependent differences in chromatin accessibility, G quadruplex distribution and sequence content. Using conditional gene excision, we show that loss of RNase H1 leads to transient growth perturbation and permanently abrogates the differences in DNA replication timing across chromosomes, as well as altering levels of aneuploidy and increasing chromosome instability in a size-dependent manner. This work provides a link between R-loop homeostasis and DNA replication timing in a eukaryotic parasite and demonstrates that orchestration of DNA replication dictates levels of genome plasticity in Leishmania.

Nanopore sequencing reveals that DNA replication compartmentalisation dictates genome stability and instability in Trypanosoma brucei

The Trypanosoma brucei genome is structurally complex. Eleven megabase-sized chromosomes each comprise a transcribed core flanked by silent subtelomeres, housing thousands of Variant Surface Glycoprotein (VSG) genes. Additionally, hundreds of sub-megabase chromosomes contain 177 bp repeats of unknown function, and VSG transcription sites localise to many telomeres. DNA replication dynamics have only been described in the megabase chromosome cores, and in the single active VSG transcription site. Using a Nanopore genome assembly, we show that megabase chromosome subtelomeres display a paucity of replication initiation events relative to the core, correlating with increased instability. In addition, replication of the active VSG transcription site is shown to originate from the telomere, likely causing targeted VSG recombination. Lastly, we provide evidence that the 177 bp repeats act as conserved DNA replication origins, explaining submegabase chromosome stability. Compartmentalized DNA replication therefore explains how T. brucei balances stable genome transmission with localised instability driving immune evasion.

Transcriptomic profiling of "brain-eating amoeba" Naegleria fowleri infection in mice: the host and the protozoa perspectives

The free-living amoeba Naegleria fowleri (NF) causes a rare but lethal parasitic meningoencephalitis (PAM) in humans. Currently, this disease lacks effective treatments and the specific molecular mechanisms that govern NF pathogenesis and host brain response remain unknown. To address some of these issues, we sought to explore naturally existing virulence diversity within environmental NF isolates. Herein, we purified two new NF environmental isolates (NF45 and NF1) and tested their in vivo virulence using experimental infection in mice. We found that NF45 was highly virulent (NF45_HV) compared with NF1 (low virulence, NF1_LV), based on in vivo amoeba growth kinetics and mouse survival. To identify underlying differences, we conducted RNA-seq and bioinformatics analyses from the infected mouse brains. Our results showed that NF1_LV and NF45_HV modulated the expression of their genes during mouse brain infection. Differentially expressed genes (DEGs) in NF1_LV were mostly involved in Translational protein, Protein-binding activity modulator, Protein modifying enzyme, while DEGs in NF45_HV were related to DNA metabolism, Cytoskeletal protein, Protein-binding activity modulator. Proteases (namely the virulence factor Cathepsin B) were upregulated in NF1_LV, while downregulated in NF45_HV. When analyzing the host response against infection by these two NF strains, enrichment analyses uncovered genes and mechanisms related to the host immune responses and nervous systems. We detected more DEGs in NF1_LV infected mice compared to NF45_HV, related to blood brain barrier leakage, immune cell recruitment, cytokine production (including IL-6, IFN-Ɣ and TNFα), inflammation of astrocytes and microglia, and oligodendrocyte and neurons degeneration. Increased expression of neuromotor-related genes such as Adam22, Cacnb4 and Zic1 (activated by NF1_LV infection) and ChAt (activated by NF45_LV infection) could explain PAM symptoms such as muscle weakness and seizures. Globally, our results showed that NF isolated from the environment can have different levels of virulence and differentially modulate their gene expression during brain infection. We also provided, for the first time, a comprehensive information for the molecular mechanisms of neuro-immune and host-pathogen interactions during PAM disease. As the host and the protozoa are strongly implicated in PAM lethality, new therapies targeting both the parasite, and the host should be considered to treat PAM infection.

Mechanisms of life cycle simplification in African trypanosomes

African trypanosomes are important parasites in sub-Saharan Africa that undergo a quorum-sensing dependent development to morphologically 'stumpy forms' in mammalian hosts to favour transmission by tsetse flies. However, some trypanosome clades have simplified their lifecycle by escaping dependence on tsetse allowing an expanded geographic range, with direct transmission between hosts achieved via blood-feeding biting flies and vampire bats (Trypanosoma brucei evansi, causing 'surra') or through sexual transmission (Trypanosoma brucei equiperdum, causing 'dourine'). Concomitantly, stumpy formation is reduced and the isolates are described as monomorphic, with infections spread widely in Africa, Asia, South America and parts of Europe. Here, using genomic analysis of distinct field isolates, we identify molecular changes that accompany the loss of the stumpy formation in monomorphic clades. Using CRISPR-mediated allelic replacement, mutations in two exemplar genes (Tb927.2.4020; Tb927.5.2580) are confirmed to reduce stumpy formation whereas another (Tb927.11.3400) is implicated in altered motility. Using laboratory selection we identify downregulation of RNA regulators as important in the initial development of monomorphism. This identifies a trajectory of events that simplify the life cycle in emergent and established monomorphic trypanosomes, with impact on disease spread, vector control strategies, geographical range and virulence.

Molecular Mechanisms of Drug Resistance in Leishmania spp

The protozoan parasite Leishmania causes leishmaniasis, a neglected tropical disease, that disproportionately affects underdeveloped countries. This disease has major health, economic, and social implications, particularly because of the limited treatment options, high cost, the severe side effects associated with available therapeutics, and the high rate of treatment failure caused by the parasites' growing resistance to current medications. In this review, we describe first the common strategies used by pathogens to develop drug resistance and then focus on the arsenal of available drugs to treat leishmaniasis, their modes of action, and the molecular mechanisms contributing to drug resistance in Leishmania spp., including the role of genomic, transcriptional, and translational control. We focus more specifically on our recent discovery of translational reprogramming as a major driver of drug resistance leading to coordinated changes in the translation of transcripts and orchestrating changes in metabolome and lipidome to support drug resistance. A thorough understanding of these mechanisms is essential to identify the key elements needed to combat resistance and improve leishmaniasis treatment methods.

Genetic Variability in Leishmaniasis-Causing Leishmania infantum in Humans and Dogs from North-East Spain

Leishmania infantum is the primary cause of visceral and cutaneous leishmaniasis in the European Mediterranean region. Subspecies-level characterization of L. infantum aids epidemiological studies by offering insights into the evolution and geographical distribution of the parasite and reservoir identity. In this study, conducted in north-east Spain, 26 DNA samples of L. infantum were analyzed, comprising 21 from 10 humans and 5 from 5 dogs. Minicircle kinetoplast DNA (kDNA) polymerase chain reaction assays using primers MC1 and MC2, followed by sequencing, were employed to assess intraspecific genetic variability. Single-nucleotide polymorphism (SNP) analysis detected seven genotypes (G1, G2, G12*-G15*, and G17*), with five being reported for the first time (*). The most prevalent was the newly described G13 (54%), while the other currently identified genotypes were predominantly found in single samples. The in silico restriction fragment length polymorphism (RFLP) method revealed five genotypes (B, F, N, P, and W), one of them previously unreported (W). Genotype B was the most prevalent (85%), comprising three SNP genotypes (G1, G2, and G13), whereas the other RFLP genotypes were associated with single SNP genotypes. These kDNA genotyping methods revealed significant intraspecific genetic diversity in L. infantum, demonstrating their suitability for fingerprinting and strain monitoring.

Structural and mechanistic insights into the function of Leishmania ribosome lacking a single pseudouridine modification

Leishmania is the causative agent of cutaneous and visceral diseases affecting millions of individuals worldwide. Pseudouridine (Ψ), the most abundant modification on rRNA, changes during the parasite life cycle. Alterations in the level of a specific Ψ in helix 69 (H69) affected ribosome function. To decipher the molecular mechanism of this phenotype, we determine the structure of ribosomes lacking the single Ψ and its parental strain at ∼2.4-3 Å resolution using cryo-EM. Our findings demonstrate the significance of a single Ψ on H69 to its structure and the importance for its interactions with helix 44 and specific tRNAs. Our study suggests that rRNA modification affects translation of mRNAs carrying codon bias due to selective accommodation of tRNAs by the ribosome. Based on the high-resolution structures, we propose a mechanism explaining how the ribosome selects specific tRNAs.

A comparative genomics approach reveals a local genetic signature of Leishmania tropica in Morocco

In Morocco, cutaneous leishmaniasis (CL) caused by Leishmania (L.) tropica is an important health problem. Despite the high incidence of CL in the country, the genomic heterogeneity of these parasites is still incompletely understood. In this study, we sequenced the genomes of 14 Moroccan isolates of L. tropica collected from confirmed cases of CL to investigate their genomic heterogeneity. Comparative genomics analyses were conducted by applying the recently established Genome Instability Pipeline (GIP), which allowed us to conduct phylogenomic and principal components analyses (PCA), and to assess genomic variations at the levels of the karyotype, gene copy number, single nucleotide polymorphisms (SNPs) and small insertions/deletions (INDELs) variants. Read-depth analyses revealed a mostly disomic karyotype, with the exception of the stable tetrasomy of chromosome 31. In contrast, we identified important gene copy number variations across all isolates, which affect known virulence genes and thus were probably selected in the field. SNP-based cluster analysis of the 14 isolates revealed a core group of 12 strains that formed a tight cluster and shared 45.1 % (87 751) of SNPs, as well as two strains (M3015, Ltr_16) that clustered separately from each other and the core group, suggesting the circulation of genetically highly diverse strains in Morocco. Phylogenetic analysis, which compared our 14 L. tropica isolates against 40 published genomes of L. tropica from a diverse array of locations, confirmed the genetic difference of our Moroccan isolates from all other isolates examined. In conclusion, our results indicate potential regional variations in SNP profiles that may differentiate Moroccan L. tropica from other L. tropica strains circulating in endemic countries in the Middle East. Our report paves the way for future research with a larger number of strains that will allow correlation of diverse phenotypes (resistance to treatments, virulence) and origins (geography, host species, year of isolation) to defined genomic signals such as gene copy number variations or SNP profiles that may represent interesting biomarker candidates.

Histone H3.1 is a chromatin-embedded redox sensor triggered by tumor cells developing adaptive phenotypic plasticity and multidrug resistance

Chromatin structure is regulated through posttranslational modifications of histone variants that modulate transcription. Although highly homologous, histone variants display unique amino acid sequences associated with specific functions. Abnormal incorporation of histone variants contributes to cancer initiation, therapy resistance, and metastasis. This study reports that, among its biologic functions, histone H3.1 serves as a chromatin redox sensor that is engaged by mitochondrial H2O2. In breast cancer cells, the oxidation of H3.1Cys96 promotes its eviction and replacement by H3.3 in specific promoters. We also report that this process facilitates the opening of silenced chromatin domains and transcriptional activation of epithelial-to-mesenchymal genes associated with cell plasticity. Scavenging nuclear H2O2 or amino acid substitution of H3.1(C96S) suppresses plasticity, restores sensitivity to chemotherapy, and induces remission of metastatic lesions. Hence, it appears that increased levels of H2O2 produced by mitochondria of breast cancer cells directly promote redox-regulated H3.1-dependent chromatin remodeling involved in chemoresistance and metastasis.

A single pseudouridine on rRNA regulates ribosome structure and function in the mammalian parasite Trypanosoma brucei

Trypanosomes are protozoan parasites that cycle between insect and mammalian hosts and are the causative agent of sleeping sickness. Here, we describe the changes of pseudouridine (Ψ) modification on rRNA in the two life stages of the parasite using four different genome-wide approaches. CRISPR-Cas9 knock-outs of all four snoRNAs guiding Ψ on helix 69 (H69) of the large rRNA subunit were lethal. A single knock-out of a snoRNA guiding Ψ530 on H69 altered the composition of the 80S monosome. These changes specifically affected the translation of only a subset of proteins. This study correlates a single site Ψ modification with changes in ribosomal protein stoichiometry, supported by a high-resolution cryo-EM structure. We propose that alteration in rRNA modifications could generate ribosomes preferentially translating state-beneficial proteins.

The adaptive roles of aneuploidy and polyclonality in Leishmania in response to environmental stress

Aneuploidy is generally considered harmful, but in some microorganisms, it can act as an adaptive mechanism against environmental stress. Here, we use Leishmania-a protozoan parasite with remarkable genome plasticity-to study the early steps of aneuploidy evolution under high drug pressure (using antimony or miltefosine as stressors). By combining single-cell genomics, lineage tracing with cellular barcodes, and longitudinal genome characterization, we reveal that aneuploidy changes under antimony pressure result from polyclonal selection of pre-existing karyotypes, complemented by further and rapid de novo alterations in chromosome copy number along evolution. In the case of miltefosine, early parasite adaptation is associated with independent point mutations in a miltefosine transporter gene, while aneuploidy changes only emerge later, upon exposure to increased drug levels. Therefore, polyclonality and genome plasticity are hallmarks of parasite adaptation, but the scenario of aneuploidy dynamics depends on the nature and strength of the environmental stress as well as on the existence of other pre-adaptive mechanisms.

Infectivity and stress tolerance traits affect community assembly of plant pathogenic fungi

Understanding how ecological communities assemble is an urgent research priority. In this study, we used a community ecology approach to examine how ecological and evolutionary processes shape biodiversity patterns of plant pathogenic fungi, Fusarium graminearum and F. asiaticum. High-throughput screening revealed that the isolates had a wide range of phenotypic variation in stress tolerance traits. Net Relatedness Index (NRI) and Nearest Taxon Index (NTI) values were computed based on stress-tolerant distance matrices. Certain local regions exhibited positive values of NRI and NTI, indicating phenotypic clustering within the fungal communities. Competition assays of the pooled strains were conducted to investigate the cause of clustering. During stress conditions and wheat colonization, only a few strains dominated the fungal communities, resulting in reduced diversity. Overall, our findings support the modern coexistence theory that abiotic stress and competition lead to phenotypic similarities among coexisting organisms by excluding large, low-competitive clades. We suggest that agricultural environments and competition for host infection lead to locally clustered communities of plant pathogenic fungi in the field.

Kinetoplast Genome of Leishmania spp. Is under Strong Purifying Selection

Instability is an intriguing characteristic of many protist genomes, and trypanosomatids are not an exception in this respect. Some regions of trypanosomatid genomes evolve fast. For instance, the trypanosomatid mitochondrial (kinetoplast) genome consists of fairly conserved maxicircle and minicircle molecules that can, nevertheless, possess high nucleotide substitution rates between closely related strains. Recent experiments have demonstrated that rapid laboratory evolution can result in the non-functionality of multiple genes of kinetoplast genomes due to the accumulation of mutations or loss of critical genomic components. An example of a loss of critical components is the reported loss of entire minicircle classes in Leishmania tarentolae during laboratory cultivation, which results in an inability to generate some correctly encoded genes. In the current work, we estimated the evolutionary rates of mitochondrial and nuclear genome regions of multiple natural Leishmania spp. We analyzed synonymous and non-synonymous substitutions and, rather unexpectedly, found that the coding regions of kinetoplast maxicircles are among the most variable regions of both genomes. In addition, we demonstrate that synonymous substitutions greatly predominate among maxicircle coding regions and that most maxicircle genes show signs of purifying selection. These results imply that maxicircles in natural Leishmania populations remain functional despite their high mutation rate.

Leishmania allelic selection during experimental sand fly infection correlates with mutational signatures of oxidative DNA damage

Trypanosomatid pathogens are transmitted by blood-feeding insects, causing devastating human infections. These parasites show important phenotypic shifts that often impact parasite pathogenicity, tissue tropism, or drug susceptibility. The evolutionary mechanisms that allow for the selection of such adaptive phenotypes remain only poorly investigated. Here, we use Leishmania donovani as a trypanosomatid model pathogen to assess parasite evolutionary adaptation during experimental sand fly infection. Comparing the genome of the parasites before and after sand fly infection revealed a strong population bottleneck effect as judged by allele frequency analysis. Apart from random genetic drift caused by the bottleneck effect, our analyses revealed haplotype and allelic changes during sand fly infection that seem under natural selection given their convergence between independent biological replicates. Our analyses further uncovered signature mutations of oxidative DNA damage in the parasite genomes after sand fly infection, suggesting that Leishmania suffers from oxidative stress inside the insect digestive tract. Our results propose a model of Leishmania genomic adaptation during sand fly infection, with oxidative DNA damage and DNA repair processes likely driving haplotype and allelic selection. The experimental and computational framework presented here provides a useful blueprint to assess evolutionary adaptation of other eukaryotic pathogens inside their insect vectors, such as Plasmodium spp, Trypanosoma brucei, and Trypanosoma cruzi.

Life in plastic, it's fantastic! How Leishmania exploit genome instability to shape gene expression

Leishmania are kinetoplastid pathogens that cause leishmaniasis, a debilitating and potentially life-threatening infection if untreated. Unusually, Leishmania regulate their gene expression largely post-transcriptionally due to the arrangement of their coding genes into polycistronic transcription units that may contain 100s of functionally unrelated genes. Yet, Leishmania are capable of rapid and responsive changes in gene expression to challenging environments, often instead correlating with dynamic changes in their genome composition, ranging from chromosome and gene copy number variations to the generation of extrachromosomal DNA and the accumulation of point mutations. Typically, such events indicate genome instability in other eukaryotes, coinciding with genetic abnormalities, but for Leishmania, exploiting these products of genome instability can provide selectable substrates to catalyse necessary gene expression changes by modifying gene copy number. Unorthodox DNA replication, DNA repair, replication stress factors and DNA repeats are recognised in Leishmania as contributors to this intrinsic instability, but how Leishmania regulate genome plasticity to enhance fitness whilst limiting toxic under- or over-expression of co-amplified and co-transcribed genes is unclear. Herein, we focus on fresh, and detailed insights that improve our understanding of genome plasticity in Leishmania. Furthermore, we discuss emerging models and factors that potentially circumvent regulatory issues arising from polycistronic transcription. Lastly, we highlight key gaps in our understanding of Leishmania genome plasticity and discuss future studies to define, in higher resolution, these complex regulatory interactions.

Genetic slippage after sex maintains diversity for parasite resistance in a natural host population

Although parasite-mediated selection is a major driver of host evolution, its influence on genetic variation for parasite resistance is not yet well understood. We monitored resistance in a large population of the planktonic crustacean Daphnia magna over 8 years, as it underwent yearly epidemics of the bacterial pathogen Pasteuria ramosa. We observed cyclic dynamics of resistance: Resistance increased throughout the epidemics, but susceptibility was restored each spring when hosts hatched from sexual resting stages. Host resting stages collected across the year showed that largely resistant host populations can produce susceptible sexual offspring. A genetic model of resistance developed for this host-parasite system, based on multiple loci and strong epistasis, is in partial agreement with our findings. Our results reveal that, despite strong selection for resistance in a natural host population, genetic slippage after sexual reproduction can be a strong factor for the maintenance of genetic diversity of host resistance.

A feedback mechanism controls rDNA copy number evolution in yeast independently of natural selection

Ribosomal DNA (rDNA) is the genetic loci that encodes rRNA in eukaryotes. It is typically arranged as tandem repeats that vary in copy number within the same species. We have recently shown that rDNA repeats copy number in the yeast Saccharomyces cerevisiae is controlled by cell volume via a feedback circuit that senses cell volume by means of the concentration of the free upstream activator factor (UAF). The UAF strongly binds the rDNA gene promoter, but is also able to repress SIR2 deacetylase gene transcription that, in turn, represses rDNA amplification. In this way, the cells with a smaller DNA copy number than what is optimal evolve to increase that copy number until they reach a number that sequestrates free UAF and provokes SIR2 derepression that, in turn, blocks rDNA amplification. Here we propose a mathematical model to show that this evolutionary process can amplify rDNA repeats independently of the selective advantage of yeast cells having bigger or smaller rDNA copy numbers. We test several variants of this process and show that it can explain the observed experimental results independently of natural selection. These results predict that an autoregulated feedback circuit may, in some instances, drive to non Darwinian deterministic evolution for a limited time period.