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

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

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

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

基于高密度SNP标记估计群体间遗传关联

联合育种的准确性受到群体间遗传关联程度的影响.本研究通过比较基于系谱数据和基因组数据计算的群体遗传关联,探究高密度SNP标记在遗传关联估计中的应用前景.本研究同时使用了模拟数据和真实数据,采用6种不同的遗传关联计算方法,包括PEVD(prediction error variance of differences)、P...     

复杂疾病关联研究中的若干问题

关联研究广泛应用于阐述心血管疾病、2型糖尿病、原发性高血压和肥胖等人类复杂疾病的遗传学基础。文中就关联研究中混杂的识别与控制、候选基因的选择、中间表型的应用、单体型分析方法的应用,以及结果的判定等问题进行了讨论。人群分层是关联研究混杂的主要来源之一。选择患者亲属做对照、基因组对照和选择遗传背景较为一致的隔离人群都可以减少混杂。候选基因的选择可以基于与疾病间的生物学联系或是该基因与疾病某已知相关基因的同源性。适当的应用中间表型和单体型分析方法可以增加关联研究有意义发现的机会。本文认为,优化研究设计、足够的样本含量、正确选择对照,结合先进的数据分析方法,关联研究必将为困扰人类的常见疾病的易感性研究发挥更大的作用。     

SCtool:识别复杂性状和疾病间遗传关联的工具

识别复杂性状和疾病间遗传关联可以提供有用的病因学见解,并有助于确定可能的因果关系的优先级。尽管已有很多工具可以实现复杂性状和疾病间遗传关联,但是某些工具代码可读性差、并且不同工具基于不同的计算机语言、工具间的串联性较差。因此,本研究基于全基因组关联研究(GWAS)数据,提出了SCtool,一个开源、跨平台和用户友好的软件工具。SCtool整合了ldsc, TwosampleMR和MR-BMA三种软件,其主要功能是基于GWAS汇总水平的数据,识别复杂性状和疾病、复杂性状和复杂性状以及疾病与疾病间的遗传相关性并探究其间潜在的因果关联。最后,使用SCtool揭示了全身性铁状态(铁蛋白,血清铁,转铁蛋白,转铁蛋白饱和度)与表观遗传时钟GrimAge之间的遗传关联。     

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

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

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

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

随机SNP在全基因组关联研究人群分层分析中的应用

在复杂疾病的全基因组关联研究中,人群分层现象会增加结果的假阳性率,因此考虑人群遗传结构、控制人群分层是很有必要的。而在人群分层研究中,使用随机选择的SNP的效果还有待进一步探讨。文章利用HapMap Phase2人群中无关个体的Affymetrix SNP 6.0芯片分型数据,在全基因组上随机均匀选择不同数量的SNP,同时利用f值和Fisher精确检验方法筛选祖先信息标记（Ancestry Informative Markers,AIMs）。然后利用HapMap Phase3中的无关个体的数据,以F-statistics和STRUCTURE分析两种方法评估所选出的不同SNP组合对人群的区分效果。研究发现,随机均匀分布于全基因组的SNP可用于识别人群内部存在的遗传结构。文章进一步提示,在全基因组关联研究中,当没有针对特定人群的AIMs时,可在全基因组上随机选择3000以上均匀分布的SNP来控制人群分层。     

后GWAS时代结直肠癌致病SNP功能机制的研究进展

结直肠癌(colorectal cancer, CRC)是受遗传与环境因素共同影响的复杂疾病,其中遗传因素发挥重要作用。至今,全基因组关联研究(genome-wide association studies, GWAS)已经发现了大量与结直肠癌风险相关的遗传变异。随之而来的后GWAS时代,越来越多的研究侧重于利用多组学数据和功能实验对潜在的致病位点进行解析。分析表明绝大多数风险单核苷酸多态性(single nucleotide polymorphism,SNP)位于非编码区,可能通过影响转录因子结合、表观遗传修饰、染色质可及性、基因组高级结构等,调控靶基因表达。本文对后GWAS时代结直肠癌致病位点的机制研究进行综述,阐述了后GWAS对于理解结直肠癌分子机制的重要意义,并探讨了结直肠癌GWAS的应用和前景,为实现GWAS成果转化提供参考。     

新冠肺炎遗传易感基因的研究进展

新型冠状病毒(severe acute respiratory syndrome coronavirus 2,SARS-CoV-2)感染人体后,个体间存在显著不同的新冠肺炎(corona virus disease 2019,COVID-19)临床症状。机体遗传因素在新冠病毒感染后的临床转归过程中发挥重要的作用。以全基因组关联研究(genome-wide association studies, GWAS)为代表的遗传关联研究方法,已成功鉴定了多个与新冠肺炎相关的易感基因,为新冠肺炎防诊治措施的研发提供了理论基础。本文综述了新冠肺炎遗传易感基因的研究进展,包括多种表型、多个人群、多种遗传变异类型的新冠肺炎全基因组关联研究以及易感基因区域的精细定位研究等,旨在为新冠肺炎遗传易感基因的后续研究提供参考。     

