Molecular Vision 2012; 18:2976-2981 <http://www.molvis.org/molvis/v18/a304>
Received 29 August 2012 | Accepted 12 December 2012 | Published 14 December 2012

The role of lysyl oxidase-like 1 DNA copy number variants in exfoliation glaucoma

Yutao Liu,1,2 Benjamin T. Whigham,1 Joshua Wheeler,1 Susan E. I. Williams,3 Robyn M. Rautenbach,4 Ari Ziskind,4 Michele Ramsay,5 Trevor R. Carmichael,3 Allison E. Ashley-Koch,1 R. Rand Allingham,2 Michael A. Hauser1,2

1Center for Human Genetics, Department of Medicine, Duke University Medical Center, Durham, NC; 2Department of Ophthalmology, Duke University Medical Center, Durham, NC; 3Division of Ophthalmology, Department of Neurosciences, University of the Witwatersrand, Johannesburg, South Africa; 4Division of Ophthalmology, Department of Surgical Sciences, Stellenbosch University, Stellenbosch, South Africa; 5Division of Human Genetics, NHLS and School of Pathology, University of Witwatersrand, Johannesburg, South Africa

Correspondence to: Michael A. Hauser, Center for Human Genetics, DUMC Box 3445 Durham, NC 27710; Phone: (919) 684-3508; FAX: (919) 684-0919; email: mike.hauser@duke.edu

Abstract

Purpose: To investigate whether DNA copy number variants (CNVs) in the lysyl oxidase-like 1 (LOXL1) gene are associated with exfoliation glaucoma (XFG) in black South Africans.

Methods: Black South African subjects with XFG and age-matched unaffected controls were recruited from the St. John Eye Hospital in Soweto (Johannesburg, South Africa) and East London Hospital Complex (Eastern Cape, South Africa) using standard clinical examination techniques. A customized array comparative genomic hybridization (aCGH) from Roche NimbleGen was designed to cover a 1.5 million base genomic region centered on the LOXL1 gene on chromosome 15. Twenty selected XFG cases were examined using this custom aCGH to identify common CNVs in the LOXL1 gene. The potential DNA copy number variants identified from aCGH were further validated using TaqMan probe-based CNV real-time PCR in a data set containing 91 XFG cases and 52 controls. The frequencies of CNVs in the LOXL1 region were compared between the XFG cases and the controls using Fisher's exact test.

Results: Several DNA CNV variants were identified in the LOXL1 genomic region using aCGH in the selected XFG cases. However, we were unable to validate these candidate CNVs using real-time PCR-based TaqMan CNV assays. There was no significant difference in the frequency of the DNA copy number variants in the LOXL1 region between the XFG cases and the controls.

Conclusions: This represents the first DNA CNV study of LOXL1 in the black South African population with XFG. Our study did not identify any significant DNA copy number alterations in the genomic region containing the LOXL1 gene. This suggests that other as yet unknown causal variants of LOXL1 or variants in other genes in linkage disequilibrium with the LOXL1 locus contribute to the genetic risk of XFG in black South Africans.

Introduction

Glaucoma is characterized by progressive loss of retinal ganglion cells, optic nerve cupping, and visual field loss [1]. Glaucoma is the second most common cause of vision loss worldwide [2,3]. Exfoliation glaucoma (XFG) is the most common identifiable secondary cause of open-angle glaucoma [4-6]. XFG is a direct result of exfoliation syndrome (XFS), a systemic condition characterized by pathological deposits of microfibrillar material within the eye and other non-ocular tissues [5-8]. Although XFG is uncommon in African Americans and virtually non-existent in West Africans [9-11], it is relatively common in black South Africans and accounts for approximately 16%–20% of glaucoma in that population [12,13]. This ancestral population can be divided into distinct linguistic groups, but they are closely related to each other genetically [14]. Recently, we reported the strong association of coding variants in the lysyl oxidase-like 1 (LOXL1) gene with XFG in this population [15,16]. Interestingly, the major allele of rs3825942 that confers increased risk of XFG in all non-African populations is associated with reduced risk in South Africans [4,15-17]. This finding, in addition to the fact that the two associated variants, R141L (rs1048661) and G153D (rs3825942) [17-21], have been shown to have no effect on the amine oxidase activity of LOXL1 protein, suggests that these LOXL1 variants tag the true functional variants [22]. These functional variants remain unknown [4], but may exert an effect through disruption of LOXL1 regulation [4,8,23].

