细菌全基因组关联研究的方法与应用

随着测序技术的发展和全基因组序列的不断积累,全基因组关联研究(genome-wide association study, GWAS)在人类复杂疾病研究中取得了丰硕成果,10余年间发现了数以万计的疾病风险因子。同样,GWAS也为探索细菌表型的遗传机制提供了新的工具。自2013年第一项细菌GWAS(bacterial GWAS, BGWAS)工作发表以来,目前已有10多项相关研究报道,分别揭示了细菌宿主适应性、耐药性及毒力等表型的遗传机制,极大加深了人们对细菌遗传、进化及传播等方面的认识。本文对目前BGWAS的研究方法、应用成果及存在的问题进行了总结,并对BGWAS的研究前景进行了展望,旨在为微生物学领域开展BGWAS研究提供参考。     

全基因组关联研究在复杂性疾病遗传研究中的应用

在过去的几年中,人们应用全基因组关联研究(genomewide association studies,GWAS)对多种人类复杂性疾病及性状进行研究,如糖尿病、肿瘤、心血管疾病、神经精神系统疾病、自身免疫性疾病等,且已经鉴定出大量与之密切相关的遗传变异,为进一步探索人类复杂性疾病的遗传特征提供重要线索。但是,由于影响复杂性疾病的因素较多,许多已发现遗传变异对疾病贡献较小,作用机制尚不清楚,现全基因组关联研究亦存在许多问题。今本文就GWAS在复杂性疾病中的应用做一综述,并就其前景做一展望。     

冠心病全基因组关联研究进展

近年来全基因组关联研究在世界范围内发展迅猛,研究者应用全基因组关联研究策略发现了一系列疾病的相关基因或变异,将疾病的基因组研究推向一个新的阶段。冠心病是一种由环境因素和遗传因素共同作用导致的复杂疾病,且是世界范围内死亡和致残的首要原因之一,世界各地的研究者应用此策略发现了候选基因关联研究未曾发现的多个冠心病相关易感区域。文章对近年来世界范围内针对冠心病的全基因组关联研究取得的重要进展进行简要总结,然后就现阶段全基因组关联研究所面临的挑战以及对未来研究的发展趋势进行分析阐述,为进一步探究冠心病的遗传机制提供指导。     

基于基因组关联数据识别阿尔茨海默病相关通路

阿尔茨海默病又称老年性痴呆,是一种复杂的中枢神经系统退行性疾病,本文选取一套阿尔茨海默病全基因组关联分析(GWAS)数据,利用Proxy Gene LD软件进行基因水平上的检验,利用Web Gestalt数据库进行遗传通路分析,识别出320个显著(P0.05)的阿尔茨海默病相关基因、8个显著的KEGG通路和41个显著的GO功能类,这些研究结果对进一步揭示阿尔茨海默病潜在的发病机制具有重要意义。     

全基因组关联研究现状

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

基于高通量测序的全基因组关联研究策略

全基因组关联研究（Genome-wide association study, GWAS）是人类复杂疾病研究的重要组成部分之一,在群体水平检测全基因组范围的遗传变异与可观测性状间的遗传关联。传统的GWAS是以芯片（Array）技术获得高密度的遗传变异,尽管硕果累累,但也存在不少问题。如：所谓的“缺失的遗传力”,即利用关联分析检测达到全基因组水平显著的遗传变异位点只能解释小部分遗传力;在某些性状上不同研究的结果一致性较弱;显著关联的遗传变异位点的功能较难解释等。高通量测序技术,也称第二代测序（Next-generation sequencing, NGS）技术,可以快速、准确地产出高通量的变异位点数据,为解决以上问题提供了可行的方案。基于NGS技术的GWAS方法（NGS-GWAS）,可在一定程度上弥补传统GWAS的不足。文章对NGS-GWAS策略和方法进行了系统性调研,提出了目前较为可行的NGS-GWAS的实施策略和方法,并对NGS-GWAS如何应用于个体化医疗（Personalized medicine, PM）进行了展望。     

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

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

肺结核全基因组关联研究进展

肺结核是由结核分枝杆菌感染引起的一类古老但仍对人类造成巨大影响的传染性疾病。到目前为止, 肺结核依然是由单一病原菌导致死亡人数最多的疾病, 并且随着耐药菌株的出现而呈现死灰复燃之势。近几年, 肺结核全基因组关联研究在世界范围内取得了阶段性成果, 发现了与肺结核相关联的遗传易感位点和区域, 使肺结核的遗传学研究进入了一个崭新的阶段, 为后续肺结核的早期和综合防治提供了重要线索。然而, 由于人群遗传结构差异和宿主/病原体相互作用, 与其他复杂疾病相比, 肺结核全基因组关联研究依旧面临重重困难, 进展缓慢。文章对不同人群肺结核全基因组关联研究及其验证进行综述, 并系统阐述了目前研究中存在的困难及可能的应对策略。     

