Two common variants in the gene encoding complement factor H (CFH), the Y402H substitution (rs1061170, c.1204C>T)(1-4) and the intronic rs1410996 SNP(5,6), explain 17% of age-related macular degeneration (AMD) liability. However, proof for the involvement of CFH, as opposed to a neighboring transcript, and knowledge of the potential mechanism of susceptibility alleles are lacking. Assuming that rare functional variants might provide mechanistic insights, we used genotype data and high-throughput sequencing to discover a rare, high-risk CFH haplotype with a c.3628C>T mutation that resulted in an R1210C substitution. This allele has been implicated previously in atypical hemolytic uremic syndrome, and it abrogates C-terminal ligand binding(7,8). Genotyping R1210C in 2,423 AMD cases and 1,122 controls demonstrated high penetrance (present in 40 cases versus 1 control, P = 7.0 × 10(-6)) and an association with a 6-year-earlier onset of disease (P = 2.3 × 10(-6)). This result suggests that loss-of-function alleles at CFH are likely to drive AMD risk. This finding represents one of the first instances in which a common complex disease variant has led to the discovery of a rare penetrant mutation.

Although genome-wide association studies have implicated many individual loci in complex diseases, identifying the exact causal alleles and the cell types within which they act remains greatly challenging. To ultimately understand disease mechanism, researchers must carefully conceive functional studies in relevant pathogenic cell types to demonstrate the cellular impact of disease-associated genetic variants. This challenge is highlighted in autoimmune diseases, such as rheumatoid arthritis, where any of a broad range of immunological cell types might potentially be impacted by genetic variation to cause disease. To this end, we developed a statistical approach to identify potentially pathogenic cell types in autoimmune diseases by using a gene-expression data set of 223 murine-sorted immune cells from the Immunological Genome Consortium. We found enrichment of transitional B cell genes in systemic lupus erythematosus (p = 5.9 × 10(-6)) and epithelial-associated stimulated dendritic cell genes in Crohn disease (p = 1.6 × 10(-5)). Finally, we demonstrated enrichment of CD4+ effector memory T cell genes within rheumatoid arthritis loci (p < 10(-6)). To further validate the role of CD4+ effector memory T cells within rheumatoid arthritis, we identified 436 loci that were not yet known to be associated with the disease but that had a statistically suggestive association in a recent genome-wide association study (GWAS) meta-analysis (p(GWAS) < 0.001). Even among these putative loci, we noted a significant enrichment for genes specifically expressed in CD4+ effector memory T cells (p = 1.25 × 10(-4)). These cell types are primary candidates for future functional studies to reveal the role of risk alleles in autoimmunity. Our approach has application in other phenotypes, outside of autoimmunity, where many loci have been discovered and high-quality cell-type-specific gene expression is available.

Advances in genotyping and sequencing technologies have revolutionized the genetics of complex disease by locating rare and common variants that influence an individual's risk for diseases, such as diabetes, cancers, and psychiatric disorders. However, to capitalize on these data for prevention and therapies requires the identification of causal alleles and a mechanistic understanding for how these variants contribute to the disease. After discussing the strategies currently used to map variants for complex diseases, this Primer explores how variants may be prioritized for follow-up functional studies and the challenges and approaches for assessing the contributions of rare and common variants to disease phenotypes.

MOTIVATION: As disease loci are rapidly discovered, an emerging challenge is to identify common pathways and biological functionality across loci. Such pathways might point to potential disease mechanisms. One strategy is to look for functionally related or interacting genes across genetic loci. Previously, we defined a statistical strategy, Gene Relationships Across Implicated Loci (GRAIL), to identify whether pair-wise gene relationships defined using PubMed text similarity are enriched across loci. Here, we have implemented VIZ-GRAIL, a software tool to display those relationships and to depict the underlying biological patterns. RESULTS: Our tool can seamlessly interact with the GRAIL web site to obtain the results of analyses and create easy to read visual displays. To most clearly display results, VIZ-GRAIL arranges genes and genetic loci to minimize intersecting pair-wise gene connections. VIZ-GRAIL can be easily applied to other types of functional connections, beyond those from GRAIL. This method should help investigators appreciate the presence of potentially important common functions across loci. AVAILABILITY: The GRAIL algorithm is implemented online at http://www.broadinstitute.org/mpg/grail/grail.php. VIZ-GRAIL source-code is at http://www.broadinstitute.org/mpg/grail/vizgrail.html.