DNA copy number variants (CNVs) have been shown to be responsible for phenotypic variations [24] as well as human genetic disorders, including primary open-angle glaucoma [25,26], normal tension glaucoma [27], and age-related macular degeneration [28-31]. CNVs can change the copy number of entire genes, or may duplicate or delete regulatory elements that affect the expression level of one or more genes [24]. Given the reported alteration of LOXL1 expression in XFS and XFG [8,20,32,33], it is important to examine the LOXL1 region for potential disease-associated CNV changes, studies that have not been performed to date in XFS or XFG.

To identify potential CNVs encompassing the LOXL1 gene, we applied a two-stage approach. First, we designed a custom comparative genomic hybridization array targeting the LOXL1 locus in a set of 20 XFG cases. Second, we used TaqMan-based real-time PCR in 91 XFG cases and 52 controls to validate candidate CNV changes. This is the first report to examine the roles of LOXL1-related CNVs in XFG, especially in black South African individuals.

Methods

Study participants

This study adhered to the tenets of the Declaration of Helsinki. The research protocol was approved by all participating universities including the University of the Witwatersrand Human Research Ethics Committee and Stellenbosch University Health Research Ethics Committee. Black Southern African participants with clinically diagnosed XFG and unaffected control subjects were recruited from the St. John Eye Hospital in Soweto (Johannesburg, South Africa) and East London Hospital Complex (Eastern Cape, South Africa), and have been previously described [15,16]. Written informed consent was obtained from all participants. Ethnic affiliation was established by the home language of participants and that of their parents and grandparents. All participants underwent a standardized detailed ophthalmic examination by an ophthalmologist (S.E.I.W. or R.M.R.). The examination included measurement of intraocular pressure (IOP) by applanation, slit-lamp biomicroscopy, gonioscopy, dilated pupil examination of the lens and fundus, and visual field testing. Subjects with XFG were defined as those with clinical evidence of exfoliation material on the pupil margin, anterior lens surface, and the presence of glaucomatous optic neuropathy and visual field loss. Gender- and ethnicity-matched South African subjects with normal IOP and normal anterior segment and optic nerve examinations were recruited as control subjects. Control subjects were older than 55 years of age to reduce the possibility of misdiagnosis.

DNA copy number analysis with array comparative genomic hybridization

Genomic DNA was extracted using a salting out procedure from nucleated cells from the venous blood samples from all subjects [34]. A total of 91 XFG cases and 52 controls were included in our study. Twenty XFG cases were selected for Roche NimbleGen (Madison, WI) array CGH analysis using custom designed arrays targeting the genomic region around the LOXL1 gene. The array CGH was performed as described previously [25]. Briefly, 500 ng high-quality undegraded genomic DNA was labeled with Cy3, and the reference sample from HapMap was labeled with Cy5. The hybridization was done according to the standard procedures recommended by the manufacturer. The arrays were scanned on an Axon4100A scanner (Molecular Devices, Sunnyvale, CA). Scanned TIFF images were processed using NimbleScan 2.6 software. The CNV calls were determined using segMNT algorithm in NimbleScan and were visualized using SignalMap 1.8 from Roche NimbleGen. Briefly, blood was obtained using peripheral venipuncture. The blood was anticoagulated and lysed. Cellular protein was salted out of solution using NaCl. Genomic DNA, still in solution, was removed, precipitated in ethanol, and resuspended in water.

Copy number variant validation and screening with real-time polymerase chain reaction

