UGT1A1基因多态性与新生儿黄疸遗传关联性的Meta分析

杨 琳 丁俊杰 周文浩

尿苷二磷酸葡萄糖醛酸转移酶1A1(UGT1A1)是胆红素在肝脏结合的关键酶，可使胆红素与葡萄糖醛酸结合形成结合胆红素,UGT1A1缺陷可导致非溶血性高未结合胆红素血症[1,2]。UGT1A1基因在启动子及编码区的不同多态性改变，可引起UGT1A1基因的表达减低或其酶的活性降低。其中在高加索人中最常见的UGT1A1基因的多态性是启动子TATA盒正常序列A(TA)6TAA中插入一个额外的TA，称为UGT1A1*28 ，即A (TA)6TAA to A (TA)7TAA[3]。在黄种人，编码区的错义突变，尤其是第一外显子UGT1A1*6 G211A(GLY71ARG)最为常见[4]。UGT1A1基因的多态性亦可见于2、4和5外显子等，但由于相关文献较少，尚不能进行Meta分析。国内外已有多项研究报道了UGT1A1基因 GLY71ARG多态性和TATA重复多态性与新生儿黄疸的关联性，但单项研究的结果并不一致，同时有研究显示，黄种人新生儿黄疸水平显著高于高加索人种，黄种人可能存在新生儿黄疸的遗传易感性[5]。为明确不同人群UGT1A1基因 GLY71ARG多态性和TATA重复多态性与新生儿黄疸的关联性，本研究检索UGT1A1基因GLY71AGR和TATA重复多态性与新生儿黄疸遗传关联性的相关文献，采用Meta分析方法进行定量综合，以期为新生儿黄疸的遗传病因学提供有力的证据。

1 方法

1.1 文献纳入和排除标准 同时符合以下条件的文献被纳入：①新生儿病理性黄疸的诊断符合目前国内外公认的定义[6～8]； ②研究设计为病例对照研究，病例组为新生儿黄疸；③文献中能提取UGT1A1基因启动子TATA和GLY71ARG基因型和(或)等位基因的频数；④纳入文献研究对象的种族不限;⑤仅纳入中英文文献。

1.2 文献检索策略

1.2.1 数据库 西文数据库：PubMed、EMBASE、Ovid、Springer和Web of science；中文数据库：中国期刊全文数据库、维普中文科技期刊数据库、万方数据库和中国生物医学文献数据库。文献检索时间均从建库至2010年2月。未手工检索灰色文献。

1.2.2 检索策略 英文检索式：(jaundice or hyperbilirubinemia)and(uridine diphosphate-glucuronosyl transferase 1A1 orUGT1A1 or UDP- glucuronosyl transferase)and(polymorphism)；中文检索式：(黄疸 or 高胆红素血症)and(尿苷二磷酸葡萄糖醛酸转移酶 orUGT1A1)and(新生儿 or 足月儿 or 婴儿)。检索语种不限。

1.3 纳入和排除文献 由杨琳按照纳入和排除标准选择文献，如遇不确定的文献与丁俊杰讨论决定。

1.4 资料提取 按照预先设定的资料提取表，由杨琳提取和录入数据。

1.4.1 研究背景和设计信息提取 包括作者、发表时间、国家、胆红素水平和种族等，以及1.5项下文献质量评价的相关信息。

1.4.2 遗传关联性信息的提取 病例组和对照组的样本量、GLY71ARG A/A、A/G和G/G基因型的频数，等位基因A和G的频数；TATA 7/7、7/6和6/6基因型的频数，等位基因7和6的频数。

