Reproduction in Trypanosomatids: Past and Present

Detecting complex infections in trypanosomatids using whole genome sequencing

BACKGROUND

Trypanosomatid parasites are a group of protozoans that cause devastating diseases that disproportionately affect developing countries. These protozoans have developed several mechanisms for adaptation to survive in the mammalian host, such as extensive expansion of multigene families enrolled in host-parasite interaction, adaptation to invade and modulate host cells, and the presence of aneuploidy and polyploidy. Two mechanisms might result in "complex" isolates, with more than two haplotypes being present in a single sample: multiplicity of infections (MOI) and polyploidy. We have developed and validated a methodology to identify multiclonal infections and polyploidy using whole genome sequencing reads, based on fluctuations in allelic read depth in heterozygous positions, which can be easily implemented in experiments sequencing genomes from one sample to larger population surveys.

RESULTS

The methodology estimates the complexity index (CI) of an isolate, and compares real samples with simulated clonal infections at individual and populational level, excluding regions with somy and gene copy number variation. It was primarily validated with simulated MOI and known polyploid isolates respectively from Leishmania and Trypanosoma cruzi. Then, the approach was used to assess the complexity of infection using genome wide SNP data from 497 trypanosomatid samples from four clades, L. donovani/L. infantum, L. braziliensis, T. cruzi and T. brucei providing an overview of multiclonal infection and polyploidy in these cultured parasites. We show that our method robustly detects complex infections in samples with at least 25x coverage, 100 heterozygous SNPs and where 5-10% of the reads correspond to the secondary clone. We find that relatively small proportions (≤ 7%) of cultured trypanosomatid isolates are complex.

CONCLUSIONS

The method can accurately identify polyploid isolates, and can identify multiclonal infections in scenarios with sufficient genome read coverage. We pack our method in a single R script that requires only a standard variant call format (VCF) file to run ( https://github.com/jaumlrc/Complex-Infections ). Our analyses indicate that multiclonality and polyploidy do occur in all clades, but not very frequently in cultured trypanosomatids. We caution that our estimates are lower bounds due to the limitations of current laboratory and bioinformatic methods.

Shining the spotlight on the neglected: new high-quality genome assemblies as a gateway to understanding the evolution of Trypanosomatidae

BACKGROUND

Protists of the family Trypanosomatidae (phylum Euglenozoa) have gained notoriety as parasites affecting humans, domestic animals, and agricultural plants. However, the true extent of the group's diversity spreads far beyond the medically and veterinary relevant species. We address several knowledge gaps in trypanosomatid research by undertaking sequencing, assembly, and analysis of genomes from previously overlooked representatives of this protistan group.

RESULTS

We assembled genomes for twenty-one trypanosomatid species, with a primary focus on insect parasites and Trypanosoma spp. parasitizing non-human hosts. The assemblies exhibit sizes consistent with previously sequenced trypanosomatid genomes, ranging from approximately 18 Mb for Obscuromonas modryi to 35 Mb for Crithidia brevicula and Zelonia costaricensis. Despite being the smallest, the genome of O. modryi has the highest content of repetitive elements, contributing nearly half of its total size. Conversely, the highest proportion of unique DNA is found in the genomes of Wallacemonas spp., with repeats accounting for less than 8% of the assembly length. The majority of examined species exhibit varying degrees of aneuploidy, with trisomy being the most frequently observed condition after disomy.

CONCLUSIONS

The genome of Obscuromonas modryi represents a very unusual, if not unique, example of evolution driven by two antidromous forces: i) increasing dependence on the host leading to genomic shrinkage and ii) expansion of repeats causing genome enlargement. The observed variation in somy within and between trypanosomatid genera suggests that these flagellates are largely predisposed to aneuploidy and, apparently, exploit it to gain a fitness advantage. High heterogeneity in the genome size, repeat content, and variation in chromosome copy numbers in the newly-sequenced species highlight the remarkable genome plasticity exhibited by trypanosomatid flagellates. These new genome assemblies are a robust foundation for future research on the genetic basis of life cycle changes and adaptation to different hosts in the family Trypanosomatidae.

Regulatory basis for reproductive flexibility in a meningitis-causing fungal pathogen

Pathogenic fungi of the genus Cryptococcus can undergo two sexual cycles, involving either bisexual diploidization (after fusion of haploid cells of different mating type) or unisexual diploidization (by autodiploidization of a single cell). Here, we construct a gene-deletion library for 111 transcription factor genes in Cryptococcus deneoformans, and explore the roles of these regulatory networks in the two reproductive modes. We show that transcription factors crucial for bisexual syngamy induce the expression of known mating determinants as well as other conserved genes of unknown function. Deletion of one of these genes, which we term FMP1, leads to defects in bisexual reproduction in C. deneoformans, its sister species Cryptococcus neoformans, and the ascomycete Neurospora crassa. Furthermore, we show that a recently evolved regulatory cascade mediates pre-meiotic unisexual autodiploidization, supporting that this reproductive process is a recent evolutionary innovation. Our findings indicate that genetic circuits with different evolutionary ages govern hallmark events distinguishing unisexual and bisexual reproduction in Cryptococcus.

