Search

Predicting individual vaccine responses remains a significant challenge due to the complexity and variability of immune processes. To address this gap, we developed immunaut , an open-source, data-driven framework implemented as an R package specifically designed for all systems vaccinologists seeking to analyze and predict immunological outcomes across diverse vaccination settings. Leveraging one of the most comprehensive live attenuated influenza vaccine (LAIV) datasets to date - 244 Gambian children enrolled in a phase 4 immunogenicity study - immunaut integrates humoral, mucosal, cellular, transcriptomic, and microbiological parameters collected before and after vaccination, providing an unprecedentedly holistic view of LAIV-induced immunity. Through advanced dimensionality reduction, clustering, and predictive modeling, immunaut identifies distinct immunophenotypic responder profiles and their underlying baseline determinants. In this study, immunaut delineated three immunophenotypes: (1) CD8 T-cell responders, marked by strong baseline mucosal immunity and extensive prior influenza virus exposure that boosts memory CD8 T-cell responses, without generating influenza virus-specific antibody responses; (2) Mucosal responders, characterized by pre-existing systemic influenza A virus immunity (specifically to H3N2) and stable epithelial integrity, leading to potent mucosal IgA expansions and subsequent seroconversion to influenza B virus; and (3) Systemic, broad influenza A virus responders, who start with relatively naive immunity and leverage greater initial viral replication to drive broad systemic antibody responses against multiple influenza A virus variants beyond those included in the LAIV vaccine. By integrating pathway-level analysis, model-derived contribution scores, and hierarchical decision rules, immunaut elucidates how distinct immunological landscapes shape each response trajectory and how key baseline features, including pre-existing immunity, mucosal preparedness, and cellular support, dictate vaccine outcomes. Collectively, these findings emphasize the power of integrative, predictive frameworks to advance precision vaccinology, and highlight immunaut as a versatile, community-available resource for optimizing immunization strategies across diverse populations and vaccine platforms., One-Sentence Summary: Using one of the most comprehensive LAIV datasets compiled to date, immunaut , an integrative machine learning framework, identifies distinct immunophenotypic responder groups shaped by baseline immune landscapes, advancing precision vaccinology and guiding more effective, personalized immunization strategies., Graphical Abstract: Immunaut , an automated framework for mapping and predicting vaccine response immunotypes. Step 1 outlines the identification of vaccine response outcomes using pre- and post-vaccination data integration across immune features, including antibodies, flu-specific T-cells, and immunophenotyping at mucosal and systemic sites. Clustering methods define the vaccine response landscape, stability, and validation through t-SNE-based visualization. Step 2 leverages an automated machine learning modeling approach, to enhance the accuracy and interpretability of vaccine response predictions, enabling stratification and targeted intervention strategies for personalized vaccine immunogenicity.

Chikungunya fever (CHIKF), caused by the Chikungunya virus (CHIKV), manifests as acute febrile illness often associated with polyarthritis and polyarthralgia. Although the acute symptoms resolve within two weeks, many patients experience prolonged joint pain and inflammation, resembling rheumatoid arthritis (RA). This study aimed to identify molecular markers related to joint pain and chronicity in CHIKV-infected individuals by analyzing blood transcriptomes using bulk RNA sequencing. B- and T-cell receptor (BCR and TCR) diversity was assessed through computational analysis of RNA-seq data, revealing a significant reduction in CDR3 diversity in CHIKV-infected individuals compared to healthy controls. This reduced diversity was associated with the upregulation of genes involved in osteoclast differentiation and activation, particularly through the RANK/RANKL signaling pathway. These findings suggest a potential link between immune dysregulation and enhanced osteoclast activity, which may contribute to the persistence of joint pain in chronic CHIKF. Targeting osteoclast-related pathways could offer therapeutic strategies for managing chronic symptoms in CHIKF patients.

DNA is an intelligent data storage medium due to its stability and high density. It has been used by nature for over 3.5 billion years. Compared with traditional methods, DNA offers better compression and physical density. DNA can retain information for thousands of years. However, challenges exist in scalability, standardization, metadata gathering, biocybersecurity, and specialized tools. Addressing these challenges is crucial for widespread implementation. Collaboration among experts, as well as keeping the future in mind, is needed to unlock the full potential of DNA data storage, which promises low energy costs, high-density storage, and long-term stability., Competing Interests: Declaration of interests The authors have no conflicts of interest to declare., (Copyright © 2023 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Introduction: The rVSVDG-ZEBOV-GP (Ervebo®) vaccine is both immunogenic and protective against Ebola. However, the vaccine can cause a broad range of transient adverse reactions, from headache to arthritis. Identifying baseline reactogenicity signatures can advance personalized vaccinology and increase our understanding of the molecular factors associated with such adverse events., Methods: In this study, we developed a machine learning approach to integrate prevaccination gene expression data with adverse events that occurred within 14 days post-vaccination., Results and Discussion: We analyzed the expression of 144 genes across 343 blood samples collected from participants of 4 phase I clinical trial cohorts: Switzerland, USA, Gabon, and Kenya. Our machine learning approach revealed 22 key genes associated with adverse events such as local reactions, fatigue, headache, myalgia, fever, chills, arthralgia, nausea, and arthritis, providing insights into potential biological mechanisms linked to vaccine reactogenicity., Competing Interests: Author IM was employed by the company Microbiotec Srl. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2023 Gonzalez Dias Carvalho, Dominguez Crespo Hirata, Mano Alves, Moscardini, do Nascimento, Costa-Martins, Sorgi, Harandi, Ferreira, Vianello, Haks, Ottenhoff, Santoro, Martinez-Murillo, Huttner, Siegrist, Medaglini and Nakaya.)

