Whole blood transcriptome profiling identifies candidate genes associated with alopecia in male giant pandas (Ailuropoda melanoleuca)
BMC Genomics volume 23, Article number: 297 (2022)
The giant panda (Ailuropoda melanoleuca) is a threatened species endemic to China. Alopecia, characterized by thinning and broken hair, mostly occurs in breeding males. Alopecia significantly affects the health and public image of the giant panda and the cause of alopecia is unclear.
Here, we researched gene expression profiles of four alopecia giant pandas and seven healthy giant pandas. All pandas were approximately ten years old and their blood samples collected during the breeding season. A total of 458 up-regulated DEGs and 211 down-regulated DEGs were identified. KEGG pathway enrichment identified that upregulated genes were enriched in the Notch signaling pathway and downregulated genes were enriched in ribosome, oxidative phosphorylation, and thermogenesis pathways. We obtained 28 hair growth-related DEGs, and identified three hub genes NOTCH1, SMAD3, and TGFB1 in PPI analysis. Five hair growth-related signaling pathways were identified with abnormal expression, these were Notch, Wnt, TGF-β, Mapk, and PI3K-Akt. The overexpression of NOTCH1 delays inner root sheath differentiation and results in hair shaft abnormalities. The delayed hair regression was associated with a significant decrease in the expression levels of TGFB1.
Our data confirmed the abnormal expression of several hair-related genes and pathways and identified alopecia candidate genes in the giant panda. Results of this study provide theoretical basis for the establishment of prevention and treatment strategies for giant pandas with alopecia.
The main structures of hair include the dermal papilla, bulge, sebaceous gland, and hair shaft . The hair shaft is produced by a hair follicle (HF) and possesses a unique ability to regenerate itself. Morphogenesis of the HF occurs in the embryo, and postnatal follicles undergo a cycle of renewal in three phases: anagen (growth), catagen (degradation) and telogen (resting) . Follicle stem cells are activated at the telogen-to-anagen transition, to initiate a new round of hair growth . Many genes and signaling pathways are involved in HF development, such as Notch-1 receptor and its three ligands, Delta-1, Jagged-1 and Jagged-2 in Notch , Wnt10b in Wnt , TGF-β2 in TGF-β , Gab1 in Mapk , and FGF16, CSF3, interleukin (IL)6, and oncostatin M in PI3K-Akt  signaling pathways. The investigation of genes and molecular signaling in the regulation of hair follicles is ongoing.
Mammalian hair has important biological functions, including sunlight reflection or absorption, heat retention, water resistance, mate attraction, skin protection and social communication . Despite these important biological functions, premature or atypical hair loss (i.e. alopecia) is common. For example, recent studies have reported that rhesus macaques , mice , Andean bears  and polar bears  displayed some extent of alopecia. Alopecia in captive female rhesus macaques occurred during gestation and parturition  and was associated with pregnancy and chronic stress . Whereas otherwise healthy captive Andean bears acquired the slowly progressive alopecia syndrome of unknown aetiology . While 28% of polar bears in the southern Beaufort Sea region had varying degrees of alopecia that was concomitant with body condition declines . Numerous factors may contribute to alopecia, ranging from naturally occurring processes (seasonality, aging) to various diseases (stress, endocrine disorders, genetic mutations, immunologic diseases, bacterial and fungal infections) [13, 14]. Alopecia is not fully understood but it is a common and complex phenomenon in captive mammals .
The endemic giant panda (Ailuropoda melanoleuca), is a national treasure in China and a flagship species within China and internationally. The unique bamboo-specialist carnivore has become a global conservation symbol due to efforts to save the species, by overcoming difficulties in captive breeding, have been captured on film and now on social media. Humans prefer animals with physiological and behavioral similarities and these bears often display behaviors similar to human children, which are frequently comedic and therefore are loved by people across the world . Akin to other captive species, the giant panda can develop alopecia and mostly occurs in breeding males in spring and summer. The alopecia, characterized by symmetrical thinning and broken hair or hair loss, is mainly concentrated on the limbs and abdomen. Occasionally pandas can become bald in certain parts of the body. Alopecia can seriously affect the health of giant pandas, but can also damage the public image of the flagship species.
At present, the causes of giant panda alopecia are unclear and might be complex. It is likely that mapping and analyzing genes involved in alopecia will assist with understanding the cause and may provide avenues for treatment. Blood is the main component of the animal immune system and plays a key role in the animal immune system. At the same time, the gene expression of blood cells that circulate through various tissues and organs can experience unique changes due to disease or injury in other tissues and cells of the body. Additionally, many genes previously thought to be expressed only in non-blood tissues are also expressed in peripheral blood cells . Therefore, the phenotype of blood cell can reflect the physiological and pathological events of other tissues, and different blood markers have been widely used as substitutes for organs to monitor their health status . Next-generation sequencing analysis has greatly improved genetic diagnostics and therefore, we aimed to mine the differentially expressed genes (DEGs) in male giant panda peripheral blood with and without alopecia and discuss the possible cause of giant panda alopecia.
Material and methods
We collected peripheral blood samples from 11 captive adult male giant pandas (9 to 11 years old) from April to June. These 11 pandas were housed in the China Conservation and Research Center for Giant Panda in Chengdu, China. Samples were taken during routine physical examinations to minimize stress and the need for additional captures (i.e. anesthesia). We divided the 11 pandas into two groups, with and without alopecia. The (with) alopecia group had four giant pandas and the normal (control) group consisted of seven pandas (Additional file 1: Table S1).
All giant pandas continued their normal captive existence after completion of the sampling and were in no way physically harmed due to our research. The China Conservation and Research Center for Giant Panda agreed to provide blood samples and they assisted with the research.
Library preparation and sequencing
Approximately 3 ml blood was collected from each individual and immediately transferred into a paxgene blood RNA tube (BD company, USA) to stabilize the RNA in the cell and avoid degradation. Total RNA was extracted using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s protocol and treated with RNase-free DNase I. Nanodrop 8000 Spectrophotometer (Thermo scientific) was used to evaluate the purity and concentration of RNA. We confirmed RNA integrity using an RNA PicoChip with Agilent 2100 Bioanalyzer (Agilent Technologies). We then used the extracted RNA samples for the cDNA synthesis and amplified double-stranded cDNA. We used kits from the Illumina Company to perform sequencing libraries construction, quality control and quantification as per manufacture's recommendations. We sequenced the cDNA library on the Illumina sequencing platform (HiSeq 2000) according to standard procedures. The Genome Analyzer Pipeline version 2.0 ran the image analysis and base-calling using default parameters. Our sequencing generated 150-bp paired-end reads.
Quantification and mapping
All raw reads were processed with adapter trimming and reads filtering by trim-galore version 0.5.0 (http://www.bioinformatics.babraham.ac.uk/projects/trim_galore/) to ensure data quality control. We used FastQC version 0.11.5 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) to generate quality reports thus ensuring that clean data were used in subsequent analyses. We built a genome index and mapped reads data on the giant panda reference genome using HISAT2 version 2.1.0 . The giant panda reference genome (v90) and reference annotation were downloaded from the Ensembl website (ftp://ftp.ensembl.org/pub/release-90/).
Calculating differentially expressed genes (DEGs)
The alignment program HISAT2 generated SAM files that stored mapping information and these files were sorted by SAMtools version 1.7 . We obtained a sample read counts file from BAM files using featureCounts version v1.6.2 , allowing quicker counting of features. Reads counts were used to calculate the expression value of transcripts per million (TPM). The data matrix of the number of reads was analyzed using principal components analysis (PCA) in R function prcomp. A PCA plot with clustering information was illustrated by R package ggbiplot (https://github.com/vqv/ggbiplot). The counts file was entered into R package edgeR  that identified genes with significant expression differences based on the Poisson model. We calculated the expression fold change (FC) and false discovery rate (FDR) to identify differentially expressed genes (DEGs) in the two groups. We set a cut-off of 1 for |log2FC| and 0.05 for FDR. The filtered genes were retained for further analysis.
Analysis of gene enrichment
Genes were annotated and analyzed using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases to obtain known functions of these gene sets. We compared our dataset with three GO categories; molecular function, biological process and cellular component. Genes were clustered to GO terms using a web server g:Profiler . We calculated the FDR using g:SCS algorithm and a threshold of 0.05 was set for FDR. The KEGG pathway map analysis was performed by KOBAS 2.0  and this analysis presented biological interpretation of higher-level systemic functions,.
Analysis of protein–protein interaction network
We input DEGs into the STRING  database of known and predicted protein–protein interactions to obtain protein–protein interaction networks. Protein–protein interaction outputs were simple texts in tabular form. We used Cytoscape  to visualize molecular interactions and a Cytoscape plugin, cytoHubba , to calculate hub genes. The EcCentricity algorithm was used to calculate hub genes. Betweeness, closeness, and degree topological algorithms were used for verification.
Real-time quantitative PCR (qRT-PCR)
Real-time quantitative PCR (qRT-PCR) was performed to confirm the expression changes of genes and six genes were selected randomly. We used the same reference gene, GAPDH , as previous studies. The real-time quantitative PCR was performed with six samples (A1, A2, A3, N1, N2, N3) because RNA of several samples had already been used. The sequences of primers were predicted using Primer-BLAST  and were synthesized at TSINGKE (listed in Additional file 2: Table S2). Each qPCR reaction system was 10 μL containing 5 μL of 2 × M5 Hiper SYBR Premix EsTaq (mei5, Beijing, China), 0.2 μL cDNA, 0.2 μL of forward and reverse primer and nuclease-free water. Reactions were performed in triplicate on Bio-Rad CFX96 Touch. We used the following program: 95 °C for 30 s, followed by 40 cycles of 95 °C/10 s, 60 °C/15 s; then 72 °C for 10 s. Gene expression levels were calculated using optimized comparative Ct (2−ΔΔCt) value method . GraphPad Prism 8.0.2 (https://www.graphpad.com/scientific-software/prism/) was used to analyze data.
Reads sequencing and processing
Clean reads data were generated from raw Illumina RNA-Seq data. A total of 79.27 Gb of paired-end clean data were generated (Additional file 1: Table S1) and the percent of Q30 was greater than 85% according to FastQC. HISAT2 mapping results indicated that sample overall alignment rates were between 86 and 90%. We conducted read summarization using program featureCounts and counts were transformed into a numerical matrix. We found that the 11 samples were divided into two groups according to a PCA normalized matrix, which agreed with the groups being divided by phenotype (Fig. 1).
Identification of DEGs
We identified 458 up-regulated DEGs and 211 down-regulated DEGs in the alopecia group compared to the normal group (Additional file 3: Table S3). Of the 458 up-regulated genes, 365 were annotated, while 140 of the 211 down-regulated DEGs had annotations. Our study was mainly based on annotated genes. We identified the most significant DEGs by sorting by log2FC value. We selected the top 50 up-regulated and down-regulated genes to draw the gene expression heatmap (Fig. 2).
In the top 10 up-regulated DEGs, SHISA7 (shisa family member 7) and ASCL5 (achaete-scute family bHLH transcription factor 5) were 10 and sixfold higher in alopecia males than normal males, respectively. Also identified in the top 10, TMEM170B (transmembrane protein 170B), IFITM5 (interferon induced transmembrane protein 5), and FZD8 (frizzled class receptor 8) are associated with Wnt signaling. Wnt signaling participates in numerous fundamental processes, such as embryonic development and hair growth.
In the top 10 down-regulated DEGs, KLK5 (kallikrein related peptidase 5) and TCEA3 (transcription elongation factor A3) were 10 and sixfold lower in alopecia males than normal males, respectively. KLK5 is highly expressed in human hair follicles and sweat glands, and is involved in keratinocyte desquamation and pigmentation. Other top 10 down-regulated DEGs included CD52 (CD52 molecule) and TRDC (T cell receptor delta constant) and these are related to the immune system. TRDC is involved in the formation of T cell receptors. RPL21 (ribosomal protein L21) was also a top 10 and is known to encode a ribosomal protein.
Gene Ontology enrichment of DEGs
We found that up-regulated DEGs were enriched in 37 GO terms, consisting of three terms in molecular function and 34 terms in biological process (Fig. 3). The three molecular function terms were transcription regulator activity (GO:0,140,110), semaphorin receptor binding (GO:0,030,215), and binding (GO:0,005,488). The 34 biological process terms included negative regulation of myotube differentiation (GO:0,010,832), regulation of cell communication (GO:0,010,646), and regulation of response to stimulus (GO:0,048,583).
The down-regulated DEGs were enriched in 78 GO terms, with 14 terms in molecular function, 28 terms in biological process, and 36 terms in cellular component (Fig. 4). The top three biological process terms were organonitrogen compound biosynthetic process (GO:1,901,566), peptide metabolic process (GO:0,006,518), and translation (GO:0,006,412). The top three cellular component terms were ribosome (GO:0,005,840), respirasome (GO:0,070,469), and mitochondrial inner membrane (GO:0,005,743). The terms in molecular function mainly included structural constituent of ribosome (GO:0,003,735), NADH dehydrogenase activity (GO:0,003,954), mRNA 5'-UTR binding (GO:0,048,027), and cytochrome-c oxidase activity (GO:0,004,129). All GO term enrichments are shown in Additional file 4: Table S4.
Pathway enrichment of DEGs
Up-regulated DEGs were significantly enriched in the Notch signaling pathway (aml04330; Fig. 5). The Notch signaling pathway controls the fate of adult hair follicular stem cells in the bulge. We found that down-regulated DEGs were enriched in four pathways (Fig. 5). These were the ribosome pathway (aml03010), thermogenesis pathway (aml04714) oxidative phosphorylation (aml00190) and ferroptosis (aml04216). The thermogenesis pathway is the child term of environmental adaptation. Oxidative phosphorylation (aml00190) is the downstream term of thermogenesis. Ferroptosis (aml04216) is a regulated form of cell death.
Expression of hair growth-associated genes
We identified 28 hair growth-associated genes with altered expression (Table 1). DUSP3, GNG5, RAC1, and TGFB1 were found to be highly down-regulated while the other genes were all up-regulated. These genes belonged to five signaling pathways, being Notch, Wnt, TGF-β, Mapk, and PI3K-Akt signaling pathways. NOTCH1 (notch receptor 1) and DTX4 (deltex E3 ubiquitin ligase 4) participate in the Notch signaling pathway. DVL1 (dishevelled segment polarity protein 1), DVL2 (dishevelled segment polarity protein 2), FZD8 (frizzled class receptor 8), NFATC2 (nuclear factor of activated T cells 2), and VANGL2 (VANGL planar cell polarity protein 2) are linked to the Wnt signaling pathway. CDKN2B (cyclin dependent kinase inhibitor 2B), SMAD3 (SMAD family member 3), SMAD7 (SMAD family member 7), and TGFB1 (transforming growth factor beta 1) are members of the TGF-β signaling pathway. DUSP4 (dual specificity phosphatase 4), DUSP7 (dual specificity phosphatase 7), IKBKG (inhibitor of nuclear factor kappa B kinase regulatory subunit gamma), and RAC1 (Rac family small GTPase 1) function within the Mapk signaling pathway. Lastly, IL6R (interleukin 6 receptor), OSM (oncostatin M), and RXRA (retinoid X receptor alpha) are associated with the PI3K-Akt signaling pathway.
Protein–protein interaction network of hair growth-associated DEGs
We converted hair growth-associated DEGs into proteins using STRING. After removing 12 isolated nodes, a total of 22 interaction edges between 16 nodes were extracted from the database. Functional enrichment of STRING local network cluster was then constructed. The genes were significantly enriched in 13 clusters, including the Notch signaling pathway (CL:18,099), the TGF-β signaling pathway (CL:17,598), and the Wnt signaling pathway (CL:17,371). We calculated the hub genes using cytoHubba and obtained a score of 36 genes using EcCentricity algorithm (Additional file 5: Table S5). The top three genes were NOTCH1, SMAD3, and TGFB1. These three genes were the candidate genes associated with alopecia in male giant pandas. We plotted the hub genes interaction network (Fig. 6).
Real-time quantitative PCR (qRT-PCR) validation
We selected six DEGs (IL4R, SH2B3, CMC1, TGFB1, RAC1, CTSC) for verification. The results of qRT-PCR indicated similar expression tendencies with transcriptome sequencing (Fig. 7). The qRT-PCR validation improves the reliability of the present study.
Hair is a unique character in mammals, and has important biological functions, such as mate attraction, skin protection, and social communication [1, 31]. Alopecia is a common phenomenon in mammal, including humans, and can seriously affect health. Here we monitored the expression of hair growth-related genes based on peripheral blood transcriptomes of giant panda. We identified 28 hair growth-related genes with altered expression in alopecia male giant pandas compared to normal male pandas.
The GO term enrichment of “ribosome” and “structural constituent of ribosome”, and KEGG pathway enrichment of “ribosome” were observed in down-regulated genes. Several studies have identified a connection between the ribosome and hair growth [32,33,34,35,36]. Lv et al. found that the ribosome pathway was the most significantly enriched pathway by skin transcriptome analysis and may be related to the development and the density of secondary wool follicles in Chinese Merino sheep . By miRNA sequencing and differentially expressed miRNAs analysis, the ribosome pathway may participate in dermal papilla cells viability and proliferation, and thereby affects hair cycling . Plowman et al. showed that ribosomal proteins were detected in dissected portions from 30 hair follicles by Gel-free proteomic analysis . In patients with severe active alopecia areata, the ribosome pathway was significantly enriched by functional annotation analysis . The ribosome may influence protein synthesis that is important for human dermal papilla cells . The oxidative phosphorylation pathway was down-regulated in giant pandas with alopecia. Oxidative phosphorylation can produce ATP for hair follicle stem cells differentiation and proliferation . Our enrichment results are consistent with a previous study where the oxidative phosphorylation pathway was reported to be related to hair follicle cycling . The thermogenesis pathway was another enriched pathway in our study, which was unusual of the hair follicle cycle activities, underlining the relationship between environmental temperature and hair growth . The present study confirmed the abnormal expression of several pathways related to hair growth.
Notch signaling pathway
The Notch pathway was found to be related to Hu sheep hair follicle growth and development when the source genes of differentially expressed circRNAs were compared between small waves and straight wool . Additionally, the results of microRNA expression profiles showed that the Notch pathway may be an important regulator in the process of coat color formation in the goat (Capra hircus) . In the current study, the Notch signaling pathway was significantly enriched in alopecia giant pandas. Notch is a local cell-signaling pathway and comprises evolutionarily conserved transmembrane receptors and has three possible functional roles in hair follicle development: lateral inhibition, boundary formation and lineage decision . In the bulge, Notch signaling is required for follicular fate selection of adult hair follicle stem cells during hair follicle morphogenesis . In the hair bulb, Notch controls the proper development of the hair shaft and inner root sheath . Notch signaling cannot only suppress TGF-β to provide an optimal matrix proliferating environment, but it can also activate Wnt family members to regulate hair follicle keratinocyte differentiation . NOTCH1 is expressed in the follicle bulb and outer root sheath . Inactivation of Notch1 in adult mice is characterized by a premature entry into catagen and leads to alopecia . The overexpression of NOTCH1 delays inner root sheath differentiation and results in hair shaft abnormalities . NOTCH1 was highly expressed and up-regulated (about fourfold) in giant pandas with alopecia. DTX4 (deltex E3 ubiquitin ligase 4) encodes Deltex protein that is a positive regulator of Notch signaling . DTX4 was also up-regulated in alopecia giant pandas. Therefore, alopecia in our giant pandas may be related to dysregulation of Notch signaling and NOTCH1 may be the candidate gene.
Wnt signaling pathway
The Wnt signaling pathway includes canonical Wnt/β-catenin, Wnt/PCP, and Wnt/Ca2+ signaling pathways . In the canonical Wnt pathway, the Dvl protein is a key component of Wnt signaling and is recruited by Frizzled receptors and LRP co-receptors, preventing the constitutive destruction of cytosolic β-catenin [41, 46, 47]. The stabilization of β-catenin may play a role in hair follicle development . Wnt signaling was the earliest signal in the development of the epidermis and has been found in the embryonic skin of mice . The overexpression of DKK1, a Wnt inhibitor, could prevent the formation of hair follicles in mice . The miR-29 family was reported to inhibit the Wnt and BMP signaling pathways causing eventual hair loss and alopecia . The Wnt was reported to be involved in the hair follicle cycle in yaks (Bos grunniens) via mRNA and lncRNA analysis [50, 51]. FZD8 (frizzled class receptor 8) as Wnt family receptor, was found to be expressed in dermal papillae cells that controlled feather regeneration in chickens . FZD8 was up-regulated in the induction of dermal fibroblasts into dermal papilla cell-like cells in nude mice, suggesting the important role of FZD8 and the Wnt signaling pathway in hair follicle regeneration . DVL1 and DVL2 as the key regulators of Wnt signaling, were up-regulated 2.3-fold and 2.1-fold, respectively. Overexpression of DVL2 in the outer root may affect hair growth and structure, resulting in short hair . Lee et al. found CXXC-type zinc finger protein 5 interacted with Dvl protein via its C-terminal Dvl-binding motif to negatively regulate hair regrowth . In the Wnt/Ca2+ signaling pathway, NFATC2 (nuclear factor of activated T cells 2) was up-regulated in giant pandas with alopecia. The inhibition of the nuclear factor of activated T cells pathway could enhance hair growth in follicular keratinocytes .In the Wnt/PCP signaling pathway, VANGL2 (VANGL planar cell polarity protein 2) could drive the anterior-directed tilt of hair follicle placodes during early stages of hair follicle morphogenesis . In vitro evidence showed that VANGL2 plays a role in hair follicle polarization and orientation . Our findings implied that the Wnt signaling pathway was the candidate pathway in giant panda alopecia syndrome.
TGF-β signaling pathway
TGF-β signaling could activate intracellular effectors such as Smad proteins and kinases via ligands (TGF-β1/2/3, activins, BMPs, GDFs) bind to receptors (BMPR, TGFBR), controlling the activities of key transcription factors that promote epithelial differentiation . TGF‐β was reported to regulate various biological processes in the skin, including the inhibition of keratinocyte proliferation and the stimulation of epithelial cell apoptosis . During the hair follicle cycle of mice and humans, TGF‐β contributes to the induction of catagen . Western blot analysis showed that the delayed hair regression was associated with a significant decrease in the expression levels of TGFB1, and high expression of TGFB1 in the epidermis resulting in the inhibition of normal skin development . Liu et al. found that chronic expression of TGFB1 led to alopecia in adult transgenic mice . The down-regulation of TGFB1 in the giant panda may imply their partial responsibility for giant panda alopecia. All-trans-retinoic acid could activate the phosphorylation of Smad2/3, inducing apoptosis of dermal papilla cells . Low levels of phosphorylated Smad2/3 delay the hair cycle in Tsukushi-null mutant mice . SMAD7 as an inhibitor of TGFβ signaling, was up-regulated in mice with alopecia , which was consistent with our findings in giant pandas with alopecia. Overexpression of Smad7 was responsible for aberrant hair follicle morphogenesis and hyperproliferation in the epidermis, leading to severe pathological alterations in transgenic mice epithelial tissues . In addition, CDKN2B (cyclin dependent kinase inhibitor 2B) is a TGFβ target gene and encodes the cell-cycle inhibitor protein, and may be involved in hair follicle cycle . Our data showed that several hair growth-related genes in TGFβ signaling were abnormally expressed in male giant pandas with alopecia.
Mapk signaling pathway
MAPK (mitogen-activated protein kinase) comprises at least three groups, ERK (extracellular signal-related kinases), JNK (Jun amino-terminal kinases), and p38, and MAPK is involved in the regulation of normal cell proliferation, survival, differentiation, and migration . Corticotropin‐releasing hormone has been reported to induce alopecia via the MAPK signaling pathway . MAPK was suggested to regulate the quiescence of hair follicle stem cells and control the hair cycle . IKBKG (inhibitor of nuclear factor kappa B kinase regulatory subunit gamma, also known as NEMO, IKKγ) is a downstream messenger of MAPK and necessary for full activation of NF-κB . Mice lacking IKKγ or IKKβ develop severe inflammatory skin diseases, suggesting a significant role of NF-κB and IKKγ . The IKBKG gene is found mutated in hypohidrotic ectodermal dysplasia, accompanied by variable degrees of alopecia . ERK signaling is involved in the proliferation of the hair matrix cell and dermal papilla cell [71, 72]. DUSP4 (dual specificity phosphatase 4) and DUSP7 (dual specificity phosphatase 7) were upregulated in giant pandas with alopecia. DUSP4 and DUSP7 are two members of mitogen-activated protein kinase phosphatases family and are involved in dephosphorylation of ERK and negatively regulate MAPK [73, 74]. JNK functions in the proliferation and differentiation of bulge hair follicle stem cells. Hyperactivation of the JNK was associated with alopecia and topical application of the JNK inhibitor reverted faster hair growth [75, 76]. RAC1 (Rac family small GTPase 1) is associated with JNK to initiate various cellular responses and maintains the differentiated state of hair follicle keratinocytes . The deletion of RAC1 resulted in the rapid depletion of stem cells in the adult mouse epidermis . RAC1 functions in keratinocyte proliferation and migration, and the inhibition of RAC1 impairs the proliferation of keratinocytes . The down-regulation of RAC1 may be the cause of alopecia in giant pandas.
PI3K-Akt signaling pathway
Inhibition of PI3K or Akt noticeably suppressed hair follicle regeneration, suggesting a role of the PI3K-Akt signaling pathway in hair follicle de novo regeneration . The activation of PI3K signaling can promote the proliferation and migration of dermal papilla cells . Cai et al. found that lncRNA5322 could promote the proliferation and differentiation of hair follicle stem cells via the PI3K‑AKT signaling pathway . Liu et al. found that in Rex rabbits ocu-miR-205 altered PI3K-Akt, Wnt, and Notch signaling and promoted hair follicle transition from the growth phase to the regression and resting phase . OSM (oncostatin M) and IL6R (interleukin 6 receptor) were significantly up-regulated in giant pandas with alopecia. OSM, belonging to IL-6 family, negatively regulates hair growth and increased expression of OSM might contribute to alopecia [83, 84]. Kwack et al. found that IL-6 was upregulated in dermal papilla cells from balding patients . RXRA (retinoid X receptor alpha) is highly expressed in skin and in the hair follicle outer root sheath . The mutant mice without RXRA developed hair follicle degeneration, then alopecia . Blocking retinoic acid signaling delayed anagen initiation of hair follicle, while increasing retinol accelerated the transition from telogen to anagen . Therefore, we conclude that the PI3K-Akt signaling pathway might also be a candidate pathway in giant pandas with alopecia.
We investigated the blood transcriptome expression profiles of four alopecia pandas and seven healthy giant pandas. A total of 458 up-regulated DEGs and 211 down-regulated DEGs were identified. We obtained 28 hair growth-related genes with altered expression, and identified three hub genes NOTCH1, SMAD3, and TGFB1 in our PPI analysis. Five hair growth-related signaling pathways with abnormal expression were identified, being Notch, Wnt, TGF-β, Mapk, and PI3K-Akt. In conclusion, giant panda alopecia may be complex, involving abnormal expression of multiple genes and pathways. Our study provides a foundation for further studies in prevention and treatment strategies for giant pandas with alopecia.
Availability of data and materials
Raw sequence data are accessible at NCBI under the BioProject accession number PRJNA747307 (https://www.ncbi.nlm.nih.gov/sra/PRJNA747307).
Differentially expressed genes
Transcripts per million
Principal components analysis
False discovery rate
Kyoto Encyclopedia of Genes and Genomes
Schneider MR, Schmidt-Ullrich R, Paus R. The hair follicle as a dynamic miniorgan. Curr Biol. 2009;19(3):R132–42.
Alonso L, Fuchs E. The hair cycle. J Cell Sci. 2006;119(3):391–3.
Powell BC, Passmore EA, Nesci A, Dunn SM. The Notch signalling pathway in hair growth. Mech Dev. 1998;78(1–2):189–92.
Wu ZY, Zhu YL, Liu HL, Liu GY, Li FC: Wnt10b promotes hair follicles growth and dermal papilla cells proliferation via Wnt/beta-Catenin signaling pathway in Rex rabbits. Biosci Rep. 2020;40(2):BSR20191248.
Song LL, Cui Y, Yu SJ, Liu PG, He JF. TGF- and HSP70 profiles during transformation of yak hair follicles from the anagen to catagen stage. J Cell Physiol. 2019;234(9):15638–46.
Ozturk OA, Pakula H, Chmielowiec J, Qi JJ, Stein S, Lan LX, Sasaki Y, Rajewsky K, Birchmeier W. Gab1 and mapk signaling are essential in the hair cycle and hair follicle stem cell quiescence. Cell Rep. 2015;13(3):561–72.
Chen Y, Fan ZM, Wang XX, Mo MH, Zeng SB, Xu RH, Wang XS, Wu YJ. PI3K/Akt signaling pathway is essential for de novo hair follicle regeneration. Stem Cell Res Ther. 2020;11(1):144.
Kroeker R, Lee GH, Bellanca RU, Thom JP, Worlein JM. Prior facility affects alopecia in adulthood for rhesus macaques. Am J Primatol. 2017;79(1):1–9.
Hayashi A, Mikami Y, Miyamoto K, Kamada N, Sato T, Mizuno S, Naganuma M, Teratani T, Aoki R, Fukuda S, et al. Intestinal dysbiosis and biotin deprivation induce alopecia through overgrowth of lactobacillus murinus in mice. Cell Rep. 2017;20(7):1513–24.
Nicolau A, Lemberger K, Mosca M, Leclerc A, Lecu A, Pin D. Clinical and histopathological aspects of an alopecia syndrome in captive Andean bears (Tremarctos ornatus). Vet Dermatol. 2018;29(3):234.
Bowen L, Miles AK, Stott J, Waters S, Atwood T. Enhanced biological processes associated with alopecia in polar bears (Ursus maritimus). Sci Total Environ. 2015;529:114–20.
Lutz CK, Menard MT, Rosenberg K, Meyer JS, Novak MA. Alopecia in rhesus macaques (Macaca mulatta): association with pregnancy and chronic stress. J Med Primatol. 2019;48(4):251–6.
Novak MA, Meyer JS. Alopecia: possible causes and treatments, particularly in captive nonhuman primates. Comp Med. 2009;59(1):18–26.
Anzai A, Wang EHC, Lee EY, Aoki V, Christiano AM. Pathomechanisms of immune-mediated alopecia. Int Immunol. 2019;31(7):439–47.
Beisner BA, Isbell LA. Factors influencing hair loss among female captive rhesus macaques (Macaca mulatta). Appl Anim Behav Sci. 2009;119(1–2):91–100.
Prokop P, Randler C. Biological predispositions and individual differences in human attitudes toward animals. R.R.N. Alves, U.P. Alburquerque (Eds.), Ethnozoology. Amsterdam: Elsevier; 2018. p. 447–66.
Liew CC, Ma J, Tang HC, Zheng R, Dempsey AA. The peripheral blood transcriptome dynamically reflects system wide biology: a potential diagnostic tool. J Lab Clin Med. 2006;147(3):126–32.
Mohr S, Liew CC. The peripheral-blood transcriptome: new insights into disease and risk assessment. Trends Mol Med. 2007;13(10):422–32.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, Genome Project Data P. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
Reimand J, Arak T, Adler P, Kolberg L, Reisberg S, Peterson H, Vilo J. g:Profiler-a web server for functional interpretation of gene lists (2016 update). Nucleic Acids Res. 2016;44(W1):W83–9.
Xie C, Mao XZ, Huang JJ, Ding Y, Wu JM, Dong S, Kong L, Gao G, Li CY, Wei LP. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39:W316–22.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.
Smoot ME, Ono K, Ruscheinski J, Wang PL, Ideker T. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011;27(3):431–2.
Chin CH, Chen SH, Wu HH, Ho CW, Ko MT, Lin CY. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 2014;84 Suppl 4:S11.
Shen HB, Li CW, He M, Huang Y, Wang J, Wang ML, Yue BS, Zhang XY. Immune profiles of male giant panda (Ailuropoda melanoleuca) during the breeding season. BMC Genomics. 2021;22(1):143.
Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden TL. Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. Bmc Bioinformatics. 2012;13:134.
Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C-T method. Nat Protoc. 2008;3(6):1101–8.
Buffoli B, Rinaldi F, Labanca M, Sorbellini E, Trink A, Guanziroli E, Rezzani R, Rodella LF. The human hair: from anatomy to physiology. Int J Dermatol. 2014;53(3):331–41.
Lv XF, Chen L, He SG, Liu CX, Han B, Liu ZL, Yusupu M, Blair H, Kenyon P, Morris S, et al. Effect of nutritional restriction on the hair follicles development and skin transcriptome of Chinese merino sheep. Animals. 2020;10(6):1058.
Zhou Q, Song YJ, Zheng QL, Han R, Cheng H. Expression profile analysis of dermal papilla cells mRNA in response to WNT10B treatment. Exp Ther Med. 2020;19(2):1017–23.
Zhang YL, Xia SZ, Wang TC, Wang SH, Yuan D, Li F, Wang X. Chi-miR-30b-5p inhibits dermal papilla cells proliferation by targetingCaMKII delta gene in cashmere goat. Bmc Genomics. 2020;21(1):430.
Plowman JE, Harland DP, Ganeshan S, Woods JL, van Shaijik B, Deb-Choudhury S, Thomas A, Clerens S, Scobie DR. The proteomics of wool fibre morphogenesis. J Struct Biol. 2015;191(3):341–51.
Sheng YY, Qi SS, Hu RM, Zhao J, Rui WL, Miao Y, Ma JW, Yang QP. Identification of blood microRNA alterations in patients with severe active alopecia areata. J Cell Biochem. 2019;120(9):14421–30.
Kim CS, Ding XL, Allmeroth K, Biggs LC, Kolenc OI, L’Hoest N, Chacon-Martinez CA, Edlich-Muth C, Giavalisco P, Quinn KP, et al. Glutamine metabolism controls stem cell fate reversibility and long-term maintenance in the hair follicle. Cell Metab. 2020;32(4):629.
Nocelli C, Cappelli K, Capomaccio S, Pascucci L, Mercati F, Pazzaglia I, Mecocci S, Antonini M, Renieri C. Shedding light on cashmere goat hair follicle biology: from morphology analyses to transcriptomic landascape. Bmc Genomics. 2020;21(1):458.
Lv XY, Chen WH, Sun W, Hussain Z, Chen L, Wang SH, Wang JY. Expression profile analysis to identify circular RNA expression signatures in hair follicle of Hu sheep lambskin. Genomics. 2020;112(6):4454–62.
Wu ZY, Fu YH, Cao JH, Yu M, Tang XH, Zhao SH. Identification of differentially expressed miRNAs between white and black hair follicles by RNA-Sequencing in the Goat (Capra hircus). Int J Mol Sci. 2014;15(6):9531–45.
Rishikaysh P, Dev K, Diaz D, Qureshi WMS, Filip S, Mokry J. Signaling involved in hair follicle morphogenesis and development. Int J Mol Sci. 2014;15(1):1647–70.
Aubin-Houzelstein G.Notch signaling and the developing hair folliclE. In: Notch Signaling in Embryology and Cancer. Edited by Reichrath J, Reichrath S, vol. 727; 2012: 142–160.
Vauclair S, Nicolas M, Barrandon Y, Radtke F. Notch1 is essential for postnatal hair follicle development and homeostasis. Dev Biol. 2005;284(1):184–93.
Uyttendaele H, Panteleyev AA, de Berker D, Tobin DJ, Christiano AM. Activation of Notch1 in the hair follicle leads to cell-fate switch and Mohawk alopecia. Differentiation. 2004;72(8):396–409.
Kojika S, Griffin JD. Notch receptors and hematopoiesis. Exp Hematol. 2001;29(9):1041–52.
Zhan T, Rindtorff N, Boutros M. Wnt signaling in cancer. Oncogene. 2017;36(11):1461–73.
Gao C, Chen YG. Dishevelled: The hub of Wnt signaling. Cell Signal. 2010;22(5):717–27.
Choi BY. Hair-growth potential of ginseng and its major metabolites: a review on its molecular mechanisms. Int J Mol Sci. 2018;19(9):2703.
Ge MX, Liu CC, Li L, Lan MM, Yu YY, Gu LJ, Su Y, Zhang K, Zhang YY, Wang TT, et al. miR-29a/b1 inhibits hair follicle stem cell lineage progression by spatiotemporally suppressing WNT and BMP signaling. Cell Reports. 2019;29(8):2489.
Bao PJ, Luo JY, Liu YB, Chu M, Ren QM, Guo X, Tang BL, Ding XZ, Qiu Q, Pan HP, et al. The seasonal development dynamics of the yak hair cycle transcriptome. BMC Genomics. 2020;21(1):355.
Zhang XL, Bao Q, Jia CJ, Li C, Chang YF, Wu XY, Liang CN, Bao PJ, Yan P. Genome-wide detection and sequence conservation analysis of long non-coding RNA during hair follicle cycle of yak. Bmc Genomics. 2020;21(1):681.
Chu QQ, Cai LY, Fu Y, Chen X, Yan ZP, Lin X, Zhou GX, Han H, Widelitz RB, Chuong CM, et al. Dkk2/Frzb in the dermal papillae regulates feather regeneration. Dev Biol. 2014;387(2):167–78.
Xie B, Chen MT, Ding PH, Lei L, Zhang X, Zhu D, Zou YJ, Deng ZL, Sun GM, Li J et al: Induction of dermal fibroblasts into dermal papilla cell-like cells in hydrogel microcapsules for enhanced hair follicle regeneration. Appl Mater Today. 2020;21.
Millar SE, Willert K, Salinas PC, Roelink H, Nusse R, Sussman DJ, Barsh GS. WNT signaling in the control of hair growth and structure. Dev Biol. 1999;207(1):133–49.
Lee SH, Seo SH, Lee DH, Pi LQ, Lee WS, Choi KY. Targeting of CXXC5 by a competing peptide stimulates hair regrowth and wound-induced hair neogenesis. J Investig Dermatol. 2017;137(11):2260–9.
Gafter-Gvili A, Sredni B, Gal R, Gafter U, Kalechman Y. Cyclosporin A-induced hair growth in mice is associated with inhibition of calcineurin-dependent activation of NFAT in follicular keratinocytes. Am J Physiol Cell Physiol. 2003;284(6):C1593–603.
Cetera M, Leybova L, Woo FW, Deans M, Devenport D. Planar cell polarity-dependent and independent functions in the emergence of tissue-scale hair follicle patterns. Dev Biol. 2017;428(1):188–203.
Devenport D, Fuchs E. Planar polarization in embryonic epidermis orchestrates global asymmetric morphogenesis of hair follicles. Nat Cell Biol. 2008;10(11):1257-U1212.
Kahata K, Dadras MS, Moustakas A. TGF-beta family signaling in epithelial differentiation and epithelial-mesenchymal transition. Cold Spring Harb Perspect Biol. 2018;10(1):a022194.
Begum S, Lee MR, Gu LJ, Hossain J, Sung CK. Exogenous stimulation with Eclipta alba promotes hair matrix keratinocyte proliferation and downregulates TGF-beta 1 expression in nude mice. Int J Mol Med. 2015;35(2):496–502.
Liu X, Alexander V, Vijayachandra K, Bhogte E, Diamond I, Glick A. Conditional epidermal expression of TGF beta 1 blocks neonatal lethality but causes a reversible hyperplasia and alopecia. Proc Natl Acad Sci USA. 2001;98(16):9139–44.
Nan WX, Li GY, Si HZ, Lou YJ, Wang DY, Guo R, Zhang HH. All-trans-retinoic acid inhibits mink hair follicle growth via inhibiting proliferation and inducing apoptosis of dermal papilla cells through TGF-beta 2/Smad2/3 pathway. Acta Histochem. 2020;122(7):151603.
Niimori D, Kawano R, Felemban A, Niimori-Kita K, Tanaka H, Ihn H, Ohta K. Tsukushi controls the hair cycle by regulating TGF-beta 1 signaling. Dev Biol. 2012;372(1):81–7.
He W, Li AG, Wang DY, Han SH, Zheng B, Goumans MJ, ten Dijke P, Wang XJ. Overexpression of Smad7 results in severe pathological alterations in multiple epithelial tissues. EMBO J. 2002;21(11):2580–90.
Lefever T, Pedersen E, Basse A, Paus R, Quondamatteo F, Stanley AC, Langbein L, Wu XW, Wehland J, Lommel S, et al. N-WASP is a novel regulator of hair-follicle cycling that controls antiproliferative TGF beta pathways. J Cell Sci. 2010;123(1):128–40.
Mishra S, Tripathi A, Chaudhari BP, Dwivedi PD, Pandey HP, Das M. Deoxynivalenol induced mouse skin cell proliferation and inflammation via MAPK pathway. Toxicol Appl Pharmacol. 2014;279(2):186–97.
Nam YJ, Lee EY, Choi EJ, Kang S, Kim J, Choi YS, Kim DH, An JH, Han I, Lee S, et al. CRH receptor antagonists from Pulsatilla chinensis prevent CRH-induced premature catagen transition in human hair follicles. J Cosmet Dermatol. 2020;19(11):3058–66.
Ma CA, Wang HY, Temmerman S, Zhao YG, Wu LM, Hornung RL, Wara D, Jain A. Dendritic cells from humans with hypomorphic mutations in IKBKG/NEMO have impaired mitogen-activated protein kinase activity. Hum Mutat. 2011;32(3):318–24.
Page A, Navarro M, Garin M, Perez P, Casanova ML, Moreno R, Jorcano JL, Cascallana JL, Bravo A, Ramirez A. IKK beta leads to an inflammatory skin disease resembling interface dermatitis. J Investig Dermatol. 2010;130(6):1598–610.
Duverger O, Morasso MI. To grow or not to grow: hair morphogenesis and human genetic hair disorders. Semin Cell Dev Biol. 2014;25:22–33.
Tang P, Wang XE, Zhang M, Huang SM, Lin CX, Yan F, Deng Y, Zhang L, Zhang L. Activin B stimulates mouse vibrissae growth and regulates cell proliferation and cell cycle progression of hair matrix cells through ERK signaling. Int J Mol Sci. 2019;20(4):853.
Huang HC, Lin H, Huang MC. Lactoferrin promotes hair growth in mice and increases dermal papilla cell proliferation through Erk/Akt and Wnt signaling pathways. Arch Dermatol Res. 2019;311(5):411–20.
Zhang R, Wang G, Zhang PF, Zhang J, Huang YX, Lu YM, Da W, Sun Q, Zhu JS. Sanguinarine inhibits growth and invasion of gastric cancer cells via regulation of the DUSP4/ERK pathway. J Cell Mol Med. 2017;21(6):1117–27.
Peng WX, Huang JG, Yang L, Gong AH, Mo YY. Linc-RoR promotes MAPK/ERK signaling and confers estrogen-independent growth of breast cancer. Mol Cancer. 2017;16:161.
Li Z, Ryu SW, Lee J, Choi K, Kim S, Choi C. Protopanaxatirol type ginsenoside Re promotes cyclic growth of hair follicles via inhibiting transforming growth factor beta signaling cascades. Biochem Biophys Res Commun. 2016;470(4):924–9.
Chovatiya GL, Sunkara RR, Roy S, Godbole SR, Waghmare SK. Context-dependent effect of sPLA(2)-IIA induced proliferation on murine hair follicle stem cells and human epithelial cancer. EBioMedicine. 2019;48:364–76.
Chrostek A, Wu XW, Quondamatteo F, Hu R, Sanecka A, Niemann C, Langbein L, Haase I, Brakebusch C. Rac1 is crucial for hair follicle integrity but is not essential for maintenance of the epidermis. Mol Cell Biol. 2006;26(18):6957–70.
Benitah SA, Frye M, Glogauer M, Watt FM. Stem cell depletion through epidermal deletion of Rac1. Science. 2005;309(5736):933–5.
Tscharntke M, Pofahl R, Chrostek-Grashoff A, Smyth N, Niessen C, Niemann C, Hartwig B, Herzog V, Klein HW, Krieg T, et al. Impaired epidermal wound healing in vivo upon inhibition or deletion of Rac1. J Cell Sci. 2007;120(8):1480–90.
Zhang HH, Su YM, Wang JT, Gao Y, Yang FH, Li GY, Shi QM. Ginsenoside Rb1 promotes the growth of mink hair follicle via PI3K/AKT/GSK-3 beta signaling pathway. Life Sci. 2019;229:210–8.
Cai BJ, Zheng YP, Ma SS, Xing Q, Wang XX, Yang B, Yin GW, Guan FX. Long non-coding RNA regulates hair follicle stem cell proliferation and differentiation through PI3K/AKT signal pathway. Mol Med Rep. 2018;17(4):5477–83.
Liu GY, Li S, Liu HL, Zhu YL, Bai LY, Sun HT, Gao SX, Jiang WX, Li FC. The functions of ocu-miR-205 in regulating hair follicle development in Rex rabbits. Bmc Dev Biol. 2020;20(1):8.
Yu M, Kissling S, Freyschmidt-Paul P, Hoffmann R, Shapiro J, McElwee KJ. Interleukin-6 cytokine family member oncostatin M is a hair-follicle-expressed factor with hair growth inhibitory properties. Exp Dermatol. 2008;17(1):12–9.
Wang ECE, Dai ZP, Ferrante AW, Drake CG, Christiano AM. A subset of TREM2(+) dermal macrophages secretes oncostatin m to maintain hair follicle stem cell quiescence and inhibit hair growth. Cell Stem Cell. 2019;24(4):654.
Kwack MH, Ahn JS, Kim MK, Kim JC, Sung YK. Dihydrotestosterone-inducible IL-6 inhibits elongation of human hair shafts by suppressing matrix cell proliferation and promotes regression of hair follicles in mice. J Investig Dermatol. 2012;132(1):43–9.
Hochfeld LM, Anhalt T, Reinbold CS, Herrera-Rivero M, Fricker N, Nothen MM, Heilmann-Heimbach S. Expression profiling and bioinformatic analyses suggest new target genes and pathways for human hair follicle related microRNAs. Bmc Dermatol. 2017;17:3.
Ghyselinck NB, Chapellier B, Calleja C, Indra AK, Li M, Messaddeq N, Mark M, Metzger D, Chambon P. Genetic dissection of retinoic acid function in epidermis physiology. Ann Dermatol Venereol. 2002;129(5):793–9.
Everts HB. Endogenous retinoids in the hair follicle and sebaceous gland. BBA-Mol Cell Biol L. 2012;1821(1):222–9.
We acknowledge Dr. Megan Price for providing language help.
This work was supported by the Open Project of Key Laboratory of State Forestry and grassland administration (KLSFGAGP2020.010) and Project of Key Laboratory of State Forestry and grassland administration (CCRCGP181913). The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
This study was carried out according to the Regulation on the Administration of Laboratory Animals (2017 Revision) published by the Ministry of Science and Technology of the People´s Republic of China. All study procedures and animal care activities were approved by the Academic and Ethics Committee of Sichuan University.
Consent for publication
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Shen, H., Li, C., He, M. et al. Whole blood transcriptome profiling identifies candidate genes associated with alopecia in male giant pandas (Ailuropoda melanoleuca). BMC Genomics 23, 297 (2022). https://doi.org/10.1186/s12864-022-08501-z
- Giant panda
- Peripheral blood
- Transcriptome profiling
- Differential gene expression