Molecular Vision 2017; 23:660-665 <http://www.molvis.org/molvis/v23/660>
Received 30 May 2017 | Accepted 26 September 2017 | Published 28 September 2017

Genome-wide linkage and association analysis of primary open-angle glaucoma endophenotypes in the Norfolk Island isolate

Elizabeth Matovinovic,1 Pik Fang Kho,1 Rodney A. Lea,1 Miles C. Benton,1 David A. Eccles,1 Larisa M. Haupt,1 Alex W. Hewitt,2,3 Justin C. Sherwin,2 David A. Mackey,3,4 Lyn R. Griffiths1

1Genomics Research Centre, Institute of Health and Biomedical Innovation, School of Biomedical Sciences, Queensland University of Technology, Brisbane, Australia; 2Centre for Eye Research Australia & Royal Victorian Eye and Ear Hospital, University of Melbourne, East Melbourne, Victoria, Australia; 3Menzies Institute for Medical Research, School of Medicine, University of Tasmania, Hobart, Australia; 4Lions Eye Institute, Centre for Ophthalmology and Visual Science, University of Western Australia, Perth, Australia

Correspondence to: Lyn Griffiths, Institute of Health and Biomedical Innovation, Queensland University of Technology 60 Musk Ave, Kelvin Grove, QLD, 4059, Australia; Phone: +617 3138 6102; FAX: +617 3138 6039; email: lyn.griffiths@qut.edu.au

Abstract

Purpose: Primary open-angle glaucoma (POAG) refers to a group of heterogeneous diseases involving optic nerve damage. Two well-established risk factors for POAG are elevated intraocular pressure (IOP) and a thinner central corneal thickness (CCT). These endophenotypes exhibit a high degree of heritability across populations. Large-scale genome-wide association studies (GWASs) of outbred populations have robustly implicated several susceptibility gene variants for both IOP and CCT. Despite this progress, a substantial amount of genetic variance remains unexplained. Population-specific variants that might be rare in outbred populations may also influence POAG endophenotypes. The Norfolk Island population is a founder-effect genetic isolate that has been well characterized for POAG endophenotypes. This population is therefore a suitable candidate for mapping new variants that influence these complex traits.

Methods: Three hundred and thirty participants from the Norfolk Island Eye Study (NIES) core pedigree provided DNA. Ocular measurements of CCT and IOP were also taken for analysis. Heritability analyses and genome-wide linkage analyses of short tandem repeats (STRs) were conducted using SOLAR. Pedigree-based GWASs of single-nucleotide polymorphisms (SNPs) were performed using the GenABEL software.

Results: CCT was the most heritable endophenotype in this cohort (h2 = 0.77, p = 6×10−6), while IOP showed a heritability of 0.39 (p = 0.008). A genome-wide linkage analysis of these POAG phenotypes identified a maximum logarithm of the odds (LOD) score of 1.9 for CCT on chromosome 20 (p = 0.0016) and 1.3 for IOP on chromosome 15 (p = 0.0072). The GWAS results revealed a study-wise significant association for IOP at rs790357, which is located within DLG2 on chr11q14.1 (p = 1.02×10−7). DLG2 is involved in neuronal signaling and development, and while it has not previously been associated with IOP, it has been associated with myopia. An analysis of 12 known SNPs for IOP showed that rs12419342 in RAPSN on chromosome 11 was nominally associated in Norfolk Island (NI; p = 0.0021). For CCT, an analysis of 26 known SNPs showed rs9938149 in BANP-ZNF469 on chromosome 16 was nominally associated in NI (p = 0.002).

Conclusions: These study results indicate that CCT and IOP exhibit a substantial degree of heritability in the NI pedigree, indicating a genetic component. A genome-wide linkage analysis of POAG endophenotypes did not reveal any major effect loci, but the GWASs did implicate several known loci, as well as a potential new locus in DLG2, suggesting a role for neuronal signaling in development in IOP and perhaps POAG. These results also highlight the need to target rarer variants via whole genome sequencing in this genetic isolate.

Introduction