Background: Public health research frequently requires the integration of information from different data sources. However, errors in the records and the high computational costs involved make linking large administrative databases using record linkage (RL) methodologies a major challenge., Methods: We present Tucuxi-BLAST, a versatile tool for probabilistic RL that utilizes a DNA-encoded approach to encrypt, analyze and link massive administrative databases. Tucuxi-BLAST encodes the identification records into DNA. BLASTn algorithm is then used to align the sequences between databases. We tested and benchmarked on a simulated database containing records for 300 million individuals and also on four large administrative databases containing real data on Brazilian patients., Results: Our method was able to overcome misspellings and typographical errors in administrative databases. In processing the RL of the largest simulated dataset (200k records), the state-of-the-art method took 5 days and 7 h to perform the RL, while Tucuxi-BLAST only took 23 h. When compared with five existing RL tools applied to a gold-standard dataset from real health-related databases, Tucuxi-BLAST had the highest accuracy and speed. By repurposing genomic tools, Tucuxi-BLAST can improve data-driven medical research and provide a fast and accurate way to link individual information across several administrative databases., Competing Interests: Helder I. Nakaya and Robson Souza are Academic Editors for PeerJ., (© 2022 Araujo et al.)

Globally, malaria and human immunodeficiency virus (HIV) are both independently associated with a massive burden of disease and death. While their co-infection has been well studied for Plasmodium falciparum, scarce data exist regarding the association of P. vivax and HIV. In this cohort study, we assessed the effect of HIV on the risk of vivax malaria infection and recurrence during a 4-year follow-up period in an endemic area of the Brazilian Amazon. For the purpose of this study, we obtained clinical information from January 2012 to December 2016 from two databases. HIV screening data were acquired from the clinical information system at the tropical hospital Fundação de Medicina Tropical Dr. Heitor Vieira Dourado (FMT-HVD). The National Malaria Surveillance database (SIVEP malaria) was utilized to identify malaria infections during a 4-year follow-up period after diagnosis of HIV. Both datasets were combined via data linkage. Between 2012 and 2016, a total of 42,121 people were screened for HIV, with 1569 testing positive (3.7%). Out of all the patients diagnosed with HIV, 198 had at least one episode of P. vivax malaria in the follow-up. In the HIV-negative group, 711 participants had at least one P. vivax malaria episode. When comparing both groups, HIV patients had a 6.48 [(5.37-7.83); P < 0.0001] (adjusted relative risk) greater chance of acquiring P. vivax malaria. Moreover, being of the male gender [ARR = 1.41 (1.17-1.71); P < 0.0001], Amerindian ethnicity [ARR = 2.77 (1.46-5.28); P < 0.0001], and a resident in a municipality of the Metropolitan region of Manaus [ARR = 1.48 (1.02-2.15); P = 0.038] were independent risk factors associated with an increased risk of clinical malaria. Education ≥ 8 years [ARR = 0.41 (0.26-0.64); P < 0.0001] and living in the urban area [ARR = 0.44 (0.24-0.80); P = 0.007] were associated to a lower risk of P. vivax malaria. A total of 28 (14.1%) and 180 (25.3%) recurrences (at least a second clinical malaria episode) were reported in the HIV-positive and HIV-negative groups, respectively. After adjusting for sex and education, HIV-positive status was associated with a tendency towards protection from P. vivax malaria recurrences [ARR = 0.55 (0.27-1.10); P = 0.090]. HIV status was not associated with hospitalizations due to P. vivax malaria. CD4 + counts and viral load were not associated with recurrences of P. vivax malaria. No significant differences were found in the distribution of parasitemia between HIV-negative and HIV-positive P. vivax malaria patients. Our results suggest that HIV-positive status is a risk factor for vivax malaria infection, which represents an additional challenge that should be addressed during elimination efforts., (© 2022. The Author(s).)

