dgfr: an R package to assess sequence diversity of gene families.

BACKGROUND: Gene families are groups of homologous genes that often have similar biological functions. These families are formed by gene duplication events throughout evolution, resulting in multiple copies of an ancestral gene. Over time, these copies can acquire mutations and structural variations, resulting in members that may vary in size, motif ordering and sequence. Multigene families have been described in a broad range of organisms, from single-celled bacteria to complex multicellular organisms, and have been linked to an array of phenomena, such as host-pathogen interactions, immune evasion and embryonic development. Despite the importance of gene families, few approaches have been developed for estimating and graphically visualizing their diversity patterns and expression profiles in genome-wide studies. RESULTS: Here, we introduce an R package named dgfr, which estimates and enables the visualization of sequence divergence within gene families, as well as the visualization of secondary data such as gene expression. The package takes as input a multi-fasta file containing the coding sequences (CDS) or amino acid sequences from a multigene family, performs a pairwise alignment among all sequences, and estimates their distance, which is subjected to dimension reduction, optimal cluster determination, and gene assignment to each cluster. The result is a dataset that allows for the visualization of sequence divergence and expression within the gene family, an approximation of the number of clusters present in the family. CONCLUSIONS: dgfr provides a way to estimate and study the diversity of gene families, as well as visualize the dispersion and secondary profile of the sequences. The dgfr package is available at https://github.com/lailaviana/dgfr under the GPL-3 license.

Effects of in ovo injection of bacterial peptides and CpG-ODN on Salmonella enterica serovar Heidelberg infection in specific pathogen-free (SPF) chicks.

RESEARCH HIGHLIGHTS: Peptides + CpG-ODN reduced SH in caeca at the first week post-infection.Administered formulations did not reduce SH-faecal excretion.Levels of intestinal IgA were similar between all groups.CpG-ODN improved some parameters associated with chick intestinal health.

Complete assembly, annotation of virulence genes and CRISPR editing of the genome of Leishmania amazonensis PH8 strain.

We report the sequencing and assembly of the PH8 strain of Leishmania amazonensis one of the etiological agents of leishmaniasis. After combining data from long Pacbio reads, short Illumina reads and synteny with the Leishmania mexicana genome, the sequence of 34 chromosomes with 8317 annotated genes was generated. Multigene families encoding three virulence factors, A2, amastins and the GP63 metalloproteases, were identified and compared to their annotation in other Leishmania species. As they have been recently recognized as virulence factors essential for disease establishment and progression of the infection, we also identified 14 genes encoding proteins involved in parasite iron and heme metabolism and compared to genes from other Trypanosomatids. To follow these studies with a genetic approach to address the role of virulence factors, we tested two CRISPR-Cas9 protocols to generate L. amazonensis knockout cell lines, using the Miltefosine transporter gene as a proof of concept.

Assessing the composition of the plasma membrane of Leishmania (Leishmania) infantum and L. (L.) amazonensis using label-free proteomics.

Protozoan parasites of the genus Leishmania are causative agents of leishmaniasis, a wide range of diseases affecting 12 million people worldwide. The species L. infantum and L. amazonensis are etiologic agents of visceral and cutaneous leishmaniasis, respectively. Most proteome analyses of Leishmania have been carried out on whole-cell extracts, but such an approach tends to underrepresent membrane-associated proteins due to their high hydrophobicity and low solubility. Considering the relevance of this category of proteins in virulence, invasiveness and the host-parasite interface, this study applied label-free proteomics to assess the plasma membrane sub-proteome of L. infantum and L. amazonensis. The number of proteins identified in L. infantum and L. amazonensis promastigotes was 1168 and 1455, respectively. After rigorous data processing and mining, 157 proteins were classified as putative plasma membrane-associated proteins, of which 56 proteins were detected in both species, six proteins were detected only in L. infantum and 39 proteins were exclusive to L. amazonensis. The quantitative analysis revealed that two proteins were more abundant in L. infantum, including the glucose transporter 2, and five proteins were more abundant in L. amazonensis. The identified proteins associated with distinct processes and functions. In this regard, proteins of L. infantum were linked to metabolic processes whereas L. amazonensis proteins were involved in signal transduction. Moreover, transmembrane transport was a significant process among the group of proteins detected in both species and members of the superfamily of ABC transporters were highly represented. Interestingly, some proteins of this family were solely detected in L. amazonensis, such as ABCA9. GP63, a well-known virulence factor, was the only GPI-anchored protein identified in the membrane preparations of both species. Finally, we found several proteins with uncharacterized functions, including differentially abundant ones, highlighting a gap in the study of Leishmania proteins. Proteins characterization could provide a better biological understanding of these parasites and deliver new possibilities regarding the discovery of therapeutic targets, drug resistance and vaccine candidates.

