0 引言
【研究意义】巴马香猪原产于广西巴马瑶族自治县,是中国优良的小型地方猪品种之一[1]。既有独特的肉用经济价值[2],又可用作良好的实验动物模型[3-4]。皮肤是动物外皮系统中最大的器官,在保护机体内环境和抵御外界环境病原菌侵入等方面起着重要屏障作用[5]。同时,猪皮是皮革产业的重要原料之一[6],也是现代医学异种皮肤移植的理想材料[7]。测量巴马香猪周身皮肤厚度,了解巴马香猪皮肤厚度的变化规律,可增加人们对巴马香猪皮肤的认知。在巴马香猪群体中验证7号染色体影响皮肤厚度的主效QTL,可为进一步在巴马香猪群体中大规模开展皮肤厚度等形态变化的分子遗传控制机理及其相关基因功能研究奠定基础,为更好地开发利用巴马香猪提供参考。【前人研究进展】巴马香猪作为人类疾病模型已在皮肤创伤移植、心血管疾病和肌肉萎缩等疾病研究中应用。例如：2006年陈俊颖和魏泓用巴马香猪皮制作无细胞真皮,可作为创伤修复的真皮支架[8];2008年刘宇以洛伐他汀为典型药物研究巴马香猪和人代谢的异同[9]。目前,涉及皮肤厚度的研究主要集中在皮肤科学的医药研究和皮肤病理治疗上,而针对皮肤厚度的分子遗传控制机理研究鲜见报道。比如,SCHMOOK 等2001 年报道了皮肤厚度对不同药物的体外渗透力的影响[10];SETO 等2010 年报道了不同方法估计人皮肤厚度的研究[11];WONG等2011 年报道了乳腺癌患者在乳房切除术后放射治疗的皮肤厚度变化情况[12];YONDEDA等2011 年报道了在成年女性中骨密度与皮肤厚度的相关性[13]。在猪中,有少量文献报道了有关猪皮肤厚度的研究,比如猪皮肤厚度与瘦肉率的关系研究[14],云南本地猪及其杂交猪的皮肤厚度变化情况[15],猪皮部位差异的基本规律及对制革工艺的影响[16]等。目前尚没有关于巴马香猪的皮肤厚度研究报道。而从分子遗传学角度研究猪皮肤厚度只有AI等2014年报道的一篇关于猪皮厚性状连锁分析和全基因组关联分析的研究。该研究在二花脸×白色杜洛克F2资源群体和苏太猪群体中,利用194个微卫星标记和猪60K SNP芯片,发现猪4、7和15号染色体存在影响皮肤厚度主效QTL,其中7号染色体~35Mb处的主效QTL效应最强[17]。【本研究切入点】以中国地方小型猪品种巴马香猪为研究素材,随机选取了50头300日龄的成年猪,测量猪只周身9个典型部位的皮肤厚度,统计皮肤厚度在不同部位和不同性别间的差异性以及不同部位间的相关性。在巴马香猪群体中进行了皮肤厚度与7号染色体34.5—36.2 Mb候选区域46个SNPs位点的关联分析,分析不同部位皮肤厚度与这个区域SNPs的关联性,根据强关联的SNPs及其物理位置上所在或临近基因,结合基因的生物学功能,推测影响皮肤厚度变化的位置候选基因。【拟解决的关键问题】全面地了解巴马香猪周身皮肤厚度的变化情况,揭示其变化规律。同时在巴马香猪群体中验证在SSC7影响皮肤厚度变化的主效QTL,探讨影响猪皮肤厚度的遗传控制机理,增强人们对猪皮肤的认知。1 材料与方法
1.1 试验动物
试验猪品种为巴马香猪,2014年购自广西巴马自治县巴马香猪繁育基地,统一饲养于江西南昌蒋巷镇国鸿生态园猪场,供猪只自由采食和饮水。购买时,巴马香猪的日龄为40—50日龄,其中公猪于30—35日龄被阉割。所有猪只饲养至（300±5）日龄,运至附近的国鸿生猪屠宰场,共屠宰了319头巴马香成年猪,其中阉割公猪160头,母猪159头,为43头父本和61头母本的后代。1.2 方法
1.2.1 表型测定 巴马香猪宰前禁水禁食24h后进行屠宰。在剔毛、去除内脏后,随机抽取50头（公猪23头,母猪27头,为17头父本和20头母本的后代）进行皮肤组织取样。用手术刀和剪刀,分别在头脸部、肩部、背部、肷部、臀部、胸部、下腹部、腋下和前管等9个部位（图1）处,取约5 cm2大小的完整皮肤组织样品。然后用手术刀垂直割取皮肤形成平整的切面,立刻用精度为0.01mm的电子游标卡尺垂直卡在皮肤切割面上进行测量,并读取、记录皮肤厚度的数据。1.2.2 DNA提取及SNPs位点判型 采集0.5—1.0 g猪耳组织,按照标准酚/氯仿法提取DNA。其A260/A280比值在1.8—2.0且琼脂糖电泳检测片段长度大于10 kb的,视为合格DNA,储存于-30℃冰箱待用。
使用Sequenom公司的MassARRAY®时间飞行质谱分型技术,进行7号染色体34.5—36.2 Mb内的46个SNPs位点的基因分型。即根据SNPs碱基分子量存在差异这一特性,通过PCR扩增得到一定数量级的包含SNPs在内的目的片段,再利用质谱仪测定不同基因型目的片段的质荷比,对比片段分子量差异,从而实现SNPs位点基因分型。
显示原图|下载原图ZIP|生成PPT图1巴马香猪皮肤取样的9个躯体部位
-->Fig. 1The 9 body sites where skin were sampled in Bamaxiang pig
-->
1.2.3 统计分析 巴马香猪周身9个部位皮肤厚度的描述性统计均使用R语言的相关统计函数进行计算分析。用boxplot绘制箱式分布图,使用shapiro.test进行正态分布检验,使用单因素方差分析aov检验不同部位皮肤厚度的总差异性,使用t.test检验两两不同部位间皮肤厚度的差异显著性,使用cor.test进行相关系数计算和显著性分析。
运用R语言SNPassoc 软件包中的association函数进行不同部位皮肤厚度与46个位点候选基因型的相关性分析。建立的广义混合线性模型如下：y=u+geno+ sex+batch+e。其中y为测定的皮肤厚度表型值,u表示皮肤厚度的平均数,geno表示SNPs的基因型,sex为性别效应,batch为批次效应,e为残差。
2 结果
2.1 皮肤厚度的描述性统计分析
精确测量了50头300日龄成年巴马香猪的头脸、肩、背、肷、臀、管、胸、下腹和腋下等9个典型部位的皮肤厚度（表1,图2）。从厚到薄依次是肷、背、肩、头脸、臀、管、胸、腋下和下腹,其中肷部和背部的皮肤厚度分别为（5.15±0.92）和（4.97±0.85） mm,下腹和腋下的皮肤厚度分别为（1.77±0.36）和（1.97±0.68）mm。背部的皮肤厚度变异度最小,为17.1%,腋下的最大,为34.8%。正态分布检验表明,头脸、肩、肷、胸、背和下腹皮厚的P值均大于0.05,符合正态分布;而腋下、管和臀的P值小于0.05,说明这个3个部位的皮肤厚度不符合正态分布。腋下的可能是由于不同性别间的皮肤厚度存在显著性差异（P = 0.021）;管部皮肤厚度在母猪和阉割公猪间的差异接近显著水平（P = 0.072）,可能是造成管部正态检验显著的原因之一;臀部皮肤厚度可能是由于测量的个体数偏少,测量个体的随机性造成的。
Table 1
表1
表1巴马香猪周身9个部位皮肤厚度描述性统计量a
Table 1Descriptive statistics for skin thickness at 9 body sites in Bamaxiang pig
| 部位 Body site | 肷部 Loin | 背部 Back | 肩部 Shoulder | 头脸 Cheek | 臀部 Rump | 管部 Trotter | 胸部 Chest | 腋下 Armpit | 下腹 Belly |
|---|---|---|---|---|---|---|---|---|---|
| 平均数 Mean | 5.15 | 4.97 | 4.33 | 3.03 | 2.87 | 2.42 | 2.03 | 1.97 | 1.77 |
| 标准差 SD | 0.92 | 0.85 | 0.82 | 0.74 | 0.73 | 0.74 | 0.40 | 0.68 | 0.36 |
| 最大值 Max | 7.06 | 6.43 | 5.84 | 4.80 | 5.46 | 4.84 | 2.91 | 3.96 | 2.54 |
| 最小值 Min | 2.81 | 3.01 | 2.79 | 1.64 | 1.78 | 1.08 | 1.27 | 1.02 | 1.11 |
| 变异系数 CV | 18.1% | 17.1% | 19.2% | 24.5% | 25.6% | 31.0% | 20.0% | 34.8% | 20.7% |
| P值bP value | 0.873 | 0.085 | 0.349 | 0.140 | 0.004 | 0.024 | 0.376 | 0.001 | 0.106 |
新窗口打开
显示原图|下载原图ZIP|生成PPT图2巴马香猪不同躯体部位的皮肤厚度
-->Fig. 2skin thickness at different body sites in Bamaxiang pig
-->
2.2 不同部位皮肤厚度的差异分析
巴马香猪9个部位皮肤厚度在整体上差异极显著,P值为2.95×10-117,表2列出了它们两两之间相互比较的差异显著性。由此可知,除了背部与肷部、头脸部与臀部、胸部与腋下以及下腹与腋下的皮肤厚度差异不显著（P>0.05）以外,其余各个部位之间的皮肤厚度均存在极显著差异（P<0.01）,其中以背部、肷部与下腹、腋下之间的皮肤厚度差异最为显著。性别可能是影响皮肤厚度的影响因素之一,因此本文以性别作为分类依据进行皮肤厚度的显著性检验。检验结果表明,阉割的公猪个体腋下皮肤厚度显著小于母猪的皮肤厚度（P = 0.021）,其它部位均不显著,但管部皮肤和下腹皮肤在母猪和阉割公猪间的厚度差异接近显著水平,P值分别为0.072和0.054（图3）。
显示原图|下载原图ZIP|生成PPT图3巴马香猪不同躯体部位皮肤厚度在不同性别间的差异
F为母猪,M为阉公猪
-->Fig. 3Difference of skin thickness between two gander groups at 9 body sites in Bamaxiang pigs
F presents females, M presents males
-->
Table 2
表2
表2巴马香猪9个部位皮肤厚度差异显著性P值
Table 2P value of t-test for skin thickness at 9 body sites in Bamaxiang pigs
| 肷部 Loin | 背部 Back | 肩部 Shoulder | 头脸 Cheek | 臀部 Rump | 管部 Trotter | 胸部 Chest | 腋下 Armpit | |
|---|---|---|---|---|---|---|---|---|
| 背部 Back | 0.293 | |||||||
| 肩部 Shoulder | 8.95×10-6 | 2.56×10-4 | ||||||
| 头脸 Cheek | 7.65×10-22 | 5.34×10-21 | 1.05×10-12 | |||||
| 臀部 Rump | 7.02×10-24 | 3.12×10-23 | 5.07×10-15 | 0.273 | ||||
| 管部 Trotter | 6.44×10-29 | 1.17×10-28 | 5.38×10-21 | 7.60×10-5 | 0.003 | |||
| 胸部 Chest | 4.09×10-32 | 3.17×10-33 | 1.43×10-27 | 1.78×10-12 | 4.98×10-10 | 0.002 | ||
| 腋下 Armpit | 4.27×10-34 | 1.44×10-34 | 1.14×10-27 | 3.46×10-11 | 6.11×10-9 | 0.002 | 0.595 | |
| 下腹 Belly | 1.36×10-33 | 7.17×10-35 | 6.75×10-30 | 1.13×10-16 | 2.40×10-14 | 5.64×10-7 | 0.001 | 0.076 |
新窗口打开
2.3 不同部位皮肤厚度的相关性分析
巴马香猪9个部位之间的皮肤厚度除了下腹部与背部、肩部和头脸部的不相关外,其余部位两两之间均呈正相关（P<0.05,表3）。肩、背和肷三个部位的皮肤厚度两两之间都呈极显著正相关,管与腋下也呈现极显著正相关,它们之间的相关系数均大于0.6,P值均远小于0.01。其余部位两两之间的皮肤厚度相关性在0.28—0.6之间,P值也都小于0.05,均呈极显著正相关或显著正相关。2.4 皮肤厚度表型与候选SNPs位点的关联分析
根据此前2014年AI等[17]的研究,在7号染色体的目标QTL区域34.5—36.2 Mb内选取和设计了46个SNPs位点,并进行MassARRAY的质谱基因分型,图4显示了这些SNPs在基因组中的物理位置。利用1.2.3方法中线性加性模型进行这46个SNPs位点的基因分型数据与9个部位皮肤厚度表型的关联分析。这9个皮肤厚度表型均与这个区域的某些SNPs存在极显著的相关,与这个区域SNPs位点的相关性从强到弱依次是：肷部>肩部>背部>腋下>臀部>胸部>头脸>管部>下腹（图4）。Table 3
表3
表3巴马香猪9个部位皮肤厚度的相关性
Table 3Correlation of skin thickness between pairwise different body sites in Bamaxiang pigs
| 肷 Loin | 背 Back | 肩 Shoulder | 头脸 Cheek | 臀 Rump | 管 Trotter | 胸 Chest | 腋下 Armpit | 下腹 Belly | |
|---|---|---|---|---|---|---|---|---|---|
| 肷 Loin | 8.94×10-9 | 3.72×10-7 | 5.49×10-5 | 5.13×10-6 | 2.57×10-4 | 4.10×10-5 | 2.09×10-4 | 0.003 | |
| 背 Back | 0.708 | 2.96×10-6 | 7.83×10-4 | 0.002 | 0.020 | 7.33×10-4 | 3.58×10-4 | 0.312 | |
| 肩 Shoulder | 0.648 | 0.607 | 0.021 | 8.06×10-4 | 0.042 | 0.002 | 1.80×10-4 | 0.342 | |
| 头脸 Cheek | 0.538 | 0.460 | 0.325 | 0.008 | 0.002 | 8.13×10-4 | 0.005 | 0.148 | |
| 臀 Rump | 0.595 | 0.420 | 0.459 | 0.370 | 0.005 | 2.89×10-4 | 0.028 | 0.005 | |
| 管 Trotter | 0.495 | 0.328 | 0.289 | 0.433 | 0.394 | 8.92×10-4 | 3.28×10-6 | 7.37×10-6 | |
| 胸 Chest | 0.546 | 0.462 | 0.427 | 0.458 | 0.492 | 0.455 | 7.53×10-4 | 0.006 | |
| 腋下 Armpit | 0.501 | 0.485 | 0.506 | 0.389 | 0.311 | 0.605 | 0.461 | 0.019 | |
| 下腹 Belly | 0.416 | 0.146 | 0.137 | 0.208 | 0.388 | 0.587 | 0.387 | 0.330 |
新窗口打开
肷部的皮肤厚度与SNP（chr7: 34856565）的相关显著性最强（Pcor = 5.15×10-6）,这个SNP同时也是肩部皮肤厚度最相关（Pcor=5.75×10-6）的位点;SNP（chr7: 35543837）是腋下（Pcor = 3.05×10-5）、臀部（Pcor = 0.010）、胸部（Pcor = 0.013）和头脸部（Pcor = 0.025）皮肤厚度的最相关位点,也是肩部皮肤厚度的次最相关位点。SNP（chr7: 35573869）则是背部皮肤厚度的最相关位点（Pcor=1.17×10-5）,SNPs（chr7: 35543837和chr7: 34856565)次之（表4）。这3个与皮肤厚度强关联的SNPs分别位于GRM4,TAF11和ANKS1A等3个基因中（表4）。
显示原图|下载原图ZIP|生成PPT图4巴马香猪9个不同部位皮肤厚度与7号染色体候选区域46个SNPs位点的关联分析
基因用黑框表示;SNPs的所在物理位置用黑短竖线标示;虚线表示P = 0.01的显著水平
-->Fig. 4Association assay between skin thickness at 9 body sites and 46 candidate SNPs on SSC7 in Bamaxiang pig
Black solid box presents gene; black vertical line indicates SNP position; dash line presents the significant level of a 0.01 probability
-->
Table 4
表4
表4关联分析中P值最高SNPs位点的信息
Table 4Information of the most significant SNPs in association assay
| 皮肤厚度性状 Traits of Skin thickness | SNP位置 SNP position | 等位基因 Alleles | 次等位基因 Minor allele | 次等位基因频率 Minor allele frequence | 校正P值a Corrected P value | 附近基因b Nearest gene |
|---|---|---|---|---|---|---|
| 肷 Loin | 34856565 | A/G | G | 0.410 | 5.15×10-6 | GRM4(0) |
| 背 Back | 35573869 | A/G | G | 0.498 | 1.17×10-5 | ANKS1A(0) |
| 肩 Shoulder | 34856565 | A/G | G | 0.410 | 5.75×10-6 | GRM4(0) |
| 头脸 Cheek | 35543837 | C/G | C | 0.473 | 0.025 | TAF11(0) |
| 臀 Rump | 35543837 | C/G | C | 0.473 | 0.010 | TAF11(0) |
| 管c Trotter | 35184667 | T/G | T | 0.444 | 0.088 | SPDEF(-29.7 Kb) |
| 胸 Chest | 35543837 | C/G | C | 0.473 | 0.013 | TAF11(0) |
| 下腹d Belly | 34676565 | C/T | C | 0.388 | 0.088 | MLN(-3.7 Kb) |
| 腋下 Armpit | 35543837 | C/G | C | 0.473 | 3.05×10-5 | TAF11(0) |
新窗口打开
3 讨论
研究猪及其各类器官的形态结构变化、生长发育和病理变化等生命规律可为人类的相关研究提供重要借鉴。猪皮肤在形态结构上跟人类皮肤相近,猪皮可作为人类皮肤替代物进行药物试验或模拟治疗[18],利用猪皮作为人类烧伤皮肤的异种植皮来源也得到深入研究[7]。为了更好地了解猪皮肤的形态构造变化和生长发育规律,本文以深入了解巴马香猪皮肤厚度的变化规律和探究其分子遗传控制机理为出发点,开展相关猪皮肤厚度的研究工作,即对中国地方小型猪种巴马香猪进行了9个典型部位皮肤厚度的表型统计分析,及其与7号染色体影响皮肤厚度主效QTL区域46个SNPs位点的关联分析。有研究表明,在哺乳动物中皮肤厚度的变化因个体或部位而异,比如在人中一般为0.5—4.0 mm[19],在马中为1.5—4.6 mm[20]。针对猪皮肤厚度变化规律的研究,有两篇较为经典的文献报道。1993年陈德端和王瑞华研究了云南当地猪和各杂交组合猪只的肩、背、臀和腹部4个部位皮肤厚度变化,发现猪皮肤厚度在躯体的变化通常为1.1—5.8 mm[15]。1994年王全杰比较全面地介绍了中国国产猪皮的厚度变化规律：整张猪皮以尾根部为底脊背为中线,呈现长把梨状三角形,厚度由底部沿着中垂线向顶角逐级递减,直到其2/3处不再明显下降,基本保持稳定[16]。笔者经过精确测量,发现巴马香猪的皮肤厚度也因个体不同或部位不同而不一样,变化范围为1.02—7.06 mm。巴马香猪周身9个部位皮肤从薄到厚依次为：下腹<腋下<胸部<管部<臀部<头脸<肩部<背部<肷部。由此可见,巴马香猪躯体背侧朝上的部位皮肤较厚,而躯体两端和肢蹄端的皮肤厚度相对背侧较薄,而躯体的腹侧朝下或朝内的部位皮肤最薄。
巴马香猪躯体各部位的皮肤厚度变化的差异非常大,尤其是背侧与腹侧皮肤相比较;不同部位间的皮肤厚度呈现不同程度的正相关,同侧皮肤厚度的相关性强或较强,不同侧皮肤厚度相关性较弱或无相关,如背部皮肤厚度与肷部、肩部的呈现极显著正相关,而背部皮肤与下腹皮肤不呈现显著相关（P>0.05）。性别可能是影响皮肤厚度的一个重要影响因素。DERRAIK等2014年研究发现成年男性大腿上的皮肤显著地厚于女性的皮肤,但在腹部的皮肤男女性的差别不显著[21]。本研究发现在大部分躯体部位的皮肤厚度在母猪和阉割公猪间没有显著差别,但腋下皮肤在两者之间存在显著差异,这也可能是造成腋下的皮肤厚度不呈正态分布的原因之一。
2014年AI等利用二花脸×白色杜洛克的F2代资源群体,使用194个微卫星标记进行了猪皮肤厚度性状的连锁分析,发现了3个主效QTL位点,分别位于7号染色体的58 cM、15号的78 cM 和4号的107 cM 上,其中7号染色体的QTL 解释了23.9%遗传变异。同时,新增苏太猪群体,使用猪60K SNP芯片,进行了两个群体猪皮肤厚度性状的全基因组水平关联分析,将验证了这3个主效QTL位点。将7 号染色体主效QTL置信区域缩短到683kb,将15 号染色体主效QTL 区域也缩短到2.2 Mb [17]。本研究则利用中国南方地方小型猪巴马香猪群体,进行了躯体9个部位皮肤厚度与7号染色体34.5—36.2Mb候选区域46个SNPs位点的关联分析。首先验证了巴马香猪群体也存在影响猪皮肤厚度的7号染色体主效QTL,为今后利用巴马香猪群体进行皮肤厚度的全基因组水平关联分析奠定了基础。其次,发现背侧（如：背部、肷部和肩部）的皮肤厚度与候选SNPs位点的关联性比躯体两端或腹侧（如：管部、胸部、臀部、下腹和头脸部）的更高。这提示以下两点：（1）在今后的皮肤厚度关联分析研究中,肩、背和肷部的皮肤取样部位可能会更为合适,在统计上造成的误差更小;（2）腹侧皮肤与背侧皮肤的皮肤厚度有可能存在某些不同的内在分子遗传控制机理,这一点值得今后更深入的研究和探讨。
在7号染色体34.5—36.2Mb的46个SNPs中,有3个与皮肤厚度关联性最强的SNPs值得进一步关注（表4）,分别是SNP（Chr7: 34856565）,SNP（Chr7: 35543837）和SNP（Chr7: 35573869）,它们分别位于GRM4,TAF11和ANKS1A基因中。GRM4基因编码一种代谢性谷氨酸受体,参与抑制神经系统中的cAMP级联反应（http://www.ncbi.nlm.nih.gov/gene/ 2914）。TAF11基因编码TF11D的小亚基,是一种转录因子,参与RNA聚合酶II的转录起始过程（http://www.ncbi.nlm.nih.gov/gene/6882）。ANKS1A基因编码一种锚蛋白重复及SAM结构域蛋白（http://www.ncbi.nlm.nih.gov/gene/23294）,也称作为奥丁（Odin）。PANDEY等2002年克隆了奥丁,发现奥丁是一种含有酪氨酸磷酸化结合域的蛋白,揭示出奥丁可能在生长因子受体信号通路起着重要的负调控作用[22]。2013年TANG等进一步证实奥丁在表皮生长因子受体的循环利用上行使着效应物功能,通过贩运作用回收或降解表皮生长因子受体的方式来稳定表皮生长因子受体在细胞中的剂量[23]。从生物学功能上看,GRM4和TAF11这两个基因似乎与皮肤生长发育并不相关,而ANKS1A基因编码的蛋白与表皮生长因子受体有非常密切的相互作用,值得深入研究。7号染色体34.5—36.2Mb区域是一个基因富含的区间,包含有ITPR3,UQCC2,LEMD2,MLN,GRM4,HMGA1,RPS10,SPDEF,PACSIN1,C6orf106,SNRPC,UHRF1BP1,TAF11,ANKS1A,TCP11,SCUBE3,TEAD3,RPL10A,FANCE和PPARD等20个基因（图4）。这些基因均可能是影响皮肤厚度的因果基因,需要进一步做相关的生物学功能实验来验证。其中,HMGA1基因也值得关注。HMGA1基因编码一种参与多个细胞过程的非组蛋白转录因子,常被乙酰化存在于细胞核内,优先结合在双链DNA的AT富集区的小沟上（http://www.ncbi.nlm.nih.gov/ gene/3159）。有文献报道,在哺乳动物中HMGA1基因的突变影响着身高和体型的变化[24-25]。据此可初步推断,ANKS1A和HMGA1基因为影响巴马香猪皮肤厚度变化的候选因果基因,但均需进行生物学功能试验的深入研究,加以证实。
4 结论
首次较全面地测量了中国小型地方猪品种巴马香猪周身皮肤厚度,揭示了巴马香猪皮肤厚度在不同部位之间的变化规律。在巴马香猪群体中验证了影响皮肤厚度的7号染色体主效QTL位点,肷部、背部和肩部皮肤厚度表型性状与候选SNPs位点关联性更强,可能适合下一步大规模深入分析。ANKS1A和HMGA1基因可能是影响皮肤厚度的候选因果基因,但需要进一步生物学功能试验证明。The authors have declared that no competing interests exist.