1.5 文献质量评价标准 依据NHI-NHGRI研究工作组2007年制定的遗传关联性研究报告规范为基础[9]，并依据文献[10] 选取其中14条标准用于遗传关联性研究的质量评价，包括：①计算把握度：对于前瞻性及回顾性研究行把握度计算，可评估Ⅱ类错误的概率；②对照组特征：对照组和病例组是否来自同一人群，并且充分描述了具备可重复性的特征；③病例组特征：人口学及临床信息是否以文字或表格的形式描述，仅提供表型的定义标准是不充分的；④连锁不平衡检测：用于多态性的选择和(或)结果的解释；⑤多态性的鉴定：根据NCBI提供的序列准确鉴定多态性，仅提供扩增引物序列和(或)仅以之前的文献为参考依据是不充分的；⑥计算错误率：是否以相同或其他方法重复了基因型检测，以计算错误率；⑦Hardy-Weinberg平衡检测：对照组是否进行了Hardy-Weinberg平衡检测或提供的数据可进行Hardy-Weinberg平衡检测；⑧盲法：对实验者和结果评价人员是否实施盲法；⑨多重检验校正：研究位点较多时，存在Ⅰ类错误的概率被扩大，是否调整α水平来控制；⑩协变量调整：在协方差分析中是否进行协变量调整；危险因素分析：是否采用回归模型对于其他相关危险因素进行回归分析，得到不同危险因素的OR值及其95%CI；人群分层分析：根据独立的多态性是否进行人群分层的评估和校正；重复性：是否在另外的独立样本中进行验证；功能验证性实验：仅提供功能实验的参考文献是不充分的。每条标准以“是”(满足此条标准)、“否”(不满足)和“不清楚”(文献未描述)评价。对进入分析的每篇文献逐条进行评价，杨琳独立评价文章质量，如遇不确定的文献与丁俊杰讨论决定。

1.6 统计学方法 采用RevMan5.0软件进行Meta分析。异质性分析采用χ2检验，P≥0.1为研究间具同质性，采用固定效应模型描述；P<0.1为研究间具异质性，采用随机效应模型表达。两组基因型和等位基因频率以OR及其95%CI表示。P<0.05为差异有统计学意义。采用漏斗图检验是否存在发表偏倚。视进入分析文献的人群分布考虑是否采用亚组分析，并根据异质性检验结果行原因分析。

2 结果

2.1 纳入文献的一般情况 共检索到相关文献284篇，其中中文文献18篇，英文文献259篇，其他语种文献7篇(法文文献3篇，韩文文献2篇，西班牙文文献2篇)。通过阅读摘要，有91篇文献经过初筛纳入，进一步阅读全文排除不满足纳入标准的文献60篇。最终22篇文献[11～32]仅入Meta分析(图1)，其中英文文献18篇，中文文献4篇。

图1 文献筛选流程图

Fig 1 Flow chart of article screening and selection process

纳入文献中病例组1 444例，对照组1 835例，纳入文献来自7个国家(中国[11～14，19，26，29]、日本[25,30～32]、泰国[17]、马来西亚[15,20,22,24,27]、印度[16]、土耳其[21,23,28]和美国[18])。文献[12,15,18,20,24,27,28,30～32]描述了纳入研究人群的种族，文献[18]详细说明了人种组成，白种人64.1%，黑种人15.7%，西班牙和葡萄牙人10.5%，亚洲人6.5%，其他人种3.3%，但数据结果未按照种族分类。16篇文献[11～13,15,16, 20～24,26～31]根据总胆红素的峰值>220～342 μmol·L-1、文献[14,25,32]根据日龄总胆红素水平、文献[17,18]根据总胆红素高于P95、文献[19]根据黄疸持续>28 d，总胆红素>100 μmol·L-1诊断新生儿黄疸。13篇文献[11～13,16,18～20,22,24,26～28,31]纳入新生儿黄疸为间接胆红素升高，9篇文献[14,15,17,21,23,25,27,30,32]未说明，其中8篇文献[14,15,17,21,23,25,30,32]除外引起直接胆红素升高的病理性因素；除文献[18,26,27]外，余19篇文献纳入对象均为不明原因黄疸。7篇文献[12,13,15,21,25,26,29]以生理性黄疸新生儿为对照组，13篇文献[14,16～20,22～24,27,28,30,31]以正常健康或无黄疸新生儿为对照组，文献[11]以正常或因肺炎、脐炎入院新生儿为对照组，文献[32]以健康成人为对照组。所有文献均采用PCR扩增目的片段，然后再经直接测序或限制性片段长度多态性检测(RFLP)或单链构像多态性检测(SSCP)、等位基因特异性核苷酸杂交法(ASO)或实时荧光定量PCR法进行UGT1A1基因位点多态性的检测。15篇文献[11～14,16,17,22,23,25,26,28～32]明确说明病例组与对照组基线具可比性，文献[19]病例组与对照组性别差异有统计学意义，文献[20]两组孕周差异有统计学意义，文献[15,18,21,24,27]未明确交代两组的基线可比情况。纳入文献的一般情况见表1。

