Snake venoms are biological weapon systems composed of secreted proteins and peptides that are used for immobilizing or killing prey. Although post-translational modifications are widely investigated because of their importance in many biological phenomena, we currently still have little understanding of how protein glycosylation impacts the variation and stability of venom proteomes. To address these issues, here we characterized the venom proteomes of seven Bothrops snakes using a shotgun proteomics strategy. Moreover, we compared the electrophoretic profiles of native and deglycosylated venoms and, in order to assess their subproteomes of glycoproteins, we identified the proteins with affinity for three lectins with different saccharide specificities and their putative glycosylation sites. As proteinases are abundant glycosylated toxins, we examined the effect of N-deglycosylation on their catalytic activities and show that the proteinases of the seven venoms were similarly affected by removal of N-glycans. Moreover, we prospected putative glycosylation sites of transcripts of a B. jararaca venom gland data set and detected toxin family related patterns of glycosylation. Based on our global analysis, we report that Bothrops venom proteomes and glycoproteomes contain a core of components that markedly define their composition, which is conserved upon evolution in parallel to other molecular markers that determine their phylogenetic classification.

Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome of Bothrops jararaca, a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A2 (PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option.

Interspecific differences in snake venom compositions can result from distinct regulatory mechanisms acting in each species. However, comparative analyses focusing on identifying regulatory elements and patterns that led to distinct venom composition are still scarce. Among venomous snakes, Bothrops cotiara and Bothrops fonsecai represent ideal models to complement our understanding of the regulatory mechanisms of venom production. These recently diverged species share a similar specialized diet, habitat, and natural history, but each presents a distinct venom phenotype. Here, we integrated data from the venom gland transcriptome and miRNome and the venom proteome of B. fonsecai and B. cotiara to better understand the regulatory mechanisms that may be acting to produce differing venom compositions. We detected not only the presence of similar toxin isoforms in both species but also distinct expression profiles of phospholipases A2 (PLA2) and some snake venom metalloproteinases (SVMPs) and snake venom serine proteinases (SVSPs) isoforms. We found evidence of modular expression regulation of several toxin isoforms implicated in venom divergence and observed correlated expression of several transcription factors. We did not find strong evidence for miRNAs shaping interspecific divergence of the venom phenotypes, but we identified a subset of toxin isoforms whose final expression may be fine-tuned by specific miRNAs. Sequence analysis on orthologous toxins showed a high rate of substitutions between PLA2s, which indicates that these toxins may be under strong positive selection or represent paralogous toxins in these species. Our results support other recent studies in suggesting that gene regulation is a principal mode of venom evolution across recent timescales, especially among species with conserved ecotypes.

DNA methylation is one of the epigenetic modifications that configures gene transcription programs. This study describes the DNA methylation profile of HIV-infected individuals with distinct characteristics related to natural and artificial viremia control. Sheared DNA from circulating mononuclear cells was subjected to target enrichment bisulfite sequencing designed to cover CpG-rich genomic regions. Gene expression was assessed through RNA-seq. Hypermethylation in virologic responders was highly distributed closer to Transcription Start Sites (p-value = 0.03). Hyper and hypomethylation levels within TSS adjacencies varied according to disease progression status (Kruskal–Wallis, p < 0.001), and specific differentially methylated regions associated genes were identified for each group. The lower the promoter methylation, the higher the gene expression in subjects undergoing virologic failure (R = − 0.82, p = 0.00068). Among the inversely correlated genes, those supporting glycolysis and its related pathways were hypomethylated and up-regulated in virologic failures. Disease progression heterogeneity was associated with distinct DNA methylation patterns in terms of rates and distribution. Methylation was associated with the expression of genes sustaining intracellular glucose metabolism in subjects undergoing antiretroviral virologic failure. Our findings highlight that DNA methylation is associated with latency, disease progression, and fundamental cellular processes.

No abstract