复杂疾病全基因组关联研究进展—— 研究设计和遗传标记

实现全基因组关联研究(Genome-wide association study, GWA)在数年前还是遗传学家们的梦想, 如今它已经变成了现实。自2005年Science杂志报道了第一项有关年龄相关性(视网膜)黄斑变性全基因组关联研究研究以来, 有关与复杂疾病的全基因组关联研究如雨后春笋般层出不穷。文中介绍了近两年来全基因组关联研究在复杂疾病研究领域内的主要发现、全基因组关联研究设计原理、遗传标记的选择、比较及相关商品信息。最后介绍了人类基因组拷贝数变异的研究进展, 总结了人类全基因组关联研究所取得成就和存在的问题, 并对全基因组关联研究未来的研究重点和要解决的问题进行了展望。     

羊重要性状全基因组关联分析研究进展

全基因组关联分析(genome-wide association study, GWAS)是一种复杂性状功能基因鉴定的分析策略,已成为挖掘畜禽重要经济性状候选基因的重要手段。随着绵羊和山羊基因组完成和公布,以及不同密度的SNP (single nucleotide polymorphism)芯片的推出并进行商业化推广,不仅大大丰富了羊标记辅助选择可利用的分子标记,而且还为开展重要性状的分子机理的探索提供了重要技术支撑。本文主要针对羊角、羊毛、羊奶、生长发育、肉质、繁殖和疾病等重要性状的GWAS研究所用的群体、主要研究方法和研究结果进行了综述,并对GWAS方法研究现状进行了归纳,以期为进一步利用GWAS进行羊的各种性状的遗传基础研究提供参考。     

全基因组关联分析的进展与反思

全基因组关联分析(genomewide association study,GWAS)是应用人类基因组中数以百万计的单核苷酸多态性(single nucleotide polymorphism,SNP)为标记进行病例-对照关联分析,以期发现影响复杂性疾病发生的遗传特征的一种新策略。近年来,随着人类基因组计划和基因组单倍体图谱计划的实施,人们已通过GWAS方法发现并鉴定了大量与人类性状或复杂性疾病关联的遗传变异,为进一步了解控制人类复杂性疾病发生的遗传特征提供了重要的线索。然而,由于造成复杂性疾病/性状的因素较多,而且GWAS研究系统较为复杂,因此目前GWAS本身亦存在诸多的问题。本文将从研究方式、研究对象、遗传标记,以及统计分析等方面,探讨GWAS的研究现状以及存在的潜在问题,并展望GWAS今后的发展方向。     

拷贝数变异的全基因组关联分析

基因组拷贝数变异(copy number variations,CNVs)是指与基因组参考序列相比,基因组中≥1 kb的DNA片段插入、缺失和/或扩增,及其互相组合衍生出的复杂变异．由于其具有分布范围广、可遗传、相对稳定和高度异质性等特点,目前认为,CNVs是一种新的可以作为疾病易感标志的基因组DNA多态性,其变异引起的基因剂量改变可以导致表型改变．最近,一种基于CNVs的新的疾病易感基因鉴定策略——CNV全基因组关联分析开始出现,这一策略和传统的基于单核苷酸多态性的关联分析具有互补性,通过认识基因组结构变异可以认识复杂疾病的分子机制和遗传基础．     

全基因组关联研究现状

在过去的5年中, 全基因组关联研究(Genome-wide association study, GWAS)方法已被证明是研究复杂疾病和性状遗传易感变异的一种有效手段。目前, 各国科学家在多种复杂疾病和性状中开展了大量的GWAS, 对肿瘤、糖尿病、心脏病、神经精神疾病、自身免疫及免疫相关疾病等复杂疾病以及一些常见性状(如身高、体重、血脂、色素等)的遗传易感基因研究取得了重大成果。截止到2010年9月11日, 运用GWAS开展了对近200种复杂疾病/性状的研究, 发现了3 000多个疾病相关的遗传变异。文章就GWAS的发展及其在复杂疾病/性状中的应用做一综述。     

基于简化基因组测序的大黄鱼耐高温性状全基因组关联分析