16篇文献[12～14,17～20,22,24～27,29～32]报道了UGT1A1基因GLY71ARG多态性与新生儿黄疸的关联性，基因型和等位基因频率信息见表2。14篇文献[11,15,16,18,19,21,23～26,28,29,31,32]报道了UGT1A1基因TATA重复多态性，基因型和等位基因频率信息见表3。

2.2 文献质量评价 4篇文献[16～18,20]根据等位基因的出现频率及把握度进行了样本量的计算；文献[27,32]不能认为对照组与病例组来自同一人群；文献[27]对照组特征描述不充分；12篇文献[12～15,18,20,24,27,28,30～32]描述了病例组国家、地区、人种学信息及临床资料；所有文献均未进行连锁不平衡分析；13篇文献[12～17,20,22,24,26,27,29,31]给出了多态性的鉴定序列或测序结果图，可以说明多态性的鉴定结果；6篇文献[13,14,20,22,26，27]用相同或不同方法进行了基因型结果的检测；文献[19]所提供的数据不能进行Hardy-Weinberg平衡检测，余文献对照组均符合Hardy-Weinberg平衡；所有文献均未描述是否采用盲法；纳入文献均未采用多个位点进行比较，不需行多重检验校正；仅文献[12,26]进行了协变量调整；7篇文献[14,16,17,19,20,24,26]进行了危险因素分析；所有文献均未在数据分析时进行人群分层；所有文献均未在其他人群中进行结果重复；文献[15]进行了功能检测，其余文献未提及。文献质量评价结果见图2。

表2 外显子1 GLY71ARG多态性基因型频率和等位基因频率

表3 启动子TATA重复多态性基因型频率和等位基因频率

Notes:NA: not available

图2 纳入22篇文献的质量评价结果

Fig 2 Quality of 22 included studies

Notes 1:power;2:controls characterization;3:case characterization;4: LD exploration;5:polymorphism identification;6:genotyping error check;7:Hardy-Weinberg equilibrium;8: blinding;9:multiple testing ;10:covariate adjusment;11:risks;12:population stratification adjustment;13:replication;14:functional study

2.3 Meta分析结果

2.3.1 人群分层分析的考虑 鉴于UGT1A1基因多态性具有人种差异，为考察不同人群GLY71ARG和TATA多态性与新生儿黄疸的关联性，根据纳入本文分析的22篇文献的人群构成情况，将中国和日本分别作为一个亚组，泰国和马来西亚同处于南亚地区，合并后作为一个亚组；按照人种学分类，印度和土耳其人群属于高加索人群，文献[18]中以高加索人群为主(64.1%)，故进行合并作为高加索人群进行亚组分析。

2.3.2UGT1A1基因GLY71ARG多态性与新生儿黄疸的关联性 16篇文献[4,12～14,17～20,22,24～27,29～31]报道了GLY71ARG多态性与新生儿黄疸的关联性，其中病例组956例，对照组1 553例。进行Egger′s回归分析，并绘制漏斗图，未见发表偏倚(图3)。

图3 GLY71ARG多态性发表偏倚的漏斗图

Fig 3 Funnel plot for publication bias of GLY71ARG polymorphism

2.3.2.1 基因型A/A+G/A与新生儿黄疸的关联性 异质性检验P=0.83，研究间具有同质性，采用固定效应模型进行分析。Meta分析结果显示，OR=2.79，95%CI：2.26～3.43，病例组与对照组基因型A/A+G/A频率差异有统计学意义，P<0.000 01。进一步按人种进行亚组分析，6篇文献[12～14,19,26,29]报道了中国人群的研究结果，Meta分析结果显示，OR=2.84，95%CI：2.14～3.76，差异有统计学意义，P<0.000 01。4篇文献[25,30～32]报道了日本人群的研究结果，Meta分析结果显示，OR=3.22，95%CI：2.03～5.11，差异有统计学意义，P<0.000 01。5篇文献[17,20,22,24,27]报道了马来西亚和泰国人群的研究结果，Meta分析结果显示，OR=2.41，95%CI：1.56～3.72，差异有统计学意义，P<0.000 1。文献[18]报道了高加索人群的研究结果，Meta分析结果显示，OR=1.98，95%CI：0.49～8.03，差异无统计学意义，P=0.34。提示A/A+G/A基因型与中国、日本、马来西亚和泰国人群新生儿黄疸具有关联性，与高加索人群新生儿黄疸无显著关联性(图4)。

