Next Article in Journal
Fine Wrinkle Improvement through Bioactive Materials That Modulate EDAR and BNC2 Gene Expression
Previous Article in Journal
Extracellular Vesicles as Delivery Vehicles for Non-Coding RNAs: Potential Biomarkers for Chronic Liver Diseases
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Complex Pattern of Gene Expression in Tissue Affected by Viperid Snake Envenoming: The Emerging Role of Autophagy-Related Genes

by
Ana Karina de Oliveira
1,
Alexandra Rucavado
2,
Teresa Escalante
2,
José María Gutiérrez
2 and
Jay W. Fox
3,*
1
Department of Pathology, University of Virginia School of Medicine, Charlottesville, VA 22908, USA
2
Instituto Clodomiro Picado, Facultad de Microbiología, Universidad de Costa Rica, San José 11501, Costa Rica
3
Department of Microbiology, Immunology and Cancer Biology, University of Virginia School of Medicine, Charlottesville, VA 22908, USA
*
Author to whom correspondence should be addressed.
Biomolecules 2024, 14(3), 278; https://doi.org/10.3390/biom14030278
Submission received: 10 January 2024 / Revised: 22 February 2024 / Accepted: 23 February 2024 / Published: 26 February 2024
(This article belongs to the Section Molecular Biomarkers)

Abstract

:
Viperid snake venoms induce severe tissue damage, characterized by the direct toxic action of venom components, i.e., phospholipases A2 (PLA2s) and metalloproteinases (SVMPs), concomitantly with the onset of endogenous inflammatory processes, in an intricate scenario of tissue alterations. Understanding the expression of relevant genes in muscle tissue will provide valuable insights into the undergoing pathological and inflammatory processes. In this study, we have used the Nanostring technology to evaluate the patterns of gene expression in mouse skeletal muscle 1 h, 6 h, and 24 h after injection of the venoms of Bothrops asper and Daboia russelii, two medically relevant species in Latin America and Asia, respectively, with somewhat different clinical manifestations. The dose of venoms injected (30 µg) induced local pathological effects and inflammation in muscle tissue. We focused our analysis on genes related to extracellular matrix (ECM) metabolism, immune system, programmed cell death, and autophagy. The results revealed a complex pattern of expression of genes. Regarding ECM metabolism and regulation, up-regulated genes included proteinase inhibitor Serpine 1, thrombospondin 1, collagens 1A1 and 4A1 (at 1 h in the case of B. asper), TIMP1, MMP-3 (at 24 h), and lysil oxidase (LOX). In contrast, collagen chains 5A3 and 5A1 were down-regulated, especially at 6 h. Transforming growth factor β (TGF-β) and several genes related to myofibroblast regulation were also up-regulated, which might be related to the development of fibrosis. Several genes related to cytokine and chemokine synthesis and regulation and NFκB signaling were also up-regulated. Our observations show a variable expression of genes associated with programmed cell death and autophagy, thus revealing a hitherto unknown role of autophagy in tissue affected by snake venoms. These results provide clues to understanding the complex pattern of gene expression in tissue affected by viperid snake venoms, which likely impacts the final pathophysiology of damaged tissue in envenomings.

1. Introduction

Envenomings by viperid snakes are characterized by severe pathological effects at the anatomical site of venom injection, associated with tissue necrosis, hemorrhage, blistering, and edema [1,2]. Such effects are primarily caused by the direct action of tissue-damaging toxins, mostly phospholipases A2 (PLA2s) and metalloproteinases (SVMPs) [3,4,5]. Concomitantly, a complex inflammatory response develops in the affected tissues, characterized by the recruitment of inflammatory cells, i.e., neutrophils and macrophages, the activation of resident cells, and the synthesis of a plethora of inflammatory mediators and damage-associated molecular patterns (DAMPs) from damaged cells and extracellular matrix [6,7]. Such endogenous response contributes to edema and clearance of damaged tissue and sets the stage for reparative and regenerative processes. At the same time, local inflammation might potentiate some of the deleterious effects caused by the venom.
Understanding this altered tissue environment has been challenging because of its complexity, its dynamics over time, and its spatial heterogeneity in various regions of the affected tissue [8]. Several studies have quantified specific inflammatory mediators in tissues or blood in experimental models of viperid envenomings and have contributed to our understanding of these phenomena (see reviews by Teixeira et al. [6,9]). However, to fully grasp the complexity of this inflammatory scenario it is necessary to introduce novel analytical tools that allow the analysis of a vast number of tissue components and inflammatory mediators. This will provide novel clues on the genesis and dynamics of inflammatory and tissue reparative and regenerative processes. To this end, the Nanostring nCounter platform is one approach to provide novel clues by allowing the analysis of the expression of up to 800 genes simultaneously, without amplification, cDNA, or library preparation [10,11]. This technology uses direct, digital quantitation of mRNA transcripts by hybridization to color-barcode-specific probes in a high-sensitivity manner [10]. The ability of probes to hybridize a complementary region of 100 nucleotides long allows this technology to measure mRNA with low abundance and sub-optimal RNA quality. In this context, analyzing RNA extracted from formalin-fixed paraffin-embedded (FFPE) tissue samples allows us to investigate the correlation between mRNA expression and the histological alterations occurring in the tissue (spatial omics).
Over the past years, our group has focused our studies on several critical aspects of the pathological and inflammatory effects induced by the venoms of two viperid species, Daboia russelii and Bothrops asper, which are medically relevant snakes in Asia and Latin America, respectively, owing to the high number of cases they inflict and to the severity of clinical manifestations. These venoms inflict envenomings with some similar clinical features, including local tissue damage, hemorrhage, coagulopathy, and hemodynamic disturbances [12,13]. At the same time, there are notable differences in their actions since D. russelii envenomings are also characterized by systemic capillary leakage syndrome, high incidence of acute kidney injury, and, in some geographical locations, neurotoxicity [12], effects, which are largely absent in B. asper envenomings. Experimentally, these two venoms induce a different pattern of vasculotoxicity since B. asper venom causes strong hemorrhage in rodent models, while D. russelii venom induces less extensive hemorrhage but a conspicuous increase in vascular permeability, associated with hemoconcentration [14,15]. Furthermore, they show variations in the patterns of increments in some inflammatory mediators [16]. Lipidomic and metabolomic analyses of plasma from mice injected with these venoms demonstrated both similarities and differences [17]. However, there is a notorious lack of information on the patterns of expression of genes that might be involved in the events occurring in tissue injected with venoms.
To further explore and compare the effects of these two venoms in the muscle tissue of mice, in this study, we have employed the Nanostring nCounter technology using the mouse Fibrosis V2 panel to follow the mRNA expression of several markers associated with extracellular matrix synthesis and degradation, innate immune system changes, and cell death pathways. Our goal was to delineate the patterns of gene expression in envenomed tissue, which may explain some of the pathological effects of envenoming and reveal novel mediators and processes involved in tissue damage and repair. Our findings highlight similarities and differences between the action of these venoms and identify trends in the time course of the expression of these markers. They also allow the identification of tissue and inflammatory markers, which further our understanding of venom-induced pathology and inflammation.

2. Materials and Methods

2.1. Venoms

The venom of D. russelii was purchased from Latoxan (Code L1132A; Lot: 015.051; Portes-lès-Valence, France). It is a pool of venom collected from several adult specimens from Pakistan. The venom of B. asper was provided by Instituto Clodomiro Picado, University of Costa Rica. It is a pool of venom from more than 20 adult specimens collected in the Pacific region of Costa Rica and kept at the Serpentarium of Instituto Clodomiro Picado. Venoms were freeze-dried and stored at −20 °C. They were dissolved in 0.12 M NaCl, 0.04 M phosphates, pH 7.2 (PBS) immediately before use.