Efficacy of sulfadiazine and pyrimetamine for treatment of experimental toxoplasmosis with strains obtained from human cases of congenital disease in Brazil.

Toxoplasmosis in South America presents great health impacts and is a topic of research interest not only because of the severity of native cases but also due to the predominant atypical genotypes of the parasite circulating in this continent. Typically, symptomatic toxoplasmosis is treated with a combination of sulfadiazine (SDZ) and pyrimethamine (PYR). However, some clinical cases present treatment failures due to an inability of the drugs to control the infection or their significant adverse effects, which can lead to treatment interruption. Although resistance/susceptibility to the aforementioned drugs has been well described for clonal strains of Toxoplasma gondii spread to the Northern Hemisphere, less is known about the South American atypical strains. In this study, the effectiveness of SDZ and PYR for the treatment of mice during acute infection with different atypical T. gondii strains was evaluated. Swiss mice were infected with seven T. gondii strains obtained from newborn patients with congenital toxoplasmosis in Brazil. The infected mice were treated with 10-640 mg/kg per day of SDZ, 3-200 mg/kg per day of PYR, or a combination of both drugs with a lower dosage. The mice were evaluated for parameters including mortality, anti-T. gondii IgG production by ELISA and the presence of brain cysts. In addition, the presence of polymorphisms in the dhps gene was verified by gene sequencing. A descriptive analysis was used to assess the association between susceptibility to SDZ and/or PYR and the genotype. The TgCTBr4 and TgCTBr17 strains (genotype 108) presented lower susceptibility to SDZ or PYR treatment. The TgCTBr1 and TgCTBr25 strains (genotype 206) presented similar susceptibility to PYR but not SDZ treatment. The TgCTBr9 strain (genotype 11) was the only strain with high susceptibility to treatment with both drugs. The TgCTBr13 strain (genotype 208) was not susceptible to treatment with the lower PYR or SDZ doses. The TgCTBR23 strain (genotype 41) was more susceptible to PYR than to SDZ treatment. However, the association of low SDZ and PYR doses showed good efficacy for the treatment of experimental toxoplasmosis with T. gondii atypical strains obtained from newborns in Brazil. A new mutation in the T. gondii dhps gene (I347M) was identified that might be associated with the SDZ low sensitivity profile observed for the TgCTBr4 and TgCTBr17 isolates.

Whole genome sequencing of Trypanosoma cruzi field isolates reveals extensive genomic variability and complex aneuploidy patterns within TcII DTU.