图4 GLY71ARG基因型A/A/+G/A与新生儿黄疸关联性的Meta分析

Fig 4 Meta analysis of the association between A/A+G/A genotype of GLY71ARG and hyperbilirubinemia in neonates

2.3.2.2 基因型A/A与新生儿黄疸的关联性 异质性检验P=0.71，研究间具有同质性，采用固定效应模型进行分析。Meta分析结果显示，OR=7.05，95%CI：4.19～11.85，病例组与对照组基因型A/A频率差异有统计学意义，P<0.000 01。进一步按人种进行亚组分析，6篇文献[12～14,19,26,29]报道了中国人群的研究结果，Meta分析结果显示，OR=6.47，95%CI：3.24～12.94，差异有统计学意义，P<0.000 01。4篇文献[25,30～32]报道了日本人群的研究结果，Meta分析结果显示，OR=3.08，95%CI：1.00～9.49，差异无统计学意义,P=0.05。5篇文献[17,20,22,24,27]报道了马来西亚和泰国人群的研究结果，Meta分析结果显示，OR=21.01，95%CI：5.21～84.79，差异有统计学意义，P<0.000 1。文献[18]报道了高加索人群的研究结果，Meta分析结果显示，OR=5.89，95%CI：0.24～145.49，差异无统计学意义，P=0.28。提示基因型A/A与中国、马来西亚和泰国人群新生儿黄疸有显著关联性，与日本和高加索人群新生儿黄疸无显著关联性(图5)。

图5 GLY71ARG基因型A/A与新生儿黄疸关联性的Meta分析

Fig 5 Meta analysis of the association between A/A genotype of GLY71ARG and hyperbilirubinemia in neonates

2.3.2.3 A等位基因与新生儿黄疸的关联性 异质性检验P=0.90，研究间具有同质性，采用固定模型进行分析。Meta分析结果显示，OR=2.78，95%CI：2.33～3.31，病例组和对照组A等位基因频率差异有统计学意义，P<0.000 01。按人种进行亚组分析，6篇文献[12～14,19,26,29]报道中国人群的研究结果，Meta分析结果显示，OR=2.82，95%CI：2.22～3.58，差异有统计学意义，P<0.000 01。4篇文献[4,25,30,31]报道日本人群的研究结果，Meta分析结果显示，OR=2.50，95%CI：1.72～3.62，差异有统计学意义，P<0.000 01。5篇文献[17,20,22,24,27]报道马来西亚和泰国人群的研究结果，Meta分析结果显示，OR=3.01，95%CI：2.07～4.37，差异有统计学意义，P<0.000 01。文献[18]报道了高加索人群的研究结果，OR=2.47，95%CI：0.66～9.25，差异无统计学意义，P=0.18。提示A等位基因与中国、日本、马来西亚和泰国人群新生儿黄疸有显著关联性，与高加索人群新生儿黄疸无显著关联性(图6)。

2.3.3UGT1A1基因TATA重复多态性与新生儿黄疸的关联性 对TATA重复多态性与新生儿黄疸遗传关联性的文献行Egger′s回归分析，并绘制漏斗图，提示存在发表偏倚(图7)。

图6 GLY71ARG A等位基因与新生儿黄疸关联性的Meta分析

图7 TATA重复多态性发表偏倚的漏斗图

Fig 7 Funnel plot for publication bias of TATA

2.3.3.1 基因型7/7+6/7与新生儿黄疸的关联性 异质性检验P=0.002，研究间具异质性，采用随机模型进行分析。结果显示，OR=0.83，95%CI：0.55～1.24，病例组与对照组基因型7/7+6/7频率差异无统计学意义，P=0.36，提示TATA 7/7+6/7基因型与新生儿黄疸无显著关联性。按人种行亚组分析，4篇文献[11,19,26,29]报道中国人群的研究结果，研究间具同质性，OR=0.59，95%CI：0.36～0.96，对照组显著高于病例组，P=0.03。3篇文献[25,31,32]报道了日本人群的研究结果，研究间具同质性，OR=0.15，95%CI：0.04～0.51，对照组显著高于病例组，P=0.003。2篇文献[15,24]报道了马来西亚人群的研究结果，研究间具异质性，OR=1.31，95%CI：0.59～2.92，差异无统计学意义，P=0.51。提示TATA基因型7/7+6/7与中国、日本和马来西亚人群新生儿黄疸无显著关联性(图8)。