Glaucoma refers to a group of heterogeneous diseases causing optic nerve fiber damage. Primary open-angle glaucoma (POAG) is the most common type of glaucoma in most populations. Elevated intraocular pressure (IOP), defined as >22 mmHg, is a principal risk factor for POAG and can cause optic nerve damage [1]. Heritability estimates (h2) for IOP range from 0.29 to 0.67 [1].

Genome-wide association studies (GWASs) have been successful in mapping genetic variants for IOP. On a large scale, a GWAS meta-analysis was performed on 18 cohorts comprising 35,296 individuals. This study showed robust associations among 12 different single-nucleotide polymorphisms (SNPs) representing nine different genes/loci. Notably, four of these genes (TMCO1, CAV7, ABCA1, and GAS7) were also associated with POAG, highlighting the key role of this endophenotype in end-stage ocular disease [2].

Linkage mapping of genes for IOP identified seven regions of interest on chromosomes 2, 5, 6, 7, 12, 15, and 19 [3]. Chromosomes 6 and 13 were implicated in another IOP linkage study [4]. Evidence for the linkage of a maximum IOP has been mapped to chromosome 10 near region 10p15-p14, which is near another locus involved in up to 17% of low-tension glaucoma pedigrees [5]. Chromosome 2p15-p16.3 is implicated in IOP [6], while 2cen-q13 is found in normal-tension glaucoma, with onset at over 50 years; chromosome 7q35-q36 has also been found to be associated with IOP [7].

A thinner central corneal thickness (CCT) is also a risk factor for POAG and has been shown to have high heritability across populations (0.65–0.95) [1]. A large GWAS of >20,000 individuals confirmed the previously identified 10 different loci and revealed the presence of 16 new loci associated with CCT [8]. Of particular interest is the FNDC3B gene on chromosome 3, which was also associated with POAG in this study and with IOP in the GWAS of Hysi et al. [2]. This suggests that variants of the FNDC3B gene may be conferring a pleiotropic effect on POAG endophenotypes and end-stage clinical diseases.

Despite this progress, a substantial amount of genetic variance remains unexplained. Population-specific variants that might be rare in outbred populations may also influence POAG endophenotypes. Founder-effect genetic isolates may offer some advantages for gene mapping because of such features as reduced genetic heterogeneity and an extensive familial structure. Linkage and association analyses of ocular phenotypes in a Dutch population isolate revealed new linkage regions for optic disc morphology at 20p13 (SIRPA and RNF24/PANK2 loci) and 2q33-q34 (IGFBP2 locus). In addition, significant and suggestive GWAS signals at the previously identified RERE, LRP1B, CDC7, TGFBR3, and ATOH7 loci were replicated in this isolate [9].

The Norfolk Island (NI) population, located off the east coast of Australia, is a genetic isolate with well-defined familial lineages dating back 12 generations to the originating founders, 12 Tahitian women and six European men. As part of the Norfolk Island Eye Study (NIES), the prevalence of chronic ocular diseases has previously been reported [10]. In general, the prevalence of blindness and visual impairment was observed to be lower than that in mainland Australia. However, the prevalence of glaucoma was higher in NI at 6% compared to 3% in mainland Australia [10]. Differences between NI pedigree and non-pedigree groups were observed for several POAG risk factor traits, suggesting a possible genetic influence some traits [11]. In an effort to understand better the genetic basis of POAG, we focused on two primary risk factor phenotypes—elevated IOP and a thinner CCT—by conducting heritability, genome-wide linkage and association studies of a large core pedigree from the NI isolate.

Methods

Subjects and phenotyping

The study sample consisted of 330 participants from the NIES. These individuals belong to a single large NI pedigree of 6,537 members dating back to the original Bounty mutineers and Tahitian founders. Pedigree construction and characteristics have been well described in previous papers [12-14]. All participants underwent a full eye examination as part of the NIES, which was designed to investigate the genetics and epidemiology of ocular disease traits in this population. Comprehensive ocular measurements for IOP and CCT were taken (see [11] for details).

