Skip to main content

Transcriptome profiling of lung immune responses potentially related to acute respiratory distress syndrome in forest musk deer



Forest musk deer is an endangered species globally. The death of captive forest musk deer can be caused by certain respiratory system diseases. Acute respiratory distress syndrome (ARDS) is a huge threat to the life of forest muck deer that breed in our department.


Lung histopathologic analysis was conducted by hematoxylin and eosin (HE) staining. The lung gene changes triggered by ARDS were examined by RNA sequencing and related bioinformatics analysis in forest musk deer. The potential functions of unigenes were investigated by NR, SwissProt KOG, GO, and KEGG annotation analyses. Vital biological processes or pathways in ARDS were examined by GO and KEGG enrichment analyses.


A total of 3265 unigenes were differentially expressed (|log2fold-change|> 2 and adjusted P value < 0.01) in lung tissues of 3 forest musk deer with ARDS compared with normal lung tissues of the non-ARDS group. These differentially expressed unigenes (DEGs) played crucial roles in immunity and defense responses to pathogens. Moreover, we identified the DEGs related to one or more of the following biological processes: lung development, immunity, and bacterial/viral/fungal infection. And six DEGs that might be involved in lung injury caused by immune dysregulation or viral/fungal infection were identified.


ARDS-mediated lung gene alterations were identified in forest musk deer. Moreover, multiple genes involved in lung development and lung defense responses to bacteria/viruses/fungi in ARDS were filtered out in forest musk deer.

Peer Review reports


Forest musk deer (Moschus berezovskii), mainly distributed in China and Vietnam, are listed as “Critically Endangered under criterion A2cd” in the International Union for Conservation of Nature (IUCN) RED LIST CATEGORY AND CRITERIA (version 3.1) [1]. Moreover, forest musk deer are classified as a first-class nationally-protected wildlife species under the Chinese Wild Animal Protection Law [2, 3]. Forest musk deer can secrete the musk, which has multiple potential pharmacological properties such as neuroprotective, antioxidative, anti-inflammatory, and anticancer effects in the traditional Asian medicine industry [4, 5]. Moreover, musk has been widely used in the international perfume industry [6, 7]. Thus, many forest musk deer were hunted and trapped in the past, which resulted in a steep decline in the population of forest musk deer.

To relieve the poaching and resource pressure on wild forest musk deer, the Chinese government has encouraged enterprises to establish artificial breeding systems for forest musk deer since the 1950s [8]. Despite much progress in the artificial breeding of forest musk deer in China, the long-term development of forest musk deer breeding is hindered by multiple diseases, including respiratory disorders [9,10,11]. In the animal farming process, acute respiratory distress syndrome (ARDS) is a major threat to the life of captive forest musk deer. ARDS, a form of acute respiratory failure with a higher mortality rate, is closely linked with multiple diseases including acute lung injury, sepsis, and pneumonia [12,13,14]. Moreover, it has been reported that infections with multiple pathogens (e.g., bacteria, viruses, fungi, and parasites) and dysfunctions of the immune system are involved in the development of pneumonia, acute lung injury, and ARDS [14,15,16]. In-depth insight into the molecular mechanisms underlying ARDS might deepen our understanding of the pathogenesis of ARDS and contribute to the better management of ARDS.

In this project, ARDS-induced lung gene alterations were examined by RNA sequencing (RNA-seq) in forest musk deer. Moreover, the major biological processes or signaling pathways regulated by these differentially expressed unigenes (DEGs) were identified. Additionally, multiple DEGs that might play vital roles in lung development, immunity, and infections of pathogens (e.g., bacteria, viruses, and fungi) were screened out.

Materials and methods

Animals and sample collection

The animal experiments were performed with the approval of the Ethics Committee of Shaanxi Institute of Zoology. Forest musk deer were raised in the Shaanxi Institute of Zoology (Xi an, China). Lung samples were collected from 3 adult forest musk deer (2♂1♀, 4.5 years old) with ARDS that died within 3 days of the onset of respiratory symptoms and 3 adult forest musk deer(♂, 3.5 years old) that died from fighting against each other. The middle lobe samples of the right lungs were collected immediately after death in the control group. Two anterior lobe samples and one middle lobe sample of the right lungs were obtained within 4 h after natural death in the ARDS group. Lung tissue samples were stored at -80 °C.

Histopathologic analysis

The histopathologic alterations of lung tissues were examined by hematoxylin and eosin (HE) staining. Briefly, paraffin-embedded lung tissue Sects. (4 μm) were prepared after fixation, dehydration, and vitrification. Next, the lung tissue sections were treated and stained using the Hematoxylin and Eosin Staining Kit (Beyotime, Shanghai, China) following the protocols of the manufacturer.

RNA extraction and library preparation

