Genetic Susceptibility to Postdiarrheal Hemolytic-Uremic Syndrome After Shiga Toxin-Producing Escherichia coli Infection: A Centers for Disease Control and Prevention FoodNet Study.

Background: Postdiarrheal hemolytic-uremic syndrome (D+HUS) following Shiga toxin-producing Escherichia coli (STEC) infection is a serious condition lacking specific treatment. Host immune dysregulation and genetic susceptibility to complement hyperactivation are implicated in non-STEC-related HUS. However, genetic susceptibility to D+HUS remains largely uncharacterized. Methods: Patients with culture-confirmed STEC diarrhea, identified through the Centers for Disease Control and Prevention FoodNet surveillance system (2007-2012), were serotyped and classified by laboratory and/or clinical criteria as having suspected, probable, or confirmed D+HUS or as controls and underwent genotyping at 200 loci linked to nondiarrheal HUS or similar pathologies. Genetic associations with D+HUS were explored by multivariable regression, with adjustment for known risk factors. Results: Of 641 enrollees with STEC O157:H7, 80 had suspected D+HUS (41 with probable and 32 with confirmed D+HUS). Twelve genes related to cytokine signaling, complement pathways, platelet function, pathogen recognition, iron transport, and endothelial function were associated with D+HUS in multivariable-adjusted analyses (P ≤ .05). Of 12 significant single-nucleotide polymorphisms (SNPs), 5 were associated with all levels of D+HUS (intergenic SNP rs10874639, TFRC rs3804141, EDN1 rs5370, GP1BA rs121908064, and B2M rs16966334), and 7 SNPs (6 non-complement related) were associated with confirmed D+HUS (all P < .05). Conclusions: Polymorphisms in many non-complement-related genes may contribute to D+HUS susceptibility. These results require replication, but they suggest novel therapeutic targets in patients with D+HUS.

Rare and common variants in extracellular matrix gene Fibrillin 2 (FBN2) are associated with macular degeneration.

Neurodegenerative diseases affecting the macula constitute a major cause of incurable vision loss and exhibit considerable clinical and genetic heterogeneity, from early-onset monogenic disease to multifactorial late-onset age-related macular degeneration (AMD). As part of our continued efforts to define genetic causes of macular degeneration, we performed whole exome sequencing in four individuals of a two-generation family with autosomal dominant maculopathy and identified a rare variant p.Glu1144Lys in Fibrillin 2 (FBN2), a glycoprotein of the elastin-rich extracellular matrix (ECM). Sanger sequencing validated the segregation of this variant in the complete pedigree, including two additional affected and one unaffected individual. Sequencing of 192 maculopathy patients revealed additional rare variants, predicted to disrupt FBN2 function. We then undertook additional studies to explore the relationship of FBN2 to macular disease. We show that FBN2 localizes to Bruch's membrane and its expression appears to be reduced in aging and AMD eyes, prompting us to examine its relationship with AMD. We detect suggestive association of a common FBN2 non-synonymous variant, rs154001 (p.Val965Ile) with AMD in 10 337 cases and 11 174 controls (OR = 1.10; P-value = 3.79 × 10(-5)). Thus, it appears that rare and common variants in a single gene--FBN2--can contribute to Mendelian and complex forms of macular degeneration. Our studies provide genetic evidence for a key role of elastin microfibers and Bruch's membrane in maintaining blood-retina homeostasis and establish the importance of studying orphan diseases for understanding more common clinical phenotypes.

Rare variants in RTEL1 are associated with familial interstitial pneumonia.

