{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,8,25]],"date-time":"2026-08-25T05:19:09Z","timestamp":1787635149603,"version":"build-2736575974"},"reference-count":55,"publisher":"MDPI AG","issue":"16","license":[{"start":{"date-parts":[[2026,8,21]],"date-time":"2026-08-21T00:00:00Z","timestamp":1787270400000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"name":"Bayannaoer Research Institute Young Scientist Project","award":["2024BYNECAU004"],"award-info":[{"award-number":["2024BYNECAU004"]}]},{"DOI":"10.13039\/501100005236","name":"Chinese Universities Scientific Fund","doi-asserted-by":"crossref","award":["2026RC017"],"award-info":[{"award-number":["2026RC017"]}],"id":[{"id":"10.13039\/501100005236","id-type":"DOI","asserted-by":"crossref"}]},{"name":"National Key Research and Development Program of China","award":["2023YFF10009"],"award-info":[{"award-number":["2023YFF10009"]}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["Int. J. Mol. Sci.","IJMS"],"abstract":"<jats:p>Skin is a highly specialized barrier organ, which serves essential barrier, immune, and sensory functions and exhibits conserved structure. In addition, each species exhibits a unique skin transcriptional profile for adapting distinct environmental conditions and physiological demands. However, these transcriptional profiles across mammals remain incompletely understood. In this study, we integrated 61 publicly available skin RNA-seq datasets from eight mammalian species: human (Homo sapiens), macaque (Macaca fascicularis), mouse (Mus musculus), sheep (Ovis aries), goat (Capra hircus), donkey (Equus asinus), rabbit (Oryctolagus cuniculus), and pig (Sus scrofa). We identified a large number (2286 to 4588) of differentially expressed genes (DEGs) based on 10,504 one-to-one orthologous genes in all pairwise species comparisons. A total of 44 conserved genes were identified by integrating tissue specificity, expression variability, and mean expression level, which were enriched in epidermal development, keratinocyte differentiation, and desmosome organization. We also identified species-specific genes and transcription factors, such as STMN1 (donkey), COL1A1 (goat), and ACTN2 (mouse), which are associated with the cell cycle, extracellular matrix, and muscle contraction pathways, respectively. WGCNA further revealed goat-associated module 2 (M2) enriched in the negative regulation of angiogenesis. Collectively, this study provides a comprehensive cross-species transcriptional profile resource for mammalian skin, providing a valuable resource for understanding the regulatory plasticity underlying skin evolution across mammals.<\/jats:p>","DOI":"10.3390\/ijms27167508","type":"journal-article","created":{"date-parts":[[2026,8,24]],"date-time":"2026-08-24T11:47:19Z","timestamp":1787572039000},"page":"7508","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":0,"title":["Multispecies Transcriptomics of Mammalian Skin Reveals Conserved and Divergent Regulatory Features"],"prefix":"10.3390","volume":"27","author":[{"ORCID":"https:\/\/orcid.org\/0000-0001-8049-5565","authenticated-orcid":false,"given":"Wen-Tian","family":"Wei","sequence":"first","affiliation":[{"name":"Frontiers Science Center for Molecular Design Breeding (MOE), State Key Laboratory of Animal Biotech Breeding, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Chen","family":"Zhou","sequence":"additional","affiliation":[{"name":"Frontiers Science Center for Molecular Design Breeding (MOE), State Key Laboratory of Animal Biotech Breeding, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Qi-Xuan","family":"Huang","sequence":"additional","affiliation":[{"name":"Frontiers Science Center for Molecular Design Breeding (MOE), State Key Laboratory of Animal Biotech Breeding, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Zhong-Hua","family":"Ning","sequence":"additional","affiliation":[{"name":"State Key Laboratory of Animal Nutrition, National Engineering Laboratory for Animal Breeding, Department of Animal Genetics and Breeding, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Yi-Lin","family":"Bai","sequence":"additional","affiliation":[{"name":"Key Laboratory of Innovative Utilization of Local Cattle and Sheep Germplasm Resources (Co-Construction by Ministry and Province), Ministry of Agriculture and Rural Affairs, Zhengzhou University, Zhengzhou 450001, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Yue-Yu","family":"Bai","sequence":"additional","affiliation":[{"name":"Key Laboratory of Innovative Utilization of Local Cattle and Sheep Germplasm Resources (Co-Construction by Ministry and Province), Ministry of Agriculture and Rural Affairs, Zhengzhou University, Zhengzhou 450001, China"},{"name":"Institute of Livestock Breeding and Reproduction, Henan University, Zhengzhou 450046, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Song-Song","family":"Xu","sequence":"additional","affiliation":[{"name":"Frontiers Science Center for Molecular Design Breeding (MOE), State Key Laboratory of Animal Biotech Breeding, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"1968","published-online":{"date-parts":[[2026,8,21]]},"reference":[{"key":"ref_1","doi-asserted-by":"crossref","first-page":"a005124","DOI":"10.1101\/cshperspect.a005124","article-title":"Cell-extracellular matrix interactions in normal and diseased skin","volume":"3","author":"Watt","year":"2011","journal-title":"Cold Spring Harb. Perspect. Biol."},{"key":"ref_2","doi-asserted-by":"crossref","first-page":"3543","DOI":"10.1002\/ar.24908","article-title":"Comparison of vertebrate skin structure at class level: A review","volume":"305","author":"Akat","year":"2022","journal-title":"Anat. Rec."},{"key":"ref_3","doi-asserted-by":"crossref","first-page":"834","DOI":"10.1038\/nature05659","article-title":"Scratching the surface of skin development","volume":"445","author":"Fuchs","year":"2007","journal-title":"Nature"},{"key":"ref_4","doi-asserted-by":"crossref","first-page":"315","DOI":"10.1002\/path.5589","article-title":"A developmental basis for the anatomical diversity of dermis in homeostasis and wound repair","volume":"253","author":"Usansky","year":"2021","journal-title":"J. Pathol."},{"key":"ref_5","doi-asserted-by":"crossref","first-page":"1179","DOI":"10.1111\/brv.12151","article-title":"Resting no more: Re-defining telogen, the maintenance stage of the hair growth cycle","volume":"90","author":"Geyfman","year":"2015","journal-title":"Biol. Rev. Camb. Philos. Soc."},{"key":"ref_6","doi-asserted-by":"crossref","first-page":"eadg9619","DOI":"10.1126\/sciadv.adg9619","article-title":"Transient agonism of the sonic hedgehog pathway triggers a permanent transition of skin appendage fate in the chicken embryo","volume":"9","author":"Cooper","year":"2023","journal-title":"Sci. Adv."},{"key":"ref_7","doi-asserted-by":"crossref","unstructured":"Aldea, D., Kokalari, B., Atsuta, Y., Dingwall, H.L., Zheng, Y., Nace, A., Cotsarelis, G., and Kamberov, Y.G. (2023). Differential modularity of the mammalian Engrailed 1 enhancer network directs sweat gland development. PLoS Genet., 19.","DOI":"10.1371\/journal.pgen.1010614"},{"key":"ref_8","doi-asserted-by":"crossref","first-page":"343","DOI":"10.1038\/nature10532","article-title":"The evolution of gene expression levels in mammalian organs","volume":"478","author":"Brawand","year":"2011","journal-title":"Nature"},{"key":"ref_9","doi-asserted-by":"crossref","first-page":"1593","DOI":"10.1126\/science.1228186","article-title":"Evolutionary dynamics of gene and isoform regulation in mammalian tissues","volume":"338","author":"Merkin","year":"2012","journal-title":"Science"},{"key":"ref_10","doi-asserted-by":"crossref","first-page":"53","DOI":"10.1101\/gr.237636.118","article-title":"A quantitative framework for characterizing the evolutionary history of mammalian gene expression","volume":"29","author":"Chen","year":"2019","journal-title":"Genome Res."},{"key":"ref_11","doi-asserted-by":"crossref","first-page":"505","DOI":"10.1038\/s41586-019-1338-5","article-title":"Gene expression across mammalian organ development","volume":"571","author":"Halbert","year":"2019","journal-title":"Nature"},{"key":"ref_12","doi-asserted-by":"crossref","first-page":"777","DOI":"10.1016\/S0029-6465(22)02423-9","article-title":"Skin anatomy, physiology, and pathophysiology","volume":"34","author":"Wysocki","year":"1999","journal-title":"Nurs. Clin. N. Am."},{"key":"ref_13","doi-asserted-by":"crossref","first-page":"207","DOI":"10.1038\/nrm2636","article-title":"Epidermal Homeostasis: A balancing act of stem cells in the skin","volume":"10","author":"Blanpain","year":"2009","journal-title":"Nat. Rev. Mol. Cell Biol."},{"key":"ref_14","doi-asserted-by":"crossref","first-page":"498","DOI":"10.1111\/apha.12231","article-title":"Skin Temperature: Its role in thermoregulation","volume":"210","author":"Romanovsky","year":"2014","journal-title":"Acta Physiol."},{"key":"ref_15","doi-asserted-by":"crossref","first-page":"1429","DOI":"10.1038\/labinvest.2015.116","article-title":"Myofibroblasts contribute to but are not necessary for wound contraction","volume":"95","author":"Ibrahim","year":"2015","journal-title":"Lab. Investig."},{"key":"ref_16","doi-asserted-by":"crossref","first-page":"4047","DOI":"10.1091\/mbc.e11-06-0582","article-title":"\u03b2-Actin specifically controls cell growth, migration, and the G-Actin pool","volume":"22","author":"Bunnell","year":"2011","journal-title":"Mol. Biol. Cell"},{"key":"ref_17","doi-asserted-by":"crossref","first-page":"E4320","DOI":"10.1073\/pnas.1519197113","article-title":"Vimentin coordinates fibroblast proliferation and keratinocyte differentiation in wound healing via TGF-\u03b2-Slug signaling","volume":"113","author":"Cheng","year":"2016","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"ref_18","doi-asserted-by":"crossref","first-page":"1775","DOI":"10.1091\/mbc.12.6.1775","article-title":"Complete cytolysis and neonatal lethality in keratin 5 knockout mice reveal its fundamental role in skin integrity and in epidermolysis bullosa simplex","volume":"12","author":"Peters","year":"2001","journal-title":"Mol. Biol. Cell"},{"key":"ref_19","doi-asserted-by":"crossref","first-page":"1076","DOI":"10.1038\/ncb1201-1076","article-title":"Desmo plakin is essential in epidermal sheet formation","volume":"3","author":"Vasioukhin","year":"2001","journal-title":"Nat. Cell Biol."},{"key":"ref_20","doi-asserted-by":"crossref","first-page":"40355","DOI":"10.1074\/jbc.M506692200","article-title":"Mechanisms of plakoglobin-dependent adhesion: Desmosome-specific functions in assembly and regulation by epidermal growth factor receptor","volume":"280","author":"Yin","year":"2005","journal-title":"J. Biol. Chem."},{"key":"ref_21","doi-asserted-by":"crossref","first-page":"3303","DOI":"10.1242\/jcs.00636","article-title":"Lack of plakophilin 1 increases keratinocyte migration and reduces desmosome sability","volume":"116","author":"South","year":"2003","journal-title":"J. Cell Sci."},{"key":"ref_22","doi-asserted-by":"crossref","first-page":"577","DOI":"10.1016\/S0092-8674(00)80868-5","article-title":"Endless Forms: The evolution of gene regulation and morphological diversity","volume":"101","author":"Carroll","year":"2000","journal-title":"Cell"},{"key":"ref_23","doi-asserted-by":"crossref","first-page":"691","DOI":"10.1016\/j.cell.2013.01.016","article-title":"Modeling recent human evolution in mice by expression of a selected EDAR variant","volume":"152","author":"Kamberov","year":"2013","journal-title":"Cell"},{"key":"ref_24","doi-asserted-by":"crossref","first-page":"184","DOI":"10.1016\/j.ydbio.2005.05.018","article-title":"Notch1 is essential for postnatal hair follicle development and homeostasis","volume":"284","author":"Vauclair","year":"2005","journal-title":"Dev. Biol."},{"key":"ref_25","doi-asserted-by":"crossref","first-page":"721","DOI":"10.1016\/j.bbrc.2014.12.039","article-title":"Trps1 deficiency inhibits the morphogenesis of secondary hair follicles via decreased noggin expression","volume":"456","author":"Sun","year":"2015","journal-title":"Biochem. Biophys. Res. Commun."},{"key":"ref_26","doi-asserted-by":"crossref","first-page":"379","DOI":"10.1242\/dev.00201","article-title":"\u2018Cyclic Alopecia\u2019 in Msx2 mutants: Defects in hair cycling and hair shaft differentiation","volume":"130","author":"Ma","year":"2003","journal-title":"Development"},{"key":"ref_27","doi-asserted-by":"crossref","first-page":"1243","DOI":"10.1038\/jid.2008.369","article-title":"Hair follicle stem cell-specific PPAR\u03b3 deletion causes scarring alopecia","volume":"129","author":"Karnik","year":"2009","journal-title":"J. Investig. Dermatol."},{"key":"ref_28","doi-asserted-by":"crossref","first-page":"861","DOI":"10.1038\/nature00820","article-title":"TNF-mediated inflammatory skin disease in mice with epidermis-specific deletion of IKK2","volume":"417","author":"Pasparakis","year":"2002","journal-title":"Nature"},{"key":"ref_29","doi-asserted-by":"crossref","first-page":"491","DOI":"10.1016\/j.stem.2016.06.020","article-title":"Dnmt3a and Dnmt3b associate with enhancers to regulate human epidermal stem cell homeostasis","volume":"19","author":"Rinaldi","year":"2016","journal-title":"Cell Stem Cell"},{"key":"ref_30","doi-asserted-by":"crossref","first-page":"5290","DOI":"10.1038\/s44318-025-00519-9","article-title":"Specific and redundant roles for Gli2 and Gli3 in establishing cell fate during murine hair follicle development","volume":"44","author":"Gozum","year":"2025","journal-title":"EMBO J."},{"key":"ref_31","doi-asserted-by":"crossref","first-page":"11135","DOI":"10.1038\/s41598-022-14441-1","article-title":"Assessment of hair and cashmere properties and their genetic background of several goat breeds in southwest China","volume":"12","author":"Saleh","year":"2022","journal-title":"Sci. Rep."},{"key":"ref_32","doi-asserted-by":"crossref","first-page":"182","DOI":"10.1016\/j.jdin.2020.08.001","article-title":"Single-Cell transcriptomics in dermatology","volume":"1","author":"Deutsch","year":"2020","journal-title":"JAAD Int."},{"key":"ref_33","doi-asserted-by":"crossref","first-page":"R132","DOI":"10.1016\/j.cub.2008.12.005","article-title":"The hair follicle as a dynamic miniorgan","volume":"19","author":"Schneider","year":"2009","journal-title":"Curr. Biol."},{"key":"ref_34","doi-asserted-by":"crossref","first-page":"409","DOI":"10.1172\/JCI11317","article-title":"Control of hair growth and follicle size by VEGF-mediated angiogenesis","volume":"107","author":"Yano","year":"2001","journal-title":"J. Clin. Investig."},{"key":"ref_35","doi-asserted-by":"crossref","first-page":"878","DOI":"10.1046\/j.1523-1747.1999.00791.x","article-title":"Hair cycle-dependent changes in adrenergic skin innervation, and hair growth modulation by adrenergic drugs","volume":"113","author":"Botchkarev","year":"1999","journal-title":"J. Investig. Dermatol."},{"key":"ref_36","doi-asserted-by":"crossref","first-page":"703","DOI":"10.1111\/exd.14155","article-title":"Hair follicle immune privilege and its collapse in alopecia areata","volume":"29","author":"Bertolini","year":"2020","journal-title":"Exp. Dermatol."},{"key":"ref_37","doi-asserted-by":"crossref","first-page":"1155","DOI":"10.1152\/physrev.00044.2003","article-title":"Melanin pigmentation in mammalian skin and its hormonal regulation","volume":"84","author":"Slominski","year":"2004","journal-title":"Physiol. Rev."},{"key":"ref_38","doi-asserted-by":"crossref","first-page":"924","DOI":"10.1038\/nm.3194","article-title":"Direct migration of follicular melanocyte stem cells to the epidermis after wounding or UVB irradiation is dependent on Mc1r signaling","volume":"19","author":"Chou","year":"2013","journal-title":"Nat. Med."},{"key":"ref_39","doi-asserted-by":"crossref","first-page":"13502","DOI":"10.1038\/s41598-017-13986-w","article-title":"Skin transcriptome reveals the intrinsic molecular mechanisms underlying hair follicle cycling in Cashmere goats under natural and shortened photoperiod conditions","volume":"7","author":"Yang","year":"2017","journal-title":"Sci. Rep."},{"key":"ref_40","doi-asserted-by":"crossref","unstructured":"Ben Abdallah, H., Daugaard, C., Sch\u00f8sler, L., Rasmussen, M.K., Kibsgaard, L., H\u00f8gsberg, T., Johansen, C., and Bertelsen, T. (2026). Transcriptomic profiling reveals distinct molecular signatures among lesion types in hidradenitis suppurativa. Front. Immunol., 16.","DOI":"10.3389\/fimmu.2025.1715474"},{"key":"ref_41","doi-asserted-by":"crossref","first-page":"eaaw7317","DOI":"10.1126\/science.aaw7317","article-title":"Conservation, acquisition, and functional impact of sex-biased gene expression in mammals","volume":"365","author":"Naqvi","year":"2019","journal-title":"Science"},{"key":"ref_42","doi-asserted-by":"crossref","first-page":"18977","DOI":"10.1038\/s41598-025-01187-9","article-title":"Transcriptome sequencing reveals the expression profiles of lncRNAs and mRNAs in goat skin tissues with different types of wool goats","volume":"15","author":"Zhao","year":"2025","journal-title":"Sci. Rep."},{"key":"ref_43","doi-asserted-by":"crossref","first-page":"23624","DOI":"10.1038\/s41598-025-08837-y","article-title":"RNA-Seq analysis identifies key genes and signaling pathways involved in androgen promotion of sebaceous gland proliferation in Hetian sheep","volume":"15","author":"Feng","year":"2025","journal-title":"Sci. Rep."},{"key":"ref_44","doi-asserted-by":"crossref","first-page":"458","DOI":"10.1111\/eva.12887","article-title":"The fine-scale genetic structure and selection signals of Chinese indigenous pigs","volume":"13","author":"Huang","year":"2020","journal-title":"Evol. Appl."},{"key":"ref_45","doi-asserted-by":"crossref","unstructured":"Zhao, B., Yu, Y., Sun, S., Cai, J., Bao, Z., Chen, Y., and Wu, X. (2024). Integration analysis of transcriptome sequencing and whole-genome resequencing reveal wool quality-associated key genes in Zhexi angora rabbits. Vet. Sci., 11.","DOI":"10.3390\/vetsci11120651"},{"key":"ref_46","doi-asserted-by":"crossref","first-page":"368","DOI":"10.1111\/age.13196","article-title":"An analysis of skin thickness in the dezhou donkey population and identification of candidate genes by RNA-Seq","volume":"53","author":"Wang","year":"2022","journal-title":"Anim. Genet."},{"key":"ref_47","doi-asserted-by":"crossref","first-page":"e70078","DOI":"10.1002\/imt2.70078","article-title":"Fastp 1.0: An ultra-fast all-round tool for FASTQ Data Quality Control and Preprocessing","volume":"4","author":"Chen","year":"2025","journal-title":"iMeta"},{"key":"ref_48","doi-asserted-by":"crossref","first-page":"15","DOI":"10.1093\/bioinformatics\/bts635","article-title":"STAR: Ultrafast universal RNA-Seq aligner","volume":"29","author":"Dobin","year":"2013","journal-title":"Bioinformatics"},{"key":"ref_49","doi-asserted-by":"crossref","first-page":"e108","DOI":"10.1093\/nar\/gkt214","article-title":"The subread aligner: Fast, accurate and scalable read mapping by seed-and-vote","volume":"41","author":"Liao","year":"2013","journal-title":"Nucleic Acids Res."},{"key":"ref_50","doi-asserted-by":"crossref","first-page":"238","DOI":"10.1186\/s13059-019-1832-y","article-title":"OrthoFinder: Phylogenetic orthology inference for comparative genomics","volume":"20","author":"Emms","year":"2019","journal-title":"Genome Biol."},{"key":"ref_51","doi-asserted-by":"crossref","first-page":"550","DOI":"10.1186\/s13059-014-0550-8","article-title":"Moderated estimation of fold change and dispersion for RNA-Seq data with DESeq2","volume":"15","author":"Love","year":"2014","journal-title":"Genome Biol."},{"key":"ref_52","doi-asserted-by":"crossref","first-page":"284","DOI":"10.1089\/omi.2011.0118","article-title":"ClusterProfiler: An R package for comparing biological themes among gene clusters","volume":"16","author":"Yu","year":"2012","journal-title":"Omics J. Integr. Biol."},{"key":"ref_53","doi-asserted-by":"crossref","first-page":"1363","DOI":"10.1101\/gr.240663.118","article-title":"Benchmark and integration of resources for the estimation of human transcription factor activities","volume":"29","author":"Holland","year":"2019","journal-title":"Genome Res."},{"key":"ref_54","doi-asserted-by":"crossref","first-page":"vbac016","DOI":"10.1093\/bioadv\/vbac016","article-title":"decoupleR: Ensemble of computational methods to infer biological activities from omics data","volume":"2","author":"Braunger","year":"2022","journal-title":"Bioinform. Adv."},{"key":"ref_55","doi-asserted-by":"crossref","unstructured":"Langfelder, P., and Horvath, S. (2008). WGCNA: An R package for weighted correlation network analysis. BMC Bioinform., 9.","DOI":"10.1186\/1471-2105-9-559"}],"container-title":["International Journal of Molecular Sciences"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/1422-0067\/27\/16\/7508\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,8,25]],"date-time":"2026-08-25T04:31:03Z","timestamp":1787632263000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/1422-0067\/27\/16\/7508"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,8,21]]},"references-count":55,"journal-issue":{"issue":"16","published-online":{"date-parts":[[2026,8]]}},"alternative-id":["ijms27167508"],"URL":"https:\/\/doi.org\/10.3390\/ijms27167508","relation":{},"ISSN":["1422-0067"],"issn-type":[{"value":"1422-0067","type":"electronic"}],"subject":[],"published":{"date-parts":[[2026,8,21]]}}}