利用Illumina HiSeqTM 2500测序平台, 对通过高温胁迫实验筛选得到的20尾耐高温和20尾不耐高温的大黄鱼(Larimichthys crocea)进行了简化基因组测序(SLAF-seq), 每个样本的平均测序深度达到10.26×, 共获得419211个高质量的群体单核苷酸多态性(SNP)位点 。利用TASSEL软件的混合线性模型(MLM)进行全基因组关联分析(GWAS), 共筛选到38个与大黄鱼耐高温性状显著相关的SNP位点(P<2.39E–08)。利用BLAST程序定位每个SNP位点在大黄鱼基因组中的位置, 并分析其周围的功能基因。结果在38个SNPs附近共找到26个已知的功能基因, 这些基因主要与细胞转录、代谢、免疫等功能相关。研究结果可为下一步大黄鱼耐高温分子机制解析及耐高温品种的选育提供参考。     

复杂疾病全基因组关联研究进展——遗传统计分析

2005年, Science杂志首次报道了有关人类年龄相关性黄斑变性的全基因组关联研究, 此后有关肥胖、2型糖尿病、冠心病、阿尔茨海默病等一系列复杂疾病的全基因组关联研究被陆续报道, 这一阶段被称为人类全基因组关联研究的第一次浪潮。文章分别介绍了全基因组关联研究统计分析的方法、软件和应用实例; 比较了关联分析中多重检验的P值调整方法, 包括Bonferroni、递减的Bonferroni校正法、模拟运算法和控制错误发现率的方法; 还讨论了人群混杂对关联分析结果可能产生的影响及原理, 以及全基因组关联研究中控制人群混杂的方法的研究进展和应用实例。在全基因组关联研究的第一次浪潮中, 应用经典的遗传统计方法发现了许多基因-表型之间的关联并且能够对这些关联做出解释, 其中包括许多基因组中的未知基因和染色体区域。然而, 全基因组关联研究的继续发展需要进一步阐述基因组内基因之间相互作用、基因-基因之间的复杂作用网络与环境因素的相互作用在复杂疾病发生中的作用, 现有的统计分析方法肯定不能满足需要, 开发更为高级的统计分析方法势在必行。最后, 文章还给出了全基因组关联研究统计分析软件的相关网站信息。     

全基因组关联研究的深度分析策略

2005年至今,全基因组关联研究(Genome-wide association study,GWAS)发现了大量复杂疾病/性状相关变异。近来,科学家们关注的焦点又集中在了如何利用GWAS数据进行深入分析,期待发现更多复杂疾病/性状的易感基因。一些新的策略和方法已经被尝试应用到复杂疾病/性状GWAS的后续研究中,例如深入分析GWAS数据;鉴定新的复杂疾病/性状易感基因/位点;国际合作和Meta分析;易感区域精细定位及测序;多种疾病共同易感基因研究;以及基因型填补,基于通路的关联分析,基因-基因、基因-环境交互作用和上位研究等。这些策略和方法的应用弥补了经典GWAS的一些不足之处,进一步推动了人类对复杂疾病/性状遗传机制的认识。文章对上述研究的策略、方法以及所面临的问题和挑战进行了综述,为读者描绘了GWAS后期工作的一个简要框架。     

基于全基因组关联的中国常见肿瘤遗传病因学研究

肿瘤是基因-环境交互作用引起的复杂性疾病.在同样的环境暴露下,不同遗传背景的个体发生肿瘤的风险有很大差异.研究肿瘤相关遗传因素对理解肿瘤发生发展乃至诊断治疗都有重要意义.近年来发展的全基因组关联研究(genome-wide association study,GWAS)可在全基因组范围内发现与复杂疾病或表型关联的遗传因素,为复杂疾病遗传学研究提供了强有力的手段.欧美研究者运用全基因组关联研究的方法,对各种常见肿瘤进行了研究,获得了重要成果.2010年以来,中国科学家在国际核心期刊发表了一系列高水平的肿瘤全基因组关联研究成果,在中国常见肿瘤的遗传病因学研究方面取得了重要进展.     

中国汉族高原肺水肿易感基因的全基因组关联研究