Likelihood-based association analysis for nuclear families and unrelated subjects with missing genotype data

Missing data occur in genetic association studies for several reasons including missing family members and uncertain haplotype phase. Maximum likelihood is a commonly used approach to accommodate missing data, but it can be difficult to apply to family-based association studies, because of possible loss of robustness to confounding by population stratification. Here a novel likelihood for nuclear families is proposed, in which distinct sets of association parameters are used to model the parental genotypes and the offspring genotypes. This approach is robust to population structure when the data are complete, and has only minor loss of robustness when there are missing data. It also allows a novel conditioning step that gives valid analysis for multiple offspring in the presence of linkage. Unrelated subjects are included by regarding them as the children of two missing parents. Simulations and theory indicate similar operating characteristics to TRANSMIT, but with no bias with missing data in the presence of linkage. In comparison with FBAT and PCPH, the proposed model is slightly less robust to population structure but has greater power to detect strong effects. In comparison to APL and MITDT, the model is more robust to stratification and can accommodate sibships of any size. The methods are implemented for binary and continuous traits in software, UNPHASED, available from the author.     

A PCA-based method for ancestral informative markers selection in structured populations

Identification of population structure can help trace population histories and identify disease genes. Structured association (SA) is a commonly used approach for population structure identification and association mapping. A major issue with SA is that its performance greatly depends on the informativeness and the numbers of ancestral informative markers (AIMs). Present major AIM selection methods mostly require prior individual ancestry information, which is usually not available or uncertain in practice. To address this potential weakness, we herein develop a novel approach for AIM selection based on principle component analysis (PCA), which does not require prior ancestry information of study subjects. Our simulation and real genetic data analysis results suggest that, with equivalent AIMs, PCA-based selected AIMs can significantly increase the accuracy of inferred individual ancestries compared with traditionally randomly selected AIMs. Our method can easily be applied to whole genome data to select a set of highly informative AIMs in population structure, which can then be used to identify potential population structure and correct possible statistical biases caused by population stratification.     

Properties of structured association approaches to detecting population stratification

OBJECTIVE: To examine the properties of the structured association approach for the detection and correction of population stratification. METHOD: A method is developed, within a latent class analysis framework, similar to the methods proposed by Satten et al. (2001) and Pritchard et al. (2000). A series of simulations illustrate the relative impact of number and type of loci, sample size and population structure. RESULTS: The ability to detect stratification and assign individuals to population strata is determined for a number of different scenarios. CONCLUSION: The results underline the importance of careful marker selection.     

Admixture-matched case-control study: a practical approach for genetic association studies in admixed populations

Case-control genetic association studies in admixed populations are known to be susceptible to genetic confounding due to population stratification. The transmission/disequilibrium test (TDT) approach can avoid this problem. However, the TDT is expensive and impractical for late-onset diseases. Case-control study designs, in which, cases and controls are matched by admixture, can be an appealing and a suitable alternative for genetic association studies in admixed populations. In this study, we applied this matching strategy when recruiting our African American participants in the Study of African American, Asthma, Genes and Environments. Group admixture in this cohort consists of 83% African ancestry and 17% European ancestry, which was consistent with reports from other studies. By carrying out several complementary analyses, our results show that there is a substructure in the cohort, but that the admixture distributions are almost identical in cases and controls, and also in cases only. We performed association tests for asthma-related traits with ancestry, and only found that FEV(1), a measure for baseline pulmonary function, was associated with ancestry after adjusting for socio-economic and environmental risk factors (P=0.01). We did not observe an excess of type I error rate in our association tests for ancestry informative markers and asthma-related phenotypes when ancestry was not adjusted in the analyses. Furthermore, using the association tests between genetic variants in a known asthma candidate gene, beta(2) adrenergic receptor (beta(2)AR) and DeltaFEF(25-75), an asthma-related phenotype, as an example, we demonstrated population stratification was not a confounder in our genetic association. Our present work demonstrates that admixture-matched case-control strategies can efficiently control population stratification confounding in admixed populations.     

Joint analysis for genome-wide association studies in family-based designs

In family-based data, association information can be partitioned into the between-family information and the within-family information. Based on this observation, Steen et al. (Nature Genetics. 2005, 683-691) proposed an interesting two-stage test for genome-wide association (GWA) studies under family-based designs which performs genomic screening and replication using the same data set. In the first stage, a screening test based on the between-family information is used to select markers. In the second stage, an association test based on the within-family information is used to test association at the selected markers. However, we learn from the results of case-control studies (Skol et al. Nature Genetics. 2006, 209-213) that this two-stage approach may be not optimal. In this article, we propose a novel two-stage joint analysis for GWA studies under family-based designs. For this joint analysis, we first propose a new screening test that is based on the between-family information and is robust to population stratification. This new screening test is used in the first stage to select markers. Then, a joint test that combines the between-family information and within-family information is used in the second stage to test association at the selected markers. By extensive simulation studies, we demonstrate that the joint analysis always results in increased power to detect genetic association and is robust to population stratification.     

Inferences from Genomic Models in Stratified Populations