2.2. Experimental Protocol

Groups of three mice of both sexes (CD-1 strain, 18–20 g) received an intramuscular injection in the right gastrocnemius of 30 µg venom of each species, dissolved in 50 µL PBS. Control mice were injected with 50 µL PBS only. At the time intervals of 1 h, 6 h, and 24 h following injection, groups of three mice were sacrificed by CO2 inhalation, and a sample of the injected gastrocnemius muscle was excised and added to 10% formalin solution in water. After 48 h fixation, routine processing of tissues was performed, followed by embedding in paraffin. Muscle tissue sections (4 µm thickness) were stained with hematoxylin and eosin or Masson’s trichrome for histological observation of tissue damage and inflammatory infiltrate. The protocols of experiments using mice were approved by the Institutional Committee for the Care and Use of Laboratory Animals (CICUA) of the University of Costa Rica (approval number CICUA-032-2020) and meet the International Guiding Principles for Biomedical Research Involving Animals (CIOMS). Mice were maintained in Tecniplast Eurostandard Type II 1264C cages (L 25 × W 40 × H 14 cm), six mice per cage. Animals were kept at 18–24 °C, 60–65% relative humidity, and a 12:12 h light–dark cycle.

2.3. RNA Extraction

A total of 21 blocks of formalin-fixed paraffin-embedded muscle mouse tissue (Control, B. asper (BA) BA-1 h, BA-6 h, BA-24 h, D. russelii (DR) DR-1 h, DR-6 h, DR-24 h) were used for RNA extraction. Four slices of 10 µm thickness from each selected block were deparaffinized using D-Limonene (Sigma, Kawasaki, Japan), digested with proteinase K, and isolated using RNeasy® Plus Mini Kit QIAGEN, according to the manufacturer’s protocol. The final volume of extracted RNA was 14 µL. RNA concentration and purity were assessed using a NanoDrop instrument. Sample concentration was measured at 260 nm and 280 nm, and the ratio of optical density 260/280 and 260/230 were used to test for protein and phenol contamination, respectively.

2.4. NanoString nCounter Analysis

Comparative analysis between samples collected from mice injected with B. asper and D. russelii venoms was conducted by Nanostring nCounter multiplex analysis [10] using the Mouse Fibrosis V2 panel (Supplementary Table S1). Nanostring probes are made with target-specific sequences and tag-specific sequences at 5′ and 3′ tailing ends (Supplementary Table S1). The RNA samples (100 ng) were incubated for 16 h at 65 °C in a hybridization buffer containing the CodeSet (reporter and capture probes). Hybridized samples were processed using the Prep Station using high sensitivity protocol, 3 h per 12-sample cartridge. The Prep Station purifies the RNA/probe complexes and places them in a cartridge where they are immobilized and aligned for data collection. Data acquisition was carried out in the NanoString nCounter Digital Analyzer with the ‘Max’ Field of View (FOV) setting to 555 images per sample in a 6 h scan per cartridge. Raw counts were normalized using the positive controls, and target genes were normalized to the five internal reference genes (Supplementary Table S2).

2.5. Data Analysis

