Genetic and genomic studies have enhanced our understanding of complex neurodegenerative diseases that exert a devastating impact on individuals and society. One such disease, age-related macular degeneration (AMD), is a major cause of progressive and debilitating visual impairment. Since the pioneering discovery in 2005 of complement factor H (CFH) as a major AMD susceptibility gene, extensive investigations have confirmed 19 additional genetic risk loci, and more are anticipated. In addition to common variants identified by now-conventional genome-wide association studies, targeted genomic sequencing and exome-chip analyses are uncovering rare variant alleles of high impact. Here, we provide a critical review of the ongoing genetic studies and of common and rare risk variants at a total of 20 susceptibility loci, which together explain 40–60% of the disease heritability but provide limited power for diagnostic testing of disease risk. Identification of these susceptibility loci has begun to untangle the complex biological pathways underlying AMD pathophysiology, pointing to new testable paradigms for treatment.
complex disease; genetic susceptibility; neurodegeneration; retina; blindness
Refractive error (RE) is a complex, multifactorial disorder characterized by a mismatch between the optical power of the eye and its axial length that causes object images to be focused off the retina. The two major subtypes of RE are myopia (nearsightedness) and hyperopia (farsightedness), which represent opposite ends of the distribution of the quantitative measure of spherical refraction. We performed a fixed effects meta-analysis of genome-wide association results of myopia and hyperopia from 9 studies of European-derived populations: AREDS, KORA, FES, OGP-Talana, MESA, RSI, RSII, RSIII and ERF. One genome-wide significant region was observed for myopia, corresponding to a previously identified myopia locus on 8q12 (p = 1.25×10−8), which has been reported by Kiefer et al. as significantly associated with myopia age at onset and Verhoeven et al. as significantly associated to mean spherical-equivalent (MSE) refractive error. We observed two genome-wide significant associations with hyperopia. These regions overlapped with loci on 15q14 (minimum p value = 9.11×10−11) and 8q12 (minimum p value 1.82×10−11) previously reported for MSE and myopia age at onset. We also used an intermarker linkage- disequilibrium-based method for calculating the effective number of tests in targeted regional replication analyses. We analyzed myopia (which represents the closest phenotype in our data to the one used by Kiefer et al.) and showed replication of 10 additional loci associated with myopia previously reported by Kiefer et al. This is the first replication of these loci using myopia as the trait under analysis. “Replication-level” association was also seen between hyperopia and 12 of Kiefer et al.'s published loci. For the loci that show evidence of association to both myopia and hyperopia, the estimated effect of the risk alleles were in opposite directions for the two traits. This suggests that these loci are important contributors to variation of refractive error across the distribution.
For automatic segmentation of optic disc and cup from color fundus photograph, we describe a fairly general energy function that can naturally fit into a global optimization framework with graph cut. Distinguished from most previous work, our energy function includes priors on the shape & location of disc & cup, the rim thickness and the geometric interaction of “disc contains cup”. These priors together with the effective optimization of graph cut enable our algorithm to generate reliable and robust solutions. Our approach is able to outperform several state-of-the-art segmentation methods, as shown by a set of experimental comparisons with manual delineations and a series of results of correlations with the assessments of a merchant-provided software from Optical Coherence Tomography (OCT) regarding several cup and disc parameters.
To identify myopia susceptibility genes influencing common myopia in 94 African-American and 36 White families.
A prospective study of families with myopia consisting of a minimum of two individuals affected with myopia.
Extended families consisting of at least two siblings affected with myopia were ascertained. A genome-wide linkage scan using 387 markers was conducted by the Center for Inherited Disease Research. Linkage analyses were conducted with parametric and nonparametric methods. Model-free linkage analysis was performed maximizing over penetrance and over dominance (that is, fitting a wide range of both dominant and recessive models).
Under the model-free analysis, the maximum two point heterogeneity logarithm of the odds score (MALOD) was 2.87 at D6S1009 in the White cohort and the maximum multipoint MALOD was 2.42 at D12S373-D12S1042 in the same cohort. The nonpara-metric linkage (NPL) maximum multipoint at D6S1035 had a P value of .005. An overall multipoint NPL score was obtained by combining NPL scores from both populations. The highest combined NPL score was observed at D20S478 with a significant P value of .008. Suggestive evidence of linkage in the White cohort mapped to a previously mapped locus on chromosome 11 at D11S1981 (NPL = 2.14; P = .02).
Suggestive evidence of linkage to myopia in both African Americans and Whites was seen on chromosome 20 and became more significant when the scores were combined for both groups. The locus on chromosome 11 independently confirms a report by Hammond and associates mapping a myopia quantitative trait locus to this region.
Refractive errors, myopia and hyperopia, are the most common causes of visual impairment worldwide. Recent advances in genetics have been utilized to identify a wealth of genetic loci believed to contain susceptibility genes for refractive error. The current genetic evidence confirms that refractive error is influenced by both common and rare variants with a significant environmental component. These studies argue that only by combining genetic knowledge with in vivo measurements of biological states will it be possible to understand the underlying biology of refractive error that will lead to novel therapeutic targets and accurate genetic predictions.
Refractive Error; Myopia; Hyperopia; Genetic Variation; Animal Models; Gene-Environment Interaction
Knowledge of individual ancestry is important for genetic association studies where population structure leads to false positive signals. Estimating individual ancestry with targeted sequence data, which constitutes the bulk of current sequence datasets, is challenging. Here, we propose a new method for accurate estimation of genetic ancestry. Our method skips genotype calling and directly analyzes sequence reads. We validate the method using simulated and empirical data and show that the method can accurately infer worldwide continental ancestry with whole genome shotgun coverage as low as 0.001X. For estimates of fine-scale ancestry within Europe, the method performs well with coverage of 0.1X. At an even finer-scale, the method improves discrimination between exome-sequenced participants originating from different provinces within Finland. Finally, we show that our method can be used to improve case-control matching in genetic association studies and reduce the risk of spurious findings due to population structure.
Visual refractive errors (REs) are complex genetic traits with a largely unknown etiology. To date, genome-wide association studies (GWASs) of moderate size have identified several novel risk markers for RE, measured here as mean spherical equivalent (MSE). We performed a GWAS using a total of 7280 samples from five cohorts: the Age-Related Eye Disease Study (AREDS); the KORA study (‘Cooperative Health Research in the Region of Augsburg’); the Framingham Eye Study (FES); the Ogliastra Genetic Park-Talana (OGP-Talana) Study and the Multiethnic Study of Atherosclerosis (MESA). Genotyping was performed on Illumina and Affymetrix platforms with additional markers imputed to the HapMap II reference panel. We identified a new genome-wide significant locus on chromosome 16 (rs10500355, P = 3.9 × 10−9) in a combined discovery and replication set (26 953 samples). This single nucleotide polymorphism (SNP) is located within the RBFOX1 gene which is a neuron-specific splicing factor regulating a wide range of alternative splicing events implicated in neuronal development and maturation, including transcription factors, other splicing factors and synaptic proteins.
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 revealed two large-effect rare variants; previously described R1210C in the CFH gene (fcase = 0.51%, fcontrol = 0.02%, OR = 23.11), and newly identified K155Q in the C3 gene (fcase = 1.06%, fcontrol = 0.39%, OR = 2.68). The variants suggest decreased inhibition of C3 by Factor H, resulting in increased activation of the alternative complement pathway, as a key component of disease biology.
A previous study of Old Order Amish families has shown association of ocular refraction with markers proximal to matrix metalloproteinase (MMP) genes MMP1 and MMP10 and intragenic to MMP2. We conducted a candidate gene replication study of association between refraction and single nucleotide polymorphisms (SNPs) within these genomic regions.
Candidate gene genetic association study.
2,000 participants drawn from the Age Related Eye Disease Study (AREDS) were chosen for genotyping. After quality control filtering, 1912 individuals were available for analysis.
Microarray genotyping was performed using the HumanOmni 2.5 bead array. SNPs originally typed in the previous Amish association study were extracted for analysis. In addition, haplotype tagging SNPs were genotyped using TaqMan assays. Quantitative trait association analyses of mean spherical equivalent refraction (MSE) were performed on 30 markers using linear regression models and an additive genetic risk model, while adjusting for age, sex, education, and population substructure. Post-hoc analyses were performed after stratifying on a dichotomous education variable. Pointwise (P-emp) and multiple-test study-wise (P-multi) significance levels were calculated empirically through permutation.
Main outcome measures
MSE was used as a quantitative measure of ocular refraction.
The mean age and ocular refraction were 68 years (SD=4.7) and +0.55 D (SD=2.14), respectively. Pointwise statistical significance was obtained for rs1939008 (P-emp=0.0326). No SNP attained statistical significance after correcting for multiple testing. In stratified analyses, multiple SNPs reached pointwise significance in the lower-education group: 2 of these were statistically significant after multiple testing correction. The two highest-ranking SNPs in Amish families (rs1939008 and rs9928731) showed pointwise P-emp<0.01 in the lower-education stratum of AREDS participants.
We show suggestive evidence of replication of an association signal for ocular refraction to a marker between MMP1 and MMP10. We also provide evidence of a gene-environment interaction between previously-reported markers and education on refractive error. Variants in MMP1- MMP10 and MMP2 regions appear to affect population variation in ocular refraction in environmental conditions less favorable for myopia development.
refraction; refractive error; myopia; association study; gene-environment interaction; matrix metalloproteinase; MMP; genetics
Refractive error is the most common eye disorder worldwide, and a prominent cause of blindness. Myopia affects over 30% of Western populations, and up to 80% of Asians. The CREAM consortium conducted genome-wide meta-analyses including 37,382 individuals from 27 studies of European ancestry, and 8,376 from 5 Asian cohorts. We identified 16 new loci for refractive error in subjects of European ancestry, of which 8 were shared with Asians. Combined analysis revealed 8 additional loci. The new loci include genes with functions in neurotransmission (GRIA4), ion channels (KCNQ5), retinoic acid metabolism (RDH5), extracellular matrix remodeling (LAMA2, BMP2), and eye development (SIX6, PRSS56). We also confirmed previously reported associations with GJD2 and RASGRF1. Risk score analysis using associated SNPs showed a tenfold increased risk of myopia for subjects with the highest genetic load. Our results, accumulated across independent multi-ethnic studies, considerably advance understanding of mechanisms involved in refractive error and myopia.
Copy number variations (CNVs) are associated with many complex diseases. Next generation sequencing data enable one to identify precise CNV breakpoints to better under the underlying molecular mechanisms and to design more efficient assays. Using the CIGAR strings of the reads, we develop a method that can identify the exact CNV breakpoints, and in cases when the breakpoints are in a repeated region, the method reports a range where the breakpoints can slide. Our method identifies the breakpoints of a CNV using both the positions and CIGAR strings of the reads that cover breakpoints of a CNV. A read with a long soft clipped part (denoted as S in CIGAR) at its 3′(right) end can be used to identify the 5′(left)-side of the breakpoints, and a read with a long S part at the 5′ end can be used to identify the breakpoint at the 3′-side. To ensure both types of reads cover the same CNV, we require the overlapped common string to include both of the soft clipped parts. When a CNV starts and ends in the same repeated regions, its breakpoints are not unique, in which case our method reports the left most positions for the breakpoints and a range within which the breakpoints can be incremented without changing the variant sequence. We have implemented the methods in a C++ package intended for the current Illumina Miseq and Hiseq platforms for both whole genome and exon-sequencing. Our simulation studies have shown that our method compares favorably with other similar methods in terms of true discovery rate, false positive rate and breakpoint accuracy. Our results from a real application have shown that the detected CNVs are consistent with zygosity and read depth information. The software package is available at http://statgene.med.upenn.edu/softprog.html.
structural variation; breakpoint; duplication; deletion; exon sequencing
Refractive error is a complex trait with multiple genetic and environmental risk factors, and is the most common cause of preventable blindness worldwide. The common nature of the trait suggests the presence of many genetic factors that individually may have modest effects. To achieve an adequate sample size to detect these common variants, large, international collaborations have formed. These consortia typically use meta-analysis to combine multiple studies from many different populations. This approach is robust to differences between populations; however, it does not compensate for the different haplotypes in each genetic background evidenced by different alleles in linkage disequilibrium with the causative variant. We used the Age-Related Eye Disease Study (AREDS) cohort to replicate published significant associations at two loci on chromosome 15 from two genome-wide association studies (GWASs). The single nucleotide polymorphisms (SNPs) that exhibited association on chromosome 15 in the original studies did not show evidence of association with refractive error in the AREDS cohort. This paper seeks to determine whether the non-replication in this AREDS sample may be due to the limited number of SNPs chosen for replication.
We selected all SNPs genotyped on the Illumina Omni2.5v1_B array or custom TaqMan assays or imputed from the GWAS data, in the region surrounding the SNPs from the Consortium for Refractive Error and Myopia study. We analyzed the SNPs for association with refractive error using standard regression methods in PLINK. The effective number of tests was calculated using the Genetic Type I Error Calculator.
Although use of the same SNPs used in the Consortium for Refractive Error and Myopia study did not show any evidence of association with refractive error in this AREDS sample, other SNPs within the candidate regions demonstrated an association with refractive error. Significant evidence of association was found using the hyperopia categorical trait, with the most significant SNPs rs1357179 on 15q14 (p=1.69×10−3) and rs7164400 on 15q25 (p=8.39×10−4), which passed the replication thresholds.
This study adds to the growing body of evidence that attempting to replicate the most significant SNPs found in one population may not be significant in another population due to differences in the linkage disequilibrium structure and/or allele frequency. This suggests that replication studies should include less significant SNPs in an associated region rather than only a few selected SNPs chosen by a significance threshold.
Age-related macular degeneration (AMD) is the most common cause of incurable visual impairment in high-income countries. Previous studies report inconsistent associations between AMD and apolipoprotein E (APOE), a lipid transport protein involved in low-density cholesterol modulation. Potential interaction between APOE and sex, and smoking status, has been reported. We present a pooled analysis (n=21,160) demonstrating associations between late AMD and APOε4 (OR=0.72 per haplotype; CI: 0.65–0.74; P=4.41×10−11) and APOε2 (OR=1.83 for homozygote carriers; CI: 1.04–3.23; P=0.04), following adjustment for age-group and sex within each study and smoking status. No evidence of interaction between APOE and sex or smoking was found. Ever smokers had significant increased risk relative to never smokers for both neovascular (OR=1.54; CI: 1.38–1.72; P=2.8×10−15) and atrophic (OR=1.38; CI: 1.18–1.61; P=3.37×10−5) AMD but not early AMD (OR=0.94; CI: 0.86–1.03; P=0.16), implicating smoking as a major contributing factor to disease progression from early signs to the visually disabling late forms. Extended haplotype analysis incorporating rs405509 did not identify additional risks beyondε2 and ε4 haplotypes. Our expanded analysis substantially improves our understanding of the association between the APOE locus and AMD. It further provides evidence supporting the role of cholesterol modulation, and low-density cholesterol specifically, in AMD disease etiology.
age-related macular degeneration; AMD; apolipoprotein E; APOE; case-control association study
Variation in the apolipoprotein E gene (APOE) has been reported to be associated with longevity in humans. The authors assessed the allelic distribution of APOE isoforms ε2, ε3, and ε4 among 10,623 participants from 15 case-control and cohort studies of age-related macular degeneration (AMD) in populations of European ancestry (study dates ranged from 1990 to 2009). The authors included only the 10,623 control subjects from these studies who were classified as having no evidence of AMD, since variation within the APOE gene has previously been associated with AMD. In an analysis stratified by study center, gender, and smoking status, there was a decreasing frequency of the APOE ε4 isoform with increasing age (χ2 for trend = 14.9 (1 df); P = 0.0001), with a concomitant increase in the ε3 isoform (χ2 for trend = 11.3 (1 df); P = 0.001). The association with age was strongest in ε4 homozygotes; the frequency of ε4 homozygosity decreased from 2.7% for participants aged 60 years or less to 0.8% for those over age 85 years, while the proportion of participants with the ε3/ε4 genotype decreased from 26.8% to 17.5% across the same age range. Gender had no significant effect on the isoform frequencies. This study provides strong support for an association of the APOE gene with human longevity.
aged; apolipoprotein E2; apolipoprotein E3; apolipoprotein E4; apolipoproteins E; longevity; meta-analysis; multicenter study
Myopia is a complex genetic disorder and a common cause of visual impairment among working age adults. Genome-wide association studies have identified susceptibility loci on chromosomes 15q14 and 15q25 in Caucasian populations of European ancestry. Here, we present a confirmation and meta-analysis study in which we assessed whether these two loci are also associated with myopia in other populations. The study population comprised 31 cohorts from the Consortium of Refractive Error and Myopia (CREAM) representing 4 different continents with 55,177 individuals; 42,845 Caucasians and 12,332 Asians. We performed a meta-analysis of 14 single nucleotide polymorphisms (SNPs) on 15q14 and 5 SNPs on 15q25 using linear regression analysis with spherical equivalent as a quantitative outcome, adjusted for age and sex. We calculated the odds ratio (OR) of myopia versus hyperopia for carriers of the top-SNP alleles using a fixed effects meta-analysis. At locus 15q14, all SNPs were significantly replicated, with the lowest P value 3.87 × 10−12 for SNP rs634990 in Caucasians, and 9.65 × 10−4 for rs8032019 in Asians. The overall meta-analysis provided P value 9.20 × 10−23 for the top SNP rs634990. The risk of myopia versus hyperopia was OR 1.88 (95 % CI 1.64, 2.16, P < 0.001) for homozygous carriers of the risk allele at the top SNP rs634990, and OR 1.33 (95 % CI 1.19, 1.49, P < 0.001) for heterozygous carriers. SNPs at locus 15q25 did not replicate significantly (P value 5.81 × 10−2 for top SNP rs939661). We conclude that common variants at chromosome 15q14 influence susceptibility for myopia in Caucasian and Asian populations world-wide.
Electronic supplementary material
The online version of this article (doi:10.1007/s00439-012-1176-0) contains supplementary material, which is available to authorized users.
Complement factor D catalyzes a critical step in the alternative complement activation pathway. The authors report a significant elevation in plasma CFD concentrations in age-related macular degeneration (AMD) patients compared with controls and a weak genetic association between CFD gene variants and AMD.
To examine the role of complement factor D (CFD) in age-related macular degeneration (AMD) by analysis of genetic association, copy number variation, and plasma CFD concentrations.
Single nucleotide polymorphisms (SNPs) in the CFD gene were genotyped and the results analyzed by binary logistic regression. CFD gene copy number was analyzed by gene copy number assay. Plasma CFD was measured by an enzyme-linked immunosorbent assay.
Genetic association was found between CFD gene SNP rs3826945 and AMD (odds ratio 1.44; P = 0.028) in a small discovery case-control series (462 cases and 325 controls) and replicated in a combined cohorts meta-analysis of 4765 cases and 2693 controls, with an odds ratio of 1.11 (P = 0.032), with the association almost confined to females. Copy number variation in the CFD gene was identified in 13 out of 640 samples examined but there was no difference in frequency between AMD cases (1.3%) and controls (2.7%). Plasma CFD concentration was measured in 751 AMD cases and 474 controls and found to be elevated in AMD cases (P = 0.00025). The odds ratio for those in the highest versus lowest quartile for plasma CFD was 1.81. The difference in plasma CFD was again almost confined to females.
CFD regulates activation of the alternative complement pathway, which is implicated in AMD pathogenesis. The authors found evidence for genetic association between a CFD gene SNP and AMD and a significant increase in plasma CFD concentration in AMD cases compared with controls, consistent with a role for CFD in AMD pathogenesis.
Quasi-least squares (QLS) is a two-stage computational approach for estimation of the correlation parameters in the framework of generalized estimating equations. We prove two general results for the class of mixed linear correlation structures: namely, that the stage one QLS estimate of the correlation parameter always exists and is feasible (yields a positive definite estimated correlation matrix) for any correlation structure, while the stage two estimator exists and is unique (and therefore consistent) with probability one, for the class of mixed linear correlation structures. Our general results justify the implementation of QLS for particular members of the class of mixed linear correlation structures that are appropriate for analysis of data from families that may vary in size and composition. We describe the familial structures and implement them in an analysis of optical spherical values in the Old Order Amish (OOA). For the OOA analysis, we show that we would suffer a substantial loss in efficiency, if the familial structures were the true structures, but were misspecified as simpler approximate structures. To help bridge the interface between Statistics and Medicine, we also provide R software so that medical researchers can implement the familial structures in a QLS analysis of their own data.
Quasi-least squares; Linear correlation structure; Mixed correlation structure; Familial data
Complement factor H shows very strong association with Age-related Macular Degeneration (AMD), and recent data suggest that multiple causal variants are associated with disease. To refine the location of the disease associated variants, we characterized in detail the structural variation at CFH and its paralogs, including two copy number polymorphisms (CNP), CNP147 and CNP148, and several rare deletions and duplications. Examination of 34 AMD-enriched extended families (N = 293) and AMD cases (White N = 4210 Indian = 134; Malay = 140) and controls (White N = 3229; Indian = 117; Malay = 2390) demonstrated that deletion CNP148 was protective against AMD, independent of SNPs at CFH. Regression analysis of seven common haplotypes showed three haplotypes, H1, H6 and H7, as conferring risk for AMD development. Being the most common haplotype H1 confers the greatest risk by increasing the odds of AMD by 2.75-fold (95% CI = [2.51, 3.01]; p = 8.31×10−109); Caucasian (H6) and Indian-specific (H7) recombinant haplotypes increase the odds of AMD by 1.85-fold (p = 3.52×10−9) and by 15.57-fold (P = 0.007), respectively. We identified a 32-kb region downstream of Y402H (rs1061170), shared by all three risk haplotypes, suggesting that this region may be critical for AMD development. Further analysis showed that two SNPs within the 32 kb block, rs1329428 and rs203687, optimally explain disease association. rs1329428 resides in 20 kb unique sequence block, but rs203687 resides in a 12 kb block that is 89% similar to a noncoding region contained in ΔCNP148. We conclude that causal variation in this region potentially encompasses both regulatory effects at single markers and copy number.
The development of refractive error is mediated by both environmental and genetic factors. We performed regression-based quantitative trait locus (QTL) linkage analysis on Ashkenazi Jewish families to identify regions in the genome responsible for ocular refraction. We measured refractive error on individuals in 49 multi-generational American families of Ashkenazi Jewish descent. The average family size was 11.1 individuals and was composed of 2.7 generations. Recruitment criteria specified that each family contain at least two myopic members. The mean spherical equivalent refractive error in the sample was −3.46D (SD=3.29) and 87% of individuals were myopic. Microsatellite genotyping with 387 markers was performed on 411 individuals. We performed multipoint regression-based linkage analysis for ocular refraction and a log transformation of the trait using the statistical package Merlin-Regress. Empirical genomewide significance levels were estimated through gene-dropping simulations by generating random genotypes at each of the 387 markers in 200 replicates of our pedigrees. Maximum LOD scores of 9.5 for ocular refraction and 8.7 for log-transformed refraction (LTR) were observed at 49.1 cM on chromosome 1p36 between markers D1S552 and D1S1622. The empirical genomewide significance levels were P=0.065 for ocular refraction and P<0.005 for LTR, providing strong evidence for linkage of refraction to this locus. The inter-marker region containing the peak spans 11 Mb and contains approximately 189 genes. Conclusion: We found genomewide significant evidence for linkage of refractive error to a novel QTL on chromosome 1p36 in an Ashkenazi Jewish population.
Refractive development is influenced by environmental and genetic factors. Genetic studies have identified several regions of linkage to ocular refraction, but none have been carried out in African-derived populations. We performed quantitative trait locus linkage analyses in African-American (AA) families to identify genomic regions responsible for refraction. We recruited 493 AA individuals in 96 families to participate in the Myopia Family Study. Genotyping of 387 microsatellite markers was performed on 398 participants. The mean refraction among genotyped individuals was −2.87 D (SD = 3.58) and myopia of at least 1 D was present in 267 (68%) participants. Multipoint, regression-based, linkage analyses were carried out on a logarithmic transformation of ocular refraction using the statistical package MERLIN-REGRESS. Empirical significance levels were determined via 4,898 whole-genome gene-dropping simulations. Linkage analyses were repeated after clustering families into two subgroups based on admixture proportions as determined by the software package STRUCTURE. Genomewide significant linkage was seen at 47 cM on chromosome 7 (logarithm of the odds ratio (LOD) = 5.87, P = 0.00005). In addition, three regions on chromosomes 2p, 3p and 10p showed suggestive evidence of linkage (LOD > 2, P < 0.005) for ocular refraction. We mapped the first quantitative trait locus for ocular refraction in an AA population to chr.7p15. Two previous studies in European-derived families reported some evidence of linkage to a nearby region, suggesting that this region may contain polymorphisms that mediate refraction across populations. The genomic region under our linkage peak spans ~17 Mb and contains ~170 genes. Further refinement of this region will be pursued in future studies.
ocular refraction; myopia; linkage; genetics; African-American
The role of polymorphisms in matrix metalloproteinase genes with variations in refractive error among Orthodox Ashkenazi Jewish and Old Order Amish families was investigated. The results showed statistically significant genetic associations of refraction with two markers in MMP regions in the Amish population, but not among the Ashkenazim.
Matrix metalloproteinases (MMPs) and tissue inhibitors of metalloproteinases (TIMPs) are involved in scleral extracellular matrix remodeling and have shown differential expression in experimental myopia. The genetic association of refractive error and polymorphisms in MMP and TIMP genes in Old Order Amish (AMISH) and Ashkenazi Jewish (ASHK) families was investigated.
Individuals from 55 AMISH and 63 ASHK families participated in the study. Ascertainment was designed to enrich the families for myopia; the mean spherical equivalent (MSE) refractive error (SD) was −1.61 (2.72) D in the AMISH, and −3.56 (3.32) D in the ASHK. One hundred forty-six common haplotype tagging SNPs covering 14 MMP and 4 TIMP genes were genotyped in 358 AMISH and 535 ASHK participants. Association analyses of MSE and the spherical component of refraction (SPH) were performed separately for the AMISH and the ASHK. Bonferroni-corrected significance thresholds and local false discovery rates were used to account for multiple testing.
After they were filtered for quality-control, 127 SNPs were included in the analyses. No polymorphisms showed statistically significant association to refraction in the ASHK (minimum P = 0.0132). In AMISH, two SNPs showed evidence of association with refractive phenotypes: rs1939008 (P = 0.00016 for SPH); and rs9928731 (P = 0.00026 for SPH). These markers were each estimated to explain <5% of the variance of SPH in the AMISH sample.
Statistically significant genetic associations of ocular refraction to polymorphisms near MMP1 and within MMP2 were identified in the AMISH but not among the ASHK families. The results suggest that the MMP1 and MMP2 genes are involved in refractive variation in the AMISH. Genetic and/or environmental heterogeneity most likely contribute to differences in association results between ethnic groups.
Despite many years of research, most of the genetic factors contributing to myopia development remain unknown. Genetic studies have pointed to a strong inherited component, but although many candidate regions have been implicated, few genes have been positively identified.
We have previously reported 2 genomewide linkage scans in a population of 63 highly aggregated Ashkenazi Jewish families that identified a locus on chromosome 22. Here we used ordered subset analysis (OSA), conditioned on non-parametric linkage to chromosome 22 to detect other chromosomal regions which had evidence of linkage to myopia in subsets of the families, but not the overall sample.
Strong evidence of linkage to a 19-cM linkage interval with a peak OSA nonparametric allele-sharing logarithm-of-odds (LOD) score of 3.14 on 20p12-q11.1 (ΔLOD=2.39, empirical p=0.029) was identified in a subset of 20 families that also exhibited strong evidence of linkage to chromosome 22. One other locus also presented with suggestive LOD scores >2.0 on chromosome 11p14-q14 and one locus on chromosome 6q22-q24 had an OSA LOD score=1.76 (ΔLOD=1.65, empirical p=0.02).
The chromosome 6 and 20 loci are entirely novel and appear linked in a subset of families whose myopia is known to be linked to chromosome 22. The chromosome 11 locus overlaps with the known Myopia-7 (MYP7, OMIM 609256) locus. Using ordered subset analysis allows us to find additional loci linked to myopia in subsets of families, and underlines the complex genetic heterogeneity of myopia even in highly aggregated families and genetically isolated populations such as the Ashkenazi Jews.
Genomewide linkage scans were performed in Caucasian (CAUC) and Old Order Amish (OOA) families to identify genomic regions containing genes responsible for refractive error control. We also performed a meta-analysis by combining these results with our previous linkage results from Ashkenazi Jewish (ASHK) and African American (AFRAM) families.
Two hundred seventy-one CAUC and 411 OOA participants (36 and 61 families, respectively) were recruited to participate in the Myopia Family Study. Recruitment criteria were designed to enrich the sample for multiplex myopic families. Genomewide, model-free, multipoint linkage analyses were performed separately for each population by using >370 microsatellite markers. Empirical significance levels were determined via gene-dropping simulations. A meta-analysis was performed by combining linkage results from the CAUC, OOA, AFRAM, and ASHK samples, and results were compared to previously reported loci for myopia and refraction.
Suggestive evidence of linkage was found at 12q24 (LOD = 4.583, P = 0.00037) and 4q21 (LOD = 2.72, P = 0.0028) in the CAUC sample and at 5qter (LOD = 3.271, P = 0.0014) in the OOA. Meta-analysis linkage results were largely driven by population-specific signals from ASHK and AFRAM families. The meta-analysis showed suggestive evidence of linkage to 4q21-22 (meta-P = 0.00214) adjacent to the previously reported MYP9 and MYP11 loci.
The results showed suggestive evidence of linkage of ocular refraction to 12q24 and 4q21 in CAUC and to 5qter in OOA families. The meta-analysis supports the view that several genes play a role in refractive development across populations. In MFS families, four broad genomic regions (on 1p, 4q, 7p, and 12q) most likely contain genes that influence ocular refraction.
Multiple studies demonstrate a strong association between three variants at chromosome 10q26 – rs10490924, del443ins54, and rs11200638 – near the age-related maculopathy susceptibility 2 (ARMS2) and high-temperature requirement factor A1 (HTRA1) genes with susceptibility to age-related macular degeneration (AMD). In different reports, the del443ins54 and rs11200638 variants are suggested to affect ARMS2 mRNA stability and/or HTRA1 mRNA expression, respectively. The goal of this study is to examine whether these AMD-associated variants alter expression levels of ARMS2 and HTRA1 in human retina samples.
Genomic DNA and total RNA were obtained from 35 human retinas (three young controls, average age=32 years; twenty aged controls, average age=72 years; and twelve AMD retinas, average age=77 years) using standard procedures. As ARMS2 exhibits higher expression in the human placenta, we also included eighteen placenta samples in our analysis. Four polymorphisms – rs2736911, rs10490924, del443ins54, and rs11200638 – were genotyped by PCR followed by sequencing. Expression of ARMS2, HTRA1 and three endogenous control genes (rRNA [rRNA], hypoxanthine phosphoribosyltransferase 1 [HPRT1], and glyceraldehyde-3-phosphate dehydrogenase [GAPDH]) was measured by real-time quantitative RT–PCR using Taqman gene expression or SYBR Green assays.
ARMS2 and HTRA1 mRNA levels did not show a significant difference in expression among the control (young and elderly) and AMD retinas. No association of del443ins54 and rs11200638 variants was detected with mRNA expression levels of ARMS2 or HTRA1 in the retina. Human placenta samples showed high variability in expression levels.
We did not find association between AMD susceptibility variants at 10q26 and steady-state expression levels of either ARMS2 or HTRA1 in the human retina.
A previous genome-wide study in Orthodox Ashkenazi Jewish pedigrees showed significant linkage of ocular refraction to a Quantitative Trait Locus (QTL) on 1p34-36.1. We carried out a fine-mapping study of this region in Orthodox Ashkenazi Jewish (ASHK) and Old Order Amish (OOA) families to confirm linkage and narrow the candidate region.
Families were recruited from ASHK and OOA American communities. The samples included: 402 individuals in 53 OOA families; and 596 members in 68 ASHK families. Families were ascertained to contain multiple myopic individuals. Genotyping of 1,367 SNPs was carried out within a 35cM (~23.9 Mb) candidate QTL region on 1p34-36. Multipoint variance components (VC) and regression-based (REG) linkage analyses were carried out separately in OOA and ASHK groups, and in a combined analysis that included all families.
Evidence of linkage of refractive error was found in both OOA (VC LOD=3.45, REG LOD=3.38 at ~59 cM) and ASHK families (VC LOD=3.12, REG LOD=4.263 at ~66 cM). Combined analyses showed three highly significant linkage peaks, separated by ~11cM (or 10 Mb), within the candidate region.
In a fine-mapping linkage study of OOA and ASHK families, we have confirmed linkage of refractive error to a QTL on 1p. The area of linkage has been narrowed down to a gene-rich region at 1p34.2-35.1 containing ~124 genes.