基于大家系设计的静息心率与常见慢性病双表型遗传度估计

郑鸿尘，薛恩慈，王雪珩，陈 曦，王斯悦，黄 辉，江 锦，叶 莺，黄春兰，周 筠，高文静，余灿清，吕 筠，吴小玲，黄小明，曹卫华，严延生，吴 涛，李立明

(1. 北京大学公共卫生学院流行病与卫生统计学系，北京 100191； 2. 福建省疾病预防控制中心地方病防治科，福州 350001； 3.福建省漳州市南靖县疾病预防控制中心卫生科，福建南靖 363600； 4. 首都医科大学附属天坛医院国家神经系统疾病临床医学研究中心，北京 100070)

静息心率是指个体在安静、无活动状态下心脏每分钟搏动的次数。静息心率不仅是心脏生理功能的重要指标，既往研究发现静息心率也是多种慢性非传染性疾病(如高血压、糖尿病、冠心病)的独立预测因素[1-3]。既往双生子研究与家系研究提示，静息心率与血压、血脂、血糖、体重指数(body mass index, BMI)等多种表型的关联基础可能受到相同的遗传因素作用[4-5]。一项欧洲人群的全基因组关联研究发现，17个与静息心率相关的单核苷酸多态性位点与心血管疾病、代谢性疾病和肾疾病有关[6]。以上研究结果提示，静息心率与多种慢性病可能存在共同的遗传基础。

遗传度是评价遗传因素对性状变异影响的指标，能够反映总体遗传效应决定性状的比例，双生子设计和家系设计则是估计遗传度的常用方法。由于家系设计基于人群抽样，其得到的遗传度估计往往更为保守。大家系设计则能够得到较为丰富的家系关系，使得遗传度估计更加稳健[7]。既往双生子研究发现，静息心率的遗传度为0.61[4]，而家系设计得到的静息心率遗传度为0.11～0.31[8]。目前，我国人群中静息心率单表型、双表型遗传度相关研究较少见。因此，本研究基于福建省南靖县募集的大家系，利用基线调查横断面数据，探索了静息心率与常见慢性病的遗传相关性。

1 资料与方法

1.1 研究对象

本研究资料来源于“中国环境流行病学特殊人群队列研究”项目[9]，研究现场为福建省漳州市南靖县土楼地区，包括2015年8月至2017年12月塔下村、曲江村、南欧村募集的张姓居民，草坂村、图美村、背岭村募集的陈姓居民。研究对象的纳入标准为：(1)所有张姓、陈姓氏族成员及其有血缘关系的亲属和配偶；(2)自愿参与并完成问卷调查、家系信息收集、生化检查、体格检查；(3)无重大躯体性、精神性疾病(如残疾、重度痴呆等)， 可以配合调查。

研究开始前获得北京大学生物医学伦理委员会审查批准(IRB00001052-14021)， 所有研究对象包括患者和健康人均签署了知情同意书。

1.2 家系构建

家系构建的依据为研究对象提供的家系信息，搜寻其亲属中被调查的研究对象，同时依据南靖县张姓、陈姓族谱进行家系信息补充并排除仅有1名调查对象的家系。

1.3 结局定义

基线调查包括问卷调查、体格检查和血生化检测。采用电子血压计测量血压和静息心率，采用皮尺测量身高、腰围、臀围，以上变量每名调查对象测量2次，若两次测量结果之差大于5个单位，则测量第3次并取后两次测量值的均值进入统计分析。

研究涉及结局的定义：(1)高血压：自报二级以上医院诊断高血压病史或两周内服用降血压药物或体检中收缩压≥140 mmHg和(或)舒张压≥90 mmHg[10]；(2)糖尿病：自报二级以上医院诊断的糖尿病史或空腹血糖≥7.0 mmol/L[11]；(3)血脂异常：自报二级以上医院诊断的高脂血症史或血生化检测总胆固醇≥6.2 mmol/L或甘油三酯≥2.3 mmol/L或高密度脂蛋白<1.0 mmol/L或低密度脂蛋白≥4.1 mmol/L[12]。