Microevolution and subspecific taxonomy of Trypanosoma cruzi

Trypanosoma cruzi, the agent of Chagas disease, is a highly polymorphic species, subdivided into 6 main evolutionary lineages or near-clades (formerly discrete typing units or DTUs). An additional near-clade (TC-bat) has recently been evidenced. This pattern is considered to be the result of predominant clonal evolution (PCE). PCE is compatible with occasional mating/hybridization, which do not break the prevalent pattern of clonal evolution, the main trait of it being the presence of Multigene Bifurcating Trees (MGBTs) at all evolutionary levels ("clonal frame"). The development of highly resolutive genetic (microsatellites*) and genomic (sequencing and multi-single nucleotide polymorphism {SNP}* typing) markers shows that PCE also operates at a microevolutionary* level within each of the near-clades ("Russian doll pattern"), in spite of occasional meiosis and hybridization events. Within each near-clade, one can evidence widespread clonal multilocus genotypes*, linkage disequilibrium*, Multigene Bifurcating Trees and lesser near-clades. The within near-clade population structure is like a miniature picture of that of the whole species, suggesting gradual rather than saltatory evolution. Additional data are required to evaluate the stability of these lesser near-clades in the long run and to evaluate the need for an adequate nomenclature for this microevolutionary level.

Impact of Genetic Diversity and Genome Plasticity of Leishmania spp. in Treatment and the Search for Novel Chemotherapeutic Targets

Leishmaniasis is one of the major public health concerns in Latin America, Africa, Asia, and Europe. The absence of vaccines for human use and the lack of effective vector control programs make chemotherapy the main strategy to control all forms of the disease. However, the high toxicity of available drugs, limited choice of therapeutic agents, and occurrence of drug-resistant parasite strains are the main challenges related to chemotherapy. Currently, only a small number of drugs are available for leishmaniasis treatment, including pentavalent antimonials (Sb^V), amphotericin B and its formulations, miltefosine, paromomycin sulphate, and pentamidine isethionate. In addition to drug toxicity, therapeutic failure of leishmaniasis is a serious concern. The occurrence of drug-resistant parasites is one of the causes of therapeutic failure and is closely related to the diversity of parasites in this genus. Owing to the enormous plasticity of the genome, resistance can occur by altering different metabolic pathways, demonstrating that resistance mechanisms are multifactorial and extremely complex. Genetic variability and genome plasticity cause not only the available drugs to have limitations, but also make the search for new drugs challenging. Here, we examined the biological characteristics of parasites that hinder drug discovery.

Leishmania and the Model of Predominant Clonal Evolution

As it is the case for other pathogenic microorganisms, the respective impact of clonality and genetic exchange on Leishmania natural populations has been the object of lively debates since the early 1980s. The predominant clonal evolution (PCE) model states that genetic exchange in these parasites’ natural populations may have a high relevance on an evolutionary scale, but is not sufficient to erase a persistent phylogenetic signal and the existence of bifurcating trees. Recent data based on high-resolution markers and genomic polymorphisms fully confirm the PCE model down to a microevolutionary level.

Analyses of Leishmania-LRV Co-Phylogenetic Patterns and Evolutionary Variability of Viral Proteins

Leishmania spp. are important pathogens causing a vector-borne disease with a broad range of clinical manifestations from self-healing ulcers to the life-threatening visceral forms. Presence of Leishmania RNA virus (LRV) confers survival advantage to these parasites by suppressing anti-leishmanial immunity in the vertebrate host. The two viral species, LRV1 and LRV2 infect species of the subgenera Viannia and Leishmania, respectively. In this work we investigated co-phylogenetic patterns of leishmaniae and their viruses on a small scale (LRV2 in L. major) and demonstrated their predominant coevolution, occasionally broken by intraspecific host switches. Our analysis of the two viral genes, encoding the capsid and RNA-dependent RNA polymerase (RDRP), revealed them to be under the pressure of purifying selection, which was considerably stronger for the former gene across the whole tree. The selective pressure also differs between the LRV clades and correlates with the frequency of interspecific host switches. In addition, using experimental (capsid) and predicted (RDRP) models we demonstrated that the evolutionary variability across the structure is strikingly different in these two viral proteins.