RATIONALE: Up to 20% of cases of idiopathic interstitial pneumonia cluster in families, comprising the syndrome of familial interstitial pneumonia (FIP); however, the genetic basis of FIP remains uncertain in most families. OBJECTIVES: To determine if new disease-causing rare genetic variants could be identified using whole-exome sequencing of affected members from FIP families, providing additional insights into disease pathogenesis. METHODS: Affected subjects from 25 kindreds were selected from an ongoing FIP registry for whole-exome sequencing from genomic DNA. Candidate rare variants were confirmed by Sanger sequencing, and cosegregation analysis was performed in families, followed by additional sequencing of affected individuals from another 163 kindreds. MEASUREMENTS AND MAIN RESULTS: We identified a potentially damaging rare variant in the gene encoding for regulator of telomere elongation helicase 1 (RTEL1) that segregated with disease and was associated with very short telomeres in peripheral blood mononuclear cells in 1 of 25 families in our original whole-exome sequencing cohort. Evaluation of affected individuals in 163 additional kindreds revealed another eight families (4.7%) with heterozygous rare variants in RTEL1 that segregated with clinical FIP. Probands and unaffected carriers of these rare variants had short telomeres (<10% for age) in peripheral blood mononuclear cells and increased T-circle formation, suggesting impaired RTEL1 function. CONCLUSIONS: Rare loss-of-function variants in RTEL1 represent a newly defined genetic predisposition for FIP, supporting the importance of telomere-related pathways in pulmonary fibrosis.

Common variants at 9p21 and 8q22 are associated with increased susceptibility to optic nerve degeneration in glaucoma.

Optic nerve degeneration caused by glaucoma is a leading cause of blindness worldwide. Patients affected by the normal-pressure form of glaucoma are more likely to harbor risk alleles for glaucoma-related optic nerve disease. We have performed a meta-analysis of two independent genome-wide association studies for primary open angle glaucoma (POAG) followed by a normal-pressure glaucoma (NPG, defined by intraocular pressure (IOP) less than 22 mmHg) subgroup analysis. The single-nucleotide polymorphisms that showed the most significant associations were tested for association with a second form of glaucoma, exfoliation-syndrome glaucoma. The overall meta-analysis of the GLAUGEN and NEIGHBOR dataset results (3,146 cases and 3,487 controls) identified significant associations between two loci and POAG: the CDKN2BAS region on 9p21 (rs2157719 [G], OR = 0.69 [95%CI 0.63-0.75], p = 1.86×10⁻¹8), and the SIX1/SIX6 region on chromosome 14q23 (rs10483727 [A], OR = 1.32 [95%CI 1.21-1.43], p = 3.87×10⁻¹¹). In sub-group analysis two loci were significantly associated with NPG: 9p21 containing the CDKN2BAS gene (rs2157719 [G], OR = 0.58 [95% CI 0.50-0.67], p = 1.17×10⁻¹²) and a probable regulatory region on 8q22 (rs284489 [G], OR = 0.62 [95% CI 0.53-0.72], p = 8.88×10⁻¹°). Both NPG loci were also nominally associated with a second type of glaucoma, exfoliation syndrome glaucoma (rs2157719 [G], OR = 0.59 [95% CI 0.41-0.87], p = 0.004 and rs284489 [G], OR = 0.76 [95% CI 0.54-1.06], p = 0.021), suggesting that these loci might contribute more generally to optic nerve degeneration in glaucoma. Because both loci influence transforming growth factor beta (TGF-beta) signaling, we performed a genomic pathway analysis that showed an association between the TGF-beta pathway and NPG (permuted p = 0.009). These results suggest that neuro-protective therapies targeting TGF-beta signaling could be effective for multiple forms of glaucoma.

Mapping of the disease locus and identification of ADAMTS10 as a candidate gene in a canine model of primary open angle glaucoma.

Primary open angle glaucoma (POAG) is a leading cause of blindness worldwide, with elevated intraocular pressure as an important risk factor. Increased resistance to outflow of aqueous humor through the trabecular meshwork causes elevated intraocular pressure, but the specific mechanisms are unknown. In this study, we used genome-wide SNP arrays to map the disease gene in a colony of Beagle dogs with inherited POAG to within a single 4 Mb locus on canine chromosome 20. The Beagle POAG locus is syntenic to a previously mapped human quantitative trait locus for intraocular pressure on human chromosome 19. Sequence capture and next-generation sequencing of the entire canine POAG locus revealed a total of 2,692 SNPs segregating with disease. Of the disease-segregating SNPs, 54 were within exons, 8 of which result in amino acid substitutions. The strongest candidate variant causes a glycine to arginine substitution in a highly conserved region of the metalloproteinase ADAMTS10. Western blotting revealed ADAMTS10 protein is preferentially expressed in the trabecular meshwork, supporting an effect of the variant specific to aqueous humor outflow. The Gly661Arg variant in ADAMTS10 found in the POAG Beagles suggests that altered processing of extracellular matrix and/or defects in microfibril structure or function may be involved in raising intraocular pressure, offering specific biochemical targets for future research and treatment strategies.