5篇文献[16,18,21,23,28]报道了高加索人群的研究结果，研究间具异质性，Meta分析结果显示，OR=1.18，95%CI：0.68～2.02，差异无统计学意义，P=0.56(图8)。进一步分析异质性产生的原因，剔除文献[16](印度人群)，研究间具同质性(P=0.78)，OR=0.97，95%CI：0.72～1.32，差异无统计学意义，P=0.86；提示TATA基因型7/7+6/7与非印度高加索人群新生儿黄疸无显著关联性；与印度人群新生儿黄疸有显著关联性，OR=3.81，95%CI：1.75～8.33，P=0.000 8(图9)。

图8 TATA基因型 7/7+6/7与新生儿黄疸关联性的Meta分析

图9 高加索人群TATA基因型 7/7+6/7与新生儿黄疸关联性的敏感性分析

Fig9 Sensitivity analysis of the association between 7/7+6/7 genotype of TATA box and hyperbilirubinemia in Caucasian

2.3.3.2 基因型7/7与新生儿黄疸的关联性 异质性检验P=0.80，研究间具有同质性，采用固定效应模型进行分析。Meta分析结果显示，OR=1.42，95%CI：0.92～2.19，病例组与对照组基因型7/7频率差异无统计学意义，P=0.11，提示TATA基因型7/7与新生儿黄疸无显著关联性。进一步按人种行亚组分析，3篇文献[11,26,29]报道了中国人群的研究结果，Meta分析结果显示，OR=1.78，95%CI：0.11～28.69，差异无统计学意义，P=0.68。3篇文献[25,31,32]报道了日本人群的研究结果，Meta分析结果显示，OR=0.38，95%CI：0.04～3.56，差异无统计学意义，P=0.40。2篇文献[15,24]报道了马来西亚人群的研究结果，Meta分析结果显示，OR=2.46，95%CI：0.46～13.06，差异无统计学意义，P=0.29。5篇文献[16,18,21,23,28]报道了高加索人群的研究结果，Meta分析结果显示，OR=1.45，95%CI：0.91～2.33，差异无统计学意义，P=0.12(图10)。提示，中国、日本、马来西亚和高加索人群TATA基因型7/7与新生儿黄疸均无显著关联性。

2.3.3.3 7等位基因与新生儿黄疸的关联性 见图11，异质性检验P=0.002，研究间具异质性，采用随机效应模型进行分析。Meta分析结果显示，OR=0.90，95%CI：0.63～1.29，病例组与对照组7等位基因频率差异无统计学意义，P=0.57。提示TATA 7等位基因与新生儿黄疸无显著关联性。进一步按人种行亚组分析。3篇文献[11,26,29]报道中国人群的研究结果，研究间具同质性，Meta分析结果显示，OR=0.65，95%CI：0.35～1.21，差异无统计学意义，P=0.18。3篇文献[25,31,32]报道了日本人群的研究结果，研究间具同质性，Meta分析结果显示，OR=0.15，95%CI：0.04～0.50，对照组显著高于病例组，P=0.002。2篇文献[15,24]报道了马来西亚人群的研究结果，Meta分析结果显示，OR=1.40，95%CI：0.59～3.29，差异无统计学意义，P=0.45(图11)。提示TATA 7等位基因与中国、日本和马来西亚人群新生儿黄疸无显著关联性。

图10 TATA基因型7/7与新生儿黄疸关联性的Meta分析

图11 TATA 7等位基因与新生儿黄疸关联性的Meta分析

Fig11 Meta analysis of the association between 7 allele of TATA box and hyperbilirubinemia in neonates

5篇文献[16,18,21,23,28]报道了高加索人群的研究，研究间具异质性，Meta分析结果显示，OR=1.17，95%CI：0.80～1.69，差异无统计学意义，P=0.42(图11)。分析异质性产生的原因，剔除文献[16](印度人群)，研究间具同质性(P=0.53)，结果显示，OR=1.07,95%CI:0.84～1.35,P=0.58，提示TATA 7等位基因与非印度高加索人群新生儿黄疸无显著关联性；与印度人群新生儿黄疸有显著关联性，OR=2.32，95%CI：1.36～3.97，P=0.002(图12)。

图12 高加索人群TATA 7等位基因与新生儿黄疸关联性的敏感性分析

Fig12 Sensitivity analysis of the association between 7 allele of TATA box and hyperbilirubinemia in Caucasian