Breastfeeding protects against mucosal infections in infants. The underlying mechanisms through which immunity develops in human milk following maternal infection with mucosal pathogens are not well understood. We simulated nasal mucosal influenza infection through live attenuated influenza vaccination (LAIV) and compared immune responses in milk to inactivated influenza vaccination (IIV). Transcriptomic analysis was performed on RNA extracted from human milk cells to evaluate differentially expressed genes and pathways on days 1 and 7 post-vaccination. Both LAIV and IIV vaccines induced influenza-specific IgA that persisted for at least 6 months. Regulation of type I interferon production, toll-like receptor, and pattern recognition receptor signaling pathways were highly upregulated in milk on day 1 following LAIV but not IIV at any time point. Upregulation of innate immunity in human milk may provide timely protection against mucosal infections until antigen-specific immunity develops in the human milk-fed infant., (© 2022. The Author(s), under exclusive licence to Society for Mucosal Immunology.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Cortez, Lima, Reis-Cunha, Bartholomeu, Villacis, Rogatto, Costa-Martins, Marchiano, do Carmo, da Silveira and Marini.)

In children lacking influenza-specific adaptive immunity, upper respiratory tract innate immune responses may influence viral replication and disease outcome. We use trivalent live attenuated influenza vaccine (LAIV) as a surrogate challenge model in children aged 24-59 months to identify pre-infection mucosal transcriptomic signatures associated with subsequent viral shedding. Upregulation of interferon signaling pathways prior to LAIV is significantly associated with lower strain-specific viral loads (VLs) at days 2 and 7. Several interferon-stimulated genes are differentially expressed in children with pre-LAIV asymptomatic respiratory viral infections and negatively correlated with LAIV VLs. Upregulation of genes enriched in macrophages, neutrophils, and eosinophils is associated with lower VLs and found more commonly in children with asymptomatic viral infections. Variability in pre-infection mucosal interferon gene expression in children may impact the course of subsequent influenza infections. This variability may be due to frequent respiratory viral infections, demonstrating the potential importance of mucosal virus-virus interactions in children., Competing Interests: J.O.M. reports compensation for consulting services with Cellarity. No other conflicts of interest are reported by the authors., (© 2021 The Author(s).)

This work reports on the whole-genome sequencing of Leishmania infantum chagasi from Honduras (Central America) and Brazil (South America).

Background: Mycobacterium tuberculosis, the etiological agent of tuberculosis, has at least four ATP-Binding Cassette (ABC) transporters dedicated to carbohydrate uptake: LpqY/SugABC, UspABC, Rv2038c-41c, and UgpAEBC. LpqY/SugABC transporter is essential for M. tuberculosis survival in vivo and potentially involved in the recycling of cell wall components. The three-dimensional structures of substrate-binding proteins (SBPs) LpqY, UspC, and UgpB were described, however, questions about how these proteins interact with the cognate transporter are still being explored. Components of these transporters, such as SBPs, show high immunogenicity and could be used for the development of diagnostic and therapeutic tools. In this work, we used a phylogenetic and structural bioinformatics approach to compare the four systems, in an attempt to predict functionally important regions., Results: Through the analysis of the putative orthologs of the carbohydrate ABC importers in species of Mycobacterium genus it was shown that Rv2038c-41c and UgpAEBC systems are restricted to pathogenic species. We showed that the components of the four ABC importers are phylogenetically separated into four groups defined by structural differences in regions that modulate the functional activity or the interaction with domain partners. The regulatory region in nucleotide-binding domains, the periplasmic interface in transmembrane domains and the ligand-binding pocket of the substrate-binding proteins define their substrates and segregation in different branches. The interface between transmembrane domains and nucleotide-binding domains show conservation of residues and charge., Conclusions: The presence of four ABC transporters in M. tuberculosis dedicated to uptake and transport of different carbohydrate sources, and the exclusivity of at least two of them being present only in pathogenic species of Mycobacterium genus, highlights their relevance in virulence and pathogenesis. The significant differences in the SBPs, not present in eukaryotes, and in the regulatory region of NBDs can be explored for the development of inhibitory drugs targeting the bacillus. The possible promiscuity of NBDs also contributes to a less specific and more comprehensive control approach., (© 2021. The Author(s).)

Background & Aims: Obesity is associated with low grade systemic inflammation and insulin resistance. Although metabolic and immunological changes may contribute to the increased risk for COVID-19 mortality in obese, little is known about the impact of obesity in the lungs of patients with COVID-19., Methods: We analyzed gene expression profiles of autopsy lungs of a cohort of 14 COVID-19 patients and 4 control individuals. Patients were divided into 3 groups according to their comorbidities: hypertension, type 2 diabetes (T2D) and obesity. We then identified the molecular alterations associated with these comorbidities in COVID-19 patients., Results: Patients with only hypertension showed higher levels of inflammatory genes and B-cell related genes when compared to those with T2D and obesity. However, the levels of IFN-gamma, IL22, and CD274 (a ligand that binds to receptor PD1) were higher in COVID-19 patients with T2D and obesity. Several metabolic- and immune-associated genes such as G6PD, LCK and IL10 were significantly induced in the lungs of the obese group., Conclusion: Our findings suggest that SARS-CoV-2 infection in the lungs may exacerbate the immune response and chronic condition in obese COVID-19 patients., Competing Interests: Declaration of competing interest None of the authors have conflicts of interest to disclose., (Copyright © 2021 European Society for Clinical Nutrition and Metabolism. Published by Elsevier Ltd. All rights reserved.)

