Genome-wide association studies have contributed extensively to the discovery of disease-associated common variants. However, the genetic contribution to complex traits is still largely difficult to interpret. We report a genome-wide association study of 2394 cases and 2393 controls for age-related macular degeneration (AMD) via whole-genome sequencing, with 46.9 million genetic variants. Our study reveals significant single-variant association signals at four loci and independent gene-based signals in CFH, C2, C3, and NRTN. Using data from the Exome Aggregation Consortium (ExAC) for a gene-based test, we demonstrate an enrichment of predicted rare loss-of-function variants in CFH, CFI, and an as-yet unreported gene in AMD, ORMDL2. Our method of using a large variant list without individual-level genotypes as an external reference provides a flexible and convenient approach to leverage the publicly available variant datasets to augment the search for rare variant associations, which can explain additional disease risk in AMD.
Genome-wide association studies (GWAS) have identified 52 independent variants at 34 genetic loci that are associated with age-related macular degeneration (AMD), the most common cause of incurable vision loss in the elderly worldwide. However, causal genes at the majority of these loci remain unknown. In this study, we performed whole exome sequencing of 264 individuals from 63 multiplex families with AMD and analyzed the data for rare protein-altering variants in candidate target genes at AMD-associated loci. Rare coding variants were identified in the CFH, PUS7, RXFP2, PHF12 and TACC2 genes in three or more families. In addition, we detected rare coding variants in the C9, SPEF2 and BCAR1 genes, which were previously suggested as likely causative genes at respective AMD susceptibility loci. Identification of rare variants in the CFH and C9 genes in our study validated previous reports of rare variants in complement pathway genes in AMD. We then extended our exome-wide analysis and identified rare protein-altering variants in 13 genes outside the AMD-GWAS loci in three or more families. Two of these genes, SCN10A and KIR2DL4, are of interest because variants in these genes also showed association with AMD in case-control cohorts, albeit not at the level of genome-wide significance. Our study presents the first large-scale, exome-wide analysis of rare variants in AMD. Further independent replications and molecular investigation of candidate target genes, reported here, would assist in gaining novel insights into mechanisms underlying AMD pathogenesis.
Andrea R. Waksmunski, Michelle Grunin, Tyler G. Kinzy, Robert P. Igo Jr, Jonathan L. Haines, and Jessica N. Cooke Bailey; for the International Age-Related Macular Degeneration Genomics Consortium Department of Genetics and Genome Sciences, Case Western Reserve University, Cleveland, Ohio, United States Cleveland Institute for Computational Biology, Case Western Reserve University, Cleveland, Ohio, United States Department of Population and Quantitative Health Sciences, Case Western Reserve University, Cleveland, Ohio, United States
Purpose: X-linked retinitis pigmentosa can manifest in female carriers with widely variable severity, whereas others remain unaffected. The contribution of X-chromosome inactivation (XCI) to phenotypic variation has been postulated but not demonstrated. Furthermore, the impact of genotype and genetic modifiers has been demonstrated in affected males but has not been well established in female carriers. The purpose of this study was to describe the scope of clinical phenotype in female carriers with mutations in RPGR and quantify the contribution of genotype, genetic modifiers, and XCI to phenotypic severity. Design: Cohort study. Participants: Seventy-seven female carriers with RPGR mutations from 41 pedigrees. Methods: Coding single nucleotide polymorphisms were sequenced in candidate genetic modifier genes encoding known RPGR-interacting proteins. X-chromosome inactivation ratios were determined in genomic DNA isolated from blood (n = 42) and saliva (n = 20) using methylation status of X-linked polymorphic repeats. These genetic data were compared with disease severity based on quantitative clinical parameters. Main Outcome Measures: Visual acuity, Humphrey visual field (HVF) results, full-field electroretinography results, and dark adaptation. Results: Most individuals at all ages were mildly affected or unaffected, whereas those who progressed to moderate or severe vision loss were older than 30 years. RPGR genotype was not associated with clinical severity. The D1264N variant in RPGRIP1L was associated with more severe disease. Skewed XCI toward inactivation of the normal RPGR allele was associated with more severe disease. The XCI ratio in both blood and saliva was a predictor of visual function as measured by HVF diameter, rod amplitude, flicker amplitude, and flicker implicit time. For carriers with extreme XCI skewing of 80:20 or more, 57% were affected severely compared with 8% for those with XCI of less than 80:20 (P = 0.002). Conclusions: Female carriers with mutations in RPGR demonstrate widely variable clinical severity. X-chromosome inactivation ratios correlate with clinical severity and may serve as a predictor of clinically significant disease. Because RPGR gene therapy trials are underway, a future imperative exists to determine which carriers require intervention and when to intervene. X-chromosome inactivation analysis may be useful for identifying candidates for early intervention. (C) 2019 by the American Academy of Ophthalmology
Retinal dystrophies are a phenotypically and genetically complex group of conditions. Because of this complexity, it can be challenging in many families to determine the inheritance based on pedigree analysis alone. Clinical examinations were performed and blood samples were collected from a North American (M1186) and a consanguineous Pakistani (PKRD168) pedigree affected with two different retinal dystrophies (RD). Based on the structure of the pedigrees, inheritance patterns in the families were difficult to determine. In one family, linkage analysis was performed with markers on X-chromosome. In the second family, whole-exome sequencing (WES) was performed. Subsequent Sanger sequencing of genes of interest was performed. Linkage and haplotype analysis localized the disease interval to a 70 Mb region on the X chromosome that encompassed RP2 and RPGR in M1186. The disease haplotype segregated with RD in all individuals except for an unaffected man (IV:3) and his affected son (V:1) in this pedigree. Subsequent analysis identified a novel RPGR mutation (p. Lys857Glu fs221X) in all affected members of M1186 except V:1. This information suggests that there is an unidentified second cause of retinitis pigmentosa (RP) within the family. A novel two-base-pair deletion (p. Tyr565Ter fsX) in CHM (choroideremia) was found to segregate with RD in PKRD168. This paper highlights the challenges of interpreting family history in families with RD and reports on the identification of novel mutations in two RD families.
Autoimmune retinopathy (AIR) causes rapidly progressive vision loss that is treatable but often is confused with other forms of retinal degeneration including retinitis pigmentosa (RP). Measurement of anti-retinal antibodies (ARA) by Western blot is a commonly used laboratory assay that supports the diagnosis yet does not reflect current disease activity. To search for better diagnostic indicators, this study was designed to compare immune biomarkers and responses toward the retinal protein, recoverin, between newly diagnosed AIR patients, slow progressing RP patients and healthy controls. All individuals had measurable anti-recoverin IgG and IgM antibodies by ELISA regardless of disease status or Western blot results. Many AIR patients had elevated anti-recoverin IgG1 levels and a strong cellular response toward recoverin dominated by IFNγ. RP patients and controls responded to recoverin with a lower IFNγ response that was balanced by IL-10 production. Both AIR and RP patients displayed lower levels of total peripheral blood mononuclear cells that were due to reductions of CD4+ TH cells. A comparison of messenger RNA (mRNA) for immune-related genes in whole blood of AIR patients versus RP patients or controls indicated lower expression of ATG5 and PTPN22 and higher expression of several genes involved in TH cell signaling/transcription and adhesion. These data indicate that an immune response toward recoverin is normal in humans, but that in AIR patients the balance shifts dramatically toward higher IFNγ production and cellular activation.
Citation for published version: Stanton, CM, Yates, JRW, den Hollander, AI, Seddon, JM, Swaroop, A, Stambolian, D, Fauser, S, Hoyng, C, Yu, Y, Atsuhiro, K, Branham, K, Othman, M, Chen, W, Kortvely, E, Chalmers, K, Hayward, C, Moore, AT, Dhillon, B, Ueffing, M & Wright, AF 2011, 'Complement factor D in age-related macular degeneration' Investigative Ophthalmology & Visual Science, vol 52, no. 12, pp. 8828-34. DOI: 10.1167/iovs.11-7933
Bicuspid aortic valve (BAV) is a heritable congenital heart defect and an important risk factor for valvulopathy and aortopathy. Here we report a genome-wide association scan of 466 BAV cases and 4,660 age, sex and ethnicity-matched controls with replication in up to 1,326 cases and 8,103 controls. We identify association with a noncoding variant 151 kb from the gene encoding the cardiac-specific transcription factor, GATA4, and near-significance for p.Ser377Gly in GATA4. GATA4 was interrupted by CRISPR-Cas9 in induced pluripotent stem cells from healthy donors. The disruption of GATA4 significantly impaired the transition from endothelial cells into mesenchymal cells, a critical step in heart valve development.
Iris Heid, Gonçalo Abecasis, Sudha Iyengar and colleagues report the results of a large genome-wide association meta-analysis of macular degeneration based on over 43,000 subjects. They identify 16 new risk loci, including some very rare coding variants. Advanced age-related macular degeneration (AMD) is the leading cause of blindness in the elderly, with limited therapeutic options. Here we report on a study of >12 million variants, including 163,714 directly genotyped, mostly rare, protein-altering variants. Analyzing 16,144 patients and 17,832 controls, we identify 52 independently associated common and rare variants (P < 5 × 10−8) distributed across 34 loci. Although wet and dry AMD subtypes exhibit predominantly shared genetics, we identify the first genetic association signal specific to wet AMD, near MMP9 (difference P value = 4.1 × 10−10). Very rare coding variants (frequency <0.1%) in CFH, CFI and TIMP3 suggest causal roles for these genes, as does a splice variant in SLC16A8. Our results support the hypothesis that rare coding variants can pinpoint causal genes within known genetic loci and illustrate that applying the approach systematically to detect new loci requires extremely large sample sizes.
BACKGROUND:Sequencing studies of exonic regions aim to identify rare variants contributing to complex traits. With high coverage and large sample size, these studies tend to apply simple variant calling algorithms. However, coverage is often heterogeneous; sites with insufficient coverage may benefit from sophisticated calling algorithms used in low-coverage sequencing studies. We evaluate the potential benefits of different calling strategies by performing a comparative analysis of variant calling methods on exonic data from 202 genes sequenced at 24x in 7,842 individuals. We call variants using individual-based, population-based and linkage disequilibrium (LD)-aware methods with stringent quality control. We measure genotype accuracy by the concordance with on-target GWAS genotypes and between 80 pairs of sequencing replicates. We validate selected singleton variants using capillary sequencing.RESULTS:Using these calling methods, we detected over 27,500 variants at the targeted exons; >57% were singletons. The singletons identified by individual-based analyses were of the highest quality. However, individual-based analyses generated more missing genotypes (4.72%) than population-based (0.47%) and LD-aware (0.17%) analyses. Moreover, individual-based genotypes were the least concordant with array-based genotypes and replicates. Population-based genotypes were less concordant than genotypes from LD-aware analyses with extended haplotypes. We reanalyzed the same dataset with a second set of callers and showed again that the individual-based caller identified more high-quality singletons than the population-based caller. We also replicated this result in a second dataset of 57 genes sequenced at 127.5x in 3,124 individuals.CONCLUSIONS:We recommend population-based analyses for high quality variant calls with few missing genotypes. With extended haplotypes, LD-aware methods generate the most accurate and complete genotypes. In addition, individual-based analyses should complement the above methods to obtain the most singleton variants.
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.
Elevated intraocular pressure (IOP) is a major risk factor for glaucoma and is influenced by genetic and environmental factors. Recent genome-wide association studies (GWAS) reported associations with IOP at TMCO1 and GAS7, and with primary open-angle glaucoma (POAG) at CDKN2B-AS1, CAV1/CAV2, and SIX1/SIX6. To identify novel genetic variants and replicate the published findings, we performed GWAS and meta-analysis of IOP in >6,000 subjects of European ancestry collected in three datasets: the NEI Glaucoma Human genetics collaBORation, GLAUcoma Genes and ENvironment study, and a subset of the Age-related Macular Degeneration-Michigan, Mayo, AREDS and Pennsylvania study. While no signal achieved genome-wide significance in individual datasets, a meta-analysis identified significant associations with IOP at TMCO1 (rs7518099-G, p = 8.0 × 10−8). Focused analyses of five loci previously reported for IOP and/or POAG, i.e., TMCO1, CDKN2B-AS1, GAS7, CAV1/CAV2, and SIX1/SIX6, revealed associations with IOP that were largely consistent across our three datasets, and replicated the previously reported associations in both effect size and direction. These results confirm the involvement of common variants in multiple genomic regions in regulating IOP and/or glaucoma risk.
Age-related macular degeneration (AMD) is a leading cause of visual impairment worldwide. Aberrant DNA methylation within the promoter of IL17RC in peripheral blood mononuclear cells has recently been reported in AMD. To validate this association, we examined DNA methylation of the IL17RC promoter in peripheral blood. First, we used Illumina Human Methylation450 Bead Arrays, a widely accepted platform for measuring global DNA methylation. Second, methylation status at multiple sites within the IL17RC promoter was determined by bisulfite pyrosequencing in two cohorts. Third, a methylation-sensitive quantitative PCR-based assay was performed on a subset of samples. In contrast to previous findings, we did not find evidence of differential methylation between AMD cases and age-matched controls. We conclude that hypomethylation within the IL17RC gene promoter in peripheral blood is not suitable for use as a clinical biomarker of AMD. This study highlights the need for considerable replication of epigenetic association studies prior to clinical application.
Importance For patients with X-linked retinitis pigmentosa and clinicians alike, phenotypic variability can be challenging because it complicates counseling regarding patients’ likely visual prognosis. Objective To evaluate the clinical findings from patients with X-linked retinitis pigmentosa with 13 distinctRPGRmutations and assess for phenotypic concordance or variability. Design Retrospective medical record review of data collected from 1985 to 2011. Setting Kellogg Eye Center, University of Michigan. Patients A total of 42 patients with X-linked retinitis pigmentosa with mutations inRPGR. Age at first visit ranged from 4 to 53 years, with follow-up ranging from 1 to 11 visits (median follow-up time, 5.5 years; range, 1.4-32.7 years, for 23 patients with >1 visit). Main Outcomes and Measures Clinical data assessed for concordance included visual acuity (VA), Goldmann visual fields (GVFs), and full-field electroretinography (ERG). Electroretinography phenotype (cone-rod vs rod-cone dysfunction) was defined by the extent of photopic vs scotopic abnormality. Qualitative GVF phenotype was determined by the GVF pattern, where central or peripheral loss suggested cone or rod dysfunction, respectively. Goldmann visual fields were also quantified and compared among patients. Results Each mutation was detected in 2 or more related or unrelated patients. Five mutations in 11 patients displayed strong concordance of VA, while 4 mutations in 16 patients revealed moderate concordance of VA. A definitive cone-rod or rod-cone ERG pattern consistent among patients was found in 6 of 13 mutations (46.2%); the remaining mutations were characterized by patients demonstrating both phenotypes or who had limited data or nonrecordable ERG values. Concordant GVF phenotypes (7 rod-cone pattern vs 4 cone-rod pattern) were seen in 11 of 13 mutations (84.6%). All 6 mutations displaying a constant ERG pattern within the mutation group revealed a GVF phenotype consistent with the ERG findings. Conclusions and Relevance While VA and ERG phenotypes are concordant in only some patients carrying identical mutations, assessment of GVF phenotypes revealed stronger phenotypic conservation. Phenotypic concordance is important for establishing proper counseling of patients diagnosed as having X-linked retinitis pigmentosa, as well as for establishing accurate patient selection and efficacy monitoring in therapeutic trials.
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.
Animal models of human disease are an invaluable component of studies aimed at understanding disease pathogenesis and therapeutic possibilities. Mutations in the gene encoding retinitis pigmentosa GTPase regulator (RPGR) are the most common cause of X-linked retinitis pigmentosa (XLRP) and are estimated to cause 20% of all retinal dystrophy cases. A majority of RPGR mutations are present in ORF15, the purine-rich terminal exon of the predominant splice-variant expressed in retina. Here we describe the genetic and phenotypic characterization of the retinal degeneration 9 (Rd9) strain of mice, a naturally occurring animal model of XLRP. Rd9 mice were found to carry a 32-base-pair duplication within ORF15 that causes a shift in the reading frame that introduces a premature-stop codon. Rpgr ORF15 transcripts, but not protein, were detected in retinas from Rd9/Y male mice that exhibited retinal pathology, including pigment loss and slowly progressing decrease in outer nuclear layer thickness. The levels of rhodopsin and transducin in rod outer segments were also decreased, and M-cone opsin appeared mislocalized within cone photoreceptors. In addition, electroretinogram (ERG) a- and b-wave amplitudes of both Rd9/Y male and Rd9/Rd9 female mice showed moderate gradual reduction that continued to 24 months of age. The presence of multiple retinal features that correlate with findings in individuals with XLRP identifies Rd9 as a valuable model for use in gaining insight into ORF15-associated disease progression and pathogenesis, as well as accelerating the development and testing of therapeutic strategies for this common form of retinal dystrophy.