Search

Abstract Many molecular mechanisms that lead to the host antibody response to COVID-19 vaccines remain largely unknown. In this study, we used serum antibody detection combined with whole blood RNA-based transcriptome analysis to investigate variability in vaccine response in healthy recipients of a booster (third) dose schedule of the mRNA BNT162b2 vaccine against COVID-19. The cohort was divided into two groups: (1) low-stable individuals, with antibody concentration anti-SARS-CoV IgG S1 below 0.4 percentile at 180 days after boosting vaccination; and (2) high-stable individuals, with antibody values greater than 0.6 percentile of the range in the same period (median 9525 [185–80,000] AU/mL). Differential gene expression, expressed single nucleotide variants and insertions/deletions, differential splicing events, and allelic imbalance were explored to broaden our understanding of the immune response sustenance. Our analysis revealed a differential expression of genes with immunological functions in individuals with low antibody titers, compared to those with higher antibody titers, underscoring the fundamental importance of the innate immune response for boosting immunity. Our findings also provide new insights into the determinants of the immune response variability to the SARS-CoV-2 mRNA vaccine booster, highlighting the significance of differential splicing regulatory mechanisms, mainly concerning HLA alleles, in delineating vaccine immunogenicity.

Abstract Background Phylogenetic gaps of public databases of reference sequences are a major obstacle for comparative genomics and management of marine resources, particularly in the Global South, where economically important fisheries and conservation flagship species often lack closely-related references. We applied target-enrichment to obtain complete mitochondrial genomes of marine ichthyofauna from the Brazilian coast selected based on economic significance, conservation status and lack of phylogenetically-close references. These included sardines (Dorosomatidae, Alosidae), mackerels (Scombridae) croakers (Sciaenidae), groupers (Epinephelidae) and snappers (Lutjanidae). Results Custom baits were designed to enrich mitochondrial DNA across a broad phylogenetic range of fishes. Sequencing generated approximately 100k reads per sample, which were assembled in a total of 70 complete mitochondrial genomes and include fifty-two new additions to GenBank, including five species with no previous mitochondrial data. Departures from the typical gene content and order occurred in only three taxa and mostly involved tRNA gene duplications. Start-codons for all genes, except Cytochrome C Oxidase subunit I (COI), were consistently ATG, whilst a wide range of stop-codons deviated from the prevailing TAA. Phylogenetic analysis confirmed assembly accuracy and revealed signs of cryptic diversification within the Mullus genus. Lineage delimitation methods using Sardinella aurita and S. brasiliensis mitochondrial genomes support a single Operational Taxonomic Unit. Conclusions Target enrichment was highly efficient, providing complete novel mitochondrial genomes with little sequencing effort. These sequences are deposited in public databases to enable subsequent studies in population genetics and adaptation of Latin American fish species and serve as a vital resource for conservation and management programs that rely on molecular data for species and genus-level identification.

Abstract Background Multisystem Inflammatory Syndrome in Children (MIS-C) is a life-threatening complication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, which manifests as a hyper inflammatory process with multiorgan involvement in predominantly healthy children in the weeks following mild or asymptomatic coronavirus disease 2019 (COVID-19). However, host monogenic predisposing factors to MIS-C remain elusive. Methods Herein, we used whole exome sequencing (WES) on 16 MIS-C Brazilian patients to identify single nucleotide/InDels variants as predisposition factors associated with MIS-C. Results We identified ten very rare variants in eight genes (FREM1, MPO, POLG, C6, C9, ABCA4, ABCC6, and BSCL2) as the most promising candidates to be related to a higher risk of MIS-C development. These variants may propitiate a less effective immune response to infection or trigger the inflammatory response or yet a delayed hyperimmune response to SARS-CoV-2. Protein–Protein Interactions (PPIs) among the products of the mutated genes revealed an integrated network, enriched for immune and inflammatory response mechanisms with some of the direct partners representing gene products previously associated with MIS-C and Kawasaki disease (KD). In addition, the PPIs direct partners are also enriched for COVID-19-related gene sets. HLA alleles prediction from WES data allowed the identification of at least one risk allele in 100% of the MIS-C patients. Conclusions This study is the first to explore host MIS-C-associated variants in a Latin American admixed population. Besides expanding the spectrum of MIS-C-associated variants, our findings highlight the relevance of using WES for characterising the genetic interindividual variability associated with COVID-19 complications and ratify the presence of overlapping/convergent mechanisms among MIS-C, KD and COVID-19, crucial for future therapeutic management.

Measurement(s) Metagenome Technology Type(s) Illumina Sequencing

Background Technological advances involving RNA-Seq and Bioinformatics allow quantifying the transcriptional levels of genes in cells, tissues, and cell lines, permitting the identification of Differentially Expressed Genes (DEGs). DESeq2 and edgeR are well-established computational tools used for this purpose and they are based upon generalized linear models (GLMs) that consider only fixed effects in modeling. However, the inclusion of random effects reduces the risk of missing potential DEGs that may be essential in the context of the biological phenomenon under investigation. The generalized linear mixed models (GLMM) can be used to include both effects. Methods We present DEGRE (Differentially Expressed Genes with Random Effects), a user-friendly tool capable of inferring DEGs where fixed and random effects on individuals are considered in the experimental design of RNA-Seq research. DEGRE preprocesses the raw matrices before fitting GLMMs on the genes and the derived regression coefficients are analyzed using the Wald statistical test. DEGRE offers the Benjamini-Hochberg or Bonferroni techniques for P-value adjustment. Results The datasets used for DEGRE assessment were simulated with known identification of DEGs. These have fixed effects, and the random effects were estimated and inserted to measure the impact of experimental designs with high biological variability. For DEGs’ inference, preprocessing effectively prepares the data and retains overdispersed genes. The biological coefficient of variation is inferred from the counting matrices to assess variability before and after the preprocessing. The DEGRE is computationally validated through its performance by the simulation of counting matrices, which have biological variability related to fixed and random effects. DEGRE also provides improved assessment measures for detecting DEGs in cases with higher biological variability. We show that the preprocessing established here effectively removes technical variation from those matrices. This tool also detects new potential candidate DEGs in the transcriptome data of patients with bipolar disorder, presenting a promising tool to detect more relevant genes. Conclusions DEGRE provides data preprocessing and applies GLMMs for DEGs’ inference. The preprocessing allows efficient remotion of genes that could impact the inference. Also, the computational and biological validation of DEGRE has shown to be promising in identifying possible DEGs in experiments derived from complex experimental designs. This tool may help handle random effects on individuals in the inference of DEGs and presents a potential for discovering new interesting DEGs for further biological investigation.