Common variants near CAV1 and CAV2 are associated with primary open-angle glaucoma in Caucasians from the USA.

Primary open-angle glaucoma (POAG) is a genetically complex common disease characterized by progressive optic nerve degeneration that results in irreversible blindness. Recently, a genome-wide association study (GWAS) for POAG in an Icelandic population identified significant associations with single nucleotide polymorphisms (SNPs) between the CAV1 and CAV2 genes on chromosome 7q31. In this study, we confirm that the identified SNPs are associated with POAG in our Caucasian US population and that specific haplotypes located in the CAV1/CAV2 intergenic region are associated with the disease. We also present data suggesting that associations with several CAV1/CAV2 SNPs are significant mostly in women.

Pitfalls of merging GWAS data: lessons learned in the eMERGE network and quality control procedures to maintain high data quality.

Genome-wide association studies (GWAS) are a useful approach in the study of the genetic components of complex phenotypes. Aside from large cohorts, GWAS have generally been limited to the study of one or a few diseases or traits. The emergence of biobanks linked to electronic medical records (EMRs) allows the efficient reuse of genetic data to yield meaningful genotype-phenotype associations for multiple phenotypes or traits. Phase I of the electronic MEdical Records and GEnomics (eMERGE-I) Network is a National Human Genome Research Institute-supported consortium composed of five sites to perform various genetic association studies using DNA repositories and EMR systems. Each eMERGE site has developed EMR-based algorithms to comprise a core set of 14 phenotypes for extraction of study samples from each site's DNA repository. Each eMERGE site selected samples for a specific phenotype, and these samples were genotyped at either the Broad Institute or at the Center for Inherited Disease Research using the Illumina Infinium BeadChip technology. In all, approximately 17,000 samples from across the five sites were genotyped. A unified quality control (QC) pipeline was developed by the eMERGE Genomics Working Group and used to ensure thorough cleaning of the data. This process includes examination of sample and marker quality and various batch effects. Upon completion of the genotyping and QC analyses for each site's primary study, eMERGE Coordinating Center merged the datasets from all five sites. This larger merged dataset reentered the established eMERGE QC pipeline. Based on lessons learned during the process, additional analyses and QC checkpoints were added to the pipeline to ensure proper merging. Here, we explore the challenges associated with combining datasets from different genotyping centers and describe the expansion to eMERGE QC pipeline for merged datasets. These additional steps will be useful as the eMERGE project expands to include additional sites in eMERGE-II, and also serve as a starting point for investigators merging multiple genotype datasets accessible through the National Center for Biotechnology Information in the database of Genotypes and Phenotypes. Our experience demonstrates that merging multiple datasets after additional QC can be an efficient use of genotype data despite new challenges that appear in the process.

Deletion of CFHR3 and CFHR1 genes in age-related macular degeneration.