BACKGROUND: Trypanosoma cruzi, the etiologic agent of Chagas disease, is currently divided into six discrete typing units (DTUs), named TcI-TcVI. TcII is among the major DTUs enrolled in human infections in South America southern cone, where it is associated with severe cardiac and digestive symptoms. Despite the importance of TcII in Chagas disease epidemiology and pathology, so far, no genome-wide comparisons of the mitochondrial and nuclear genomes of TcII field isolates have been performed to track the variability and evolution of this DTU in endemic regions. RESULTS: In the present work, we have sequenced and compared the whole nuclear and mitochondrial genomes of seven TcII strains isolated from chagasic patients from the central and northeastern regions of Minas Gerais, Brazil, revealing an extensive genetic variability within this DTU. A comparison of the phylogeny based on the nuclear or mitochondrial genomes revealed that the majority of branches were shared by both sequences. The subtle divergences in the branches are probably consequence of mitochondrial introgression events between TcII strains. Two T. cruzi strains isolated from patients living in the central region of Minas Gerais, S15 and S162a, were clustered in the nuclear and mitochondrial phylogeny analysis. These two strains were isolated from the other five by the Espinhaço Mountains, a geographic barrier that could have restricted the traffic of insect vectors during T. cruzi evolution in the Minas Gerais state. Finally, the presence of aneuploidies was evaluated, revealing that all seven TcII strains have a different pattern of chromosomal duplication/loss. CONCLUSIONS: Analysis of genomic variability and aneuploidies suggests that there is significant genomic variability within Minas Gerais TcII strains, which could be exploited by the parasite to allow rapid selection of favorable phenotypes. Also, the aneuploidy patterns vary among T. cruzi strains and does not correlate with the nuclear phylogeny, suggesting that chromosomal duplication/loss are recent and frequent events in the parasite evolution.

Gene and Chromosomal Copy Number Variations as an Adaptive Mechanism Towards a Parasitic Lifestyle in Trypanosomatids.

Trypanosomatids are a group of kinetoplastid parasites including some of great public health importance, causing debilitating and life-long lasting diseases that affect more than 24 million people worldwide. Among the trypanosomatids, Trypanosoma cruzi, Trypanosoma brucei and species from the Leishmania genus are the most well studied parasites, due to their high prevalence in human infections. These parasites have an extreme genomic and phenotypic variability, with a massive expansion in the copy number of species-specific multigene families enrolled in host-parasite interactions that mediate cellular invasion and immune evasion processes. As most trypanosomatids are heteroxenous, and therefore their lifecycles involve the transition between different hosts, these parasites have developed several strategies to ensure a rapid adaptation to changing environments. Among these strategies, a rapid shift in the repertoire of expressed genes, genetic variability and genome plasticity are key mechanisms. Trypanosomatid genomes are organized into large directional gene clusters that are transcribed polycistronically, where genes derived from the same polycistron may have very distinct mRNA levels. This particular mode of transcription implies that the control of gene expression operates mainly at post-transcriptional level. In this sense, gene duplications/losses were already associated with changes in mRNA levels in these parasites. Gene duplications also allow the generation of sequence variability, as the newly formed copy can diverge without loss of function of the original copy. Recently, aneuploidies have been shown to occur in several Leishmania species and T. cruzi strains. Although aneuploidies are usually associated with debilitating phenotypes in superior eukaryotes, recent data shows that it could also provide increased fitness in stress conditions and generate drug resistance in unicellular eukaryotes. In this review, we will focus on gene and chromosomal copy number variations and their relevance to the evolution of trypanosomatid parasites.

Chromosomal copy number variation reveals differential levels of genomic plasticity in distinct Trypanosoma cruzi strains.

BACKGROUND: Trypanosoma cruzi, the etiologic agent of Chagas disease, is currently divided into six discrete typing units (DTUs), named TcI-TcVI. CL Brener, the reference strain of the T. cruzi genome project, is a hybrid with a genome assembled into 41 putative chromosomes. Gene copy number variation (CNV) is well documented as an important mechanism to enhance gene expression and variability in T. cruzi. Chromosomal CNV (CCNV) is another level of gene CNV in which whole blocks of genes are expanded simultaneously. Although the T. cruzi karyotype is not well defined, several studies have demonstrated a significant variation in the size and content of chromosomes between different T. cruzi strains. Despite these studies, the extent of diversity in CCNV among T. cruzi strains based on a read depth coverage analysis has not been determined. RESULTS: We identify the CCNV in T. cruzi strains from the TcI, TcII and TcIII DTUs, by analyzing the depth coverage of short reads from these strains using the 41 CL Brener chromosomes as reference. This study led to the identification of a broader extent of CCNV in T. cruzi than was previously speculated. The TcI DTU strains have very few aneuploidies, while the strains from TcII and TcIII DTUs present a high degree of chromosomal expansions. Chromosome 31, which is the only chromosome that is supernumerary in all six T. cruzi samples evaluated in this study, is enriched with genes related to glycosylation pathways, highlighting the importance of glycosylation to parasite survival. CONCLUSIONS: Increased gene copy number due to chromosome amplification may contribute to alterations in gene expression, which represents a strategy that may be crucial for parasites that mainly depend on post-transcriptional mechanisms to control gene expression.

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.