Background: The progression and severity of COVID-19 vary significantly in the population. While the hallmarks of SARS-CoV-2 and severe COVID-19 within routine laboratory parameters are emerging, the impact of sex and age on these profiles is still unknown., Methods: A multidimensional analysis was performed involving millions of records of laboratory parameters and diagnostic tests for 178 887 individuals from Brazil, of whom 33 266 tested positive for SARS-CoV-2. Analyzed data included those relating to complete blood cell count, electrolytes, metabolites, arterial blood gases, enzymes, hormones, cancer biomarkers, and others., Findings: COVID-19 induced similar alterations in laboratory parameters in males and females. CRP and ferritin were increased, especially in older men with COVID-19, whereas abnormal liver function tests were common across several age groups, except for young women. Low peripheral blood basophils and eosinophils were more common in the elderly with COVID-19. Both male and female COVID-19 patients admitted to intensive care units displayed alterations in the coagulation system, and higher values for neutrophils, CRP, and lactate dehydrogenase., Conclusions: Our study uncovered the laboratory profiles of a large cohort of COVID-19 patients, which formed the basis of discrepancies influenced by aging and biological sex. These profiles directly linked COVID-19 disease presentation to an intricate interplay between sex, age, and immune activation., (Copyright © 2021 The Authors. Published by Elsevier Ltd.. All rights reserved.)

The etiological agent of Chagas disease, Trypanosoma cruzi, is a complex of seven genetic subdivisions termed discrete typing units (DTUs), TcI-TcVI and Tcbat. The relevance of T. cruzi genetic diversity to the variable clinical course of the disease, virulence, pathogenicity, drug resistance, transmission cycles and ecological distribution requires understanding the parasite origin and population structure. In this study, we introduce the PhyloQuant approach to infer the evolutionary relationships between organisms based on differential mass spectrometry-based quantitative features. In particular, large scale quantitative bottom-up proteomics features (MS1, iBAQ and LFQ) were analyzed using maximum parsimony, showing a correlation between T. cruzi DTUs and closely related trypanosomes' protein expression and sequence-based clustering. Character mapping enabled the identification of synapomorphies, herein the proteins and their respective expression profiles that differentiate T. cruzi DTUs and trypanosome species. The distance matrices based on phylogenetics and PhyloQuant clustering showed statistically significant correlation highlighting the complementarity between the two strategies. Moreover, PhyloQuant allows the identification of differentially regulated and strain/DTU/species-specific proteins, and has potential application in the identification of specific biomarkers and candidate therapeutic targets.

Retrotransposon Hot Spot (RHS) is the most abundant gene family in Trypanosoma cruzi, with unknown function in this parasite. The aim of this work was to shed light on the organization and expression of RHS in T. cruzi. The diversity of the RHS protein family in T. cruzi was demonstrated by phylogenetic and recombination analyses. Transcribed sequences carrying the RHS domain were classified into ten distinct groups of monophyletic origin. We identified numerous recombination events among the RHS and traced the origins of the donors and target sequences. The transcribed RHS genes have a mosaic structure that may contain fragments of different RHS inserted in the target sequence. About 30% of RHS sequences are located in the subtelomere, a region very susceptible to recombination. The evolution of the RHS family has been marked by many events, including gene duplication by unequal mitotic crossing-over, homologous, as well as ectopic recombination, and gene conversion. The expression of RHS was analyzed by immunofluorescence and immunoblotting using anti-RHS antibodies. RHS proteins are evenly distributed in the nuclear region of T. cruzi replicative forms (amastigote and epimastigote), suggesting that they could be involved in the control of the chromatin structure and gene expression, as has been proposed for T. brucei .

The structure and function of the immune system is governed by complex networks of interactions between cells and molecular components. Vaccination perturbs these networks, triggering specific pathways to induce cellular and humoral immunity. Systems vaccinology studies have generated vast data sets describing the genes related to vaccination, motivating the use of new approaches to identify patterns within the data. Here, we describe a framework called Network Vaccinology to explore the structure and function of biological networks responsible for vaccine-induced immunity. We demonstrate how the principles of graph theory can be used to identify modules of genes, proteins, and metabolites that are associated with innate and adaptive immune responses. Network vaccinology can be used to assess specific and shared molecular mechanisms of different types of vaccines, adjuvants, and routes of administration to direct rational design of the next generation of vaccines., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

