 |
|
How to Interpret a Genome-wide Association Study
Thomas A. Pearson, MD, MPH, PhD;
Teri A. Manolio, MD, PhD
JAMA. 2008;299(11):1335-1344.
ABSTRACT
| |
Genome-wide association (GWA) studies use high-throughput genotyping technologies to assay hundreds of thousands of single-nucleotide polymorphisms (SNPs) and relate them to clinical conditions and measurable traits. Since 2005, nearly 100 loci for as many as 40 common diseases and traits have been identified and replicated in GWA studies, many in genes not previously suspected of having a role in the disease under study, and some in genomic regions containing no known genes. GWA studies are an important advance in discovering genetic variants influencing disease but also have important limitations, including their potential for false-positive and false-negative results and for biases related to selection of study participants and genotyping errors. Although these studies are clearly many steps removed from actual clinical use, and specific applications of GWA findings in prevention and treatment are actively being pursued, at present these studies mainly represent a valuable discovery tool for examining genomic function and clarifying pathophysiologic mechanisms. This article describes the design, interpretation, application, and limitations of GWA studies for clinicians and scientists for whom this evolving science may have great relevance.
INTRODUCTION
In the past 2 years, there has been a dramatic increase in genomic discoveries involvng complex, non-Mendelian diseases, with nearly 100 loci for as many as 40 common diseases robustly identified and replicated in genome-wide association (GWA) studies (T.A.M.; unpublished data, 2008). These studies use high-throughput genotyping technologies to assay hundreds of thousands of the most common form of genetic variant, the single-nucleotide polymorphism (SNP), and relate these variants to diseases or health-related traits.1 Nearly 12 million unique human SNPs have been assigned a reference SNP (rs) number in the National Center for Biotechnology Information's dbSNP database2 and characterized as to specific alleles (alternate forms of the SNP), summary allele frequencies, and other genomic information.3
The GWA approach is revolutionary because it permits interrogation of the entire human genome at levels of resolution previously unattainable, in thousands of unrelated individuals, unconstrained by prior hypotheses regarding genetic associations with disease.4 However, the GWA approach can also be problematic because the massive number of statistical tests performed presents an unprecedented potential for false-positive results, leading to new stringency in acceptable levels of statistical significance and requirements for replication of findings.5
The genome-wide, nonhypothesis-driven nature of GWA studies represents an important step beyond candidate gene studies, in which the high cost of genotyping had limited the number of variants assayed to several hundred at most. This required careful selection of variants to be studied, often based on imperfect understanding of the biologic pathways relating genes to disease.6 Many such associations failed to be replicated in subsequent studies,7-8 leading to calls for all genetic association reports to include documented replication of findings as a prerequisite for publication.9-10
For non-Mendelian conditions, GWA studies also represent a valuable advance over family-based linkage studies, in which multiply affected families are arduously assembled and inheritance patterns are related to several hundred markers throughout the genome. Family-based linkage studies, although successful in identifying genes of large effect in Mendelian diseases such as cystic fibrosis and neurofibromatosis, have had more limited success in common diseases like atherosclerosis and asthma.11 Major limitations of linkage studies are relatively low power for complex disorders influenced by multiple genes, and the large size of the chromosomal regions shared among family members (often comprising hundreds of genes), in whom it can be difficult to narrow the linkage signal sufficiently to identify a causative gene.
GWA studies build on the valuable lessons learned from candidate gene and family linkage studies, as well as the expanding knowledge of the relationships among SNP variants generated by the International HapMap Project,12-13 to capture the great majority of common genetic differences among individuals and relate them to health and disease. These studies not only represent a powerful new tool for identification of genes influencing common diseases, but also use new terminologies (Box 1), apply new models, and present new challenges in interpretation. GWA studies rely on the "common disease, common variant" hypothesis, which suggests that genetic influences on many common diseases will be at least partly attributable to a limited number of allelic variants present in more than 1% to 5% of the population.14 Many important disease-causing variants may be rarer than this and are unlikely to be detected with this approach.
| Box 1. Terms Frequently Used in Genome-wide Association Studies
- Alleles
- Alternate forms of a gene or chromosomal locus that differ in DNA sequence
- Candidate gene
- A gene believed to influence expression of complex phenotypes due to known biological and/or physiological properties of its products,or to its location near a region of association or linkage
- Copy number variants
- Stretches of genomic sequence of roughly 1 kb to 3 Mb in size that are deleted or are duplicated in varying numbers
- False discovery rate 59-60
- Proportion of significant associations that are actually false positives
- False-positive report probability 61
- Probability that the null hypothesis is true, given a statistically significant finding
- Functional studies
- Investigations of the role or mechanism of a genetic variant in causation of a disease or trait
- Gene-environment interactions
- Modification of gene-disease associations in the presence of environmental factors
- Genome-wide association study
- Any study of genetic variation across the entire human genome designed to identify genetic association with observable traits or the presence or absence of a disease, usually referring to studies with genetic marker density of 100 000 or more to represent a large proportion of variation in the human genome
- Genotyping call rate
- Proportion of samples or SNPs for which a specific allele SNP can be reliably identified by a genotyping method
- Haplotype
- A group of specific alleles at neighboring genes or markers that tend to be inherited together
- HapMap 12-13
- Genome-wide database of patterns of common human genetic sequence variation among multiple ancestral population samples
- Hardy Weinberg equilibrium
- Population distribution of 2 alleles (with frequencies p and q) such that the distribution is stable from generation to generation and genotypes occur at frequencies of p2, 2pq, and q2 for the major allele homozygote, heterozygote, and minor allele homozygote, respectively
- Linkage disequilibrium
- Association between 2 alleles located near each other on a chromosome,such that they are inherited together more frequently than expected by chance
- Mendelian disease
- Condition caused almost entirely by a single major gene, such as cystic fibrosis or Huntington's disease, in which disease is manifested in only 1 (recessive) or 2 (dominant) of the 3 possible genotype groups
- Minor allele
- The allele of a biallelic polymorphism that is less frequent in the study population
- Minor allele frequency
- Proportion of the less common of 2 alleles in a population (with 2 alleles carried by each person at each autosomal locus) ranging from less than 1% to less than 50%
- Modest effect
- Association between a gene variant and disease or trait that is statistically significant but carries a small odds ratio (usually <1.5)
- Non-Mendelian disease (also "common" or "complex" disease)
- Condition influenced by multiple genes and environmental factors and not showing Mendelian inheritance patterns
- Nonsynonymous SNP
- A polymorphism that results in a change in the amino acid sequence of a protein (and therefore may affect the function of the protein)
- Platform
- Arrays or chips on which high-throughput genotyping is performed
- Polymorphic
- A gene or site with multiple allelic forms. The term polymorphism usually implies a minor allele frequency of at least 1%
- Population attributable risk
- Proportion of a disease or trait in the population that is due to a specific cause, such as a genetic variant
- Population stratification (also "population structure")
- A form of confounding in genetic association studies caused by genetic differences between cases and controls unrelated to disease but due to sampling them from populations of different ancestries
- Power
- A statistical term for the probability of identifying a difference between 2 groups in a study when a difference truly exists
- Single-nucleotide polymorphism
- Most common form of genetic variation in the genome, in which a single-base substitution has created 2 forms of a DNA sequence that differ by a single nucleotide
- Tag SNP
- A readily measured SNP that is in strong linkage disequilibrium with multiple other SNPs so that it can serve as a proxy for these SNPs on large-scale genotyping platforms
- Trio
- Genetic study design including an affected offspring and both parents
Abbreviation: SNP, single-nucleotide polymorphism.
RETURN TO TEXT |
|
Although GWA discovery studies provide important clues to genomic function and pathophysiologic mechanisms, they are as yet many steps removed from actual clinical application. Nonetheless, they have gained considerable media attention and have the potential for generating queries from patients about whether to get tested for the "new gene for disease X" based on the latest report. In this article, we describe the design, interpretation, application, and limitations of GWA studies for clinicians and scientists for whom this evolving science may have great relevance.
Overview of GWA Studies
A GWA study is defined by the National Institutes of Health as a study of common genetic variation across the entire human genome designed to identify genetic associations with observable traits.15 Although family linkage studies and studies comprising tens of thousands of gene-based SNPs also assay genetic variation across the genome,16 the National Institutes of Health definition requires sufficient density and selection of genetic markers to capture a large proportion of the common variants in the study population, measured in enough individuals to provide sufficient power to detect variants of modest effect.
The present discussion focuses on studies attempting to assay at least 100 000 SNPs selected to serve as proxies for the largest possible number of SNPs.12 The typical GWA study has 4 parts: (1) selection of a large number of individuals with the disease or trait of interest and a suitable comparison group; (2) DNA isolation, genotyping, and data review to ensure high genotyping quality; (3) statistical tests for associations between the SNPs passing quality thresholds and the disease/trait; and (4) replication of identified associations in an independent population sample or examination of functional implications experimentally.
Most of the roughly 100 GWA studies published by the end of 2007 were designed to identify SNPs associated with common diseases. However, the technique can also be used to identify genetic variants related to quantitative traits such as height17 or electrocardiographic QT interval,18 and to rank the relative importance of previously identified susceptibility genes, such as APOE*e4 in Alzheimer disease19 and CARD15 and IL23R in Crohn disease.20
GWA studies can also demonstrate gene-gene interactions, or modification of the association of one genetic variant by another, as with GAB2 and APOE in Alzheimer disease,21 and can detect high-risk haplotypes or combinations of multiple SNPs within a single gene, as in exfoliation glaucoma22 and atrial fibrillation.23 These studies have also been used to identify SNPs associated with gene expression, either as confirmation of a phenotypic association, such as asthma and ORMDL3 expression,24 or more globally.25 Thus, GWA studies have broader applications than those solely involving discovery of individual SNPs associated with discrete disease end points.
Study Designs Used in GWA
By far the most frequently used GWA study design to date has been the case-control design, in which allele frequencies in patients with the disease of interest are compared to those in a disease-free comparison group. These studies are often easier and less expensive to conduct than studies using other designs, especially if sufficient numbers of case and control participants can be assembled rapidly. This design also carries the most assumptions, which if not met, can lead to substantial biases and spurious associations (Table 1). The most important of these biases involve the selected, often unrepresentative nature of the study case participants, who are typically sampled from clinical sources and thus may not include fatal, mild, or silent cases not coming to clinical attention; and the lack of comparability of case and control participants, who may differ in important ways that could be related both to genetic risk factors and to disease outcomes.26
|
|
|
|
Table 1. Study Designs Used in Genome-wide Association Studies
|
|
|
If well-established principles of epidemiologic design are followed, case-control studies can produce valid results that, especially for rare diseases, may not be obtainable in any other way. However, genetic association studies using case-control methodologies have often not always adhered to these principles. The often sharply abbreviated descriptions of case and control participants and lack of comparison of key characteristics in GWA reports27 can make evaluation of potential biases and replication of findings quite difficult.28
The trio design includes the affected case participant and both of his or her parents.29 Phenotypic assessment (classification of affected status) is performed only in the offspring and only affected offspring are included, but genotyping is performed in all 3 trio members. The frequency with which an allele is transmitted to an affected offspring from heterozygous parents is then estimated.29 Under the null hypothesis of no association with disease, the transmission frequency for each allele of a given SNP will be 50%, but alleles associated with the disease will be transmitted in excess to the affected case individual. Because the trio design studies allele transmission from parents to offspring, it is not susceptible to population stratification, or genetic differences between case and control participants unrelated to disease but due to sampling them from populations of different ancestry.30 A significant challenge of the trio design in GWA studies is its sensitivity to even small degrees of genotyping error,4, 31 which can distort transmission proportions between parents and offspring, especially for uncommon alleles. Therefore, standards for genotyping quality in trio studies may need to be more stringent than for other designs.
Cohort studies involve collecting extensive baseline information in a large number of individuals who are then observed to assess the incidence of disease in subgroups defined by genetic variants. Although cohort studies are typically more expensive and take longer to conduct than case-control studies, they often include study participants who are more representative than clinical series of the population from which they are drawn, and they typically include a vast array of health-related characteristics and exposures for which genetic associations can be sought.17-18 For these reasons, genome-wide genotyping has recently been added to cohort studies such as the Framingham Heart Study32 and the Women's Health Study.33
Many GWA studies use multistage designs to reduce the number of false-positive results while minimizing the number of costly genome-wide scans performed and retaining statistical power.4 Genome-wide scans are typically performed on an initial group of case and control participants and then a smaller number of associated SNPs is replicated in a second or third group of case and control participants (Table 2). Some studies begin with small numbers of participants in the initial scan but carry forward large numbers of SNPs to minimize false-negative results.34 Other studies begin with more participants but carry forward a smaller proportion of associated SNPs.35 Optimal proportions of study participants and SNPs in each phase have yet to be determined,36 but carrying forward a small proportion (<5%) of stage 1 SNPs will often mean limiting the associations ultimately identified to those having a relatively large effect.37
|
|
|
|
Table 2. Examples of Multistage Designs in Genome-wide Association Studiesa
|
|
|
Selection of Study Participants
Many genetic studies, whether GWA or otherwise, focus on case participants more likely to have a genetic basis for their disease, such as early-onset cases or those with multiple affected relatives. Misclassification of case participants can markedly reduce study power and bias study results toward no association, particularly when large numbers of unaffected individuals are misclassified as affected. For diseases that are difficult to diagnose reliably, ensuring that cases are truly affected (as by invasive testing or imaging), is probably more important than ensuring generalizability, although the limitations on diagnostic reliability and generalizability should be clearly described so that clinicians can judge the relevance to their patients.
The control participants should be drawn from the same population as the case participants and should be at risk to develop the disease and be detected in the study. Inclusion of women as controls in genetic association studies of diseases limited to men, for example, is problematic in that this approach adds individuals to the control group who had no chance of developing the disease (but might have done so had they also inherited a Y chromosome), thus mixing the controls with possible latent cases. This artificially reduces the differences in allele frequencies between cases and controls and limits the ability of the study to detect a true difference (ie, reduces study power).
If the disease is common, such as coronary heart disease or hypertension in the United States, efforts should be made to ensure that the controls are truly disease free. Some studies address this by using super-controls or persons at high risk but without even early evidence of disease, such as persons with diabetes of long duration but without microalbuminuria in a study of diabetic nephropathy.38 The success of recent GWA studies using control groups of questionable representativeness due to volunteer bias, such as the blood donor cohort in the Wellcome Trust Case-Control Consortium,39 suggests that initial identification of SNPs associated with disease may be robust to these biases, especially given subsequent evidence of replication of these associations in studies using more traditional control groups.40-42
Of more concern may be the risk of false-negative findings, as many biases tend to reduce the magnitude of observed associations toward the null. Use of convenience controls such as blood donors, however, may also be problematic in examining potential modification of genetic associations by environmental exposures and sociocultural factors, and in the identification of less strongly associated SNPs.
A key component in articles reporting results in the epidemiology literature of observational study is an initial table comparing relevant characteristics of those with and without disease, allowing assessment of comparability and generalizability of the 2 groups. Such comparisons are infrequent in GWA studies,28 but they are important because common diseases are typically influenced by multiple environmental (as well as genetic) factors. Important differences should be adjusted for in the analysis if possible, to avoid the risk of identifying genetic associations not with the disease of interest but with a confounding factor, such as smoking43 or obesity.44
Confounding due to population stratification (also called population structure) has been cited as a major threat to the validity of genetic association studies, but its true importance is a matter of debate.45-46 When variations occur in allele frequency between population subgroups, such as those defined by ethnicity or geographic origin, that in turn differ in their risk for disease, GWA studies may then falsely identify the subgroup-associated genes as related to disease.30 Population structure should be assessed and reported in GWA studies, typically by examining the distribution of test statistics generated from the thousands of association tests performed (eg, the  2 test) and assessing their deviation from the null distribution (that expected under the null hypothesis of no SNP associated with the trait) in a quantile-quantile or "Q-Q," plot (Figure 1). In these plots, observed association statistics or calculated P values for each SNP are ranked in order from smallest to largest and plotted against the values expected had they been sampled from a distribution of known form (such as the 2 distribution).39 Deviations from the diagonal identity line suggest that either the assumed distribution is incorrect or that the sample contains values arising in some other manner, as by a true association.39
Since the underlying assumption in GWA studies is that the vast majority of assayed SNPs are not associated with the trait, strong deviations from the null suggest either a very highly associated and heavily genotyped locus (Figure 1, A), or significant differences in population structure (Figure 1, B). Several effective statistical methods are available to correct for population structure and are a standard component of rigorous GWA analyses.28, 30
Genotyping and Quality Control in GWA Studies
GWA studies rely on the typically strong associations among SNPs located near each other on a chromosome, which tend to be inherited together more often than expected by chance.50 This nonrandom association is called linkage disequilibrium; alleles of SNPs in high linkage disequilibrium are almost always inherited together and can serve as proxies for each other. Their correlation with each other in the population is measured by the r2 statistic, which is the proportion of variation of one SNP explained by the other, and ranges from 0 (no association) to 1 (perfect correlation).
Genomic coverage of GWA genotyping platforms (arrays or chips on which genotyping is performed) is often estimated by the percent of common SNPs having an r2 of 0.8 or greater with at least 1 SNP on the platform.13 Genotyping platforms comprising 500 000 to 1 000 000 SNPs have been estimated to capture 67% to 89% of common SNP variation in populations of European and Asian ancestry and 46% to 66% of variation in individuals of recent African ancestry.13 Higher density platforms now also include probes for copy number variants that are not well tagged by SNPs. Copy number variants, in which stretches of genomic sequence are deleted or are duplicated in varying numbers, have gained increasing attention because of their apparent ubiquity and potential dosage effect on gene expression.51 Newer genotyping platforms are increasingly being focused on capturing copy number variants, but other structural variants such as insertions, deletions, and inversions, remain difficult to assay.52
GWA studies frequently identify associations with multiple SNPs in a chromosomal region and display the association statistics by their genomic location on a portion of a chromosome (Figure 2). For ease of display, association statistics are typically shown as the - log10 of the P value (the probability of the observed association arising by chance alone), so that P = .01 would be plotted as "2" on the y-axis and P = 10– 7 as "7." Such displays also often plot a matrix of r2 values for each pair of SNPs in the region, with larger r2 values more intensely shaded. These plots can be used to identify linkage disequilibrium blocks containing SNPs associated with disease, allowing estimation of the independence of the SNP associations observed.53
|
|
|
|
Figure 2. Associations in the IL23R Gene Region Identified by a Genome-wide Association Study of Inflammatory Bowel Disease
Genome-wide association studies frequently identify associations with many highly correlated single-nucleotide polymorphisms (SNPs) in a chromosomal region, due in part to linkage disequilibrium, among the SNPs. This can make it difficult to determine which SNP within a group is likely to be the causative or functional variant. A, Genomic locations of 2 genes, the interleukin 23 receptor (IL23R) and the interleukin 12 receptor, beta-2 (IL12Rb2), and a hypothetical protein, NM_001013674, between positions 62700000 and 67580000 of the short arm of chromosome 1 at region 1p31, are shown. B, The –log10 P values for association with inflammatory bowel disease are plotted for each SNP genotyped in the region; those reaching a prespecified value of –log10 of 7 or greater are presumed to show association with disease. Several strong associations, at –log10 P values or greater, are seen in the region just telomeric of position approximately 67400000 and extending just centromeric of position approximate 67450000. C, Pairwise linkage disequilibrium estimates between SNPs (measured as r2 ) are plotted for the region. Higher r2 values are indicated by darker shading. The region contains 4 "triangles" or "blocks" of linkage disequilibrium, 2 on either side of position 67400000 in the IL23R gene, another in the hypothetical protein telomeric of IL23R, and a fourth in the IL12RB2 gene at the centromeric end of the region. The 2 IL23R linkage disequilibrium regions each contain SNPs associated with inflammatory bowel disease, while the IL12RB2 region does not. Reproduced with permission from Duerr et al.53
|
|
|
Genotyping errors, especially if occurring differentially between cases and controls, are an important cause of spurious associations and must be diligently sought and corrected.54 A number of quality control features should be applied both on a per-sample and a per-SNP basis. Checks on sample identity to avoid sample mix-ups should be described and a minimum rate of successfully genotyped SNPs per sample (usually 80%-90% of SNPs attempted) should be reported. Once samples failing these thresholds are removed, individual SNPs across the remaining samples are subjected to further checks or filters for probable genotyping errors, including: (1) the proportion of samples for which a SNP can be measured (the SNP call rate, typically >95%); (2) the minor allele frequency (often >1%, as rarer SNPs are difficult to measure reliably; (3) severe violations of Hardy-Weinberg equilibrium; (4) Mendelian inheritance errors in trio studies; and (5) concordance rates in duplicate samples (typically >99.5%).
Additional checks on genotyping quality should include careful visual inspection of genotype cluster plots, or intensity values generated by the genotyping assay to ensure that the strongest associations do not merely reflect genotyping artifact.28, 39 Genotyping the most strongly associated SNPs should also be confirmed using a different method.28 Associations with any known "positive controls," such as TCF7L2 in type 2 diabetes mellitus55 or HLA-DRB1 in rheumatoid arthritis,47 should be reported to increase confidence in the consistency of findings with prior reports.
Analysis and Presentation of GWA Results
Associations with the 2 alleles of each SNP are tested in a relatively straightforward manner by comparing the frequency of each allele in cases and controls (Table 3). Because each individual carries 2 copies of each autosomal SNP, the frequency of each of 3 possible genotypes can also be compared (Table 3). Exploratory analyses may also include testing of different genetic models (dominant, recessive, or additive), although additive models, in which each copy of the allele is assumed to increase risk by the same amount, tend to be the most common (T.A.M.; unpublished data, 2008). Odds ratios of disease associated with the risk allele or genotype(s) can then be calculated and are typically modest, often in the range of 1.2 to 1.3. Many studies also calculate population attributable risk, classically defined as the proportion of disease in the population associated with a given risk factor (in this case, a genetic variant).57
|
|
|
|
Table 3. Association of Alleles and Genotypes of rs6983267 on Chromosome 8q24 With Colorectal Cancera
|
|
|
Such estimates are nearly always inflated because odds ratios overestimate relative risks (especially for common diseases58) needed for population attributable risk calculations, and because odds ratios and allele frequencies in published reports have wide confidence intervals so that those selected by exceeding a specified threshold for statistical significance tend to be biased upwards, an effect of ascertainment known as the "winner's curse."59 This exaggerated initial estimate of the odds ratio often leads to replication studies that lack sufficient sample size and power to replicate the association because larger samples are needed to detect smaller odds ratios.
Complexity in analysis emerges due to the multiple testing carried out in GWA studies, in that the association tests shown in Table 2 are repeated for each of the 100 000 to more than 1 million SNPs assayed (Figure 3). At the conventional P < .05 level of significance, an association study of 1 million SNPs will show 50 000 SNPs to be "associated" with disease, almost all falsely positive and due to chance alone. The most common manner of dealing with this problem is to reduce the false-positive rate by applying the Bonferroni correction, in which the conventional P value is divided by the number of tests performed.60 A 1 million SNP survey would thus use a threshold of P < .05/106, or 5 x 10– 8, to identify associations unlikely to have occurred by chance. This correction has been criticized as overly conservative because it assumes independent associations of each SNP with disease even though individual SNPs are known to be correlated to some degree due to linkage disequilibrium.
|
|
|
|
Figure 3. Genome-wide Association Findings in Rheumatoid Arthritis
Genome-wide association studies assume a priori hypotheses about candidate genes or regions that might be associated with disease; rather, they test single-nucleotide polymorphisms (SNPs) throughout the genome for possible evidence of genetic susceptibility. Associations plotted as –log10 P values for a genome-wide association study in 1522 cases with rheumatoid arthritis and 1850 controls, showing single data points for SNPs with P < 10–4 (lower horizontal red line) for 22 autosomes and the X chromosome. The predefined level of significance, at 5 x 10–8 is shown with a horizontal blue line. SNPs at PTPN22 on chromosome 1, the major histocompatibility comples (MHC) on chromosome 6, and the TRAF1-C5 locus on chromosome 9 exceed this threshold. Reproduced with permission from Plenge et al.47
|
|
|
Other approaches have been proposed, including estimation of the false discovery rate or proportion of significant associations that are actually false positive associations,61-62 false-positive report probability, or probability that the null hypothesis is true given a statistically significant finding,63 and estimation of Bayes factors that incorporate the prior probability of association based on characteristics of the disease or the specific SNP.39 To date, Bonferroni correction has generally been the most commonly used correction for multiple comparisons in GWA reports (T.A.M.; unpublished data, 2008).
Replication and Functional Studies
Given the major challenge of separating the many false-positive associations from the few true-positive associations with disease in GWA studies, an important strategy has been replication of results in independent samples.28 This is typically included in a single GWA report as part of a multistage design34-35 or may be reported separately.39, 64 Consensus criteria for replication have recently been published and include study of the same or very similar phenotype and population, and demonstration of a similar magnitude of effect and significance (in the same genetic model and same direction) for the same SNP and the same allele as the initial report.28 Replication is usually first attempted in studies as similar as possible to the initial report, but then may be extended to related phenotypes (such as fat mass in addition to obesity44), different populations (such as West Africans in addition to Icelanders65), or different study designs53 to refine and extend the initial findings and increase confidence in verity.
Lack of reproducibility of genetic associations has been frequently observed and has been varyingly attributed to population stratification, phenotype differences, selection biases, genotyping errors, and other factors.28, 66 At present, the best way of resolving these inconsistencies appears to be additional replication studies with larger sample sizes, although this may not be feasible for rare conditions or for associations identified in unique populations.28
Identification of a robustly replicating SNP-disease association is a crucial first step in identifying disease-causing genetic variants and developing suitable treatments, but it is only a first step. Association studies essentially identify a genomic location related to disease but provide little information on gene function unless SNPs with predictable effects on gene expression or the transcribed product happened to be identified. Few of the associations identified to date have involved genes previously suspected of being related to the disease under study, and some have been in genomic locations harboring no known genes.27, 67 Examination of known SNPs in high linkage disequilibrium with the associated SNP may identify variants with plausible biologic effects, or sequencing of a suitable surrounding interval may be undertaken to identify rarer variants with more obvious functional implications. Tissue samples or cell lines can be examined for expression of the gene variant. Other functional studies may include genetic manipulations in cell or animal models, such as knockouts or knock-ins.68
|