Abstract Background X-linked agammaglobulinemia (XLA) is an Inborn Errors of Immunity (IEI) characterized by pan-hypogammaglobulinemia and low numbers of B lymphocytes due to mutations in BTK gene. Usually, XLA patients are not susceptible to respiratory tract infections by viruses and do not present interstitial lung disease (ILD) such as bronchiolitis obliterans (BO) as a consequence of acute or chronic bacterial infections of the respiratory tract. Although many pathogenic variants have already been described in XLA, the heterogeneous clinical presentations in affected patients suggest a more complex genetic landscape underlying this disorder. Case presentation We report two pediatric cases from male siblings with X-Linked Agammaglobulinemia and bronchiolitis obliterans, a phenotype not often observed in XLA phenotype. The whole-exome sequencing (WES) analysis showed a rare hemizygous missense variant NM_000061.2(BTK):c.1751G>A(p.Gly584Glu) in BTK gene of both patients. We also identified a gain-of-function mutation in TGFβ1 (rs1800471) previously associated with transforming growth factor-beta1 production, fibrotic lung disease, and graft fibrosis after lung transplantation. TGFβ1 plays a key role in the regulation of immune processes and inflammatory response associated with pulmonary impairment. Conclusions Our report illustrates a possible role for WES in patients with known inborn errors of immunity, but uncommon clinical presentations, providing a personalized understanding of genetic basis, with possible implications in the identification of potential treatments, and prognosis for patients and their families.

A new variant of concern for SARS-CoV-2, Omicron (B.1.1.529), was designated by the World Health Organization on November 26, 2021. This study analyzed the viral genome sequencing data of 108 samples collected from patients infected with Omicron. First, we found that the enrichment efficiency of viral nucleic acids was reduced due to mutations in the region where the primers anneal to. Second, the Omicron variant possesses an excessive number of mutations compared to other variants circulating at the same time (median: 62 vs. 45), especially in the Spike gene. Mutations in the Spike gene confer alterations in 32 amino acid residues, more than those observed in other SARS-CoV-2 variants. Moreover, a large number of nonsynonymous mutations occur in the codons for the amino acid residues located on the surface of the Spike protein, which could potentially affect the replication, infectivity, and antigenicity of SARS-CoV-2. Third, there are 53 mutations between the Omicron variant and its closest sequences available in public databases. Many of these mutations were rarely observed in public databases and had a low mutation rate. In addition, the linkage disequilibrium between these mutations was low, with a limited number of mutations concurrently observed in the same genome, suggesting that the Omicron variant would be in a different evolutionary branch from the currently prevalent variants. To improve our ability to detect and track the source of new variants rapidly, it is imperative to further strengthen genomic surveillance and data sharing globally in a timely manner.

IntroductionCell entry of SARS-CoV-2 causes genome-wide disruption of the transcriptional profiles of genes and biological pathways involved in the pathogenesis of COVID-19. Expression allelic imbalance is characterized by a deviation from the Mendelian expected 1:1 expression ratio and is an important source of allele-specific heterogeneity. Expression allelic imbalance can be measured by allele-specific expression analysis (ASE) across heterozygous informative expressed single nucleotide variants (eSNVs). ASE reflects many regulatory biological phenomena that can be assessed by combining genome and transcriptome information. ASE contributes to the interindividual variability associated with the disease. We aim to estimate the transcriptome-wide impact of SARS-CoV-2 infection by analyzing eSNVs.MethodsWe compared ASE profiles in the human lung cell lines Calu-3, A459, and H522 before and after infection with SARS-CoV-2 using RNA-Seq experiments.ResultsWe identified 34 differential ASE (DASE) sites in 13 genes (HLA-A, HLA-B, HLA-C, BRD2, EHD2, GFM2, GSPT1, HAVCR1, MAT2A, NQO2, SUPT6H, TNFRSF11A, UMPS), all of which are enriched in protein binding functions and play a role in COVID-19. Most DASE sites were assigned to the MHC class I locus and were predominantly upregulated upon infection. DASE sites in the MHC class I locus also occur in iPSC-derived airway epithelium basal cells infected with SARS-CoV-2. Using an RNA-Seq haplotype reconstruction approach, we found DASE sites and adjacent eSNVs in phase (i.e., predicted on the same DNA strand), demonstrating differential haplotype expression upon infection. We found a bias towards the expression of the HLA alleles with a higher binding affinity to SARS-CoV-2 epitopes.DiscussionIndependent of gene expression compensation, SARS-CoV-2 infection of human lung cell lines induces transcriptional allelic switching at the MHC loci. This suggests a response mechanism to SARS-CoV-2 infection that swaps HLA alleles with poor epitope binding affinity, an expectation supported by publicly available proteome data.

Abstract Background The rice weevil Sitophilus oryzae is one of the most important agricultural pests, causing extensive damage to cereal in fields and to stored grains. S. oryzae has an intracellular symbiotic relationship (endosymbiosis) with the Gram-negative bacterium Sodalis pierantonius and is a valuable model to decipher host-symbiont molecular interactions. Results We sequenced the Sitophilus oryzae genome using a combination of short and long reads to produce the best assembly for a Curculionidae species to date. We show that S. oryzae has undergone successive bursts of transposable element (TE) amplification, representing 72% of the genome. In addition, we show that many TE families are transcriptionally active, and changes in their expression are associated with insect endosymbiotic state. S. oryzae has undergone a high gene expansion rate, when compared to other beetles. Reconstruction of host-symbiont metabolic networks revealed that, despite its recent association with cereal weevils (30 kyear), S. pierantonius relies on the host for several amino acids and nucleotides to survive and to produce vitamins and essential amino acids required for insect development and cuticle biosynthesis. Conclusions Here we present the genome of an agricultural pest beetle, which may act as a foundation for pest control. In addition, S. oryzae may be a useful model for endosymbiosis, and studying TE evolution and regulation, along with the impact of TEs on eukaryotic genomes.

Abstract Background The Rhizobiales (Proteobacteria) order is an abundant and diverse group of microorganisms, being extensively studied for its lifestyle based on the association with plants, animals, and humans. New studies have demonstrated that the last common ancestor (LCA) of Rhizobiales had a free-living lifestyle, but the phylogenetic and metabolism characterization of basal lineages remains unclear. Here, we used a high-resolution phylogenomic approach to test the monophyly of the Aestuariivirgaceae family, a new taxonomic group of Rhizobiales. Furthermore, a deep metabolic investigation provided an overview of the main functional traits that can be associated with its lifestyle. We hypothesized that the presence of pathways (e.g., Glycolysis/Gluconeogenesis) and the absence of pathogenic genes would be associated with a free-living lifestyle in Aestuariivirgaceae. Results Using high-resolution phylogenomics approaches, our results revealed a clear separation of Aestuariivirgaceae into a distinct clade of other Rhizobiales family, suggesting a basal split early group and corroborate the monophyly of this group. A deep functional annotation indicated a metabolic versatility, which includes putative genes related to sugar degradation and aerobic respiration. Furthermore, many of these traits could reflect a basal metabolism and adaptations of Rhizobiales, as such the presence of Glycolysis/Gluconeogenesis pathway and the absence of pathogenicity genes, suggesting a free-living lifestyle in the Aestuariivirgaceae members. Conclusions Aestuariivirgaceae (Rhizobiales) family is a monophyletic taxon of the Rhizobiales with a free-living lifestyle and a versatile metabolism that allows these microorganisms to survive in the most diverse microbiomes, demonstrating their adaptability to living in systems with different conditions, such as extremely cold environments to tropical rivers.

ABSTRACT The epidemiology of antimicrobial resistance (AMR) is complex, with multiple interfaces (human-animal-environment). In this context, One Health surveillance is essential for understanding the distribution of microorganisms and antimicrobial resistance genes (ARGs). This report describes a multicentric study undertaken to evaluate the bacterial communities and resistomes of food-producing animals (cattle, poultry, and swine) and healthy humans sampled simultaneously from five Brazilian regions. Metagenomic analysis showed that a total of 21,029 unique species were identified in 107 rectal swabs collected from distinct hosts, the highest numbers of which belonged to the domain Bacteria, mainly Ruminiclostridium spp. and Bacteroides spp., and the order Enterobacterales. We detected 405 ARGs for 12 distinct antimicrobial classes. Genes encoding antibiotic-modifying enzymes were the most frequent, followed by genes related to target alteration and efflux systems. Interestingly, carbapenemase-encoding genes such as blaAIM-1, blaCAM-1, blaGIM-2, and blaHMB-1 were identified in distinct hosts. Our results revealed that, in general, the bacterial communities from humans were present in isolated clusters, except for the Northeastern region, where an overlap of the bacterial species from humans and food-producing animals was observed. Additionally, a large resistome was observed among all analyzed hosts, with emphasis on the presence of carbapenemase-encoding genes not previously reported in Latin America. IMPORTANCE Humans and food production animals have been reported to be important reservoirs of antimicrobial resistance (AMR) genes (ARGs). The frequency of these multidrug-resistant (MDR) bacteria tends to be higher in low- and middle-income countries (LMICs), due mainly to a lack of public health policies. Although studies on AMR in humans or animals have been carried out in Brazil, this is the first multicenter study that simultaneously collected rectal swabs from humans and food-producing animals for metagenomics. Our results indicate high microbial diversity among all analyzed hosts, and several ARGs for different antimicrobial classes were also found. As far as we know, we have detected for the first time ARGs encoding carbapenemases, such as blaAIM-1, blaCAM-1, blaGIM-2, and blaHMB-1, in Latin America. Thus, our results support the importance of metagenomics as a tool to track the colonization of food-producing animals and humans by antimicrobial-resistant bacteria. In addition, a network surveillance system called GUARANI, created for this study, is ready to be expanded and to collect additional data.

Abstract Background Cestoda is a class of endoparasitic worms in the flatworm phylum (Platyhelminthes). During the course of their evolution cestodes have evolved some interesting aspects, such as their increased reproductive capacity. In this sense, they have serial repetition of their reproductive organs in the adult stage, which is often associated with external segmentation in a developmental process called strobilation. However, the molecular basis of strobilation is poorly understood. To assess this issue, an evolutionary comparative study among strobilated and non-strobilated flatworm species was conducted to identify genes and proteins related to the strobilation process. Results We compared the genomic content of 10 parasitic platyhelminth species; five from cestode species, representing strobilated parasitic platyhelminths, and five from trematode species, representing non-strobilated parasitic platyhelminths. This dataset was used to identify 1813 genes with orthologues that are present in all cestode (strobilated) species, but absent from at least one trematode (non-strobilated) species. Development-related genes, along with genes of unknown function (UF), were then selected based on their transcriptional profiles, resulting in a total of 34 genes that were differentially expressed between the larval (pre-strobilation) and adult (strobilated) stages in at least one cestode species. These 34 genes were then assumed to be strobilation related; they included 12 encoding proteins of known function, with 6 related to the Wnt, TGF-β/BMP, or G-protein coupled receptor signaling pathways; and 22 encoding UF proteins. In order to assign function to at least some of the UF genes/proteins, a global gene co-expression analysis was performed for the cestode species Echinococcus multilocularis. This resulted in eight UF genes/proteins being predicted as related to developmental, reproductive, vesicle transport, or signaling processes. Conclusions Overall, the described in silico data provided evidence of the involvement of 34 genes/proteins and at least 3 developmental pathways in the cestode strobilation process. These results highlight on the molecular mechanisms and evolution of the cestode strobilation process, and point to several interesting proteins as potential developmental markers and/or targets for the development of novel antihelminthic drugs.

Objectives: Enterobacter hormaechei is an important causative agent of severe infections in critically ill patients. Aminoglycosides are among the main antibiotics for the treatment of E. hormaechei infections, however the development of antimicrobial resistance is an increasing problem. RmtG is a 16S rRNA methyltransferase, a class of enzymes conferring high-level resistance to clinically relevant aminoglycosides. The aim of this study was to characterise the full genetic context of plasmids harbouring the rmtG gene in two aminoglycoside-resistant E. hormaechei isolated in Brazil. Methods: ThermtG-harbouring plasmids were transferred to an Escherichia coli J53 recipient strain and were fully sequenced using a MiSeq sequencing system. Complete genome assemblies were accomplished using a combination of Newbler v.3.0, SPAdes 3.10.0 and phrap/cross_match programs. Plasmid sequences were annotated using RAST server and were then manually curated using BLAST databases and ISfinder. Easyfig 2.0 was used to map and compare regions of interest containing rmtG in both plasmids. Results: Both isolates carried thermtG gene on an IncA/C plasmid of ˜152 kb and ˜235 kb, respectively, associated with a Tn3 transposon. The plasmids contain a transfer region as well as genes involved in plasmid stability and resistance to β-lactams, sulfonamides and quaternary ammonium compounds. One of the plasmids also carried the mrk operon encoding type 3 fimbriae. Conclusion: This first detection ofrmtG in E. hormaechei supports the ability for horizontal transfer. The location in complex genetic platforms carried by Tn3 transposons in IncA/C plasmids may facilitate dissemination to other Gram-negative pathogens, further limiting treatment options.