The nitro-heterocyclic compound benznidazole (BZ) is the first-line drug for the treatment of Chagas disease, caused by the protozoan Trypanosoma cruzi. However, therapeutic failures are common for reasons that include the influences of parasite and host genetics, the effects of toxicity on adherence to treatment, and difficulties in demonstrating parasitological cure. To obtain information on the origin of the resistance to BZ and eliminate from the scenery the participation of the host, initially we mapped the susceptibility to the drug in thirteen species of seven genera of the family Trypanosomatidae. We verified that all Trypanosoma species are sensitive to low concentrations of the drug (IC 50 2.7 to 25 µM) while Non-Trypanosoma species are highly resistant to these concentrations. The two groups of parasites correspond to the major phylogenetic lineages of trypanosomatids. Next, we searched in the trypanosomatid genome databases homologs of two type-I nitroreductases (NTR-1 and OYE) and an ABC transporter (ABCG1) that have been associated with BZ resistance in T. cruzi. The predicted proteins were characterized regarding domains and used for phylogenetic analyses. Homologous NTR-1 genes were found in all trypanosomatids investigated and the structural characteristics of the enzyme suggest that it may be functional. OYE genes were absent in BZ-sensitive African trypanosomes, which excludes the participation of this enzyme in BZ bio-activation. Two copies of ABCG1 genes were observed in most BZ resistant species, while Trypanosoma species exhibit only one copy per haploid genome. Functional studies are required to verify the involvement of these genes in BZ resistance. In addition, since multiple mechanisms can contribute to BZ susceptibility, our study poses a range of organisms highly resistant to BZ in which these aspects can be investigated. Preliminary studies on BZ uptake indicate marked differences between BZ-sensitive and BZ-resistant species., (Copyright © 2019 Elsevier B.V. All rights reserved.)

Background: Trypanosoma conorhini and Trypanosoma rangeli, like Trypanosoma cruzi, are kinetoplastid protist parasites of mammals displaying divergent hosts, geographic ranges and lifestyles. Largely nonpathogenic T. rangeli and T. conorhini represent clades that are phylogenetically closely related to the T. cruzi and T. cruzi-like taxa and provide insights into the evolution of pathogenicity in those parasites. T. rangeli, like T. cruzi is endemic in many Latin American countries, whereas T. conorhini is tropicopolitan. T. rangeli and T. conorhini are exclusively extracellular, while T. cruzi has an intracellular stage in the mammalian host., Results: Here we provide the first comprehensive sequence analysis of T. rangeli AM80 and T. conorhini 025E, and provide a comparison of their genomes to those of T. cruzi G and T. cruzi CL, respectively members of T. cruzi lineages TcI and TcVI. We report de novo assembled genome sequences of the low-virulent T. cruzi G, T. rangeli AM80, and T. conorhini 025E ranging from ~ 21-25 Mbp, with ~ 10,000 to 13,000 genes, and for the highly virulent and hybrid T. cruzi CL we present a ~ 65 Mbp in-house assembled haplotyped genome with ~ 12,500 genes per haplotype. Single copy orthologs of the two T. cruzi strains exhibited ~ 97% amino acid identity, and ~ 78% identity to proteins of T. rangeli or T. conorhini. Proteins of the latter two organisms exhibited ~ 84% identity. T. cruzi CL exhibited the highest heterozygosity. T. rangeli and T. conorhini displayed greater metabolic capabilities for utilization of complex carbohydrates, and contained fewer retrotransposons and multigene family copies, i.e. trans-sialidases, mucins, DGF-1, and MASP, compared to T. cruzi., Conclusions: Our analyses of the T. rangeli and T. conorhini genomes closely reflected their phylogenetic proximity to the T. cruzi clade, and were largely consistent with their divergent life cycles. Our results provide a greater context for understanding the life cycles, host range expansion, immunity evasion, and pathogenesis of these trypanosomatids.

Small Heat-Shock Proteins (sHSPs) and other proteins bearing alpha-crystallin domains (ACD) participate in defense against heat and oxidative stress and play important roles in cell cycle, cytoskeleton dynamics, and immunological and pathological mechanisms in eukaryotes. However, little is known about these proteins in early-diverging lineages of protists such as the kinetoplastids. Here, ACD-like proteins (ACDp) were investigated in genomes of 61 species of 12 kinetoplastid genera, including Trypanosoma spp. (23 species of mammals, reptiles and frogs), Leishmania spp. (mammals and lizards), trypanosomatids of insects, Phytomonas spp. of plants, and bodonids. Comparison of ACDps based on domain architecture, predicted tertiary structure, phylogeny and genome organization reveals a kinetoplastid evolutionarily conserved repertoire, which diversified prior to trypanosomatid adaptation to parasitic life. We identified 9 ACDp orthologs classified in 8 families of TryACD: four previously recognized (HSP20, Tryp23A, Tryp23B and ATOM69), and four characterized for the first time in kinetoplastids (TryACDP, TrySGT1, TryDYX1C1 and TryNudC). A single copy of each ortholog was identified in each genome alongside TryNudC1/TrypNudC2 homologs and, overall, ACDPs were under strong selection pressures at main phylogenetic lineages. Transcripts of all ACDPs were identified across the life stages of T. cruzi, T. brucei and Leishmania spp., but proteomic profiles suggested that most ACDPs may be species- and stage-regulated. Our findings establish the basis for functional studies, and provided evolutionary and structural support for an underestimated repertoire of ACDps in the kinetoplastids., Competing Interests: The authors have declared that no competing interests exist.