Unaccounted population stratification can lead to spurious associations in genome-wide association studies (GWAS) and in this context several methods have been proposed to deal with this problem. An alternative line of research uses whole-genome random regression (WGRR) models that fit all markers simultaneously. Important objectives in WGRR studies are to estimate the proportion of variance accounted for by the markers, the effect of individual markers, prediction of genetic values for complex traits, and prediction of genetic risk of diseases. Proposals to account for stratification in this context are unsatisfactory. Here we address this problem and describe a reparameterization of a WGRR model, based on an eigenvalue decomposition, for simultaneous inference of parameters and unobserved population structure. This allows estimation of genomic parameters with and without inclusion of marker-derived eigenvectors that account for stratification. The method is illustrated with grain yield in wheat typed for 1279 genetic markers, and with height, HDL cholesterol and systolic blood pressure from the British 1958 cohort study typed for 1 million SNP genotypes. Both sets of data show signs of population structure but with different consequences on inferences. The method is compared to an advocated approach consisting of including eigenvectors as fixed-effect covariates in a WGRR model. We show that this approach, used in the context of WGRR models, is ill posed and illustrate the advantages of the proposed model. In summary, our method permits a unified approach to the study of population structure and inference of parameters, is computationally efficient, and is easy to implement.     

Matching strategies for genetic association studies in structured populations

Association studies in populations that are genetically heterogeneous can yield large numbers of spurious associations if population subgroups are unequally represented among cases and controls. This problem is particularly acute for studies involving pooled genotyping of very large numbers of single-nucleotide-polymorphism (SNP) markers, because most methods for analysis of association in structured populations require individual genotyping data. In this study, we present several strategies for matching case and control pools to have similar genetic compositions, based on ancestry information inferred from genotype data for approximately 300 SNPs tiled on an oligonucleotide-based genotyping array. We also discuss methods for measuring the impact of population stratification on an association study. Results for an admixed population and a phenotype strongly confounded with ancestry show that these simple matching strategies can effectively mitigate the impact of population stratification.     

Quantitative similarity-based association tests using population samples

Although genetic association studies using unrelated individuals may be subject to bias caused by population stratification, alternative methods that are robust to population stratification, such as family-based association designs, may be less powerful. Furthermore, it is often more feasible and less expensive to collect unrelated individuals. Recently, several statistical methods have been proposed for case-control association tests in a structured population; these methods may be robust to population stratification. In the present study, we propose a quantitative similarity-based association test (QSAT) to identify association between a candidate marker and a quantitative trait of interest, through use of unrelated individuals. For the QSAT, we first determine whether two individuals are from the same subpopulation or from different subpopulations, using genotype data at a set of independent markers. We then perform an association test between the candidate marker and the quantitative trait, through incorporation of such information. Simulation results based on either coalescent models or empirical population genetics data show that the QSAT has a correct type I error rate in the presence of population stratification and that the power of the QSAT is higher than that of family-based association designs.     

The bias introduced by population stratification in IBD based linkage analysis

The lack of replication of model-free linkage analyses performed on complex diseases raises questions about the robustness of these methods to various biases. The confounding effect of population stratification on a genetic association study has long been recognized in the genetic epidemiology community. Because the estimation of the number of alleles shared identical by descent (IBD) does not depend on the marker allele frequency when founders of families are observed, model-free linkage analysis is usually thought to be robust to population stratification. However, for common complex diseases, the genotypes of founders are often unobserved and therefore population stratification has the potential to impair model-free linkage analysis. Here, we demonstrate that, when some or all of the founder genotypes are missing, population stratification can introduce deleterious effects on various model-free linkage methods or designs. For an affected sib pair design, it can cause excess false-positive discoveries even when the trait distribution is homogeneous among subpopulations. After incorporating a control group of discordant sib pairs or for a quantitative trait, two circumstances must be met for population stratification to be a confounder: the distributions for both the marker and the trait must be heterogeneous among subpopulations. When this occurs, the bias can result in either a liberal, and hence invalid, test or a conservative test. Bias can be eliminated or alleviated by inclusion of founders' or other family members' genotype data. When this is not possible, new methods need to be developed to be robust to population stratification.     

A unified association analysis approach for family and unrelated samples correcting for stratification

There are two common designs for association mapping of complex diseases: case-control and family-based designs. A case-control sample is more powerful to detect genetic effects than a family-based sample that contains the same numbers of affected and unaffected persons, although additional markers may be required to control for spurious association. When family and unrelated samples are available, statistical analyses are often performed in the family and unrelated samples separately, conditioning on parental information for the former, thus resulting in reduced power. In this report, we propose a unified approach that can incorporate both family and case-control samples and, provided the additional markers are available, at the same time corrects for population stratification. We apply the principal components of a marker matrix to adjust for the effect of population stratification. This unified approach makes it unnecessary to perform a conditional analysis of the family data and is more powerful than the separate analyses of unrelated and family samples, or a meta-analysis performed by combining the results of the usual separate analyses. This property is demonstrated in both a variety of simulation models and empirical data. The proposed approach can be equally applied to the analysis of both qualitative and quantitative traits.