高原肺水肿(High-altitude pulmonary edema, HAPE)是一种特发于高原低氧环境的肺水肿, 是遗传和环境因素共同作用的结果。为了寻找与中国汉族高原肺水肿相关的单核苷酸多态性(Single nucleotide polymorphism, SNP)位点及易感基因, 文章利用Affymetrix SNP Array 6.0芯片, 对2010年5月至2012年7月在青海省玉树地区执行援建任务时来自平原地区的40例HAPE患者和33例健康对照进行全基因组SNP分型, 通过PLINK软件对芯片结果进行全基因组关联分析(Genome-wide association study, GWAS), 筛选出在病例组和对照组中间有显著差异(P < 10E-7)的SNP位点57个, 通过对57个SNP位点附近74个基因进行GO与Pathway富集分析, 发现这些基因与“前列腺素代谢”、“四烯酸代谢”、“氮代谢”显著相关(adjust P < 0.05), 以上代谢过程与HAPE病理生理机制相关。结果表明, 高原肺水肿受遗传多态性影响, 与多个基因以及位点相关。     

Genetic predisposition for sudden cardiac death in myocardial ischaemia: the Arrhythmia Genetics in the NEtherlandS study

Sudden cardiac death from ventricular fibrillation during myocardial infarction is a leading cause of total and cardiovascular mortality. This multifactorial, complex condition clusters in families, suggesting a substantial genetic cause. We carried out a genomewide association study (GWAS) for sudden cardiac death, in the AGNES (Arrhythmia Genetics in the Netherlands) population, consisting of patients with (cases) and without (controls) ventricular fibrillation during a first ST-elevation myocardial infarction. The most significant association was found at chromosome 21q21 (rs2824292; odds ratio = 1.78, 95% CI 1.47–2.13, P = 3.3 × 10⁻¹⁰), 98 kb proximal of the CXADR gene, encoding the Coxsackie and adenovirus receptor. This locus has not previously been implicated in arrhythmia susceptibility. Further research on the mechanism of this locus will ultimately provide novel insight into arrhythmia mechanisms in this condition.     

多基因遗传病基因研究的策略和方法

基因在决定个体表型方面起着决定性的作用。虽然单基因疾病的致病基因的克隆工作取得了显著的进展,但对于多基因疾病来说,仍然存在许多问题,同时也是巨大的挑战。本文综述了多基因疾病的遗传特点和多基因疾病易感基因识别、分离和克隆的一般步骤和存在的问题,介绍了人类基因组计划在此过程中的作用和单核苷酸多态性的应用前景,提出 了最有可能克隆出多基因疾病易感基因的策略和方法。     

Genetic architecture of type 2 diabetes

Genome-wide association studies (GWAS) have identified over 70 loci associated with type 2 diabetes (T2D). Most genetic variants associated with T2D are common variants with modest effects on T2D and are shared with major ancestry groups. To what extent the genetic component of T2D can be explained by common variants relies upon the shape of the genetic architecture of T2D. Fine mapping utilizing populations with different patterns of linkage disequilibrium and functional annotation derived from experiments in relevant tissues are mandatory to track down causal variants responsible for the pathogenesis of T2D.     

中国北方汉族人群核呼吸因子-1基因+141G/T单核苷酸多态与冠心病发病的关联研究

目的:本研究旨在探讨IRF-1基因+141 G/T单核苷酸多态位点与中国北方汉族人群冠心病发病的相关关系。方法:本研究采用聚合酶链反应-限制性片段长度多态性对经过冠脉造影证实的冠状动脉有一条主要分支狭窄大于70%的675例冠心病患者和经过冠状动脉造影证实冠状动脉狭窄小于20%或完全正常的636例对照患者进行检验检,分析核呼吸因子IRF-1基因+141G/T单核苷酸多态位点的基因型和等位基因频率在两组间的分布情况。结果:核呼吸因子IRF-1基因+141 G/T单核苷酸多态位点三种基因型(GG型,GT型和TT型)在中国北方汉族人群冠心病组的分布频率分别为53.8%,36.2%和10.1%,在对照组的分布频率分别为45.6%,46.2%和8.2%,核呼吸因子IRF-1基因+141 G/T单核苷酸多态位点的基因型和等位基因频率分布在对照组和冠心病组之间存在统计学差异(P0.05)。Logistic回归分别校正冠心病的其他危险因素性别、年龄、体重指数、吸烟、高血压、高脂血症、糖尿病等后,核呼吸因子IRF-1基因+141 G/T单核苷酸多态位点与中国北方汉族人群的冠心病的发病存在相关关系(P0.05)。结论:核呼吸因子IRF-1基因+141 G/T单核苷酸多态与中国北方汉族人群冠心病的发病存在相关关系,IRF-1基因+141 G/T多态可能是中国北方汉族人群冠心病发病的独立危险因子。     