Trypanosoma cruzi , the etiologic agent of Chagas disease, cycles through different life stages characterized by defined molecular traits associated with the proliferative or differentiation state. In particular, T. cruzi epimastigotes are the replicative forms that colonize the intestine of the Triatomine insect vector before entering the stationary phase that is crucial for differentiation into metacyclic trypomastigotes, which are the infective forms of mammalian hosts. The transition from proliferative exponential phase to quiescent stationary phase represents an important step that recapitulates the early molecular events of metacyclogenesis, opening new possibilities for understanding this process. In this study, we report a quantitative shotgun proteomic analysis of the T. cruzi epimastigote in the exponential and stationary growth phases. More than 3000 proteins were detected and quantified, highlighting the regulation of proteins involved in different subcellular compartments. Ribosomal proteins were upregulated in the exponential phase, supporting the higher replication rate of this growth phase. Autophagy-related proteins were upregulated in the stationary growth phase, indicating the onset of the metacyclogenesis process. Moreover, this study reports the regulation of N-terminally acetylated proteins during growth phase transitioning, adding a new layer of regulation to this process. Taken together, this study reports a proteome-wide rewiring during T. cruzi transit from the replicative exponential phase to the stationary growth phase, which is the preparatory phase for differentiation.

Trypanosoma rangeli and Trypanosoma cruzi are generalist trypanosomes sharing a wide range of mammalian hosts; they are transmitted by triatomine bugs, and are the only trypanosomes infecting humans in the Neotropics. Their origins, phylogenetic relationships, and emergence as human parasites have long been subjects of interest. In the present study, taxon-rich analyses (20 trypanosome species from bats and terrestrial mammals) using ssrRNA, glycosomal glyceraldehyde-3-phosphate dehydrogenase (gGAPDH), heat shock protein-70 (HSP70) and Spliced Leader RNA sequences, and multilocus phylogenetic analyses using 11 single copy genes from 15 selected trypanosomes, provide increased resolution of relationships between species and clades, strongly supporting two main sister lineages: lineage Schizotrypanum, comprising T. cruzi and bat-restricted trypanosomes, and Tra[Tve-Tco] formed by T. rangeli, Trypanosoma vespertilionis and Trypanosoma conorhini clades. Tve comprises European T. vespertilionis and African T. vespertilionis-like of bats and bat cimicids characterised in the present study and Trypanosoma sp. Hoch reported in monkeys and herein detected in bats. Tco included the triatomine-transmitted tropicopolitan T. conorhini from rats and the African NanDoum1 trypanosome of civet (carnivore). Consistent with their very close relationships, Tra[Tve-Tco] species shared highly similar Spliced Leader RNA structures that were highly divergent from those of Schizotrypanum. In a plausible evolutionary scenario, a bat trypanosome transmitted by cimicids gave origin to the deeply rooted Tra[Tve-Tco] and Schizotrypanum lineages, and bat trypanosomes of diverse genetic backgrounds jumped to new hosts. A long and independent evolutionary history of T. rangeli more related to Old World trypanosomes from bats, rats, monkeys and civets than to Schizotrypanum spp., and the adaptation of these distantly related trypanosomes to different niches of shared mammals and vectors, is consistent with the marked differences in transmission routes, life-cycles and host-parasite interactions, resulting in T. cruzi (but not T. rangeli) being pathogenic to humans., (Copyright © 2018 Australian Society for Parasitology. Published by Elsevier Ltd. All rights reserved.)

Background: Trypanosoma (Duttonella) vivax is a major pathogen of livestock in Africa and South America (SA), and genetic studies limited to small sampling suggest greater diversity in East Africa (EA) compared to both West Africa (WA) and SA., Methods: Multidimensional scaling and phylogenetic analyses of 112 sequences of the glycosomal glyceraldehyde phosphate dehydrogenase (gGAPDH) gene and 263 sequences of the internal transcribed spacer of rDNA (ITS rDNA) were performed to compare trypanosomes from tsetse flies from Gorongosa National Park and Niassa National Reserve of Mozambique (MZ), wild ungulates and livestock from EA, and livestock isolates from WA and SA., Results: Multidimensional scaling (MDS) supported Tvv (T. vivax) and TvL (T. vivax-like) evolutionary lineages: 1) Tvv comprises two main groups, TvvA/B (all SA and WA isolates plus some isolates from EA) and TvvC/D (exclusively from EA). The network revealed five ITS-genotypes within Tvv: Tvv1 (WA/EA isolates), Tvv2 (SA) and Tvv3-5 (EA). EA genotypes of Tvv ranged from highly related to largely different from WA/SA genotypes. 2) TvL comprises two gGAPDH-groups formed exclusively by EA sequences, TvLA (Tanzania/Kenya) and TvLB-D (MZ). This lineage contains more than 11 ITS-genotypes, seven forming the lineage TvL-Gorongosa that diverged from T. vivax Y486 enough to be identified as another species of the subgenus Duttonella. While gGAPDH sequences were fundamental for classification at the subgenus, major evolutionary lineages and species levels, ITS rDNA sequences permitted identification of known and novel genotypes., Conclusions: Our results corroborate a remarkable diversity of Duttonella trypanosomes in EA, especially in wildlife conservation areas, compared to the moderate diversity in WA. Surveys in wilderness areas in WA may reveal greater diversity. Biogeographical and phylogenetic data point to EA as the place of origin, diversification and spread of Duttonella trypanosomes across Africa, providing relevant insights towards the understanding of T. vivax evolutionary history.

