# Linkage analysis of high myopia susceptibility locus in 26 families

Sandrine Paget, Sophie Julia, Zulma Vitezica, Vincent Soler, François Malecaze, Patrick Calvas

► To cite this version:

Sandrine Paget, Sophie Julia, Zulma Vitezica, Vincent Soler, François Malecaze, et al.. Linkage analysis of high myopia susceptibility locus in 26 families. Molecular Vision, 2008, 14 (291-95), pp.2566-2574.
hal-02663955

HAL Id: hal-02663955

https://hal.inrae.fr/hal-02663955

Submitted on 31 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Linkage analysis of high myopia susceptibility locus in 26 families

Sandrine Paget,<sup>1,2</sup> Sophie Julia,<sup>2,3</sup> Zulma G. Vitezica,<sup>1,2</sup> Vincent Soler,<sup>1,4</sup> François Malecaze,<sup>1,2,4</sup>
Patrick Calvas<sup>1,2,3</sup>

<sup>1</sup>Inserm, U563, Centre de Physiopathologie de Toulouse Purpan, Toulouse, France; <sup>2</sup>Université Toulouse III Paul-Sabatier, UMRs563, Toulouse, France; <sup>3</sup>CHU de Toulouse, Hôpital Purpan, Service de Génétique Médicale, Toulouse, France; <sup>4</sup>CHU de Toulouse, Hôpital Purpan, Service d'Ophtalmologie, Toulouse, France

**Purpose:** We conducted a linkage analysis in high myopia families to replicate suggestive results from chromosome 7q36 using a model of autosomal dominant inheritance and genetic heterogeneity. We also performed a genome-wide scan to identify novel loci.

**Methods:** Twenty-six families, with at least two high-myopic subjects (ie. refractive value in the less affected eye of -5 diopters) in each family, were included. Phenotypic examination included standard autorefractometry, ultrasonographic eye length measurement, and clinical confirmation of the non-syndromic character of the refractive disorder. Nine families were collected de novo including 136 available members of whom 34 were highly myopic subjects. Twenty new subjects were added in 5 of the 17 remaining families. A total of 233 subjects were submitted to a genome scan using ABI linkage mapping set LMSv2-MD-10, additional markers in all regions where preliminary LOD scores were greater than 1.5 were used. Multipoint parametric and non-parametric analyses were conducted with the software packages Genehunter 2.0 and Merlin 1.0.1. Two autosomal recessive, two autosomal dominant, and four autosomal additive models were used in the parametric linkage analyses.

**Results:** No linkage was found using the subset of nine newly collected families. Study of the entire population of 26 families with a parametric model did not yield a significant LOD score (>3), even for the previously suggestive locus on 7q36. A non-parametric model demonstrated significant linkage to chromosome 7p15 in the entire population (Z-NPL=4.07, p=0.00002). The interval is 7.81 centiMorgans (cM) between markers D7S2458 and D7S2515.

**Conclusions:** The significant interval reported here needs confirmation in other cohorts. Among possible susceptibility genes in the interval, certain candidates are likely to be involved in eye growth and development.

Myopia is currently divided into low to moderate myopia (refractive values between -0.5 and -5 diopters [D]) and high myopia (beyond -5 D). The prevalence of myopia varies moderately in Western countries, ranging from 16% in Australia and 18% in the Netherlands, to an average value of 25% in the USA in adults between 40 and 80 years old [1]. The Asian population seems to be more affected than Western populations. The prevalence of myopia ranges from 16% in Australia and 18% in the Netherlands to 25% in the United States in adults aged 40-80 years [1] to much higher values in Eastern Asian countries. Over 38% of urban Singaporean Chinese adults [2] and up to 80% of teenagers (16-18 years old) in urban Taiwan [3] are affected. For high myopia, the prevalence is 4.5% in populations of Western European origin [1] as compared to the 8%-9% [2,4] observed in Eastern Asian adults over the age of 40. Laser refractive surgery as a myopia-related cost was estimated to be 4.6 billion dollars for the United States alone in 1990 [5]. Stambolian et al. [6] estimated that by 2005 this cost had doubled.

Increasing prevalence and associated health costs [5,6] make myopia an important public health problem [7,8]. High myopia predisposes patients to premature cataracts and an increased risk for retinal detachment, glaucoma, macular degeneration, and blindness [9]. It is one of the major causes of legal blindness worldwide [10,11]. No one certain cause of myopia has yet been identified. Several family and twin studies have shown the role of genetics in the etiology of myopia [12-16]. Different inheritance models have been proposed [14,17]. It is mostly thought that myopia results from interactions between genetic and environmental factors. For example, close vision and near-work activities increase myopia prevalence (for review, see Gilmartin [18] and Morgan and Rose [19]). Different chromosomal localizations for high myopia have been reported in linkage analyses over the last decade. To date, syndromic and isolated forms of X-linked high myopia have been respectively associated with Xq28 [20] and Xq23-25 [21,22]. Larger autosomal regions have been reported for non-syndromic high myopia. Some were identified in single large families [23-26] or in subsets of families [27-30]. Low and moderate myopia has also been linked to different regions [31-34]. Candidate genes have also been suggested [35-39]. For instance, the 5' region of *HGF* (hepatocyte growth factor) in mouse [40] and man may