Investigators have linked rare copy number variation (CNVs) to neuropsychiatric diseases, such as schizophrenia. One hypothesis is that CNV events cause disease by affecting genes with specific brain functions. Under these circumstances, we expect that CNV events in cases should impact brain-function genes more frequently than those events in controls. Previous publications have applied "pathway" analyses to genes within neuropsychiatric case CNVs to show enrichment for brain-functions. While such analyses have been suggestive, they often have not rigorously compared the rates of CNVs impacting genes with brain function in cases to controls, and therefore do not address important confounders such as the large size of brain genes and overall differences in rates and sizes of CNVs. To demonstrate the potential impact of confounders, we genotyped rare CNV events in 2,415 unaffected controls with Affymetrix 6.0; we then applied standard pathway analyses using four sets of brain-function genes and observed an apparently highly significant enrichment for each set. The enrichment is simply driven by the large size of brain-function genes. Instead, we propose a case-control statistical test, cnv-enrichment-test, to compare the rate of CNVs impacting specific gene sets in cases versus controls. With simulations, we demonstrate that cnv-enrichment-test is robust to case-control differences in CNV size, CNV rate, and systematic differences in gene size. Finally, we apply cnv-enrichment-test to rare CNV events published by the International Schizophrenia Consortium (ISC). This approach reveals nominal evidence of case-association in neuronal-activity and the learning gene sets, but not the other two examined gene sets. The neuronal-activity genes have been associated in a separate set of schizophrenia cases and controls; however, testing in independent samples is necessary to definitively confirm this association. Our method is implemented in the PLINK software package.

Obesity is globally prevalent and highly heritable, but its underlying genetic factors remain largely elusive. To identify genetic loci for obesity susceptibility, we examined associations between body mass index and ∼ 2.8 million SNPs in up to 123,865 individuals with targeted follow up of 42 SNPs in up to 125,931 additional individuals. We confirmed 14 known obesity susceptibility loci and identified 18 new loci associated with body mass index (P < 5 × 10⁻⁸), one of which includes a copy number variant near GPRC5B. Some loci (at MC4R, POMC, SH2B1 and BDNF) map near key hypothalamic regulators of energy balance, and one of these loci is near GIPR, an incretin receptor. Furthermore, genes in other newly associated loci may provide new insights into human body weight regulation.

