《Frontiers in Cell and Developmental Biology》：Deciphering the cellular landscape and genetic underpinnings of fiber diameter determined by dermal papilla cells in fine-wool sheep
编辑推荐：
羊毛纤维直径（Fiber Diameter, FD）是细毛羊的关键经济性状，其生长过程与毛囊（Hair Follicles, HFs）的发育及周期性生长密切相关。因此，全面了解调控毛囊发育的分子机制具有重要意义。本研究对同一饲养条件下但表现出不同纤维直径的六只
羊毛纤维直径（Fiber Diameter, FD）是细毛羊的关键经济性状，其生长过程与毛囊（Hair Follicles, HFs）的发育及周期性生长密切相关。因此，全面了解调控毛囊发育的分子机制具有重要意义。本研究对同一饲养条件下但表现出不同纤维直径的六只半同胞细毛羊的皮肤组织进行了单细胞RNA测序（scRNA-seq），成功构建了毛囊微环境的转录图谱。研究人员共捕获了59,732个高质量细胞，并利用已知标记基因将其系统注释为14种不同的细胞类型。每种细胞类型均表现出独特的基因表达谱。对三个实验组细胞比例的分析显示，与另外两组相比，超细毛组中真皮乳头细胞（Dermal Papilla Cells, DPCs）、毛囊干细胞（Hair Follicle Stem Cells, HFSCs）和外根鞘细胞（Outer Root Sheath cells, ORS）的比例更高。细胞间通讯分析鉴定了多个参与毛囊调控的配体-受体对，表明DPCs可能在调节毛囊生长发育中发挥重要作用。伪时间轨迹分析重建了表皮谱系的发育动态。此外，研究人员发现DPCs可能存在三种不同的状态以执行其生理功能。为了阐明纤维直径的分子机制，研究人员识别并分析了超细、中细和细毛组中各细胞类型的差异表达基因（Differentially Expressed Genes, DEGs），并对在DPCs中显著上调的CRABP1进行了细胞水平验证。结果显示，过表达该基因能显著上调PCNA和CTNNB1的表达，抑制SFRP2和BMP2的表达，进而促进DPCs的增殖。本研究成功构建了周岁细毛羊皮肤组织的单细胞图谱，揭示了DPCs在毛囊发育中的潜在重要作用，并发现了纤维直径的潜在分子机制。这些发现为羊毛囊发育生物学提供了新见解，并为超细毛羊的分子靶向精准育种提供了宝贵参考。
论文解读
研究背景与意义
羊毛品质主要取决于作为皮肤重要附属器官的毛囊（HFs）。毛囊分子调控机制具有高度的保守性和复杂性。羊毛纤维直径（FD）是决定羊毛柔软度和舒适度等关键经济性状的指标。随着生活水平提高，市场对超细羊毛的需求日益增长，培育超细毛羊成为行业焦点。然而，目前关于细毛羊纤维直径变异的分子遗传机制鲜有报道。尽管已有研究通过全基因组关联分析和转录组学筛选了部分候选基因，但在单细胞分辨率下解析毛囊微环境异质性及其对纤维直径影响的机制仍属空白。因此，研究人员旨在通过单细胞RNA测序（scRNA-seq）技术，系统分析阿尔卑斯细毛羊皮肤组织的细胞异质性，揭示调控毛囊发育进而影响羊毛纤维直径的潜在分子机制，为分子靶向精准育种提供理论依据。该研究发表于《Frontiers in Cell and Developmental Biology》。
关键技术方法
研究人员选取了来自同一父系、饲养条件相同的6只一岁阿尔卑斯母羊，按平均纤维直径分为超细毛组、中细毛组和细毛组，每组2只。从肩胛骨后缘采集皮肤组织样本，利用10x Genomics平台进行单细胞RNA测序（scRNA-seq），构建转录图谱。通过生物信息学分析，包括使用Seurat包进行数据标准化、降维和聚类识别，Monocle进行伪时间轨迹分析，以及CellChat进行细胞间通讯网络分析。体外实验部分，研究人员分离培养了真皮乳头细胞（DPCs），通过构建过表达和干扰载体进行细胞转染，并利用实时荧光定量PCR（RT-qPCR）、EDU增殖实验和CCK-8法检测基因表达对细胞增殖的影响及相关信号通路的变化。
研究结果
3.1 细毛羊皮肤组织细胞类型鉴定
研究人员通过对55,048个高质量细胞的UMAP降维聚类分析，结合经典标记基因，成功将细胞注释为14种不同类型，涵盖上皮细胞、真皮细胞、内皮细胞和免疫细胞等谱系。比较不同纤维直径组的细胞比例发现，超细毛组中真皮乳头细胞（DPCs）、外根鞘细胞（ORS）和毛囊干细胞（HFSCs）的比例显著高于其他两组，而内根鞘细胞（IMCs）和内皮细胞比例较低。免疫组化验证了标记基因的表达模式，证实了细胞类型注释的可靠性。
3.2 免疫组化验证标记基因
通过对皮肤组织结构进行苏木精-伊红（HE）染色观察，确认了表皮、真皮及毛囊的基本结构。免疫组化染色进一步验证了特定标记蛋白的定位：KRT15在表皮细胞中强阳性，VIM在真皮细胞中阳性，KRT35在毛干（HS）中阳性，KRT71在内根鞘（IRS）中阳性，HOXC13在基质（Mx）中阳性，TAGLN在外根鞘（ORS）中强阳性。这些结果与单细胞数据分析高度一致。
3.3 表皮-真皮相互作用共同调控毛囊发育
利用CellChat分析细胞间通讯网络，研究人员发现所有细胞类型之间存在广泛且复杂的配体-受体相互作用，其中上皮细胞占据主导地位。真皮乳头细胞（DPCs）和毛囊干细胞（HFSCs）是活跃的信号发送者。DPCs主要通过SEMA3、RA和IGFBP等信号通路传递信号，并在BMP、FGF、TGF-β、PDGF、DHEA和SLIT等经典毛囊发育通路中充当主要的信号发送源，表明其在调控毛囊发育中具有核心地位。
3.4 表皮分化过程重建
针对毛囊表皮细胞谱系（包括基质Mx、内根鞘IRS和毛干HS），研究人员利用伪时间分析重建了其分化轨迹。结果显示，轨迹起始于Mx细胞，随后分化为两个分支：命运1主要为HS，命运2主要为IRS。这一轨迹清晰地反映了毛囊发育过程中，干细胞先激活分化为基质细胞，再进一步向毛干和内根鞘分化的连续动态过程，验证了预期的分化路径。
3.5 真皮乳头细胞分化过程重建
鉴于DPCs的关键作用，研究人员对其进行了亚群细分，最终鉴定出三个不同的亚群（DP1、DP2、DP3）。KEGG富集分析揭示了功能异质性：DP1主要参与感知微环境和整合促生存信号；DP2主要涉及细胞骨架重塑和细胞间粘附；DP3则富含Hedgehog、Wnt和TGF-β等信号通路，暗示其在调节毛囊发育中起重要作用。伪时间分析表明，DPCs从原始的DP3状态开始，分别向DP1和DP2状态分化，呈现连续的命运决定过程。
3.6 与羊毛纤维直径相关的组间差异表达基因
为探究影响纤维直径的分子机制，研究人员比较了不同组别间的差异表达基因（DEGs）。在超细毛组与细毛组的DPCs比较中，鉴定出147个显著上下调基因。其中，CRABP1在超细组中显著上调，鉴于其在视黄酸转运及小鼠毛发曲率调节中的作用，研究人员推测其可能通过调节毛囊形态影响纤维直径。此外，LGR5、COL1A1、KRT10以及MYOC等基因也被鉴定为与纤维直径潜在相关的关键候选基因。
3.7 CRABP1对真皮乳头细胞增殖的影响
3.7.1 真皮乳头细胞的分离纯化与鉴定
研究人员通过组织块法成功分离培养出DPCs，细胞呈梭形或不规则形态。免疫荧光染色证实细胞高表达经典的DPCs标记蛋白α-SMA和SOX2，验证了分离细胞的纯度与身份，为后续功能实验奠定基础。
3.7.2 过表达CRABP1对真皮乳头细胞增殖的影响
功能实验显示，过表达CRABP1能显著促进DPCs的增殖能力。机制研究表明，过表达CRABP1会上调增殖相关基因PCNA和CTNNB1的表达，同时抑制SFRP2和BMP2的表达。这表明CRABP1通过调控Wnt/β-catenin和BMP信号通路，正向调控DPCs的增殖活性。
3.7.3 敲低CRABP1对真皮乳头细胞增殖的影响
与之相反，敲低CRABP1则显著抑制了DPCs的增殖。在分子水平上，敲低CRABP1导致PCNA和CTNNB1表达下调，而SFRP2和BMP2表达上调。这进一步证实了CRABP1在促进DPCs增殖及维持毛囊生长活性中的关键作用。
讨论与结论总结
在讨论部分，研究人员指出，本研究利用scRNA-seq技术构建了高分辨率的细毛羊皮肤组织单细胞图谱，弥补了该物种在毛囊微环境异质性研究上的空白。研究发现DPCs在超细毛组中比例更高，且通过FGF、BMP等多条信号通路与表皮细胞进行密切的双向通讯，发挥了“信号枢纽”的作用。特别地，研究人员首次在细毛羊中揭示了DPCs存在三种功能状态（DP1、DP2、DP3），其中DP3处于原始状态，具有多向分化潜能。研究还重点探讨了CRABP1的分子机制，证实其通过抑制SFRP2和BMP2，激活PCNA和CTNNB1，从而促进DPCs增殖，最终影响羊毛纤维直径。此外，研究鉴定出的MYOC等其他候选基因为后续研究提供了新方向。尽管样本量有限，但这项初步探索为理解毛囊发育生物学及实现超细毛羊的分子设计育种提供了重要的理论基础和候选基因资源。
结论
综上所述，本研究利用scRNA-seq技术探究了细毛羊皮肤组织内的细胞异质性。研究鉴定了14种不同的细胞类型，描绘了毛囊发育和分化的细胞特征，阐明了细胞间通讯网络，并进一步揭示了调控纤维直径的潜在分子机制。这些发现共同为毛囊生理学和繁殖的生物学研究以及细毛羊育种实践提供了新颖的见解。