To validate the CNVs identified from the array CGH experiments, we selected three candidate CNVs in the LOXL1 genomic region. We used TaqMan Copy Number Assays from ABI (Applied Biosystems Inc., Carlsbad, CA) on the ABI Prism 7900HT Sequence Detection System for validation and further screening in all XFG cases and controls. A Vic-labeled Copy Number Assay for RNase P was selected as an internal control as it performed in the same reaction with gene-specific assays. Four Copy Number Assays were ordered from ABI. Hs01538855_cn, HS01338991_cn, and HS00942919_cn assays to target the exon 1 coding region, which would be used to validate any CNVs overlapping with the LOXL1 exon 1. The Hs03894435_cn assay targets intron 6 of the LOXL1 gene, which could be used to validate any CNVs affecting the LOXL1 intron 6. Each sample was assayed with four replicates by using 10 ng DNA in each reaction in a 384-well format. The CNV calls were generated with SDS software and CopyCaller from ABI (Applied Biosystems). A known CEPH sample was used as a reference for a copy number of 2. This CEPH sample had been previously confirmed with a single nucleotide polymorphism microarray to contain two copies of candidate CNV regions. To make CNV calls in CopyCaller software, a confidence score of greater than 0.95 was required with four replicates.

Results

The cases in this study had exfoliation-related open-angle glaucoma [15,16]. Our data set included 91 black South Africans with XFG and 52 black South African controls. The mean age at enrollment for the cases and controls was 75.01±8.98 and 74.5±8.49 years, respectively. The percentage of women in our data set was 41% in the cases and 56% in the controls.

We designed a custom array CGH based on the CGH 4X72k array from Roche NimbleGen. The size of the area examined in the LOXL1 genomic region was 1.5 Mb (chr15:71,250,002–72,749,772, Hg18) and was centered on the LOXL1 gene. More than 42,000 probes were used in the analysis (Figure 1). These designed probes ranged from 50 to 72 bp to optimize their thermodynamic properties. The mean interval between adjacent probes was 27 bp. No probes were designed within genomic regions with repeated sequences, represented by the gaps in Figure 1. After quality control and data analysis, we identified three potential DNA copy number variations in the LOXL1 region (Table 1), which met the following criteria: 1) at least ten probes per candidate CNV, 2) genomic size ≥1 kb, and 3) good-quality metrics of array CGH hybridization. Two CNV candidates, CNV1 and CNV3, overlapped with exon 1 and intron 1 of the LOXL1 gene while CNV2 overlapped exons 3–7 of the LOXL1 gene (Figure 2).

To validate the DNA copy number changes, we selected four predesigned DNA copy number assays, whose relative genomic locations are shown in Figure 2. The first three assays (Hs01538855_cn, Hs01338991_cn, and Hs00942919_cn) targeted the first exon and were used to validate CNV1 and CNV3 identified in the aCGH experiment. The assay Hs03894435_cn targeted intron 6 and was used to validate CNV2 identified in that experiment. TaqMan-based real-time PCR reactions in quadruplicate failed to validate any of the three identified DNA duplication events. Instead, we found that a single control sample had a DNA duplication (Figure 3). This DNA duplication was confirmed with all four copy number assays (data not shown), indicating that this specific DNA duplication in this unaffected individual covered the entire LOXL1 genomic region. The frequency of this duplication was not significantly different between the XFG cases and the controls (Fisher’s exact test, two-tailed p value=0.36).

Discussion

This study is the first study of LOXL1 copy number variants in exfoliation glaucoma conducted in black South African subjects. We applied a comprehensive approach using array CGH and copy number real-time PCR assays. We failed to identify LOXL1 CNVs that were significantly different between the XFG cases and the controls.

Although our initial CGH array analysis identified three potential CNVs in the XFG cases, these CNVs were not validated with TaqMan probe-based copy number real-time PCR assays. The lack of validation might be related to the DNA probes selected in the CGH array. Many of the probes on the CGH array overlapped since the median interval between adjacent probes was 27 bp and the length of the probes ranged from 50 bp to 72 bp. This probe overlap may have interfered with the DNA hybridization and led to false-positive CNV findings, which might be the main reason for our failure to validate the identified CNVs. The CNV-specific real-time PCR assays were performed with four replicates for every sample with an internal reference assay targeting the RNase P gene. This robust design enabled more reliable and replicable real-time PCR results. Therefore, validation of CNV-specific real-time PCR was considered the main criterion in our study. This study underscores the importance of validation when studying CNVs. The array CGH covers a relatively large genomic region (>1 kb at least) while real-time PCR-based CNV assays target a relatively small genomic region (<100 bp).