Next generation genome-wide association tool: design and coverage of a high-throughput European-optimized SNP array

The success of genome-wide association studies has paralleled the development of efficient genotyping technologies. We describe the development of a next-generation microarray based on the new highly-efficient Affymetrix Axiom genotyping technology that we are using to genotype individuals of European ancestry from the Kaiser Permanente Research Program on Genes, Environment and Health (RPGEH). The array contains 674,517 SNPs, and provides excellent genome-wide as well as gene-based and candidate-SNP coverage. Coverage was calculated using an approach based on imputation and cross validation. Preliminary results for the first 80,301 saliva-derived DNA samples from the RPGEH demonstrate very high quality genotypes, with sample success rates above 94% and over 98% of successful samples having SNP call rates exceeding 98%. At steady state, we have produced 462 million genotypes per week for each Axiom system. The new array provides a valuable addition to the repertoire of tools for large scale genome-wide association studies.     

Whole genome association study identifies regions of the bovine genome and biological pathways involved in carcass trait performance in Holstein-Friesian cattle

Background

Four traits related to carcass performance have been identified as economically important in beef production: carcass weight, carcass fat, carcass conformation of progeny and cull cow carcass weight. Although Holstein-Friesian cattle are primarily utilized for milk production, they are also an important source of meat for beef production and export. Because of this, there is great interest in understanding the underlying genomic structure influencing these traits. Several genome-wide association studies have identified regions of the bovine genome associated with growth or carcass traits, however, little is known about the mechanisms or underlying biological pathways involved. This study aims to detect regions of the bovine genome associated with carcass performance traits (employing a panel of 54,001 SNPs) using measures of genetic merit (as predicted transmitting abilities) for 5,705 Irish Holstein-Friesian animals. Candidate genes and biological pathways were then identified for each trait under investigation.

Results

Following adjustment for false discovery (q-value < 0.05), 479 quantitative trait loci (QTL) were associated with at least one of the four carcass traits using a single SNP regression approach. Using a Bayesian approach, 46 QTL were associated (posterior probability > 0.5) with at least one of the four traits. In total, 557 unique bovine genes, which mapped to 426 human orthologs, were within 500kbs of QTL found associated with a trait using the Bayesian approach. Using this information, 24 significantly over-represented pathways were identified across all traits. The most significantly over-represented biological pathway was the peroxisome proliferator-activated receptor (PPAR) signaling pathway.

Conclusions

A large number of genomic regions putatively associated with bovine carcass traits were detected using two different statistical approaches. Notably, several significant associations were detected in close proximity to genes with a known role in animal growth such as glucagon and leptin. Several biological pathways, including PPAR signaling, were shown to be involved in various aspects of bovine carcass performance. These core genes and biological processes may form the foundation for further investigation to identify causative mutations involved in each trait. Results reported here support previous findings suggesting conservation of key biological processes involved in growth and metabolism.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-837) contains supplementary material, which is available to authorized users.     

血管紧张素转换酶基因多态性与冠心病发病的相关关系

