{{short description|Statistical method used in genetic association studies}} '''Genomic control (GC)''' is a statistical method that is used to control for the confounding effects of population stratification in genetic association studies. The method was originally outlined by Bernie Devlin and Kathryn Roeder in a 1999 paper.<ref>{{Cite journal |last1=Devlin |first1=Bernie |author-link=Bernie Devlin |last2=Roeder |first2=Kathryn |author-link2=Kathryn Roeder |date=1999 |title=Genomic Control for Association Studies |journal=Biometrics |language=en |volume=55 |issue=4 |pages=997–1004 |doi=10.1111/j.0006-341X.1999.00997.x |pmid=11315092 |issn=1541-0420|citeseerx=10.1.1.420.1751 |s2cid=6297807 }}</ref> It involves using a set of anonymous genetic markers to estimate the effect of population structure on the distribution of the chi-square statistic. The distribution of the chi-square statistics for a given allele that is suspected to be associated with a given trait can then be compared to the distribution of the same statistics for an allele that is expected not to be related to the trait.<ref>{{Cite journal |last1=Donnelly |first1=Peter |author-link=Peter Donnelly |last2=Phillips |first2=Michael S. |last3=Cardon |first3=Lon R. |last4=Marchini |first4=Jonathan |author-link4=Jonathan Marchini |date=May 2004 |title=The effects of human population structure on large genetic association studies |journal=Nature Genetics |language=en |volume=36 |issue=5 |pages=512–517 |doi=10.1038/ng1337 |pmid=15052271 |issn=1546-1718|doi-access=free }}</ref><ref>{{Cite journal |last1=Altshuler |first1=David |author-link=David Altshuler (physician) |last2=Hirschhorn |first2=Joel N. |author-link2=Joel N. Hirschhorn |last3=Henderson |first3=Brian |last4=Sklar |first4=Pamela |author-link4=Pamela Sklar |last5=Lander |first5=Eric S. |author-link5=Eric Lander |last6=Kolonel |first6=Laurence N. |last7=Petryshen |first7=Tracey L. |last8=Pato |first8=Michele T. |last9=Pato |first9=Carlos N. |date=April 2004 |title=Assessing the impact of population stratification on genetic association studies |journal=Nature Genetics |language=en |volume=36 |issue=4 |pages=388–393 |doi=10.1038/ng1333 |pmid=15052270 |issn=1546-1718|doi-access=free }}</ref> The method is supposed to involve the use of markers that are not linked to the marker being tested for a possible association.<ref>{{Cite journal |last1=Krawczak |first1=Michael |last2=Dempfle |first2=Astrid |last3=Lieb |first3=Wolfgang |last4=Freitag-Wolf |first4=Sandra |last5=Yadav |first5=Pankaj |date=2015-10-01 |title=Allowing for population stratification in case-only studies of gene–environment interaction, using genomic control |journal=Human Genetics |language=en |volume=134 |issue=10 |pages=1117–1125 |doi=10.1007/s00439-015-1593-y |pmid=26297539 |s2cid=18146948 |issn=1432-1203}}</ref> In theory, it takes advantage of the tendency of population structure to cause overdispersion of test statistics in association analyses.<ref>{{Cite journal |last1=Devlin |first1=Bernie |last2=Roeder |first2=Kathryn |last3=Wasserman |first3=Larry |s2cid=11547174 |author-link3=Larry A. Wasserman |date=2001-11-01 |title=Genomic Control, a New Approach to Genetic-Based Association Studies |journal=Theoretical Population Biology |language=en |volume=60 |issue=3 |pages=155–166 |doi=10.1006/tpbi.2001.1542 |pmid=11855950 |issn=0040-5809|url=https://figshare.com/articles/journal_contribution/6586649 }}</ref> The genomic control method is as robust as family-based designs, despite being applied to population-based data.<ref>{{Cite journal |last1=Roeder |first1=Kathryn |last2=Devlin |first2=B. |last3=Bacanu |first3=Silviu-Alin |date=2000-06-01 |title=The Power of Genomic Control |url= |journal=The American Journal of Human Genetics |language=English |volume=66 |issue=6 |pages=1933–1944 |doi=10.1086/302929 |issn=1537-6605 |pmc=1378064 |pmid=10801388}}</ref> It has the potential to lead to a decrease in statistical power to detect a true association, and it may also fail to eliminate the biasing effects of population stratification.<ref>{{Cite journal |last1=Greenberg |first1=David A. |last2=Zhang |first2=Junying |last3=Shmulewitz |first3=Dvora |date=2004 |title=Case-Control Association Studies in Mixed Populations: Correcting Using Genomic Control |journal=Human Heredity |language=english |volume=58 |issue=3–4 |pages=145–153 |doi=10.1159/000083541 |issn=1423-0062 |pmid=15812171|s2cid=24635575 }}</ref> A more robust form of the genomic control method can be performed by expressing the association being studied as two Cochran–Armitage trend tests, and then applying the method to each test separately.<ref>{{Cite journal |last1=Gastwirth |first1=Joseph L. |last2=Freidlin |first2=Boris |last3=Zheng |first3=Gang |date=2006-02-01 |title=Robust Genomic Control for Association Studies |url= |journal=The American Journal of Human Genetics |language=English |volume=78 |issue=2 |pages=350–356 |doi=10.1086/500054 |issn=1537-6605 |pmc=1380242 |pmid=16400614}}</ref>
The assumption of population homogeneity in association studies, especially case-control studies, can easily be violated and can lead to both type I and type II errors. It is therefore important for the models used in the study to compensate for the population structure. The problem in case control studies is that if there is a genetic involvement in the disease, the case population is more likely to be related than the individuals in the control population. This means that the assumption of independence of observations is violated. Often this will lead to an overestimation of the significance of an association but it depends on the way the sample was chosen. If, coincidentally, there is a higher allele frequency in a subpopulation of the cases, you will find association with any trait that is more prevalent in the case population.<ref>{{cite journal | vauthors = Lander ES, Schork NJ | title = Genetic dissection of complex traits | journal = Science | volume = 265 | issue = 5181 | pages = 2037–48 | date = September 1994 | pmid = 8091226 | doi = 10.1126/science.8091226 | bibcode = 1994Sci...265.2037L }}</ref> This kind of spurious association increases as the sample population grows so the problem should be of special concern in large scale association studies when loci only cause relatively small effects on the trait. A method that in some cases can compensate for the above described problems has been developed by Devlin and Roeder (1999).<ref name="devlin_roeder1999">{{cite journal | vauthors = Devlin B, Roeder K | title = Genomic control for association studies | journal = Biometrics | volume = 55 | issue = 4 | pages = 997–1004 | date = December 1999 | pmid = 11315092 | doi = 10.1111/j.0006-341X.1999.00997.x | s2cid = 6297807 }}</ref> It uses both a frequentist and a Bayesian approach (the latter being appropriate when dealing with a large number of candidate genes).
The frequentist way of correcting for population structure works by using markers that are not linked with the trait in question to correct for any inflation of the statistic caused by population structure. The method was first developed for binary traits but has since been generalized for quantitative ones.<ref>{{cite journal | vauthors = Bacanu SA, Devlin B, Roeder K | title = Association studies for quantitative traits in structured populations | journal = Genetic Epidemiology | volume = 22 | issue = 1 | pages = 78–93 | date = January 2002 | pmid = 11754475 | doi = 10.1002/gepi.1045 | s2cid = 3053350 | url = https://figshare.com/articles/journal_contribution/6586442 }}</ref> For the binary one, which applies to finding genetic differences between the case and control populations, Devlin and Roeder (1999) use Armitage's trend test
: <math> Y^2=\frac{N(N(r_1+2r_2)-R(n_1+2n_2))^2}{R(N-R)(N(n_1 + 4n_2) - (n_1 + 2n_2)^2)} </math>
and the <math>\chi^2</math> test for allelic frequencies
: <math> \chi^2\sim X_A^2 = \frac{2N (2N(r_1 + 2r_2) - R(n_1 + 2n_2))^2} {4R(N - R) (2N(n_1 + 2n_2) - (n_1 + 2n_2)^2)} </math>
{| class="wikitable" |- ! Alleles ! aa ! Aa ! AA ! total |- | Case | ''r''<sub>0</sub> | ''r''<sub>1</sub> | ''r''<sub>2</sub> | ''R'' |- | Control | ''s''<sub>0</sub> | ''s''<sub>1</sub> | ''s''<sub>2</sub> | ''S'' |- | total | ''n''<sub>0</sub> | ''n''<sub>1</sub> | ''n''<sub>2</sub> | ''N'' |}
If the population is in Hardy–Weinberg equilibrium the two statistics are approximately equal. Under the null hypothesis of no population stratification the trend test is asymptotic <math>\chi^2</math> distribution with one degree of freedom. The idea is that the statistic is inflated by a factor <math>\lambda</math> so that <math>Y^2\sim\lambda\chi_1^2</math> where <math>\lambda</math> depends on the effect of stratification. The above method rests upon the assumptions that the inflation factor <math>\lambda</math> is constant, which means that the loci should have roughly equal mutation rates, should not be under different selection in the two populations, and the amount of Hardy–Weinberg disequilibrium measured in Wright's coefficient of inbreeding ''F'' should not differ between the different loci. The last of these is of greatest concern. If the effect of the stratification is similar across the different loci <math>\lambda</math> can be estimated from the unlinked markers
: <math>\widehat{\lambda}= \frac{\operatorname{median}(Y_1^2,Y_2^2,\ldots, Y_L^2)}{0.456} </math>
where ''L'' is the number of unlinked markers. The denominator is derived from the gamma distribution as a robust estimator of <math>\lambda</math>. Other estimators have been suggested, for example, Reich and Goldstein<ref>{{cite journal | vauthors = Reich DE, Goldstein DB | title = Detecting association in a case-control study while correcting for population stratification | journal = Genetic Epidemiology | volume = 20 | issue = 1 | pages = 4–16 | date = January 2001 | pmid = 11119293 | doi = 10.1002/1098-2272(200101)20:1<4::AID-GEPI2>3.0.CO;2-T | s2cid = 17480622 }}</ref> suggested using the mean of the statistics instead. This is not the only way to estimate <math>\lambda</math> but according to Bacanu et al.<ref>{{cite journal | vauthors = Bacanu SA, Devlin B, Roeder K | title = The power of genomic control | journal = American Journal of Human Genetics | volume = 66 | issue = 6 | pages = 1933–44 | date = June 2000 | pmid = 10801388 | pmc = 1378064 | doi = 10.1086/302929 }}</ref> it is an appropriate estimate even if some of the unlinked markers are actually in disequilibrium with a disease causing locus or are themselves associated with the disease. Under the null hypothesis and when correcting for stratification using ''L'' unlinked genes, <math>Y^2</math> is approximately <math>\chi^2_1</math> distributed. With this correction the overall type I error rate should be approximately equal to <math>\alpha</math> even when the population is stratified. Devlin and Roeder (1999)<ref name="devlin_roeder1999"/> mostly considered the situation where <math>\alpha=0.05</math> gives a 95% confidence level and not smaller p-values. Marchini et al. (2004)<ref>{{cite journal | vauthors = Marchini J, Cardon LR, Phillips MS, Donnelly P | title = The effects of human population structure on large genetic association studies | journal = Nature Genetics | volume = 36 | issue = 5 | pages = 512–7 | date = May 2004 | pmid = 15052271 | doi = 10.1038/ng1337 | s2cid = 11694537 | doi-access = free }}</ref> demonstrates by simulation that genomic control can lead to an anti-conservative p-value if this value is very small and the two populations (case and control) are extremely distinct. This was especially a problem if the number of unlinked markers were in the order 50−100. This can result in false positives (at that significance level).
==References== {{Reflist}}
Category:Statistical genetics Category:Biostatistical methods