Molecular Vision 2013; 19:1238-1246 <http://www.molvis.org/molvis/v19/1238>
Received 06 November 2012 | Accepted 04 June 2013 | Published 06 June 2013

Interrogation of the platelet-derived growth factor receptor alpha locus and corneal astigmatism in Australians of Northern European ancestry: Results of a genome-wide association study

Seyhan Yazar,1 Aniket Mishra,2 Wei Ang,3 Lisa S. Kearns,4 Jenny A. Mountain,5 Craig Pennell,3 Grant W. Montgomery,2 Terri L. Young,6,7 Christopher J. Hammond,8 Stuart Macgregor,2 David A. Mackey,1,4 Alex W. Hewitt1,4

The first two authors contributed equally to this work

1Centre for Ophthalmology and Visual Science, University of Western Australia, Lions Eye Institute, Perth, Australia; 2Queensland Institute of Medical Research, Brisbane, Australia; 3School of Women’s and Infants’ Health, University of Western Australia, Perth, Australia; 4Centre for Eye Research Australia, University of Melbourne, Royal Victorian Eye and Ear Hospital, Melbourne, Australia; 5Telethon Institute for Child Health Research, Centre for Child Health Research, University of Western Australia, Perth, Australia; 6Center for Human Genetics, Duke University Medical Center, Durham, NC; 7Duke Eye Center, Duke University, Durham, NC; 8Department of Twin Research and Genetic Epidemiology, King’s College London, St. Thomas’ Hospital, London, UK

Correspondence to: Seyhan Yazar, Lions Eye Institute, 2 Verdun St, Nedlands, Perth, Australia, 6009; Phone: +61 8 9381 0759; FAX: +61 8 9381 0700; email: seyhanyazar@lei.org.au

Abstract

Purpose: Corneal astigmatism is a common eye disorder characterized by irregularities in corneal curvature. Recently, the rs7677751 single nucleotide polymorphism (SNP) at the platelet-derived growth factor receptor alpha (PDGFRA) locus was found to be associated with corneal astigmatism in people of Asian ancestry. In the present study, we sought to replicate this finding and identify other genetic markers of corneal astigmatism in an Australian population of Northern European ancestry.

Methods: Data from two cohorts were included in this study. The first cohort consisted of 1,013 individuals who were part of the Western Australian Pregnancy Cohort (Raine) Study: 20-year follow-up Eye Study. The second cohort comprised 1,788 individuals of 857 twin families who were recruited through the Twins Eye Study in Tasmania and the Brisbane Adolescent Twin Study. Corneal astigmatism was calculated as the absolute difference between the keratometry readings in two meridians, and genotype data were extracted from genome-wide arrays. Initially, each cohort was analyzed separately, before being combined for meta- and subsequent genome-wide pathway analysis.

Results: Following meta-analysis, SNP rs7677751 at the PDGFRA locus had a combined p=0.32. No variant was found to be statistically significantly associated with corneal astigmatism at the genome-wide level (p<5.0×10−8). The SNP with strongest association was rs1164064 (p=1.86×10−6) on chromosome 3q13. Gene-based pathway analysis identified a significant association between the Gene Ontology “segmentation” (GO:0035282) pathway, corrected p=0.009.

Conclusions: Our data suggest that the PDGFRA locus does not transfer a major risk of corneal astigmatism in people of Northern European ancestry. Better-powered studies are required to validate the novel putative findings of our study.

Introduction

The majority of light refraction occurs at the air-tear film/cornea interface, as light enters the eye. Consequently, irregularities of this surface, manifesting as corneal astigmatism, can cause dramatic reductions in visual acuity. Approximately, 40% of participants from a Singapore Chinese population were astigmatic (as defined by cylindrical autorefraction readings >0.5 diopter [D]) [1,2]. In separate studies using identical definitions, more than 50% of rural Asian Indian and Persian populations were found to have astigmatism [3,4]. Interestingly, with marginally higher astigmatic groupings (either ≥0.75D or ≥1.00D), the age-adjusted prevalence of astigmatism was reported to be just over 35% in Caucasian populations from Australia and the United States [5,6].