<footnote>Correspondence to: Patrick Calvas, INSERM, U563 - CPTP, Bâtiment B, BP 3028, Place du Dr Baylac, 31024, Toulouse cedex 3, France; Phone: +33(0)5 61 77 90 79; FAX: +33(0)5 61 77 90 73; email: calvas.p@chu-toulouse.fr</footnote>
contain a polymorphism associated with early-onset extreme myopia (beyond −10 D) in the Han Chinese population [41].

Our previous studies analyzed the inheritance of high myopia, suggested an autosomal dominant model with weak penetrance [42], and found a suggestive linkage to high myopia on chromosome 7q36 [28]. To confirm the previous linkage peak and to look for other loci involved in high myopia, a novel study was conducted in an extended population. We were able to collect nine additional families. In this study, our goal is to perform a genome-wide scan (GWS) on these nine new families as well as the combined linkage analysis on the full sample of families.

## METHODS

**Subjects:** Volunteers received information about the study in agreement with the Helsinki Declaration principle. Participants were included according to the French laws governing participation in biomedical research and with approval of the local ethics committee. Each participant provided written, informed consent. Families of high myopia were enrolled through a proband affected with non-syndromic high myopia (refractive value, $RV \le -5$ D, axial length>26 mm). The family was retained for the study if it included at least two high myopic members and if at least three members on two generations agreed to participate. Exclusion criteria were myopia in prematurity retinopathy, Marfan's syndrome, Stickler's syndrome, Wagner's disease, retinitis pigmentosa, corneal dystrophy, keratoconus, myopia secondary to cataract, or any developmental genetic disorder. In addition, families with unilateral high myopia and families with X-linked compatible high myopia were excluded.

Twenty-six families were included in this study (nine new families and 17 previously studied families [28]). The new families (16, 18, 20–26) were constituted by 136 persons (34 high myopic subjects). Among the 17 previously studied families, five (1, 4, 10, 12, and 14) had a total of 20 new participants with one to eleven additional individuals per family. All the members in these families under 20 years of age at the time of the first evaluation [28] were phenotypically re-evaluated, and it was found that a total of five persons in four families (individual 15 [ID15] in family 5 [F5], ID4 in F11, ID6 in F17, ID4 and ID7 in F19) who had moderate myopia became high myopic. Age at recruitment was between 5 and 95 years old.

**Clinical examination:** Refractive values (RVs) were measured using standard aut refractometry after dilation. In patients under 16 years old, three instillations of 1% cyclopentolate at 0 min, 5 min, and 10 min were performed and RVs were measured 45–60 min after the last instillation. In patients over 16 years old, dilation was induced with 1% tropicamide at 0 min, 5 min, and 10 min, and the RV measured 30–45 min later. Tropicamide (1%) was instilled every 5 min three times, and then measurement is performed. Axial length

was evaluated by A-scan ultrasonography. A total of three readings were taken for each eye, and the average value was recorded. Astigmatism was assessed by autorefractometer. The included patients had minimal astigmatism, indicating that bias was avoided. Each proband had a complete ocular examination including visual acuity, intraocular pressure, and fundus, to confirm that the high myopia was a primary finding and in completion of a general physical examination to avoid recruitment of syndromic myopia patients. For the probands' relatives, ophthalmologists of their choice were asked to examine subjects according to the standard protocol above. In the case of cataract surgery, only the ocular measurement before surgery was considered.

**Genotyping:** DNA extraction from venous blood was performed according to standard phenol-chloroform extraction procedure [43]. For genotyping, four DNA size references were added to standardize size-calling between runs according to the manufacturer's procedure (Applied Biosystems, Foster City, CA). The average spacing distance between markers was 9.1 centiMorgan (cM), and the average heterozygosity was 0.78. The reduction of the detected chromosomal locus was made by genotyping of highly heterozygous microsatellites as per the UniSTS (NCBI) and GDB databases.

Amplifications of microsatellites were performed in 10 µl total volume, with 1X GoTaq PCR buffer (Promega, Madison, WI), 2.0–2.5 mM MgCl₂, 200 µM each dNTP, 250 nM specific primer to each marker, 0.28 U/µl GoTaq DNA polymerase (Promega), and 25 ng of genomic DNA. Fluorescent labeled amplification products were electrophoresed on an ABI PRISM® 3100 Genetic Analyzer with GeneScan 500HD ROX standard size label (Applied Biosystems) and analyzed with Genescan 3.5 software (Applied Biosystems). Alleles were analyzed with Genotyper 3.6 program (Applied Biosystems).

**Statistical analyses:** Familial relationship inconsistencies and genotyping error checking and cleaning were performed with the Merlin 1.0.1 program [44,45]. Marker allele frequencies were estimated from the founders of the pedigrees. Sex-average genetic distances were taken from the Marshfield Center for Medical Genetics. The order consistency was compared to the NCBI physical map (Build 36.2). The effects of ethnic origin, sex, and age were compared between affected and unaffected individuals using a Fisher's exact and student's t-tests with a level of significance below the p value of 0.01. In this study, RV is a spherical equivalent (the correlation between RV and spherical equivalent was 0.76 and 0.85 for the right and the left eye, respectively). The RV correlation between right and left eye was equal to 0.95. The RV mean was used for status qualification and as the trait in quantitative trait loci (QTL) analyses. Subjects with a RV mean of less than or equal to −5 D were considered to be affected. The others persons were classified as unaffected.
<table><caption>TABLE 1. PARAMETRIC MODELS USED IN THE PARAMETRIC MULTIPOINT GENOME-WIDE LINKAGE ANALYSIS.</caption><thead><tr><th>Phenocopy rates and penetrances</th><th>DD</th><th>Dd</th><th>dd</th><th>Model</th></tr></thead><tbody><tr><td rowspan="2">Autosomal recessive models</td><td>0</td><td>0</td><td>0.58</td><td>Model 1</td></tr><tr><td>0</td><td>0</td><td>0.9</td><td>Model 2</td></tr><tr><td rowspan="2">Autosomal dominant models</td><td>0</td><td>0.58</td><td>0.58</td><td>Model 3</td></tr><tr><td>0</td><td>0.9</td><td>0.9</td><td>Model 4</td></tr><tr><td rowspan="4">Autosomal additive models</td><td>0.1</td><td>0.58</td><td>0.58</td><td>Model 5</td></tr><tr><td>0.1</td><td>0.9</td><td>0.9</td><td>Model 6</td></tr><tr><td>0.2</td><td>0.58</td><td>0.58</td><td>Model 7</td></tr><tr><td>0.2</td><td>0.9</td><td>0.9</td><td>Model 8</td></tr></tbody></table>

Based on the postulation of a single, two-allele gene in which “d” would be the disease-causing allele, eight parametric models have been tested. Models one and two assume an autosomal recessive mode of inheritance. Models three and four assume an autosomal dominant mode of inheritance. Models six to eight assume autosomal additive transmission as per Chen et al. [30]. In all models, penetrances of the genotype of 0.58 and 0.90 were used according to those reported respectively by Naiglin et al. [28] or Young et al. [23,27]. In addition, models five to eight consider ten or twenty percent phenocopy rates.

<table><caption>TABLE 2. DEMOGRAPHIC CHARACTERISTICS.</caption><thead><tr><th>Demographic characteristics</th><th></th></tr></thead><tbody><tr><td>Number of families analyzed</td><td>26</td></tr><tr><td>Number of individuals</td><td>347</td></tr><tr><td>Total number of affected individuals</td><td>98</td></tr><tr><td>Total number of low and moderate myopia</td><td>47</td></tr><tr><td>Total number of individuals genotyped</td><td>233</td></tr><tr><td>Average number of generations (range)</td><td>3.35 (2 to 4)</td></tr><tr><td>Average number of individuals per family (range)</td><td>13 (5 to 24)</td></tr><tr><td>Average number of affected individuals per family (range)</td><td>4 (2 to 10)</td></tr><tr><td>Average number of genotyped individuals per family (range)</td><td>9 (3 to 22)</td></tr><tr><td>Average age, in years, of examined individuals (range)</td><td>37 (5 to 95)</td></tr><tr><td>Average spherical equivalence (range)</td><td>-4.84±6.19 (-25 to 4.5)</td></tr></tbody></table>

A summary of the population characteristics is given including average data of the ophthalmological evaluations. Affected subjects are high myopes with refractive errors beyond -5 diopters (D). Low and moderate myopes have refractive errors between -0.5 D and -5 D. The mean values and the range of age variation (in years) and spherical equivalent (in diopters) are given

Multipoint parametric linkage (PL) and non-parametric linkage (NPL) analyses were performed under multiple models with Genehunter 2.0 [46] and Merlin 1.0.1 programs [44]. A total of eight parametric models with different inheritance model, genotype penetrance, and phenocopy rates were tested with a frequency of the susceptibility allele of 0.013 (Table 1). For each model, the hypothesis of genetic heterogeneity was also tested, and Heterogeneity Log of the Odd (HLOD) was calculated.

Myopia was also treated as a quantitative trait. The RV mean was transformed to reach normality [47]. A QTL analysis was conducted with Merlin 1.0.1 using the regression option. The sex, age, and ethnic origin were included as covariates.

## RESULTS

*Population characteristics:* A total of 26 families and their DNAs were collected from all over France. Twenty-five families were of French origin and one family (of 13 members) was of Algerian origin. The entire cohort included 347 individuals. The population characteristics are shown in Table 2 and Figure 1.

Clinical data were collected for the 233 genotyped individuals. Among them, 45% were men and 55% women. Age distribution was similar in both genders (78% of men and 75% of women were older than 20 years). Eighty-one adults over 20 years old (32 men and 49 women) and 17 patients under 20 years old (6 boys and 11 girls) had high myopia. The unaffected group was formed by 97 adults over the age of 20 years (49 men and 48 women) and 38 individuals under the age of 20 years (16 boys and 14 girls). No significant effect
<figure><img src="image_3.png" /><figcaption>Figure 1. Genealogical trees of the 26 pedigrees. Black squares and circles denote subjects affected with high myopia (refractive value [RV] beyond -5 diopters [D]). Half-black squares and circles denote individuals with a RV between -0.5 and -5 D. The asterisks denote genotyped individuals.</figcaption></figure>
of sex (p=0.02), age (p=0.03), or ethnic origin (p=1) between
affected and unaffected persons was observed.

*Linkage analysis of the nine new families:* No significant or suggestive linkage was observed in parametric linkage (PL) and non-parametric linkage (NPL) analyses of the nine novel families in any model tested. No significant (p-value less than or equal to 0.000049) or suggestive (p-value greater than or equal to 0.000048 but less than or equal to 0.0017) [48] linkage.

*Linkage analyses in 26 families:* The PL analysis did not show any suggestive or significant linkage to any locus. Results of the NPL analysis are summarized in Table 3 and Figure 2. A significant NPL score of 3.74 (p=0.00002) and Z-NPL score of 4.07 (p=0.00002) was observed for markers D7S529 and D7S516 (Table 3 and Table 4). The region identified covers 7.81 cM between D7S2458 and D7S2515 (Table 4 and Figure 3). A weak signal was observed on chromosome 1p31, 6q15, and 9q21-q22 (Table 3).

*Definition of high myopia phenotype:* The definition of high myopia threshold is to some extent arbitrary. Thus, we examined the influence of the threshold used to define high myopia on the non-parametric results. Moderate myopic persons (with RV between -3 and -5 D) were considered affected in PL and NPL analysis. A suggestive linkage in 7p15.2-15.3 region was observed when the myopic threshold was equal to -6 D (p=0.0006), -3 D (p=0.00003), -2 D (p=0.00004), and -1.5 D (p=0.0006).

*QTL analysis on the 26 families:* The RV mean for the entire cohort was considered for a QTL analysis but did not give any significant or suggestive linkage. Merlin 1.0.1 gave an estimation of heritability of 10% for the RV mean. The axial

<figure><img src="image_4.png" /><figcaption>Figure 3. Non-parametric multipoint LOD score in 26 French families with myopia less than or equal to -5D for chromosome 7. LOD scores were plotted by Merlin 1.0.1 against marker distance given in centiMorgans (cM).</figcaption></figure>

length was only available for 122 subjects with clinical
examination, thus this trait was not used in the QTL analysis.

DISCUSSION

The aim of this study was to confirm previously reported
linkage positions. The results of the genome wide scan (GWS)
linkage study on the nine new families or on all 26 families
together do not support the suggestive linkage to 7q36 region.
The inclusion of new families and the enlargement of
previously recruited families considerably modified the
structure of the cohort. Although these data do not formally
exclude the linkage of a high myopia locus on chromosome
7q36, it is likely that other myopia susceptibility loci play a
more important role in the families that we studied. Our
inability to confirm the previously identified position
probably reflects the inclusion of more families with a larger
sample size and of additional markers. Genotyping errors in
the samples might vary between the study conducted by
Naiglin et al. [28] and the current study. Such errors can
influence linkage results by inflating recombination fraction
estimates in linkage analysis. Replication of GWS studies for
complex traits is difficult [48,49].

Another interesting finding is on chromosome 7p15 because linkage signal was detected in NPL analysis even when applying different definitions of the myopic phenotype. However, no significant linkage result in this region was found in any of the PL models used even under the hypothesis of genetic heterogeneity. Even if the NPL approach is robust in the face of uncertainty about the mode of inheritance [48, 49], the new interval we have identified at 7p15.2-p15.3 needs to be confirmed. However, the same locus has also been highlighted in a Ashkenazi Jewish cohort [32]. Klein et al. [50] have also found a nearby region at 7p21.

The qualification of quantitative RVs as myopia, emmetropia, or hyperopia is to some extent arbitrary. A quantitative trait analysis was performed. This approach is generally more informative and powerful because continuous data on refractive errors (e.g., RV or axial length) provide more information than discrete trait data (e.g., affected and unaffected). In our cohort, the QTL analysis did not confirm the NPL result. The weak signal of the NPL analyses of binary traits may be due to low power. This paper highlights the complex heredity of high myopia and the great difficulty of replicating mapping results.

Sixty-seven genes or putative coding sequences are contained in the significant 7p15 region detected in this study (about 7.81 cM long). Candidate genes could be similar in function or structure to genes in other loci intervening in high myopia. Neurotransmitters can modulate experimental myopia [51]. A neuropeptide gene (NPVF) lies in our candidate region on 7p15.2-p15.3. The molecule encoded by this gene belongs to the family of neuropeptides with the Arg-Phe-amide motif at their C-termini (RF-amide peptides)
<table><caption>TABLE 3. MAXIMUM NON-PARAMETRIC MULTIPOINT LINKAGE ANALYSIS RESULTS.</caption><thead><tr><th>Chromosome</th><th>Number of markers on chromosome</th><th>Marker at the maximum NPL score</th><th>Position (cM)</th><th>Z-NPL</th><th>p value</th><th>NPL</th><th>p value</th></tr></thead><tbody><tr><td>cytogenetic band</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>1p31</td><td>31</td><td>D1S218</td><td>95.31</td><td>1.74</td><td>0.04</td><td>1.38</td><td>0.006</td></tr><tr><td>2p24</td><td>30</td><td>D2S305</td><td>38.87</td><td>1.18</td><td>0.12</td><td>0.95</td><td>0.02</td></tr><tr><td>3p21.1</td><td>23</td><td>D3S1277</td><td>61.52</td><td>0.94</td><td>0.2</td><td>0.65</td><td>0.04</td></tr><tr><td>4q34.1</td><td>22</td><td>D4S1539</td><td>176.19</td><td>1.30</td><td>0.1</td><td>0.46</td><td>0.07</td></tr><tr><td>5p15.1-p14.3</td><td>22</td><td>D5S416</td><td>28.76</td><td>1.84</td><td>0.03</td><td>1.11</td><td>0.012</td></tr><tr><td>6q15</td><td>33</td><td>D6S462</td><td>99.01</td><td>3.07</td><td>0.0011</td><td>1.37</td><td>0.006</td></tr><tr><td>7p15</td><td>49</td><td>D7S529</td><td>39.82</td><td>4.07</td><td>2.10<sup>-5</sup></td><td>3.74</td><td>2.10<sup>-5</sup></td></tr><tr><td>8p22</td><td>14</td><td>D8S549</td><td>31.73</td><td>0.83</td><td>0.2</td><td>0.46</td><td>0.07</td></tr><tr><td>9q21-q22</td><td>20</td><td>D9S287</td><td>103.42</td><td>2.32</td><td>0.01</td><td>1.20</td><td>0.009</td></tr><tr><td>10q26</td><td>20</td><td>D10S587</td><td>147.57</td><td>1.20</td><td>0.12</td><td>0.42</td><td>0.08</td></tr><tr><td>11p13</td><td>18</td><td>D11S935</td><td>45.94</td><td>1.21</td><td>0.11</td><td>0.29</td><td>0.13</td></tr><tr><td>12p13.2-q24.1</td><td>19</td><td>D12S83</td><td>75.17</td><td>0.31</td><td>0.4</td><td>0.03</td><td>0.4</td></tr><tr><td>13q14</td><td>14</td><td>D13S263</td><td>38.32</td><td>1.72</td><td>0.04</td><td>0.86</td><td>0.02</td></tr><tr><td>14q24.3</td><td>14</td><td>D14S74</td><td>87.36</td><td>0.41</td><td>0.3</td><td>0.13</td><td>0.2</td></tr><tr><td>15q23</td><td>14</td><td>D15S131</td><td>71.28</td><td>1.65</td><td>0.05</td><td>0.39</td><td>0.09</td></tr><tr><td>16q23.3</td><td>13</td><td>D16S3091</td><td>111.12</td><td>1.19</td><td>0.12</td><td>0.54</td><td>0.06</td></tr><tr><td>17q25.1</td><td>28</td><td>D17S1807</td><td>99.21</td><td>1.38</td><td>0.08</td><td>0.53</td><td>0.06</td></tr><tr><td>18q11.2-q12.1</td><td>14</td><td>D18S478</td><td>52.86</td><td>0.42</td><td>0.3</td><td>0.11</td><td>0.2</td></tr><tr><td>19q13.3</td><td>12</td><td>D19S902</td><td>72.72</td><td>0.63</td><td>0.3</td><td>0.09</td><td>0.3</td></tr><tr><td>20q13.13</td><td>13</td><td>D20S196</td><td>75.01</td><td>1.87</td><td>0.03</td><td>0.85</td><td>0.02</td></tr><tr><td>21q22</td><td>5</td><td>D21S266</td><td>45.87</td><td>0.13</td><td>0.4</td><td>0.01</td><td>0.4</td></tr><tr><td>Chromosome22</td><td>7</td><td>No positive hit</td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

An overview of the data obtained by the whole genome scan is provided with the number of markers used on each chromosome, marker chromosomal locations and their positions in centiMorgans (cM) on the genetic map. Markers were selected for each chromosome by the maximum non-parametric Z and LOD scores, which are associated with their p values.
<figure><img src="image_5.png" /><figcaption>Figure 2. Non-parametric multipoint linkage analysis results for the 22 autosomes. Non-parametric multipoint linkage analysis results across the 22 autosomes. Genome-wide scan LOD scores of the entire population have been plotted against chromosomal location.</figcaption></figure>

<table><caption>TABLE 4. NON-PARAMETRIC MULTIPOINT LINKAGE ANALYSIS FOR REGION 7p15.</caption><thead><tr><th>Marker</th><th>Z NPL</th><th>p value</th><th>NPL score</th><th>p value</th></tr></thead><tbody><tr><td>D7S493</td><td>2.97</td><td>0.0015</td><td>1.68</td><td>0.003</td></tr><tr><td>D7S2458</td><td>3.5</td><td>0.0002</td><td>2.6</td><td>3x10<sup>-4</sup>*</td></tr><tr><td>D7S629</td><td>3.53</td><td>0.0002</td><td>3.05</td><td>9x10<sup>-5</sup>*</td></tr><tr><td>D7S673</td><td>3.53</td><td>0.0002</td><td>3.06</td><td>9x10<sup>-5</sup>*</td></tr><tr><td>D7S529</td><td>4.07</td><td>0.00002</td><td>3.74</td><td>2x10<sup>-5</sup>**</td></tr><tr><td>D7S516</td><td>3.97</td><td>0.00004</td><td>3.59</td><td>2x10<sup>-5</sup>**</td></tr><tr><td>D7S2515</td><td>3.31</td><td>0.0005</td><td>2.9</td><td>1.3x10<sup>-4</sup>*</td></tr><tr><td>D7S2496</td><td>2.73</td><td>0.003</td><td>2.9</td><td>1.1x10<sup>-4</sup>*</td></tr><tr><td>D7S632</td><td>2.22</td><td>0.013</td><td>1.42</td><td>0.005</td></tr></tbody></table>

Individual non-parametric Z and LOD scores are given with their p-values, for markers in the 7p15 region linked to high myopia. The asterisk indicates suggestive values and the double asterisk indicates significant values.

involved in various neurotransmission/neuromodulation processes and muscle contraction control. It is expressed in the rat and human central nervous system [52] and in the retina [53]. Interestingly, it can inhibit GABAergic neurotransmission [54]. Gamma-aminobutyric acid (GABA) is an inhibitory neurotransmitter in adults but plays a neurotrophic role in the embryo. Stone et al. [54] have demonstrated that GABA receptor agonists and antagonists affect eye growth and modulate experimental myopia in chicken. The ability of NPVF to interact with the GABAergic system would be in accordance with both the strong effect of neurotransmitter modulation in experimental myopia and the role of heredity, which makes it a candidate gene for high myopia susceptibility. Similarly, other neurotransmitters that

can be involved in experimental myopia such as neuropeptide Y (NPY) are encoded by genes present within the interval [55]. Recently, Han et al. [41] demonstrated the association of a hepatocyte growth factor (HGF) polymorphism with high myopia. The overexpression of SNX15 leads to decreased cleavage of both insulin and HGF receptors [56], thereby highlighting a link between HGF and the sorting nexin family. Interestingly, our interval includes a member of the SNX family, SNX10. Further studies will be needed to determine putative association between such genes and the myopia phenotype within these families.

## ACKNOWLEDGMENTS

This work was supported by grants from the Association Retina France, the grant ‘Berthe Fouassier’ administrated by
The 'Fondation de France' (2005004235), and the Rotary Club. The authors want to thank the patients and their family, Dr. Heather Etchevers for her advice and constructive discussion, Dr. David Levy for his participation in patient recruitment, Dr. Frederic Lecerf and Dr. Matthias Titeux for their help in bioinformatics, and Mrs. Jennifer Jampoc for her translation aid.

## REFERENCES

1. Kempen JH, Mitchell P, Lee KE, Tielsch JM, Broman AT, Taylor HR, Ikram MK, Congdon NG, O'Colmain BJ, Eye Diseases Prevalence Research Group. The prevalence of refractive errors among adults in the United States, Western Europe, and Australia. Arch Ophthalmol 2004; 122:495-505. [PMID: 15078666]

2. Wong TY, Foster PJ, Hee J, Ng TP, Tielsch JM, Chew SJ, Johnson GJ, Seah SK. Prevalence and risk factors for refractive errors in adult Chinese in Singapore. Invest Ophthalmol Vis Sci 2000; 41:2486-94. [PMID: 10937558]

3. Lin LL, Shih YF, Hsiao CK, Chen CJ. Prevalence of myopia in Taiwanese schoolchildren: 1983 to 2000. Ann Acad Med Singapore 2004; 33:27-33. [PMID: 15008558]

4. Iwase A, Araie M, Tomidokoro A, Yamamoto T, Shimizu H, Kitazawa Y. Prevalence and causes of low vision and blindness in a Japanese adult population: the Tajimi Study. Ophthalmology 2006; 113:1354-62. [PMID: 16877074]

5. Javitt JC, Chiang YP. The socioeconomic aspects of laser refractive surgery. Arch Ophthalmol 1994; 112:1526-30. [PMID: 7993206]

6. Stambolian D, Ciner EB, Reider LC, Moy C, Dana D, Owens R, Schlifka M, Holmes T, Ibay G, Bailey-Wilson JE. Genome-wide scan for myopia in the Old Order Amish. Am J Ophthalmol 2005; 140:469-76. [PMID: 16084785]

7. Tay MT, Au Eong KG, Ng CY, Lim MK. Myopia and educational attainment in 421,116 young Singaporean males. Ann Acad Med Singapore 1992; 21:785-91. [PMID: 1295418]

8. Rose KA, Morgan IG, Smith W, Mitchell P. High heritability of myopia does not preclude rapid changes in prevalence. Clin Experiment Ophthalmol 2002; 30:168-72. [PMID: 12010207]

9. Congdon NG, Friedman DS, Lietman T. Important causes of visual impairment in the world today. JAMA 2003; 290:2057-60. [PMID: 14559961]

10. Munier A, Gunning T, Kenny D, O'Keefe M. Causes of blindness in the adult population of the Republic of Ireland. Br J Ophthalmol 1998; 82:630-3. [PMID: 9797662]

11. Dandona L, Dandona R, Srinivas M, Giridhar P, Vilas K, Prasad MN, John RK, McCarty CA, Rao GN. Blindness in the Indian state of Andhra Pradesh. Invest Ophthalmol Vis Sci 2001; 42:908-16. [PMID: 11274066]

12. Sorsby A, Fraser GR. Statistical Note on the Components of Ocular Refraction in Twins. J Med Genet 1964; 5:57-9. [PMID: 14205985]

13. Chen CJ, Cohen BH, Diamond EL. Genetic and environmental effects on the development of myopia in Chinese twin children. Ophthalmic Paediatr Genet 1985; 6:353-9. [PMID: 4069597]

14. Hammond CJ, Snieder H, Gilbert CE, Spector TD. Genes and environment in refractive error: the twin eye study. Invest Ophthalmol Vis Sci 2001; 42:1232-6. [PMID: 11328732]

15. Farbrother JE, Kirov G, Owen MJ, Guggenheim JA. Family aggregation of high myopia: estimation of the sibling recurrence risk ratio. Invest Ophthalmol Vis Sci 2004; 45:2873-8. [PMID: 15326097]

16. Dirani M, Chamberlain M, Garoufalis P, Chen C, Guymer RH, Baird PN. Refractive errors in twin studies. Twin Res Hum Genet 2006; 9:566-72. [PMID: 16899164]

17. Goss DA, Hampton MJ, Wickham MG. Selected review on genetic factors in myopia. J Am Optom Assoc 1988; 59:875-84. [PMID: 3068279]

18. Gilmartin B. Myopia: precedents for research in the twenty-first century. Clin Experiment Ophthalmol 2004; 32:305-24. [PMID: 15180846]

19. Morgan I, Rose K. How genetic is school myopia? Prog Retin Eye Res 2005; 24:1-38. [PMID: 15555525]

20. Schwartz M, Haim M, Skarsholm D. X-linked myopia: Bornholm eye disease. Linkage to DNA markers on the distal part of Xq. Clin Genet 1990; 38:281-6. [PMID: 1980096]

21. Zhang Q, Guo X, Xiao X, Jia X, Li S, Hejtmancik JF. Novel locus for X linked recessive high myopia maps to Xq23-q25 but outside MYP1. J Med Genet 2006; 43:e20. [PMID: 16648373]

22. Zhang Q, Li S, Xiao X, Jia X, Guo X. Confirmation of a genetic locus for X-linked recessive high myopia outside MYP1. J Hum Genet 2007; 52:469-72. [PMID: 17351708]

23. Young TL, Ronan SM, Alvear AB, Wildenberg SC, Oetting WS, Atwood LD, Wilkin DJ, King RA. A second locus for familial high myopia maps to chromosome 12q. Am J Hum Genet 1998; 63:1419-24. [PMID: 19792869]

24. Paluru P, Ronan SM, Heon E, Devoto M, Wildenberg SC, Scavello G, Holleschau A, Mäkitie O, Cole WG, King RA, Young TL. New locus for autosomal dominant high myopia maps to the long arm of chromosome 17. Invest Ophthalmol Vis Sci 2003; 44:1830-6. [PMID: 12714612]

25. Zhang Q, Guo X, Xiao X, Jia X, Li S, Hejtmancik J. A new locus for autosomal dominant high myopia maps to 4q22-q27 between D4S1578 and D4S1612. Mol Vis 2005; 11:554-60. [PMID: 16052171]

26. Nallasamy S, Paluru PC, Devoto M, Wasserman NF, Zhou J, Young TL. Genetic linkage study of high-grade myopia in a Hutterite population from South Dakota. Mol Vis 2007; 13:229-36. [PMID: 17327828]

27. Young TL, Ronan SM, Drahozal LA, Wildenberg SC, Alvear AB, Oetting WS, Atwood LD, Wilkin DJ, King RA. Evidence that a locus for familial high myopia maps to chromosome 18p. Am J Hum Genet 1998; 63:109-19. [PMID: 9634508]

28. Naiglin L, Gazagne C, Dallongeville F, Thalamas C, Idder A, Rascol O, Malecaze F, Calvas P. A genome wide scan for familial high myopia suggests a novel locus on chromosome 7q36. J Med Genet 2002; 39:118-24. [PMID: 11836361]

29. Paluru PC, Nallasamy S, Devoto M, Rappaport EF, Young TL. Identification of a novel locus on 2q for autosomal dominant high-grade myopia. Invest Ophthalmol Vis Sci 2005; 46:2300-7. [PMID: 15980214]

30. Chen CY, Stankovich J, Scurrah KJ, Garoufalis P, Dirani M, Pertile KK, Richardson AJ, Baird PN. Linkage replication of
the MYP12 locus in common myopia. Invest Ophthalmol Vis Sci 2007; 48:4433-9. [PMID: 17898262]

31. Stambolian D, Ibay G, Reider L, Dana D, Moy C, Schlifka M, Holmes T, Ciner E, Bailey-Wilson JE. Genomewide linkage scan for myopia susceptibility loci among Ashkenazi Jewish families shows evidence of linkage on chromosome 22q12. Am J Hum Genet 2004; 75:448-59. [PMID: 15273935]

32. Stambolian D, Ibay G, Reider L, Dana D, Moy C, Schlifka M, Holmes TN, Ciner E, Bailey-Wilson JE. Genome-wide scan of additional Jewish families confirms linkage of a myopia susceptibility locus to chromosome 22q12. Mol Vis 2006; 12:1499-505. [PMID: 17167407]

33. Wojciechowski R, Moy C, Ciner E, Ibay G, Reider L, Bailey-Wilson JE, Stambolian D. Genomewide scan in Ashkenazi Jewish families demonstrates evidence of linkage of ocular refraction to a QTL on chromosome 1p36. Hum Genet 2006; 119:389-99. [PMID: 16501916]

34. Hammond CJ, Andrew T, Mak YT, Spector TD. A susceptibility locus for myopia in the normal population is linked to the PAX6 gene region on chromosome 11: a genomewide scan of dizygotic twins. Am J Hum Genet 2004; 75:294-304. [PMID: 15307048]

35. Scavello GS, Paluru PC, Ganter WR, Young TL. Sequence variants in the transforming growth beta-induced factor (TGIF) gene are not associated with high myopia. Invest Ophthalmol Vis Sci 2004; 45:2091-7. [PMID: 15223781]

36. Young TL, Deeb SS, Ronan SM, Dewan AT, Alvear AB, Scavello GS, Paluru PC, Brott MS, Hayashi T, Holleschau AM, Benegas N, Schwartz M, Atwood LD, Oetting WS, Rosenberg T, Motulsky AG, King RA. X-linked high myopia associated with cone dysfunction. Arch Ophthalmol 2004; 122:897-908. [PMID: 15197065]

