Search

Abstract 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.

ABSTRACT 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.

Visceral leishmaniasis is an important global health problem with an estimated of 50,000 to 90,000 new cases per year. VL is the most serious form of leishmaniasis as it can be fatal in 95% of the cases if it remains untreated. VL is a particularly acute problem in Brazil which contributed with 97% of all cases reported in 2020 in the Americas. In this country, VL affects mainly the poorest people in both urban and rural areas and continues to have a high mortality rate estimated around 8.15%. Here, we performed a temporal parasite population study using whole genome sequence data from a set of 34 canine isolates sampled in 2008, 2012 and 2015 from a re-emergent focus in Southeastern Brazil. Our study found the presence of two distinct sexual subpopulations that corresponded to two isolation periods. These subpopulations diverged hundreds of years ago with no apparent gene flow between them suggesting a process of rapid replacement during a two-year period. Sequence comparisons and analysis of nucleotide diversity also showed evidence of balancing selection acting on transport-related genes and antigenic families. To our knowledge this is the first population genomic study showing a turn-over of parasite populations in an endemic region for leishmaniasis. The complexity and rapid adaptability of these parasites pose new challenges to control activities and demand more integrated approaches to understand this disease in New World foci.

ABSTRACT 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.

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.

Trypanosoma cruzi, a zoonotic kinetoplastid protozoan parasite, is the causative agent of American trypanosomiasis (Chagas disease). Having a very plastic, repetitive and complex genome, the parasite displays a highly diverse repertoire of surface molecules, with pivotal roles in cell invasion, immune evasion and pathogenesis. Before 2016, the complexity of the genomic regions containing these genes impaired the assembly of a genome at chromosomal level, making it impossible to study the structure and function of the several thousand repetitive genes encoding the surface molecules of the parasite. We here describe the genome assembly of the Sylvio X10/1 genome sequence, which since 2016 has been used as a reference genome sequence for T. cruzi clade I (TcI), produced using high coverage PacBio single-molecule sequencing. It was used to analyze deep Illumina sequence data from 34 T. cruzi TcI isolates and clones from different geographic locations, sample sources and clinical outcomes. Resolution of the surface molecule gene distribution showed the unusual duality in the organization of the parasite genome, a synteny of the core genomic region with related protozoa flanked by unique and highly plastic multigene family clusters encoding surface antigens. The presence of abundant interspersed retrotransposons in these multigene family clusters suggests that these elements are involved in a recombination mechanism for the generation of antigenic variation and evasion of the host immune response on these TcI strains. The comparative genomic analysis of the cohort of TcI strains revealed multiple cases of such recombination events involving surface molecule genes and has provided new insights into T. cruzi population structure.

The public health system is extremely dependent on the use of vaccines to immunize the population from a series of infectious and dangerous diseases, preventing the system from collapsing and millions of people dying every year. However, to develop these vaccines and effectively monitor these diseases, it is necessary to use accurate diagnostic methods capable of identifying highly immunogenic regions within a given pathogenic protein. Existing experimental methods are expensive, time-consuming, and require arduous laboratory work, as they require the screening of a large number of potential candidate epitopes, making the methods extremely laborious, especially for application to larger microorganisms. In the last decades, researchers have developed in silico prediction methods, based on machine learning, to identify these markers, to drastically reduce the list of potential candidate epitopes for experimental tests, and, consequently, to reduce the laborious task associated with their mapping. Despite these efforts, the tools and methods still have low accuracy, slow diagnosis, and offline training. Thus, we develop a method to predict B-cell linear epitopes which are based on a Fuzzy-ARTMAP neural network architecture, called BepFAMN (B Epitope Prediction Fuzzy ARTMAP Artificial Neural Network). This was trained using a linear averaging scheme on 15 properties that include an amino acid ratio scale and a set of 14 physicochemical scales. The database used was obtained from the IEDB website, from which the amino acid sequences with the annotations of their positive and negative epitopes were taken. To train and validate the knowledge models, five-fold cross-validation and competition techniques were used. The BepiPred-2.0 database, an independent database, was used for the tests. In our experiment, the validation dataset reached sensitivity = 91.50%, specificity = 91.49%, accuracy = 91.49%, MCC = 0.83, and an area under the curve (AUC) ROC of approximately 0.9289. The result in the testing dataset achieves a significant improvement, with sensitivity = 81.87%, specificity = 74.75%, accuracy = 78.27%, MCC = 0.56, and AOC = 0.7831. These achieved values demonstrate that BepFAMN outperforms all other linear B-cell epitope prediction tools currently used. In addition, the architecture provides mechanisms for online training, which allow the user to find a new B-cell linear epitope, and to improve the model without need to re-train itself with the whole dataset. This fact contributes to a considerable reduction in the number of potential linear epitopes to be experimentally validated, reducing laboratory time and accelerating the development of diagnostic tests, vaccines, and immunotherapeutic approaches.