Despite much work, the etiology of astigmatism remains poorly understood. Nonetheless, genetic and environmental factors have been suggested to have important roles in its development. Using a classical twin study, Hammond and colleagues reported that dominant and additive genetic effects accounted for approximately 46% to 79% of the phenotypic variance in corneal astigmatism [7]. Similarly, Dirani et al. found a heritability of 50% and 60% in men and women, respectively [8], while Grjibovski and colleagues calculated an overall heritability of 63% for corneal astigmatism [9].

A limited number of genome-wide association studies (GWASs) investigating corneal parameters have been conducted. Since central corneal thickness (CCT) has been found to be one of the most heritable human traits, the best studied corneal trait to date is CCT [10]. In the first published GWAS for CCT, the zinc finger 469 locus on chromosome 16q24 was identified in Australian and UK twin cohorts [11], and subsequently confirmed in other populations [12,13]. Additional quantitative trait loci for CCT have been identified on chromosomes 1p34 collagen type VIII alpha 2 gene (COL8A2), 6q14, 7q11, 9p23, 9q34 collagen type V alpha 1 gene (COL5A1), and 15q26 [12-14].

Corneal curvature (CC) is another important biometric feature, which has been recently interrogated at the genetic level. The first GWAS for keratometry (as defined as mean CC) in populations of Chinese, Malay, and Indian ancestry identified associated variants near the FK506 binding protein 12-rapamycin associated protein 1 (FRAP1) and platelet-derived growth factor receptor alpha (PDGFRA) genes [15]. Interestingly, we have recently shown that the association with CC and the PDGFRA locus is transferrable to people of Caucasian ancestry [16]. In a separate study, Fan and colleagues found that this locus was also associated with corneal astigmatism in a Singaporean Asian population [17]. The purpose of our present study was to investigate the role of the variants near PDGFRA on corneal astigmatism in people of Northern European ancestry. We also present results from a genome-wide meta-analysis for corneal astigmatism in more than 2,700 people.

Methods

Ethical approval

This study was conducted in accordance with the Declaration of Helsinki, and informed consent was obtained from all adult participants and at least one parent of the child participants before examination. Approval for this study was obtained from the Human Research Ethics Committees of the University of Western Australia, University of Tasmania, Royal Victorian Eye and Ear Hospital, and Queensland Institute of Medical Research.

Sample populations

A total of 1,013 (51.3% male) unrelated individuals from the REHS and 1,788 (56.7% female) individuals of 857 twin families who were recruited through the TEST and the BATS were included in the analysis. Demographic and phenotypic characteristics of these cohorts are shown in Table 1. Two Australian cohorts of Northern European ancestry were included in this study. In both studies, corneal astigmatism was calculated as the absolute difference between horizontal and vertical keratometry readings. An inverse normal transformation was applied to the average corneal astigmatism of both eyes and used for analysis. Participants who had a pterygium or had previously undergone ocular surgery were excluded from analysis.

Western Australian Pregnancy Cohort (Raine) Study: 20-year follow-up Eye Study

The first cohort comprised of participants who are enrolled in the Raine Study [18]. At the 20-year follow up, these individuals were invited to participate in the Raine Eye Health Study (REHS) and undertake a comprehensive eye examination [19]. As part of the examination, corneal curvature was measured in horizontal and vertical meridians with IOLMaster V:5 (Carl Zeiss Meditec AG, Jena, Germany). Three consecutive measurements of corneal curvature within 0.3D within each meridian were recorded with careful alignment and focus. DNA samples and consent for GWASs were available from the previous assessments. Genotype data were generated using the genome-wide Illumina 660 Quad Array at the Centre for Applied Genomics (Toronto, Ontario, Canada). As part of quality control, we investigated any individuals who were related with π >0.1875 (second- or third-degree relatives) and excluded individuals with a higher proportion of missing data. We also excluded people who had a high degree of missing genotyping data (>3%). The data were filtered for a Hardy–Weinberg equilibrium p value >5.7×10−7, single nucleotide polymorphism (SNP) call rate >95%, and a minor allele frequency >0.01. We conducted principal component analysis (PCA) and constructed the first five principal components for a subset of 42,888 SNPs that were not in linkage disequilibrium (LD) with each other using the EIGENSTRAT program [20]. We also performed the GWAS imputation of 22 autosomes in the MACH v1.0.16 software using the CEU samples from HapMap phase2 build 36 release 22. A linear regression model in R with a PLINK [21] interface was used to determine associations between SNPs and corneal astigmatism. The model was adjusted for age, sex, and the first two principal components that accounted for the population stratification.