37. Liang CL, Wang HS, Hung KS, Hsi E, Sun A, Kuo YH, Juo SH. Evaluation of MMP3 and TIMP1 as candidate genes for high myopia in young Taiwanese men. Am J Ophthalmol 2006; 142:518-20. [PMID: 16935611]

38. Inamori Y, Ota M, Inoko H, Okada E, Nishizaki R, Shiota T, Mok J, Oka A, Ohno S, Mizuki N. The COL1A1 gene and high myopia susceptibility in Japanese. Hum Genet 2007; 122:151-7. [PMID: 17557158]

39. Liang CL, Hung KS, Tsai YY, Chang W, Wang HS, Juo SH. Systematic assessment of the tagging polymorphisms of the COL1A1 gene for high myopia. J Hum Genet 2007; 52:374-7. [PMID: 17273809]

40. Zhou G, Williams RW. Eye1 and Eye2: gene loci that modulate eye size, lens weight, and retinal area in the mouse. Invest Ophthalmol Vis Sci 1999; 40:817-25. [PMID: 10102277]

41. Han W, Yap MK, Wang J, Yip SP. Family-based association analysis of hepatocyte growth factor (HGF) gene polymorphisms in high myopia. Invest Ophthalmol Vis Sci 2006; 47:2291-9. [PMID: 16723436]

42. Naiglin L, Clayton J, Gazagne C, Dallongeville F, Malecaze F, Calvas P. Familial high myopia: evidence of an autosomal dominant mode of inheritance and genetic heterogeneity. Ann Genet 1999; 42:140-6. [PMID: 10526656]

43. Sambrook JFE, Maniatis T. Molecular cloning, a laboratory manual. New York: Cold Spring Harbor Laboratory Press; 1989.

44. Abecasis GR, Cherny SS, Cookson WO, Cardon LR. Merlin–rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet 2002; 30:97-101. [PMID: 11731797]

45. Wigginton JE, Abecasis GR. PEDSTATS: descriptive statistics, graphics and quality assessment for gene mapping data. Bioinformatics 2005; 21:3445-7. [PMID: 15947021]

46. Kruglyak L, Daly MJ, Reeve-Daly MP, Lander ES. Parametric and nonparametric linkage analysis: a unified multipoint approach. Am J Hum Genet 1996; 58:1347-63. [PMID: 8651312]

47. Yalcin B, Willis-Owen SA, Fullerton J, Meesaq A, Deacon RM, Rawlins JN, Copley RR, Morris AP, Flint J, Mott R. Genetic dissection of a behavioral quantitative trait locus shows that Rgs2 modulates anxiety in mice. Nat Genet 2004; 36:1197-202. [PMID: 15489855]

48. Altmuller J, Palmer LJ, Fischer G, Scherb H, Wjst M. Genomewide scans of complex human diseases: true linkage is hard to find. Am J Hum Genet 2001; 69:936-50. [PMID: 11565063]

49. Risch N, Merikangas K. The future of genetic studies of complex human diseases. Science 1996; 273:1516-7. [PMID: 8801636]

50. Klein AP, Duggal P, Lee KE, Klein R, Bailey-Wilson JE, Klein BE. Confirmation of linkage to ocular refraction on chromosome 22q and identification of a novel linkage region on 1q. Arch Ophthalmol 2007; 125:80-5. [PMID: 17210856]

51. Stone RA, Liu J, Sugimoto R, Capehart C, Zhu X, Pendrak K. GABA, experimental myopia, and ocular growth in chick. Invest Ophthalmol Vis Sci 2003; 44:3933-46. [PMID: 12939312]

52. Schulz HL, Stoehr H, White K, van Driel MA, Hoyng CB, Cremers F, Weber BH. Genomic structure and assessment of the retinally expressed RFamide-related peptide gene in dominant cystoid macular dystrophy. Mol Vis 2002; 8:67-71. [PMID: 11951088]

53. Jhamandas JH, Simonin F, Bourguignon JJ, Harris KH. Neuropeptide FF and neuropeptide VF inhibit GABAergic neurotransmission in parvocellular neurons of the rat hypothalamic paraventricular nucleus. Am J Physiol Regul Integr Comp Physiol 2007; 292:R1872-80. [PMID: 17289819]

54. Stone RA, Liu J, Sugimoto R, Capehart C, Zhu X, Pendrak K. GABA, experimental myopia, and ocular growth in chick. Invest Ophthalmol Vis Sci 2003; 44:3933-46. [PMID: 12939312]

55. Gong S, Zheng C, Doughty ML, Losos K, Didkovsky N, Schambra UB, Nowak NJ, Joyner A, Leblanc G, Hatten ME, Heintz N. A gene expression atlas of the central nervous system based on bacterial artificial chromosomes. Nature 2003; 425:917-25. [PMID: 14586460]

56. Phillips SA, Barr VA, Haft DH, Taylor SI, Haft CR. Identification and characterization of SNX15, a novel sorting nexin involved in protein trafficking. J Biol Chem 2001; 276:5074-84. [PMID: 11085978]

The print version of this article was created on 26 December 2008. This reflects all typographical corrections and errata to the article through that date. Details of any changes may be found in the online version of the article.