OBJECTIVE: A recent meta-analysis of published genome-wide association studies (GWAS) in populations of European descent reported novel associations of markers mapping to the CD40, CCL21 and CDK6 genes with rheumatoid arthritis (RA) susceptibility while a large-scale, case-control association study in a Japanese population identified association with multiple single nucleotide polymorphisms (SNPs) in the CD244 gene. The aim of the current study was to validate these potential RA susceptibility markers in a UK population. METHODS: A total of 4 SNPs (rs4810485 in CD40, rs2812378 in CCL21, rs42041 in CDK6 and rs6682654 in CD244) were genotyped in a UK cohort comprising 3962 UK patients with RA and 3531 healthy controls using the Sequenom iPlex platform. Genotype counts in patients and controls were analysed with the chi(2) test using Stata. RESULTS: Association to the CD40 gene was robustly replicated (p=2 x 10(-4), OR 0.86, 95% CI 0.79 to 0.93) and modest evidence was found for association with the CCL21 locus (p=0.04, OR 1.08, 95% CI 1.01 to 1.16). However, there was no evidence for association of rs42041 (CDK6) and rs6682654 (CD244) with RA susceptibility in this UK population. Following a meta-analysis including the original data, association to CD40 was confirmed (p=7.8 x 10(-8), OR 0.87 (95% CI 0.83 to 0.92). CONCLUSION: In this large UK cohort, strong association of the CD40 gene with susceptibility to RA was found, and weaker evidence for association with RA in the CCL21 locus.

BACKGROUND: Recent discoveries of risk alleles have made it possible to define genetic risk profiles for patients with rheumatoid arthritis (RA). This study examined whether a cumulative score based on 22 validated genetic risk alleles for seropositive RA would identify high-risk, asymptomatic individuals who might benefit from preventive interventions. METHODS: Eight human leucocyte antigen (HLA) alleles and 14 single-nucleotide polymorphisms representing 13 validated RA risk loci were genotyped among 289 white seropositive cases and 481 controls from the US Nurses' Health Studies (NHS) and 629 white cyclic-citrullinated peptide antibody-positive cases and 623 controls from the Swedish Epidemiologic Investigation of Rheumatoid Arthritis (EIRA). A weighted genetic risk score (GRS) was created, in which the weight for each risk allele is the log of the published odds ratio (OR). Logistic regression was used to study associations with incident RA. Area under the curve (AUC) statistics were compared from a clinical-only model and clinical plus genetic model in each cohort. RESULTS: Patients with GRS >1.25 SD of the mean had a significantly higher OR of seropositive RA in both NHS (OR=2.9, 95%CI 1.8 to 4.6) and EIRA (OR 3.4, 95% CI 2.3 to 5.0) referent to the population average. In NHS, the AUC for a clinical model was 0.57 and for a clinical plus genetic model was 0.66, and in EIRA was 0.63 and 0.75, respectively. CONCLUSION: The combination of 22 risk alleles into a weighted GRS significantly stratifies individuals for RA risk beyond clinical risk factors alone. Given the low incidence of RA, the clinical utility of a weighted GRS is limited in the general population.

Advanced age-related macular degeneration (AMD) is the leading cause of late onset blindness. We present results of a genome-wide association study of 979 advanced AMD cases and 1,709 controls using the Affymetrix 6.0 platform with replication in seven additional cohorts (totaling 5,789 unrelated cases and 4,234 unrelated controls). We also present a comprehensive analysis of copy-number variations and polymorphisms for AMD. Our discovery data implicated the association between AMD and a variant in the hepatic lipase gene (LIPC) in the high-density lipoprotein cholesterol (HDL) pathway (discovery P = 4.53e-05 for rs493258). Our LIPC association was strongest for a functional promoter variant, rs10468017, (P = 1.34e-08), that influences LIPC expression and serum HDL levels with a protective effect of the minor T allele (HDL increasing) for advanced wet and dry AMD. The association we found with LIPC was corroborated by the Michigan/Penn/Mayo genome-wide association study; the locus near the tissue inhibitor of metalloproteinase 3 was corroborated by our replication cohort for rs9621532 with P = 3.71e-09. We observed weaker associations with other HDL loci (ABCA1, P = 9.73e-04; cholesterylester transfer protein, P = 1.41e-03; FADS1-3, P = 2.69e-02). Based on a lack of consistent association between HDL increasing alleles and AMD risk, the LIPC association may not be the result of an effect on HDL levels, but it could represent a pleiotropic effect of the same functional component. Results implicate different biologic pathways than previously reported and provide new avenues for prevention and treatment of AMD.

We undertook a meta-analysis of six Crohn's disease genome-wide association studies (GWAS) comprising 6,333 affected individuals (cases) and 15,056 controls and followed up the top association signals in 15,694 cases, 14,026 controls and 414 parent-offspring trios. We identified 30 new susceptibility loci meeting genome-wide significance (P < 5 × 10⁻⁸). A series of in silico analyses highlighted particular genes within these loci and, together with manual curation, implicated functionally interesting candidate genes including SMAD3, ERAP2, IL10, IL2RA, TYK2, FUT2, DNMT3A, DENND1B, BACH2 and TAGAP. Combined with previously confirmed loci, these results identify 71 distinct loci with genome-wide significant evidence for association with Crohn's disease.

Most common human traits and diseases have a polygenic pattern of inheritance: DNA sequence variants at many genetic loci influence the phenotype. Genome-wide association (GWA) studies have identified more than 600 variants associated with human traits, but these typically explain small fractions of phenotypic variation, raising questions about the use of further studies. Here, using 183,727 individuals, we show that hundreds of genetic variants, in at least 180 loci, influence adult height, a highly heritable and classic polygenic trait. The large number of loci reveals patterns with important implications for genetic studies of common human diseases and traits. First, the 180 loci are not random, but instead are enriched for genes that are connected in biological pathways (P = 0.016) and that underlie skeletal growth defects (P < 0.001). Second, the likely causal gene is often located near the most strongly associated variant: in 13 of 21 loci containing a known skeletal growth gene, that gene was closest to the associated variant. Third, at least 19 loci have multiple independently associated variants, suggesting that allelic heterogeneity is a frequent feature of polygenic traits, that comprehensive explorations of already-discovered loci should discover additional variants and that an appreciable fraction of associated loci may have been identified. Fourth, associated variants are enriched for likely functional effects on genes, being over-represented among variants that alter amino-acid structure of proteins and expression levels of nearby genes. Our data explain approximately 10% of the phenotypic variation in height, and we estimate that unidentified common variants of similar effect sizes would increase this figure to approximately 16% of phenotypic variation (approximately 20% of heritable variation). Although additional approaches are needed to dissect the genetic architecture of polygenic human traits fully, our findings indicate that GWA studies can identify large numbers of loci that implicate biologically relevant genes and pathways.

A powerful way to discover key genes with causal roles in oncogenesis is to identify genomic regions that undergo frequent alteration in human cancers. Here we present high-resolution analyses of somatic copy-number alterations (SCNAs) from 3,131 cancer specimens, belonging largely to 26 histological types. We identify 158 regions of focal SCNA that are altered at significant frequency across several cancer types, of which 122 cannot be explained by the presence of a known cancer target gene located within these regions. Several gene families are enriched among these regions of focal SCNA, including the BCL2 family of apoptosis regulators and the NF-kappaBeta pathway. We show that cancer cells containing amplifications surrounding the MCL1 and BCL2L1 anti-apoptotic genes depend on the expression of these genes for survival. Finally, we demonstrate that a large majority of SCNAs identified in individual cancer types are present in several cancer types.

By combining genome-wide association data from 8,130 individuals with type 2 diabetes (T2D) and 38,987 controls of European descent and following up previously unidentified meta-analysis signals in a further 34,412 cases and 59,925 controls, we identified 12 new T2D association signals with combined P<5x10(-8). These include a second independent signal at the KCNQ1 locus; the first report, to our knowledge, of an X-chromosomal association (near DUSP9); and a further instance of overlap between loci implicated in monogenic and multifactorial forms of diabetes (at HNF1A). The identified loci affect both beta-cell function and insulin action, and, overall, T2D association signals show evidence of enrichment for genes involved in cell cycle regulation. We also show that a high proportion of T2D susceptibility loci harbor independent association signals influencing apparently unrelated complex traits.

OBJECTIVE: Electronic medical records (EMRs) are a rich data source for discovery research but are underutilized due to the difficulty of extracting highly accurate clinical data. We assessed whether a classification algorithm incorporating narrative EMR data (typed physician notes) more accurately classifies subjects with rheumatoid arthritis (RA) compared with an algorithm using codified EMR data alone. METHODS: Subjects with > or =1 International Classification of Diseases, Ninth Revision RA code (714.xx) or who had anti-cyclic citrullinated peptide (anti-CCP) checked in the EMR of 2 large academic centers were included in an "RA Mart" (n = 29,432). For all 29,432 subjects, we extracted narrative (using natural language processing) and codified RA clinical information. In a training set of 96 RA and 404 non-RA cases from the RA Mart classified by medical record review, we used narrative and codified data to develop classification algorithms using logistic regression. These algorithms were applied to the entire RA Mart. We calculated and compared the positive predictive value (PPV) of these algorithms by reviewing the records of an additional 400 subjects classified as having RA by the algorithms. RESULTS: A complete algorithm (narrative and codified data) classified RA subjects with a significantly higher PPV of 94% than an algorithm with codified data alone (PPV of 88%). Characteristics of the RA cohort identified by the complete algorithm were comparable to existing RA cohorts (80% women, 63% anti-CCP positive, and 59% positive for erosions). CONCLUSION: We demonstrate the ability to utilize complete EMR data to define an RA cohort with a PPV of 94%, which was superior to an algorithm using codified data alone.

To identify new genetic risk factors for rheumatoid arthritis, we conducted a genome-wide association study meta-analysis of 5,539 autoantibody-positive individuals with rheumatoid arthritis (cases) and 20,169 controls of European descent, followed by replication in an independent set of 6,768 rheumatoid arthritis cases and 8,806 controls. Of 34 SNPs selected for replication, 7 new rheumatoid arthritis risk alleles were identified at genome-wide significance (P < 5 x 10(-8)) in an analysis of all 41,282 samples. The associated SNPs are near genes of known immune function, including IL6ST, SPRED2, RBPJ, CCR6, IRF5 and PXK. We also refined associations at two established rheumatoid arthritis risk loci (IL2RA and CCL21) and confirmed the association at AFF3. These new associations bring the total number of confirmed rheumatoid arthritis risk loci to 31 among individuals of European ancestry. An additional 11 SNPs replicated at P < 0.05, many of which are validated autoimmune risk alleles, suggesting that most represent genuine rheumatoid arthritis risk alleles.

OBJECTIVE: Anti-tumor necrosis factor alpha (anti-TNF) therapy is a mainstay of treatment in rheumatoid arthritis (RA). The aim of the present study was to test established RA genetic risk factors to determine whether the same alleles also influence the response to anti-TNF therapy. METHODS: A total of 1,283 RA patients receiving etanercept, infliximab, or adalimumab therapy were studied from among an international collaborative consortium of 9 different RA cohorts. The primary end point compared RA patients with a good treatment response according to the European League Against Rheumatism (EULAR) response criteria (n = 505) with RA patients considered to be nonresponders (n = 316). The secondary end point was the change from baseline in the level of disease activity according to the Disease Activity Score in 28 joints (triangle upDAS28). Clinical factors such as age, sex, and concomitant medications were tested as possible correlates of treatment response. Thirty-one single-nucleotide polymorphisms (SNPs) associated with the risk of RA were genotyped and tested for any association with treatment response, using univariate and multivariate logistic regression models. RESULTS: Of the 31 RA-associated risk alleles, a SNP at the PTPRC (also known as CD45) gene locus (rs10919563) was associated with the primary end point, a EULAR good response versus no response (odds ratio [OR] 0.55, P = 0.0001 in the multivariate model). Similar results were obtained using the secondary end point, the triangle upDAS28 (P = 0.0002). There was suggestive evidence of a stronger association in autoantibody-positive patients with RA (OR 0.55, 95% confidence interval [95% CI] 0.39-0.76) as compared with autoantibody-negative patients (OR 0.90, 95% CI 0.41-1.99). CONCLUSION: Statistically significant associations were observed between the response to anti-TNF therapy and an RA risk allele at the PTPRC gene locus. Additional studies will be required to replicate this finding in additional patient collections.

PURPOSE OF REVIEW: To review the recently discovered genetic risk loci in rheumatoid arthritis (RA), the pathways they implicate, and the genetic architecture of RA. RECENT FINDINGS: Since 2008 investigators have identified many common genetic variants that confer disease risk through single nucleotide polymorphism genotyping studies; the list of variants will no doubt continue to expand at a rapid rate as genotyping technologies evolve and case-control sample collections continue to grow. In aggregate, these variants implicate pathways leading to NF-kappaB (nuclear factor kappa-light-chain-enhancer of activated B cells) activation, the interluekin-2 signaling pathway, and T-cell activation. SUMMARY: Although the effect of any individual variant is modest and even in aggregate considerably less than that of the major histocompatability complex, discovery of recent risk variants suggests immunological processes that are involved in disease pathogenesis.

The purpose of this article is to place these genetic discoveries in the context of current and future therapeutic strategies for patients with RA. More specifically, this article focuses on (1) a brief overview of genetic studies, (2) human genetics as an approach to identify the Achilles heel of disease pathways, (3) humans as the model organism for functional studies of human mutations, (4) pharmacogenetic studies to gain insight into the mechanism of action of drugs, and (5) next-generation patient registries to enable large-scale genotype-phenotype studies.

OBJECTIVE: The severity of joint destruction in rheumatoid arthritis (RA) is highly variable from patient to patient and is influenced by genetic factors. Genome-wide association studies have enormously boosted the field of the genetics of RA susceptibility, but risk loci for RA severity remain poorly defined. A recent meta-analysis of genome-wide association studies identified 6 genetic regions for susceptibility to autoantibody-positive RA: CD40, KIF5A/PIP4K2C, CDK6, CCL21, PRKCQ, and MMEL1/TNFRSF14. The purpose of this study was to investigate whether these newly described genetic regions are associated with the rate of joint destruction. METHODS: RA patients enrolled in the Leiden Early Arthritis Clinic were studied (n=563). Yearly radiographs were scored using the Sharp/van der Heijde method (median followup 5 years; maximum followup 9 years). The rate of joint destruction between genotype groups was compared using a linear mixed model, correcting for age, sex, and treatment strategies. A total of 393 anti-citrullinated protein antibody (ACPA)-positive RA patients from the North American Rheumatoid Arthritis Consortium (NARAC) who had radiographic data available were used for the replication study. RESULTS: The TT and CC/CG genotypes of 2 single-nucleotide polymorphisms, rs4810485 (CD40) and rs42041 (CDK6), respectively, were associated with a higher rate of joint destruction in ACPA-positive RA patients (P=0.003 and P=0.012, respectively), with rs4810485 being significant after Bonferroni correction for multiple testing. The association of the CD40 minor allele with the rate of radiographic progression was replicated in the NARAC cohort (P=0.021). CONCLUSION: A polymorphism in the CD40 locus is associated with the rate of joint destruction in patients with ACPA-positive RA. Our findings provide one of the first non-HLA-related genetic severity factors that has been replicated.