Twins Eye Study in Tasmania and Brisbane Adolescent Twin Study

The second cohort comprised participants from the Twins Eye Study in Tasmania (TEST) and the Brisbane Adolescent Twin Study (BATS) [22]. In both studies, corneal curvature was measured using a Humphrey-598 Automatic Refractor/Keratometer (Carl Zeiss Meditec, Inc., Miami, FL), and there was no significant difference between the measurements of the right and left eyes (Student t test, p value=0.24). In the BATS and the TEST, DNA was obtained from either saliva or peripheral blood samples. Blood was collected in tubes containing ethylenediaminetetraacetic acid and saliva samples were collected using an Oragene saliva kit (DNA Genotek, Inc., Kanata, ON, Canada). The extracted DNA from these samples was genotyped on the Illumina HumanHap 610W Quad Arrays (Illumina Inc., San Diego, CA). The majority of the BATS samples were genotyped at deCODE Genetics (Sturlugata 8; Reykjavik, Iceland) as part of a larger project. All TEST samples and a small proportion of the BATS samples (50) were genotyped at the Centre for Inherited Disease Research (CIDR; Baltimore, MD). As outlined previously, genotype data were excluded if they did not satisfy a Hardy–Weinberg equilibrium test p value ≥10−6, SNP call rate >95%, Illumina BeadStudio GenCall score ≥0.7, or a minor allele frequency ≥1% [23].

Ancestral outliers were corrected with PCA using the smartpca program from v3.0 of EIGENSOFT [20]. The Australian twin data were compared with all populations in HapMap phase 3 and a collection of five other GenomEUTWIN populations [24,25]. When the outliers were identified and filtered, only PC1 (the difference between the African population and others) and PC2 (the difference between the East Asian population and others) with the highest eigenvalues were used. We calculated the mean and standard deviation of the ancestral relation of the collective European population for reference PC1 and PC2 scores. Any individual who fell away from the mean by >6 times the standard deviation on PC1 and PC2 were removed. Considering the sensitivity of imputation toward missingness and SNP distribution, we conducted imputation using 469,117 common SNPs from the genotyping data present in HapMap CEU I+II data (release 22, build 36). This imputation was performed using the MACH v1.0.16b and mimimac packages [26,27], which generated association statistics for 2,543,887. These SNPs further underwent quality control with the following criteria: Hardy–Weinberg equilibrium test p value ≥10−6, a minor allele frequency ≥1%, and Rsq score >0.3. A total of 2,428,106 SNPs passed the filtering step and were used for further analysis. The association of these SNPs with corneal astigmatism was tested using the –fastAssoc option in MERLIN [28]. The association model was adjusted for age and sex.

In Table 2, the quality control details of the genotyping in both studies are outlined. The PCA of both population structures is shown in Appendix 1.

Joint cohort analysis

Meta-analysis of the data from two cohorts was conducted using the β-coefficients method of the METAL program [29]. Only the common SNPs imputed in both cohorts (n about 2.5 million) were included in the meta-analysis. Regional associations were generated using SNAP [30].

Pathway analysis was undertaken using Pathway-VEGAS, an extension of the recently developed gene-based analysis tool Versatile Gene-based Association Study (VEGAS) program [31]. We selected pathways from the Gene Ontology (GO) database if the pathway size ranged in 10 to 1,000 genes, which resulted in 4,628 for further analysis.

To perform pathway analysis with Pathway-VEGAS, we first conducted a gene-based test on the summary data generated from the meta-analysis. To include most regulatory effects, each gene region was defined as being 50 kb up- and downstream of a gene. VEGAS calculated the gene-based test statistics by incorporating the effects of all SNPs in the gene region by correcting the linkage disequilibrium between the SNPs through a simulation approach for the multivariate normal distribution. Since the participants in our sample are European descendents, we used the linkage disequilibrium pattern from the HapMap2 CEU reference sample. Pathway p values were calculated by summing the χ2 test statistics of the respective gene derived from the VEGAS p values. These summarized p values were compared with 500,000 simulations where the summarized χ2 test statistics of randomly drawn genes depending on the pathway size to calculate the empirical p values of the pathway. To avoid adverse effects due to clustered genes, we considered only one gene from each cluster of genes, chosen randomly, and dropped others if the distance between them was <500 kb.