Genotyping

Genotypes for all 330 individuals were obtained for 400 microsatellite (STR) markers with an average heterozygosity level of 76% spaced at 10cM. Chromosomal marker maps were obtained from the Marshfield Centre for Medical Genetics. Cytobands for markers were obtained from the University of California Santa Cruz (UCSC) Genome Browser.

Individuals were also genotyped for a genome-wide panel of SNPs using Illumina Infinium High Density (HD) Human610-Quad BeadChip version 1, according to the manufacturer’s instructions. Samples were scanned on the Illumina Bead Array 500GX Reader. Raw data were obtained using the Illumina Bead Scan image data acquisition software (version 2.3.0.13). Raw data from the Illumina data files were then SNP genotyped in R using the CRLMM package [15]. Genotype data then underwent initial quality control routines using PLINK [16]. SNPs were filtered based on a minor allele frequency of >0.01, a call rate of >0.95, and a Hardy–Weinberg equilibrium testing p value of >10−5. After this initial quality control, 590,603 SNPs were exported from PLINK and imported into the CRAN package GenABEL [17]. Further filtering (including Mendelian inheritance violations and sex-checking based on available X and Y markers) in GenABEL led to the reduction of the SNP set to a total of ~480,000; this included the removal of both X and Y chromosome SNPs after gender checking, as well as the removal of mitochondrial and XY SNPs.

Statistical analyses

Heritability analyses were conducted using a variance components-based methodology implemented in the Sequential Oligonucleotide Linkage Analysis Routines (SOLAR) version 4.0.6 software package [18]. Heritability estimates (h2) were calculated as the ratio of the trait variance explained by additive polygenic effects to the total phenotypic variance of the trait [19]. The applied polygenic model assumes an infinite number of genetic factors, each with a small additive effect contributing to the trait variance. Estimates were screened for covariate effects of sex, age, age*sex, age^2, age^2*sex, height, ultraviolet autofluorescence, glaucoma status, visual acuity pinhole, sphere, cylinder, axis, kvalue-v, kvalue axis, CCT, ACD, IOP, pterygium, eye color, visual acuity, axial length, outdoors, and kvalue-h. Glaucoma was treated as a continuous variable using “EnableDiscrete 0” to circumvent convergence errors due to discrete trait modeling in SOLAR. Covariates with p values less than or equal to 0.10 were retained in the final model. SOLAR requires kurtosis to be <0.8 and a standard deviation of <0.5 to proceed with the analysis. The SOLAR bivariate analysis was used to calculate genotypic, environmental, and phenotypic correlations between metrics. Genome-wide linkage mapping of adjusted eye metric traits was conducted with SOLAR using a multipoint analysis.

A pedigree-based GWAS analysis of all heritable endophenotypes was batched using custom R scripts in GenABEL [17]. A correction was made for the relatedness and admixture inherent in the NI population using the polygenic model with age and sex interactions, as well as the genetic structure —the top two genomic principal components of the complete SNP set, as calculated by KING [20]. For an association analysis, the mmscore function implemented in GenABEL was used. This function represents a mixed-model approximation analysis of the association between a trait and genetic polymorphism(s) and it is specifically designed for association testing in samples of related individuals. This allows for per-SNP association testing using a mixed-model polygenic approach. After adjusting for linkage disequilibrium (LD) among SNPs using Meff [21], the study-wide significance was set at p = 1.84×10−7. This significance threshold has been used in previous GWAS publications of NI [22,23]. GWAS Manhattan plots were generated for each trait association using a custom modified version of the GenABEL plot.scan.gwaa function (for all Manhattan plots, see Supplemental Material).

Previous studies of CVD risk have determined the current trimmed pedigree structure and size had >80% power to detect the heritability of phenotypes whose variation is partially attributable to additive genetic effects and to detect major effect loci via linkages and associations as that are statistically significant at accepted thresholds [24,25]. Evidence that the major effect loci exist for complex quantitative traits in this NI pedigree has also been published [22,23].

Results