Total RNA was extracted from diseased lung tissues (Pne_1, Pne_2, and Pne_3) (ARDS group) and normal lung tissues (Ctrl_1, Ctrl_2, and Ctrl_3) (no-ARDS group) using the Trizol reagent (Thermo Scientific, Waltham, MA, USA) according to the protocols of the manufacturer. RNA concentration and purity were examined by NanoDrop 2000 spectrophotometer (Thermo Scientific). RNA purity and integrity were tested by gel electrophoresis on a 1% agarose gel. Also, the concentration and integrity of RNA were further analyzed using the Agilent RNA Nano 6000 Assay Kit on Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). Next, high-quality RNA (RNA Integrity number > 7.0; 1 µg/sample) was used in the RNA sample preparations. The cDNA libraries were constructed using the NEBNextUltra RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA) following the manufacturer’s instructions and index codes were added. Briefly, mRNA with poly-A structure was enriched using poly-T oligo-attached magnetic beads and then fragmented using divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5X). Then, mRNA was reversely transcribed into first-strand cDNA using M-MuLV Reverse Transcriptase (RNase H) and random hexamer primer. Next, second-strand cDNA was synthesized using DNA Polymerase I and RNase H. Subsequently, remaining overhangs of cDNA were converted into blunt ends and NEBNext Adaptor with hairpin loop structure were ligated with adenylated DNAs (3’ends). The cDNA library fragments (250 ~ 300 bp) were purified with the AMPure XP system (Beckman Coulter, Brea, CA, USA). Next, cDNA was co-incubated with USER Enzyme (NEB) for 15 min at 37 °C. After enzyme inactivation for 5 min at 95 °C, the PCR reaction was performed using Phusion High-Fidelity DNA polymerase, Universal PCR primers, and Index Primer. Next, PCR products were purified using the AMPure XP system (Beckman Coulter), and the quality of the cDNA library was evaluated on the Agilent Bioanalyzer 2100 system (Agilent Technologies). The index-coded samples were clustered on a cBot Cluster Generation System using TruSeq PE Cluster Kit v3-cBot-HS (Illumia) according to the manufacturer’s instructions. After cluster analysis, the libraries were sequenced on the Illumina Novaseq platform with paired-end 150 bp reads.

Data analysis