The black South African subjects in this study are uniquely valuable and confer extraordinary power for studying XFG [15,16]. We have reported a strong association (p=5.2 × 10−13) between XFG disease status and the A allele of the variant rs3825942 despite a small data set [15,16]. This study examined the potential role of LOXL1 copy number variants in this special population. Our study indicated a limited contribution of the potential common LOXL1 copy number variants to XFG in black South African individuals. It will be necessary to replicate our negative findings in a larger sample of black South Africans and other populations.

This study was limited by the relatively small sample size. The 20 selected cases in the array CGH experiment were used to identify common CNVs (≥2.5% of allele frequency) in the LOXL1 genomic region. Our approach may have missed potential CNVs with lower allele frequency. Our real-time PCR-based experimental finding was also limited by the relatively small number of cases and controls from South Africa. Previously, we reported the significant genetic risk associated with LOXL1 coding variants rs1048661 and rs3825942 in the black South African population with odds ratio 23.2 and 0.092, respectively [16]. Assuming an CNV allele frequency of 0.05, we calculated using Quanto software [35] that our case-control data set has more than 99% power to detect either risk. More XFG cases and controls from South Africa could be integrated to improve the statistical power in the future.

In summary, we conducted a two-stage DNA copy number study in the LOXL1 genomic region with XFG in a black South African population. We did not detect a significant contribution of common CNV in LOXL1 to the increased risk of XFG. This study does not rule out the possible contribution of CNVs in other populations with XFS or XFG.

Acknowledgments

The authors thank all the study participants at the St. John Eye Hospital and the East London Hospital Complex; Fahmida Essop, Punita Pitamber, and the staff of the Molecular Genetics Laboratory at the NHLS in Braamfontein, for their assistance with DNA extractions; Dr. Soraya Bardien of Stellenbosch University for DNA extraction; and the study staff at the Duke Center for Human Genetics for their assistance. This publication was made possible (in part) by a grant from the Carnegie Corporation of New York. This work is also supported by National Institutes of Health (NIH) grants R01EY013315 (MAH), R01EY019126 (MAH), R03EY014939 (RRA), and R01EY015543 (RRA), the Glaucoma Research Foundation (YL), the National Glaucoma Research of American Health Assistance Foundation (YL), The Glaucoma Foundation (YL). This research was also supported by Duke University’s CTSA grant TL1RR024126 from NCRR/NIH. The statements made and views expressed are, however, solely the responsibility of the authors.