3 讨论

3.1 纳入文献的质量 目前对于RCT研究的Meta分析可采用Cochrane Handbook for Systematic Reviews of Interventions 5.0.2 的质量评价标准对纳入文献的质量进行评价；对于诊断性研究的Meta分析可采用Whiting 等[33]制订的 QUADAS(quality assessment of diagnostic accuracy studies)14条标准评价文献质量。但对于遗传关联性的研究，目前尚无统一和公认的质量评价体系；加强性报告遗传学关联研究(strengthening the reporting of genetic association studies，STREGA)写作清单是目前国际公认的撰写遗传关联性研究论文的规范，但其条目较多，评价费时，不适用于评价遗传关联性研究的质量。本研究采用的评价方法是根据NHI-NHGRI研究工作组2007年制定的遗传关联性研究报告规范(发表于Nature)总结的关于遗传关联性研究的14条质量评价方法，可操作性强，可客观和全面评价遗传关联性的研究质量。

本研究文献质量评价结果显示：①研究设计和一般情况描述：仅4/22篇文献进行了样本量计算；1/22篇文献对照组特征描述不充分；10/22篇文献未描述病例组国家、地区、人种学信息及临床资料。②实验方法：13/22篇文献给出了多态性的鉴定序列或测序结果图；仅5/22篇文献用相同或不同方法进行了基因型结果的检测；所有文献均未描述是否采用盲法。③统计学方法：所有文献均未进行连锁不平衡分析；仅1/22篇文献所提供的数据不能进行Hardy-Weinberg平衡检测；2/22篇文献进行了协变量调整；7/22篇文献进行了危险因素分析。④其他：所有文献均未在其他人群中进行结果重复；仅1/22篇文献进行了功能检测。提示，本Meta分析纳入的文献在研究设计、实验方法、统计学方法等方面与国际一流杂志对遗传关联性研究的要求相比存在一定的不足，文献质量总体上偏低。

3.2UGT1A1基因GLY71ARG和TATA多态性与新生儿黄疸的关联性 本Meta分析结果显示，中国、日本、马来西亚和泰国人群GLY71ARG基因型A/A+G/A、等位基因A频率在病例组与对照组间差异均有统计学意义，且研究间均具同质性。提示UGT1A1基因第1外显子nt211为纯合或杂合A，与中国、日本、马来西亚和泰国人群新生儿黄疸存在显著关联性。中国和日本人群TATA基因型7/7+6/7频率在病例组与对照组间差异有统计学意义(病例组显著低于对照组)，马来西亚人群基因型7/7+6/7、7/7及等位基因7频率在病例组与对照组间差异无统计学意义；且研究间均具同质性，提示TATA重复多态性与新生儿黄疸无显著关联性。

高加索人群GLY71ARG基因型A/A+G/A、A/A及等位基因A频率在病例组与对照组差异均无统计学意义，提示GLY71ARG多态性与高加索人群新生儿黄疸无显著关联性；但纳入文献仅1篇，样本量较局限，因此尚不能下此结论。高加索人群TATA基因型7/7+6/7和等位基因7频率的Meta分析存在显著的异质性，剔除印度人群后，各研究间均具同质性；对于非印度高加索人群基因型7/7+6/7、7/7及等位基因7频率在病例组与对照组间差异无统计学意义，提示与新生儿黄疸无关联性；印度人群仅1篇文献，虽然Meta分析结果提示TATA重复多态性与新生儿黄疸有关联性，但结论不可靠，需进一步扩大样本量后明确。

3.3 关于遗传关联性研究的建议 ①实验设计：在研究设计最初应进行样本量的计算，明确把握度；全面清晰阐述病例组人种、临床相关资料；对照组与病例组来源于同一人群，也应具备相同人种及临床资料；病例组与对照组间各项资料的匹配情况应说明。②实验方法：给出基因型检测的序列或测序结果图；用相同或者不同方法进行结果验证；对实验人员及结果评价人员应采用盲法。③统计学方法：应进行连锁不平衡及Hardy-Weinberg平衡检测；在数据分析时进行人群分层；对于多个位点比较采用多重检测校正，调整α水平；采用协变量调整；进行危险因素的回归分析，给出OR值及其95%CI。④其他：在其他人群中进行验证，并进行功能学的研究。