Age-related macular degeneration (AMD) impairs vision for approximately 7.5 million Americans. Both susceptibility variants and protective haplotypes in the complement factor H (CFH) gene modulate risk for AMD. Recently, deletion of the 'CFH-related' genes CFHR1 and CFHR3 was found to be segregating with a particular CFH haplotype, which reduced the risk of AMD. We tested the deletion for association in a Caucasian population of 780 cases and 265 controls and examined its effect in the context of known AMD risk factors. The deletion did not segregate perfectly with any one SNP, as previously suggested. CFH haplotype P2 was the most frequent haplotype in deletion homozygotes (47%), and the majority (14/16) of these individuals were homozygous for the non-risk allele of Y402H. Overall, deletion homozygosity was significantly more frequent in controls than cases (2.6% controls, 0.8% cases, P = 0.025, OR = 0.29, 95% CI = 0.10-0.86). After controlling for age, Y402H, smoking and LOC387715 A69S, the protective effect of the deletion was no longer statistically significant (P = 0.27). However, using a CFH haplotype that all deletion homozygotes share as a surrogate for the deletion, this marker remained modestly associated with AMD after adjustment for known risk factors (OR = 0.63, 95% CI 0.39-1.04, P = 0.07). Therefore, deletion of CFHR1 and CFHR3 may account for a small portion of the protection from AMD associated with particular haplotypes in CFH. The presence of protective haplotypes in CFH that do not carry the deletion, suggests that other protective variants in this region have yet to be discovered.

C3 R102G polymorphism increases risk of age-related macular degeneration.

Inflammation has long been suspected to play a role in the pathogenesis of age-related macular degeneration (AMD). Association of variants in the complement factor H (CFH) and complement factor B (CFB) genes has targeted the search for additional loci to the alternative complement cascade, of which C3 is a major component. Two non-synonymous coding polymorphisms within C3, R102G and L314P, have previously been strongly associated with increased risk. These variants are in strong linkage disequilibrium (LD), making the contribution of this locus to AMD even more difficult to ascertain. We sought to determine whether the C3 association resulted primarily from only one of these two variants or from a combined effect of both in 223 families and an independent dataset of 701 cases and 286 unrelated controls. The C3 polymorphisms were in strong LD (r(2) = 0.85), and both were associated in the family-based and case-control datasets (R102G genoPDT P = 0.02, case-control genotypic P = 0.004; L314P genoPDT P = 0.001, case-control genotypic P = 0.04). In conditional analyses in the case-control dataset, R102G remained associated with disease in the L314P risk allele carriers (P = 0.01), but there was no effect of L314P in the R102G risk allele carriers (P = 0.2). After adjusting for age, smoking, CFH Y402H, LOC387715 A69S, and CFB R32Q, the effect of R102G remained strong [P = 0.015, odds ratio = 1.55, 95% confidence interval 1.09 to 2.21, adjusted PAR(population attributable risk) = 0.17]. Therefore, while the strong LD between R102G and L314P makes it difficult to disentangle their individual effects on disease risk, the R102G polymorphism acting alone provides the best model for disease in our data.

Haplotypes spanning the complement factor H gene are protective against age-related macular degeneration.

PURPOSE: Age-related macular degeneration (AMD) is a devastating disorder that adversely affects the quality of life of nearly 2 million Americans who have advanced forms of the disease. Besides the well-known risk imparted by carrying the Y402H variant in the complement factor H (CFH) gene on chromosome 1, recent evidence of the existence of protective haplotypes spanning CFH has been reported. METHODS: The haplo.stats program was used to test for association of the protective haplotypes after adjusting for age in the dataset of 584 sporadic cases and 248 control samples. Logistic regression modeling and likelihood ratio tests were used to investigate an interaction between a particular haplotype and smoking status. The HBAT option of FBAT was used to confirm the associations in an independent dataset of 201 families. RESULTS: Two protective (P) haplotypes in a family-based dataset (P1 = CAATTTAG, P = 0.021; and P2 = CGGCTTAG, P = 0.018) were identified for the first time. Age-adjusted score statistics provided support for these protective haplotypes in the case-control dataset (P1 frequency in cases approximately 13%, in controls approximately 20%, P = 0.001; P2 frequency in cases approximately 5%, in controls approximately 8%, P = 0.077). There was also tentative evidence of an interaction between one of the protective haplotypes and cigarette smoking (P = 0.04 likelihood ratio test for P2-smoking interaction). CONCLUSIONS: Replication of the association between the protective haplotypes and decreased AMD susceptibility provides increased evidence that these associations have biological meaning. The suggestion of a haplotype-smoking interaction adds to the growing body of evidence that smoking is an important environmental covariate in AMD that should be considered in genetic studies. Identification of the protective variant(s) carried within these haplotypes is critical for understanding the etiology of AMD.