The tropism of Zika virus (ZIKV) has been described in the nervous system, blood, placenta, thymus, and skeletal muscle. We investigated the mechanisms of skeletal muscle susceptibility to ZIKV using an in vitro model of human skeletal muscle myogenesis, in which myoblasts differentiate into myotubes. Myoblasts were permissive to ZIKV infection, generating productive viral particles, while myotubes controlled ZIKV replication. To investigate the underlying mechanisms, we used gene expression profiling. First, we assessed gene changes in myotubes compared with myoblasts in the model without infection. As expected, we observed an increase in genes and pathways related to the contractile muscle system in the myotubes, a reduction in processes linked to proliferation, migration and cytokine production, among others, confirming the myogenic capacity of our system in vitro. A comparison between non-infected and infected myoblasts revealed more than 500 differentially expressed genes (DEGs). In contrast, infected myotubes showed almost 2,000 DEGs, among which we detected genes and pathways highly or exclusively expressed in myotubes, including those related to antiviral and innate immune responses. Such gene modulation could explain our findings showing that ZIKV also invades myotubes but does not replicate in these differentiated cells. In conclusion, we showed that ZIKV largely (but differentially) disrupts gene expression in human myoblasts and myotubes. Identifying genes involved in myotube resistance can shed light on potential antiviral mechanisms against ZIKV infection.

BACKGROUND During routine Coronavirus disease 2019 (COVID-19) diagnosis, an unusually high viral load was detected by reverse transcription real-time polymerase chain reaction (RT-qPCR) in a nasopharyngeal swab sample collected from a patient with respiratory and neurological symptoms who rapidly succumbed to the disease. Therefore we sought to characterise the infection. OBJECTIVES We aimed to determine and characterise the etiological agent responsible for the poor outcome. METHODS Classical virological methods, such as plaque assay and plaque reduction neutralisation test combined with amplicon-based sequencing, as well as a viral metagenomic approach, were performed to characterise the etiological agents of the infection. FINDINGS Plaque assay revealed two distinct plaque phenotypes, suggesting either the presence of two severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) strains or a productive coinfection of two different species of virus. Amplicon-based sequencing did not support the presence of any SARS-CoV-2 genetic variants that would explain the high viral load and suggested the presence of a single SARS-CoV-2 strain. Nonetheless, the viral metagenomic analysis revealed that Coronaviridae and Herpesviridae were the predominant virus families within the sample. This finding was confirmed by a plaque reduction neutralisation test and PCR. MAIN CONCLUSIONS We characterised a productive coinfection of SARS-CoV-2 and Herpes simplex virus 1 (HSV-1) in a patient with severe symptoms that succumbed to the disease. Although we cannot establish the causal relationship between the coinfection and the severity of the clinical case, this work serves as a warning for future studies focused on the interplay between SARS-CoV-2 and HSV-1 coinfection and COVID-19 severity.

Predicting the physical or functional associations through protein-protein interactions (PPIs) represents an integral approach for inferring novel protein functions and discovering new drug targets during repositioning analysis. Recent advances in high-throughput data generation and multi-omics techniques have enabled large-scale PPI predictions, thus promoting several computational methods based on different levels of biological evidence. However, integrating multiple results and strategies to optimize, extract interaction features automatically and scale up the entire PPI prediction process is still challenging. Most procedures do not offer an in-silico validation process to evaluate the predicted PPIs. In this context, this paper presents the PredPrIn scientific workflow that enables PPI prediction based on multiple lines of evidence, including the structure, sequence, and functional annotation categories, by combining boosting and stacking machine learning techniques. We also present a pipeline (PPIVPro) for the validation process based on cellular co-localization filtering and a focused search of PPI evidence on scientific publications. Thus, our combined approach provides means to extensive scale training or prediction of new PPIs and a strategy to evaluate the prediction quality. PredPrIn and PPIVPro are publicly available at https://github.com/YasCoMa/predprin and https://github.com/YasCoMa/ppi_validation_process.

Congenital Zika Syndrome (CZS) is a critical illness with a wide range of severity caused by Zika virus (ZIKV) infection during pregnancy. Life-threatening neurodevelopmental dysfunctions are among the most common phenotypes observed in affected newborns. Risk factors that contribute to susceptibility and response to ZIKV infection may be related to the virus itself, the environment, and maternal genetic background. Nevertheless, the newborn's genetic contribution to the critical illness is still not elucidated. Here, we aimed to identify possible genetic variants as well as relevant biological pathways that might be associated with CZS phenotypes. For this purpose, we performed a whole-exome sequencing in 40 children born to women with confirmed exposure to ZIKV during pregnancy. We investigated the occurrence of rare harmful single-nucleotide variants (SNVs) possibly associated with inborn errors in genes ontologically related to CZS phenotypes. Moreover, an exome-wide association analysis was also performed using a case-control design (29 CZS cases and 11 controls), for both common and rare variants. Five out of the 29 CZS patients harbored known pathogenic variants likely to contribute to mild to severe manifestations observed. Approximately, 30% of affected individuals carried at least one pathogenic or likely pathogenic SNV in genes candidates to play a role in CZS. Our common variant association analysis detected a suggestive protective effect of the rs2076469 in DISP3 gene (p-value: 1.39 x 10-5). The IL12RB2 gene (p-value: 2.18x10-11) also showed an unusual distribution of nonsynonymous rare SNVs in control samples. Finally, genes harboring harmful variants are involved in processes related to CZS phenotypes such as neurological development and immunity. Therefore, both rare and common variations may be likely to contribute as the underlying genetic cause of CZS susceptibility. The variations and pathways identified in this study may also have implications for the development of therapeutic strategies in the future.