全基因组关联研究通路分析方法现状

全基因组关联研究(genome-wide association study, GWAS)自2005年首次发表以来已不断增进人们对疾病遗传机制的认识,结合系统生物学并改进统计分析方法是对GWAS数据进行深度挖掘的重要途径。通路分析(pathway analysis)将GWAS所检测的遗传变异根据一定的生物学含义组合为集合进行分析,有利于发现对疾病单独效应小却在通路中相互关联的遗传变异,更有利于进行生物学解释。当前通路分析在GWAS数据上已有较为广泛的应用并取得初步成果。与此同时,通路分析的统计方法仍在不断发展。本文旨在介绍现有直接以SNP为对象的GWAS通路分析算法,根据方法中是否采用核函数分为非核算法和核算法两大类,其中非核算法主要包括基因功能富集分析(gene set enrichment analysis, GSEA)和分层贝叶斯优取(hierarchical Bayes prioritization, HBP),核算法包括线性核(linear kernel, LIN)、状态认证核(identity-by-status kernel, IBS)和尺度不变核(powered exponential kernel)。通过介绍这些方法的计算原理和优缺点,以期为新算法的构建提供更好的思路,为GWAS领域研究方法的选择提供参考。     

创伤弧菌毒力相关因子的全基因组关联分析研究

【目的】创伤弧菌是致死率最高的弧菌物种,但目前尚无在全基因组层面挖掘毒力相关因子的研究。本研究以创伤弧菌分离来源(临床和环境)作为不同表型,通过与260株基因组序列进行关联分析,挖掘毒力相关因子,从而进一步了解创伤弧菌致病因素。【方法】对139株创伤弧菌分离株进行高通量测序,获取其全基因组序列;与公共数据库已公开发表的121株基因组整合,使用pyseer软件进行全基因组关联分析,对与不同分离来源显著相关的基因进行注释和解读。【结果】共发现11个基因与临床分离株显著相关,其中9个是本研究新发现的创伤弧菌潜在毒力相关因子。【结论】本研究使用群体基因组学和统计遗传学方法,在全基因组范围扫描挖掘了创伤弧菌毒力相关因子,为深入揭示该物种致病机制、设计新的疫苗和治疗靶点提供了重要依据。     

Genome-Wide Association Study of Partial Resistance to P. sojae in Wild Soybeans from Heilongjiang Province,China

Phytophthora root rot (PRR) is a destructive disease of soybeans (Glycine max (L.) Merr) caused by Phytophthora sojae (P. sojae). The most effective way to prevent the disease is growing resistant or tolerant varieties. Partial resistance provides a more durable resistance against the pathogen compared to complete resistance. Wild soybean (Glycine soja Sieb. & Zucc.) seems to be an extraordinarily important gene pool for soybean improvement due to its high level of genetic variation. In this study, 242 wild soybean germplasms originating from different regions of Heilongjiang province were used to identify resistance genes to P. sojae race 1 using a genome-wide association study (GWAS). A total of nine significant SNPs were detected, repeatedly associated with P. sojae resistance and located on chromosomes 1, 10, 12, 15, 17, 19 and 20. Among them, seven favorable allelic variations associated with P. sojae resistance were evaluated by a t-test. Eight candidate genes were predicted to explore the mechanistic hypotheses of partial resistance, including Glysoja.19G051583, which encodes an LRR receptor-like serine/threonine protein kinase protein, Glysoja.19G051581, which encodes a receptor-like cytosolic serine/threonine protein kinase protein. These findings will provide additional insights into the genetic architecture of P. sojae resistance in a large sample of wild soybeans and P. sojae-resistant breeding through marker-assisted selection.     

Penalized Multimarker vs. Single-Marker Regression Methods for Genome-Wide Association Studies of Quantitative Traits