We screened the genomes of a broad panel of kinetoplastid protists for genes encoding proteins associated with the RNA interference (RNAi) system using probes from the Argonaute (AGO1), Dicer1 (DCL1), and Dicer2 (DCL2) genes of Leishmania brasiliensis and Crithidia fasciculata. We identified homologs for all the three of these genes in the genomes of a subset of these organisms. However, several of these organisms lacked evidence for any of these genes, while others lacked only DCL2. The open reading frames encoding these putative proteins were structurally analyzed in silico. The alignments indicated that the genes are homologous with a high degree of confidence, and three-dimensional structural models strongly supported a functional relationship to previously characterized AGO1, DCL1, and DCL2 proteins. Phylogenetic analysis of these putative proteins showed that these genes, when present, evolved in parallel with other nuclear genes, arguing that the RNAi system genes share a common progenitor, likely across all Kinetoplastea. In addition, the genome segments bearing these genes are highly conserved and syntenic, even among those taxa in which they are absent. However, taxa in which these genes are apparently absent represent several widely divergent branches of kinetoplastids, arguing that these genes were independently lost at least six times in the evolutionary history of these organisms. The mechanisms responsible for the apparent coordinate loss of these RNAi system genes independently in several lineages of kinetoplastids, while being maintained in other related lineages, are currently unknown.

Background: Proline racemase (PRAC) enzymes of Trypanosoma cruzi (TcPRAC), the agent of Chagas disease, and Trypanosoma vivax (TvPRAC), the agent of livestock trypanosomosis, have been implicated in the B-cells polyclonal activation contributing to immunosuppression and the evasion of host defences. The similarity to prokaryotic PRAC and the absence in Trypanosoma brucei and Trypanosoma congolense have raised many questions about the origin, evolution, and functions of trypanosome PRAC (TryPRAC) enzymes., Findings: We identified TryPRAC homologs as single copy genes per haploid genome in 12 of 15 Trypanosoma species, including T. cruzi and T. cruzi marinkellei, T. dionisii, T. erneyi, T. rangeli, T. conorhini and T. lewisi, all parasites of mammals. Polymorphisms in TcPRAC genes matched T. cruzi genotypes: TcI-TcIV and Tcbat have unique genes, while the hybrids TcV and TcVI contain TcPRACA and TcPRACB from parental TcII and TcIII, respectively. PRAC homologs were identified in trypanosomes from anurans, snakes, crocodiles, lizards, and birds. Most trypanosomes have intact PRAC genes. T. rangeli possesses only pseudogenes, maybe in the process of being lost. T. brucei, T. congolense and their allied species, except the more distantly related T. vivax, have completely lost PRAC genes., Conclusions: The genealogy of TryPRAC homologs supports an evolutionary history congruent with the Trypanosoma phylogeny. This finding, together with the synteny of PRAC loci, the relationships with prokaryotic PRAC inferred by taxon-rich phylogenetic analysis, and the absence in trypanosomatids of any other genera or in bodonids or euglenids suggest that a common ancestor of Trypanosoma gained PRAC gene by a single and ancient horizontal gene transfer (HGT) from a Firmicutes bacterium more closely related to Gemella and other species of Bacilli than to Clostridium as previously suggested. Our broad phylogenetic study allowed investigation of TryPRAC evolution over long and short timescales. TryPRAC genes diverged to become species-specific and genotype-specific for T. cruzi and T. rangeli, with resulting genealogies congruent with those obtained using vertically inherited genes. The inventory of TryPRAC genes described here is the first step toward the understanding of the roles of PRAC enzymes in trypanosomes differing in life cycles, virulence, and infection and immune evasion strategies.