Abstract Background The Metarhizium genus harbors important entomopathogenic fungi. These species have been widely explored as biological control agents, and strategies to improve the fungal virulence are under investigation. Thus, the interaction between Metarhizium species and susceptible hosts have been explored employing different methods in order to characterize putative virulence determinants. However, the impact of epigenetic modulation on the infection cycle of Metarhizium is still an open topic. Among the different epigenetic modifications, DNA methylation of cytosine bases is an important mechanism to control gene expression in several organisms. To better understand if DNA methylation can govern Metarhizium-host interactions, the genome-wide DNA methylation profile of Metarhizium anisopliae was explored in two conditions: tick mimicked infection and a saprophytic-like control. Results Using a genome wide DNA methylation profile based on bisulfite sequencing (BS-Seq), approximately 0.60% of the total cytosines were methylated in saprophytic-like condition, which was lower than the DNA methylation level (0.89%) in tick mimicked infection condition. A total of 670 mRNA genes were found to be putatively methylated, with 390 mRNA genes uniquely methylated in the tick mimicked infection condition. GO terms linked to response to stimuli, cell wall morphogenesis, cytoskeleton morphogenesis and secondary metabolism biosynthesis were over-represented in the tick mimicked infection condition, suggesting that energy metabolism is directed towards the regulation of genes associated with infection. However, recognized virulence determinants known to be expressed at distinct infection steps, such as the destruxin backbone gene and the collagen-like protein gene Mcl1, were found methylated, suggesting that a dynamic pattern of methylation could be found during the infectious process. These results were further endorsed employing RT-qPCR from cultures treated or not with the DNA methyltransferase inhibitor 5-Azacytidine. Conclusions The set of genes here analyzed focused on secondary metabolites associated genes, known to be involved in several processes, including virulence. The BS-Seq pipeline and RT-qPCR analysis employing 5-Azacytidine led to identification of methylated virulence genes in M. anisopliae. The results provided evidences that DNA methylation in M. anisopliae comprises another layer of gene expression regulation, suggesting a main role of DNA methylation regulating putative virulence determinants during M. anisopliae infection cycle.

Abstract Background Typhoid fever, caused by Salmonella Typhi, follows a fecal-oral transmission route and is a major global public health concern, especially in developing countries like Bangladesh. Increasing emergence of antimicrobial resistance (AMR) is a serious issue; the list of treatments for typhoid fever is ever-decreasing. In addition to IncHI1-type plasmids, Salmonella genomic island (SGI) 11 has been reported to carry AMR genes. Although reports suggest a recent reduction in multidrug resistance (MDR) in the Indian subcontinent, the corresponding genomic changes in the background are unknown. Results Here, we assembled and annotated complete closed chromosomes and plasmids for 73 S. Typhi isolates using short-length Illumina reads. S. Typhi had an open pan-genome, and the core genome was smaller than previously reported. Considering AMR genes, we identified five variants of SGI11, including the previously reported reference sequence. Five plasmids were identified, including the new plasmids pK91 and pK43; pK43and pHCM2 were not related to AMR. The pHCM1, pPRJEB21992 and pK91 plasmids carried AMR genes and, along with the SGI11 variants, were responsible for resistance phenotypes. pK91 also contained qnr genes, conferred high ciprofloxacin resistance and was related to the H58-sublineage Bdq, which shows the same phenotype. The presence of plasmids (pHCM1 and pK91) and SGI11 were linked to two H58-lineages, Ia and Bd. Loss of plasmids and integration of resistance genes in genomic islands could contribute to the fitness advantage of lineage Ia isolates. Conclusions Such events may explain why lineage Ia is globally widespread, while the Bd lineage is locally restricted. Further studies are required to understand how these S. Typhi AMR elements spread and generate new variants. Preventive measures such as vaccination programs should also be considered in endemic countries; such initiatives could potentially reduce the spread of AMR.

Arboviruses pose a major threat throughout the world and represent a great burden in tropical countries of South America. Although generally associated with moderate febrile illness, in more severe cases they can lead to neurological outcomes, such as encephalitis, Guillain-Barré syndrome, and Congenital Syndromes. In this context astrocytes play a central role in production of inflammatory cytokines, regulation of extracellular matrix, and control of glutamate driven neurotoxicity in the central nervous system. Here, we presented a comprehensive genome-wide transcriptome analysis of human primary astrocytes infected with Chikungunya, Mayaro, Oropouche, or Zika viruses. Analyses of differentially expressed genes (DEGs), pathway enrichment, and interactomes have shown that Alphaviruses up-regulated genes related to elastic fiber formation and N-glycosylation of glycoproteins, with down-regulation of cell cycle and DNA stability and chromosome maintenance genes. In contrast, Oropouche virus up-regulated cell cycle and DNA maintenance and condensation pathways while down-regulated extracellular matrix, collagen metabolism, glutamate and ion transporters pathways. Zika virus infection only up-regulated eukaryotic translation machinery while down-regulated interferon pathways. Reactome and integration analysis revealed a common signature in down-regulation of innate immune response, antiviral response, and inflammatory cytokines associated to interferon pathway for all arboviruses tested. Validation of interferon stimulated genes by reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) corroborated our transcriptome findings. Altogether, our results showed a co-evolution in the mechanisms involved in the escape of arboviruses to antiviral immune response mediated by the interferon (IFN) pathway.

Chikungunya virus (CHIKV) is a re-emergent arbovirus that causes a disease characterized primarily by fever, rash and severe persistent polyarthralgia, although

The symbiosis in trypanosomatids is a mutualistic relationship characterized by extensive metabolic exchanges between the bacterium and the protozoan. The symbiotic bacterium can complete host essential metabolic pathways, such as those for heme, amino acid, and vitamin production. Experimental assays indicate that the symbiont acquires phospholipids from the host trypanosomatid, especially phosphatidylcholine, which is often present in bacteria that have a close association with eukaryotic cells. In this work, an in-silico study was performed to find genes involved in the glycerophospholipid (GPL) production of Symbiont Harboring Trypanosomatids (SHTs) and their respective bacteria, also extending the search for trypanosomatids that naturally do not have symbionts. Results showed that most genes for GPL synthesis are only present in the SHT. The bacterium has an exclusive sequence related to phosphatidylglycerol production and contains genes for phosphatidic acid production, which may enhance SHT phosphatidic acid production. Phylogenetic data did not indicate gene transfers from the bacterium to the SHT nucleus, proposing that enzymes participating in GPL route have eukaryotic characteristics. Taken together, our data indicate that, differently from other metabolic pathways described so far, the symbiont contributes little to the production of GPLs and acquires most of these molecules from the SHT.