The data from genome-wide association studies (GWAS) in humans are still predominantly analyzed using single-marker association methods. As an alternative to single-marker analysis (SMA), all or subsets of markers can be tested simultaneously. This approach requires a form of penalized regression (PR) as the number of SNPs is much larger than the sample size. Here we review PR methods in the context of GWAS, extend them to perform penalty parameter and SNP selection by false discovery rate (FDR) control, and assess their performance in comparison with SMA. PR methods were compared with SMA, using realistically simulated GWAS data with a continuous phenotype and real data. Based on these comparisons our analytic FDR criterion may currently be the best approach to SNP selection using PR for GWAS. We found that PR with FDR control provides substantially more power than SMA with genome-wide type-I error control but somewhat less power than SMA with Benjamini–Hochberg FDR control (SMA-BH). PR with FDR-based penalty parameter selection controlled the FDR somewhat conservatively while SMA-BH may not achieve FDR control in all situations. Differences among PR methods seem quite small when the focus is on SNP selection with FDR control. Incorporating linkage disequilibrium into the penalization by adapting penalties developed for covariates measured on graphs can improve power but also generate more false positives or wider regions for follow-up. We recommend the elastic net with a mixing weight for the Lasso penalty near 0.5 as the best method.     

小麦根部耐盐性状全基因组关联分析

为了解小麦耐盐相关性状的遗传机理,挖掘与小麦耐盐性显著相关的SNP位点及候选基因,本研究利用浓度200 mmol/L的NaCl溶液和正常营养液对全国300份小麦品种(系)进行耐盐性试验,并利用小麦90 K芯片对分布于小麦全基因组的16650个SNP,采用Q+K关联混合模型对小麦最长根长、根干重、根鲜重、根平均直径、根尖...     

普通菜豆镰孢菌枯萎病抗性种质资源筛选及全基因组关联分析

普通菜豆镰孢菌枯萎病是严重制约菜豆(Phaseolus vulgaris)产量的主要病害之一。采用下胚轴双孔注射法对601份普通菜豆种质资源进行枯萎病抗性鉴定, 共筛选出4份高抗材料。在此基础上, 基于分布在全基因组上的3 765 456个单核苷酸多态性(SNP)标记, 进行全基因组关联分析, 以P<1×10^-5为阈值。结果检测到57个显著关联的SNP位点, 分布于1、2、6、8和11号染色体上; 共获得8个显著关联区域, 其中位于1号染色体上的区域1包含SNP最多(48个), 最显著SNP P值为2.18E-07。在8个显著关联区域中, 共检测到186个基因, 其中157个基因有注释信息, 编码过氧化物酶、抗病蛋白、转录因子和蛋白激酶等。结合KEGG富集分析和序列同源性比对, 鉴定出9个候选基因可能与抗性相关。     

Association study highlights the influence of ELOVL fatty acid elongase 6 gene region on backfat fatty acid composition in Large White pig breed

Dietary fatty acid (FA) composition has an impact on human health. There is an increasing request from consumers for healthier food and pork industry must respond to it without worsening performance and the technological properties of pork products. The inclusion of genetic markers for carcass FA composition in pig selection schemes could be a useful tool to reach the right balance between unsaturated and saturated FAs to satisfy market demands. With the aim of finding genomic regions associated with porcine backfat FA composition, a genome-wide association study was performed on 798 Italian Large White pigs genotyped using Illumina PorcineSNP60 k. The strongest associations with backfat contents of palmitic, palmitoleic, oleic, medium-chain and long-chain FAs were found for the Sus scrofa chromosome (SSC) 8 region located at 119 to 122 Mb, where the gene ELOVL FA elongase 6 is mapped. Palmitic, palmitoleic, stearic and oleic acid contents were also found associated with SSC14, in particular with the genomic region at 121 to 124 Mb, where stearoyl-CoA desaturase Δ9 gene lies. On the other hand, the genomic regions associated with backfat contents of arachidic, arachidonic, n-6 and n-3 FAs showed to harbour mainly genes involved in dietary lipids and carbohydrates digestion, absorption and utilisation. To our knowledge, this is the first study performed in Large White pigs identifying markers and genomic regions associated with backfat FA composition. The results validate in Large White some associations previously detected in other pig breeds and indicate the involvement of distinct metabolic pathways in the deposition pattern of essential and non-essential FAs.     

A Foundation for Provitamin A Biofortification of Maize: Genome-Wide Association and Genomic Prediction Models of Carotenoid Levels