3.4 本文的局限性和不足之处 ①纳入文献的质量偏低，可能与本文所采用的证据评价标准晚于各研究及要求标准较高有关；②对于遗传关联性研究，虽然本文纳入的总体样本量病例组和对照组均有1 000余例，但样本量仍较局限；③纳入的部分文献未详细描述病例组和对照组的特征，因此在人群分层分析时无法完全按照种族分析；④关于GLY71ARG多态性研究的文献中，涉及高加索人种及其他人种的文献较少，符合纳入标准的仅1篇，故而影响了这部分结果的可靠性；⑤对于中英文文献，由于语言限制，未进行评价，仅纳入中英文文献，可能存在发表偏倚和语种偏倚；⑥在信息采集、文献评价中主要以1名评价员进行，缺乏相互检验。

3.5 结论 现有证据表明，UGT1A1基因GLY71ARG多态性与中国、日本、马来西亚和泰国人群新生儿黄疸有显著关联性；启动子TATA重复多态性可能与印度人的新生儿黄疸有关联性，与中国、日本、马来西亚和泰国、非印度高加索人群新生儿黄疸无显著关联性。

致谢 新疆医科大学公共卫生学院严卫丽教授给予本文方法学的指导。

[1]Tukey RH, Strassburg CP.Human UDP-glucuronosyltransfera-ses: metabolism, expression, and disease.Annu Rev Pharmacol Toxicol, 2000,40:581-616

[2]Lin Z, Fontaine J, Watchko JF.Coexpression of gene polymorphisms involved in bilirubin production and metabolism.Pediatrics, 2008 ,122(1):156-162

[3]Kaplan M, Hammerman C, Maisels MJ.Bilirubin genetics for the nongeneticist: hereditary defects of neonatal bilirubin conjugation.Pediatrics, 2003, 111 (4 p t 1) : 886-893

[4]Sun G, Wu M, Cao J, et al.Cord blood bilirubin level in relation to bilirubin UDP-glucuronosyltransferase gene missense allele in Chinese neonates.Acta Paediatr, 2007,96(11):1622-1625

[5]Setia S, Villaveces A, Dhillon P, et al.Neonatal jaundice in Asian, white, and mixed-race infants.Arch Pediatr Adolesc Med, 2002,156(3):276-279

[6]中华医学会儿科学分会新生儿学组. 新生儿黄疸干预推荐方案.Chin J Pediatr(中华儿科杂志),2001,39 (3):185-187

[7]American Academy of Pediatrics, Subcommittee on Neonatal Hyperbilirubinemia. Neonatal jaundice and kernicterus. Pediatrics,2001,108(3):763-765

[8]Bhutani VK, Johnson LH, Jeffrey Maisels M,et al.Kernicterus: epidemiological strategies for its prevention through systems-based approaches.J Perinatol,2004,24(10):650-662

[9]Chanock SJ, Manolio T, Boehnke M, et al.Replicating genotype-phenotype associations.Nature,2007,447(7145):655-660

[10]Flores C, Pino-Yanes Mdel M, Villar J.A quality assessment of genetic association studies supporting susceptibility and outcome in acute lung injury.Crit Care, 2008,12(5):130

[11]Fu WP(傅雯萍), Zhao X, Qu YM.Effect of UDP-glucoronosyl transferase transferase 1A1 gene TATA mutation on the development of breast-milk jaundice. Acta Academiae Medicinae CPAF(武警医学院学报),2009,18(11):926-928

[12]Tian GY(田桂英), Xu FS, Zhu FX, et al.Association of Gly71Arg mutation in gene of bilirubin uridine 5′-diphosphate-glucuronosyl transferase and neonatal jaundice. J Appl Clin Pediatr(实用儿科临床杂志),2008,23(2):129-130

[13]Zhong DN(钟丹妮), Liu Y, Lin WX. The relationship between gene mutation of UGT1A1 and neonatal hyperbilirubinemia in Guangxi region. Journal of Guangxi Medical University(广西医科大学学报),2008,25(6):838-840

[14]Gao ZY(高宗燕). UGT1A1基因突变对广西新生儿黄疸发病的影响.广西医科大学.2008

[15]Yusoff S, Takeuchi A, Ashi C, et al. A polymorphic mutation, c.-3279T > G, in the UGT1A1 promoter is a risk factor for neonatal jaundice in the Malay Population. Pediatr Res, 2010,67(4):401-406