Transcriptomic analysis of benznidazole-resistant and susceptible Trypanosoma cruzi populations.

BACKGROUND: Chagas disease (CD), caused by the parasite Trypanosoma cruzi, is a serious public health concern in Latin America. Nifurtimox and benznidazole (BZ), the only two drugs currently approved for the treatment of CD, have very low efficacies in the chronic phase of the disease and several toxic side effects. Trypanosoma cruzi strains that are naturally resistant to both drugs have been reported. We performed a comparative transcriptomic analysis of wild-type and BZ-resistant T. cruzi populations using high-throughput RNA sequencing to elucidate the metabolic pathways related to clinical drug resistance and identify promising molecular targets for the development of new drugs for treating CD. METHODS: All complementary DNA (cDNA) libraries were constructed from the epimastigote forms of each line, sequenced and analysed using the Prinseq and Trimmomatic tools for the quality analysis, STAR as the aligner for mapping the reads against the reference genome (T. cruzi Dm28c-2018), the Bioconductor package EdgeR for statistical analysis of differential expression and the Python-based library GOATools for the functional enrichment analysis. RESULTS: The analytical pipeline with an adjusted P-value of < 0.05 and fold-change > 1.5 identified 1819 transcripts that were differentially expressed (DE) between wild-type and BZ-resistant T. cruzi populations. Of these, 1522 (83.7%) presented functional annotations and 297 (16.2%) were assigned as hypothetical proteins. In total, 1067 transcripts were upregulated and 752 were downregulated in the BZ-resistant T. cruzi population. Functional enrichment analysis of the DE transcripts identified 10 and 111 functional categories enriched for the up- and downregulated transcripts, respectively. Through functional analysis we identified several biological processes potentially associated with the BZ-resistant phenotype: cellular amino acid metabolic processes, translation, proteolysis, protein phosphorylation, RNA modification, DNA repair, generation of precursor metabolites and energy, oxidation-reduction processes, protein folding, purine nucleotide metabolic processes and lipid biosynthetic processes. CONCLUSIONS: The transcriptomic profile of T. cruzi revealed a robust set of genes from different metabolic pathways associated with the BZ-resistant phenotype, proving that T. cruzi resistance mechanisms are multifactorial and complex. Biological processes associated with parasite drug resistance include antioxidant defenses and RNA processing. The identified transcripts, such as ascorbate peroxidase (APX) and iron superoxide dismutase (Fe-SOD), provide important information on the resistant phenotype. These DE transcripts can be further evaluated as molecular targets for new drugs against CD.

Complete Genome Sequences of Avian Metapneumovirus Subtype B Vaccine Strains from Brazil.

Avian metapneumovirus (aMPV) causes a highly contagious upper respiratory and reproductive disease in chickens, turkeys, and ducks. Here, complete genome sequences of aMPV-B vaccine strains BR/1890/E1/19 (PL21, Nemovac; Boehringer Ingelheim Animal Health, Brazil) and BR/1891/E2/19 (1062; Hipraviar, France) were sequenced and compared with the pathogenic field strain VCO3/60616.

Antigenic diversity of MASP gene family of Trypanosoma cruzi.