Genome-wide and Ordered-Subset linkage analyses provide support for autism loci on 17q and 19p with evidence of phenotypic and interlocus genetic correlates.

BACKGROUND: Autism is a neurobehavioral spectrum of phenotypes characterized by deficits in the development of language and social relationships and patterns of repetitive, rigid and compulsive behaviors. Twin and family studies point to a significant genetic etiology, and several groups have performed genomic linkage screens to identify susceptibility loci. METHODS: We performed a genome-wide linkage screen in 158 combined Tufts, Vanderbilt and AGRE (Autism Genetics Research Exchange) multiplex autism families using parametric and nonparametric methods with a categorical autism diagnosis to identify loci of main effect. Hypothesizing interdependence of genetic risk factors prompted us to perform exploratory studies applying the Ordered-Subset Analysis (OSA) approach using LOD scores as the trait covariate for ranking families. We employed OSA to test for interlocus correlations between loci with LOD scores > or =1.5, and empirically determined significance of linkage in optimal OSA subsets using permutation testing. Exploring phenotypic correlates as the basis for linkage increases involved comparison of mean scores for quantitative trait-based subsets of autism between optimal subsets and the remaining families. RESULTS: A genome-wide screen for autism loci identified the best evidence for linkage to 17q11.2 and 19p13, with maximum multipoint heterogeneity LOD scores of 2.9 and 2.6, respectively. Suggestive linkage (LOD scores > or =1.5) at other loci included 3p, 6q, 7q, 12p, and 16p. OSA revealed positive correlations of linkage between the 19p locus and 17q, between 19p and 6q, and between 7q and 5p. While potential phenotypic correlates for these findings were not identified for the chromosome 7/5 combination, differences indicating more rapid achievement of "developmental milestones" was apparent in the chromosome 19 OSA-defined subsets for 17q and 6q. OSA was used to test the hypothesis that 19p linkage involved more rapid achievement of these milestones and it revealed significantly increased LOD* scores at 19p13. CONCLUSIONS: Our results further support 19p13 as harboring an autism susceptibility locus, confirm other linkage findings at 17q11.2, and demonstrate the need to analyze more discreet trait-based subsets of complex phenotypes to improve ability to detect genetic effects.

Screening for UGT1A1 Genotype in Study A5257 Would Have Markedly Reduced Premature Discontinuation of Atazanavir for Hyperbilirubinemia.

Background. Some patients are not prescribed atazanavir because of concern about possible jaundice. Atazanavir-associated hyperbilirubinemia correlates with UGT1A1 rs887829 genotype. We examined bilirubin-related discontinuation of atazanavir in participants from AIDS Clinical Trials Group Study A5257. Methods. Discriminatory properties of UGT1A1 T/T genotype for predicting bilirubin-related atazanavir discontinuation through 96 weeks after antiretroviral initiation were estimated. Results. Genetic analyses involved 1450 participants, including 481 who initiated randomized atazanavir/ritonavir. Positive predictive values of rs887829 T/T for bilirubin-related discontinuation of atazanavir (with 95% confidence intervals [CIs]) were 20% (CI, 9%-36%) in Black, 60% (CI, 32%-84%) in White, and 29% (CI, 8%-58%) in Hispanic participants; negative predictive values were 97% (CI, 93%-99%), 95% (CI, 90%-98%), and 97% (CI, 90%-100%), respectively. Conclusions. Bilirubin-related discontinuation of atazanavir was rare in participants not homozygous for rs887829 T/T, regardless of race or ethnicity. We hypothesize that the higher rate of discontinuation among White participants homozygous for rs887829 T/T may reflect differences in physical manifestations of jaundice by race and ethnicity. Selective avoidance of atazanavir initiation among individuals with T/T genotypes would markedly reduce the likelihood of bilirubin-related discontinuation of atazanavir while allowing atazanavir to be prescribed to the majority of individuals. This genetic association will also affect atazanavir/cobicistat.