1.4 统计学分析

对连续性变量采用Kolmogorov-Smirnov检验其正态性，若符合正态分布，以均数±标准差描述集中趋势和离散趋势。对不符合正态分布的变量，经逆正态变换后纳入分析。对分类变量采用频率描述变量的水平，描述性分析采用Stata 软件13.1版本。连续型变量的单表型、双表型遗传度分析采用最大似然估计的方差组分模型进行，分类变量的单表型、双表型遗传度分析采用易患性-阈值模型进行[13-14]。为避免年龄、性别及其交互作用可能对结果造成的偏倚，根据既往研究经验，所有分析中均调整性别、年龄、性别×年龄、年龄2、性别×年龄2作为协变量，遗传度估计会在扣除协变量效应之后进行。单表型遗传度定义为加性遗传组分占表型变异的比例，双表型遗传度定义为两表型加性遗传组分相关系数。数据分析使用SOLAR软件(版本 8.1.1)，P<0.05认为差异有统计学意义。

2 结果

2.1 基本人口学特征

本研究共纳入研究对象1 563人(表1)，平均年龄57.2岁，男性占39.4%，全部为汉族。研究对象的高血压患病率为 49.2%， 糖尿病患病率为10.0%，血脂异常的患病率为45.2%。男性各种慢性病患病率均高于女性。本研究构建的家系包括1个七代家系，2个五代家系，23个四代家系，186个三代家系，240个两代家系(表2)。

表1 研究对象基本人口学特征

BMI, body mass index; SBP, systolic blood pressure; DBP, diastolic blood pressure.

表2 家系构建基本信息

2.2 单变量遗传度估计

基于本研究构建的大家系，估计得到静息心率遗传度为0.263(P<0.001)，且协变量对其贡献较小。高血压遗传度为0.404(P=0.001)，血脂异常遗传度为0.799(P<0.001)，其余表型遗传度见表3。

表3 静息心率及常见慢性病单表型遗传度

a, covariates include gender, age, age×gender, age2, gender×age2; b, inverse normal transformation was applied for the variables not following normal distribution. SBP, systolic blood pressure; DBP, diastolic blood pressure; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol;CI, confidence interval.

2.3 双变量遗传度估计

基于本研究构建的大家系，双表型遗传度分析显示，静息心率与高血压、糖尿病、舒张压、血糖、甘油三酯均存在表型相关性，其中静息心率与血糖的表型相关性最高，为0.203(P<0.05)。静息心率与血糖和糖尿病具有较强的遗传相关性，分别为0.485(P<0.05)和0.795(P<0.05)。

3 讨论

慢性非传染性疾病是造成我国居民健康寿命年损失的主要因素。2002年至2017年，中国高血压患病率增长约2.5倍[15-16]，达到44.7%，糖尿病的患病率也持续增长达到了 9.1%[17]。为深入探索常见慢性病的病因，本研究基于福建南靖县土楼地区大家系，分析发现静息心率与高血压、糖尿病、血脂相关指标均存在表型相关性，其中糖尿病患病状态和血糖值与静息心率存在较强的遗传相关性，提示血糖水平和静息心率可能受到共同的遗传因素影响，这为后续遗传致病因素的深入探索提供了重要依据，静息心率相关的遗传危险因素可能同时影响糖尿病的发病风险。