Heritability analyses showed that CCT was the most heritable trait, with an unadjusted h2 of 0.85 (p = 1.5×10−6). Adjusted heritability for CCT decreased to h2 = 0.77 (p = 5.7×10−6), with the inclusion of IOP, pterygium, and sphere as significant covariates. The adjusted IOP heritability was 0.39 (p = 0.008) with significant covariates anterior chamber depth (ACD) and ultraviolet autofluorescence, a biomarker for sun exposure accounting for 1% of the trait variation. Overall, the heritability values for CCT and IOP in the NI population fell within the range of heritability estimates among worldwide populations [1]. Importantly, these results indicate a substantial influence of genetic factors on these traits in the NI pedigree.

A genome-wide linkage and association analysis was performed on the core NI pedigree to identify loci for CCT and IOP. Table 1 summarizes logarithm of the odds (LOD) scores for the traits adjusted for significant covariates. There were no statistically significant linkage peaks (i.e., LOD ≥ 3). The highest linkage peak was for CCT on chromosome 20 (LOD = 1.9, p = 0.0016). Table 2 summarizes the GWAS results. For IOP, there was a statistically significant association at rs790357 located in DLG2 on chromosome 11q14.1 (odds ratio [OR] = 1.70, p = 1.02×10−7). This top-ranking SNP was part of a cluster of seven SNPs forming an LD block spanning a 16.6-kb region (Pall< p = 1.00×10−6). There were also several suggestive hits on chromosomes 10 and 20 for IOP. For CCT, there were no association peaks that met the genome-wide level of statistical significance, but there were several suggestive hits, with the most significant being for rs13095933 on chromosome 3 (p = 2.86×10−6). We also examined the previously associated index SNPs for IOP (n = 12) and CCT (n = 26), as reported in the large-scale GWASs of Hysi et al. (2014) and Lu et al. (2013), respectively [2,8]. These results showed that for IOP, only one SNP—rs12419342 in RAPSN on chromosome 11—was statistically significant at the nominal level (p = 0.0021). For CCT, the only statistically significant SNP was rs9938149 in BANP-ZNF469 on chromosome 16 (p = 0.002).

Discussion

Herein, we examined the genetic isolate of NI and performed heritability and genome-wide linkage and association studies of phenotypes associated with POAG, CCT, and IOP. Overall, the heritability estimates were consistent with previous values, and the indicated genetic factors influence these traits. A linkage analysis for IOP showed that the highest LOD Score occurred on chromosome 15 (LOD 1.3), where a peak was also found in the Beaver Dam Eye Study [4]. Linkage results for CCT did not overlap any previously published peaks. Although a previous study indicated a linkage to chromosome 3 for IOP [7], in this study, chromosome 3 was associated with CCT.

A GWAS of ocular endophenotypes in this study identified a statistically significant association at a SNP representing a 16.6-kb locus harboring the Disks large homolog 2 (DLG2) gene. DLG2, also known as channel-associated protein of synapse-110 (chapsyn-110) or postsynaptic density protein 93 (PSD-93), is involved in neuronal signaling and development. Specifically, the protein forms a heterodimer with a related family member that may interact at postsynaptic sites to form a multimeric scaffold for the clustering of receptors, ion channels, and associated signaling proteins. DLG2 is expressed in many neuronal and immunological tissues and the retina [26]. DLG2 variants have not previously been associated with IOP or POAG. However, a GWAS performed by Kiefer et al. identified a significant and replicated association of a DLG2 SNP (rs2155413) and myopia in a large cohort of >400,000 individuals [26]. Myopia is associated with intra-ocular pressure and is a risk factor for POAG [26].

Our GWAS also revealed several suggestive minor effect peaks in the NI population that may represent new loci requiring replication. Surprisingly, association testing of the known SNPs for IOP and CCT did not indicate confirmation in NI for most loci. This may reflect the modest sample size of NI, a different genetic architecture for these traits in NI, or both. In addition, the fact that CCT does not exhibit a major effect locus indicates it is of a greater than expected complexity (e.g., genetic heterogeneity) in the NI isolate.