NanoString gene expression profiling was conducted using the mouse nCounter Fibrosis V2 Panel (NanoString Technologies, Seattle, WA, USA), which contains probes specific to 760 endogenous genes and 10 housekeeping genes (ACAD9, ARMH3, CNOT10, GUSB, MTMR14, NOL7, NUBBP1, PGK1, PPIA, RPLP0) described in Supplementary Table S1.
Data were analyzed using ROSALIND (https://rosalind.bio/, accessed on 1 September 2023), with a hyper-scale cloud developed by ROSALIND, Inc. (San Diego, CA, USA), accessed on 5 March 2023. The quality control of samples was verified using spike-ins of positive and negative probes and housekeeping expression levels (Supplementary Figure S1). Normalization, fold changes, and p-values were calculated using criteria provided by NanoString. ROSALIND follows the nCounter Advanced Analysis protocol of dividing counts within a lane by the geometric mean of the normalizer probes from the same lane. The geNorm algorithm was used to select normalizer probes in the Bioconductor package NormqPCR, removing candidate housekeepers with the least stable expression relative to other candidates (Supplementary Table S2).
Fold changes and p-values were calculated using the “Fast” method, as described in the nCounter Advanced Analysis 2.0 User Manual (MAN-10030-03). Differential gene expression analysis was calculated using the generalized linear model (GLM) that was developed by NanoString for analysis of count data, assuming a negative binomial distribution. Adjusted p-values are calculated using the Benjamini-Hochberg False Discovery Rate (FDR) methodology (Supplementary Table S3).

2.6. Enrichment Pathway Analysis

The enrichment pathways were analyzed using the Gene Set Analysis (GSA) module from NanoString. GSA summarizes the change in regulation within each defined gene set relative to the baseline, as described in the manufacturer’s manual (MAN-10030-03). The values calculated are the Global Significance Score (also called Undirected Global Significance Score), which measures the overall differential expression of the selected gene set relative to selected patient populations, ignoring whether each gene is up- or down-regulated (Supplementary Table S4).
Another variant calculated was the Directed Global Significance scores, which measure the extent to which a given gene set is up- or down-regulated relative to a given covariate. It is calculated similarly to the undirected global significance score, but it takes the sign of the t-statistics into account (Supplementary Table S4).

3. Results and Discussion

3.1. Histopathological Observations

The dose of venoms used in this study (30 µg) was selected on the basis of previous studies and corresponds to a dose that induces overt pathological and inflammatory effects in the muscle tissue of mice [14,16]. In agreement with these previous reports [14,16], histological analysis of muscle tissue injected with venoms of B. asper and D. russelii revealed a different pattern of local tissue damage. B. asper venom induced myonecrosis and hemorrhage, evidenced by the presence of hypercontracted muscle fibers and abundant erythrocytes in the interstitial space. In contrast, tissue from mice injected with D. russelii venom showed muscle fiber necrosis, but hemorrhage was largely absent. In both cases, an inflammatory infiltrate was observed in the damaged tissue (Figure 1)

3.2. ECM Synthesis, Degradation, and Modification

A variable pattern of expression of genes associated with ECM synthesis, degradation, and modification was observed at the three time intervals in samples from muscle tissue injected with B. asper and D. russelii venoms, with quantitative differences in the fold change expression between venoms and between time intervals (Figure 2 and Figure 3). Up-regulated genes, as compared to control samples, included the proteinase inhibitor Serpine 1, thrombospondin 1, collagens 1A1 and 4A1 (at 1 h in the case of B. asper), TIMP1, MMP-3 (at 24 h), lysil oxidase (LOX), and integrin subunit α5, among others. As a general trend, changes were more pronounced in samples collected from B. asper-injected mice, especially at 1 and 6 h (Figure 2 and Figure 3). In contrast, genes coding for MMP-2 and collagen 5A3 and 5A1 were down-regulated, especially at 6 h, compared to controls injected with PBS only (Figure 2 and Figure 3).
These findings reveal a variable pattern of gene regulation, with some genes coding for fibrillar (collagen I) and basement membrane (collagen IV) collagens being up-regulated, while chains of collagen V were down-regulated. Expression of MMPs also evidenced up- and down-regulation, whereas two proteinase inhibitors were up-regulated, indicating a process of modulation of ECM degradation. Fragments of several serpins and TIMPS have been described in exudates collected from mice injected with venom and tissue-damaging toxins from B. asper [18,19], implying a control over ECM degradation by these proteinase inhibitors. Our observations reveal a variable pattern of expression of genes coding for ECM-associated proteins within the previously described pattern of ECM degradation, synthesis, and modulation [18,19,20,21,22]. The consequences in the overall tissue damage of this differential expression of MMPs and ECM components remain to be investigated. Despite some qualitatively similar patterns observed in gene expression between the two venoms, tissue from B. asper-injected venom showed a higher quantitative expression of most of these proteins. This is likely related to the fact that B. asper venom has higher proteolytic and hemorrhagic activities than D. russelii venom, in agreement with the higher content of SVMPs in B. asper venom [23,24,25], thus causing more extensive ECM degradation and turnover. It is necessary to further investigate the actual degradation patterns of ECM proteins in these in vivo models, to assess which components are degraded at different time intervals. Our findings suggest a higher ECM turnover in muscle injected with B. asper venom.
Thrombospondin 1 is a counter-adhesive ECM protein, which is also found in platelets. It plays many roles and was detected by proteomics in the exudate from B. asper venom-injected mice [19]. Thrombospondin 1 is known to bind collagens [26] and to activate TGFβ [27], which is a pleiotropic cytokine that, among its many actions, promotes fibrosis [28]. The Nanostring analysis of the TGFβ pathway revealed up-regulation of the genes coding for TGFβ and TGFβ receptor at some time intervals, as well as LRRC 32, a key regulator of TGFβ (Figure 4 and Figure 5) [29]. The enzyme lysyl oxidase was also up-regulated. It is related to the metabolism of collagen and has been described in a variety of fibrotic diseases [30]. Taken together, these observations highlight a scenario of up-regulation of fibrosis-promoting genes.
We also investigated gene expression related to myofibroblast regulation since these cells play a key role in wound contraction but may also stimulate fibrosis in a variety of pathologies [31]. The genes coding for chemokine interferon-γ-inducible protein-10 (CXCL10), integrin subunit α5, and lysyl oxidase, and connective tissue growth factor (CTGF), all related to the myofibroblast regulation pathways, were up-regulated in tissue affected by the two venoms (Figure 2 and Figure 3), reflecting stimuli for myofibroblast proliferation and collagen production. Our findings agree with previous studies on muscle repair and regeneration after myonecrosis induced by B. asper venom showing a deficient regeneration associated with increased fibrosis, which is very likely a consequence of fibroblast proliferation and fibrillar collagen synthesis [32,33,34]. Thus, our data provide clues to understanding the imbalance between myoblast and fibroblast proliferation in muscle tissue affected by viperid venoms, which may shed light on possible ways to modulate this imbalance to favor muscle regeneration instead of fibrosis.

3.3. Innate Immune System

Gene expression related to the pathways of NF-κB, cytokine signaling, and Toll-like receptors were evaluated. A variable pattern of up-regulation was observed in many genes in tissue samples from mice injected with these venoms, although there was a higher expression of several genes associated with these pathways in the case of B. asper venom. The highest level of most of these transcripts in the case of B. asper was observed at 6 h (Figure 4 and Figure 5). Some up-regulated genes from these pathways with the two venoms included IL1β, IL1R1, CD14, Myd88, TGFβ, TGFβR, and the chemokines CCL2 and CXCL10. (Figure 4 and Figure 5).
Up-regulation of IL1β agrees with previous observations of elevated levels of this cytokine in muscle tissue and peritoneal fluid lavage from mice injected with several Bothrops sp venoms, including B. asper [35,36,37]. IL1β is a pleiotropic cytokine with a variety of proinflammatory effects, including neutrophil, T cell, B cell, and endothelial cell activation, monocyte differentiation, and TH17 response, by high expression of IL1β in macrophages and neutrophils inducing an increase in CD45 cells infiltration [38,39]. In chronic inflammation, IL1β is related to angiogenesis and immunosuppression of macrophages characterized by the expression of M2 macrophages [40]. Interestingly, IL1β is also involved in wound healing delay, promotes fibrosis, and impairs re-epithelialization [39,41].
CCL2 (monocyte chemoattractant protein-1, MCP-1) is a potent chemoattractant for monocytes/macrophages but also attracts other cells such as T cells, B cells, NK cells, and DC cells [42]. In addition to its chemoattractant effect, it impacts leukocyte function, induces autophagy, macrophage polarization according to the context of the microenvironment, and cytokine secretion, including TNF-α, IL1β, IL6, IL10, and IL12 [42,43,44]. CXCL10 mediates the activation and recruitment of leukocytes such as T cells, monocytes, and NK cells. It is secreted by several cells in response to IFN-γ. CXCL10 has pleiotropic effects, including inflammatory process, inhibition of neovascularization in early wound healing, and immune fibroblast activation [45,46]. Overall, the observed patterns of expression of cytokine and chemokine genes underscore an inflammatory landscape associated with leukocyte infiltration, activation of immune cells, synthesis of a plethora of mediators, and other inflammatory events, in agreement with the pro-inflammatory action of viperid snake venoms [6].
It is of little surprise that we identified the overexpression of myeloid cell markers CD14 and Myd88. The former is a membrane-associated protein that binds the complex of bacterial lipopolysaccharide (LPS) and LPS-binding protein and is expressed in monocytes, macrophages, and dendritic cells [47,48]. Myd88 is an adapter in immune cells that has a pivotal role in innate immunity through Toll-like receptors, acting as an adaptor molecule that relays signals from outside the cell to intracellular proteins [47,49].
Innate immunity is likely to play a role in the inflammation in tissues affected by snake venom toxins by detecting molecular patterns associated with venom components, i.e., venom-associated molecular patterns (VAMPs) [50], and with fragments of cells and extracellular matrix, i.e., damage-associated molecular patterns (DAMPs) released as a consequence of the tissue-damaging activity of venoms [7]. It was, therefore, of interest to assess the expression of genes related to TLR signaling. The observed up-regulation of TLR4 in the case of B. asper venom at 6 h (Figure 4) is relevant as this receptor has been associated with inflammation induced by exudates collected from tissue affected by B. asper venom [7] and by other venoms and toxins [51,52]. Thus, the up-regulation of genes related to TLR signaling may reflect an ongoing process of innate immune system response to the action of these venoms in muscle tissue. Noteworthy, there were differences in the expression of these genes between tissues injected with these venoms (Figure 5). Such differences might be related to previously described variations in the inflammatory response after injection of these two venoms in mice [16].

3.4. Programmed Cell Death and Autophagy

Genes associated with programmed cell death and autophagy were also analyzed (Figure 6 and Figure 7), as it is likely that these phenomena develop in tissue affected by snake venoms and toxins [34,53]. Bothrops sp. venom components have been shown to induce apoptosis and autophagy of various cell types in culture [54,55,56,57]. The genes BIRC3 and BCL2L1 were up-regulated in tissue injected with the two venoms. These genes are associated with programmed cell death. BIRC3 encodes a member of inhibitors of apoptosis (IAP) proteins that inhibit apoptosis by binding to TNFα-associated factors TRAF1 and TRAF2 and interfere with the activation of caspases [58]. BCL2L1 encodes for a protein of the BCL-2 family. These proteins are located in the mitochondrial outer membrane, regulating its permeabilization, and members of this family of proteins exert pro- and anti-apoptotic functions. Moreover, XIAP, an inhibitor of apoptosis [59], is down-regulated in the case of B. asper venom. In tissue injected with this venom, there is down-regulation of CASP7, a proapoptotic caspase (Figure 7). The observed pattern of expression of programmed cell death-related genes in our model calls for further studies aimed at understanding their role in venom-induced pathology.
Little is known about the role of autophagy in tissues affected by snake venoms. Autophagy plays a key role in tissue homeostasis and is regulated by an abundant set of genes [60]. In vitro, some snake toxins have been described to induce autophagy [56], but whether this phenomenon occurs in vivo in experimental envenoming is largely unknown. With some exceptions, most genes associated with autophagy were down-regulated with the two venoms, although a few were up-regulated, with quantitative differences between venoms (Figure 6 and Figure 7). Among the genes associated with autophagy and up-regulated in muscle injected with B. asper venom are CTSL and PPP2CA (Figure 7). The former encodes cathepsin L, a lysosomal cysteine proteinase that has been associated with several pathologies and is responsible for myofibrillar degradation in myopathies [61]. PPP2CA gene encodes the protein phosphatase 2 catalytic subunit and is implicated in the negative control of cell growth and division [62].
The gene LAMTOR2 is down-regulated in tissue injected with both venoms. It codes for a member of the Ragulator/LAMTOR complex known to regulate mTOR and has a role as a late endosomal/lysosomal adaptor [63,64]. Deficiencies in this gene affect innate and adaptive immunity, having an impact on endosomal and lysosomal biogenesis [63]. The gene IGF1R, a receptor of insulin-like growth factor, is up-regulated at 6 h in both venoms (Figure 7). Its inhibition is known to reduce autophagy [65], thus, its expression would stimulate this cell death process. The gene VAMP8 is down-regulated in tissue injected with B. asper venom. It codes for a SNARE protein located in the lysosomal membrane that plays a key role in the fusion of lysosomes with autophagosomes [66]. On the other hand, the genes ATG2B and IRS1 are down-regulated in tissue affected by both venoms, especially in the case of B. asper (Figure 7). The former is related to the autophagic flow [67]. Moreover, several genes related to the subunits of AMP-activated protein kinase, i.e., PRKAG1, PRKACA, PRKAB2, are down-regulated in tissue affected by the venoms. These genes are involved in responses to various forms of cellular stress and are positive regulators of autophagy [68]. Taken together, our findings reveal a trend toward an expression pattern that results in the inhibition of autophagy. Since this process is an adaptation of cells to various forms of stress, its inhibition may have implications for the overall development of venom-induced tissue damage and repair.

4. Conclusions

Our findings reveal a complex and intricate pattern of gene expression in muscle injected with B. asper and D. russelii venoms, with similarities and differences between the venoms and variations along the time course of the experiment. Rapid up-regulation and down-regulation of many genes occurred starting in the first hour after the onset of envenoming, underscoring a rapid response of muscle tissue to the venoms. We focused our analysis on genes related to ECM synthesis and degradation and myofibroblast regulation, as well as on genes related to the innate immune response, programmed cell death, and autophagy. Up-regulated genes included various genes related to collagen metabolism and myofibroblast proliferation, probably associated with the onset of fibrosis and tissue repair. A variety of genes related to the immune system were also up-regulated, suggesting activation of innate immune pathways as a consequence of tissue-damaging effects of the venoms. Moreover, alterations also occurred in genes related to programmed cell death and autophagy, underscoring that these processes may be involved in venom-induced tissue damage and repair. In particular, our findings on the expression of genes related to autophagy provide novel clues on this hitherto unexplored aspect of venom-induced pathology. These findings pave the way for future experiments in which drugs that modulate apoptosis and autophagy can be used to assess the effects of these cell death pathways in the overall outcome of venom-induced tissue damage. Further studies are needed to link these acute gene expression events with the development of medium- and long-term phenomena taking place in the tissues, associated with repair, regeneration, and chronic inflammation in snakebite envenoming.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/biom14030278/s1. Figure S1. Quality Control and Normalized Data. (A) Principal Component Analysis (PCA) showed a clear separation of venom toxins time point group can be seen in the PC1 vs. PC2 results, meaning that changes in this variable cause clear, consistent changes in the data. The analysis was made using https://www.bioinformatics.com.cn/srplot, accessed on 23 September 2023 (B–D) The quality control of samples was verified using spike-ins of positive and negative probes, and housekeeping expression levels. The graphics were generated using the Rosalind cloud platform. Table S1. Mouse Fibrosis Gene List and Probe Annotation. Table S2. Raw Counts and Log2 Normalized Counts. Data acquisition is carried out in the NanoString nCounter Digital Analyzer with the ‘Max’ Field of View (FOV) setting to 555 images per sample in a 6 h scan per cartridge. Raw counts were normalized using the positive controls, and target genes were normalized to the five internal reference genes. Table S3. Differential analysis—The gene expression was analyzed by comparing each time point of venom injection to the control group (PBS). Differential gene expression analysis was calculated using the generalized linear model (GLM) assuming a negative binomial distribution. Table S4. GSA Pathway Analysis. The enrichment pathways were analyzed using the Gene Set Analysis (GSA) module from NanoString on the Rosalind cloud platform.

Author Contributions

Conceptualization: A.K.d.O., A.R., T.E., J.M.G. and J.W.F. Methodology: A.K.d.O., A.R., T.E. and J.W.F. Resources: J.W.F. Data curation: A.K.d.O. Writing—original draft: A.K.d.O., J.M.G. and J.W.F. Writing—review and editing: A.K.d.O., A.R., T.E., J.M.G. and J.W.F. Project administration: J.W.F. Funding acquisition: A.R. and J.W.F. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Office of Research Core Administration, University of Virginia School of Medicine, Virginia, USA, and by Vicerrectoría de Investigación, Universidad de Costa Rica.

Institutional Review Board Statement

The protocols of experiments using mice were approved by the Institutional Committee for the Care and Use of Laboratory Animals (CICUA) of the University of Costa Rica (approval number CICUA-032-2020).

Data Availability Statement

Data are available in the GeoAccession number GSE248215 at: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE248215, accessed on 1 September 2023.

Acknowledgments

We acknowledge the use of the University of Virginia School of Medicine Spatial Biology Core Facility, RRID: SCR_023281; University of Virginia School of Medicine Research Histology Core Facility.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Warrell, D.A. Snake bite. Lancet 2010, 375, 77–88. [Google Scholar] [CrossRef]
  2. Gutiérrez, J.M.; Calvete, J.J.; Habib, A.G.; Harrison, R.A.; Williams, D.J.; Warrell, D.A. Snakebite envenoming. Nat. Rev. Dis. Prim. 2017, 3, 17063. [Google Scholar] [CrossRef]
  3. Gutiérrez, J.M.; Ownby, C.L. Skeletal muscle degeneration induced by venom phospholipases A 2: Insights into the mechanisms of local and systemic myotoxicity. Toxicon 2003, 42, 915–931. [Google Scholar] [CrossRef]
  4. Gutiérrez, J.; Rucavado, A.; Escalante, T. Snake venom metalloproteinases. Biological roles and participation in the pathophysiology of envenomation. In Handbook of Venoms and Toxins of Reptiles, 1st ed.; Mackessy, S.P., Ed.; CRC Press: Boca Raton, FL, USA, 2009; pp. 115–138. [Google Scholar]
  5. Fox, J.W.; Serrano, S.M.T. Structural considerations of the snake venom metalloproteinases, key members of the M12 reprolysin family of metalloproteinases. Toxicon 2005, 45, 969–985. [Google Scholar] [CrossRef]
  6. Teixeira, C.; Moreira, V.; Gutiérrez, J.M. Venoms. In Inflammation: From Molecular and Cellular Mechanisms to the Clinic; Wiley: Weinheim, Germany, 2017; pp. 99–128. [Google Scholar] [CrossRef]
  7. Rucavado, A.; Nicolau, C.A.; Escalante, T.; Kim, J.; Herrera, C.; Gutiérrez, J.M.; Fox, J.W. Viperid Envenomation Wound Exudate Contributes to Increased Vascular Permeability via a DAMPs/TLR-4 Mediated Pathway. Toxins 2016, 8, 349. [Google Scholar] [CrossRef]
  8. de Oliveira, A.K.; Pramoonjago, P.; Rucavado, A.; Moskaluk, C.; Silva, D.T.; Escalante, T.; Gutiérrez, J.M.; Fox, J.W. Mapping the Immune Cell Microenvironment with Spatial Profiling in Muscle Tissue Injected with the Venom of Daboia russelii. Toxins 2023, 15, 208. [Google Scholar] [CrossRef]
  9. Teixeira, C.; Cury, Y.; Moreira, V.; Picolo, G.; Chaves, F. Inflammation induced by Bothrops asper venom. Toxicon 2009, 54, 67–76. [Google Scholar] [CrossRef]
  10. Geiss, G.K.; Bumgarner, R.E.; Birditt, B.; Dahl, T.; Dowidar, N.; Dunaway, D.L.; Fell, H.P.; Ferree, S.; George, R.D.; Grogan, T.; et al. Direct multiplexed measurement of gene expression with color-coded probe pairs. Nat. Biotechnol. 2008, 26, 317–325. [Google Scholar] [CrossRef]
  11. Vider, J.; Croaker, A.; Cox, A.J.; Raymond, E.; Rogers, R.; Adamson, S.; Doyle, M.; O’Brien, B.; Cripps, A.W.; West, N.P. Comparison of skin biopsy sample processing and storage methods on high dimensional immune gene expression using the Nanostring nCounter system. Diagn. Pathol. 2020, 15, 57. [Google Scholar] [CrossRef]
  12. Warrell, D.A. Clinical Toxicology of Snakebite in Asia. In Handbook of Clinical Toxicology of Animal Venoms and Poisons; CRC Press: Boca Raton, FL, USA, 2017; pp. 493–594. [Google Scholar] [CrossRef]
  13. Otero-Patiño, R. Epidemiological, clinical and therapeutic aspects of Bothrops asper bites. Toxicon 2009, 54, 998–1011. [Google Scholar] [CrossRef]
  14. Rucavado, A.; Escalante, T.; Camacho, E.; Gutiérrez, J.M.; Fox, J.W. Systemic vascular leakage induced in mice by Russell’s viper venom from Pakistan. Sci. Rep. 2018, 8, 16088. [Google Scholar] [CrossRef] [PubMed]
  15. Lingam, T.M.C.; Tan, K.Y.; Tan, C.H. Capillary leak syndrome induced by the venoms of Russell’s Vipers (Daboia russelii and Daboia siamensis) from eight locales and neutralization of the differential toxicity by three snake antivenoms. Comp. Biochem. Physiol. C Toxicol. Pharmacol. 2021, 250, 109186. [Google Scholar] [CrossRef]
  16. Rucavado, A.; Escalante, T.; Kalogeropoulos, K.; Camacho, E.; Gutiérrez, J.M.; Fox, J.W. Analysis of wound exudates reveals differences in the patterns of tissue damage and inflammation induced by the venoms of Daboia russelii and Bothrops asper in mice. Toxicon 2020, 186, 94–104. [Google Scholar] [CrossRef]
  17. Wase, N.; Gutiérrez, J.M.; Rucavado, A.; Fox, J.W. Longitudinal Metabolomics and Lipidomics Analyses Reveal Alterations Associated with Envenoming by Bothrops asper and Daboia russelii in an Experimental Murine Model. Toxins 2022, 14, 657. [Google Scholar] [CrossRef] [PubMed]
  18. Escalante, T.; Rucavado, A.; Pinto, A.F.M.; Terra, R.M.S.; Gutiérrez, J.M.; Fox, J.W. Wound exudate as a proteomic window to reveal different mechanisms of tissue damage by snake venom toxins. J. Proteome Res. 2009, 8, 5120–5131. [Google Scholar] [CrossRef]
  19. Rucavado, A.; Escalante, T.; Shannon, J.; Gutiérrez, J.M.; Fox, J.W. Proteomics of wound exudate in snake venom-induced pathology: Search for biomarkers to assess tissue damage and therapeutic success. J. Proteome Res. 2011, 10, 1987–2005. [Google Scholar] [CrossRef]
  20. Escalante, T.; Rucavado, A.; Fox, J.W.; Gutiérrez, J.M. Key events in microvascular damage induced by snake venom hemorrhagic metalloproteinases. J. Proteom. 2011, 74, 1781–1794. [Google Scholar] [CrossRef]
  21. Paes Leme, A.F.; Sherman, N.E.; Smalley, D.M.; Sizukusa, L.O.; Oliveira, A.K.; Menezes, M.C.; Fox, J.W.; Serrano, S.M.T. Hemorrhagic activity of HF3, a snake venom metalloproteinase: Insights from the proteomic analysis of mouse skin and blood plasma. J. Proteome Res. 2012, 11, 279–291. [Google Scholar] [CrossRef]
  22. Freitas-de-sousa, L.A.; Colombini, M.; Lopes-Ferreira, M.; Serrano, S.M.T.; Moura-da-silva, A.M. Insights into the Mechanisms Involved in Strong Hemorrhage and Dermonecrosis Induced by Atroxlysin-Ia, a PI-Class Snake Venom Metalloproteinase. Toxins 2017, 9, 239. [Google Scholar] [CrossRef]
  23. Alape-Girón, A.; Sanz, L.; Escolano, J.; Flores-Díaz, M.; Madrigal, M.; Sasa, M.; Calvete, J.J. Snake venomics of the lancehead pitviper Bothrops asper: Geographic, individual, and ontogenetic variations. J. Proteome Res. 2008, 7, 3556–3571. [Google Scholar] [CrossRef] [PubMed]
  24. Pla, D.; Sanz, L.; Quesada-Bernat, S.; Villalta, M.; Baal, J.; Chowdhury, M.A.W.; León, G.; Gutiérrez, J.M.; Kuch, U.; Calvete, J.J. Phylovenomics of Daboia russelii across the Indian subcontinent. Bioactivities and comparative in vivo neutralization and in vitro third-generation antivenomics of antivenoms against venoms from India, Bangladesh and Sri Lanka. J. Proteom. 2019, 207, 103443. [Google Scholar] [CrossRef] [PubMed]
  25. Senji Laxme, R.R.; Khochare, S.; Attarde, S.; Suranse, V.; Iyer, A.; Casewell, N.R.; Whitaker, R.; Martin, G.; Sunagar, K. Biogeographic venom variation in Russell’s viper (Daboia russelii) and the preclinical inefficacy of antivenom therapy in snakebite hotspots. PLoS Negl. Trop. Dis. 2021, 15, e0009247. [Google Scholar] [CrossRef] [PubMed]
  26. Galvin, N.J.; Vance, P.M.; Dixit, V.M.; Fink, B.; Frazier, W.A. Interaction of human thrombospondin with types I-V collagen: Direct binding and electron microscopy. J. Cell Biol. 1987, 104, 1413–1422. [Google Scholar] [CrossRef] [PubMed]
  27. Murphy-Ullrich, J.E.; Suto, M.J. Thrombospondin-1 regulation of latent TGF-β activation: A therapeutic target for fibrotic disease. Matrix Biol. 2018, 68–69, 28–43. [Google Scholar] [CrossRef]
  28. Lodyga, M.; Hinz, B. TGF-β1—A truly transforming growth factor in fibrosis and immunity. Semin. Cell Dev. Biol. 2020, 101, 123–139. [Google Scholar] [CrossRef] [PubMed]
  29. Harel, T.; Levy-Lahad, E.; Daana, M.; Mechoulam, H.; Horowitz-Cederboim, S.; Gur, M.; Meiner, V.; Elpeleg, O. Homozygous stop-gain variant in LRRC32, encoding a TGFβ receptor, associated with cleft palate, proliferative retinopathy, and developmental delay. Eur. J. Hum. Genet. 2019, 27, 1315–1319. [Google Scholar] [CrossRef]
  30. Poe, A.; Martinez Yus, M.; Wang, H.; Santhanam, L. Lysyl oxidase like-2 in fibrosis and cardiovascular disease. Am. J. Physiol. Cell Physiol. 2023, 325, C694–C707. [Google Scholar] [CrossRef]
  31. Tai, Y.; Woods, E.L.; Dally, J.; Kong, D.; Steadman, R.; Moseley, R.; Midgley, A.C. Myofibroblasts: Function, Formation, and Scope of Molecular Therapies for Skin Fibrosis. Biomolecules 2021, 11, 1095. [Google Scholar] [CrossRef]
  32. Gutiérrez, J.; Ownby, C.L.; Odell, G.V. Skeletal muscle regeneration after myonecrosis induced by crude venom and a myotoxin from the snake Bothrops asper (Fer-de-Lance). Toxicon 1984, 22, 719–731. [Google Scholar] [CrossRef]
  33. Arce, V.; Brenes, F.; Gutiérrez, J.M.; Arce, D.V. Degenerative and regenerative changes in murine skeletal muscle after injection of venom from the snake Bothrops asper: A histochemical and immunocytochemical study. Int. J. Exp. Pathol. 1991, 72, 211. [Google Scholar]
  34. Hernández, R.; Cabalceta, C.; Saravia-Otten, P.; Chaves, A.; Gutiérrez, J.M.; Rucavado, A. Poor regenerative outcome after skeletal muscle necrosis induced by Bothrops asper venom: Alterations in microvasculature and nerves. PLoS ONE 2011, 6, e19834. [Google Scholar] [CrossRef]
  35. Petricevich, V.L.; Teixeira, C.F.P.; Tambourgi, D.V.; Gutiérrez, J.M. Increments in serum cytokine and nitric oxide levels in mice injected with Bothrops asper and Bothrops jararaca snake venoms. Toxicon 2000, 38, 1253–1266. [Google Scholar] [CrossRef]
  36. Wanderley, C.W.S.; Silva, C.M.S.; Wong, D.V.T.; Ximenes, R.M.; Morelo, D.F.C.; Cosker, F.; Aragao, K.S.; Fernandes, C.; Palheta-Junior, R.C.; Havt, A.; et al. Bothrops jararacussu snake venom-induces a local inflammatory response in prostanoid-and neuthrophil-dependent manner. Toxicon 2014, 90, 134–147. [Google Scholar] [CrossRef] [PubMed]
  37. Chaves, F.; Teixeira, C.F.P.; Gutiérrez, J.M. Role of TNF-α, IL-1β and IL-6 in the local tissue damage induced by Bothrops asper snake venom: An experimental assessment in mice. Toxicon 2005, 45, 171–178. [Google Scholar] [CrossRef] [PubMed]
  38. Kaneko, N.; Kurata, M.; Yamamoto, T.; Morikawa, S.; Masumoto, J. The role of interleukin-1 in general pathology. Inflamm. Regen. 2019, 39, 12. [Google Scholar] [CrossRef]
  39. Pyrillou, K.; Burzynski, L.C.; Clarke, M.C.H. Alternative Pathways of IL-1 Activation, and Its Role in Health and Disease. Front. Immunol. 2020, 11, 613170. [Google Scholar] [CrossRef]
  40. Mantovani, A. Cancer: An infernal triangle. Nature 2007, 448, 547–548. [Google Scholar] [CrossRef] [PubMed]
  41. Weinheimer-Haus, E.M.; Mirza, R.E.; Koh, T.J. Nod-like receptor protein-3 inflammasome plays an important role during early stages of wound healing. PLoS ONE 2015, 10, e0119106. [Google Scholar] [CrossRef]
  42. Gschwandtner, M.; Derler, R.; Midwood, K.S. More Than Just Attractive: How CCL2 Influences Myeloid Cell Behavior Beyond Chemotaxis. Front. Immunol. 2019, 10, 2759. [Google Scholar] [CrossRef]
  43. Martinez, F.O.; Sica, A.; Mantovani, A.; Locati, M. Macrophage activation and polarization. Front. Biosci. 2008, 13, 453–461. [Google Scholar] [CrossRef]
  44. Mantovani, A.; Sica, A.; Sozzani, S.; Allavena, P.; Vecchi, A.; Locati, M. The chemokine system in diverse forms of macrophage activation and polarization. Trends Immunol. 2004, 25, 677–686. [Google Scholar] [CrossRef] [PubMed]
  45. Vazirinejad, R.; Ahmadi, Z.; Arababadi, M.K.; Hassanshahi, G.; Kennedy, D. The biological functions, structure and sources of CXCL10 and its outstanding part in the pathophysiology of multiple sclerosis. Neuroimmunomodulation 2014, 21, 322–330. [Google Scholar] [CrossRef]
  46. Ridiandries, A.; Tan, J.T.M.; Bursill, C.A. The Role of Chemokines in Wound Healing. Int. J. Mol. Sci. 2018, 19, 3217. [Google Scholar] [CrossRef]
  47. Jiang, Z.; Georgel, P.; Du, X.; Shamel, L.; Sovath, S.; Mudd, S.; Huber, M.; Kalis, C.; Keck, S.; Galanos, C.; et al. CD14 is required for MyD88-independent LPS signaling. Nat. Immunol. 2005, 6, 565–570. [Google Scholar] [CrossRef] [PubMed]
  48. Zamani, F.; Shahneh, F.Z.; Aghebati-Maleki, L.; Baradaran, B. Induction of CD14 Expression and Differentiation to Monocytes or Mature Macrophages in Promyelocytic Cell Lines: New Approach. Adv. Pharm. Bull. 2013, 3, 329–332. [Google Scholar] [CrossRef]
  49. Sharygin, D.; Koniaris, L.G.; Wells, C.; Zimmers, T.A.; Hamidi, T. Role of CD14 in human disease. Immunology 2023, 169, 260–270. [Google Scholar] [CrossRef] [PubMed]
  50. Zoccal, K.F.; Bitencourt, C.D.S.; Paula-Silva, F.W.G.; Sorgi, C.A.; De Castro Figueiredo Bordon, K.; Arantes, E.C.; Faccioli, L.H. TLR2, TLR4 and CD14 recognize venom-associated molecular patterns from Tityus serrulatus to induce macrophage-derived inflammatory mediators. PLoS ONE 2014, 9, e88174. [Google Scholar] [CrossRef]
  51. Ikenohuchi, Y.J.; Silva, M.D.S.; Rego, C.M.A.; Francisco, A.F.; da Silva Setúbal, S.; Ferreira e Ferreira, A.A.; Boeno, C.N.; Santana, H.M.; Felipin, K.P.; de Lima, A.M.; et al. A C-type lectin induces NLRP3 inflammasome activation via TLR4 interaction in human peripheral blood mononuclear cells. Cell. Mol. Life Sci. 2023, 80, 188. [Google Scholar] [CrossRef]
  52. Fontana, B.C.; Soares, A.M.; Zuliani, J.P.; Gonçalves, G.M. Role of Toll-like receptors in local effects in a model of experimental envenoming induced by Bothrops jararacussu snake venom and by two phospholipases A2. Toxicon 2022, 214, 145–154. [Google Scholar] [CrossRef]
  53. Jiménez, N.; Escalante, T.; Gutiérrez, J.M.; Rucavado, A. Skin pathology induced by snake venom metalloproteinase: Acute damage, revascularization, and re-epithelization in a mouse ear model. J. Investig. Dermatol. 2008, 128, 2421–2428. [Google Scholar] [CrossRef]
  54. Mora, R.; Valverde, B.; Díaz, C.; Lomonte, B.; Gutiérrez, J.M. A Lys49 phospholipase A(2) homologue from Bothrops asper snake venom induces proliferation, apoptosis and necrosis in a lymphoblastoid cell line. Toxicon 2005, 45, 651–660. [Google Scholar] [CrossRef]
  55. Bezerra, P.H.A.; Ferreira, I.M.; Franceschi, B.T.; Bianchini, F.; Ambrósio, L.; Cintra, A.C.O.; Sampaio, S.V.; De Castro, F.A.; Torqueti, M.R. BthTX-I from Bothrops jararacussu induces apoptosis in human breast cancer cell lines and decreases cancer stem cell subpopulation. J. Venom. Anim. Toxins Incl. Trop. Dis. 2019, 25, e20190010. [Google Scholar] [CrossRef]
  56. Costal-Oliveira, F.; Stransky, S.; Guerra-Duarte, C.; Naves de Souza, D.L.; Vivas-Ruiz, D.E.; Yarlequé, A.; Sanchez, E.F.; Chávez-Olórtegui, C.; Braga, V.M.M. L-amino acid oxidase from Bothrops atrox snake venom triggers autophagy, apoptosis and necrosis in normal human keratinocytes. Sci. Rep. 2019, 9, 781. [Google Scholar] [CrossRef] [PubMed]
  57. Burin, S.M.; Da Costa Cacemiro, M.; Cominal, J.G.; De Grandis, R.A.; Machado, A.R.T.; Donaires, F.S.; Cintra, A.C.O.; Ambrosio, L.; Antunes, L.M.G.; Sampaio, S.V.; et al. Bothrops moojeni L-amino acid oxidase induces apoptosis and epigenetic modulation on Bcr-Abl+ cells. J. Venom. Anim. Toxins Incl. Trop. Dis. 2020, 26, e20200123. [Google Scholar] [CrossRef] [PubMed]
  58. Saleem, M.; Qadir, M.I.; Perveen, N.; Ahmad, B.; Saleem, U.; Irshad, T. Inhibitors of apoptotic proteins: New targets for anticancer therapy. Chem. Biol. Drug Des. 2013, 82, 243–251. [Google Scholar] [CrossRef] [PubMed]
  59. Wolf, P. Inhibitor of apoptosis proteins as therapeutic targets in bladder cancer. Front. Oncol. 2023, 13, 1124600. [Google Scholar] [CrossRef]
  60. Levine, B.; Kroemer, G. Biological Functions of Autophagy Genes: A Disease Perspective. Cell 2019, 176, 11–42. [Google Scholar] [CrossRef] [PubMed]
  61. Komamura, K.; Shirotani-Ikejima, H.; Tatsumi, R.; Tsujita-Kuroda, Y.; Kitakaze, M.; Miyatake, K.; Sunagawa, K.; Miyata, T. Differential gene expression in the rat skeletal and heart muscle in glucocorticoid-induced myopathy: Analysis by microarray. Cardiovasc. Drugs Ther. 2003, 17, 303–310. [Google Scholar] [CrossRef] [PubMed]
  62. Janssens, V.; Goris, J. Protein phosphatase 2A: A highly regulated family of serine/threonine phosphatases implicated in cell growth and signalling. Biochem. J. 2001, 353, 417–439. [Google Scholar] [CrossRef]
  63. Scheffler, J.M.; Sparber, F.; Tripp, C.H.; Herrmann, C.; Humenberger, A.; Blitz, J.; Romani, N.; Stoitzner, P.; Huber, L.A. LAMTOR2 regulates dendritic cell homeostasis through FLT3-dependent mTOR signalling. Nat. Commun. 2014, 5, 5138. [Google Scholar] [CrossRef]
  64. Łyszkiewicz, M.; Kotlarz, D.; Ziȩtara, N.; Brandes, G.; Diestelhorst, J.; Glage, S.; Hobeika, E.; Reth, M.; Huber, L.A.; Krueger, A.; et al. LAMTOR2 (p14) Controls B Cell Differentiation by Orchestrating Endosomal BCR Trafficking. Front. Immunol. 2019, 10, 497. [Google Scholar] [CrossRef]
  65. Renna, M.; Bento, C.F.; Fleming, A.; Menzies, F.M.; Siddiqi, F.H.; Ravikumar, B.; Puri, C.; Garcia-Arencibia, M.; Sadiq, O.; Corrochano, S.; et al. IGF-1 receptor antagonism inhibits autophagy. Hum. Mol. Genet. 2013, 22, 4528–4544. [Google Scholar] [CrossRef] [PubMed]
  66. Chen, P.C.; Kuo, Y.C.; Chuong, C.M.; Huang, Y.H. Niche Modulation of IGF-1R Signaling: Its Role in Stem Cell Pluripotency, Cancer Reprogramming, and Therapeutic Applications. Front. Cell Dev. Biol. 2021, 8, 625943. [Google Scholar] [CrossRef] [PubMed]
  67. Tang, Z.; Takahashi, Y.; Chen, C.; Liu, Y.; He, H.; Tsotakos, N.; Serfass, J.M.; Gebru, M.T.; Chen, H.; Young, M.M.; et al. Atg2A/B deficiency switches cytoprotective autophagy to non-canonical caspase-8 activation and apoptosis. Cell Death Differ. 2017, 24, 2127–2138. [Google Scholar] [CrossRef] [PubMed]
  68. Wang, S.; Li, H.; Yuan, M.; Fan, H.; Cai, Z. Role of AMPK in autophagy. Front. Physiol. 2022, 13, 2479. [Google Scholar] [CrossRef]
Figure 1. Light micrographs of sections of skeletal muscle from mice injected in the right gastrocnemius muscle with (A) PBS, (B) 30 µg B. asper venom dissolved in 50 µL PBS, or (C) 30 µg D. russelii venom dissolved in 50 µL PBS. Mice were sacrificed 24 h after injection, and the gastrocnemius muscles were excised, placed in formalin fixative, and routinely processed for embedding in paraffin (see Section 2 for details). Muscle tissue sections (4 µm thickness) were stained with Mallory’s trichrome for histological observation of tissue damage and inflammatory infiltrate. Both venoms induced muscle fiber necrosis (arrows), but only B. asper venom induced hemorrhage, as evidenced by the presence of erythrocytes in the interstitial space (*). An inflammatory infiltrate of leukocytes is present in envenomed tissue. Bars correspond to 100 µm.
Figure 1. Light micrographs of sections of skeletal muscle from mice injected in the right gastrocnemius muscle with (A) PBS, (B) 30 µg B. asper venom dissolved in 50 µL PBS, or (C) 30 µg D. russelii venom dissolved in 50 µL PBS. Mice were sacrificed 24 h after injection, and the gastrocnemius muscles were excised, placed in formalin fixative, and routinely processed for embedding in paraffin (see Section 2 for details). Muscle tissue sections (4 µm thickness) were stained with Mallory’s trichrome for histological observation of tissue damage and inflammatory infiltrate. Both venoms induced muscle fiber necrosis (arrows), but only B. asper venom induced hemorrhage, as evidenced by the presence of erythrocytes in the interstitial space (*). An inflammatory infiltrate of leukocytes is present in envenomed tissue. Bars correspond to 100 µm.
Biomolecules 14 00278 g001
Figure 2. NanoString GSA pathway analysis of extracellular matrix synthesis (a), extracellular matrix degradation (b), collagen biosynthesis and modification (c), and myofibroblast regulation (d) genes alterations in muscle tissue injected with venoms of B. asper and D. russelii, as compared to controls injected with PBS. The Directed Global Significance scores are represented by up-regulated (red) and down-regulated (blue) genes at different time points in relation to the control samples. The p-value and Log2 FoldChange are shown to the right. Global Significance Score > 1.0. Graphic created in the Rosalind software (https://rosalind.bio/ version 3.38.6.1, accessed on 5 March 2023).
Figure 2. NanoString GSA pathway analysis of extracellular matrix synthesis (a), extracellular matrix degradation (b), collagen biosynthesis and modification (c), and myofibroblast regulation (d) genes alterations in muscle tissue injected with venoms of B. asper and D. russelii, as compared to controls injected with PBS. The Directed Global Significance scores are represented by up-regulated (red) and down-regulated (blue) genes at different time points in relation to the control samples. The p-value and Log2 FoldChange are shown to the right. Global Significance Score > 1.0. Graphic created in the Rosalind software (https://rosalind.bio/ version 3.38.6.1, accessed on 5 March 2023).
Biomolecules 14 00278 g002aBiomolecules 14 00278 g002b
Figure 3. Heat maps of the extracellular matrix and myofibroblast regulation genes expression in muscle tissue injected with PBS, B. asper venom, or D. russelii venom. Heat maps were made using Log2 of transcript counts of the up- and down-regulated significant genes in each pathway using Gene Set Analysis (GSA). Each row represents the gene monitoring over time after venom injection. Euclidean distance and complete linkage methods were applied. Hierarchical clustering was performed with SRPLOT (accessed on 1 January 2023 at https://www.bioinformatics.com.cn).
Figure 3. Heat maps of the extracellular matrix and myofibroblast regulation genes expression in muscle tissue injected with PBS, B. asper venom, or D. russelii venom. Heat maps were made using Log2 of transcript counts of the up- and down-regulated significant genes in each pathway using Gene Set Analysis (GSA). Each row represents the gene monitoring over time after venom injection. Euclidean distance and complete linkage methods were applied. Hierarchical clustering was performed with SRPLOT (accessed on 1 January 2023 at https://www.bioinformatics.com.cn).
Biomolecules 14 00278 g003
Figure 4. NanoString GSA pathway analysis of cytokine signaling (a), Toll-like receptors (b), NK-κB (c), and TGFβ (d) genes in muscle tissue injected with venoms of B. asper and D. russelii, as compared to controls injected with PBS. The Directed Global Significance scores are represented by up-regulated (red) and down-regulated (blue) genes at different time points in relation to the control samples. The p-value and Log2 FoldChange are shown to the right. Global Significance Score > 1.0. Graphic created in the Rosalind software (https://rosalind.bio/ version 3.38.6.1, accessed on 5 March 2023).
Figure 4. NanoString GSA pathway analysis of cytokine signaling (a), Toll-like receptors (b), NK-κB (c), and TGFβ (d) genes in muscle tissue injected with venoms of B. asper and D. russelii, as compared to controls injected with PBS. The Directed Global Significance scores are represented by up-regulated (red) and down-regulated (blue) genes at different time points in relation to the control samples. The p-value and Log2 FoldChange are shown to the right. Global Significance Score > 1.0. Graphic created in the Rosalind software (https://rosalind.bio/ version 3.38.6.1, accessed on 5 March 2023).
Biomolecules 14 00278 g004aBiomolecules 14 00278 g004b
Figure 5. Heat maps of the cytokine signaling, Toll-like receptors, NK-κB, and TGFβgenes alterations in samples from mice injected with either PBS (control), B. asper, or D. russelii venoms. The heat maps were made using Log2 of transcript counts of the up- and down-regulated significant genes in each pathway using the Gene Set Analysis (GSA). Each row represents the gene monitoring over time after venom injection. Euclidean distance and complete linkage methods were applied. Hierarchical clustering was performed with SRPLOT (accessed on 1 September 2023 at https://www.bioinformatics.com.cn).
Figure 5. Heat maps of the cytokine signaling, Toll-like receptors, NK-κB, and TGFβgenes alterations in samples from mice injected with either PBS (control), B. asper, or D. russelii venoms. The heat maps were made using Log2 of transcript counts of the up- and down-regulated significant genes in each pathway using the Gene Set Analysis (GSA). Each row represents the gene monitoring over time after venom injection. Euclidean distance and complete linkage methods were applied. Hierarchical clustering was performed with SRPLOT (accessed on 1 September 2023 at https://www.bioinformatics.com.cn).
Biomolecules 14 00278 g005
Figure 6. NanoString GSA pathway analysis of programmed cell death and autophagy gene alterations in muscle tissue injected with venoms of B. asper and D. russelii, as compared to controls injected with PBS. NanoString GSA pathway analysis. The Directed Global Significance scores are represented by up-regulated (red) and down-regulated (blue) genes at different time points in relation to the control samples. The p-value and Log2 FoldChange are shown to the right. Global Significance Score > 1.0. Graphic created in the Rosalind software (https://rosalind.bio/ version 3.38.6.1, accessed on 5 March 2023).
Figure 6. NanoString GSA pathway analysis of programmed cell death and autophagy gene alterations in muscle tissue injected with venoms of B. asper and D. russelii, as compared to controls injected with PBS. NanoString GSA pathway analysis. The Directed Global Significance scores are represented by up-regulated (red) and down-regulated (blue) genes at different time points in relation to the control samples. The p-value and Log2 FoldChange are shown to the right. Global Significance Score > 1.0. Graphic created in the Rosalind software (https://rosalind.bio/ version 3.38.6.1, accessed on 5 March 2023).
Biomolecules 14 00278 g006
Figure 7. Heat maps of the programmed cell death and autophagy gene alterations in samples from mice injected with either PBS (control), B. asper, or D. russelii venoms. The heat maps were made using Log2 of transcript counts of the up-and down-regulated significant genes in each pathway using the Gene Set Analysis (GSA). Each row represents the gene monitoring over time after venom injection. Euclidean distance and complete linkage methods were applied. Hierarchical clustering is performed with SRPLOT (accessed on 1 September 2023 at https://www.bioinformatics.com.cn).
Figure 7. Heat maps of the programmed cell death and autophagy gene alterations in samples from mice injected with either PBS (control), B. asper, or D. russelii venoms. The heat maps were made using Log2 of transcript counts of the up-and down-regulated significant genes in each pathway using the Gene Set Analysis (GSA). Each row represents the gene monitoring over time after venom injection. Euclidean distance and complete linkage methods were applied. Hierarchical clustering is performed with SRPLOT (accessed on 1 September 2023 at https://www.bioinformatics.com.cn).
Biomolecules 14 00278 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Oliveira, A.K.d.; Rucavado, A.; Escalante, T.; Gutiérrez, J.M.; Fox, J.W. A Complex Pattern of Gene Expression in Tissue Affected by Viperid Snake Envenoming: The Emerging Role of Autophagy-Related Genes. Biomolecules 2024, 14, 278. https://doi.org/10.3390/biom14030278

AMA Style

Oliveira AKd, Rucavado A, Escalante T, Gutiérrez JM, Fox JW. A Complex Pattern of Gene Expression in Tissue Affected by Viperid Snake Envenoming: The Emerging Role of Autophagy-Related Genes. Biomolecules. 2024; 14(3):278. https://doi.org/10.3390/biom14030278

Chicago/Turabian Style

Oliveira, Ana Karina de, Alexandra Rucavado, Teresa Escalante, José María Gutiérrez, and Jay W. Fox. 2024. "A Complex Pattern of Gene Expression in Tissue Affected by Viperid Snake Envenoming: The Emerging Role of Autophagy-Related Genes" Biomolecules 14, no. 3: 278. https://doi.org/10.3390/biom14030278

APA Style

Oliveira, A. K. d., Rucavado, A., Escalante, T., Gutiérrez, J. M., & Fox, J. W. (2024). A Complex Pattern of Gene Expression in Tissue Affected by Viperid Snake Envenoming: The Emerging Role of Autophagy-Related Genes. Biomolecules, 14(3), 278. https://doi.org/10.3390/biom14030278

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop