- Shaanxi Key Laboratory for Animal Conservation, Shaanxi Institute of Zoology, Xi'an, China
Pneumonia can seriously threaten the life of forest musk deer (FMD, an endangered species). To gain a comprehensive understanding of pneumonia pathogenesis in FMD, iTRAQ-based proteomics analysis was performed in diseased (Pne group) lung tissues of FMD that died of pneumonia and normal lung tissues (Ctrl group) of FMD that died from fighting against each other. Results showed that 355 proteins were differentially expressed (fold change ≥ 1.2 and adjusted P-value < 0.05) in Pne vs. Ctrl. GO/KEGG annotation and enrichment analyses showed that dysregulated proteins might play vital roles in bacterial infection and immunity. Given the close association between bacterial infection and pneumonia, 32 dysregulated proteins related to Staphylococcus aureus infection, bacterial invasion of epithelial cells, and pathogenic Escherichia coli infection were screened out. Among these 32 proteins, 13 proteins were mapped to the bovine genome. Given the close phylogenetic relationships of FMD and bovine, the protein-protein interaction networks of the above-mentioned 13 proteins were constructed by the String database. Based on the node degree analysis, 5 potential key proteins related to pneumonia-related bacterial infection in FMD were filtered out. Moreover, 85 dysregulated proteins related to the immune system process were identified given the tight connection between immune dysregulation and pneumonia pathogenesis. Additionally, 12 proteins that might function as crucial players in pneumonia-related immune response in FMD were screened out using the same experimental strategies described above. In conclusion, some vital proteins, biological processes, and pathways in pneumonia development were identified in FMD.
Introduction
Forest musk deer (FMD, Moschus berezovskii) is a type of mammal that mainly lives in the alpine forests in China and Vietnam (1, 2). The populations of FMD are sharply declined since the 1950s due to habitat destruction/degradation and massive hunting for their musk (Moschus) (3, 4). Musk is the dried secretion from the musk sac gland of male musk deer, such as Moschus berezovskii Flerov, Moschus moschiferus Linnaeus, and Moschus sifanicus Przewalski (5, 6). Musk is a superior component in perfume and is believed to have potential therapeutic values for multiple diseases such as cancers, strokes, and heart diseases in the traditional Asian medicine industry (6, 7). This species has been listed in the First-Class National Protected Animal List of China and is protected under the Chinese Wild Animal Protection Law (4, 8). Moreover, FMD is the major species of musk deer that can be reared artificially in special farms under the support of the Chinese government, which contributes to the growth of the population, reduction of poaching behaviors, and better utilization of FMD resources (9, 10).
The increase in the population of FMD is also limited by some fatal diseases including pneumonia (11–13). Pneumonia can be caused by multiple pathogens including bacteria (14–16). However, previous research on pneumonia mainly focused on the isolation, identification, and genome analysis of pathogens in FMD (12, 17, 18). To reduce the risk and harm of pneumonia for the health and life of FMD, it is imperative to have an in-depth insight into the molecular mechanisms underlying pneumonia development and identify key molecules or pathways related to pneumonia pathogenesis.
Recently, mass spectrometry (MS)-based proteomics has attracted much attention from researchers because proteins are responsible for most biological functions and proteomics can simultaneously capture and quantify thousands of proteins rather than RNAs in a cost-effective manner (19, 20). Isobaric tag for relative and absolute quantification (iTRAQ), an isotope labeling strategy, has been widely used in proteomics studies needing relative quantification due to the multiple advantages such as multiplexing capacity, reproducibility, easy operation, and flexibility (21–23). The combination of iTRAQ and MS-based proteomics technologies and bioinformatics analytical methods have emerged as a powerful strategy for identifying vital proteins related to disease pathogenesis, comprehensively understanding protein roles and basic biological functions, and deciphering complicated molecular mechanisms underlying disease development in multiple animals (24–26).
However, to our knowledge, there is no proteomics data to explore the pathogenesis of pneumonia in FMD to date. To build up a general and comprehensive understanding of pneumonia pathogenesis, the iTRAQ-based LC-MS/MS technique was used to explore the proteomics alterations in diseased lung tissues of FMD that died of pneumonia than in normal lung tissues of FMD died from fighting against each other. Moreover, some genes/proteins, biological processes, and signaling pathways that might play vital roles in pneumonia progression were screened out based on differential expression, annotation, enrichment, and protein-protein interaction analysis.
Materials and methods
Animal samples
Forest musk deer are reared in the Shaanxi Institute of Zoology (China). The Animal breeding area (34.210832°N, 106.902117°E) is located in Fengxian, Southwest of BaoJi City, Shaanxi Province, China, a region of Qinling mountain at an altitude of 1,500 m. Diseased lung tissues were obtained from 3 adult FMD (2♂1♀, 4.5 years old) that died of pneumonia. Normal lung tissues were obtained from 3 adult FMD (♂, 3.5 years old) that died from fighting against each other. Tissue samples with a weight of no <200 mg were stored at −80°C. The study was approved by the Academic Committee of Shaanxi Institute of Zoology with Ethical Approval No.: 20210327001.
Histological analysis
The tissues mentioned above were fixed in a PBS buffered formaldehyde solution for 48 h. After routine dehydration and transparency, sectioned at a thickness of 4 μm and stained with Eosin Staining Kit (Beyotime, Shanghai, China) following the protocols of the manufacturer, and examined by light microscopy.
ITRAQ-based proteomics analysis
ITRAQ-based proteomics analysis was performed as previously described (27, 28). The detailed experimental procedures of proteomics analysis including sample preparation, iTRAQ labeling and fractionation, LC-MS/MS analysis, and data analysis were shown in Supplementary file 1. Briefly, tissue samples were ground to a fine powder in liquid nitrogen and then lysed using the protein lysis buffer [7M Urea/4% SDS/2M Thiourea/40 mM Tris-HCl (pH 8.5)] supplemented with 2 mM EDTA and 1mM phenylmethylsulfonyl fluoride. Samples were labeled using the iTRAQ Reagent-8 plex Multiplex Kit (SCIEX, Framingham, MA, USA) according to the protocols of the manufacturer, in which only 6 channels were used in our project. The information of sample-corresponding channels was shown in Table 1. Pne and Ctrl groups represented diseased and normal lung tissue groups, respectively. LC-MS/MS analysis was carried out on TripleTOF 5600+ mass spectrometry (SCIEX) coupled with an EksigentnanoLC system (SCIEX). Raw data analysis was performed using the Protein Pilot Software (version 4.5, SCIEX). The raw MS/MS file data were searched against the PR1-19060015_pep. fasta (containing 24,352 sequences). Proteins were regarded to be significantly differentially expressed when fold-change ≥ 1.2 and adjusted P-value < 0.05.
Bioinformatics and annotations
To determine the biological and functional properties of all the identified proteins, the identified protein sequences were mapped with those in the Swiss-Prot database using BLASTP. In addition, a homology search was performed for the differentially expressed protein sequences using a localized NCBI blastp program against the NCBI non-redundant protein (NR) animal database. Moreover, the GO and KEGG annotation information of matched proteins was extracted. GO and KEGG pathway enrichment analysis was performed using the hypergeometric test. GO and KEGG pathway terms were considered to be significantly enriched at a P-value < 0.05. Protein-protein interaction (PPI) networks were constructed using the STRING database (version: 11.5) (https://cn.string-db.org/).
qRT-PCR analysis for gene expression
Ten mRNAs were randomly selected for expression analysis by qRT-PCR to validate the data. The primer sequences are listed in Table 2. The GAPDH gene was used as the internal control. The total RNA was extracted from the Lung tissues with the RNAiso plus reagent (Takara, Dalian, China) following the manufacturer's protocols. The qRT-PCR was performed using SYBR Premix ExTaq (TaKaRa, Dalian, China) and a Thermal Cycler CFX96 Real Time-PCR detection system (Bio-Rad, Hercules, CA, USA) with the following parameters: 95 °C for 60 s; 40 cycles at 95 °C for 15 s; 60 °C for 30 s; and 72 °C 10 s. The concentration and purity of total RNA were measured using a GE Nanovue™ Spectrophotometer (GE Health care Biosciences, Pittsburgh, USA). cDNA was synthesized using the SYBR Prime Script™ RT Master Mix (Perfect Real Time) Kit (Takara, Dalian, China). The relative expression of each gene was calculated with the 2–ΔΔCt method. There were three biological sample replicates, and each biological sample replicates included three technical replicates.
Results
Histological observation of lung tissue
Histological analysis showed that the alveolar cavity has inflammatory cell exudates and the alveolar wall capillary hyperemia in the pneumonia group (Figures 1C,D). Numerous broken neutrophils were exuded from the alveolar cavity and obvious bleeding was noticed in the pneumonia group (Figures 1E,F). And, the most notable pathological changes were interstitial pneumonia and hemorrhagic pneumonia in the pneumonia group (Figures 1G,H). Moreover, red blood cell, inflammatory cell, and fibrin exudate were present in the alveolar lumen, and the lung interstitium was widened in the pneumonia group (Figures 1G,H).
Figure 1. Histopathological changes in lungs of forest musk deer. (A,B) normal lung tissue in the control group, HE staining of normal tissues was performed (5×). The black arrow indicates the alveolar cavity, the red arrow points to the alveolar wall. (C,D) Inflammatory cells (black arrow) and exudate can be seen in the alveolar cavity (red arrow), and the capillaries in the alveolar wall are hyperemic (yellow arrow). (E,F) Numerous broken neutrophils (yellow arrow) exuded from the alveolar cavity, bleeding obvious (black arrow). (G,H) Red blood cell (yellow arrow), inflammatory cell (black arrow) and fibrin exudates (red arrow) present in alveolar lumen, and lung interstitium is widened (blue arrow). Bars: (B,D,F,H) 50 μm, Bars: (A,C,E,G) 200 μm, (magnification, 5.0×, 20.0×).
Identification of differentially expressed proteins
In our proteomics analysis, 355 proteins (169 down-regulated and 186 up-regulated) were found to be differentially expressed (up-regulated ratio ≥ 1.2 or down-regulated ratio ≤ 0.83; adjusted P-value < 0.05) in the diseased lung tissues of FMDs who died of pneumonia compared to the normal lung tissue group (Figure 2; Supplementary Table 1). The volcano plot of differentially expressed proteins was shown in Figure 2. Among these differentially expressed proteins, 158 proteins were annotated in the bovine Swiss-Prot database.
Figure 2. The volcano plot of differentially expressed and unchanged proteins. Red dots represent the up-regulated proteins and green dots denote the down-regulated proteins.
GO and KEGG annotation analysis of differentially expressed proteins
To screen out key proteins related to the pathogenesis of pneumonia, the sequences of differentially expressed proteins were compared against the NCBI NR database using the NCBI-BLAST. These differentially expressed proteins were also annotated by comparisons against the GO and KEGG databases. Based on the principle of sequence similarity, proteins with similar sequences have similar functions. GO annotation analysis revealed that most of the down-regulated and up-regulated proteins were involved in the regulation of biological processes such as cellular process, metabolic process, biological regulation, response to stimulus, and cellular component organization or biogenesis (Figure 3A; Supplementary Table 2). Also, many differentially expressed proteins were implicated in the immune system process, death, locomotion, cell proliferation, biological adhesion, and growth (Figure 3A; Supplementary Table 2). Moreover, KEGG pathway annotation analysis showed that most up-regulated and down-regulated proteins played crucial roles in the pathways related to focal adhesion, phagosome, microbial metabolism in diverse environments, leukocyte transendothelial migration, bacterial invasion of epithelial cells, endocytosis, Staphylococcus aureus infection, and pathogenic Escherichia coli infection (Supplementary Tables 3, 4). The statistics of the top 20 KEGG pathways of up-regulated and down-regulated proteins were shown in Figure 3B.
Figure 3. GO and KEGG annotation analysis of up-regulated and down-regulated proteins. (A) The percentage of up-regulated and down-regulated proteins in each GO term. (B) The percentage of up-regulated and down-regulated proteins in the top 20 annotated KEGG pathway terms.
GO and KEGG enrichment analysis of differentially expressed proteins
GO enrichment analysis showed that differentially expressed proteins were significantly enriched in biological processes such as acute-phase response, leukocyte adhesion, leukocyte migration, phagocytosis, regulation of tumor necrosis factor biosynthetic process, defense response to Gram-negative bacterium, regulation of locomotion, receptor-mediated endocytosis, cell structure disassembly during apoptosis, defense response to fungus (Supplementary Table 5). The top 20 GO biological process terms that were significantly enriched by differentially expressed genes were displayed in Figure 4A. KEGG enrichment analysis disclosed that differentially expressed proteins were significantly enriched in pathways related to Staphylococcus aureus infection, focal adhesion, complement and coagulation cascades, phagosome, antigen processing and presentation, bacterial invasion of epithelial cells, and pathogenic Escherichia coli infection (Supplementary Table 6). The top 20 KEGG pathways that were significantly enriched by the differentially expressed proteins were shown in Figure 4B.
Figure 4. GO and KEGG enrichment analysis of differentially expressed proteins. (A) The top 20 GO biological process terms that were significantly enriched by the differentially expressed proteins. (B) The top 20 KEGG pathways that were significantly enriched by the differentially expressed proteins.
Screening and PPI network construction of dysregulated proteins related to bacterial infection
Both KEGG annotation and enrichment outcomes suggested that pathways related to Staphylococcus aureus infection, bacterial invasion of epithelial cells, and pathogenic Escherichia coli infection might play vital roles in pneumonia progression. Given the close association between bacterial infection and pneumonia pathogenesis, dysregulated proteins in the above-mentioned pathways were filtered out based on KEGG annotation analysis. The information on these proteins was shown in Supplementary Table 7. As presented in Supplementary Table 7, 13 (9 down-regulated and 4 up-regulated), 14 (7 down-regulated and 7 up-regulated), or 12 (4 down-regulated and 8 up-regulated) differentially expressed proteins were identified to be implicated in Staphylococcus aureus infection, bacterial invasion of epithelial cells, or pathogenic Escherichia coli infection, respectively. The above-mentioned differentially expressed proteins related to bacterial infection (total number: 32) were integrated into Supplementary Table 8. The STRING database has been widely used to construct the PPI network and identify hub proteins in previous studies (29, 30). Prior phylogenetic tree analysis showed that FMD was a member of the suborder Ruminantia and order Artiodactyla with close phylogenetic relationships with four members of the family Bovidae (sheep, yak, cattle, and Tibetan antelope) (31). Also, a recent study showed that most FMD unigenes that were identified by de novo assembly of heart and musk gland transcriptomes were homologous with bovine genes (32). Given the close phylogenetic relationships of FMD and bovine, the PPI networks of filtered proteins were constructed based on the information of the organism species Bos taurus (bovine). Among filtered 32 differentially expressed proteins related to bacterial infection, 13 proteins were annotated in the bovine Swiss-Prot database (Supplementary Table 8). The PPI networks of the 13 proteins were constructed and displayed in Figure 4 (organism: Bos taurus, combined interaction score ≥ 0.4) (the solitary proteins had been removed from the network). The detailed interaction information of these 13 proteins was shown in Supplementary Table 9. The node degrees (number of interacted proteins) of proteins in the PPI networks (Figure 5) were exhibited in Supplementary Table 10. The node degree can be used to identify hub proteins in the PPI networks (33, 34). Results suggested that 5 proteins with greater node degrees [i.e. catenin beta-1 (CTNNB1), integrin beta-1 (ITGB1), catenin alpha-1 (CTNNA1), dynamin-2 (DNM2), Keratin, type I cytoskeletal 19 (KRT19)] might function as crucial players in pneumonia-related bacterial infection in FMD.
Figure 5. The PPI networks of differentially expressed proteins (species: bovine) related to bacterial infection. The darkturquoise and pink lines represented the interactions that were identified from curated databases and experimentally determined, respectively. The green, red, and blue lines represented the predicted interactions based on gene neighborhood, gene fusions, and gene co-occurrence relationships, respectively. The yellow-green, black, and purple lines denoted the potential interactions that were identified by textmining, co-expression, and protein homology. The full names of proteins in this figure were shown as follows: C3, complement C3; CFH, complement factor H, KRT19, keratin, type I cytoskeletal 19; KRT25, keratin, type I cytoskeletal 25; DNM2: dynamin-2; ILK: integrin-linked protein kinase; CLTC, clathrin heavy chain 1; CTNNA1, catenin alpha-1; ITGB1, integrin beta-1; CTNNB1, catenin beta-1; BCAM, basal cell adhesion molecule; TUBA4A, tubulin alpha-4A chain; EZR, ezrin.
Screening and PPI network construction of dysregulated proteins related to immunity
It has been reported that the pathogenesis of pneumonia is closely linked with the dysfunction of the immune system (35, 36). In this project, 53 down-regulated and 32 up-regulated proteins that were implicated in the immune system process were screened out based on the GO annotation analysis. These 85 proteins related to the immune system process were shown in Supplementary Table 11. Among these 85 proteins, 49 proteins that were mapped to the bovine genome were screened out (Supplementary Table 11). Next, the PPI networks of these 49 proteins were established and presented in Figure 6 (organism: Bos taurus, combined interaction score ≥ 0.4) (the solitary proteins have been removed from the network). The detailed protein-protein interaction information and node degrees of the above-mentioned 49 proteins in the PPI networks were displayed in Supplementary Tables 12, 13, respectively. The outcomes suggested that CTNNB1, ITGB1, Annexin A5 (ANXA5), calreticulin (CALR), prothrombin (F2), matrix metalloproteinase-9 (MMP9), platelet endothelial cell adhesion molecule (PECAM1), thrombospondin-1 (THBS1), heat shock protein HSP 90-beta (HSP90AB1), endoplasmin (HSP90B1), integrin alpha-3 (ITGA3), and moesin (MSN) might be the hub proteins in the PPI networks because they had greater node degrees. In other words, these proteins might play vital roles in the immune response related to pneumonia in FMD.
Figure 6. The PPI networks of differentially expressed proteins (species: bovine) related to immunity. The full names of the proteins in this figure were presented as below: CALR, calreticulin; LGALS9, galectin-9; C3, complement C3; EPB42, protein 4.2; HSP90B1, endoplasmin; CORO1A, coronin-1A; ANXA1, annexin A1; MMP9, matrix metalloproteinase-9; AP3B1,AP-3 complex subunit beta-1; BPI, bactericidal permeability-increasing protein; THBS1, thrombospondin-1; ANG1, angiogenin-1; S100A12, protein S100-A12; LBP, lipopolysaccharide-binding protein; ANXA5, annexin A5; S100A8, protein S100-A8; S100A9, protein S100-A9; HSP90AB1, heat shock protein HSP 90-beta; DHX9,ATP-dependent RNA helicase A; CAPZA1, f-actin-capping protein subunit alpha-1; GPI, glucose-6-phosphate isomerase; CLTC, clathrin heavy chain 1; DNM2,dynamin-2;F2, prothrombin; THBD, thrombomodulin (fragment); CFH, complement factor H; EEF2, elongation factor 2; BPIFB1, BPI fold-containing family B member 1; PDIA3, protein disulfide-isomerase A3; FLOT2, flotillin-2; UBA52, ubiquitin-60S ribosomal protein L40; MSN, moesin; PIGR, polymeric immunoglobulin receptor; ITGB1, integrin beta-1; LGMN, legumain; COLEC12, collectin-12; CTNNB1, catenin beta-1; SFXN1, sideroflexin-1; BCAM, basal cell adhesion molecule; HMOX1, heme oxygenase 1; ANXA3, annexin A3; MYO1C, unconventional myosin-Ic; CTSD, cathepsin D; APOA1, apolipoprotein A-I; DPP4, dipeptidyl peptidase 4; ITGA3, integrin alpha-3; EZR, ezrin; PECAM1, platelet endothelial cell adhesion molecule; FTH1, ferritin heavy chain.
Validation of differentially expressed proteins by qRT-PCR analysis
The expression patterns determined by qRT-PCR were consistent with those obtained by iTRAQ, with 90% agreement between the qRT-PCR and iTRAQ results (Figure 7). This result indicated that the differential proteomic analysis outcomes in this study were reliable.
Figure 7. Quantitative RT-PCR analyses of gene expression in normal and diseased tissues. Quantitative expression patterns of genes, which was calculated based on Ct value normalized against the housekeeping GAPDH gene. The *, **, and *** symbols indicate the P-value < 0.05, P-value < 0.01, and P-value < 0.001 respectively.
Discussion
In this project, a total of 355 differentially expressed (up-regulated ratio ≥ 1.2 or down-regulated ratio ≤ 0.83; adjusted P-value < 0.05) proteins were identified in the diseased lung tissues of FMDs who died of pneumonia vs. the normal group. Among these dysregulated proteins, 158 proteins were mapped to the bovine genome.
Moreover, KEGG pathway annotation analysis showed that 9 pathways (i.e., metabolic pathways, focal adhesion, pathways in cancer, regulation of actin cytoskeleton, amoebiasis, microbial metabolism in diverse environments, tight junction, leukocyte transendothelial migration, bacterial invasion of epithelial cells) were shared in the top 20 KEGG pathways of up-regulated and down-regulated proteins. Among the 46 metabolic pathways-related proteins (28 up-regulated and 18 down-regulated), alpha-enolase (ENO1), neutrophil gelatinase-associated lipocalin (LCN2), and Acetyl-CoA acetyltransferase (ACAT1) have been found to be related to pathogens-induced pneumonia. For instance, ENO1 facilitated lipopolysaccharide (LPS)-driven monocyte recruitment to the acutely inflamed lung, and ENO1 was highly expressed in blood monocyte cell surface and alveolar mononuclear cells of patients with pneumonia (37). LCN2 had a potential protective effect against Escherichia coli-induced pneumonia (38). LCN2 knockout notably improved the susceptibility of mice to Acinetobacter baumannii pneumonia (39). LCN2 hindered the clearance of pneumococcal pneumonia and exacerbated pneumococcal pneumonia in mice and humans (40). ACAT1 expression was notably increased in THP-1-derived macrophages following the infection of Chlamydia pneumonia (41). The inhibition of ACAT1 weakened pulmonary inflammation and inhibited macrophage activation in bleomycin-induced acute lung injury (42). Among the 27 focal adhesion-related dysregulated proteins (18 up-regulated and 9 down-regulated), thrombospondin-1 (THBS1) and caveolin-1 (CAV1) have been reported to be associated with pneumonia. For example, THBS1 loss promoted the clearance of lung Klebsiella pneumonia, decreased lung inflammation burden, and enhanced the innate immune responses against Klebsiella pneumonia infection (43). CAV1 depletion reduced mouse survival rate, enhanced bacterial burdens, facilitated bacterial dissemination, and potentiated pro-inflammatory responses in mice infected with Klebsiella pneumonia (44). In other words, CAV1 enhanced the resistance of mice to Klebsiella pneumonia infection (44). Among 20 proteins related to the regulation of actin cytoskeleton, Rho guanine nucleotide exchange factor 1 (Arhgef1) and myosin light chain kinase (MYLK) have been demonstrated to be related to lung inflammation. For instance, Brown et al. demonstrated that Arhgef1 knockout mice presented decreased airway hyperreactivity and lung inflammation (45). The intravenous injection of MYLK peptide inhibitor reduced lipopolysaccharide-induced lung inflammation in mice (46). Also, 19 (9 down-regulated and 10 up-regulated) and 23 (9 down-regulated and 14 up-regulated) dysregulated proteins in the pneumonia group vs. the control group were identified to be implicated in the amoebiasis and cancer pathways, respectively. As mentioned above, the amoebiasis and cancer pathways were shared in the top 10 KEGG pathways of up-regulated and down-regulated proteins. Although amoebiasis and cancer might be irrelevant or not the study's subject, we supposed that the immunology activity for pneumonia had the same or like response activity in amoebiasis and cancer due to the central roles of pneumonia-induced dysregulated proteins in amoebiasis and cancer. Also, prior studies of some proteins related to amoebiasis and cancer preliminarily validated our speculation. For example, among amoebiasis pathway-related proteins, integrin beta-2 (ITGB2), heat shock protein beta-1 (HSPB1), LCN2, leukocyte elastase inhibitor (SERPINB1), laminin subunit alpha-4 (LAMA4), and fibronectin (FN1) have been found to be correlated with immunity. Wang et al. demonstrated that ITGB2 depletion in combination with CXCR7 and PDGFB knockdown markedly suppressed Chlamydia pneumonia entry into human cells (47). Also, ITGB2 has been identified as an immune-related gene (48, 49). HSPB1 inhibitor J2 reduced lung inflammation (50). Epinephelus coioides HSPB1 was a negative regulator in Singapore grouper iridovirus (SGIV)-induced innate immune response and apoptosis (51). LCN2 not only plays a vital role in antibacterial infection but also functions as a crucial player in the immune response to pathogenic inflammatory stimuli (52, 53). SERPINB1 loss increased the susceptibility of mice to pulmonary bacterial and viral infections (54, 55). Also, SERPINB1 controlled neutrophil survival and homeostatic expansion of IL-17+ γδ and CD4+ Th17 cells (56, 57). LAMA4 deficient mice presented impaired recruitment of neutrophils, monocytes, and lymphocytes to inflammatory loci relative to wild-type mice (58). FN1 also has been found to be involved in the regulation of innate immune response and to be correlated with immune infiltrates in cancers (59–61). Among the cancer-related proteins, signal transducer and activator of transcription 3 (STAT3) has been well documented to be inflammation and immunity (62–64). Also, STAT3 served as a positive regulator of pneumonia induced by influenza virus H1N1 (65), Agiostrongylus cantonensis (66), and Mycoplasma pneumonia (67). Combined with these data, we supposed that dysregulated proteins related to amoebiasis and cancer might also function as crucial players in pneumonia, immune, and inflammation, suggesting the same or like immune response activity between pneumonia and amoebiasis or cancer.
Both KEGG pathway annotation and enrichment analyses showed that dysregulated proteins played vital roles in pathways related to bacterial infection (e.g., bacterial invasion of epithelial cells, Staphylococcus aureus infection, and pathogenic Escherichia coli infection) and immunity (focal adhesion, phagosome, and complement and coagulation cascades), suggesting that pathways related to bacterial infection and immunity might be closely linked with the development of pneumonia in FMD. Moreover, it has been reported that bacteria including Staphylococcus aureus and Escherichia coli are common risk factors for pneumonia in FMD (11, 12, 17, 68). Moreover, multiple bacterial pathogens, such as Leclercia spp., Stenotrophononas maltophila, Staphylococcus aureus, and Staphylococcus sciur, have been identified in bovine pneumonia (69, 70). Thus, differentially expressed proteins in the pathways of Staphylococcus aureus infection, bacterial invasion of epithelial cells, and pathogenic Escherichia coli infection were screened out. After integration, a total of 32 dysregulated proteins were identified to be implicated in bacterial infection. Among these 32 proteins, 13 proteins, whose sequences were aligned onto the bovine genome, were screened out for further exploration given the close genetic relationships between FMD and bovine. Based on the PPI and node degree analyses of these 13 proteins, we supposed that 5 proteins (CTNNB1, ITGB1, CTNNA1, DNM2, and KRT19) might play crucial roles in bacteria-related pneumonia in FMD. CTNNB1 (protein name: β-catenin) and CTNNA1 are two vital players in the Wnt signaling pathway (71, 72). Wnt/β-catenin signaling has been reported to be a target of bacterial virulence factors (73) and a vital player in lung development and lung diseases (74–77). Additionally, Chen et al. demonstrated that morusin could mitigate mycoplasma pneumonia by inhibiting the Wnt/β-catenin signaling pathway in mice lung tissues (78). ITGB1, also named β1-integrin, hindered bacterial clearance and facilitated bacterial infection in cystic fibrosis airway cells and cystic fibrosis mice (79).
Given the close correlation between immune system dysfunction and pneumonia development, 85 differentially expressed proteins (53 down-regulated and 32 up-regulated) that were implicated in the immune system process were filtered out based on GO annotation analysis. Among these 85 proteins, the sequences of 49 proteins were mapped to the bovine genome. PPI and node degree analysis of these 49 proteins suggested that CTNNB1, ITGB1, ANXA5, CALR, F2, MMP9, PECAM1, THBS1, HSP90AB1, HSP90B1, ITGA3, and MSN might be the hub proteins in the pneumonia-related immune responses in FMD. Some of these proteins have been found to be implicated in pneumonia, lung inflammation, and lung injury. For instance, CALR blockade alleviated acute lung injury (ALI), reduced pro-inflammatory cytokine expression, and inhibited neutrophil and T cell infiltration in bronchoalveolar lavage and lung tissues in lipopolysaccharide (LPS)-induced ALI mouse model (80). MMP9 loss facilitated pulmonary cell death and aggravated lung injury in an interleukin-1β (IL-1β)-induced lung injury mouse model (81). MMP9 acted as a potentially protective factor against Streptococcus pneumonia infection (82, 83). PECAM1, an endothelial cell adhesion molecule, played a potentially protective role in lung injury and acute respiratory distress syndrome (84, 85). THBS1 also has been found to be implicated in the pathogenesis of gram-positive bacteria and the development of lung injury (86, 87). For instance, THBS1 loss reduced mouse survival rate, increased lung bacterial burden and lung microvascular permeability, impaired host defense against Pseudomonas aeruginosa (P. aeruginosa), and potentiated inflammatory injury during P. aeruginosa acute intrapulmonary infection (87), while P. aeruginosa is a common pathogen of pneumonia in FMD (17). HSP90B1 depletion reduced the phagocytic capacity of macrophages against Klebsiella pneumonia (K. pneumonia) (a common gram-negative bacteria that can cause pneumonia) and inhibited pro-inflammatory mediator release in alveolar and peritoneal macrophages treated with LPS derived from K. pneumonia or heat-killed K. pneumonia (88). Moreover, HSP90B1 loss in macrophages led to the increase of mouse lung K. pneumonia loads and are duction in mouse survival rate during K. pneumonia (88). CTNNA1 and CTNNB1 are two members of the catenin family (89). Multiple members of catenin family including CTNNA1 and CTNNB1 have been identified to be implicated in immune responses (90–92). For example, CTNNB1 activation enhanced the inflammatory activity of alveolar macrophages and facilitated acute host morbidity in a murine influenza pneumonia model (93). Integrins are crucial players in cell development, cell adhesion, pathogen clearance, inflammation, and immune responses (35, 36). ITGB1 and ITGA3 are two integrin family subunits (94, 95). It has been reported that ITGB1 mediated the entry of coronavirus severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) (96) and the conditional depletion of ITGB1 in type 2 alveolar epithelial cells could trigger emphysema, epithelial dysfunction, increased efferocytosis and pulmonary macrophage infiltration, and widespread lung inflammation in mice (97). Also, Li et al. suggested that ITGA3 was implicated in the infiltration of 6 immune cells (i.e., B cells, CD8 T cells, CD4 T cells, macrophages, neutrophils, and dendritic cells) in breast cancer (98).
Conclusions
Taken together, our proteomics analysis revealed that 355 proteins were differentially expressed in diseased lung tissues of FMD that died of pneumonia compared to the normal control group. KEGG annotation and enrichment analysis showed that these dysregulated proteins mainly be associated with bacterial infection and immunity. Moreover, we further screened out the dysregulated proteins related to bacterial infection (n = 32) and immunity (n = 85). Some key proteins in pneumonia-related bacterial infection and immunity were identified based on PPI and node degree analyses in the FMD. This is the first study to investigate the lung proteomics alterations caused by pneumonia in FMD, which can deepen our understanding of the molecular mechanisms of pneumonia in this rare species. Additionally, the identification of some pathways and proteins that might play vital roles in pneumonia development might contribute to the better management of pneumonia and reduction of mortality rate in FMD. However, only 6 FMD with 3 FMD in each group were used due to the rareness of this species and the difficulty in the acquisition of their organs.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material. The mass spectrometry proteomics data have been deposited via the iProX partner repository (https://www.iprox.cn/page/home.html) with the dataset identifier PXD031240.
Author contributions
JT conceived this study, designed, and supervised the experiments. LS, FL, CY, and KB performed the experiments, conducted data analysis, prepared figures and tables. JT wrote the manuscript. YW modified the manuscript. All authors reviewed and approved the manuscript.
Funding
This research was supported by Science and Technology Program of Shaanxi Academy of Science, China (Program No. 2021K-37) and Shaanxi Province Forestry Science and Technology Innovation Plan (SXLK2021-0222).
Acknowledgments
We would like to express our gratitude to all those who have helped us during the writing of this thesis. Also, we appreciate the Wuhan GeneCreate Biological Engineering Co., Ltd. for Proteomics sequencing.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2022.1012276/full#supplementary-material
References
1. Wang Y, Harris R. Moschus berezovskii. The IUCN Red List of Threatened Species. 2015:e.T13894A103431781. (2015). doi: 10.2305/IUCN.UK.2015-4.RLTS.T13894A61976926.en
2. Jie H, Zhang P, Xu Z, Mishra SK, Lei M, Zeng D, et al. microRNA and other small RNA sequence profiling across six tissues of Chinese forest musk deer (Moschus berezovskii). Biomed Res Int. (2019) 2019:4370704. doi: 10.1155/2019/4370704
3. Jie H, Zheng CL, Wang JM, Feng XL, Zeng DJ, Zhao GJ. [Research progress on molecular genetics of forest musk deer]. Zhongguo Zhong Yao Za Zhi. (2015) 40:4319–23.
4. Yang Q, Meng X, Xia L, Feng Z. Conservation status and causes of decline of musk deer (Moschus spp.) in China. Biol Conserv. (2003) 109:333–42. doi: 10.1016/S0006-3207(02)00159-3
5. Lv S, Lei Z, Yan G, Shah SA, Ahmed S, Sun T. Chemical compositions and pharmacological activities of natural musk (Moschus) and artificial musk: a review. J Ethnopharmacol. (2022) 284:114799. doi: 10.1016/j.jep.2021.114799
6. Liu K, Xie L, Deng M, Zhang X, Luo J, Li X. Zoology, chemical composition, pharmacology, quality control and future perspective of Musk (Moschus): a review. Chin Med. (2021) 16:46. doi: 10.1186/s13020-021-00457-8
7. Wang J, Xing H, Qin X, Ren Q, Yang J, Li L. Pharmacological effects and mechanisms of muscone. J Ethnopharmacol. (2020) 262:113120. doi: 10.1016/j.jep.2020.113120
8. Guan T-L, Zeng B, Peng Q-K, Yue B-S, Zou F-D. Microsatellite analysis of the genetic structure of captive forest musk deer populations and its implication for conservation. Biochem Syst Ecol. (2009) 37:166–73. doi: 10.1016/j.bse.2009.04.001
9. Xiuxiang M, Caiquan Z, Jinchu H, Cao L, Zhibin M, Jinchao F, et al. Musk deer farming in China. Anim Sci. (2006) 82:1–6. doi: 10.1079/ASC200516
10. He L. Li L-h, Wang W-x, Liu G, Liu S-q, Liu W-h, et al. Welfare of farmed musk deer: changes in the biological characteristics of musk deer in farming environments. Appl Anim Behav Sci. (2014) 156:1–5. doi: 10.1016/j.applanim.2014.03.011
11. Qiao J, Wu X, Su L. A review of mainly affected on musk-deer diseases: purulent, respiratory system and parasitic diseases. J Econ Anim. (2009) 13:104–7.
12. Tian Q, Zhou X, Cheng J, Luo Y, Dai L, Zhao W, et al. Genome sequence of lung pathogenic Escherichia coli O78, a chimeric strain isolated from pneumonia forest musk deer. Genes Genomics. (2017) 39:805–15. doi: 10.1007/s13258-017-0545-4
13. Zhao W, Yu D, Cheng JG, Luo Y, Wang Y, Yao XP, et al. Affinity of Streptococcal G protein to forest musk deer (Moschus berezovskii) serum immunoglobulin G. J Wildl Dis. (2020) 56:684–6. doi: 10.7589/2019-09-223
14. Yan M, Yan Q, Yang G. The mass diseases of captive musk deer. J Econ Anim. (2016) 20:112–7. doi: 10.13326/j.jea.2016.112
15. Luo Y, Kang J, Cheng J, Zou L, Li B, Wang C, et al. Physicochemical properties of musk deer pneumonia and purulent disease viruses. Animal Husbandry and Feed Science. (2010) 2:37–40. doi: 10.19578/j.cnki.ahfs.2010.05.014
16. Karikalan M, Ram H, Pathak S, Banerjee PS, Chandra MS, Pawdeand AM, et al. Pneumonia due to protostrongylid lung worms in wild himalayan bluesheep (Pseudoisnayaur). Indian J Vet Pathol. (2020) 44:181–3. doi: 10.5958/0973-970X.2020.00036.X
17. Zhao W, Ren Z, Luo Y. Metagenomics analysis of the gut microbiome in healthy and bacterial pneumonia forest musk deer. Genes Genomics. (2021) 43:43–53. doi: 10.1007/s13258-020-01029-0
18. Zhao W, Tian Q, Luo Y, Wang Y, Yang ZX, Yao XP, et al. Isolation, identification, and genome analysis of lung pathogenic klebsiella pneumoniae (LPKP) in forest musk deer. J Zoo Wildl Med. (2017) 48:1039–48. doi: 10.1638/2016-0241.1
19. McArdle AJ, Menikou S. What is proteomics? Arch Dis Child Educ Pract Ed. (2021) 106:178–81. doi: 10.1136/archdischild-2019-317434
20. Noor Z, Ahn SB, Baker MS, Ranganathan S, Mohamedali A. Mass spectrometry-based protein identification in proteomics-a review. Brief Bioinform. (2021) 22:1620–38. doi: 10.1093/bib/bbz163
21. Evans C, Noirel J, Ow SY, Salim M, Pereira-Medrano AG, Couto N, et al. An insight into iTRAQ: where do we stand now? Anal Bioanal Chem. (2012) 404:1011–27. doi: 10.1007/s00216-012-5918-6
22. Rauniyar N, Yates JR 3rd. Isobaric labeling-based relative quantification in shotgun proteomics. J Proteome Res. (2014) 13:5293–309. doi: 10.1021/pr500880b
23. Chen X, Sun Y, Zhang T, Shu L, Roepstorff P, Yang F. Quantitative proteomics using isobaric labeling: a practical guide. Genomics Proteomics Bioinformatics. (2022) 19:689–706. doi: 10.1016/j.gpb.2021.08.012
24. Li Y, Guo T, Wang X, Ni W, Hu R, Cui Y, et al. ITRAQ-based quantitative proteomics reveals the proteome profiles of MDBK cells infected with bovine viral diarrhea virus. Virol J. (2021) 18:119. doi: 10.1186/s12985-021-01592-2
25. Miller I, de Almeida AM, Eckersall PD. Across the great divide: proteomics becoming an essential tool for animal and veterinary sciences. J Proteomics. (2021) 241:104225. doi: 10.1016/j.jprot.2021.104225
26. Zhang Z, Yu J, Wang P, Lin L, Liu R, Zeng R, et al. iTRAQ-based proteomic profiling reveals protein alterations after traumatic brain injury and supports thyroxine as a potential treatment. Mol Brain. (2021) 14:25. doi: 10.1186/s13041-021-00739-0
27. Bian J, Liang M, Ding S, Wang L, Ni W, Xiong S, et al. iTRAQ-based high-throughput proteomics analysis reveals alterations of plasma proteins in patients infected with human bocavirus. PLoS ONE. (2019) 14:e0225261. doi: 10.1371/journal.pone.0225261
28. Tang J, Yu Y, Zheng H, Yin L, Sun M, Wang W, et al. ITRAQ-based quantitative proteomic analysis of Cynops orientalis limb regeneration. BMC Genomics. (2017) 18:750. doi: 10.1186/s12864-017-4125-4
29. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. (2019) 47:D607–13. doi: 10.1093/nar/gky1131
30. Szklarczyk D, Gable AL, Nastou KC, Lyon D, Kirsch R, Pyysalo S, et al. The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. (2021) 49:D605–d612. doi: 10.1093/nar/gkaa1074
31. Fan Z, Li W, Jin J, Cui K, Yan C, Peng C, et al. The draft genome sequence of forest musk deer (Moschus berezovskii). GigaScience. (2018) 7:giy038. doi: 10.1093/gigascience/giy038
32. Xu Z, Jie H, Chen B, Gaur U, Wu N, Gao J, et al. Illumina-based de novo transcriptome sequencing and analysis of Chinese forest musk deer. J Genet. (2017) 96:1033–40. doi: 10.1007/s12041-017-0872-x
33. Fox AD, Hescott BJ, Blumer AC, Slonim DK. Connectedness of PPI network neighborhoods identifies regulatory hub proteins. Bioinformatics. (2011) 27:1135–42. doi: 10.1093/bioinformatics/btr099
34. Vallabhajosyula RR, Chakravarti D, Lutfeali S, Ray A, Raval A. Identifying hubs in protein interaction networks. PLoS ONE. (2009) 4:e5344. doi: 10.1371/journal.pone.0005344
35. Quinton LJ, Walkey AJ, Mizgerd JP. Integrative physiology of pneumonia. Physiol Rev. (2018) 98:1417–64. doi: 10.1152/physrev.00032.2017
36. Mizgerd JP. Pathogenesis of severe pneumonia: advances and knowledge gaps. Curr Opin Pulm Med. (2017) 23:193–7. doi: 10.1097/MCP.0000000000000365
37. Wygrecka M, Marsh LM, Morty RE, Henneke I, Guenther A, Lohmeyer J, et al. Enolase-1 promotes plasminogen-mediated recruitment of monocytes to the acutely inflamed lung. Blood. (2009) 113:5588–98. doi: 10.1182/blood-2008-08-170837
38. Wu H, Santoni-Rugiu E, Ralfkiaer E, Porse BT, Moser C, Høiby N, et al. Lipocalin 2 is protective against E. coli pneumonia. Respir Res. (2010) 11:96. doi: 10.1186/1465-9921-11-96
39. Sheldon JR, Himmel LE, Kunkle DE, Monteith AJ, Maloney KN, Skaar EP. Lipocalin-2 is an essential component of the innate immune response to Acinetobacter baumannii infection. PLoS Pathog. (2022) 18:e1010809. doi: 10.1371/journal.ppat.1010809
40. Warszawska JM, Gawish R, Sharif O, Sigel S, Doninger B, Lakovits K, et al. Lipocalin 2 deactivates macrophages and worsens pneumococcal pneumonia outcomes. J Clin Invest. (2013) 123:3363–72. doi: 10.1172/JCI67911
41. Liu W, He P, Cheng B, Mei CL, Wang YF, Wan JJ. Chlamydia pneumoniae disturbs cholesterol homeostasis in human THP-1 macrophages via JNK-PPARγ dependent signal transduction pathways. Microbes Infect. (2010) 12:1226–35. doi: 10.1016/j.micinf.2010.09.004
42. Stevenson ER, Wilkinson ML, Abramova E, Guo C, Gow AJ. Intratracheal administration of Acyl Coenzyme A Acyltransferase-1 inhibitor K-604 reduces pulmonary inflammation following bleomycin-induced lung injury. J Pharmacol Exp Ther. (2022) 382:356–65. doi: 10.1124/jpet.122.001284
43. Zhao Y, Olonisakin TF, Xiong Z, Hulver M, Sayeed S, Yu MT, et al. Thrombospondin-1 restrains neutrophil granule serine protease function and regulates the innate immune response during Klebsiella pneumoniae infection. Mucosal Immunol. (2015) 8:896–905. doi: 10.1038/mi.2014.120
44. Guo Q, Shen N, Yuan K, Li J, Wu H, Zeng Y 3rd, Bansal AK, Singh BB, Gao H Caveolin-1 plays a critical role in host immunity against Klebsiella pneumoniae by regulating STAT5 and Akt activity. Eur J Immunol. (2012) 42:1500–11. doi: 10.1002/eji.201142051
45. Brown JP, Taube C, Miyahara N, Koya T, Pelanda R, Gelfand EW, et al. Arhgef1 is required by T cells for the development of airway hyperreactivity and inflammation. Am J Respir Crit Care Med. (2007) 176:10–9. doi: 10.1164/rccm.200702-270OC
46. Mirzapoiazova T, Moitra J, Moreno-Vinasco L, Sammani S, Turner JR, Chiang ET, et al. Non-muscle myosin light chain kinase isoform is a viable molecular target in acute inflammatory lung injury. Am J Respir Cell Mol Biol. (2011) 44:40–52. doi: 10.1165/rcmb.2009-0197OC
47. Wang A, Johnston SC, Chou J, Dean D. A systemic network for Chlamydia pneumoniae entry into human cells. J Bacteriol. (2010) 192:2809–15. doi: 10.1128/JB.01462-09
48. Wei J, Huang XJ, Huang Y, Xiong MY, Yao XY, Huang ZN Li SN, et al. Key immune-related gene ITGB2 as a prognostic signature for acute myeloid leukemia. Ann Transl Med. (2021) 9:1386. doi: 10.21037/atm-21-3641
49. Altorki T, Muller W, Brass A, Cruickshank S. The role of β(2) integrin in dendritic cell migration during infection. BMC Immunol. (2021) 22:2. doi: 10.1186/s12865-020-00394-5
50. Oh A, Jeon S, Jeong MG, Kim HK, Kang J, Lee YS, et al. HSPB1 inhibitor J2 attenuates lung inflammation through direct modulation of Ym1 production and paracrine signaling. Biomed Pharmacother. (2021) 143:112225. doi: 10.1016/j.biopha.2021.112225
51. Li PH, Cai YJ, Zhu XL, Yang JD, Yang SQ, Huang W, et al. Epinephelus coioides Hsp27 negatively regulates innate immune response and apoptosis induced by Singapore grouper iridovirus (SGIV) infection. Fish Shellfish Immunol. (2022) 120:470–80. doi: 10.1016/j.fsi.2021.12.016
52. Guardado S, Ojeda-Juárez D, Kaul M, Nordgren TM. Comprehensive review of lipocalin 2-mediated effects in lung inflammation. Am J Physiol Lung Cell Mol Physiol. (2021) 321:L726–l733. doi: 10.1152/ajplung.00080.2021
53. Xiao X, Yeoh BS, Vijay-Kumar M. Lipocalin 2: An Emerging Player in Iron Homeostasis and Inflammation. Annu Rev Nutr. (2017) 37:103–30. doi: 10.1146/annurev-nutr-071816-064559
54. Benarafa C, Priebe GP, Remold-O'Donnell E. The neutrophil serine protease inhibitor serpinb1 preserves lung defense functions in Pseudomonas aeruginosa infection. J Exp Med. (2007) 204:1901–9. doi: 10.1084/jem.20070494
55. Farley K, Stolley JM, Zhao P, Cooley J, Remold-O'Donnell E. A serpinB1 regulatory mechanism is essential for restricting neutrophil extracellular trap generation. J Immunol. (2012) 189:4574–81. doi: 10.4049/jimmunol.1201167
56. Zhao P, Hou L, Farley K, Sundrud MS, Remold-O'Donnell E. SerpinB1 regulates homeostatic expansion of IL-17+ γδ and CD4+ Th17 cells. J Leukoc Biol. (2014) 95:521–30. doi: 10.1189/jlb.0613331
57. Benarafa C, LeCuyer TE, Baumann M, Stolley JM, Cremona TP. Remold-O'Donnell E: SerpinB1 protects the mature neutrophil reserve in the bone marrow. J Leukoc Biol. (2011) 90:21–9. doi: 10.1189/jlb.0810461
58. Kenne E, Soehnlein O, Genové G, Rotzius P, Eriksson EE, Lindbom L. Immune cell recruitment to inflammatory loci is impaired in mice deficient in basement membrane protein laminin alpha4. J Leukoc Biol. (2010) 88:523–8. doi: 10.1189/jlb.0110043
59. Zheng M, Ambesi A, McKeown-Longo PJ. Role of TLR4 receptor complex in the regulation of the innate immune response by fibronectin. Cells. (2020) 9:216. doi: 10.3390/cells9010216
60. Kelsh R, You R, Horzempa C, Zheng M, McKeown-Longo PJ. Regulation of the innate immune response by fibronectin: synergism between the III-1 and EDA domains. PLoS ONE. (2014) 9:e102974. doi: 10.1371/journal.pone.0102974
61. Zhang XX, Luo JH, Wu LQ. FN1 overexpression is correlated with unfavorable prognosis and immune infiltrates in breast cancer. Front Genet. (2022) 13:913659. doi: 10.3389/fgene.2022.913659
62. Wang Y, Shen Y, Wang S, Shen Q, Zhou X. The role of STAT3 in leading the crosstalk between human cancers and the immune system. Cancer Lett. (2018) 415:117–28. doi: 10.1016/j.canlet.2017.12.003
63. Chen W, Dai X, Chen Y, Tian F, Zhang Y, Zhang Q, et al. Significance of STAT3 in immune infiltration and drug response in cancer. Biomolecules. (2020) 10:834. doi: 10.3390/biom10060834
64. Hu YS, Han X, Liu XH. STAT3: a potential drug target for tumor and inflammation. Curr Top Med Chem. (2019) 19:1305–17. doi: 10.2174/1568026619666190620145052
65. Guo L, Wang Q, Zhang D. MicroRNA-4485 ameliorates severe influenza pneumonia via inhibition of the STAT3/PI3K/AKT signaling pathway. Oncol Lett. (2020) 20:215. doi: 10.3892/ol.2020.12078
66. Zhou H, Lu Y, Wei H, Chen Y, Limpanon Y, Dekumyoy P, et al. Stat3/IL-6 signaling mediates sustained pneumonia induced by Agiostrongylus cantonensis. PLoS Negl Trop Dis. (2022) 16:e0010461. doi: 10.1371/journal.pntd.0010461
67. Chen Y, Dong S, Tian L, Chen H, Chen J, He C. Combination of azithromycin and methylprednisolone alleviates Mycoplasma pneumoniae induced pneumonia by regulating miR-499a-5p/STAT3 axis. Exp Ther Med. (2022) 24:578. doi: 10.3892/etm.2022.11515
68. Shaukat A, Hanif S, Shaukat I, Shukat R, Rajput SA, Jiang K, et al. Upregulated-gene expression of pro-inflammatory cytokines, oxidative stress and apoptotic markers through inflammatory, oxidative and apoptosis mediated signaling pathways in Bovine Pneumonia. Microb Pathog. (2021) 155:104935. doi: 10.1016/j.micpath.2021.104935
69. Choudhary M, Choudhary BK, Chandra Ghosh R, Bhoyar S, Chaudhari S, Barbuddhe SB. Cultivable microbiota and pulmonary lesions in polymicrobial bovine pneumonia. Microb Pathog. (2019) 134:103577. doi: 10.1016/j.micpath.2019.103577
70. Panciera RJ, Confer AW. Pathogenesis and pathology of bovine pneumonia. Vet Clin North Am Food Anim Pract. (2010) 26:191–214. doi: 10.1016/j.cvfa.2010.04.001
71. van der Wal T, van Amerongen R. Walking the tight wire between cell adhesion and WNT signalling: a balancing act for β-catenin. Open Biol. (2020) 10:200267. doi: 10.1098/rsob.200267
72. Jesinghaus M, Konukiewitz B, Foersch S, Stenzinger A, Steiger K. Appendiceal goblet cell carcinoids and adenocarcinomas ex-goblet cell carcinoid are genetically distinct from primary colorectal-type adenocarcinoma of the appendix. Mod Pathol. (2018) 31:829–39. doi: 10.1038/modpathol.2017.184
73. Silva-García O, Valdez-Alarcón JJ, Baizabal-Aguirre VM. Wnt/β-catenin signaling as a molecular target by pathogenic bacteria. Front Immunol. (2019) 10:2135. doi: 10.3389/fimmu.2019.02135
74. Hussain M, Xu C, Lu M, Wu X, Tang L, Wu X. Wnt/β-catenin signaling links embryonic lung development and asthmatic airway remodeling. Biochim Biophys Acta Mol Basis Dis. (2017) 1863:3226–42. doi: 10.1016/j.bbadis.2017.08.031
75. Raslan AA, Yoon JK. WNT signaling in lung repair and regeneration. Mol Cells. (2020) 43:774–83. doi: 10.14348/molcells.2020.0059
76. Baarsma HA, Königshoff M. 'WNT-er is coming': WNT signalling in chronic lung diseases. Thorax. (2017) 72:746–59. doi: 10.1136/thoraxjnl-2016-209753
77. Königshoff M, Eickelberg O. WNT signaling in lung disease: a failure or a regeneration signal? Am J Respir Cell Mol Biol. (2010) 42:21–31. doi: 10.1165/rcmb.2008-0485TR
78. Chen C, Wang J, Chen J, Zhou L, Wang H, Chen J, et al. Morusin alleviates mycoplasma pneumonia via the inhibition of Wnt/β-catenin and NF-κB signaling. Biosci Rep. (2019) 39:BSR20190190. doi: 10.1042/BSR20190190
79. Grassmé H, Henry B, Ziobro R, Becker KA, Riethmüller J, Gardner A, et al. β1-integrin accumulates in cystic fibrosis luminal airway epithelial membranes and decreases sphingosine, promoting bacterial infections. Cell Host Microbe. (2017) 21:707–18.e708. doi: 10.1016/j.chom.2017.05.001
80. Jiang Z, Chen Z, Hu L, Qiu L, Zhu L. Calreticulin blockade attenuates murine acute lung injury by inducing polarization of M2 subtype macrophages. Front Immunol. (2020) 11:11. doi: 10.3389/fimmu.2020.00011
81. Lukkarinen H, Hogmalm A, Lappalainen U, Bry K. Matrix metalloproteinase-9 deficiency worsens lung injury in a model of bronchopulmonary dysplasia. Am J Respir Cell Mol Biol. (2009) 41:59–68. doi: 10.1165/rcmb.2008-0179OC
82. Böttcher T, Spreer A, Azeh I, Nau R, Gerber J. Matrix metalloproteinase-9 deficiency impairs host defense mechanisms against Streptococcus pneumoniae in a mouse model of bacterial meningitis. Neurosci Lett. (2003) 338:201–4. doi: 10.1016/S0304-3940(02)01406-4
83. Hong JS, Greenlee KJ, Pitchumani R, Lee SH, Song LZ, Shan M, et al. Dual protective mechanisms of matrix metalloproteinases 2 and 9 in immune defense against Streptococcus pneumoniae. J Immunol. (2011) 186:6427–36. doi: 10.4049/jimmunol.1003449
84. Villar J, Zhang H, Slutsky AS. Lung repair and regeneration in ARDS: role of PECAM1 and Wnt signaling. Chest. (2019) 155:587–94. doi: 10.1016/j.chest.2018.10.022
85. Shi J, Hu CL, Gao YF, Liao XX, Xu H. The relationship between platelet endothelial cell adhesion molecule-1 and paraquat-induced lung injury in rabbits. World J Emerg Med. (2012) 3:60–4. doi: 10.5847/wjem.j.issn.1920-8642.2012.01.011
86. Binsker U, Kohler TP, Hammerschmidt S. Contribution of human Thrombospondin-1 to the pathogenesis of gram-positive bacteria. J Innate Immun. (2019) 11:303–15. doi: 10.1159/000496033
87. Qu Y, Olonisakin T, Bain W, Zupetic J, Brown R, Hulver M, et al. Thrombospondin-1 protects against pathogen-induced lung injury by limiting extracellular matrix proteolysis. JCI Insight. (2018) 3:e96914. doi: 10.1172/jci.insight.96914
88. Anas AA, de Vos AF, Hoogendijk AJ, van Lieshout MH, van Heijst JW, Florquin S, et al. Endoplasmic reticulum chaperone gp96 in macrophages is essential for protective immunity during Gram-negative pneumonia. J Pathol. (2016) 238:74–84. doi: 10.1002/path.4637
89. Smith A, Bourdeau I, Wang J, Bondy CA. Expression of Catenin family members CTNNA1, CTNNA2, CTNNB1 and JUP in the primate prefrontal cortex and hippocampus. Brain Res Mol Brain Res. (2005) 135:225–31. doi: 10.1016/j.molbrainres.2004.12.025
90. Chignalia AZ, Vogel SM, Reynolds AB, Mehta D, Dull RO, Minshall RD, et al. p120-catenin expressed in alveolar type II cells is essential for the regulation of lung innate immune response. Am J Pathol. (2015) 185:1251–63. doi: 10.1016/j.ajpath.2015.01.022
91. Pai SG, Carneiro BA, Mota JM, Costa R, Leite CA, Barroso-Sousa R, et al. Wnt/beta-catenin pathway: modulating anticancer immune response. J Hematol Oncol. (2017) 10:101. doi: 10.1186/s13045-017-0471-6
92. Tong J, Ji X, Zhang H, Xiong B, Cui D, Jiang L. The analysis of the ubiquitylomic responses to streptococcus agalactiae infection in bovine mammary gland epithelial cells. J Inflamm Res. (2022) 15:4331–43. doi: 10.2147/JIR.S368779
93. Zhu B, Wu Y, Huang S, Zhang R, Son YM, Li C, et al. Uncoupling of macrophage inflammation from self-renewal modulates host recovery from respiratory viral infection. Immunity. (2021) 54:1200–18.e1209. doi: 10.1016/j.immuni.2021.04.001
94. Liu M, Zhang Y, Yang J, Cui X, Zhou Z, Zhan H, et al. ZIP4 Increases expression of transcription factor ZEB1 to promote integrin α3β1 signaling and inhibit expression of the gemcitabine transporter ENT1 in pancreatic cancer cells. Gastroenterology. (2020) 158:679–92.e671. doi: 10.1053/j.gastro.2019.10.038
95. Mautone L, Ferravante C, Tortora A, Tarallo R, Giurato G, Weisz A, et al. Higher integrin Alpha 3 beta1 expression in papillary thyroid cancer is associated with worst outcome. Cancers. (2021) 13:2937. doi: 10.3390/cancers13122937
96. Zhang Y, Shang L, Zhang J, Liu Y, Jin C, Zhao Y, et al. An antibody-based proximity labeling map reveals mechanisms of SARS-CoV-2 inhibition of antiviral immunity. Cell Chem Biol. (2022) 29:5–18.e16. doi: 10.1016/j.chembiol.2021.10.008
97. Plosa EJ, Benjamin JT, Sucre JM, Gulleman PM, Gleaves LA, Han W, et al. β1 Integrin regulates adult lung alveolar epithelial cell inflammation. JCI Insight. (2020) 5:129259. doi: 10.1172/jci.insight.129259
Keywords: forest musk deer, pneumonia, proteomics, bacterium, immunity
Citation: Tang J, Suo L, Li F, Yang C, Bian K and Wang Y (2022) ITRAQ-based quantitative proteomics analysis of forest musk deer with pneumonia. Front. Vet. Sci. 9:1012276. doi: 10.3389/fvets.2022.1012276
Received: 05 August 2022; Accepted: 05 October 2022;
Published: 26 October 2022.
Edited by:
Nelida Rodriguez-Osorio, Universidad de la República, UruguayReviewed by:
Mauricio Corredor, University of Antioquia, ColombiaWuzi Dong, Northwest A&F University, China
Copyright © 2022 Tang, Suo, Li, Yang, Bian and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jie Tang, eWF5YTE4NCYjeDAwMDQwO3hhYi5hYy5jbg==