|
|
||||||||
1From the School of Optometry and Vision Sciences, Cardiff University, Cardiff, Wales, United Kingdom; the 2Department of Psychological Medicine, University of Wales College of Medicine, Cardiff, Wales, United Kingdom; and the 3Department of Genetics and Biometry, Roslin Institute, Roslin, United Kingdom.
| Abstract |
|---|
|
|
|---|
METHODS. The families comprised 245 subjects with phenotypic information and DNA available, of whom 170 were classified as affected. Subjects were genotyped for microsatellite markers spanning
40cM regions on 18p (MYP2), 12q (MYP3) and 17q, together with markers flanking COL2A1, COL11A1, and FBN1. Two-point linkage analyses were performed using the same disease gene segregation model as was used in the original publications, followed by nonparametric and multipoint analyses using Genehunter (http://linkage.rockefeller.edu/soft/gh/ provided in the public domain by Rockefeller University, New York, NY), with additional maximization over the parameter
, the proportion of linked families.
RESULTS. Evidence of linkage was found for the MYP3 locus on 12q (two-point Zmax = 2.54, P = 0.0003 and multipoint hLOD = 1.08 at
= 0.24, P = 0.023 for marker D12S332; nonparametric linkage [NPL] = 1.49, P = 0.07 for marker D12S1607). For the 17q locus there was weak evidence of excess allele sharing and linkage under a recessive model (NPL = 1.34, P = 0.09 for marker D17S956; two-point hLOD = 1.24 at
= 0.30 for marker D17S1795; multipoint hLOD = 1.24 at
= 0.17, P = 0.014 for marker at 77.68 cM, between markers D17S956 and D17S1853). No significant linkage was found to the MYP2 locus on 18p, or to the COL2A1, COL11A1, and FBN1 genes.
CONCLUSIONS. These results suggest that the MYP3 locus on 12q could be responsible for high myopia in approximately 25% of the U.K. families showing apparent autosomal dominant transmission, but that the loci on 18p and 17q are less common causes. Thus, additional loci for high myopia are likely to be the cause of the majority of cases of high myopia in the United Kingdom.
Dichotomizing on the basis of a refractive error of 6.00 D of myopia in the most minus meridian of the right eye (RE MMM), Young et al.8 identified a locus for fully penetrant high myopia in a large family of Chinese descent (two-point LOD score 3.90 for marker D18S1138 at 18p11.31). Six of the other seven families described in this report were also consistent with linkage to this locus, which has since been named MYP2. A follow-up study in Hong Kong10 reported evidence of linkage to this locus in 5 of the 15 participating families, and an association study11 showed that the MYP2 locus was an important contributor to the inheritance of high myopia in an isolated community in Sardinia. A second locus9 for fully penetrant nonsyndromic high myopia, MYP3 on 12q, was found in a single, large, German-Italian family (two-point LOD score 3.85 for markers D12S1706 and D12S327 at 12q21-23). More recently, a third locus has been identified on 17q in a single large family of English-Canadian descent (two-point LOD score 3.17 for marker D17S1604 on 17q21-22, dichotomization criterion 5.00 D RE MMM)7 and there has been suggestive evidence for a fourth locus on 7q in a study of 23 small families in France.6 Linkage analysis in 53 nuclear families of juvenile-onset myopes by Mutti et al.12 found no evidence of involvement of the MYP2 or MYP3 loci (dichotomization criterion 0.75 D RE least minus meridian with onset before 17 years of age).
Therefore, it is known from molecular genetic studies that at least three high-penetrance loci exist for autosomal dominant high myopia. However, it is not clear to what extent these genes of large effect are responsible for causing high myopia in the general population. To address this question, we assessed the extent of genetic linkage to the 18p, 12q, and 17q loci in 51 U.K. pedigrees in which at least two siblings had myopia of greater than 6.00 D. These results have been presented in abstract form (Farbrother JE, et al. IOVS 2003;44:ARVO E-Abstract 4780).
| Methods |
|---|
|
|
|---|
Pedigrees containing at least two siblings with myopia of
6.00 D in the least-minus meridian of both eyes were recruited for the study with the help of optometrists distributed across the United Kingdom and ophthalmologists at selected U.K. hospitals. (Seven of the families recruited from optometric practice were ascertained as part of a contemporaneous investigation into the genetic epidemiology of high myopia13 .) All families meeting these criteria who agreed to participate in the study were included. Six of the 51 families were bilineal (that is, high myopia was present in more than one of the unrelated founder or married-in individuals). These families were not excluded, since simulation studies suggest that bilineality has little effect on the mean maximum LOD score and the estimate of the recombination fraction.14 The ocular history and current spectacle or contact lens prescription was obtained from each subjects optometrist. For subjects who had undergone cataract removal, presurgery spectacle prescriptions were obtained. The family of any subject with high myopia as part of a known syndrome or associated with a systemic or ocular condition was excluded from the study. The 51 families comprised 306 individuals, with phenotypic information available for 254. According to the criteria used by Young et al.,8 9 173 subjects were classed as affected. DNA was available for 245 subjects, of whom 170 were affected. The subjects phenotype information is presented in Supplemental Table S1 at www.iovs.org/cgi/content/full/45/9/2879/DC1. Pedigree diagrams are presented in Supplemental Figure S1 at www.iovs.org/cgi/content/full/45/9/2879/DC1.
DNA was obtained by postal collection of saline mouthwashes.15 Microsatellite repeat sequences were selected for genotyping on the basis of their genetic location within and adjacent to the previously linked regions on 18p, 12q, and 17q.7 8 9 In addition, markers flanking the genes at which causal mutations have been identified for Stickler syndrome type I (COL2A1),16 Stickler syndrome type II (COL11A1),17 and Marfan syndrome (FBN1)18 were genotyped to exclude linkage to these regions, under the assumption that these syndromes are the most prevalent systemic causes of high myopia.19 Physical locations were established for all but one of the markers.
For genotyping, PCR reaction volumes of 12 µL comprised 1x PCR buffer (Promega, Madison, WI), 200 µM each dNTP (Amersham-Pharmacia, Piscataway, NJ), 1.5 mM MgCl2, 0.6 or 0.3 µM of each primer (Sigma-Genosys, Woodland, TX), 0.025 U Taq DNA polymerase (Promega) and
14 ng template genomic DNA. One primer was 5' labeled with a fluorescent dye (FAM, TET, or HEX). Microsatellites were PCR amplified individually, multiplexed along with TAMRA 350 size standards (Perkin Elmer, Boston, MA) and electrophoresed on a 10% denaturing polyacrylamide gel (LongRanger acrylamide; BMA Bioproducts, Lowell, MA) on an automated sequencer and analyzed with a model 373 DNA sequencer running Genescan 3.0 software (Applied Biosystems, Foster City, CA). A sample of DNA pooled from 83 control subjects was included on each plate to standardize the allele size calling between gels. Alleles were sized manually using the peak labeling function of the software program (Genotyper 1.1; ABI).
Allele frequencies were calculated from the founder and married-in individuals (n = 87). Parametric, two-point linkage analysis was performed using MLINK.20 21 The disease gene segregation models used by Young et al.8 9 22 were used for the parametric analyses. For the regions on chromosome 18p and 12q, this consisted of a 100% penetrant, dominant disease gene with frequency 0.0133 and no phenocopies. For the region on 17q and the Marfan and Stickler loci, this consisted of a 90% penetrant, dominant disease gene, with a frequency of 0.0133 and no phenocopies. The presence of heterogeneity in the participating families was considered likely, and hence the A test23 24 was performed using the Table program to maximize the two-point LOD scores over the additional parameter
, the proportion of linked families. GeneHunter 2.125 was used to perform nonparametric linkage analyses and calculate multipoint LOD scores. The data for all markers were checked for Mendelian errors with PedCheck,26 and unlikely double recombinants were checked for maximum-likelihood haplotypes generated by Genehunter 2.1. Statistical significance for the multipoint heterogeneity LOD scores and power were estimated from analysis of simulated marker data generated for the pedigrees by SLINK27 28 for an eight-allele marker with equal allele frequencies at
= 0.0 and
= 0.00, respectively.
After the planned linkage analyses had been performed (i.e., using linkage models that replicated those used for the discovery of the 18p, 12q, and 17q loci), a series of exploratory linkage analyses were undertaken to examine the effect of analyzing refractive error as a quantitative trait. Refractive error in the least minus meridian (RE LMM) was chosen as the quantitative trait of interest, rather than mean spherical equivalent (MSE) or RE MMM, to avoid potential bias from cylinder powers.29 Refractions were transformed using a Box-Cox type function30 {RETRANS = sign(RE LMM) x [(|RE LMM| + 1)0.415 1] x 5} which has the property of removing the kurtosis and most of the skew from the refractive error distribution of the general population. Linkage analysis was performed using Merlin31 (for standard quantitative trait locus [QTL] or variance components analysis) or Merlin-Regress, which implements an extension of the Haseman-Elston procedure.32 The correlation in refractive error between eyes was rs = 0.95 (P < 0.0001). Quantitative analyses repeated with the average refractive error in the two eyes used as the quantitative trait instead of that in just the right eye produced very similar results.
| Results |
|---|
|
|
|---|
|
|
|
0.5 would be expected to occur 87% of the time if 60% of the families were linked and 68% of the time if 40% of the families were linked.
MYP3 Locus on 12q
For the MYP3 locus, the maximum two-point LOD score (Zmax) was 2.54 at
= 0.15 for D12S332 (asymptotic P = 0.0003; Table 2 ). The multipoint analysis results in this region were lower than the two-point results, but remained positive. The maximum multipoint hLOD was 1.08, for
= 24% for marker D12S332 (2-df, asymptotic P = 0.041; empiric estimate P = 0.023). The maximum multipoint nonparametric linkage (NPL) was 1.49 for marker D12S1607 (P = 0.07). Haplotype assessment demonstrated that 21 of the 42 families were consistent with linkage, but of these only one pedigree (number 24) had a multipoint LOD exceeding 1.0. Recombination events in the families defined a linked region between the markers D12S348 and D12S1605.
|
= 0.30 for D17S1795, maximum multipoint hLOD = 1.24 between markers D17S956 and D17S1853,
= 0.17; 2-df, asymptotic P = 0.029; empiric estimate P = 0.014).
|
|
Analysis of Myopia as a Quantitative Trait
After assessing the evidence for linkage to the three known myopia loci with linkage models that replicated those used to identify the loci originally, we performed a series of exploratory linkage analyses in which refractive error was examined as a quantitative trait. There was no evidence of linkage to any of the three loci, when using variance components analysis (Merlin-VC) or an extension of the Haseman-Elston regression method (Merlin-Regress). Using the QTL analysis method implemented in Merlin, there was no evidence of linkage at the 18p (Zmax = 0.76 at 0 cM, P = 0.20) or 12q (Zmax = 1.10 at 111.87 cM, P = 0.14) loci, but there was suggestive evidence of linkage at the 17q locus (Zmax = 2.42 at 68.44 cM, P = 0.008). Because of the exploratory nature of these quantitative analyses and the multiple testing that this entailed, we consider the absolute values of these Z-scores and their associated probabilities to be less meaningful than their relative levels.
| Discussion |
|---|
|
|
|---|
3 in the absence of locus heterogeneity was 100% in all models. However, the power to detect a locus by using linkage analysis inevitably decreased as the proportion of linked families (
) decreased. For our families, the power to detect an hLOD
3 exceeded 80% when 80% of families were linked, but decreased to 20% to 25% when only 40% of families were linked, depending on the exact linkage model used. Only five families (pedigrees 24, 30, 32, 33, and 50) were of sufficient size to produce a LOD score
1 in isolation. Hence, although consistency with linkage can be described in the case of individual families, this is not conclusive evidence of inheritance due to the locus in question.
MYP2 Locus on 18p
Of the 42 families consistent with fully penetrant autosomal dominant transmission, linkage to MYP2 (telomeric to D18S52)33 could be excluded in 21, in which maximum multipoint LOD scores were negative and haplotype analysis was inconsistent with linkage. Coupled with the absence of significant linkage and excess allele-sharing across the entire data set, the findings strongly suggest that the MYP2 locus is not responsible for most cases of fully penetrant autosomal dominant high myopia in the United Kingdom. Because seven of the eight pedigrees examined by Young et al.8 in their original report were linked to MYP2, it could be that the etiology of high myopia in the United Kingdom is different from that in the United States. Alternatively, the high proportion of pedigrees linked to the MYP2 locus in the cohort studied by Young et al. could have occurred by chance. Consequently, fully penetrant disease allele(s) at the MYP2 locus may be rarer in many population groups than originally seemed the case.
MYP3 Locus on 12q
The multipoint linkage results for the MYP3 locus just reached statistical significance when considered in isolation (P = 0.023). However, once multiple testing is taken into account, the possibility that this result is a false positive becomes more likely. If seeking to replicate linkage to three separate loci is considered as performing three independent statistical tests, the results did not reach statistical significance (empiric estimate P = 0.07). If the weak linkage signal is not a false-positive result, then this study suggests that the MYP3 locus could be a relatively common cause of familial high myopia in the United Kingdom. Our estimate of the proportion of linked families,
= 0.24 (95% confidence interval; 0.020.48), suggests that approximately one quarter of seemingly autosomal dominant families are linked to the MYP3 locus. The 95% confidence interval suggests that it is unlikely that either none of the families, or more than half of the families are linked to this locus. However, even in view of the wide confidence interval, our estimate of the proportion of linked families must be considered tentative, not least since the accuracy of the maximum likelihood estimate of
is dependent on the accuracy of the linkage model chosen (and this is subject to inevitable uncertainty regarding the criteria used to define affectation status, penetrance, and disease gene frequency) as well as the characteristics of the alternative unlinked loci.34 35 Assessment of haplotypes demonstrated that 21 of the 42 families were consistent with linkage, but of these only one pedigree (number 24) had a multipoint LOD exceeding 1.0. Recombination events in these families suggest the disease gene is harbored within the narrowed interval (
15 cM) between markers D12S348 and D12S1605.
The difference between the two-point and multipoint results for MYP3 is explained in part by the extra identity-by-descent (IBD) information available in the analysis of multiple markers. In pedigree 33 in particular, additional data for the centromeric markers reduced the likelihood that all the affected individuals shared the region between D12S348 and D12S1605. Lower multipoint results could also have occurred as a consequence of the constraint of
, in combination with heterogeneity in the data set.36 The data were rigorously checked for genotyping and marker map distance errors to exclude these as a potential cause of the two-point versus multipoint LOD score disparity.
None of the myopes among the 80 unaffected subjects included in the linkage analysis were younger than 18 years at the time of the eye examination on which their classification was based, and, as such, their status could be determined with relative confidence37 (it should be noted that 5 of these 80 unaffected subjects were not genotyped). However, because of the potential for the measurement of spectacle refractions to vary from optometrist to optometrist, the effect of changing the affectation threshold was tested. Increasing the stringency of the threshold by 0.50 D increased the multipoint hLOD by 49% and the NPL statistic by 32%, whereas reducing the stringency by 0.50 D reduced the multipoint hLOD by 15% but increased the NPL statistic by 5%. However, the maxima remained close to marker D12S332 and the proportion of linked families stayed between 20% and 30%.
17q Locus
The suggestion of excess allele sharing found using nonparametric analysis, even though linkage was absent in the original, dominant parametric model,22 prompted us to investigate linkage under a recessive model for this locus. When considered in isolation, the evidence for linkage with this recessive model reached statistical significance (P = 0.014) and explained the excess allele sharing. After accounting for multiple testing, the result would be of only borderline statistical significance. Yet the fact that, under a quantitative trait model, the linkage signal for this locus was the highest of the three suggests that it could either be the cause of the high myopia in a small number of the families studied or could act as a susceptibility gene in a larger proportion of them. It is conceivable that allelic variants in the same gene at the locus on 17q reported by Paluru et al.7 and Young et al.22 could give rise to either dominant or recessive inheritance in different families. Alternatively, our results would be consistent with the presence of a recessive or additive locus for high myopia close to the previously identified high penetrance dominant locus.
Analysis of Myopia as a Quantitative Trait
The highly selective ascertainment scheme adopted in this study meant that the degree of phenotypic variability in refractive error in these pedigrees was limited. This meant that only one of the three QTL analysis methods tested, the nonparametric QTL analysis option of Merlin, was strictly appropriate for assessing linkage. (Both the variance components and the regression-based analysis methods are designed for traits with a normal distribution and may be sensitive to deviations from normality. In addition, variance components analysis is better suited to nonselectively ascertained cohorts or cohorts ascertained using selection from both tails of the distribution, rather than the one-sided selection scheme used in the present study.38 Finally, the regression analysis approach is highly sensitive to misspecification of trait parameters.32 ) Using Merlin-QTL, all three loci showed at least a weak linkage signal, with that from the 17q locus being the strongest. We interpret these quantitative analyses as backing up the results of the main linkage analyses, in suggesting that whereas the three loci examined may contribute to high myopia in these U.K. pedigrees, they must either act only in a minority of families, or harbor susceptibility genes, rather than high-penetrance myopia genes.
Implications for Future Linkage Studies
The 51 U.K. pedigrees examined herein were ascertained on the basis of an affected sibling pair, which would be expected to enrich the sample for high penetrance autosomal dominant disease genes. Yet the three high myopia loci for which convincing evidence of linkage has so far been obtained, all of which show apparent high-penetrance autosomal dominant transmission, appeared to be the cause of high myopia in only a minority of cases. This suggests either that additional high myopia loci remain to be discovered and that together these loci account for much of the high myopia in the general population or that most cases of high myopia are not genetic in origin. Our genetic epidemiology investigation of high myopia in the United Kingdom13 reached the conclusion that high-penetrance dominant inheritance of high myopia is the exception rather than the rule. Thus, finding additional loci may be difficult, requiring either the fortuitous discovery of further large families suitable for parametric linkage analysis (especially consanguineous pedigrees for mapping recessive high myopia genes) or very large collections of small pedigrees suitable for mapping "common genes of small effect." The results from a recent linkage study in a clinically unselected sample of dizygotic twins in which refractive error was analyzed as a quantitative trait were very encouraging (Hammond, et al. IOVS 2003;44:ARVO E-Abstract 2015). Such approaches should be complementary to those using extreme proband selection schemes like the affected sib-pair design, to determine whether the same genes or pathways underlie both high and low myopia.
| Conclusions |
|---|
|
|
|---|
| Electronic Database Information |
|---|
|
|
|---|
| Note Added in Proof |
|---|
|
|
|---|
| Acknowledgements |
|---|
| Footnotes |
|---|
Submitted for publication October 21, 2003; revised April 13, 2004; accepted May 10, 2004.
Disclosure: J.E. Farbrother, None; G. Kirov, None; M.J. Owen, None; R. Pong-Wong, None; C.S. Haley, None; J.A. Guggenheim, None
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked "advertisement" in accordance with 18 U.S.C.
1734 solely to indicate this fact.
Corresponding author: Jeremy A. Guggenheim, School of Optometry and Vision Sciences, Cardiff University, Redwood Building, King Edward VII Avenue, Cardiff CF10 3NB, UK; guggenheim{at}cardiff.ac.uk.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. T. Pardue, A. E. Faulkner, A. Fernandes, H. Yin, F. Schaeffel, R. W. Williams, N. Pozdeyev, and P. M. Iuvone High Susceptibility to Experimental Myopia in a Mouse Model with a Retinal ON Pathway Defect Invest. Ophthalmol. Vis. Sci., February 1, 2008; 49(2): 706 - 712. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. K. Pertile, M. Schache, F. M. A. Islam, C. Y. Chen, M. Dirani, P. Mitchell, and P. N. Baird Assessment of TGIF as a Candidate Gene for Myopia Invest. Ophthalmol. Vis. Sci., January 1, 2008; 49(1): 49 - 54. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Schache, A. J. Richardson, K. K. Pertile, M. Dirani, K. Scurrah, and P. N. Baird Genetic Mapping of Myopia Susceptibility Loci Invest. Ophthalmol. Vis. Sci., November 1, 2007; 48(11): 4924 - 4929. [Abstract] [Full Text] [PDF] |
||||
![]() |
Q Zhang, X Guo, X Xiao, X Jia, S Li, and J F Hejtmancik Novel locus for X linked recessive high myopia maps to Xq23-q25 but outside MYP1. J. Med. Genet., May 1, 2006; 43(5): e20 - e20. [Abstract] [Full Text] [PDF] |
||||
![]() |
X.-J. Qin, T. H. Margrain, C. H. To, N. Bromham, and J. A. Guggenheim Anisometropia Is Independently Associated with Both Spherical and Cylindrical Ametropia Invest. Ophthalmol. Vis. Sci., November 1, 2005; 46(11): 4024 - 4031. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |