Purpose
To investigate the clinical and electrophysiologic natural history of Stargardt disease and correlate with the genotype.
Design
Cohort study of 59 patients.
Methods
Clinical history, examination, and electrophysiologic assessment were undertaken in a longitudinal survey. Patients were classified into 3 groups based on electrophysiologic findings, as previously published: Group 1 had dysfunction confined to the macula; Group 2 had macular and generalized cone system dysfunction; and Group 3 had macular and both generalized cone and rod system dysfunction. At baseline, there were 27 patients in Group 1, 17 in Group 2, and 15 in Group 3. Amplitude reduction of >50% in the relevant electroretinogram (ERG) component or a peak time shift of >3 ms for the 30 Hz flicker ERG or bright flash a-wave was considered clinically significant ERG deterioration. Molecular screening of ABCA4 was undertaken.
Results
The mean age at baseline was 31.7 years, with the mean follow-up interval being 10.5 years. A total of 22% of patients from Group 1 showed ERG group transition during follow-up, with 11% progressing to Group 2 and 11% to Group 3. Forty-seven percent of patients in Group 2 progressed to Group 3. There was clinically significant ERG deterioration in 54% of all subjects: 22% of Group 1, 65% of Group 2, and 100% of Group 3. At least 1 disease-causing ABCA4 variant was identified in 47 patients.
Conclusions
All patients with initial rod ERG involvement demonstrated clinically significant electrophysiologic deterioration; only 20% of patients with normal full-field ERGs at baseline showed clinically significant progression. Such data assist counseling by providing more accurate prognostic information and are also highly relevant in the design, patient selection, and monitoring of potential therapeutic interventions.
Stargardt disease is one of the most common inherited retinal disorders, with a prevalence of 1 in 10 000. It is inherited as an autosomal recessive trait. Most cases present with central visual loss and there is typically macular atrophy with yellow-white flecks at the posterior pole, which are at the level of the retinal pigment epithelium (RPE). Autofluorescence (AF) imaging and fluorescein angiography can be helpful in confirming the diagnosis. The age of onset is usually in the early teens, but there is wide variation, with a later age of onset being associated with a better visual prognosis.
Since the discovery of ABCA4 variants underlying Stargardt disease, multiple studies have described the wide phenotypic variability in ABCA4 -associated retinopathy. There is also extensive allelic heterogeneity, with more than 600 sequence variations having been reported to date in the ABCA4 gene. These 2 features make comprehensive genotype/phenotype correlations challenging. A previous cross-sectional study of 63 patients with Stargardt disease classified subjects into 3 functional electroretinogram (ERG) phenotypes: Group 1: dysfunction confined to the macula; Group 2: macular and generalized cone ERG abnormalities; and Group 3: macular and both generalized cone and rod ERG abnormalities. Differences in rod or cone function between groups could not be explained by differences in age of onset or duration of disease. It was thereby concluded that these 3 groups may represent distinct phenotypic subtypes of Stargardt disease and it was suggested, based on the cross-sectional data, that patients in Group 1 were likely to have a more favorable prognosis.
The purpose of the present study was to determine whether longitudinal data from a cohort of Stargardt disease patients support the value of full-field ERG to visual prognosis previously suggested by cross-sectional data. We have assessed the progression of Stargardt disease by repeated clinical and electrophysiologic examinations over time and probed whether the initial phenotype predicts long-term prognosis.
Patients and Methods
A cohort of 59 patients with a clinical diagnosis of Stargardt disease and a minimum of 7 years of follow-up were ascertained at Moorfields Eye Hospital. All patients were first diagnosed between 1997 and 2000, with the latest examinations performed between 2009 and 2011. The baseline clinical and electrophysiologic characteristics of 33 of these 59 patients have been previously reported. The panel included 5 sibships (4 sibling pairs and 1 set of 3 siblings). Informed consent was obtained from all participants. Blood samples were taken from all individuals for DNA extraction. The protocol of the study adhered to the tenets of the Declaration of Helsinki and was approved by the Ethics Committee of Moorfields Eye Hospital.
Clinical Assessment
Fifty-nine patients were assessed on at least 2 occasions, with the first and most recent visits taken as the baseline and “follow-up” examinations, respectively, for the purposes of data analysis. A full medical history was obtained and a comprehensive ophthalmologic examination performed for all patients. The age of onset was defined as the age at which visual loss was first noted by the patient. The duration of the disease was calculated as the difference between age at onset and age at the baseline examination when an electrophysiologic assessment was obtained. The interval of observation was determined by the difference between the age at baseline and the age at the most recent electrophysiologic examination. Clinical assessment included best-corrected Snellen visual acuity (converted to equivalent logarithm of minimal angle of resolution [logMAR] visual acuity for the purpose of data analysis), dilated ophthalmoscopy, and color fundus photography.
Electrophysiology
All patients underwent electrophysiologic assessment, to include full-field ERG and pattern electroretinography (PERG), incorporating the minimum standards of the International Society for Clinical Electrophysiology of Vision (ISCEV). ERG examination was comprehensive and included: (1) dark-adapted dim flash 0.01 candela second (cd·s)/m 2 (dark-adapted 0.01); (2) dark-adapted bright flash 11.0 cd·s/m 2 (dark-adapted 11.0); (3) light-adapted 3.0 cd·s/m 2 30 Hz flicker ERG (light-adapted 30 Hz); and (4) light-adapted 3.0 cd·s/m 2 at 2 Hz (light-adapted 3.0). All recordings were performed with gold-foil recording electrodes with reference electrodes at the ipsilateral outer canthi.
The patient data were compared against those of 16 healthy subjects younger than 50 years and 19 subjects older than 50 years, to maintain consistency with the original cross-sectional study. PERGs were compared against those from 28 normal subjects, with N95 peak time not being used for interpretation because of its accepted variability. The limits of ERG normality were defined for all the components of the ERG and PERG as the mean value ±2 standard deviations ( Supplemental Tables 1 and 2 , available at AJO.com ). The threshold values for the minimum amplitude/maximum peak time for subjects younger than 50 years were defined as 135 μV/107 ms (dark-adapted 0.01), 250 μV/13 ms and 320 μV/56 ms (dark-adapted 11.0 a- and b-wave, respectively), 70 μV/27 ms (light-adapted 30 Hz), and 30 μV/15 ms and 95 μV/32 ms (light-adapted 3.0 a- and b-wave, respectively); and for patients older than 50 years as 30 μV/117 ms (dark-adapted 0.01), 105 μV/16 ms and 235 μV/57 ms (dark-adapted 11.0 a- and b-wave, respectively), 50 μV/29 ms (light-adapted 30 Hz); and 15 μV/16 ms and 90 μV/32 ms (light-adapted 3.0 a- and b-wave, respectively). The threshold values for the PERG P50 minimum amplitude/maximum peak time were defined as 2.1 μV/58.5 ms.
All the components of the ERG and PERG from each eye were taken into account when classifying patients into 1 of the 3 ERG groups at baseline and follow-up. Group 1 was defined as PERG abnormality with normal ERGs. In Group 2, there was PERG abnormality and abnormal cone function (assessed with light-adapted 30 Hz and light-adapted 3.0) on ERG. In Group 3, there was additional rod ERG abnormality (assessed using dark-adapted 0.01 and dark-adapted 11.0). The overall classification was based on the more severe eye in the small number of patients with different ERG groups between eyes. The data obtained at follow-up were compared with those at baseline. Concordance for ERG group between siblings was defined as siblings having the same ERG group classification both at baseline and at follow-up.
Amplitude reduction was calculated as the difference between amplitude at baseline and at follow-up. The percentage reduction in amplitudes was obtained by dividing the amplitude reduction by baseline amplitude. A yearly amplitude reduction and a yearly percentage reduction were calculated by dividing the amplitude reduction or the percentage reduction by the follow-up time. A yearly peak time shift (difference between peak time at baseline and at follow-up) was also calculated by dividing by the follow-up time.
An amplitude reduction of over 50% in any ERG component and/or a peak time shift of over 3 ms for the light-adapted 30 Hz ERG or dark-adapted 11.0 ERG a-wave were considered evidence of clinically significant ERG deterioration/progression. Patients were thereby classified into 2 subsets: those with clinically significant ERG deterioration and those without significant ERG deterioration (stable ERG).
Mutation Screening
Mutation analysis was performed using the single-stranded conformation polymorphism (SSCP) strategy of the whole coding region of ABCA4 in 33 subjects and the arrayed primer extension (APEX) microarray (ABCR400 chip; Asper Ophthalmics, Tartu, Estonia) for previously reported variants in 27 patients. Direct Sanger sequencing was done in siblings of probands and parents, when available, to confirm segregation of alleles, as well as in 8 subjects either to confirm putative novel variants or where the variants found with SSCP and APEX differed ( Supplemental Table 3 , available at AJO.com ).
Non-null variants were analyzed using 2 software prediction programs: SIFT (Sorting Intolerant from Tolerant; http://sift.jcvi.org/ ) and PolyPhen2 ( http://genetics.bwh.harvard.edu/pph/index.html ). All variants were also analyzed for their effect on splicing using the Human Splicing Finder program, version 2.4.1 ( http://www.umd.be/HSF/ ). All variants were compared with variants in the Exome Variant Server, NHLBI Exome Sequencing Project, Seattle, Washington, USA ( http://snp.gs.washington.edu/EVS/ ).
Each patient was classified into 4 mutually exclusive genotype groups on the basis of the molecular analysis: (A) patients with at least 1 null variant, (B) subjects with 2 or more non-null variants, (C) individuals with 1 non-null variant, and (D) patients with no detectable variants. Null variants were those that would be expected to affect splicing, or to introduce a premature truncating codon in the protein if translated. The term “variants” for the purpose of this study includes those sequence changes previously shown to be enriched in Stargardt patients from prior studies, or for very rare variants, those not found at an allele frequency greater than 0.1% on the exome variant database (Accessed March 1, 2012).
Statistical Analysis
Statistical analysis has been undertaken using data from only 1 eye in each subject. For the 57 patients with the same ERG grouping in both eyes, the eye used for analysis was selected according to the Random Integer Generator ( http://www.random.org/ ). For the 2 patients (Patients 26 and 48) with a different ERG group in each eye, the eye with the more severe ERG grouping (ie, more generalized retinal dysfunction) was selected for analysis.
The Mann-Whitney U test was used to explore whether differences observed between patients with clinically significant electrophysiologic deterioration and those without were statistically significant with regard to age of onset, duration of disease, age at baseline, the interval of observation, logMAR visual acuity at baseline, logMAR visual acuity reduction (defined as the difference between visual acuity at baseline and at follow-up), and yearly percentage amplitude reduction and yearly peak time shift in both the light-adapted 11.0 a-wave and light-adapted 30 Hz.
The Kruskal-Wallis test with Steel-Dwass multiple comparisons was performed to compare the 3 baseline ERG groups (ERG Group 1, 2, and 3) and the 3 genotype groups (genotype A, B, and C) for the 10 aforementioned parameters. Where evidence was found of a difference between these groups, all pairwise comparisons were made.
The association between genotype group classification and baseline ERG group classification was tested using the Goodman-Kruskal gamma, a measure of association for ordered categories ranging between −1 and +1 for perfect negative or positive association, respectively. P values less than .05 were considered to indicate statistical significance.
All analyses were conducted using MedCalc statistical software version 9.2.1.0 (MedCalc Software, Ostend, Belgium) and Excel Tokei 2010 (Social Survey Research Information Co Ltd, Tokyo, Japan).
Results
Clinical Findings
Fifty-nine patients, 31 female (52%, 31/59) and 28 male (48%, 28/59), were included in the study. All complained of central visual loss with a median age of onset of 20.8 years (range, 5-48 years) and a median duration of disease of 10.9 years (range, 0-31 years). The median ages at baseline and at follow-up were 31.7 and 42.2 years (range, 8-64 and 20-73 years), respectively. The mean follow-up interval was 10.5 years (range, 7-13 years). Seven patients (12%, 7/59) presented before 16 years of age and 52 (88%, 52/59) presented after age 16 years. The median logMAR visual acuities (VA) at baseline and at follow-up were 0.93 (range, 0.0-2.0) and 1.22 (range, 0.0-3.0), respectively, with a median logMAR VA reduction during the follow-up interval of 0.29 (range, −0.78-2.0). The clinical findings are summarized in Table 1 and the eye selected for data analysis is shown in Supplemental Table 3 (available at AJO.com ).
Pt | Onset (y) | Age (y) | logMAR VA | Variants Identified a | ||
---|---|---|---|---|---|---|
BL | FU | BL | FU | |||
1 | 16 | 17 | 26 | 0.0/1.0 | 0.0/0.48 | c.768G>T / p.Gly863Ala / p.Arg943Gln |
2 | 15 | 17 | 25 | 0.78/0.78 | 1.0/1.0 | p. Arg1443His |
3 | 11 | 18 | 27 | 0.78/1.0 | 1.0/1.0 | p.Trp439* / p.Gly863Ala / p.Leu1970Phe |
4 | 19 | 21 | 32 | 0.78/0.78 | 1.0/1.0 | p.Leu2027Phe |
5 | 10 | 22 | 30 | 0.48/0.48 | 1.0/0.78 | p.Gly863Ala / p.Arg943Gln / c.5461-10 T>C |
6 | 18 | 26 | 37 | 0.78/1.0 | 1.0/1.0 | p.Pro1380Phe |
7 | 25 | 28 | 40 | 0.78/1.0 | 1.3/0.78 | ND |
8 | 24 | 29 | 38 | 1.0/0.78 | 1.0/1.0 | p.Phe418Ser / p.Leu2027Phe |
9 | 24 | 31 | 44 | 1.0/1.0 | 1.3/1.0 | c.4253+5 G>T / p.Gly1507Arg |
10 | 26 | 32 | 44 | 0.78/0.78 | 1.0/1.0 | p.Cys1490Tyr / p.Arg2030Gln |
11 | 31 | 34 | 46 | 0.18/0.3 | 0.6/0.7 | ND |
12 | 17 | 35 | 47 | 1.0/1.0 | 1.0/1.0 | p.Asn96His |
13 | 23 | 35 | 45 | 1.0/0.3 | 1.0/0.48 | p.Gly1513Profs*1554 |
14 | 33 | 37 | 48 | 0.18/1.48 | 1.0/1.3 | ND |
15 | 38 | 40 | 51 | 0.18/0.78 | 1.0/1.0 | p.Arg2107His |
16 | 42 | 43 | 53 | 0.0/0.0 | 1.0/1.0 | ND |
17 | 22 | 48 | 59 | 1.0/1.0 | 1.0/1.0 | p.Cys54Tyr |
18 | 20 | 49 | 59 | 1.0/0.6 | 1.0/1.0 | p.Pro1380Leu / p.Gly1961Glu |
19 | 35 | 50 | 61 | 1.0/0.3 | 1.0/1.0 | p.Arg1108Cys |
20 | 25 | 56 | 67 | 1.3/0.18 | 1.0/1.0 | p.Trp439* / p.Gly863Ala |
21 | 48 | 59 | 71 | 1.0/0.78 | 1.0/1.0 | p. Ile156 Val / p. Cys1455Arg / p. Phe1839Ser |
22 | 21 | 22 | 31 | 0.3/1.0 | 1.0/1.0 | p.Arg2107His |
23 | 21 | 23 | 33 | 1.0/1.0 | 1.0/1.0 | p.Gly863Ala |
24 | 48 | 64 | 73 | 0.0/1.0 | 0.18/3.0 | p.Tyr1652* |
25 | 17 | 19 | 29 | 0.78/0.3 | 1.0/1.0 | c.5461-10 T>C |
26 | 17 | 21 | 33 | 1.0/0.78 | 1.0/1.0 | ND |
27 | 27 | 53 | 66 | 1.78/1.78 | 1.3/1.0 | p.Ser1071Cysfs*1084 |
28 | 5 | 14 | 21 | 0.78/0.78 | 1.0/1.0 | p.Arg408* / p.Val675lle |
29 | 9 | 15 | 27 | 1.08/1.08 | 1.0/1.0 | p.Cys2150Tyr |
30 | 14 | 24 | 32 | 1.0/0.78 | 1.0/1.0 | ND |
31 | 18 | 28 | 39 | 1.0/1.0 | 1.0/1.0 | p.Gly863Ala / p.Arg1108Cys / p.Arg943Gln |
32 | 14 | 29 | 37 | 1.0/1.0 | 1.0/1.0 | p.Arg653Cys / p.Arg2030Gln |
33 | 19 | 29 | 40 | 1.0/1.0 | 1.0/1.08 | ND |
34 | 34 | 40 | 49 | 0.3/0.48 | 1.0/1.0 | p.Gly863Ala / p.Glu1087Lys |
35 | 25 | 43 | 54 | 1.0/1.0 | 1.0/1.0 | p.Cys54Tyr / p.Gly863Ala |
36 | 38 | 60 | 69 | 1.0/1.0 | 1.3/1.08 | p.Val931Met / c.5461-10 T>C |
37 | 10 | 11 | 20 | 1.0/0.78 | 1.3/1.3 | p.Pro1380Leu |
38 | 10 | 15 | 23 | 1.0/1.0 | 1.3/1.3 | p.Ser1071Cysfs*1084 / p.Pro1380Leu |
39 | 24 | 25 | 38 | 1.56/0.3 | 2.0/2.0 | c.5461-10 T>C / c.5714+5 G>A |
40 | 18 | 26 | 36 | 1.3/1.3 | 2.0/1.3 | ND |
41 | 32 | 33 | 45 | 0.48/0.48 | 1.0/1.0 | ND |
42 | 32 | 35 | 46 | 1.3/0.0 | 3.0/1.0 | p.Cys54Tyr |
43 | 30 | 35 | 45 | 0.48/0.48 | 2.0/1.3 | ND |
44 | 15 | 41 | 49 | 1.3/1.3 | 2.0/1.3 | p.Asn965Ser |
45 | 8 | 8 | 20 | 0.78/0.78 | 1.0/1.0 | p.Thr1019Met |
46 | 10 | 11 | 23 | 1.0/1.0 | 1.0/1.0 | p.Thr1019Met |
47 | 8 | 12 | 24 | 2.0/1.56 | 1.78/1.48 | p.Cys2150Tyr |
48 | 17 | 18 | 26 | 1.0/0.78 | 1.3/1.0 | c.5461-10 T>C / p.Leu2027Phe |
49 | 8 | 21 | 33 | 1.3/1.3 | 2.0/2.0 | p.Asp574Aspfs*582 |
50 | 8 | 27 | 39 | 2.0/1.56 | 1.78/1.48 | c.5461-10 T>C |
51 | 24 | 31 | 43 | 1.18/1.18 | 1.08/1.3 | p.Arg1640Trp / p.Leu2027Phe |
52 | 11 | 31 | 42 | 1.3/1.3 | 2.0/2.0 | p.Arg1108His |
53 | 5 | 32 | 43 | 2.0/2.0 | 2.0/2.0 | c.5461-10 T>C / p.Cys2150Tyr |
54 | 5 | 32 | 43 | 2.0/2.0 | 2.0/2.0 | c.5461-10 T>C / p.Cys2150Tyr |
55 | 7 | 36 | 47 | 1.3/1.3 | 3.0/1.3 | c.5461-10 T>C / p.Cys2150Tyr |
56 | 13 | 39 | 50 | 1.25/1.56 | 3.0/3.0 | ND |
57 | 23 | 42 | 52 | 1.56/1.0 | 1.0/1.0 | p.Leu747Cysfs*787 |
58 | 40 | 43 | 54 | 0.18/0.18 | 0.78/0.78 | ND |
59 | 23 | 54 | 65 | 0.78/1.0 | 1.0/1.0 | p.Ile156Val |
At baseline, there were 27 patients (46%, 27/59) in Group 1, 17 (29%, 17/59) in Group 2, and 15 (25%, 15/59) in Group 3, compared at follow-up to 21 patients (36%, 21/59) in Group 1, 12 (20%, 12/59) in Group 2, and 26 (44%, 26/59) in Group 3 ( Table 2 ). The median age of onset for each baseline ERG group was 24.9 years in Group 1, 20.4 years in Group 2, and 14.0 years in Group 3. The median age (years) at examination/logMAR visual acuity at baseline and follow-up for each baseline ERG group was 34.4/0.78 and 45.0/1.00, respectively, in Group 1; 29.6/1.00 and 39.4/1.00, respectively, in Group 2; and 29.1/1.25 and 40.3/1.30, respectively, in Group 3 ( Table 3 ).
Electrophysiologic Group a at Baseline b | Electrophysiologic Group a at Follow-up b | ||
---|---|---|---|
Group 1 | Group 2 | Group 3 | |
Group 1 (n = 27, 6 ) | 21 | 3 ( 3 ) | 3 ( 3 ) |
Group 2 (n = 17, 11 ) | 9 ( 3 ) | 8 ( 8 ) | |
Group 3 (n = 15, 15 ) | 15 ( 15 ) | ||
Total (n = 59, 32 ) | 21 | 12 ( 6 ) | 26 ( 26 ) |
a Patients were classified into 3 groups based on electrophysiologic findings: Group 1 had dysfunction confined to the macula; Group 2 had macular and generalized cone system dysfunction; Group 3 had macular and both generalized cone and rod system dysfunction.
b Numbers in bold show the numbers of patients who demonstrated electrophysiologic evidence of deterioration. An amplitude reduction of over 50% in any electrophysiologic component and/or a peak time shift of over 3 ms for the light-adapted 30 Hz electroretinogram or dark-adapted 11.0 electroretinogram a-wave were considered evidence of significant electrophysiologic deterioration.
Median Age of Onset (y) | Median Age | Median logMAR Visual Acuity | ||||
---|---|---|---|---|---|---|
BL | FL | BL | FL | |||
Baseline electrophysiologic group | Group 1 (n = 27) | 24.9 | 34.4 | 45.0 | 0.78 | 1.00 |
Group 2 (n = 17) | 20.4 | 29.6 | 39.4 | 1.00 | 1.00 | |
Group 3 (n = 15) | 14.0 | 29.1 | 40.3 | 1.25 | 1.30 | |
Evidence of clinically significant electrophysiologic deterioration a | Stable (n = 27) | 23.4 | 33.5 | 43.8 | 0.78 | 1.00 |
Significant deterioration (n = 32) | 18.7 | 30.1 | 40.8 | 1.00 | 1.19 | |
Genotype grouping b | Genotype A (n = 19) | 17.6 | 32.6 | 42.1 | 1.08 | 1.39 |
Genotype B (n = 10) | 22.3 | 35.7 | 48.2 | 0.84 | 0.94 | |
Genotype C (n = 18) | 20.0 | 27.8 | 38.4 | 0.90 | 1.20 | |
Genotype D (n = 12) | 26.1 | 32.7 | 43.5 | 0.69 | 1.19 | |
Total (n = 59) | 20.8 | 31.7 | 42.2 | 0.93 | 1.22 |
a The subset without evidence of significant deterioration is described as “Stable.”
b Each patient was classified into 4 mutually exclusive genotype groups on the basis of the molecular analysis: (A) patients with at least 1 null variant, (B) subjects with 2 or more non-null variants, (C) individuals with 1 non-null variant, and (D) patients with no detectable variants.
Color fundus photographs of eyes in 3 representative cases (Patients 17, 42, and 53) are shown in Figure 1 ; their respective electrophysiologic traces appear in Figure 2 . Patient 17 showed no ERG group transition (Group 1 at baseline and Group 1 at follow-up). ERG transition from Group 2 to Group 3, with clinically significant ERG deterioration, was demonstrated in Patient 42. Patient 53 was in ERG Group 3 at baseline and had evidence of clinically significant ERG deterioration.
Electrophysiologic Findings
The electrophysiologic findings are summarized in Supplemental Table 4 (available at AJO.com ). PERG P50 components were undetectable (93%, 51/55) or moderately reduced (7%, 4/55; Patients 16, 24, 42, and 55) at baseline, in keeping with severe or moderately severe macular dysfunction; and were undetectable in 53 individuals (96%, 53/55) or moderately reduced in 2 patients (4%, 2/55; Patients 16 and 24) at follow-up. There were no available PERG data both at baseline and at follow-up in 2 subjects (Patients 7 and 23), and no available baseline PERGs in 2 further individuals (Patients 45 and 46), who had undetectable PERGs at follow-up.
Complete ERG data sets were available at baseline and follow-up, with few exceptions ( Supplemental Table 4 ). The dark-adapted 0.01 and dark-adapted 11.0 ERGs were abnormal in 11 and 15 patients (20%, 11/54 and 25%, 15/59), respectively, at baseline, and in 22 and 24 subjects (36%, 22/59 and 41%, 24/59), respectively, at follow-up. All those with abnormal dark-adapted 0.01 ERGs had abnormal light-adapted 30 Hz and light-adapted 3.0 ERGs. Three out of 4 patients (Patients 53-56) with undetectable dark-adapted 0.01 responses at follow-up had undetectable light-adapted ERGs at baseline and at follow-up.
Light-adapted 30 Hz and light-adapted 3.0 ERGs were abnormal in 29 and 26 patients (49%, 29/59, and 45%, 26/58), respectively, at baseline; and in 38 and 36 subjects (64%, 38/59 and 61%, 36/59), respectively, at follow-up. An abnormal light-adapted 3.0 ERG was the only baseline ERG abnormality in 2 patients (Patients 29 and 41); isolated light-adapted 30 Hz ERG abnormality occurred in another 4 subjects (Patients 28, 30, 42, and 48). All 6 of these patients showed abnormal responses in both light-adapted tests at follow-up. Isolated light-adapted 30 Hz ERG abnormality occurred in another 2 patients at follow-up.
Four out of 5 sibships were concordant (the same ERG group) both at baseline and at follow-up (Patients 11 and 14; 40 and 42; 45 and 46; 53-55). Two siblings from 1 family had discordant ERG groups, with 1 sibling in Group 3 at baseline and follow-up and the other sibling in Group 2 at baseline and follow-up (Patients 47 and 29) ( Supplemental Table 4 ).
The clinical features of each baseline group are summarized in Table 3 and Figure 3 . There was a statistically significant difference between Groups 1 and 3 and between Groups 2 and 3 in terms of onset of disease ( Supplemental Table 5 , available at AJO.com ). There was also a statistically significant difference in logMAR VA between Groups 1 and 3 and between Groups 2 and 3. No statistically significant difference was seen between groups with respect to age at baseline, duration of disease, and follow-up interval. Mean yearly electrophysiologic progression within each baseline ERG group with respect to dark-adapted 11.0 a-wave and light-adapted 30 Hz is summarized in Table 4 and Figure 3 . Statistical analysis revealed a significant difference between Groups 1 and 3 and between Groups 2 and 3 in terms of yearly amplitude reduction of dark-adapted 11.0 a-wave ( Supplemental Table 5 ). There was also a statistically significant difference in light-adapted 30 Hz yearly peak time shift between Groups 1 and 3. No significant difference was seen between groups with respect to amplitude reduction in light-adapted 30 Hz.
Dark-Adapted 11.0 A-wave | Light-Adapted 30 Hz | |||||
---|---|---|---|---|---|---|
Amplitude Reduction (μV/y) | Percentage Reduction (%/y) | Peak Time Shift (ms/y) | Amplitude Reduction (μV/y) | Percentage Reduction (%/y) | Peak Time Shift (ms/y) | |
Group 1 (n = 27) | 5.5 | 1.7 | 0.10 | 2.7 | 2.2 | 0.14 |
Group 2 (n = 17) | 4.5 | 1.5 | 0.09 | 1.1 | 1.7 | 0.19 |
Group 3 (n = 15) | 4.9 | 3.6 | 0.18 | 1.5 | 3.1 | 0.32 |
Stable (n = 27) | 3.9 | 1.2 | 0.04 | 2.2 | 1.9 | 0.07 |
Electrophysiologic Deterioration (n = 32) | 6.0 | 2.9 | 0.18 | 1.7 | 2.7 | 0.31 |
Genotype A (n = 19) | 6.5 | 3.0 | 0.14 | 2.3 | 3.0 | 0.23 |
Genotype B (n = 10) | 2.3 | 0.5 | −0.01 | 1.4 | 0.9 | 0.12 |
Genotype C (n = 18) | 5.4 | 2.1 | 0.16 | 2.4 | 3.1 | 0.33 |
Genotype D (n = 12) | 4.3 | 2.1 | 0.09 | 1.1 | 0.9 | −0.04 |
Total (n = 59) | 5.1 | 2.1 | 0.11 | 1.9 | 2.3 | 0.19 |
a A yearly amplitude reduction and a yearly percentage reduction were calculated by dividing the amplitude reduction or the percentage reduction by the follow-up time. A yearly peak time shift (difference between peak time at baseline and follow-up) was also calculated by dividing by the follow-up time.
Thirty-two patients showed evidence of clinically significant electrophysiologic deterioration ( Table 2 and Supplemental Table 4 ). Twenty-one subjects showed a greater than 50% amplitude reduction and 26 patients had more than a 3 ms peak time shift ( Supplemental Table 4 ). The clinical findings were compared between the subset of patients with evidence of ERG progression and those without (stable ERG) ( Table 3 and Figure 4 ). There was a statistically significant difference between the 2 subsets in terms of age of onset and logMAR VA at baseline ( Supplemental Table 5 and Figure 4 ). There were no statistically significant differences between the 2 subsets with respect to age at baseline, duration of disease, interval of follow-up, and reduction in logMAR VA ( Supplemental Table 5 and Figure 4 ).
There was clinically significant deterioration of ERG parameters in 22% (6/27) of patients in ERG Group 1, 65% (11/17) in Group 2, and 100% (15/15) in Group 3 ( Table 2 ). Patients with a Group 1 ERG phenotype both at baseline and at follow-up did not show significant electrophysiologic deterioration (78%, 21/27), with the Group 1 subjects (22%, 6/27) who did show ERG progression all moving to either Group 2 or Group 3 in equal proportions. Mean yearly electrophysiologic progression was compared between patients with and without clinically significant ERG deterioration ( Table 4 and Figure 4 ). Statistical analysis revealed a significant difference in terms of both amplitude reduction and peak time shift of dark-adapted 11.0 a-wave ( Supplemental Table 5 and Figure 4 ). There was also a statistically significant difference in light-adapted 30 Hz peak time shift. No significant difference was seen with respect to rate of amplitude reduction in light-adapted 30 Hz ( Supplemental Table 5 ).
Molecular Genetics
Likely disease-causing variants in ABCA4 were detected in 47 out of 59 patients, with 2 or more variants identified in 22 patients and 1 variant in 25 subjects ( Table 1 and Supplemental Table 6 , available at AJO.com ). Nineteen patients had at least 1 null variant, 10 subjects had 2 or more non-null variants, 18 individuals were identified with 1 non-null variant, and 12 patients had no detectable variants. Detailed results, including in silico analysis to assist in the prediction of pathogenicity of the variants, are shown in Supplemental Table 7 (available at AJO.com ).
Thirty-eight different variants were found in 47 patients: 11 null mutations with 3 predicted to affect splicing, and 27 non-null variants ( Supplemental Tables 6 and 7 ). Eighteen patients harbored at least 1 null variant, with a single subject having 2 null mutations. Thirty-two of these 38 variants have been previously reported and 6 are putative novel mutations: (1) c.1317G>A, p.Trp439*, (2) c.2103G>A, p.Val675lle, (3) c.2239delC, p.Leu747Cysfs*787, (4) c.4363C>T, p.Cys1455Arg, (5) c.4519G>A, p.Gly1507Arg, and (6) c.5516T>C, p.Phe1839Ser ( Supplemental Tables 6 and 7 ). At least 1 variant was identified in 22 patients (81%, 22/27) in ERG Group 1 at baseline, 12 (71%, 12/17) in Group 2, and 13 (87%, 13/15) in Group 3. At least 1 null variant was found in 8 patients (30%, 8/27) in ERG Group 1 at baseline, 4 (24%, 4/17) in Group 2, and 7 (47%, 7/15) in Group 3 ( Supplemental Table 6 and Supplemental Figure 1 , available at AJO.com ).
Genotype-Phenotype Correlations
Clinical features at baseline and electrophysiologic progression in dark-adapted 11.0 a-wave and light-adapted 30 Hz of each genotype group are summarized in Tables 3 and 4 . There was no statistically significant association identified between the severity of genotype and the extent of electrophysiologic dysfunction on the basis of baseline ERG grouping (γ = −0.126), although patients with 2 or more non-null variants (genotype B group) less frequently had rod ERG involvement ( Table 5 and Supplemental Figure 1 ).
Genotype A | Genotype B | Genotype C | Genotype D | |
---|---|---|---|---|
Group 1 (n = 27) | 8 | 5 | 9 | 5 |
Group 2 (n = 17) | 4 | 4 | 4 | 5 |
Group 3 (n = 15) | 7 | 1 | 5 | 2 |
Stable (n = 27) | 6 | 9 | 7 | 5 |
Electrophysiologic deterioration (n = 32) a | 13 | 1 | 11 | 7 |
Total (n = 59) | 19 | 10 | 18 | 12 |
a The subset without evidence of significant deterioration is described as “Stable.”
The distribution of patients with clinically significant electrophysiologic deterioration in each genotype group is shown in Table 5 and Supplemental Figure 2 (available at AJO.com ). Statistical analysis revealed a significant difference between genotype groups A and B and between genotype groups A and C in terms of age of onset. There was also a statistically significant difference between genotype groups A and B with respect to yearly amplitude reduction of dark-adapted 11.0 a-wave and light-adapted 30 Hz yearly peak time shift ( Supplemental Table 5 ). No statistically significant difference was seen between genotype groups and the other ERG parameters ( Supplemental Table 5 ).
Interestingly, 8 of the 9 patients harboring the variant c.5461-10 T>C (Patients 5, 25, 36, 39, 48, 50, 53-55) had clinically significant ERG progression. All 3 unrelated patients (1, 5, and 31) harboring p.Arg943Gln also had p.Gly863Ala, suggesting linkage disequilibrium of these 2 substitutions, with none of these subjects having clinically significant ERG deterioration.
Results
Clinical Findings
Fifty-nine patients, 31 female (52%, 31/59) and 28 male (48%, 28/59), were included in the study. All complained of central visual loss with a median age of onset of 20.8 years (range, 5-48 years) and a median duration of disease of 10.9 years (range, 0-31 years). The median ages at baseline and at follow-up were 31.7 and 42.2 years (range, 8-64 and 20-73 years), respectively. The mean follow-up interval was 10.5 years (range, 7-13 years). Seven patients (12%, 7/59) presented before 16 years of age and 52 (88%, 52/59) presented after age 16 years. The median logMAR visual acuities (VA) at baseline and at follow-up were 0.93 (range, 0.0-2.0) and 1.22 (range, 0.0-3.0), respectively, with a median logMAR VA reduction during the follow-up interval of 0.29 (range, −0.78-2.0). The clinical findings are summarized in Table 1 and the eye selected for data analysis is shown in Supplemental Table 3 (available at AJO.com ).
Pt | Onset (y) | Age (y) | logMAR VA | Variants Identified a | ||
---|---|---|---|---|---|---|
BL | FU | BL | FU | |||
1 | 16 | 17 | 26 | 0.0/1.0 | 0.0/0.48 | c.768G>T / p.Gly863Ala / p.Arg943Gln |
2 | 15 | 17 | 25 | 0.78/0.78 | 1.0/1.0 | p. Arg1443His |
3 | 11 | 18 | 27 | 0.78/1.0 | 1.0/1.0 | p.Trp439* / p.Gly863Ala / p.Leu1970Phe |
4 | 19 | 21 | 32 | 0.78/0.78 | 1.0/1.0 | p.Leu2027Phe |
5 | 10 | 22 | 30 | 0.48/0.48 | 1.0/0.78 | p.Gly863Ala / p.Arg943Gln / c.5461-10 T>C |
6 | 18 | 26 | 37 | 0.78/1.0 | 1.0/1.0 | p.Pro1380Phe |
7 | 25 | 28 | 40 | 0.78/1.0 | 1.3/0.78 | ND |
8 | 24 | 29 | 38 | 1.0/0.78 | 1.0/1.0 | p.Phe418Ser / p.Leu2027Phe |
9 | 24 | 31 | 44 | 1.0/1.0 | 1.3/1.0 | c.4253+5 G>T / p.Gly1507Arg |
10 | 26 | 32 | 44 | 0.78/0.78 | 1.0/1.0 | p.Cys1490Tyr / p.Arg2030Gln |
11 | 31 | 34 | 46 | 0.18/0.3 | 0.6/0.7 | ND |
12 | 17 | 35 | 47 | 1.0/1.0 | 1.0/1.0 | p.Asn96His |
13 | 23 | 35 | 45 | 1.0/0.3 | 1.0/0.48 | p.Gly1513Profs*1554 |
14 | 33 | 37 | 48 | 0.18/1.48 | 1.0/1.3 | ND |
15 | 38 | 40 | 51 | 0.18/0.78 | 1.0/1.0 | p.Arg2107His |
16 | 42 | 43 | 53 | 0.0/0.0 | 1.0/1.0 | ND |
17 | 22 | 48 | 59 | 1.0/1.0 | 1.0/1.0 | p.Cys54Tyr |
18 | 20 | 49 | 59 | 1.0/0.6 | 1.0/1.0 | p.Pro1380Leu / p.Gly1961Glu |
19 | 35 | 50 | 61 | 1.0/0.3 | 1.0/1.0 | p.Arg1108Cys |
20 | 25 | 56 | 67 | 1.3/0.18 | 1.0/1.0 | p.Trp439* / p.Gly863Ala |
21 | 48 | 59 | 71 | 1.0/0.78 | 1.0/1.0 | p. Ile156 Val / p. Cys1455Arg / p. Phe1839Ser |
22 | 21 | 22 | 31 | 0.3/1.0 | 1.0/1.0 | p.Arg2107His |
23 | 21 | 23 | 33 | 1.0/1.0 | 1.0/1.0 | p.Gly863Ala |
24 | 48 | 64 | 73 | 0.0/1.0 | 0.18/3.0 | p.Tyr1652* |
25 | 17 | 19 | 29 | 0.78/0.3 | 1.0/1.0 | c.5461-10 T>C |
26 | 17 | 21 | 33 | 1.0/0.78 | 1.0/1.0 | ND |
27 | 27 | 53 | 66 | 1.78/1.78 | 1.3/1.0 | p.Ser1071Cysfs*1084 |
28 | 5 | 14 | 21 | 0.78/0.78 | 1.0/1.0 | p.Arg408* / p.Val675lle |
29 | 9 | 15 | 27 | 1.08/1.08 | 1.0/1.0 | p.Cys2150Tyr |
30 | 14 | 24 | 32 | 1.0/0.78 | 1.0/1.0 | ND |
31 | 18 | 28 | 39 | 1.0/1.0 | 1.0/1.0 | p.Gly863Ala / p.Arg1108Cys / p.Arg943Gln |
32 | 14 | 29 | 37 | 1.0/1.0 | 1.0/1.0 | p.Arg653Cys / p.Arg2030Gln |
33 | 19 | 29 | 40 | 1.0/1.0 | 1.0/1.08 | ND |
34 | 34 | 40 | 49 | 0.3/0.48 | 1.0/1.0 | p.Gly863Ala / p.Glu1087Lys |
35 | 25 | 43 | 54 | 1.0/1.0 | 1.0/1.0 | p.Cys54Tyr / p.Gly863Ala |
36 | 38 | 60 | 69 | 1.0/1.0 | 1.3/1.08 | p.Val931Met / c.5461-10 T>C |
37 | 10 | 11 | 20 | 1.0/0.78 | 1.3/1.3 | p.Pro1380Leu |
38 | 10 | 15 | 23 | 1.0/1.0 | 1.3/1.3 | p.Ser1071Cysfs*1084 / p.Pro1380Leu |
39 | 24 | 25 | 38 | 1.56/0.3 | 2.0/2.0 | c.5461-10 T>C / c.5714+5 G>A |
40 | 18 | 26 | 36 | 1.3/1.3 | 2.0/1.3 | ND |
41 | 32 | 33 | 45 | 0.48/0.48 | 1.0/1.0 | ND |
42 | 32 | 35 | 46 | 1.3/0.0 | 3.0/1.0 | p.Cys54Tyr |
43 | 30 | 35 | 45 | 0.48/0.48 | 2.0/1.3 | ND |
44 | 15 | 41 | 49 | 1.3/1.3 | 2.0/1.3 | p.Asn965Ser |
45 | 8 | 8 | 20 | 0.78/0.78 | 1.0/1.0 | p.Thr1019Met |
46 | 10 | 11 | 23 | 1.0/1.0 | 1.0/1.0 | p.Thr1019Met |
47 | 8 | 12 | 24 | 2.0/1.56 | 1.78/1.48 | p.Cys2150Tyr |
48 | 17 | 18 | 26 | 1.0/0.78 | 1.3/1.0 | c.5461-10 T>C / p.Leu2027Phe |
49 | 8 | 21 | 33 | 1.3/1.3 | 2.0/2.0 | p.Asp574Aspfs*582 |
50 | 8 | 27 | 39 | 2.0/1.56 | 1.78/1.48 | c.5461-10 T>C |
51 | 24 | 31 | 43 | 1.18/1.18 | 1.08/1.3 | p.Arg1640Trp / p.Leu2027Phe |
52 | 11 | 31 | 42 | 1.3/1.3 | 2.0/2.0 | p.Arg1108His |
53 | 5 | 32 | 43 | 2.0/2.0 | 2.0/2.0 | c.5461-10 T>C / p.Cys2150Tyr |
54 | 5 | 32 | 43 | 2.0/2.0 | 2.0/2.0 | c.5461-10 T>C / p.Cys2150Tyr |
55 | 7 | 36 | 47 | 1.3/1.3 | 3.0/1.3 | c.5461-10 T>C / p.Cys2150Tyr |
56 | 13 | 39 | 50 | 1.25/1.56 | 3.0/3.0 | ND |
57 | 23 | 42 | 52 | 1.56/1.0 | 1.0/1.0 | p.Leu747Cysfs*787 |
58 | 40 | 43 | 54 | 0.18/0.18 | 0.78/0.78 | ND |
59 | 23 | 54 | 65 | 0.78/1.0 | 1.0/1.0 | p.Ile156Val |