Abstract 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.

Dengue is currently one of the most important arbovirus infections worldwide. Early diagnosis is important for disease outcome, particularly for those afflicted with the severe forms of infection. The goal of this work was to identify conserved and polymorphic linear B-cell Dengue virus (DENV) epitopes that could be used for diagnostic purposes. To this end, we aligned the predicted viral proteome of the four DENV serotype and performed in silico B-cell epitope mapping. We developed a script in Perl integrating alignment and prediction information to identify potential serotype-specific epitopes. We excluded epitopes that were similarly present in the yellow fever and zika viruses’ proteomes. A total of 15 polymorphic and nine conserved peptides among DENV serotypes were selected. Peptides were spotted on cellulose membranes and tested against sera from rabbits that were monoinfected with each DENV serotype. Although serotype-specific peptides failed to recognize any sera, three conserved peptides were recognized by all anti-dengue sera and were included on an ELISA test employing a well-characterized human sera bank. Of the three peptides, one was able to efficiently identify sera from all four DENV serotypes and to discriminate them from Zika virus positive sera.

Trypanosoma cruzi is the etiologic agent of Chagas disease, a neglected tropical disease that affects millions of people mainly in Latin America. To establish a life-long infection, T. cruzi must subvert the vertebrate host’s immune system, using strategies that can be traced to the parasite’s life cycle. Once inside the vertebrate host, metacyclic trypomastigotes rapidly invade a wide variety of nucleated host cells in a membrane-bound compartment known as the parasitophorous vacuole, which fuses to lysosomes, originating the phagolysosome. In this compartment, the parasite relies on a complex network of antioxidant enzymes to shield itself from lysosomal oxygen and nitrogen reactive species. Lysosomal acidification of the parasitophorous vacuole is an important factor that allows trypomastigote escape from the extremely oxidative environment of the phagolysosome to the cytoplasm, where it differentiates into amastigote forms. In the cytosol of infected macrophages, oxidative stress instead of being detrimental to the parasite, favors amastigote burden, which then differentiates into bloodstream trypomastigotes. Trypomastigotes released in the bloodstream upon the rupture of the host cell membrane express surface molecules, such as calreticulin and GP160 proteins, which disrupt initial and key components of the complement pathway, while others such as GPI-mucins stimulate immunoregulatory receptors, delaying the progression of a protective immune response. After an immunologically silent entry at the early phase of infection, T. cruzi elicits polyclonal B cell activation, hypergammaglobulinemia, and unspecific anti-T. cruzi antibodies, which are inefficient in controlling the infection. Additionally, the co-expression of several related but not identical epitopes derived from trypomastigote surface proteins delays the generation of T. cruzi-specific neutralizing antibodies. Later in the infection, the establishment of an anti-T. cruzi CD8+ immune response focused on the parasite’s immunodominant epitopes controls parasitemia and tissue infection, but fails to completely eliminate the parasite. This outcome is not detrimental to the parasite, as it reduces host mortality and maintains the parasite infectivity toward the insect vectors.

To improve transfection efficiency in Trypanosoma cruzi, we developed a new electroporation protocol and expression vectors which use luciferase and green and red fluorescent proteins as reporter genes. In transient transfections, the electroporation conditions reported here resulted in luciferase expression 100 times higher than the levels obtained with previously described protocols. To verify whether sequences containing different trans-splicing signals influence reporter gene expression, we compared DNA fragments corresponding to 5' untranslated plus intergenic (5' UTR plus Ig) regions from GAPDH, TcP2ß, a- and ß- tubulin and amastin genes. Vectors containing sequences derived from the first four genes presented similar efficiencies and resulted in luciferase expression in transiently transfected epimastigotes that was up to 10 times higher than that for a control vector. In contrast, the amastin 5' UTR plus Ig resulted in lower levels of reporter gene expression. We also constructed a vector containing an expression cassette designed to be targeted to the tubulin locus of the parasite. [ABSTRACT FROM AUTHOR]

The availability of Trypanosomatid genomic data in public databases has opened myriad experimental possibilities that have contributed to a more comprehensive understanding of the biology of these parasites and their interactions with hosts. In this review, after brief remarks on the history of the Trypanosoma cruzi and Leishmania genome initiatives, we present an overview of the relevant contributions of genomics, transcriptomics and functional genomics, discussing the primary obstacles, challenges, relevant achievements and future perspectives of these technologies.