[16]Agrawal SK, Kumar P, Rathi R, et al.UGT1A1 gene polymorphisms in North Indian neonates presenting with unconjugated hyperbilirubinemia.Pediatr Res,2009,65(6):675-680

[17]Prachukthum S, Nunnarumit P, Pienvichit P,et al.Genetic polymorphisms in Thai neonates with hyperbilirubinemia.Acta Paediatr, 2009,98(7):1106-1110

[18]Watchko JF, Lin Z, Clark RH, et al.Complex multifactorial nature of significant hyperbilirubinemia in neonates.Pediatrics, 2009 ,124(5):868-877

[19]Chang PF, Lin YC, Liu K, et al.Prolonged unconjugated hyperbiliriubinemia in breast-fed male infants with a mutation of uridine diphosphate-glucuronosyl transferase.J Pediatr, 2009,155(6):860-863

[20]Boo NY, Wong FL, Wang MK, et al.Homozygous variant of UGT1A1 gene mutation and severe neonatal hyperbilirubinemia.Pediatr Int, 2009,51(4):488-493

[21]Muslu N, Turhan AB, Eskandari G, et al.The frequency of UDP-glucuronosyltransferase 1A1 promoter region (TA)7 polymorphism in newborns and it′s relation with jaundice.J Trop Pediatr,2007 ,53(1):64-68

[22]Wong FL, Wang MK, Boo NY, et al.Rapid detection of the UGT1A1 single nucleotide polymorphism G211A using real-time PCR with Taqman minor groove binder probes.J Clin Lab Anal,2007,21(3):167-172

[23]Babaoglu MO, Yigit S, Aynacioglu AS, et al.Neonatal jaundice and bilirubin UDP-glucuronosyl transferase 1A1 gene polymorphism in Turkish patients.Basic Clin Pharmacol Toxicol, 2006,98(4):377-380

[24]Yusoff S, Van Rostenberghe H, Yusoff N, et al. Frequencies of A(TA)(7)TAA, G71R, and G493R mutations of the UGT1A1 gene in the Malaysian population.Biol Neonate,2006,89(3):171-176

[25]Kanai M, Kijima K, Shirahata E, et al.Neonatal hyperbilirubinemia and the bilirubin uridine diphosphate-glucuronosyltransferase gene: the common -3263T > G mutation of phenobarbital response enhancer module is not associated with the neonatal hyperbilirubinemia in Japanese.Pediatr Int, 2005,47(2):137-141

[26]Huang MJ, Kua KE, Teng HC, et al.Risk factors for severe hyperbilirubinemia in neonates.Pediatr Res, 2004,56(5):682-689

[27]Sutomo R, Talib NA, Yusoff NM, et al.Screening for G71Rmutation of the UGT1A1 gene in the Javanese-Indonesian and Malay-Malaysian populations.Pediatr Int, 2004 ,46(5):565-569

[28]Ulgenalp A, Duman N, Schaefer FV, et al.Analyses of polymorphism for UGT1*1 exon 1 promoter in neonates with pathologic and prolonged jaundice.Biol Neonate, 2003,83(4):258-262

[29]Huang CS, Chang PF, Huang MJ, et al.Relationship between bilirubin UDP-glucuronosyl transferase 1A1 gene and neonatal hyperbilirubinemia.Pediatr Res,2002,52(4):601-605

[30]Yamamoto A, Nishio H, Waku S, et al.Gly71Arg mutation of the bilirubin UDP-glucuronosyltransferase 1A1 gene is associated with neonatal hyperbilirubinemia in the Japanese population.Kobe J Med Sci, 2002,48(3-4):73-77

[31]Maruo Y, Nishizawa K, Sato H, et al.Association of neonatal hyperbilirubinemia with bilirubin UDP-glucuronosyltransferase polymorphism.Pediatrics, 1999,103(6 Pt 1):1224-1227

[32]Akaba K, Kimura T, Sasaki A, et al.Neonatal hyperbilirubinemia and mutation of the bilirubin uridine diphosphate-glucuronosyltransferase gene: a common missense mutation among Japanese, Koreans and Chinese.Biochem Mol Biol Int, 1998,46(1):21-26

[33]Whiting P, Rutjes AW, Reitsma JB, et al.The development of QUADAS: a tool for the quality assessment of studies of diagnostic accuracy included in systematic reviews.BMC Med Res Methodol,2003 ,3:25