References

  1. Wiggs JL. Genetic etiologies of glaucoma. Arch Ophthalmol. 2007; 125:30-7. [PMID: 17210849]
  2. Quigley HA, Broman AT. The number of people with glaucoma worldwide in 2010 and 2020. Br J Ophthalmol. 2006; 90:262-7. [PMID: 16488940]
  3. Liu Y, Allingham RR. Genetics of Glaucoma. Encyclopeida of Life Sciences (ELS): John Wiley & Sons 2010
  4. Liu Y, Allingham RR. Molecular genetics in glaucoma. Exp Eye Res. 2011; 93:331-9. [PMID: 21871452]
  5. Ritch R. The management of exfoliative glaucoma. Prog Brain Res. 2008; 173:211-24. [PMID: 18929111]
  6. Schlötzer-Schrehardt U, Naumann GO. Ocular and systemic pseudoexfoliation syndrome. Am J Ophthalmol. 2006; 141:921-37. [PMID: 16678509]
  7. Ritch R, Schlotzer-Schrehardt U. Exfoliation (pseudoexfoliation) syndrome: toward a new understanding. Proceedings of the First International Think Tank. Acta Ophthalmol Scand. 2001; 79:213-7. [PMID: 11284770]
  8. Schlötzer-Schrehardt U. Molecular pathology of pseudoexfoliation syndrome/glaucoma–new insights from LOXL1 gene associations. Exp Eye Res. 2009; 88:776-85. [PMID: 18809397]
  9. Ntim-Amponsah CT, Amoaku WM, Ofosu-Amaah S, Ewusi RK, Idirisuriya-Khair R, Nyatepe-Coo E, Adu-Darko M. Prevalence of glaucoma in an African population. Eye (Lond). 2004; 18:491-7. [PMID: 15131680]
  10. Cashwell LF, , Jr Shields MB. Exfoliation syndrome. Prevalence in a southeastern United States population. Arch Ophthalmol. 1988; 106:335-6. [PMID: 3345150]
  11. Ball SF. Exfoliation syndrome prevalence in the glaucoma population of South Louisiana. Acta Ophthalmol Suppl. 1988; 184:93-8. [PMID: 2853928]
  12. Rotchford AP, Kirwan JF, Muller MA, Johnson GJ, Roux P. Temba glaucoma study: a population-based cross-sectional survey in urban South Africa. Ophthalmology. 2003; 110:376-82. [PMID: 12578784]
  13. Rotchford AP, Johnson GJ. Glaucoma in Zulus: a population-based cross-sectional survey in a rural district in South Africa. Arch Ophthalmol. 2002; 120:471-8. [PMID: 11934321]
  14. Tishkoff SA, Reed FA, Friedlaender FR, Ehret C, Ranciaro A, Froment A, Hirbo JB, Awomoyi AA, Bodo JM, Doumbo O, Ibrahim M, Juma AT, Kotze MJ, Lema G, Moore JH, Mortensen H, Nyambo TB, Omar SA, Powell K, Pretorius GS, Smith MW, Thera MA, Wambebe C, Weber JL, Williams SM. The genetic structure and history of Africans and African Americans. Science. 2009; 324:1035-44. [PMID: 19407144]
  15. Rautenbach RM, Bardien S, Harvey J, Ziskind A. An investigation into LOXL1 variants in black South African individuals with exfoliation syndrome. Arch Ophthalmol. 2011; 129:206-10. [PMID: 21320968]
  16. Williams SE, Whigham BT, Liu Y, Carmichael TR, Qin X, Schmidt S, Ramsay M, Hauser MA, Allingham RR. Major LOXL1 risk allele is reversed in exfoliation glaucoma in a black South African population. Mol Vis. 2010; 16:705-12. [PMID: 20431720]
  17. Chen H, Chen LJ, Zhang M, Gong W, Tam PO, Lam DS, Pang CP. Ethnicity-based subgroup meta-analysis of the association of LOXL1 polymorphisms with glaucoma. Mol Vis. 2010; 16:167-77. [PMID: 20142848]
  18. Thorleifsson G, Magnusson KP, Sulem P, Walters GB, Gudbjartsson DF, Stefansson H, Jonsson T, Jonasdottir A, Stefansdottir G, Masson G, Hardarson GA, Petursson H, Arnarsson A, Motallebipour M, Wallerman O, Wadelius C, Gulcher JR, Thorsteinsdottir U, Kong A, Jonasson F, Stefansson K. Common sequence variants in the LOXL1 gene confer susceptibility to exfoliation glaucoma. Science. 2007; 317:1397-400. [PMID: 17690259]
  19. Fingert JH, Alward WL, Kwon YH, Wang K, Streb LM, Sheffield VC, Stone EM. LOXL1 mutations are associated with exfoliation syndrome in patients from the midwestern United States. Am J Ophthalmol. 2007; 144:974-5. [PMID: 18036875]
  20. Hewitt AW, Sharma S, Burdon KP, Wang JJ, Baird PN, Dimasi DP, Mackey DA, Mitchell P, Craig JE. Ancestral LOXL1 variants are associated with pseudoexfoliation in Caucasian Australians but with markedly lower penetrance than in Nordic people. Hum Mol Genet. 2008; 17:710-6. [PMID: 18037624]
  21. Pasutto F, Krumbiegel M, Mardin CY, Paoli D, Lammer R, Weber BH, Kruse FE, Schlotzer-Schrehardt U, Reis A. Association of LOXL1 common sequence variants in German and Italian patients with pseudoexfoliation syndrome and pseudoexfoliation glaucoma. Invest Ophthalmol Vis Sci. 2008; 49:1459-63. [PMID: 18385063]
  22. Kim S, Kim Y. Variations in LOXL1 associated with exfoliation glaucoma do not affect amine oxidase activity. Mol Vis. 2012; 18:265-70. [PMID: 22328822]
  23. Lee RK. The molecular pathophysiology of pseudoexfoliation glaucoma. Curr Opin Ophthalmol. 2008; 19:95-101. [PMID: 18301281]
  24. Zhang F, Gu W, Hurles ME, Lupski JR. Copy number variation in human health, disease, and evolution. Annu Rev Genomics Hum Genet. 2009; 10:451-81. [PMID: 19715442]
  25. Liu Y, Gibson J, Wheeler J, Kwee LC, Santiago-Turla CM, Akafo SK, Lichter PR, Gaasterland DE, Moroi SE, Challa P, Herndon LW, Girkin CA, Budenz DL, Richards JE, Allingham RR, Hauser MA. GALC deletions increase the risk of primary open-angle glaucoma: the role of Mendelian variants in complex disease. PLoS ONE. 2011; 6:e27134 [PMID: 22073273]
  26. Davis LK, Meyer KJ, Schindler EI, Beck JS, Rudd DS, Grundstad AJ, Scheetz TE, Braun TA, Fingert JH, Alward WL, Kwon YH, Folk JC, Russell SR, Wassink TH, Sheffield VC, Stone EM. Copy number variations and primary open-angle glaucoma. Invest Ophthalmol Vis Sci. 2011; 52:7122-33. [PMID: 21310917]
  27. Fingert JH, Robin AL, Stone JL, Roos BR, Davis LK, Scheetz TE, Bennett SR, Wassink TH, Kwon YH, Alward WL, Mullins RF, Sheffield VC, Stone EM. Copy number variations on chromosome 12q14 in patients with normal tension glaucoma. Hum Mol Genet. 2011; 20:2482-94. [PMID: 21447600]
  28. Meyer KJ, Davis LK, Schindler EI, Beck JS, Rudd DS, Grundstad AJ, Scheetz TE, Braun TA, Fingert JH, Alward WL, Kwon YH, Folk JC, Russell SR, Wassink TH, Stone EM, Sheffield VC. Genome-wide analysis of copy number variants in age-related macular degeneration. Hum Genet. 2011; 129:91-100. [PMID: 20981449]
  29. Kubista KE, Tosakulwong N, Wu Y, Ryu E, Roeder JL, Hecker LA, Baratz KH, Brown WL, Edwards AO. Copy number variation in the complement factor H-related genes and age-related macular degeneration. Mol Vis. 2011; 17:2080-92. [PMID: 21850184]
  30. Liu MM, Agron E, Chew E, Meyerle C, Ferris FL, , 3rd Chan CC, Tuo J. Copy number variations in candidate genes in neovascular age-related macular degeneration. Invest Ophthalmol Vis Sci. 2011; 52:3129-35. [PMID: 21273533]
  31. Sawitzke J, Im KM, Kostiha B, Dean M, Gold B. Association assessment of copy number polymorphism and risk of age-related macular degeneration. Ophthalmology. 2011; 118:2442-6. [PMID: 21856016]
  32. Khan TT, Li G, Navarro ID, Kastury RD, Zeil CJ, Semchyshyn TM, Moya FJ, Epstein DL, Gonzalez P, Challa P. LOXL1 expression in lens capsule tissue specimens from individuals with pseudoexfoliation syndrome and glaucoma. Mol Vis. 2010; 16:2236-41. [PMID: 21139690]
  33. Schlötzer-Schrehardt U, Pasutto F, Sommer P, Hornstra I, Kruse FE, Naumann GO, Reis A, Zenkel M. Genotype-correlated expression of lysyl oxidase-like 1 in ocular tissues of patients with pseudoexfoliation syndrome/glaucoma and normal patients. Am J Pathol. 2008; 173:1724-35. [PMID: 18974306]
  34. Miller SA, Dykes DD, Polesky HF. A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988; 16:1215 [PMID: 3344216]
  35. Gauderman WJ. Sample size requirements for matched case-control studies of gene-environment interaction. Stat Med. 2002; 21:35-50. [PMID: 11782049]