This study indicates that the POAG endophenotypes CCT and IOP exhibit a substantial degree of heritability in the NI pedigree, indicating a genetic component. A genome-wide linkage analysis of these traits did not reveal any major effect loci, but a GWAS did implicate several known loci, as well as a potential new locus in DLG2, suggesting a role for neuronal signaling and development in IOP and perhaps POAG. Given the substantial degree of heritability, further genetic studies are warranted, and future endeavors implementing approaches such as whole genome sequencing may help identify important, rare, or private variants of these ocular traits.

Acknowledgments

The authors acknowledge the Ian Potter Foundation, John and Jenny Corbett and the National Health and Medical Research Foundation (GNT1058806) for funding contributions to this research. The authors acknowledge the Ian Potter Foundation, Corbett and the National Health and Medical Research Foundation (no. 10588096) for funding contributions to this research.

References

  1. Sanfilippo PG, Medland SE, Hewitt AW, Kearns LS, Ruddle JB, Sun C, Hammond CJ, Young TL, Martin NG, Mackey DA. Ophthalmic phenotypes and the representativeness of twin data for the general population. Invest Ophthalmol Vis Sci. 2011; 52:5565-72. [PMID: 21498610]
  2. Hysi PG, Cheng CY, Springelkamp H, Macgregor S, Bailey JN, Wojciechowski R, Vitart V, Nag A, Hewitt AW, Höhn R, Venturini C, Mirshahi A, Ramdas WD, Thorleifsson G, Vithana E, Khor CC, Stefansson AB, Liao J, Haines JL, Amin N, Wang YX, Wild PS, Ozel AB, Li JZ, Fleck BW, Zeller T, Staffieri SE, Teo YY, Cuellar-Partida G, Luo X, Allingham RR, Richards JE, Senft A, Karssen LC, Zheng Y, Bellenguez C, Xu L, Iglesias AI, Wilson JF, Kang JH, van Leeuwen EM, Jonsson V, Thorsteinsdottir U, Despriet DD, Ennis S, Moroi SE, Martin NG, Jansonius NM, Yazar S, Tai ES, Amouyel P, Kirwan J, van Koolwijk LM, Hauser MA, Jonasson F, Leo P, Loomis SJ, Fogarty R, Rivadeneira F, Kearns L, Lackner KJ, de Jong PT, Simpson CL, Pennell CE, Oostra BA, Uitterlinden AG, Saw SM, Lotery AJ, Bailey-Wilson JE, Hofman A, Vingerling JR, Maubaret C, Pfeiffer N, Wolfs RC, Lemij HG, Young TL, Pasquale LR, Delcourt C, Spector TD, Klaver CC, Small KS, Burdon KP, Stefansson K, Wong TY, BMES GWAS Group. NEIGHBORHOOD Consortium; Wellcome Trust Case Control Consortium 2, Viswanathan A, Mackey DA, Craig JE, Wiggs JL, van Duijn CM, Hammond CJ, Aung T. Genome-wide analysis of multi-ancestry cohorts identifies new loci influencing intraocular pressure and susceptibility to glaucoma. Nat Genet. 2014; 46:1126-30. [PMID: 25173106]
  3. Duggal P, Klein AP, Lee KE, Klein R, Klein BE, Bailey-Wilson JE. Identification of novel genetic loci for intraocular pressure: a genomewide scan of the Beaver Dam Eye Study. Arch Ophthalmol. 2007; 125:74-9. [PMID: 17210855]
  4. Duggal P, Klein AP, Lee KE, Iyengar SK, Klein R, Bailey-Wilson JE, Klein BE. A genetic contribution to intraocular pressure: the beaver dam eye study. Invest Ophthalmol Vis Sci. 2005; 46:555-60. [PMID: 15671282]
  5. Charlesworth JC, Dyer TD, Stankovich JM, Blangero J, Mackey DA, Craig JE, Green CM, Foote SJ, Baird PN, Sale MM. Linkage to 10q22 for maximum intraocular pressure and 1p32 for maximum cup-to-disc ratio in an extended primary open-angleglaucoma pedigree. Invest Ophthalmol Vis Sci. 2005; 46:3723-9. [PMID: 16186355]
  6. Suriyapperuma SP, Child A, Desai T, Brice G, Kerr A, Crick RP, Sarfarazi M. A new locus (GLC1H) for adult-onset primary open-angle glaucoma maps to the 2p15-p16 region. Arch Ophthalmol. 2007; 125:86-92. [PMID: 17210857]
  7. Wirtz MK, Samples JR, Rust K, Lie J, Nordling L, Schilling K, Acott TS, Kramer PL. GLC1F, a new primary open-angle glaucoma locus, maps to 7q35-q36. Arch Ophthalmol. 1999; 117:237-41. [PMID: 10037570]
  8. Lu Y, Vitart V, Burdon KP, Khor CC, Bykhovskaya Y, Mirshahi A, Hewitt AW, Koehn D, Hysi PG, Ramdas WD, Zeller T, Vithana EN, Cornes BK, Tay WT, Tai ES, Cheng CY, Liu J, Foo JN, Saw SM, Thorleifsson G, Stefansson K, Dimasi DP, Mills RA, Mountain J, Ang W, Hoehn R, Verhoeven VJ, Grus F, Wolfs R, Castagne R, Lackner KJ, Springelkamp H, Yang J, Jonasson F, Leung DY, Chen LJ, Tham CC, Rudan I, Vatavuk Z, Hayward C, Gibson J, Cree AJ, MacLeod A, Ennis S, Polasek O, Campbell H, Wilson JF, Viswanathan AC, Fleck B, Li X, Siscovick D, Taylor KD, Rotter JI, Yazar S, Ulmer M, Li J, Yaspan BL, Ozel AB, Richards JE, Moroi SE, Haines JL, Kang JH, Pasquale LR, Allingham RR, Ashley-Koch A, NEIGHBOR Consortium. Mitchell P, Wang JJ, Wright AF, Pennell C, Spector TD, Young TL, Klaver CC, Martin NG, Montgomery GW, Anderson MG, Aung T, Willoughby CE, Wiggs JL, Pang CP, Thorsteinsdottir U, Lotery AJ, Hammond CJ, van Duijn CM, Hauser MA, Rabinowitz YS, Pfeiffer N, Mackey DA, Craig JE, Macgregor S, Wong TY. Genome-wide association analyses identify multiple loci associated with central corneal thickness and keratoconus. Nat Genet. 2013; 45:155-63. [PMID: 23291589]
  9. Axenovich T, Zorkoltseva I, Belonogova N, van Koolwijk LM, Borodin P, Kirichenko A, Babenko V, Ramdas WD, Amin N, Despriet DD, Vingerling JR, Lemij HG, Oostra BA, Klaver CC, Aulchenko Y, van Duijn CM. Linkage and association analyses of glaucoma related traits in a large pedigree from a Dutch genetically isolated population. J Med Genet. 2011; 48:802-9. [PMID: 22058429]
  10. Sherwin JC, Kearns LS, Hewitt AW, Ma Y, Kelly J, Griffiths LR, Mackey DA. Prevalence of chronic ocular diseases in a genetic isolate: the Norfolk Island Eye Study (NIES). Ophthalmic Epidemiol. 2011; 18:61-71. [PMID: 21401413]
  11. Mackey DA, Sherwin JC, Kearns LS, Ma Y, Kelly J, Chu BS, Macmillan R, Barbour JM, Wilkinson CH, Matovinovic E, Cox HC, Bellis C, Lea RA, Quinlan S, Griffiths LR, Hewitt AW. The Norfolk Island Eye Study (NIES): rationale, methodology and distribution of ocular biometry (biometry of the bounty). Twin Res Hum Genet. 2011; 14:42-52. [PMID: 21314255]
  12. Macgregor S, Bellis C, Lea RA, Cox H, Dyer T, Blangero J, Visscher PM, Griffiths LR. Legacy of mutiny on the Bounty: founder effect and admixture on Norfolk Island. Eur J Hum Genet. 2010; 18:67-72. [PMID: 19584896]
  13. Bellis C, Cox HC, Dyer TD, Charlesworth JC, Begley KN, Quinlan S, Lea RA, Heath SC, Blangero J, Griffiths LR. Linkage mapping of CVD risk traits in the isolated Norfolk Island population. Hum Genet. 2008; 124:543-52. [PMID: 18975005]
  14. McEvoy BP, Zhao ZZ, Macgregor S, Bellis C, Lea RA, Cox H, Montgomery GW, Griffiths LR, Visscher PM. European and Polynesian admixture in the Norfolk Island population. Heredity (Edinb). 2010; 105:229-34. [PMID: 19997123]
  15. Scharpf RB, Irizarry RA, Ritchie ME, Carvalho B, Ruczinski I. Using the R Package crlmm for Genotyping and Copy Number Estimation. J Stat Softw. 2011; 40:1-32. [PMID: 22523482]
  16. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007; 81:559-75. [PMID: 17701901]
  17. Aulchenko YS, Ripke S, Isaacs A, van Duijn CM. GenABEL: an R library for genome-wide association analysis. Bioinformatics. 2007; 23:1294-6. [PMID: 17384015]
  18. Almasy L, Blangero J. Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet. 1998; 62:1198-211. [PMID: 9545414]
  19. Dyer TD, Blangero J, Williams JT, Göring HH, Mahaney MC. The effect of pedigree complexity on quantitative trait linkage analysis. Genet Epidemiol. 2001; 21Suppl 1:S236-43. [PMID: 11793675]
  20. Manichaikul A, Mychaleckyj JC, Rich SS, Daly K, Sale M, Chen WM. Robust relationship inference in genome-wide association studies. Bioinformatics. 2010; 26:2867-73. [PMID: 20926424]
  21. Benton MC, Lea RA, Macartney-Coxson D, Hanna M, Eccles DA, Carless MA, Chambers GK, Bellis C, Goring HH, Curran JE, Harper JL, Gibson G, Blangero J, Griffiths LR. A Phenomic Scan of the Norfolk Island Genetic Isolate Identifies a Major Pleiotropic Effect Locus Associated with Metabolic and Renal Disorder Markers. PLoS Genet. 2015; 11:e1005593 [PMID: 26474483]
  22. Benton MC, Lea RA, Macartney-Coxson D, Carless MA, Göring HH, Bellis C, Hanna M, Eccles D, Chambers GK, Curran JE, Harper JL, Blangero J, Griffiths LR. Mapping eQTLs in the Norfolk Island genetic isolate identifies candidate genes for CVD risk traits. Am J Hum Genet. 2013; 93:1087-99. [PMID: 24314549]
  23. Cox HC, Bellis C, Lea RA, Quinlan S, Hughes R, Dyer T, Charlesworth J, Blangero J, Griffiths LR. Principal component and linkage analysis of cardiovascular risk traits in the Norfolk isolate. Hum Hered. 2009; 68:55-64. [PMID: 19339786]
  24. Bellis C, Hughes RM, Begley KN, Quinlan S, Lea RA, Heath SC, Blangero J, Griffiths LR. Phenotypical characterisation of the isolated norfolk island population focusing on epidemiological indicators of cardiovascular disease. Hum Hered. 2005; 60:211-9. [PMID: 16391489]
  25. Kiefer AK, Tung JY, Do CB, Hinds DA, Mountain JL, Francke U, Eriksson N. Genome-wide analysis points to roles for extracellular matrix remodeling, the visual cycle, and neuronal development in myopia. PLoS Genet. 2013; 9:e1003299 [PMID: 23468642]
  26. Tham YC, Aung T, Fan Q, Saw SM, Siantar RG, Wong TY, Cheng CY. Joint Effects of Intraocular Pressure and Myopia on Risk of Primary Open-Angle Glaucoma: The Singapore Epidemiology of Eye Diseases Study. Sci Rep. 2016; 6:19320 [PMID: 26758554]