Genome-Wide Association Study of Serum Creatinine Levels during Vancomycin Therapy.

Vancomycin, a commonly used antibiotic, can be nephrotoxic. Known risk factors such as age, creatinine clearance, vancomycin dose / dosing interval, and concurrent nephrotoxic medications fail to accurately predict nephrotoxicity. To identify potential genomic risk factors, we performed a genome-wide association study (GWAS) of serum creatinine levels while on vancomycin in 489 European American individuals and validated findings in three independent cohorts totaling 439 European American individuals. In primary analyses, the chromosome 6q22.31 locus was associated with increased serum creatinine levels while on vancomycin therapy (most significant variant rs2789047, risk allele A, ß = -0.06, p = 1.1 x 10(-7)). SNPs in this region had consistent directions of effect in the validation cohorts, with a meta-p of 1.1 x 10(-7). Variation in this region on chromosome 6, which includes the genes TBC1D32/C6orf170 and GJA1 (encoding connexin43), may modulate risk of vancomycin-induced kidney injury.

Examination of NRCAM, LRRN3, KIAA0716, and LAMB1 as autism candidate genes.

BACKGROUND: A substantial body of research supports a genetic involvement in autism. Furthermore, results from various genomic screens implicate a region on chromosome 7q31 as harboring an autism susceptibility variant. We previously narrowed this 34 cM region to a 3 cM critical region (located between D7S496 and D7S2418) using the Collaborative Linkage Study of Autism (CLSA) chromosome 7 linked families. This interval encompasses about 4.5 Mb of genomic DNA and encodes over fifty known and predicted genes. Four candidate genes (NRCAM, LRRN3, KIAA0716, and LAMB1) in this region were chosen for examination based on their proximity to the marker most consistently cosegregating with autism in these families (D7S1817), their tissue expression patterns, and likely biological relevance to autism. METHODS: Thirty-six intronic and exonic single nucleotide polymorphisms (SNPs) and one microsatellite marker within and around these four candidate genes were genotyped in 30 chromosome 7q31 linked families. Multiple SNPs were used to provide as complete coverage as possible since linkage disequilibrium can vary dramatically across even very short distances within a gene. Analyses of these data used the Pedigree Disequilibrium Test for single markers and a multilocus likelihood ratio test. RESULTS: As expected, linkage disequilibrium occurred within each of these genes but we did not observe significant LD across genes. None of the polymorphisms in NRCAM, LRRN3, or KIAA0716 gave p < 0.05 suggesting that none of these genes is associated with autism susceptibility in this subset of chromosome 7-linked families. However, with LAMB1, the allelic association analysis revealed suggestive evidence for a positive association, including one individual SNP (p = 0.02) and three separate two-SNP haplotypes across the gene (p = 0.007, 0.012, and 0.012). CONCLUSIONS: NRCAM, LRRN3, KIAA0716 are unlikely to be involved in autism. There is some evidence that variation in or near the LAMB1 gene may be involved in autism.

The NEIGHBOR consortium primary open-angle glaucoma genome-wide association study: rationale, study design, and clinical variables.

Primary open-angle glaucoma (POAG) is a common disease with complex inheritance. The identification of genes predisposing to POAG is an important step toward the development of novel gene-based methods of diagnosis and treatment. Genome-wide association studies (GWAS) have successfully identified genes contributing to complex traits such as POAG however, such studies frequently require very large sample sizes, and thus, collaborations and consortia have been of critical importance for the GWAS approach. In this report we describe the formation of the NEIGHBOR consortium, the harmonized case control definitions used for a POAG GWAS, the clinical features of the cases and controls, and the rationale for the GWAS study design.

CDKN2B-AS1 genotype-glaucoma feature correlations in primary open-angle glaucoma patients from the United States.