Abstract Background Bacterial and Archaeal communities have a complex, symbiotic role in crude oil bioremediation. Their biosurfactants and degradation enzymes have been in the spotlight, mainly due to the awareness of ecosystem pollution caused by crude oil accidents and their use. Initially, the scientific community studied the role of individual microbial species by characterizing and optimizing their biosurfactant and oil degradation genes, studying their individual distribution. However, with the advances in genomics, in particular with the use of New-Generation-Sequencing and Metagenomics, it is now possible to have a macro view of the complex pathways related to the symbiotic degradation of hydrocarbons and surfactant production. It is now possible, although more challenging, to obtain the DNA information of an entire microbial community before automatically characterizing it. By characterizing and understanding the interconnected role of microorganisms and the role of degradation and biosurfactant genes in an ecosystem, it becomes possible to develop new biotechnological approaches for bioremediation use. This paper analyzes 46 different metagenome samples, spanning 20 biomes from different geographies obtained from different research projects. Results A metagenomics bioinformatics pipeline, focused on the biodegradation and biosurfactant-production pathways, genes and organisms, was applied. Our main results show that: (1) surfactation and degradation are correlated events, and therefore should be studied together; (2) terrestrial biomes present more degradation genes, especially cyclic compounds, and less surfactation genes, when compared to water biomes; and (3) latitude has a significant influence on the diversity of genes involved in biodegradation and biosurfactant production. This suggests that microbiomes found near the equator are richer in genes that have a role in these processes and thus have a higher biotechnological potential. Conclusion In this work we have focused on the biogeographical distribution of hydrocarbon degrading and biosurfactant producing genes. Our principle results can be seen as an important step forward in the application of bioremediation techniques, by considering the biostimulation, optimization or manipulation of a starting microbial consortia from the areas with higher degradation and biosurfactant producing genetic diversity.

RASopathies are a group of rare genetic diseases caused by germline mutations in genes involved in the RAS–mitogen-activated protein kinase (RAS-MAPK) pathway. Whole-exome sequencing (WES) is a powerful approach for identifying new variants in coding and noncoding DNA sequences, including miRNAs. miRNAs are fine-tuning negative regulators of gene expression. The presence of variants in miRNAs could lead to malfunctions of regulation, resulting in diseases. Here, we identified 41 variants in mature miRNAs through WES analysis in five patients with previous clinical diagnosis of RASopathies syndromes. The pathways, biological processes, and diseases that were over-represented among the target genes of the mature miRNAs harboring variants included the RAS, MAPK, RAP1, and PIK3-Akt signaling pathways, neuronal differentiation, neurogenesis and nervous system development, congenital cardiac defects (hypertrophic cardiomyopathy, dilated cardiomyopathy, and arrhythmogenic right ventricular cardiomyopathy), and the phenotypes and syndromes of RASopathies (Noonan syndrome, Legius syndrome, Costello syndrome, Cafe au lait spots multiple, subaortic stenosis, pulmonary valve stenosis, and LEOPARD syndrome). Furthermore, eight selected variants in nine mature miRNAs (hsa-miR-1304, hsa-miR-146a, hsa-miR-196a2, hsa-miR-499a/hsa-miR-499b, hsa-miR-449b, hsa-miR-548l, hsa-miR-575, and hsa-miR-593) may have caused alterations in the secondary structures of miRNA precursor. Selected miRNAs containing variants such as hsa-miR-146a-3p, hsa-miR-196a-3p, hsa-miR-548l, hsa-miR-449b-5p, hsa-miR-575, and hsa-miR499a-3p could regulate classical genes associated with Rasopathies and RAS-MAPK pathways, contributing to modify the expression pattern of miRNAs in patients. RT-qPCR expression analysis revealed four differentially expressed miRNAs that were downregulated: miRNA-146a-3p in P1, P2, P3, P4, and P5, miR-1304-3p in P2, P3, P4, and P5, miR-196a2-3p in P3, and miR-499b-5p in P1. miR-499a-3p was upregulated in P1, P3, and P5. These results indicate that miRNAs show different expression patterns when these variants are present in patients. Therefore, this study characterized the role of miRNAs harboring variants related to RASopathies for the first time and indicated the possible implications of these variants for phenotypes of RASopathies such as congenital cardiac defects and cardio-cerebrovascular diseases. The expression and existence of miRNA variants may be used in the study of biomarkers of the RASopathies.

Abstract Background The emergence of multidrug-resistant Klebsiella pneumoniae is a major public health concern. Many K. pneumoniae infections can only be treated when resorting to last-line drugs such as polymyxin B (PB). However, resistance to this antibiotic is also observed, although insufficient information is described on its mode of action as well as the mechanisms used by resistant bacteria to evade its effects. We aimed to study PB resistance and the influence of abiotic stresses in a clinical K. pneumoniae strain using whole transcriptome profiling. Results We sequenced 12 cDNA libraries of K. pneumoniae Kp13 bacteria, from two biological replicates of the original strain Kp13 (Kp13) and five derivative strains: induced high-level PB resistance in acidic pH (Kp13pH), magnesium deprivation (Kp13Mg), high concentrations of calcium (Kp13Ca) and iron (Kp13Fe), and a control condition with PB (Kp13PolB). Our results show the involvement of multiple regulatory loci that differentially respond to each condition as well as a shared gene expression response elicited by PB treatment, and indicate the participation of two-regulatory components such as ArcA-ArcB, which could be involved in re-routing the K. pneumoniae metabolism following PB treatment. Modules of co-expressed genes could be determined, which correlated to growth in acid stress and PB exposure. We hypothesize that polymyxin B induces metabolic shifts in K. pneumoniae that could relate to surviving against the action of this antibiotic. Conclusions We obtained whole transcriptome data for K. pneumoniae under different environmental conditions and PB treatment. Our results supports the notion that the K. pneumoniae response to PB exposure goes beyond damaged membrane reconstruction and involves recruitment of multiple gene modules and intracellular targets.

Abstract Background The described species from the Metarhizium genus are cosmopolitan fungi that infect arthropod hosts. Interestingly, while some species infect a wide range of hosts (host-generalists), other species infect only a few arthropods (host-specialists). This singular evolutionary trait permits unique comparisons to determine how pathogens and virulence determinants emerge. Among the several virulence determinants that have been described, secondary metabolites (SMs) are suggested to play essential roles during fungal infection. Despite progress in the study of pathogen-host relationships, the majority of genes related to SM production in Metarhizium spp. are uncharacterized, and little is known about their genomic organization, expression and regulation. To better understand how infection conditions may affect SM production in Metarhizium anisopliae, we have performed a deep survey and description of SM biosynthetic gene clusters (BGCs) in M. anisopliae, analyzed RNA-seq data from fungi grown on cattle-tick cuticles, evaluated the differential expression of BGCs, and assessed conservation among the Metarhizium genus. Furthermore, our analysis extended to the construction of a phylogeny for the following three BGCs: a tropolone/citrinin-related compound (MaPKS1), a pseurotin-related compound (MaNRPS-PKS2), and a putative helvolic acid (MaTERP1). Results Among 73 BGCs identified in M. anisopliae, 20 % were up-regulated during initial tick cuticle infection and presumably possess virulence-related roles. These up-regulated BGCs include known clusters, such as destruxin, NG39x and ferricrocin, together with putative helvolic acid and, pseurotin and tropolone/citrinin-related compound clusters as well as uncharacterized clusters. Furthermore, several previously characterized and putative BGCs were silent or down-regulated in initial infection conditions, indicating minor participation over the course of infection. Interestingly, several up-regulated BGCs were not conserved in host-specialist species from the Metarhizium genus, indicating differences in the metabolic strategies employed by generalist and specialist species to overcome and kill their host. These differences in metabolic potential may have been partially shaped by horizontal gene transfer (HGT) events, as our phylogenetic analysis provided evidence that the putative helvolic acid cluster in Metarhizium spp. originated from an HGT event. Conclusions Several unknown BGCs are described, and aspects of their organization, regulation and origin are discussed, providing further support for the impact of SM on the Metarhizium genus lifestyle and infection process.