Results

No loci in the TEST/BATS or Raine populations attained genome-wide significance (p<5×10−8). Additionally, following meta-analysis on >2.5M overlapping genotyped and imputed SNPs, no locus reached the level of genome-wide significance (Figure 1 and Figure 2). Eleven loci had a nominal threshold of suggestive significance (p<1×10−5). Table 3 shows details regarding the ten most significant SNPs following the meta-analysis.

To identify genes associated with any known pathways, we tested the VEGAS results using pathways defined in the GO database. In our analysis, the top-ranking pathways were segmentation (GO:0035282) and embryonic pattern specification (GO:0009880; Table 4). Genes involved in differentiation of mesoderm (mesogenin 1 [MSGN1], mesenchyme homeobox 1 [MEOX1], mesenchyme homeobox 2 [MEOX2], teratocarcinoma-derived growth factor 1 [TDGF1]) and anterior and posterior axis formation (homeobox D8 [HOXD8], homeobox A2 [HOXA2], homeobox B6 [HOXB6]) were common in both pathways.

We found no evidence for replication of the PDGFRA locus (Figure 3). In our cohorts, the previously reported top SNP in this region (rs7677751) was not significantly associated with corneal astigmatism (beta=–0.0423, standard deviation error=0.0423; p=0.32). The minor allele frequency of rs7677751 was 0.133 and 0.123 in the Raine study and the TEST/BATS, respectively. The top SNP at this locus was rs6821576 (p=0.003).

Discussion

Understanding the molecular mechanism of cornea-related disease is useful for developing novel corrective or therapeutic strategies. Recently, through meta-analysis of five Singaporean cohorts, Fan and colleagues reported a statistically significant association between a variant (rs7677751) at the PDGFRA locus on chromosome 4q12 and corneal astigmatism [17]. In our present study, we found no strong evidence for transfer of risk for corneal astigmatism of this locus in two Australian cohorts of Northern European ancestry. Our results suggest there are underlying genetic differences between populations, which may account for differing development and prevalence of corneal astigmatism.

Although no SNP in our study was significantly associated with corneal astigmatism at the genome-wide level, we identified several putative loci, which reached a suggestive level of significance (p<1×10−5). Our strongest signal on meta-analysis (rs1151008) was located on chromosome 12p11. This SNP is approximately 100 kb centromeric to the antagonist of the mitotic exit network 1 homolog (AMN1) gene. AMN1 has been shown to be important in resetting the cell cycle [32], and the AMN1 domain of the lysine-specific demethylase 2A gene appears to inhibit keratinocyte growth in vitro [33].

Our second strongest-association (rs1164064) was on chromosome 3q13, near the developmental pluripotency-associated 4 (DPPA4) and DPPA2 genes (Figure 3). These genes have important roles in stem cell generation and are rapidly downregulated during cellular or fetal differentiation [34]. Given that the cornea is of ectodermal origin, DPPA4 regulates differentiation of embryonic stem cells into a primitive ectoderm lineage [35].

Following gene-based pathway analysis, we found that genes involved with segmentation and embryonic pattern specification were associated with the development of corneal astigmatism. In vertebrates, the periocular mesenchymal cells migrate into the cornea giving rise to cornea stroma during embryogenesis [36]. Interestingly, the MSGN1, MEOX1, MEOX2, and TDGF1 genes identified in our pathway analysis are involved in differentiation of mesoderm. Additionally, some of the genes in these pathways are part of the HOX family, which included the HOX8.1 gene that was demonstrated to be expressed during murine ocular development [37].

It is somewhat surprising that, despite our reasonable power, we were unable to replicate the association of corneal astigmatism and the rs7677751 variant [17]. We also failed to identify any locus associated with this trait at the genome-wide significance level. Our results suggest that in dissecting the genetic architecture of corneal astigmatism in people of Northern European ancestry, no major single locus will predominant, similar to other complex quantitative traits [38]. Clearly, larger, better-powered cohorts are required to intimately dissect the genetic etiology of this biometric trait.