PURPOSE: To assess the association between single nucleotide polymorphisms (SNPs) of the gene region containing cyclin-dependent kinase inhibitor 2B antisense noncoding RNA (CDKN2B-AS1) and glaucoma features among primary open-angle glaucoma (POAG) patients. DESIGN: Retrospective observational case series. METHODS: We studied associations between 10 CDKN2B-AS1 SNPs and glaucoma features among 976 POAG cases from the Glaucoma Genes and Environment (GLAUGEN) study and 1971 cases from the National Eye Institute Glaucoma Human Genetics Collaboration (NEIGHBOR) consortium. For each patient, we chose the feature from the eye with the higher value. We created cohort-specific multivariable models for glaucoma features and then meta-analyzed the results. RESULTS: For 9 of the 10 protective CDKN2B-AS1 SNPs with minor alleles associated with reduced disease risk (eg, the G allele at rs2157719), POAG patients carrying these minor alleles had smaller cup-to-disc ratio (0.05 units smaller per G allele at diagnosis; 95% CI: -0.08, -0.03; P = 6.23E-05) despite having higher intraocular pressure (IOP) (0.70 mm Hg higher per G allele at DNA collection; 95% CI: 0.40, 1.00; P = 5.45E-06). For the 1 adverse rs3217992 SNP with minor allele A associated with increased disease risk, POAG patients with A alleles had larger cup-to-disc ratio (0.05 units larger per A allele at diagnosis; 95% CI: 0.02, 0.07; P = 4.74E-04) despite having lower IOP (-0.57 mm Hg per A allele at DNA collection; 95% CI: -0.84, -0.29; P = 6.55E-05). CONCLUSION: Alleles of CDKN2B-AS1 SNPs, which influence risk of developing POAG, also modulate optic nerve degeneration among POAG patients, underscoring the role of CDKN2B-AS1 in POAG.

Mechanistic phenotypes: an aggregative phenotyping strategy to identify disease mechanisms using GWAS data.