Trypanosoma congolense is the most important agent of nagana, a wasting livestock trypanosomosis in sub-Saharan Africa. This species is a complex of three subgroups (Savannah, Forest and Kilifi) that differ in virulence, pathogenicity, drug resistance, vectors, and geographical distribution. Congopain, the major Cathepsin L-like cysteine protease (CP2) of T. congolense, has been extensively investigated as a pathogenic factor and target for drugs and vaccines, but knowledge about this enzyme is mostly restricted to the reference strain IL3000, which belongs to the Savannah subgroup. In this work we compared sequences of congopain genes from IL3000 genome database and isolates of the three subgroups of T. congolense. Results demonstrated that the congopain genes diverged into three subclades consistent with the three subgroups within T. congolense. Laboratory and field isolates of Savannah exhibited a highly polymorphic repertoire both inter- and intra-isolates: sequences sharing the archetypical catalytic triad clustered into SAV1-SAV3 groups, whereas polymorphic sequences that, in general, exhibited unusual catalytic triad (variants) assigned to SAV4 or not assigned to any group. Congopain homologous genes from Forest and Kilifi isolates showed, respectively, moderate and limited diversity. In the phylogenetic tree based on congopain and homologues, Savannah was closer to Forest than to Kilifi. All T. congolense subgroup nested into a single clade, which together with the sister clade formed by homologues from Trypanosoma simiae and Trypanosoma godfreyi formed a clade supporting the subgenus Nannomonas. A single PCR targeting congopain sequences was developed for the diagnosis of T. congolense isolates of the three subgroups. Our findings demonstrated that congopain genes are valuable targets for the diagnosis, genotyping, and phylogenetic and taxonomic inferences among T. congolense isolates and other members of the subgenus Nannomonas., (Copyright © 2014 Elsevier B.V. All rights reserved.)

Background: Trypanosomatids of the genera Angomonas and Strigomonas live in a mutualistic association characterized by extensive metabolic cooperation with obligate endosymbiotic Betaproteobacteria. However, the role played by the symbiont has been more guessed by indirect means than evidenced. Symbiont-harboring trypanosomatids, in contrast to their counterparts lacking symbionts, exhibit lower nutritional requirements and are autotrophic for essential amino acids. To evidence the symbiont's contributions to this autotrophy, entire genomes of symbionts and trypanosomatids with and without symbionts were sequenced here., Results: Analyses of the essential amino acid pathways revealed that most biosynthetic routes are in the symbiont genome. By contrast, the host trypanosomatid genome contains fewer genes, about half of which originated from different bacterial groups, perhaps only one of which (ornithine cyclodeaminase, EC:4.3.1.12) derived from the symbiont. Nutritional, enzymatic, and genomic data were jointly analyzed to construct an integrated view of essential amino acid metabolism in symbiont-harboring trypanosomatids. This comprehensive analysis showed perfect concordance among all these data, and revealed that the symbiont contains genes for enzymes that complete essential biosynthetic routes for the host amino acid production, thus explaining the low requirement for these elements in symbiont-harboring trypanosomatids. Phylogenetic analyses show that the cooperation between symbionts and their hosts is complemented by multiple horizontal gene transfers, from bacterial lineages to trypanosomatids, that occurred several times in the course of their evolution. Transfers occur preferentially in parts of the pathways that are missing from other eukaryotes., Conclusion: We have herein uncovered the genetic and evolutionary bases of essential amino acid biosynthesis in several trypanosomatids with and without endosymbionts, explaining and complementing decades of experimental results. We uncovered the remarkable plasticity in essential amino acid biosynthesis pathway evolution in these protozoans, demonstrating heavy influence of horizontal gene transfer events, from Bacteria to trypanosomatid nuclei, in the evolution of these pathways.

The congenital transmission of Chagas disease is associated with an increase in parasitemia during pregnancy, maternal and fetal immunity, and populations of Trypanosoma cruzi. In this study, the biological behavior of TcI and TcV (isolated from a human congenital case) strains and their potential for experimental congenital transmission were evaluated in female BALB/C mice. Parasitemia was estimated by fresh blood examination, semiquantitative microhematocrit, and hemoculture, while congenital transmission was evaluated by culture in the liver infusion tryptose medium and by polymerase chain reaction (PCR) of the pups' tissues on postnatal day 7 and of the pups' blood sample at 30 days after birth. Infection was detected in 100 % of the females. Both strains showed subpatent parasitemia, which was higher for TcV infection. The presence of amastigote nest was detected only in an animal infected with TcI. The inflammatory process was more frequent (p = 0.001) in the tissues of the animals infected with TcV (58.6 %) than TcI (31.1 %). The fertility rates of females mated after 35 days postinfection were similar (90 % for TcV, 88.9 % for TcI; p = 0.938). Parasitemia did not change during pregnancy. The average number of pups/female was greater (p = 0.03) in mice with TcV infection (8.30) than in those with TcI infection (4.78). Congenital transmission was detected exclusively by PCR in 50.9 % of the pups, 46.6 % for TcV and 58.1 % for TcI. The PCR positivity for TcI was higher in the blood than in the tissue (p = 0.003). These results demonstrate the T. cruzi experimental congenital infection associated with subpatent maternal parasitemia of TcI and TcV.