本研究单表型遗传度分析显示，静息心率的遗传度为0.263，接近既往家系研究的结果[8]，但低于双生子研究得出的遗传度[4]。结果提示静息心率为可遗传性状，受到加性遗传变异的影响，但遗传因素对静息心率的决定程度较低。本研究结果显示高血压、收缩压、舒张压的遗传度分别为0.404、0.379和0.306，其中高血压和收缩压除受遗传因素决定外，受年龄性别影响也较大。既往双生子研究和家系研究显示，高血压遗传度为0.3～0.5[18]，收缩压遗传度为0.34[19]，舒张压遗传度为0.368[20]，与本研究结果相近。本研究发现，血糖水平受到遗传因素的影响，而本分析中糖尿病的患者及家系样本有限，通过采用易患性-阈值模型，未发现糖尿病具有统计学意义的遗传度，需进一步扩大家系样本后验证和评估糖尿病的遗传度。本研究还发现血脂异常、总胆固醇、甘油三酯、高密度脂蛋白、低密度脂蛋白受遗传因素影响较大，其遗传度分别为0.799、0.451、0.568、0.852、0.387，与既往遗传度研究结果一致[5]。本研究发现，血脂异常作为综合多项血脂水平的指标具有较高的遗传度，但具体血脂相关指标受到遗传因素影响的程度不一，推测影响不同血脂指标的遗传危险因素可能存在一定差异。

表4 静息心率及常见慢性病双表型遗传度

a, testing the difference between the environmental correlation and 0; b, testing the difference between the genetic correlation and 0； c, testing the difference between the genetic correlation and 1 or -1; d, testing the difference between the phenotypic correlation and 0; e, inverse normal transformation was applied for the variables not following normal distribution. All variables were adjusted for gender, age, age×gender, age2, gender×age2. SBP, systolic blood pressure; DBP, diastolic blood pressure; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; CI, confidence interval.

双表型遗传度是共同遗传因素导致的相关性占两个表型总相关性的比例。双表型遗传度的估计是将两个性状的表型相关系数分解为遗传相关系数和环境相关系数，其中，表型相关系数即表型相关性，表示性状的总体相关性；遗传相关系数即遗传相关性，表示影响两个表型的遗传物质的相似程度；环境相关系数即环境相关性，表示影响两个表型的环境因素的相似程度。既往双表型遗传度研究显示，静息心率与舒张压和人体压力感受器敏感性具有相关性，但主要受到共同的环境因素影响[14]。本研究发现，静息心率与高血压、糖尿病、血脂水平存在表型相关性，其中静息心率与血糖水平相关性最高，并且，静息心率与糖尿病和血糖水平具有较强的遗传相关性，提示其可能受共同的遗传物质控制。既往多个国家和地区的前瞻性队列研究显示，高静息心率显著增加糖尿病的发生风险，即静息心率与糖尿病存在关联[21-22]。而全基因组关联研究和转录组学研究发现，静息心率与糖尿病存在共同的遗传位点，由多基因评分得到其遗传相关性为0.22[23-24]，该结果小于本研究得到的遗传相关性，提示可能尚存在未被发现的同时影响静息心率和糖尿病的遗传危险因素。而静息心率与糖尿病相关的机制可能是：高静息心率时交感神经活性增加,胰岛素抵抗，从而增加2型糖尿病的发病风险[25]。同时，2型糖尿病会引发一系列代谢系统紊乱，如肥胖、血压变化等，导致静息心率升高[26]。

福建土楼家系队列研究是我国仅有的几项家系队列研究之一,由于该地区居民以氏族聚集性生活，极适宜于收集家系资料和开展遗传流行病学研究。本研究利用该大家系资料进行遗传度估计，亲属关系丰富，因此遗传度估计结果更为可靠。同时本研究纳入1 563名研究对象，是我国目前规模较大的多代大家系研究，使得本研究得到的遗传度估计较为稳健。

本研究同时存在一定的局限性，本研究已构建的家系包含2 818名研究对象，后续工作尚需纳入未参与研究的家系对象，扩充得到更完整的家系信息。

综上所述，本研究提示静息心率与常见慢性疾病存在相关性，并且静息心率与糖尿病具有遗传相关性，提示二者可能存在共同的遗传基础，后续可继续深入探索影响静息心率和糖尿病的遗传危险因素。