A single mutation can alter cellular and global homeostatic mechanisms and give rise to multiple clinical diseases. We hypothesized that these disease mechanisms could be identified using low minor allele frequency (MAF<0.1) non-synonymous SNPs (nsSNPs) associated with "mechanistic phenotypes", comprised of collections of related diagnoses. We studied two mechanistic phenotypes: (1) thrombosis, evaluated in a population of 1,655 African Americans; and (2) four groupings of cancer diagnoses, evaluated in 3,009 white European Americans. We tested associations between nsSNPs represented on GWAS platforms and mechanistic phenotypes ascertained from electronic medical records (EMRs), and sought enrichment in functional ontologies across the top-ranked associations. We used a two-step analytic approach whereby nsSNPs were first sorted by the strength of their association with a phenotype. We tested associations using two reverse genetic models and standard additive and recessive models. In the second step, we employed a hypothesis-free ontological enrichment analysis using the sorted nsSNPs to identify functional mechanisms underlying the diagnoses comprising the mechanistic phenotypes. The thrombosis phenotype was solely associated with ontologies related to blood coagulation (Fisher's p = 0.0001, FDR p = 0.03), driven by the F5, P2RY12 and F2RL2 genes. For the cancer phenotypes, the reverse genetics models were enriched in DNA repair functions (p = 2×10-5, FDR p = 0.03) (POLG/FANCI, SLX4/FANCP, XRCC1, BRCA1, FANCA, CHD1L) while the additive model showed enrichment related to chromatid segregation (p = 4×10-6, FDR p = 0.005) (KIF25, PINX1). We were able to replicate nsSNP associations for POLG/FANCI, BRCA1, FANCA and CHD1L in independent data sets. Mechanism-oriented phenotyping using collections of EMR-derived diagnoses can elucidate fundamental disease mechanisms.

Identification of a rare coding variant in complement 3 associated with age-related macular degeneration.

Macular degeneration is a common cause of blindness in the elderly. To identify rare coding variants associated with a large increase in risk of age-related macular degeneration (AMD), we sequenced 2,335 cases and 789 controls in 10 candidate loci (57 genes). To increase power, we augmented our control set with ancestry-matched exome-sequenced controls. An analysis of coding variation in 2,268 AMD cases and 2,268 ancestry-matched controls identified 2 large-effect rare variants: previously described p.Arg1210Cys encoded in the CFH gene (case frequency (fcase) = 0.51%; control frequency (fcontrol) = 0.02%; odds ratio (OR) = 23.11) and newly identified p.Lys155Gln encoded in the C3 gene (fcase = 1.06%; fcontrol = 0.39%; OR = 2.68). The variants suggest decreased inhibition of C3 by complement factor H, resulting in increased activation of the alternative complement pathway, as a key component of disease biology.

Using genetic variation and environmental risk factor data to identify individuals at high risk for age-related macular degeneration.

A major goal of personalized medicine is to pre-symptomatically identify individuals at high risk for disease using knowledge of each individual's particular genetic profile and constellation of environmental risk factors. With the identification of several well-replicated risk factors for age-related macular degeneration (AMD), the leading cause of legal blindness in older adults, this previously unreachable goal is beginning to seem less elusive. However, recently developed algorithms have either been much less accurate than expected, given the strong effects of the identified risk factors, or have not been applied to independent datasets, leaving unknown how well they would perform in the population at large. We sought to increase accuracy by using novel modeling strategies, including multifactor dimensionality reduction (MDR) and grammatical evolution of neural networks (GENN), in addition to the traditional logistic regression approach. Furthermore, we rigorously designed and tested our models in three distinct datasets: a Vanderbilt-Miami (VM) clinic-based case-control dataset, a VM family dataset, and the population-based Age-related Maculopathy Ancillary (ARMA) Study cohort. Using a consensus approach to combine the results from logistic regression and GENN models, our algorithm was successful in differentiating between high- and low-risk groups (sensitivity 77.0%, specificity 74.1%). In the ARMA cohort, the positive and negative predictive values were 63.3% and 70.7%, respectively. We expect that future efforts to refine this algorithm by increasing the sample size available for model building, including novel susceptibility factors as they are discovered, and by calibrating the model for diverse populations will improve accuracy.

Dissection of chromosome 16p12 linkage peak suggests a possible role for CACNG3 variants in age-related macular degeneration susceptibility.

PURPOSE: Age-related macular degeneration (AMD) is a complex disorder of the retina, characterized by drusen, geographic atrophy, and choroidal neovascularization. Cigarette smoking and the genetic variants CFH Y402H, ARMS2 A69S, CFB R32Q, and C3 R102G have been strongly and consistently associated with AMD. Multiple linkage studies have found evidence suggestive of another AMD locus on chromosome 16p12 but the gene responsible has yet to be identified. METHODS: In the initial phase of the study, single-nucleotide polymorphisms (SNPs) across chromosome 16 were examined for linkage and/or association in 575 Caucasian individuals from 148 multiplex and 77 singleton families. Additional variants were tested in an independent dataset of unrelated cases and controls. According to these results, in combination with gene expression data and biological knowledge, five genes were selected for further study: CACNG3, HS3ST4, IL4R, Q7Z6F8, and ITGAM. RESULTS: After genotyping additional tagging SNPs across each gene, the strongest evidence for linkage and association was found within CACNG3 (rs757200 nonparametric LOD* = 3.3, APL (association in the presence of linkage) P = 0.06, and rs2238498 MQLS (modified quasi-likelihood score) P = 0.006 in the families; rs2283550 P = 1.3 × 10(-6), and rs4787924 P = 0.002 in the case-control dataset). After adjusting for known AMD risk factors, rs2283550 remained strongly associated (P = 2.4 × 10(-4)). Furthermore, the association signal at rs4787924 was replicated in an independent dataset (P = 0.035) and in a joint analysis of all the data (P = 0.001). CONCLUSIONS: These results suggest that CACNG3 is the best candidate for an AMD risk gene within the 16p12 linkage peak. More studies are needed to confirm this association and clarify the role of the gene in AMD pathogenesis.