The quality of raw data was evaluated using the FASTQCv0.10.1 software (parameters: extract/f fastq/t 18). Reads with less than or equal to 2 mismatches against the primers/adaptors and low-quality reads with more than 10% non-ATCG (unknown) bases or 50% of bases with the Q score ≤ 5 were removed using the NGS QC TOOLKIT v2.3.3software as previously described [17,18,19]. Also, the statistical information of raw reads/bases, clean reads/bases, valid bases (clean reads/raw reads *100%), Q20, Q30, and GC content was collected. The reference genome sequences of forest musk deer were downloaded from the ENA database ( Next, clean reads after quality control were aligned to the reference genome of forest musk deer using the Hisat2v2.1.0 software(parameters: fr/dta/ufflinks/S) and assembled into unigenes using the StringTiev1.3.3b software (default parameters).

Unigene annotation

The annotation analysis of unigenes was performed by comparison against the NCBI non-redundant (NR), SwissProt, clusters of orthologous groups for eukaryotic complete genomes (KOG), gene ontology (GO), and Kyoto encyclopedia of genes and genomes (KEGG) databases using the Basic Local Alignment Search Tool (BLAST) algorithm [20]. Forest musk deer unigenes with significant sequence homology to the genes/proteins in these databases (E value < 1E-5) were obtained and the annotation information was retrieved. GO and KEGG annotation analysis was performed using the Blast2GO software or KAAS online tool, respectively. Blast2GO and KAAS can be used for the automatic functional annotation of DNA or protein sequence data in non-model species with no available annotation information [21,22,23].

Differential expression, enrichment, and Venn analysis

Differential expression analysis was performed using the DESeq2 R package. The ARDS versus non-ARDS model was used in the DEseq2 analysis. Unigenes with |log2fold change|> 2 and adjusted P value < 0.01 were identified to be differentially expressed. Additionally, GO and KEGG enrichment analyses were conducted and the significance was examined by a hypergeometric distribution test. Terms with a P value < 0.05 were considered to be significantly enriched. Venn analysis was used to compare lists by jvenn online tool (, which can display the specific or shared elements of different lists.


Anatomical and histopathologic analysis of lung tissues of forest musk deer

As presented in Fig. 1A-C, the lung tissues of forest musk deer with ARDS presented notable clinical-pathological alterations such as purulence and hemorrhage. HE staining of normal lung tissue shows normal alveoli and alveolar walls in Fig. 1D. Histopathologic analysis showed that there were copious pus cells (Fig. 1E-G) and exudates (Fig. 1E) in the alveolar cavity and noticeable congestion in the pulmonary interstitium (Fig. 1G) and alveolar wall (Fig. 1E) in the lung tissues of forest musk deer with ARDS. Moreover, the bronchi were filled with pus cells in the ARDS group (Fig. 1G).

Fig. 1
figure 1

Anatomical and histopathologic analysis of forest musk deer lung tissues. (A-C) The clinical-pathological alterations of lung tissues of 3 forest musk deer that died of ARDS. (D) HE staining image of normal lung tissue. Black arrow: alveoli; Red arrow: alveolar wall. (E-G) HE staining images of lung tissues of forest musk deer with ARDS. The black arrow in Fig. 1E-1G represents the pus cells in the alveolar space. The red arrow in Fig. 1E indicates the exudates in the alveolar space. The red arrow in Fig. 1F denotes the alveolar wall capillary congestion. The red arrow in Fig. 1G indicates pulmonary interstitial congestion. The yellow arrow in Fig. 1G indicates a large number of pus cells in the bronchi

Illumina sequencing and assembly

To identify genes that might play crucial roles in ARDS-induced lung injury, the lung gene changes triggered by ARDS were investigated by RNA-seq in forest musk deer. As shown in Table S1, the Q30 values (percentage of bases with the sequencing error rate of less than 0.1%) were more than 94%, and the GC percentages were approximately 50%, suggesting the production of high-quality sequencing data. Next, the clean reads were assembled de novo into unigenes using StringTie software. The length and GC content distribution patterns of unigenes are presented in Fig. 2A and B.

Fig. 2
figure 2

Statistical analysis of unigenes. (A) Length distribution patterns ofunigenes. (B) GC content frequency distribution patterns of unigenes

Functional annotation analysis

Although the reference genome of forest musk deer is currently available, the functional annotation information is relatively lacking. To obtain in-depth insight into the functions of the forest musk deer genes, the unigenes were annotated using the NCBI NR, SWISSPROT, KOG, KEGG, and GO databases. The results showed that a total of 40,851 unigenes were annotated by NR, SWISSPROT, KOG, KEGG, and GO databases. Among these unigenes, 37,056 (90.71%), 30,950 (75.76%), 24,836 (60.80%), 11,740 (28.74%), and27679 (67.76%) were mapped to the NR, SWISSPROT, KOG, KEGG, and GO databases, respectively (Table S2). NCBI NR alignment results showed that 19.07% of forest musk deer unigenes matched the genome sequences of Bostaurus (Fig. 3A). Additionally, 12.23%, 11%, 9.72%, and9.09% of the BLAST hits were identified within the reference protein databases of Bubalus bubalis, Capra hircus, Bosmutus, and Ovis arie smusimon, respectively (Fig. 3A). GO annotation analysis revealed that these unigenes were involved in the regulation of multiple essential biological processes such as immune system processes, cellular component organization or biogenesis, response to stimulus, growth, and cell killing (Fig. 3B). KEGG annotation analysis showed that these unigeness were divided into 4 categories (i.e. environmental information processing, cellular processes, metabolism, and genetic information processing) (Fig. 3C).

Fig. 3
figure 3

Species and functional annotations of unigenes. (A) Species distribution of unigenes annotated by the NR database. (B) GO annotation analysis for unigenes. (C) KEGG annotation analysis for unigenes

Identification of DEGs and related GO/KEGG enrichment analysis

Differential expression analysis revealed that 1805 unigenes were markedly up-regulated and 1460 unigenes were notably down-regulated in diseased lung tissues of forest musk deer compared to the normal lung tissue group (Table S3).The heatmap of these DEGs was shown in Fig. 4A. GO enrichment analysis showed that these DEGs were implicated in a host of biological processes including innate/adaptive immune responses, defense responses to pathogens, metabolism, and cell/tissue homeostasis (Table S4). Additionally, multiple KEGG pathways (e.g., cytokine-cytokine receptor interaction, viral protein interaction with cytokine and cytokine receptor, cell adhesion molecules, JAK-STAT/TNF/NF-kappa B signaling pathways, and phagosome) were significantly enriched by these DEGs (Table S5), while these pathways have been well documented to be closely associated with innate and adaptive immune responses. The top 20 KEGG pathways are presented in Fig. 4B [20].

Fig. 4
figure 4

Identification of DEGs and related GO/KEGG enrichment analysis. (A) Heatmap of DEGsin diseased lung tissues of forest musk deer compared to the normal lung tissue group.The R package v3.6.3 and ComplexHeatmapR package v2.2.0 were used to create the heatmap. (B) Top 20 KEGG pathways enriched by DEGs [20]

Screening of DEGs related to lung injury, immunity, and pathogen infection

Next, DEGs related to lung development (Table S6), immunity (Table S7), or bacterial (Table S8)/viral (Table S9)/fungal (Table S10) infection were screened out from Table S4 by searching for items containing the keywords “lung”, “immune/immunoglobulin/immunity”, “bacterial/bacterium”, “virus/viral” or “fungiform/fungus/fungal”. The aberrant expression of these unigenes might be closely linked with the pathogenesis of ARDS in forest musk deer. Next, Venn analysis of the DEG sets in Tables S6-S10 was performed to identify the DEGs that might play vital roles in the above-mentioned two or more biological processes. The outcome of the Venn analysis is presented in Fig. 5 and Table S11. The results showed that 2, 2, 2 genes related to lung development might also be implicated in immunity (Table S6|Table S7), viral infection (Table S6|Table S9), and fungal infection (Table S6|Table S10). Additionally, 35 unigenes, 30 unigenes, and 1 unigenes were found to be involved in the immune responses to bacteria (Table S7|Table S8), viruses (Table S7|Table S9), and fungi (Table S7|Table S10), respectively. Moreover, 4 unigenes (TCONS_00011905, TCONS_00013874, TCONS_00032940, and TCONS_00030832) might be associated with infections in both bacteria and viruses (Table S8|Table S9). Additionally, 3 unigenes (TCONS_00016874, TCONS_00022939, and TCONS_00022940) related to both bacterial and fungal infections were identified (Table S8|Table S10). Additionally, there were 13 common genes in Tables S7-S9 (Table S7|Table S8|Table S9), which might play vital roles in the immune defenses against both bacteria and viruses. Additionally, 7 common unigenes were identified in Table S7, S8, and S10 (Table S7|Table S8|Table S10). These unigenes (Table S7|Table S8|Table S10) might be vital players in the immune responses to both bacteria and fungi. Furthermore, there was 1 common unigenes (TCONS_00028682) in Tables S7-S10 (Table S7|Table S8|Table S9|Table S10), suggesting that TCONS_00028682 (lactotransferrin) might be involved in the immune responses to bacteria, viruses, and fungi.

Fig. 5
figure 5

Venn analytical outcomes of unigene sets in Tables S6-S10. The upper section shows the Venn diagram of unigenes in Tables S6-S10. The middle histogram presents the number of unigenes in Tables S6-S10.The chart below displays the number of specific (1) elements of one list or shared elements by 2, 3, … lists. For instance, a total of 332 unigenes were specific to one list. In addition, 79 unigenes, 20unigenes, and 1 unigenewere shared by 2, 3, and 4 of the 5 lists, respectively


In this study, our RNA-seq analysis showed that 3265 unigenes were differentially expressed (1805 up-regulated, 1460 down-regulated) in lung tissues of forest musk deer with ARDS relative to the non-ARDS group. Moreover, we identified 29, 303,105, 99, and 18 DEGs that were potentially associated with lung development (Table S6), immunity (Table S7), bacterial infection (Table S8), viral infection (Table S9), and fungal infection (Table S10), respectively. Moreover, we further filtered out 2(TCONS_00023532 and TCONS_00023537), 2(TCONS_00040651, and TCONS_00040652), and 2(TCONS_00013967 and TCONS_00013968) DEGs that might play vital roles in lung injury related to immune system dysfunction, viral infection, or fungal infection, respectively. Among these genes, sonic hedgehog (SHH, corresponding unigenes IDs: TCONS_00013967 and TCONS_00013968) has been implicated in immunity, acute lung injury, and fungal infection. For instance, Karnam et al. demonstrated that fungi could activate the SHH signaling pathway, while SHH downstream genes (5/15-lipoxygenases (LO) or cyclooxygenase (COX)-2) could inhibit or promote proinflammatory cytokine production [24]. Chen et al. demonstrated that the inhibition of SHH signaling by its inhibitor cyclopamine could aggravate lung tissue damage and increase the thickness of alveolar septa in lipopolysaccharide (LPS, Gram-negative bacterium cell wall component)-induced acute lung injury mice [25]. Our RNA-seq data showed that SHH expression was notably down-regulated in diseased lung tissues of forest musk deer with ARDS compared with the normal group. These data suggested the close association of aberrant SHH expression and ARDS-related lung injury in forest musk deer.

Moreover, we identified 13 DEGs (e.g.,NLRC5, NLRP3 (NACHT, LRR and PYD domains-containing protein 3), tyrosine-protein kinase FGR, B-cell leukemia-3 (Bcl-3), neutrophil gelatinase-associated lipocalin (NGAL), and secretory leukocyte protease inhibitor (SLPI) that might play key roles in the immune defenses against both bacteria and viruses (Table S7|Table S8|Table S9). NLRC5 and NLRP3, members of the nucleotide-binding domain leucine-rich repeat-containing (NLR) family, have been identified as crucial players in inflammation, immunity, and host defenses against viruses and bacteria [26,27,28]. Moreover, prior studies showed that NLRP3 deficiency could markedly alleviate lung injury, ameliorate lung vascular permeability, reduce neutrophil infiltration into the lung, and inhibit pro-inflammatory cytokine secretion in an LPS-induced ALI mouse model [29]. Previous reports also pointed out that the overexpression of tyrosine-protein kinase FGR could weaken host defenses against viruses, bacteria, and parasites [30,31,32]. Bcl-3, a member of the IκB (inhibitor of NF-κB) family, functions as a vital player in immunity and NF-κB activity regulation [33, 34]. Kreisel et al. showed that Bcl-3 could attenuate acute inflammatory lung injury in mice [35]. NGAL plays a central role inantimicrobial immune responses, cell growth, metastasis, apoptosis, and differentiation [36,37,38]. Moreover, lipocalin 2 (Lcn2) (the mouse homolog of human NGAL) expression has been found to be markedly increased in septic mice with acute lung injury compared to those without acute lung injury [39]. SLPI, one of the first lines of defense against multiple bacteria (e.g.,E. coli, Staphylococcus aureus, Mycobacterium bovis, Staphylococcus epidermidis, Mycobacterium tuberculosis, and Pseudomonas aeruginosa) and fungi (e.g., Aspergillus fumigatus and Candida albicans), can protect the host from excessive inflammatory responses/tissue injury and initiate tissue/cell repair/regenerative processes [40]. Additionally, some studies revealed that SLPI plays a crucial role in protecting lung tissues from inflammatory and proteolytic damage and microbial infection [41, 42]. These data further suggested the vital roles of these genes in the immune defenses against pathogens.

Moreover, GO analysis suggested that tyrosine-protein kinase SYK, peptidoglycan recognition protein 1 (PGRP1), cathelicidin-6, cathelicidin-4, and chromogranin-A might serve as crucial players in the immune responses to both bacteria and fungi. Spleen tyrosine kinase (Syk) is involved in multiple biological processes such as innate and adaptive immunity, cellular adhesion, and vascular development [43, 44]. Moreover, previous studies showed that the inhibition of Syk reduced airway inflammation and facilitated mouse survival in an LPS-ALI mouse model [45].The repeated challenge of the fungal pathogen Aspergillus fumigates or β-glucan (the main cell wall component of various fungi) could lead to serious lung injury and the inhibition of Syk ameliorated β-glucan-induced inflammation, Th17/Treg imbalance, and lung tissue damage in BALB/C mice [46]. Peptidoglycan recognition protein 1 (PGRP1) belongs to the family of peptidoglycan recognition proteins, which not only have microbicidal activities against Gram-positive and Gram-negative bacteria and fungi but can also participate in immune and inflammatory responses [47,48,49]. Cathelicidins are a major family of host defense peptides that function as crucial players in inflammatory and immune responses and defenses against infections by pathogens including viruses, Gram-positive and Gram-negative bacteria, and fungi [50, 51]. Moreover, the expression levels of cathelicidin-3L2,cathelicidin-5,cathelicidin-6, cathelicidin-7, and cathelicidin-9 were notably increased in lung tissues of forest musk deer who died of purulent diseases compared with healthy individuals, while cathelicidin-4 was expressed at low levels in lung tissues of purulent individuals [52]. Our RNA-seq data showed that cathelicidin-3, cathelicidin-4, cathelicidin-5, cathelicidin-6, and cathelicidin-7 were all highly expressed in lung tissues of forest musk deer that died of ARDS, suggesting the close association of these cathelicidin family genes with lung injury and immune responses to bacteria and fungi. Lactotransferrin (lactoferrin), a glycoprotein widely present in most body fluids of mammals, possesses multiple pharmacological properties such as antiviral, antiparasitic, antifungal, anti-inflammatory, antimicrobial,and immunity-regulatory effects [53, 54]. Consistent with these data, we found that lactotransferrin was closely linked with immunity and pathogen (bacterium, virus, and fungus) infection.

Taken together, 3265 DEGs were identified in diseased lung tissues of forest musk deer that died of ARDS compared with normal lung tissues. This is the first study to disclose the lung gene alterations caused by ARDS in forest musk deer. Moreover, a total of 6 potential unigenes related to lung damage caused by immune abnormality (Table S6|Table S7) or viral (Table S6|Table S9)/fungal (Table S6|Table S10) infection were identified, which might contribute to the better management of lung injury induced by immune system dysfunction or pathogen (virus/fungus) infection. Additionally, we identified some unigenes that might play vital roles in the immune responses to bacteria (Table S7|Table S8), viruses (Table S7|Table S9), or fungi (Table S7|Table S10). Additionally, multiple unigenes that were involved in immune defenses against coinfection with multiple types of pathogens (bacteria + viruses: Table S7|Table S8|Table S9; bacteria + fungi: Table S7|Table S8|Table S10; bacteria + viruses + fungi: Table S7|Table S8|Table S9|Table S10) were screened out. Furthermore, DEGs in response to the single infection or dual infection of pathogens (bacteria, viruses, or fungi) were filtered out. We supposed that these unigenes might have potential value in the elimination of pathogens and improvement of resistance to bacterial/viral/fungal infection in forest musk deer. Moreover, our study might contribute to the better management of lethal acute respiratory diseases in forest musk deer. Both our data and previous studies showed that pathogen infection and immune dysregulation can trigger lung injury and result in forest musk deer death. Thus, it is imperative to maintain good breeding conditions and improve the immunity of forest musk deer. Treatment of an appropriate amount of immune enhancer might improve the immune and defensive abilities of forest musk deer against pathogens and benefit the long-term development of the forest musk deer breeding system.

Availability of data and materials

The raw data of RNA-seq has been deposited to the GEO ( =) with the number of GSE191259. The data that support the findings of this study are available from the corresponding author upon reasonable request.


  1. Wang Y, Harris R. Moschus berezovskii.(errata version published in 2016) The IUCN Red List of Threatened Species 2015: e. T13894A103431781(downloaded on 09 January 2017) 2015.

  2. 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(3):333–42.

    Article  Google Scholar 

  3. 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(3):166–73.

    Article  CAS  Google Scholar 

  4. Wang J, Xing H, Qin X, Ren Q, Yang J, Li L. Pharmacological effects and mechanisms of muscone. J Ethnopharmacol. 2020;262:113120.

    Article  CAS  PubMed  Google Scholar 

  5. 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. Chinese medicine. 2021;16(1):46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Meng X, Gong B, Ma G, Xiang L. Quantified analyses of musk deer farming in China: a tool for sustainable musk production and ex situ conservation. Asian Australas J Anim Sci. 2011;24(10):1473–82.

    Article  Google Scholar 

  7. Shrestha M. Animal welfare in the musk deer. Appl Anim Behav Sci. 1998;59(1–3):245–50.

    Article  Google Scholar 

  8. Fan J, Zheng X, Wang H, Qi H, Jiang B, Qiao M, Zhou J, Bu S: Analysis of genetic diversity and population structure in three forest musk deer captive populations with different origins. G3 (Bethesda). 2019;9(4):1037–1044.

  9. 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(2):104–7.

    Google Scholar 

  10. Li Y, Hu X, Yang S, Zhou J, Zhang T, Qi L, Sun X, Fan M, Xu S, Cha M, et al. Comparative analysis of the gut microbiota composition between captive and wild forest musk deer. Front Microbiol. 2017;8:1705.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Kang F, Liu Z, Zhang X, Wang C, Deng K. Preliminary studies on causes of death of Alpine musk deer(Moschus chrysogaster) at the Xinglongshan Musk Deer Farm. Gansu Acta Theriologica Sinica. 2008;28(4):430–3.

    Google Scholar 

  12. Banavasi H, Nguyen P, Osman H, Soubani AO. Management of ARDS - What Works and What Does Not. Am J Med Sci. 2021;362(1):13–23.

    Article  PubMed  Google Scholar 

  13. Huppert LA, Matthay MA, Ware LB. Pathogenesis of acute respiratory distress syndrome. Semin Respir Crit Care Med. 2019;40(1):31–9.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Umbrello M, Formenti P, Bolgiaghi L, Chiumello D. Current concepts of ARDS: a narrative review. Int J Mol Sci. 2016;18(1):64.

    Article  PubMed Central  Google Scholar 

  15. Kumar V. Pulmonary innate immune response determines the outcome of inflammation during pneumonia and sepsis-associated acute lung injury. Front Immunol. 2020;11:1722.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Hooven TA, Polin RA. Pneumonia. Semin Fetal Neonatal Med. 2017;22(4):206–13.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Zhao Z, Meng F, Wang W, Wang Z, Zhang C, Jiang T. Comprehensive RNA-seq transcriptomic profiling in the malignant progression of gliomas. Sci Data. 2017;4:170024.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Jenkins BH, Maguire F, Leonard G, Eaton JD, West S, Housden BE, Milner DS, Richards TA. Characterization of the RNA-interference pathway as a tool for reverse genetic analysis in the nascent phototrophic endosymbiosis, Paramecium bursaria. R Soc Open Sci. 2021;8(4):210140.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Macmanes MD. On the optimal trimming of high-throughput mRNA sequence data. Front Genet. 2014;5:13.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49:545–51.

    Article  Google Scholar 

  21. Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talón M, Dopazo J, Conesa A. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008;36(10):3420–35.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.

    Article  CAS  PubMed  Google Scholar 

  23. Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M. KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007;35(Web Server issue):W182-185.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Karnam A, Holla S, Balaji KN. Sonic hedgehog-responsive lipoxygenases and cyclooxygenase-2 modulate Dectin-1-induced inflammatory cytokines. Mol Immunol. 2015;68(2 Pt A):280–9.

    Article  CAS  PubMed  Google Scholar 

  25. Chen X, Jin Y, Hou X, Liu F, Wang Y. Sonic hedgehog signaling: evidence for Its protective role in endotoxin induced acute lung injury in mouse model. PLoS ONE. 2015;10(11):e0140886.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Lupfer C, Kanneganti TD. The expanding role of NLRs in antiviral immunity. Immunol Rev. 2013;255(1):13–24.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Brodsky IE, Monack D. NLR-mediated control of inflammasome assembly in the host response against bacterial pathogens. Semin Immunol. 2009;21(4):199–207.

    Article  CAS  PubMed  Google Scholar 

  28. Yao Y, Qian Y. Expression regulation and function of NLRC5. Protein Cell. 2013;4(3):168–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Peng Y, Wu Q, Tang H, Chen J, Wu Q, Yuan X, Xiong S. NLRP3 Regulated CXCL12 Expression in Acute Neutrophilic Lung Injury. J Inflamm Res. 2020;13:377–86.

  30. Nelson MP, Metz AE, Li S, Lowell CA, Steele C. The absence of Hck, Fgr, and Lyn tyrosine kinases augments lung innate immune responses to Pneumocystis murina. Infect Immun. 2009;77(5):1790–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Liu S, Chen S, Li X, Wu S, Zhang Q, Jin Q, Hu L, Zhou R, Yu Z, Meng F, et al. Lck/Hck/Fgr-Mediated Tyrosine Phosphorylation Negatively Regulates TBK1 to Restrain Innate Antiviral Responses. Cell Host Microbe. 2017;21(6):754-768.e755.

    Article  CAS  PubMed  Google Scholar 

  32. Lowell CA, Berton G. Resistance to endotoxic shock and reduced neutrophil migration in mice deficient for the Src-family kinases Hck and Fgr. Proc Natl Acad Sci USA. 1998;95(13):7580–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Herrington FD, Nibbs RJ. Regulation of the adaptive immune response by the IκB family protein Bcl-3. Cells. 2016;5(2):14.

    Article  PubMed Central  Google Scholar 

  34. Palmer S, Chen YH. Bcl-3, a multifaceted modulator of NF-kappaB-mediated gene transcription. Immunol Res. 2008;42(1–3):210–8.

    Article  CAS  PubMed  Google Scholar 

  35. Kreisel D, Sugimoto S, Tietjens J, Zhu J, Yamamoto S, Krupnick AS, Carmody RJ, Gelman AE. Bcl3 prevents acute inflammatory lung injury in mice by restraining emergency granulopoiesis. J Clin Investig. 2011;121(1):265–76.

    Article  CAS  PubMed  Google Scholar 

  36. Schmidt-Ott KM, Mori K, Li JY, Kalandadze A, Cohen DJ, Devarajan P, Barasch J. Dual action of neutrophil gelatinase-associated lipocalin. J Am Soc Nephrol. 2007;18(2):407–13.

    Article  CAS  PubMed  Google Scholar 

  37. Nasioudis D, Witkin SS. Neutrophil gelatinase-associated lipocalin and innate immune responses to bacterial infections. Med Microbiol Immunol. 2015;204(4):471–9.

    Article  CAS  PubMed  Google Scholar 

  38. Chakraborty S, Kaur S, Guha S, Batra SK. The multifaceted roles of neutrophil gelatinase associated lipocalin (NGAL) in inflammation and cancer. Biochem Biophys Acta. 2012;1826(1):129–69.

    CAS  PubMed  Google Scholar 

  39. Zeng G, Jia CW, Liu J, Guo SB. Lipocalin-2 test in distinguishing acute lung injury cases from septic mice without acute lung injury. Chin Med Sci J. 2014;29(2):65–77 (Chung-kuo i hsueh k&apos;o hsueh tsa chih).

    Article  CAS  PubMed  Google Scholar 

  40. Majchrzak-Gorecka M, Majewski P, Grygier B, Murzyn K, Cichy J. Secretory leukocyte protease inhibitor (SLPI), a multifunctional protein in the host defense response. Cytokine Growth Factor Rev. 2016;28:79–93.

    Article  CAS  PubMed  Google Scholar 

  41. McKiernan PJ, McElvaney NG, Greene CM. SLPI and inflammatory lung disease in females. Biochem Soc Trans. 2011;39(5):1421–6.

    Article  CAS  PubMed  Google Scholar 

  42. Zani ML, Tanga A, Saidi A, Serrano H, Dallet-Choisy S, Baranger K, Moreau T. SLPI and trappin-2 as therapeutic agents to target airway serine proteases in inflammatory lung diseases: current and future directions. Biochem Soc Trans. 2011;39(5):1441–6.

    Article  CAS  PubMed  Google Scholar 

  43. Alhazmi A. Spleen tyrosine kinase as a target therapy for pseudomonas aeruginosa infection. J Innate Immun. 2018;10(4):255–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Mócsai A, Ruland J, Tybulewicz VL. The SYK tyrosine kinase: a crucial player in diverse biological functions. Nat Rev Immunol. 2010;10(6):387–402.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Nadeem A, Ahmad SF, Al-Harbi NO, Al-Harbi MM, Ibrahim KE, Kundu S, Attia SM, Alanazi WA, AlSharari SD. Inhibition of spleen tyrosine kinase signaling protects against acute lung injury through blockade of NADPH oxidase and IL-17A in neutrophils and γδ T cells respectively in mice. Int Immunopharmacol. 2019;68:39–47.

    Article  CAS  PubMed  Google Scholar 

  46. Patel D, Gaikwad S, Challagundla N, Nivsarkar M, Agrawal-Rajput R. Spleen tyrosine kinase inhibition ameliorates airway inflammation through modulation of NLRP3 inflammosome and Th17/Treg axis. Int Immunopharmacol. 2018;54:375–84.

    Article  CAS  PubMed  Google Scholar 

  47. Dziarski R, Gupta D. How innate immunity proteins kill bacteria and why they are not prone to resistance. Curr Genet. 2018; 64(1):125–9.

  48. Dziarski R, Gupta D. The peptidoglycan recognition proteins (PGRPs). Genome Biol. 2006;7(8):232.

    Article  PubMed  PubMed Central  Google Scholar 

  49. De Marzi MC, Todone M, Ganem MB, Wang Q, Mariuzza RA, Fernández MM, Malchiodi EL. Peptidoglycan recognition protein-peptidoglycan complexes increase monocyte/macrophage activation and enhance the inflammatory response. Immunology. 2015;145(3):429–42.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Alford MA, Baquir B, Santana FL, Haney EF, Hancock REW. Cathelicidin Host defense peptides and inflammatory signaling: striking a balance. Front Microbiol. 2020;11:1902.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Agier J, Efenberger M, Brzezińska-Błaszczyk E. Cathelicidin impact on inflammatory cells. Cent-Eur J Immunol. 2015;40(2):225–35.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Zhang L, Jie H, Xiao Y, Zhou C, Lyu W, Bai W. Genomic identification and expression analysis of the cathelicidin gene family of the forest musk deer. Animals. 2019;9(8):481. an open access journal from MDPI

    Article  PubMed Central  Google Scholar 

  53. Albar AH, Almehdar HA, Uversky VN, Redwan EM. Structural heterogeneity and multifunctionality of lactoferrin. Curr Protein Pept Sci. 2014;15(8):778–97.

    Article  CAS  PubMed  Google Scholar 

  54. Hao L, Shan Q, Wei J, Ma F, Sun P. Lactoferrin: major physiological functions and applications. Curr Protein Pept Sci. 2019;20(2):139–44.

    Article  CAS  PubMed  Google Scholar 

Download references


I would like to express my gratitude to all those who have helped me during the writing of this thesis. I gratefully acknowledge the help of the Science and Technology Projects of Shaanxi Academy of Science that funded this research. And, we thank the Wuhan GeneCreate Biological Engineering Co., Ltd. for providing transcriptome sequencing.


This research was supported by Science and Technology Projects of Shaanxi Academy of Science, China (Grant No. 2018 K-01, Grant No. 2021 K-37) and Foundation of Shaanxi Academy of Forestry, China(SXLK2021-0222).

Author information

Authors and Affiliations



Jie Tang designed and performed the experiments and wrote the manuscript. Chao Yang, Lijuan Suo, and Feiran Li were involved in sample collection and the implementation of this experiment. Kun Bian and Yan Wang revised the manuscript. Yan Wang participated in the design of this research. All authors have read and approved the final manuscript.

Corresponding author

Correspondence to Chao Yang.

Ethics declarations

Ethics approval and consent to participate

No animals were anesthetized or euthanized as part of this study. All sample collections and study protocol were approved by the Ethics Committee of Shaanxi Institute of Zoology and conducted following the ARRIVE 2.0 guidelines.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no financial conflicts of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Tang, J., Suo, L., Li, F. et al. Transcriptome profiling of lung immune responses potentially related to acute respiratory distress syndrome in forest musk deer. BMC Genomics 23, 701 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Forest musk deer
  • RNA sequencing
  • Lung
  • Immunity
  • Bacterium
  • Virus
  • Fungus