Trypanosoma cruzi, the etiological agent of Chagas disease (CD), is a heterogeneous species with high genetic and phenotypic diversity. MASP is the second largest multigene family of T. cruzi. The high degree of polymorphism of the family associated with its location at the surface of infective forms of T. cruzi suggests that MASP participates in mechanisms of host-parasite interaction. In this work, MASP members were divided into 7 subgroups based on protein sequence similarity, and one representative member from each subgroup was chosen to be expressed recombinantly. Immunogenicity of recombinant MASP proteins (rMASP) was investigated using different sera panels from T. cruzi infected mice. To mimic a natural condition in which different MASP members are expressed at the same time in the parasite population, a multiplex bead-based flow cytometry assay was also standardized. Results showed that rMASPs are poorly recognized by sera from mice infected with Colombiana strain, whereas sera from mice infected with CL Brener and Y display high reactivity against the majority of rMASPs tested. Flow cytometry showed that MASP recognition profile changes 10 days after infection. Also, multiplex assay suggests that MASP M1 and M2 are more immunogenic than the other MASP members evaluated that may play an immunodominant role during infection.

Trypanosoma cruzi Genomic Variability: Array Comparative Genomic Hybridization Analysis of Clone and Parental Strain.

Trypanosoma cruzi, the etiological agent of Chagas disease, exhibits extensive inter- and intrastrain genetic diversity. As we have previously described, there are some genetic differences between the parental G strain and its clone D11, which was isolated by the limiting dilution method and infection of cultured mammalian cells. Electrophoretic karyotyping and Southern blot hybridization of chromosomal bands with specific markers revealed chromosome length polymorphisms of small size with additional chromosomal bands in clone D11 and the maintenance of large syntenic groups. Both G strain and clone D11 belong to the T. cruzi lineage TcI. Here, we designed intraspecific array-based comparative genomic hybridization (aCGH) to identify chromosomal regions harboring copy-number variations between clone D11 and the G strain. DNA losses were more extensive than DNA gains in clone D11. Most alterations were flanked by repeated sequences from multigene families that could be involved in the duplication and deletion events. Several rearrangements were detected by chromoblot hybridization and confirmed by aCGH. We have integrated the information of genomic sequence data obtained by aCGH to the electrophoretic karyotype, allowing the reconstruction of possible recombination events that could have generated the karyotype of clone D11. These rearrangements may be explained by unequal crossing over between sister or homologous chromatids mediated by flanking repeated sequences and unequal homologous recombination via break-induced replication. The genomic changes detected by aCGH suggest the presence of a dynamic genome that responds to environmental stress by varying the number of gene copies and generating segmental aneuploidy.

Disruption of multiple copies of the Prostaglandin F2alpha synthase gene affects oxidative stress response and infectivity in Trypanosoma cruzi.

Chagas disease, caused by the protozoan Trypanosoma cruzi, is a serious chronic parasitic disease, currently treated with Nifurtimox (NFX) and Benznidazole (BZ). In addition to high toxicity, these drugs have low healing efficacy, especially in the chronic phase of the disease. The existence of drug-resistant T. cruzi strains and the occurrence of cross-resistance between BZ and NFX have also been described. In this context, it is urgent to study the metabolism of these drugs in T. cruzi, to better understand the mechanisms of resistance. Prostaglandin F2α synthase (PGFS) is an enzyme that has been correlated with parasite resistance to BZ, but the mechanism by which resistance occurs is still unclear. Our results show that the genome of the CL Brener clone of T. cruzi, contains five PGFS sequences and three potential pseudogenes. Using CRISPR/Cas9 we generated knockout cell lines in which all PGFS sequences were disrupted, as shown by PCR and western blotting analyses. The PGFS deletion did not alter the growth of the parasites or their susceptibility to BZ and NFX when compared to wild-type (WT) parasites. Interestingly, NTR-1 transcripts were shown to be upregulated in ΔPGFS mutants. Furthermore, the ΔPGFS parasites were 1.6 to 1.7-fold less tolerant to oxidative stress generated by menadione, presented lower levels of lipid bodies than the control parasites during the stationary phase, and were less infective than control parasites.