BackgroundProphages play a significant role in prokaryotic evolution, often altering the function of the cell that they infect via transfer of new genes e.g., virulence or antibiotic resistance factors, inactivation of existing genes or by modifying gene expression. Recently, phage therapy has gathered renewed interest as a promising alternative to control bacterial infections. Cataloging the repertoire of prophages in large collections of species' genomes is an important initial step in understanding their evolution and potential therapeutic utility. However, current widely-used tools for identifying prophages within bacterial genome sequences are mainly web-based, can have long response times, and do not scale to keep pace with the many thousands of genomes currently being sequenced routinely.MethodologyIn this work, we present ProphET, an easy to install prophage predictor to be used in Linux operation system, without the constraints associated with a web-based tool. ProphET predictions rely on similarity searches against a database of prophage genes, taking as input a bacterial genome sequence in FASTA format and its corresponding gene annotation in GFF. ProphET identifies prophages in three steps: similarity search, calculation of the density of prophage genes, and edge refinement. ProphET performance was evaluated and compared with other phage predictors based on a set of 54 bacterial genomes containing 267 manually annotated prophages.Findings and conclusionsProphET identifies prophages in bacterial genomes with high precision and offers a fast, highly scalable alternative to widely-used web-based applications for prophage detection.

Trypanosoma cruzi, the protozoan that causes Chagas disease, has a complex life cycle involving several morphologically and biochemically distinct stages that establish intricate interactions with various insect and mammalian hosts. It has also a heterogeneous population structure comprising strains with distinct properties such as virulence, sensitivity to drugs, antigenic profile and tissue tropism. We present a comparative transcriptome analysis of two cloned T. cruzi strains that display contrasting virulence phenotypes in animal models of infection: CL Brener is a virulent clone and CL-14 is a clone that is neither infective nor pathogenic in in vivo models of infection. Gene expression analysis of trypomastigotes and intracellular amastigotes harvested at 60 and 96 hours post-infection (hpi) of human fibroblasts revealed large differences that reflect the parasite's adaptation to distinct environments during the infection of mammalian cells, including changes in energy sources, oxidative stress responses, cell cycle control and cell surface components. While extensive transcriptome remodeling was observed when trypomastigotes of both strains were compared to 60 hpi amastigotes, differences in gene expression were much less pronounced when 96 hpi amastigotes and trypomastigotes of CL Brener were compared. In contrast, the differentiation of the avirulent CL-14 from 96 hpi amastigotes to extracellular trypomastigotes was associated with considerable changes in gene expression, particularly in gene families encoding surface proteins such as trans-sialidases, mucins and the mucin associated surface proteins (MASPs). Thus, our comparative transcriptome analysis indicates that the avirulent phenotype of CL-14 may be due, at least in part, to a reduced or delayed expression of genes encoding surface proteins that are associated with the transition of amastigotes to trypomastigotes, an essential step in the establishment of the infection in the mammalian host. Confirming the role of members of the trans-sialidase family of surface proteins for parasite differentiation, transfected CL-14 constitutively expressing a trans-sialidase gene displayed faster kinetics of trypomastigote release in the supernatant of infected cells compared to wild type CL-14.

The protozoan parasite Trypanosoma cruzi is the etiologic agent of Chagas disease, a highly debilitating human pathology that affects millions of people in the Americas. The sequencing of this parasite's genome reveals that trans-sialidase/trans-sialidase-like (TcS), a polymorphic protein family known to be involved in several aspects of T. cruzi biology, is the largest T. cruzi gene family, encoding more than 1,400 genes. Despite the fact that four TcS groups are well characterized and only one of the groups contains active trans-sialidases, all members of the family are annotated in the T. cruzi genome database as trans-sialidase. After performing sequence clustering analysis with all TcS complete genes, we identified four additional groups, demonstrating that the TcS family is even more heterogeneous than previously thought. Interestingly, members of distinct TcS groups show distinctive patterns of chromosome localization. Members of the TcSgroupII, which harbor proteins involved in host cell attachment/invasion, are preferentially located in subtelomeric regions, whereas members of the largest and new TcSgroupV have internal chromosomal locations. Real-time RT-PCR confirms the expression of genes derived from new groups and shows that the pattern of expression is not similar within and between groups. We also performed B-cell epitope prediction on the family and constructed a TcS specific peptide array, which was screened with sera from T. cruzi-infected mice. We demonstrated that all seven groups represented in the array are antigenic. A highly reactive peptide occurs in sixty TcS proteins including members of two new groups and may contribute to the known cross-reactivity of T. cruzi epitopes during infection. Taken together, our results contribute to a better understanding of the real complexity of the TcS family and open new avenues for investigating novel roles of this family during T. cruzi infection.