Abstract Background Bacterial non-coding RNAs act by base-pairing as regulatory elements in crucial biological processes. We performed the identification of trans-encoded small RNAs (sRNA) from the genomes of Mycoplama hyopneumoniae, Mycoplasma flocculare and Mycoplasma hyorhinis, which are Mycoplasma species that have been identified in the porcine respiratory system. Results A total of 47, 15 and 11 putative sRNAs were predicted in M. hyopneumoniae, M. flocculare and M. hyorhinis, respectively. A comparative genomic analysis revealed the presence of species or lineage specific sRNA candidates. Furthermore, the expression profile of some M. hyopneumoniae sRNAs was determined by a reverse transcription amplification approach, in three different culture conditions. All tested sRNAs were transcribed in at least one condition. A detailed investigation revealed a differential expression profile for two M. hyopneumoniae sRNAs in response to oxidative and heat shock stress conditions, suggesting that their expression is influenced by environmental signals. Moreover, we analyzed sRNA-mRNA hybrids and accessed putative target genes for the novel sRNA candidates. The majority of the sRNAs showed interaction with multiple target genes, some of which could be linked to pathogenesis and cell homeostasis activity. Conclusion This study contributes to our knowledge of Mycoplasma sRNAs and their response to environmental changes. Furthermore, the mRNA target prediction provides a perspective for the characterization and comprehension of the function of the sRNA regulatory mechanisms.

Restriction and Modification (R-M) systems are present in all Mycoplasma species sequenced so far. The presence of these genes poses barriers to gene transfer and could protect the cell against phage infections. The number and types of R-M genes between different Mycoplasma species are variable, which is characteristic of a polymorphism. The majority of the CDSs code for Type III R-M systems and particularly for methyltransferase enzymes, which suggests that functions other than the protection against the invasion of heterologous DNA may exist. A possible function of these enzymes could be the protection against the invasion of other but similar R-M systems. In Mycoplasma hyopneumoniae strain J, three of the putative methyltransferase genes were clustered in a region forming a genomic island. Many R-M CDSs were mapped in the vicinity of transposable elements suggesting an association between these genes and reinforcing the idea of R-M systems as mobile selfish DNA. Also, many R-M genes present repeats within their coding sequences, indicating that their expression is under the control of phase variation mechanisms. Altogether, these data suggest that R-M systems are a remarkable characteristic of Mycoplasma species and are probably involved in the adaptation of these bacteria to different environmental conditions.

Bacterial cell division has been studied mainly in model systems such as Escherichia coli and Bacillus subtilis, where it is described as a complex process with the participation of a group of proteins which assemble into a multiprotein complex called the septal ring. Mycoplasmas are cell wall-less bacteria presenting a reduced genome. Thus, it was important to compare their genomes to analyze putative genes involved in cell division processes. The division and cell wall (dcw) cluster, which in E. coli and B. subtilis is composed of 16 and 17 genes, respectively, is represented by only three to four genes in mycoplasmas. Even the most conserved protein, FtsZ, is not present in all mycoplasma genomes analyzed so far. A model for the FtsZ protein from Mycoplasma hyopneumoniae and Mycoplasma synoviae has been constructed. The conserved residues, essential for GTP/GDP binding, are present in FtsZ from both species. A strong conservation of hydrophobic amino acid patterns is observed, and is probably necessary for the structural stability of the protein when active. M. synoviae FtsZ presents an extended amino acid sequence at the C-terminal portion of the protein, which may participate in interactions with other still unknown proteins crucial for the cell division process.

Mollicutes are cell wall-less bacteria with a genome characterized by its small size. Chromosomal rearrangements help these organisms evade host immune surveillance and hence cause disease. Our goal was to determine genes shared by Mollicutes genomes using the bidirectional best hit methodology. The twelve studied Mollicutes share 210 genes, most of which (> 60%) fall into the following COG categories: translation, ribosomal structure and biogenesis; DNA replication, recombination and repair; nucleotide transport and metabolism and energy production and conversion. Thirty Mollicute-specific genes were identified, 22 of them previously described as essential genes in Mycoplasma genitalium.

Trypanosoma rangeli is a hemoflagellate protozoan parasite infecting humans and other wild and domestic mammals across Central and South America. It does not cause human disease, but it can be mistaken for the etiologic agent of Chagas disease, Trypanosoma cruzi. We have sequenced the T. rangeli genome to provide new tools for elucidating the distinct and intriguing biology of this species and the key pathways related to interaction with its arthropod and mammalian hosts.The T. rangeli haploid genome is ∼ 24 Mb in length, and is the smallest and least repetitive trypanosomatid genome sequenced thus far. This parasite genome has shorter subtelomeric sequences compared to those of T. cruzi and T. brucei; displays intraspecific karyotype variability and lacks minichromosomes. Of the predicted 7,613 protein coding sequences, functional annotations could be determined for 2,415, while 5,043 are hypothetical proteins, some with evidence of protein expression. 7,101 genes (93%) are shared with other trypanosomatids that infect humans. An ortholog of the dcl2 gene involved in the T. brucei RNAi pathway was found in T. rangeli, but the RNAi machinery is non-functional since the other genes in this pathway are pseudogenized. T. rangeli is highly susceptible to oxidative stress, a phenotype that may be explained by a smaller number of anti-oxidant defense enzymes and heat-shock proteins.Phylogenetic comparison of nuclear and mitochondrial genes indicates that T. rangeli and T. cruzi are equidistant from T. brucei. In addition to revealing new aspects of trypanosome co-evolution within the vertebrate and invertebrate hosts, comparative genomic analysis with pathogenic trypanosomatids provides valuable new information that can be further explored with the aim of developing better diagnostic tools and/or therapeutic targets.