Efforts are underway for development of crops with improved levels of provitamin A carotenoids to help combat dietary vitamin A deficiency. As a global staple crop with considerable variation in kernel carotenoid composition, maize (Zea mays L.) could have a widespread impact. We performed a genome-wide association study (GWAS) of quantified seed carotenoids across a panel of maize inbreds ranging from light yellow to dark orange in grain color to identify some of the key genes controlling maize grain carotenoid composition. Significant associations at the genome-wide level were detected within the coding regions of zep1 and lut1, carotenoid biosynthetic genes not previously shown to impact grain carotenoid composition in association studies, as well as within previously associated lcyE and crtRB1 genes. We leveraged existing biochemical and genomic information to identify 58 a priori candidate genes relevant to the biosynthesis and retention of carotenoids in maize to test in a pathway-level analysis. This revealed dxs2 and lut5, genes not previously associated with kernel carotenoids. In genomic prediction models, use of markers that targeted a small set of quantitative trait loci associated with carotenoid levels in prior linkage studies were as effective as genome-wide markers for predicting carotenoid traits. Based on GWAS, pathway-level analysis, and genomic prediction studies, we outline a flexible strategy involving use of a small number of genes that can be selected for rapid conversion of elite white grain germplasm, with minimal amounts of carotenoids, to orange grain versions containing high levels of provitamin A.     

Unraveling the Genetic Etiology of Adult Antisocial Behavior: A Genome-Wide Association Study

Crime poses a major burden for society. The heterogeneous nature of criminal behavior makes it difficult to unravel its causes. Relatively little research has been conducted on the genetic influences of criminal behavior. The few twin and adoption studies that have been undertaken suggest that about half of the variance in antisocial behavior can be explained by genetic factors. In order to identify the specific common genetic variants underlying this behavior, we conduct the first genome-wide association study (GWAS) on adult antisocial behavior. Our sample comprised a community sample of 4816 individuals who had completed a self-report questionnaire. No genetic polymorphisms reached genome-wide significance for association with adult antisocial behavior. In addition, none of the traditional candidate genes can be confirmed in our study. While not genome-wide significant, the gene with the strongest association (p-value = 8.7×10⁻⁵) was DYRK1A, a gene previously related to abnormal brain development and mental retardation. Future studies should use larger, more homogeneous samples to disentangle the etiology of antisocial behavior. Biosocial criminological research allows a more empirically grounded understanding of criminal behavior, which could ultimately inform and improve current treatment strategies.     

Insights into Pancreatic Cancer Etiology from Pathway Analysis of Genome-Wide Association Study Data

Background

Pancreatic cancer is the fourth leading cause of cancer death in the U.S. and the etiology of this highly lethal disease has not been well defined. To identify genetic susceptibility factors for pancreatic cancer, we conducted pathway analysis of genome-wide association study (GWAS) data in 3,141 pancreatic cancer patients and 3,367 controls with European ancestry.

Methods

Using the gene set ridge regression in association studies (GRASS) method, we analyzed 197 pathways identified from the Kyoto Encyclopedia of Genes and Genomes database. We used the logistic kernel machine (LKM) test to identify major contributing genes to each pathway. We conducted functional enrichment analysis of the most significant genes (P<0.01) using the Database for Annotation, Visualization, and Integrated Discovery (DAVID).

Results

Two pathways were significantly associated with risk of pancreatic cancer after adjusting for multiple comparisons (P<0.00025) and in replication testing: neuroactive ligand-receptor interaction, (Ps<0.00002), and the olfactory transduction pathway (P = 0.0001). LKM test identified four genes that were significantly associated with risk of pancreatic cancer after Bonferroni correction (P<1×10⁻⁵): ABO, HNF1A, OR13C4, and SHH. Functional enrichment analysis using DAVID consistently found the G protein-coupled receptor signaling pathway (including both neuroactive ligand-receptor interaction and olfactory transduction pathways) to be the most significant pathway for pancreatic cancer risk in this study population.

Conclusion

These novel findings provide new perspectives on genetic susceptibility to and molecular mechanisms of pancreatic cancer.     

A Powerful Pathway-Based Adaptive Test for Genetic Association with Common or Rare Variants

In spite of the success of genome-wide association studies (GWASs), only a small proportion of heritability for each complex trait has been explained by identified genetic variants, mainly SNPs. Likely reasons include genetic heterogeneity (i.e., multiple causal genetic variants) and small effect sizes of causal variants, for which pathway analysis has been proposed as a promising alternative to the standard single-SNP-based analysis. A pathway contains a set of functionally related genes, each of which includes multiple SNPs. Here we propose a pathway-based test that is adaptive at both the gene and SNP levels, thus maintaining high power across a wide range of situations with varying numbers of the genes and SNPs associated with a trait. The proposed method is applicable to both common variants and rare variants and can incorporate biological knowledge on SNPs and genes to boost statistical power. We use extensively simulated data and a WTCCC GWAS dataset to compare our proposal with several existing pathway-based and SNP-set-based tests, demonstrating its promising performance and its potential use in practice.