In summary, we found no strong evidence for replication or transferability of the previously reported association between the rs7677751 variant, at the PDGFRA locus, and corneal astigmatism in our Australian cohorts of Northern European ancestry. We identified several putative loci, which clearly require replication in ongoing genetic or functional studies.

Appendix 1.

Acknowledgments

The authors thank all study participants. Funding sources include: The University of Western Australia (UWA), The Telethon Institute for Child Health Research, Raine Medical Research Foundation, Infant’s Research Foundation, Curtin University, Australian National Health and Medical Research Council Project Grants 1021105 and 350415, LEI, the Australian Foundation for the Prevention of Blindness, Alcon Research Institute, R01EY018246, Center for Inherited Diseases Research, the Clifford Craig Medical Research Trust, Ophthalmic Research Institute of Australia (ORIA), American Health Assistance Foundation, Peggy and Leslie Cranbourne Foundation, Foundation for Children, Jack Brockhoff Foundation, the Pfizer Australia Senior Research Fellowship, and the Victorian Government of Australia.

References

  1. 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]
  2. Wu HM, Seet B, Yap EP, Saw SM, Lim TH, Chia KS. Does education explain ethnic differences in myopia prevalence? A population-based study of young adult males in Singapore. Optom Vis Sci. 2001; 78:234-9. [PMID: 11349931]
  3. Raju P, Ramesh SV, Arvind H, George R, Baskaran M, Paul PG, Kumaramanickavel G, McCarty C, Vijaya L. Prevalence of refractive errors in a rural South Indian population. Invest Ophthalmol Vis Sci. 2004; 45:4268-72. [PMID: 15557431]
  4. Hashemi H, Hatef E, Fotouhi A, Mohammad K. Astigmatism and its determinants in the Tehran population: the Tehran eye study. Ophthalmic Epidemiol. 2005; 12:373-81. [PMID: 16283989]
  5. Attebo K, Ivers RQ, Mitchell P. Refractive errors in an older population: the Blue Mountains Eye Study. Ophthalmology. 1999; 106:1066-72. [PMID: 10366072]
  6. Vitale S, Ellwein L, Cotch MF, Ferris FL, Sperduto R. Prevalence of refractive error in the United States, 1999–2004. Arch Ophthalmol. 2008; 126:1111-9. [PMID: 18695106]
  7. 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]
  8. Dirani M, Islam A, Shekar SN, Baird PN. Dominant genetic effects on corneal astigmatism: the genes in myopia (GEM) twin study. Invest Ophthalmol Vis Sci. 2008; 49:1339-44. [PMID: 18385047]
  9. Grjibovski AM, Magnus P, Midelfart A, Harris JR. Epidemiology and Heritability of Astigmatism in Norwegian Twins: An Analysis of Self-Reported Data. Ophthalmic Epidemiol. 2006; 13:245-52. [PMID: 16877283]
  10. Toh T, Liew SHM, MacKinnon JR, Hewitt AW, Poulsen JL, Spector TD, Gilbert CE, Craig JE, Hammond CJ, Mackey DA. Central corneal thickness is highly heritable: the twin eye studies. Invest Ophthalmol Vis Sci. 2005; 46:3718-22. [PMID: 16186354]
  11. Lu Y, Dimasi DP, Hysi PG, Hewitt AW, Burdon KP, Toh T, Ruddle JB, Li Y-J, Mitchell P, Healey PR, Montgomery GW, Hansell N, Spector TD, Martin NG, Young TL, Hammond CJ, Macgregor S, Craig JE, Mackey DA. Common Genetic Variants near the Brittle Cornea Syndrome Locus ZNF469 Influence the Blinding Disease Risk Factor Central Corneal Thickness. PLoS Genet. 2010; 6:e1000947 [PMID: 20485516]
  12. Vitart V, Bencić G, Hayward C, Skunca Herman J, Huffman J, Campbell S, Bućan K, Navarro P, Gunjaca G, Marin J, Zgaga L, Kolcić I, Polasek O, Kirin M, Hastie ND, Wilson JF, Rudan I, Campbell H, Vatavuk Z, Fleck B, Wright A. New loci associated with central cornea thickness include COL5A1, AKAP13 and AVGR8. Hum Mol Genet. 2010; 19:4304-11. [PMID: 20719862]
  13. Vithana EN, Aung T, Khor CC, Cornes BK, Tay W-T, Sim X, Lavanya R, Wu R, Zheng Y, Hibberd ML, Chia K-S, Seielstad M, Goh L-K, Saw S-M, Tai ES, Wong TY. Collagen-related genes influence the glaucoma risk factor, central corneal thickness. Hum Mol Genet. 2011; 20:649-58. [PMID: 21098505]
  14. Cornes BK, Khor CC, Nongpiur ME, Xu L, Tay W-T, Zheng Y, Lavanya R, Li Y, Wu R, Sim X, Wang Y-X, Chen P, Teo Y-Y, Chia K-S, Seielstad M, Liu J, Hibberd ML, Cheng C-Y, Saw S-M, Tai ES, Jonas JB, Vithana EN, Wong TY, Aung T. Identification of four novel variants that influence central corneal thickness in multi-ethnic Asian populations. Hum Mol Genet. 2012; 21:437-45. [PMID: 21984434]
  15. Han S, Chen P, Fan Q, Khor CC, Sim X, Tay W-T, Ong RT-H, Suo C, Goh L-K, Lavanya R, Zheng Y, Wu R, Seielstad M, Vithana E, Liu J, Chia K-S, Lee JJ-M, Tai ES, Wong TY, Aung T, Teo Y-Y, Saw S-M. Association of variants in FRAP1 and PDGFRA with corneal curvature in Asian populations from Singapore. Hum Mol Genet. 2011; 20:3693-8. [PMID: 21665993]
  16. Mishra A, Yazar S, Hewitt AW, Mountain JA, Ang W, Pennell CE, Martin NG, Montgomery GW, Hammond CJ, Young TL, Macgregor S, Mackey DA. Genetic variants near PDGFRA are associated with corneal curvature in Australians. Invest Ophthalmol Vis Sci. 2012; 53:7131-6. [PMID: 22969067]
  17. Fan Q, Zhou X, Khor CC, Cheng C-Y, Goh L-K, Sim X, Tay W-T, Li Y-J, Ong RT-H, Suo C, Cornes B, Ikram MK, Chia K-S, Seielstad M, Liu J, Vithana E, Young TL, Tai ES, Wong TY, Aung T, Teo Y-Y, Saw S-M. Genome-Wide Meta-Analysis of Five Asian Cohorts Identifies PDGFRA as a Susceptibility Locus for Corneal Astigmatism. Wiggs J, editor. PLoS Genet. 2011;7(12):e1002402.
  18. Macdonald W, Newnham J, Gurrin L. Effect of frequent prenatal ultrasound on birthweight: follow up at 1 year of age. Lancet. 1996; 348:482 [PMID: 8709817]
  19. Yazar S, Forward H, McKnight CM, Tan A, Soloshenko A, Oates SK, Ang W, Sherwin JC, Wood D, Mountain JA, Pennell CE, Hewitt AW, Mackey DA. Raine Eye Health Study: Design, Methodology and Baseline Prevalence of Ophthalmic Disease in a Birth-cohort Study of Young Adults. Ophthalmic Genet. 2013; [PMID: 23301674]
  20. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006; 38:904-9. [PMID: 16862161]
  21. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, de Bakker PIW, 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]
  22. Mackey DA, MacKinnon JR, Brown SA, Kearns LS, Ruddle JB, Sanfilippo PG, Sun C, Hammond CJ, Young TL, Martin NG, Hewitt AW. Twins Eye Study in Tasmania (TEST): Rationale and Methodology to Recruit and Examine Twins. Twin Res Hum Genet. 2009; 12:441-54. [PMID: 19803772]
  23. Macgregor S, Hewitt AW, Hysi PG, Ruddle JB, Medland SE, Henders AK, Gordon SD, Andrew T, McEvoy B, Sanfilippo PG, Carbonaro F, Tah V, Li YJ, Bennett SL, Craig JE, Montgomery GW, Tran-Viet KN, Brown NL, Spector TD, Martin NG, Young TL, Hammond CJ, Mackey DA. Genome-wide association identifies ATOH7 as a major gene determining human optic disc size. Hum Mol Genet. 2010; 19:2716-24. [PMID: 20395239]
  24. Peltonen L. GenomEUtwin. GenomEUtwin: a strategy to identify genetic influences on health and disease. Twin Res. 2003; 6:354-60.
    Internet
    [PMID: 14624719]
  25. McEvoy BP, Montgomery GW, McRae AF, Ripatti S, Perola M, Spector TD, Cherkas L, Ahmadi KR, Boomsma D, Willemsen G, Hottenga JJ, Pedersen NL, Magnusson PKE, Kyvik KO, Christensen K, Kaprio J, Heikkila K, Palotie A, Widen E, Muilu J, Syvanen AC, Liljedahl U, Hardiman O, Cronin S, Peltonen L, Martin NG, Visscher PM. Geographical structure and differential natural selection among North European populations. Genome Res. 2009; 19:804-14. [PMID: 19265028]
  26. 1000 Genomes Project Consortium. A map of human genome variation from population-scale sequencing. Nature. 2010; 467:1061-73. [PMID: 20981092]
  27. Li Y, Willer CJ, Ding J, Scheet P, Abecasis GR. MaCH: using sequence and genotype data to estimate haplotypes and unobserved genotypes. Genet Epidemiol. 2010; 34:816-34. [PMID: 21058334]
  28. 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]
  29. Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010; 26:2190-1. [PMID: 20616382]
  30. Johnson AD, Handsaker RE, Pulit SL, Nizzari MM, O'Donnell CJ, de Bakker PIW. SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics. 2008; 24:2938-9. [PMID: 18974171]
  31. Liu JZ, McRae AF, Nyholt DR, Medland SE, Wray NR, Brown KM, AMFS Investigators. Hayward NK, Montgomery GW, Visscher PM, Martin NG, Macgregor S. A versatile gene-based test for genome-wide association studies. Am J Hum Genet. 2010; 87:139-45. [PMID: 20598278]
  32. Wang Y, Shirogane T, Liu D, Harper JW, Elledge SJ. Exit from exit: resetting the cell cycle through Amn1 inhibition of G protein signaling. Cell. 2003; 112:697-709. [PMID: 12628189]
  33. Iuchi S, Green H. Lysine-specific demethylase 2A (KDM2A) normalizes human embryonic stem cell derived keratinocytes. Proc Natl Acad Sci USA. 2012; 109:9442-7. [PMID: 22635273]
  34. Maldonado-Saldivia J, van den Bergen J, Krouskos M, Gilchrist M, Lee C, Li R, Sinclair AH, Surani MA, Western PS. Dppa2 and Dppa4 are closely linked SAP motif genes restricted to pluripotent cells and the germ line. Stem Cells. 2007; 25:19-28. [PMID: 16990585]
  35. Masaki H, Nishida T, Kitajima S, Asahina K, Teraoka H. Developmental pluripotency-associated 4 (DPPA4) localized in active chromatin inhibits mouse embryonic stem cell differentiation into a primitive ectoderm lineage. J Biol Chem. 2007; 282:33034-42. [PMID: 17855347]
  36. Graw J. Chapter Ten - Eye Development. Current Topics in Developmental Biology - Volume 90. Elsevier Inc; 2010. pages 343–86.
  37. Monaghan AP, Davidson DR, Sime C, Graham E, Baldock R, Bhattacharya SS, Hill RE. The Msh-like homeobox genes define domains in the developing vertebrate eye. Development. 1991; 112:1053-61. [PMID: 1682128]
  38. Lettre G, Jackson AU, Gieger C, Schumacher FR, Berndt SI, Sanna S, Eyheramendy S, Voight BF, Butler JL, Guiducci C, Illig T, Hackett R, Heid IM, Jacobs KB, Lyssenko V, Uda M. Diabetes Genetics Initiative, FUSION, KORA, Prostate, Lung Colorectal and Ovarian Cancer Screening Trial, Nurses' Health Study, SardiNIA, Boehnke M, Chanock SJ, Groop LC, Hu FB, Isomaa B, Kraft P, Peltonen L, Salomaa V, Schlessinger D, Hunter DJ, Hayes RB, Abecasis GR, Wichmann H-E, Mohlke KL, Hirschhorn JN. Identification of ten loci associated with height highlights new biological pathways in human growth. Nature Publishing Group.. 2008; 40:584-91.