BACKGROUND: The genus Paracoccidioides comprises human thermal dimorphic fungi, which cause paracoccidioidomycosis (PCM), an important mycosis in Latin America. Adaptation to environmental conditions is key to fungal survival during human host infection. The adaptability of carbon metabolism is a vital fitness attribute during pathogenesis. METHODOLOGY/PRINCIPAL FINDINGS: The fungal pathogen Paracoccidioides spp. is exposed to numerous adverse conditions, such as nutrient deprivation, in the human host. In this study, a comprehensive response of Paracoccidioides, Pb01, under carbon starvation was investigated using high-resolution transcriptomic (RNAseq) and proteomic (NanoUPLC-MSE) approaches. A total of 1,063 transcripts and 421 proteins were differentially regulated, providing a global view of metabolic reprogramming during carbon starvation. The main changes were those related to cells shifting to gluconeogenesis and ethanol production, supported by the degradation of amino acids and fatty acids and by the modulation of the glyoxylate and tricarboxylic cycles. This proposed carbon flow hypothesis was supported by gene and protein expression profiles assessed using qRT-PCR and western blot analysis, respectively, as well as using enzymatic, cell dry weight and fungus-macrophage interaction assays. The carbon source provides a survival advantage to Paracoccidioides inside macrophages. CONCLUSIONS/SIGNIFICANCE: For a complete understanding of the physiological processes in an organism, the integration of approaches addressing different levels of regulation is important. To the best of our knowledge, this report presents the first description of the responses of Paracoccidioides spp. to host-like conditions using large-scale expression approaches. The alternative metabolic pathways that could be adopted by the organism during carbon starvation can be important for a better understanding of the fungal adaptation to the host, because systems for detecting and responding to carbon sources play a major role in adaptation and persistence in the host niche.

Although magnetotactic bacteria (MTB) are ubiquitous in aquatic habitats, they are still considered fastidious microorganisms with regard to growth and cultivation with only a relatively low number of axenic cultures available to date. Here, we report the first axenic culture of an MTB isolated in the Southern Hemisphere (Itaipu Lagoon in Rio de Janeiro, Brazil). Cells of this new isolate are coccoid to ovoid in morphology and grow microaerophilically in semi-solid medium containing an oxygen concentration ([O2]) gradient either under chemoorganoheterotrophic or chemolithoautotrophic conditions. Each cell contains a single chain of approximately 10 elongated cuboctahedral magnetite (Fe3O4) magnetosomes. Phylogenetic analysis based on the 16S rRNA gene sequence shows that the coccoid MTB isolated in this study represents a new genus in the Alphaproteobacteria; the name Magnetofaba australis strain IT-1 is proposed. Preliminary genomic data obtained by pyrosequencing shows that Mf. australis strain IT-1 contains a genomic region with genes involved in biomineralization similar to those found in the most closely related magnetotactic cocci Magnetococcus marinus strain MC-1. However, organization of the magnetosome genes differs from Mc. marinus.

Fungal chitin metabolism involves diverse processes such as metabolically active cell wall maintenance, basic nutrition, and different aspects of virulence. Chitinases are enzymes belonging to the glycoside hydrolase family 18 (GH18) and 19 (GH19) and are responsible for the hydrolysis of β-1,4-linkages in chitin. This linear homopolymer of N-acetyl-β-D-glucosamine is an essential constituent of fungal cell walls and arthropod exoskeletons. Several chitinases have been directly implicated in structural, morphogenetic, autolytic and nutritional activities of fungal cells. In the entomopathogen Metarhizium anisopliae, chitinases are also involved in virulence. Filamentous fungi genomes exhibit a higher number of chitinase-coding genes than bacteria or yeasts. The survey performed in the M. anisopliae genome has successfully identified 24 genes belonging to glycoside hydrolase family 18, including three previously experimentally determined chitinase-coding genes named chit1, chi2 and chi3. These putative chitinases were classified based on domain organization and phylogenetic analysis into the previously described A, B and C chitinase subgroups, and into a new subgroup D. Moreover, three GH18 proteins could be classified as putative endo-N-acetyl-β-D-glucosaminidases, enzymes that are associated with deglycosylation and were therefore assigned to a new subgroup E. The transcriptional profile of the GH18 genes was evaluated by qPCR with RNA extracted from eight culture conditions, representing different stages of development or different nutritional states. The transcripts from the GH18 genes were detected in at least one of the different M. anisopliae developmental stages, thus validating the proposed genes. Moreover, not all members from the same chitinase subgroup presented equal patterns of transcript expression under the eight distinct conditions studied. The determination of M. anisopliae chitinases and ENGases and a more detailed study concerning the enzymes' roles in morphological or nutritional functions will allow comprehensive insights into the chitinolytic potential of this highly infective entomopathogenic fungus.

Endosymbiont-bearing trypanosomatids have been considered excellent models for the study of cell evolution because the host protozoan co-evolves with an intracellular bacterium in a mutualistic relationship. Such protozoa inhabit a single invertebrate host during their entire life cycle and exhibit special characteristics that group them in a particular phylogenetic cluster of the Trypanosomatidae family, thus classified as monoxenics. In an effort to better understand such symbiotic association, we used DNA pyrosequencing and a reference-guided assembly to generate reads that predicted 16,960 and 12,162 open reading frames (ORFs) in two symbiont-bearing trypanosomatids, Angomonas deanei (previously named as Crithidia deanei) and Strigomonas culicis (first known as Blastocrithidia culicis), respectively. Identification of each ORF was based primarily on TriTrypDB using tblastn, and each ORF was confirmed by employing getorf from EMBOSS and Newbler 2.6 when necessary. The monoxenic organisms revealed conserved housekeeping functions when compared to other trypanosomatids, especially compared with Leishmania major. However, major differences were found in ORFs corresponding to the cytoskeleton, the kinetoplast, and the paraflagellar structure. The monoxenic organisms also contain a large number of genes for cytosolic calpain-like and surface gp63 metalloproteases and a reduced number of compartmentalized cysteine proteases in comparison to other TriTryp organisms, reflecting adaptations to the presence of the symbiont. The assembled bacterial endosymbiont sequences exhibit a high A+T content with a total of 787 and 769 ORFs for the Angomonas deanei and Strigomonas culicis endosymbionts, respectively, and indicate that these organisms hold a common ancestor related to the Alcaligenaceae family. Importantly, both symbionts contain enzymes that complement essential host cell biosynthetic pathways, such as those for amino acid, lipid and purine/pyrimidine metabolism. These findings increase our understanding of the intricate symbiotic relationship between the bacterium and the trypanosomatid host and provide clues to better understand eukaryotic cell evolution.