Parasite Genotype Is a Major Predictor of Mortality from Visceral Leishmaniasis.

Visceral leishmaniasis (VL) is a potentially fatal disease caused mainly by Leishmania infantum in South America and Leishmania donovani in Asia and Africa. Disease outcomes have been associated with patient genotype, nutrition, age, sex, comorbidities, and coinfections. In this study, we examine the effects of parasite genetic variation on VL disease severity in Brazil. We collected and sequenced the genomes of 109 L. infantum isolates from patients in northeastern Brazil and retrieved matching patient clinical data from medical records, including mortality, sex, HIV coinfection, and laboratory data (creatinine, hemoglobin, and leukocyte and platelet counts). We identified genetic differences between parasite isolates, including single nucleotide polymorphisms (SNPs), small insertions/deletions (indels), and variations in genic, intergenic, and chromosome copy numbers (copy number variants [CNVs]). To describe associations between the parasite genotypes and clinical outcomes, we applied quantitative genetics methods of heritability and genome-wide association studies (GWAS), treating clinical outcomes as traits that may be influenced by parasite genotype. Multiple aspects of the genetic analysis indicate that parasite genotype affects clinical outcomes. We estimate that parasite genotype explains 83% chance of mortality (narrow-sense heritability [h2] = 0.83 ± 0.17) and has a significant relationship with patient sex (h2 = 0.60 ± 0.27). Impacts of parasite genotype on other clinical traits are lower (h2 ≤ 0.34). GWAS analysis identified multiple parasite genetic loci that were significantly associated with clinical outcomes; 17 CNVs were significantly associated with mortality, two with creatinine, and one with bacterial coinfection, jaundice, and HIV coinfection, and two SNPs/indels and six CNVs were associated with age, jaundice, HIV and bacterial coinfections, creatinine, and/or bleeding sites. Parasite genotype is an important factor in VL disease severity in Brazil. Our analysis indicates that specific genetic differences between parasites act as virulence factors, enhancing risks of severe disease and mortality. More detailed understanding of these virulence factors could be exploited for novel therapies. IMPORTANCE Multiple factors contribute to the risk of mortality from visceral leishmaniasis (VL), including, patient genotype, comorbidities, and nutrition. Many of these factors are influenced by socioeconomic biases. Our work suggests that the virulence of the infecting parasite is an important risk factor for mortality. We pinpoint some specific genomic markers that are associated with mortality, which can lead to a greater understanding of the molecular mechanisms that cause severe VL disease, to the identification of genetic markers for virulent parasites, and to the development of drug and vaccine therapies.

Accessing the Variability of Multicopy Genes in Complex Genomes using Unassembled Next-Generation Sequencing Reads: The Case of Trypanosoma cruzi Multigene Families.

Repetitive elements cause assembly fragmentation in complex eukaryotic genomes, limiting the study of their variability. The genome of Trypanosoma cruzi, the parasite that causes Chagas disease, has a high repetitive content, including multigene families. Although many T. cruzi multigene families encode surface proteins that play pivotal roles in host-parasite interactions, their variability is currently underestimated, as their high repetitive content results in collapsed gene variants. To estimate sequence variability and copy number variation of multigene families, we developed a read-based approach that is independent of gene-specific read mapping and de novo assembly. This methodology was used to estimate the copy number and variability of MASP, TcMUC, and Trans-Sialidase (TS), the three largest T. cruzi multigene families, in 36 strains, including members of all six parasite discrete typing units (DTUs). We found that these three families present a specific pattern of variability and copy number among the distinct parasite DTUs. Inter-DTU hybrid strains presented a higher variability of these families, suggesting that maintaining a larger content of their members could be advantageous. In addition, in a chronic murine model and chronic Chagasic human patients, the immune response was focused on TS antigens, suggesting that targeting TS conserved sequences could be a potential avenue to improve diagnosis and vaccine design against Chagas disease. Finally, the proposed approach can be applied to study multicopy genes in any organism, opening new avenues to access sequence variability in complex genomes. IMPORTANCE Sequences that have several copies in a genome, such as multicopy-gene families, mobile elements, and microsatellites, are among the most challenging genomic segments to study. They are frequently underestimated in genome assemblies, hampering the correct assessment of these important players in genome evolution and adaptation. Here, we developed a new methodology to estimate variability and copy numbers of repetitive genomic regions and employed it to characterize the T. cruzi multigene families MASP, TcMUC, and transsialidase (TS), which are important virulence factors in this parasite. We showed that multigene families vary in sequence and content among the parasite's lineages, whereas hybrid strains have a higher sequence variability that could be advantageous to the parasite's survivability. By identifying conserved sequences within multigene families, we showed that the mammalian host immune response toward these multigene families is usually focused on the TS multigene family. These TS conserved and immunogenic peptides can be explored in future works as diagnostic targets or vaccine candidates for Chagas disease. Finally, this methodology can be easily applied to any organism of interest, which will aid in our understanding of complex genomic regions.