目的:从分子遗传学角度分析血管紧张素转换酶(ACE)基因I/D多态与中国北方汉族人群中冠心病发病的相关关系。方法:本研究收集在沈阳军区总医院心内科住院的行冠脉造影检查的病例为研究对象,冠脉动脉照影检查显示冠状动脉主支狭窄程度大于等于70%的入选为冠心病组,冠状动脉照影检查显示冠状动脉主支狭窄程度小于20%的为对照组,共入选568名冠心病患者以及性别与年龄相匹配的529名对照个体为研究对象,利用测序的方法分析检测血管紧张素转换酶(ACE)基因I/D多态在冠心病组与对照组中的频率分布情况。结果:血管紧张素转换酶(ACE)基因I/D多态基因型频率符合Hardy-Weinberg定律。血管紧张素转换酶(ACE)基因I/D多态(II型,ID型和DD型)在我们入选的568例冠心病组分布频率分别为50.3%,30.3%和19.4%,而在我们入选的524例对照组中的分布频率为57.7%,31.2%和11.1%,研究发现血管紧张素转换酶(ACE)基因I/D多态可能是中国北方汉族人群冠心病发病的独立危险因素(P0.05);利用多元回归分析发现,在调整了冠心病的其他危险因素后,血管紧张素转换酶(ACE)基因I/D多态的变化仍然是中国北方汉族人群冠心病发病一个独立危险因素,可以预测中国北方汉族人群冠心病的发生。结论:在中国北方汉族人群中,血管紧张素转换酶(ACE)基因I/D多态可能是冠心病发病的独立危险因素,在临床上可以早期预判冠心病的发生。     

中国汉族人群CXCR2基因＋1235C/T多态与急性冠脉综合征的关联研究

目的：炎症反应在动脉粥样斑块变化的病理过程中发挥着重要的作用。本研究探讨CXCR2基因＋1235 C/T单核苷酸多态与中国汉族人群急性冠脉综合征发病的相关关系。方法：本研究采用聚合酶链反应-限制性片段长度多态性方法对675例急性冠脉综合征的患者和636例对照组进行检测,分析CXCR2基因＋1235 C/T单核苷酸多态的基因型和等位基因频率的分布情况,同时收集济南军区总医院心内科经冠脉造影证实为阳性的急性冠脉综合征患者360例及对照者360例,对上述关联分析的结果进行复制实验的印证。结果：CXCR2基因＋1235 C/T单核苷酸多态三种基因型（CC型,CT型和TT型）在急性冠脉综合征组分布频率分别为39.3%,45.3%和15.1%,在对照组分别为41.7%,47.2%和11.1%,CXCR2基因＋1235 C/T基因型和等位基因频率对照组和急性冠脉综合征组之间存在统计学差异（P〈0.05）。Logistic回归校正性别、年龄、体重指数、吸烟、高血压、高脂血症、糖尿病等冠心病的易患因素后,CXCR2基因＋1235 C/T多态与急性冠脉综合征的发病存在相关关系（P〈0.05）。结论：CXCR2基因＋1235 C/T多态与急性冠脉综合征发病存在相关关系,CXCR2基因＋1235 C/T多态可能是中国汉族人群急性冠脉综合征发病的独立危险因子。     

Single-nucleotide polymorphic variants of human glutathione transferase T1-1 differ in stability and functional properties

We have previously expressed hexa-histidine-tagged human glutathione transferase GST T1-1 at very high levels in an Escherichia colilacZ mutagenicity assay strain. Ethylene dibromide (EDB), which is activated by GST T1-1, produces a potent response in the mutation assay. We have now constructed and expressed two SNP variants of wild-type GST T1-1:D141N and E173K. The EDB activation activities of both variant enzymes, as measured by the lacZ mutagenicity assay, are greatly reduced The D141N variant behaved similarly to the wild-type enzyme, in terms of expression level and specific activities for conjugation of glutathione with 1,2-epoxy-3-(p-nitrophenoxy)propane (EPNP), ethylene diiodide (EDI), and 4-nitrobenzyl chloride (NBCl), and for peroxidative detoxication of cumene hydroperoxide (CuOOH). In contrast, variant E173K is poorly expressed, has no detectable activity with EPNP, NBCl, or CuOOH, and has EDI activity much lower than that of the wild-type enzyme. The circular dichroism (CD) thermal denaturation profiles of the wild-type protein and variant D141N show a sharp two-state transition between native and denatured states. Variant E173K showed a very different profile, consistent with improper or incomplete protein folding. Our results show that SNP variants can give rise to GSTT1-1 proteins with significantly altered properties.