Disruption of Active Trans-Sialidase Genes Impairs Egress from Mammalian Host Cells and Generates Highly Attenuated Trypanosoma cruzi Parasites.

Trans-sialidases (TS) are unusual enzymes present on the surface of Trypanosoma cruzi, the causative agent of Chagas disease. Encoded by the largest gene family in the T. cruzi genome, only few members of the TS family have catalytic activity. Active trans-sialidases (aTS) are responsible for transferring sialic acid from host glycoconjugates to mucins, also present on the parasite surface. The existence of several copies of TS genes has impaired the use of reverse genetics to study this highly polymorphic gene family. Using CRISPR-Cas9, we generated aTS knockout cell lines displaying undetectable levels of TS activity, as shown by sialylation assays and labeling with antibodies that recognize sialic acid-containing mucins. In vitro infection assays showed that disruption of aTS genes does not affect the parasite's capacity to invade cells or to escape from the parasitophorous vacuole but resulted in impaired differentiation of amastigotes into trypomastigotes and parasite egress from the cell. When inoculated into mice, aTS mutants were unable to establish infection even in the highly susceptible gamma interferon (IFN-γ) knockout mice. Mice immunized with aTS mutants were fully protected against a challenge infection with the virulent T. cruzi Y strain. Altogether, our results confirmed the role of aTS as a T. cruzi virulence factor and indicated that aTS play a major role during the late stages of intracellular development and parasite egress. Notably, mutants lacking TS activity are completely avirulent in animal models of infection and may be used as a live attenuated vaccine against Chagas disease. IMPORTANCE Trypanosoma cruzi is the causative agent of Chagas disease, a neglected tropical disease that affects approximately 6 to 8 million people and for which there is no effective treatment or vaccine. The parasite expresses a family of surface proteins, named trans-sialidases, responsible for transferring sialic acid from host glycoconjugates to parasite mucins. Although recognized as a main virulence factor, the multiple roles of these proteins during infection have not yet been fully characterized, mainly because the presence of several copies of aTS genes has impaired their study using reverse genetics. By applying CRISPR-Cas9, we generated aTS knockout parasites and showed that, although aTS parasite mutants were able to infect cells in vitro, they have an impaired capacity to egress from the infected cell. Importantly, aTS mutants lost the ability to cause infection in vivo but provided full protection against a challenge infection with a virulent strain.

Experimental Selection of Paromomycin Resistance in Leishmania donovani Amastigotes Induces Variable Genomic Polymorphisms.

The relatively high post-treatment relapse rates of paromomycin (PMM) in visceral leishmaniasis treatment and the swift emergence of experimental drug resistance challenge its broad application and urge for rational use and monitoring of resistance. However, no causal molecular mechanisms to Leishmania PMM resistance have been identified so far. To gain insights into potential resistance mechanisms, twelve experimentally selected Leishmania donovani clonal lines and the non-cloned preselection population, with variable degrees of PMM resistance, were subjected to whole genome sequencing. To identify genomic variations potentially associated with resistance, SNPs, Indels, chromosomal somy and gene copy number variations were compared between the different parasite lines. A total of 11 short nucleotide variations and the copy number alterations in 39 genes were correlated to PMM resistance. Some of the identified genes are involved in transcription, translation and protein turn-over (transcription elongation factor-like protein, RNA-binding protein, ribosomal protein L1a, 60S ribosomal protein L6, eukaryotic translation initiation factor 4E-1, proteasome regulatory non-ATP-ase subunit 3), virulence (major surface protease gp63, protein-tyrosine phosphatase 1-like protein), mitochondrial function (ADP/ATP mitochondrial carrier-like protein), signaling (phosphatidylinositol 3-related kinase, protein kinase putative and protein-tyrosine phosphatase 1-like protein) and vesicular trafficking (ras-related protein RAB1). These results indicate that, in Leishmania, the aminoglycoside PMM affects protein translational processes and underlines the complex and probably multifactorial origin of resistance.

Targeted Deletion of Centrin in Leishmania braziliensis Using CRISPR-Cas9-Based Editing.

Leishmania braziliensis is the main causative agent of Tegumentary Leishmaniasis in the Americas. However, difficulties related to genome manipulation, experimental infection, and parasite growth have so far limited studies with this species. CRISPR-Cas9-based technology has made genome editing more accessible, and here we have successfully employed the LeishGEdit approach to attenuate L. braziliensis. We generated a transgenic cell line expressing Cas9 and T7 RNA polymerase, which was employed for the targeted deletion of centrin, a calcium-binding cytoskeletal protein involved in the centrosome duplication in eukaryotes. Centrin-deficient Leishmania exhibit growth arrest at the amastigote stage. Whole-genome sequencing of centrin-deficient L. braziliensis (LbCen-/- ) did not indicate the presence of off-target mutations. In vitro, the growth rates of LbCen-/- and wild-type promastigotes were similar, but axenic and intracellular LbCen-/- amastigotes showed a multinucleated phenotype with impaired survival following macrophage infection. Upon inoculation into BALB/c mice, LbCen-/- were detected at an early time point but failed to induce lesion formation, contrary to control animals, infected with wild-type L. braziliensis. A significantly lower parasite burden was also observed in mice inoculated with LbCen-/- , differently from control mice. Given that centrin-deficient Leishmania sp. have become candidates for vaccine development, we propose that LbCen-/- can be further explored for the purposes of immunoprophylaxis against American Tegumentary Leishmaniasis.

The gene repertoire of the main cysteine protease of Trypanosoma cruzi, cruzipain, reveals four sub-types with distinct active sites.

Cruzipains are the main papain-like cysteine proteases of Trypanosoma cruzi, the protozoan parasite that causes Chagas disease. Encoded by a multigenic family, previous studies have estimated the presence of dozens of copies spread over multiple chromosomes in different parasite strains. Here, we describe the complete gene repertoire of cruzipain in three parasite strains, their genomic organization, and expression pattern throughout the parasite life cycle. Furthermore, we have analyzed primary sequence variations among distinct family members as well as structural differences between the main groups of cruzipains. Based on phylogenetic inferences and residue positions crucial for enzyme function and specificity, we propose the classification of cruzipains into two families (I and II), whose genes are distributed in two or three separate clusters in the parasite genome, according with the strain. Family I comprises nearly identical copies to the previously characterized cruzipain 1/cruzain, whereas Family II encompasses three structurally distinct sub-types, named cruzipain 2, cruzipain 3, and cruzipain 4. RNA-seq data derived from the CL Brener strain indicates that Family I genes are mainly expressed by epimastigotes, whereas trypomastigotes mainly express Family II genes. Significant differences in the active sites among the enzyme sub-types were also identified, which may play a role in their substrate selectivity and impact their inhibition by small molecules.