Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Association analyses of rare variants identify two genes associated with refractive error

  • Karina Patasova ,

    Contributed equally to this work with: Karina Patasova, Annechien E. G. Haarman, Anthony M. Musolf, Jeremy A. Guggenheim, Chris J. Hammond, Pirro G. Hysi

    Roles Formal analysis, Writing – original draft

    Affiliations Department of Ophthalmology, King’s College London, London, United Kingdom, Department of Twins Research and Genetic Epidemiology, King’s College London, London, United Kingdom

  • Annechien E. G. Haarman ,

    Contributed equally to this work with: Karina Patasova, Annechien E. G. Haarman, Anthony M. Musolf, Jeremy A. Guggenheim, Chris J. Hammond, Pirro G. Hysi

    Roles Validation, Writing – review & editing

    Affiliations Department of Ophthalmology, Erasmus Medical Center, Rotterdam, The Netherlands, Department of Epidemiology, Erasmus MC, Rotterdam, The Netherlands

  • Anthony M. Musolf ,

    Contributed equally to this work with: Karina Patasova, Annechien E. G. Haarman, Anthony M. Musolf, Jeremy A. Guggenheim, Chris J. Hammond, Pirro G. Hysi

    Roles Validation, Writing – review & editing

    Affiliation Computational and Statistical Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Baltimore, Maryland, United States of America

  • Omar A. Mahroo,

    Roles Writing – review & editing

    Affiliations Department of Ophthalmology, King’s College London, London, United Kingdom, Department of Twins Research and Genetic Epidemiology, King’s College London, London, United Kingdom, NIHR Biomedical Research Centre at Moorfields Eye Hospital NHS Foundation Trust and the UCL Institute of Ophthalmology, London, United Kingdom, Department of Ophthalmology, St Thomas’ Hospital, Guys and St ’Thomas’ NHS Foundation Trust, London, United Kingdom, Institute of Ophthalmology, University College London, London, United Kingdom

  • Jugnoo S. Rahi,

    Roles Conceptualization, Investigation, Writing – review & editing

    Affiliations UCL Great Ormond Street Hospital Institute of Child Health, London, United Kingdom, Ulverscroft Vision Research Group, University College London, London, United Kingdom

  • Mario Falchi,

    Roles Investigation, Methodology, Supervision

    Affiliation Department of Twins Research and Genetic Epidemiology, King’s College London, London, United Kingdom

  • Virginie J. M. Verhoeven,

    Roles Writing – review & editing

    Affiliations Department of Ophthalmology, Erasmus Medical Center, Rotterdam, The Netherlands, Department of Epidemiology, Erasmus MC, Rotterdam, The Netherlands

  • Joan E. Bailey-Wilson,

    Roles Investigation, Methodology, Writing – review & editing

    Affiliation Computational and Statistical Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Baltimore, Maryland, United States of America

  • Caroline C. W. Klaver,

    Roles Supervision, Writing – review & editing

    Affiliations Department of Ophthalmology, Erasmus Medical Center, Rotterdam, The Netherlands, Department of Epidemiology, Erasmus MC, Rotterdam, The Netherlands, Department of Ophthalmology, Radboud University Medical Center, Nijmegen, the Netherlands

  • Priya Duggal,

    Roles Writing – review & editing

    Affiliation Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, United States of America

  • Alison Klein,

    Roles Resources, Writing – review & editing

    Affiliations Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, United States of America, Department of Oncology, Sidney Kimmel Comprehensive Cancer Center at Johns Hopkins, Johns Hopkins University, Baltimore, Maryland, United States of America, Department of Pathology, Johns Hopkins School of Medicine, Johns Hopkins University, Baltimore, Maryland, United States of America

  • Jeremy A. Guggenheim ,

    Contributed equally to this work with: Karina Patasova, Annechien E. G. Haarman, Anthony M. Musolf, Jeremy A. Guggenheim, Chris J. Hammond, Pirro G. Hysi

    Roles Investigation, Supervision, Validation, Writing – review & editing

    Affiliation School of Optometry & Vision Sciences, Cardiff University, Cardiff, United Kingdom

  • Chris J. Hammond ,

    Contributed equally to this work with: Karina Patasova, Annechien E. G. Haarman, Anthony M. Musolf, Jeremy A. Guggenheim, Chris J. Hammond, Pirro G. Hysi

    Roles Conceptualization, Investigation, Methodology, Resources, Writing – review & editing

    Affiliations Department of Ophthalmology, King’s College London, London, United Kingdom, Department of Twins Research and Genetic Epidemiology, King’s College London, London, United Kingdom

  • Pirro G. Hysi ,

    Contributed equally to this work with: Karina Patasova, Annechien E. G. Haarman, Anthony M. Musolf, Jeremy A. Guggenheim, Chris J. Hammond, Pirro G. Hysi

    Roles Conceptualization, Investigation, Methodology, Resources, Writing – review & editing

    pirro.hysi@kcl.ac.uk

    Affiliations Department of Ophthalmology, King’s College London, London, United Kingdom, Department of Twins Research and Genetic Epidemiology, King’s College London, London, United Kingdom, UCL Great Ormond Street Hospital Institute of Child Health, London, United Kingdom

  • the CREAM Consortium; the UK Biobank Eye; Vision Consortium

    Membership of the CREAM Consortium and the UK Biobank Eye and Vision Consortium is provided in the Acknowledgments.

Abstract

Purpose

Genetic variants identified through population-based genome-wide studies are generally of high frequency, exerting their action in the central part of the refractive error spectrum. However, the power to identify associations with variants of lower minor allele frequency is greatly reduced, requiring considerable sample sizes. Here we aim to assess the impact of rare variants on genetic variation of refractive errors in a very large general population cohort.

Methods

Genetic association analyses of non-cyclopaedic autorefraction calculated as mean spherical equivalent (SPHE) used whole-exome sequence genotypic information from 50,893 unrelated participants in the UK Biobank of European ancestry. Gene-based analyses tested for association with SPHE using an optimised SNP-set kernel association test (SKAT-O) restricted to rare variants (minor allele frequency < 1%) within protein-coding regions of the genome. All models were adjusted for age, sex and common lead variants within the same locus reported by previous genome-wide association studies. Potentially causal markers driving association at significant loci were elucidated using sensitivity analyses by sequentially dropping the most associated variants from gene-based analyses.

Results

We found strong statistical evidence for association of SPHE with the SIX6 (p-value = 2.15 x 10−10, or Bonferroni-Corrected p = 4.41x10-06) and the CRX gene (p-value = 6.65 x 10−08, or Bonferroni-Corrected p = 0.001). The SIX6 gene codes for a transcription factor believed to be critical to the eye, retina and optic disc development and morphology, while CRX regulates photoreceptor specification and expression of over 700 genes in the retina. These novel associations suggest an important role of genes involved in eye morphogenesis in refractive error.

Conclusion

The results of our study support previous research highlighting the importance of rare variants to the genetic risk of refractive error. We explain some of the origins of the genetic signals seen in GWAS but also report for the first time a completely novel association with the CRX gene.

Introduction

Myopia is a common eye disorder characterised by an imbalance between different refractive components of the eye and axial length [1]. The prevalence rates of myopia and its related complications are on the rise in South East Asia and rapidly increasing in Europe and the US [2]. Both environmental and genetic factors play a role in the pathophysiology of refractive errors (RE). Refractive errors, especially high and pathological myopia, are leading causes of preventable vision loss worldwide and sources of significant ocular complications [1]. Previous studies have shown that individuals with myopia greater than 6 dioptres are at increased risk of other eye conditions and are more susceptible to several sight-threatening complications [1].

Recent genome-wide studies discovered hundreds of distinct loci harbouring genes involved in refraction development [3]. Cumulatively those variants accounted for approximately 18% of total heritability [3], while twin studies estimate the RE heritability between 50%-90% [4]. The missing heritability in refractive error GWAS could be attributed to several causes, such as confounding arising from linkage disequilibrium, statistical power limitations, epistasis and heritability explained by rare genetic variants that are usually not identifiable by traditional genetic association studies. Statistical power to detect associations with individual genetic variants is proportional to the magnitude of the risk they individually confer, but more crucially, to their frequency in the population.

Here we aim to evaluate the impact of lower frequency genetic variation on refractive error in the UK Biobank sub-sample of 50,893 unrelated individuals of European descent.

Materials and methods

Study population

The UK Biobank is a cohort of 500,000 volunteers for whom extensive demographic, phenotypic and biomarker data is available. Approximately 23% of individuals (N = 117,279) participated in comprehensive eye examination, including assessment of non-cycloplegic refractive error represented by mean spherical equivalent [5]. The measurements of non-cycloplegic autorefraction were performed by Tomey RC 5000 device (Tomey Corp., Nagoya, Japan) for each eye separately, and the mean spherical equivalent of the two eyes was taken.

The study was conducted with the approval of the North-West Research Ethics Committee (ref. 06/MRE08/65), following the principles of the Declaration of Helsinki. Study participants provided written informed consent.

Full genotypic data was available for all participants who were included in our analyses. Based on a pre-computed Principal Component Analysis, we ascertained individuals’ ancestry identity by state. Only individuals who were of full and homogenenous European ancestry were analysed; in cases where familial relatedness was observed (defined as a PI_HAT > 0.06), only one individual for each pair was included in the analyses.

Whole exome sequencing data

Approximately 200,629 participants from the UK Biobank cohort were selected for the second tranche of whole-exome sequencing. Individuals who had more complete data, such as baseline measurements, MRI imaging, hospital episodes and primary care records, were prioritised for the sequencing [6]. Although the whole-exome study sample was enriched for clinical outcomes and availability of physical measures, it remained largely representative of the general UK Biobank cohort in terms of demographic characteristics and composition. Exome sequencing was performed as described previously [6]. The panel targeted 39 Mb of the human genome and covered 19,396 genes on autosomal and sex chromosomes, including 4,735,722 variants within the targeted regions, comprising 1,229,303 coding synonymous, 2,498,947 nonsynonymous and 231,631 potential loss of function variants within at least one coding transcript. In addition to targeted regions, some variation in exome adjacent regions was also captured–precisely 9,693,526 nucleotide and indel variants. About 98% of the sequenced coding variants had an allele frequency below 1%. To avert issues related to population structure, we restricted our analyses to a homogeneous sample of European ancestry, as ascertained through a principal component analysis of the directly genotyped variants in the sample, as described before [3].

The CREAM consortium dataset

CREAM (Consortium for Refractive Error and Myopia) was established in 2011 as a collaboration between studies with data on refractive error which had performed genome-wide association analysis based on SNP arrays. For the current study, we included 10 participating studies with available exome chip data. These studies included: Singapore Indian Eye Study (SINDI), Age-Related Eye Study (AREDS), Rotterdam Study I (RSI), Erasmus Rucphen Family Study (ERF), Estonian Genome Center of the University of Tartu (EGCUT), Finnish Twin Study on Aging (FITSA), Ogliastra, Croatia-Korcula, TwinsUK, Raine eye health study (REHS) and Beaver Dam eye study (BDES).

The phenotypes for all individuals participating in any of the CREAM cohorts were assessed through methodologies similar to those used for the UK Bobank participants [7]. Specifically, refractive error in the CREAM-participating cohorts was measured using autorefraction. All cohorts had been genotyped on either the Illumina HumanExome-12 v 1.0 or v 1.1 array. All cohorts were jointly recalled to obtain a larger sample size of rare variants (here defined as variants with a minor allele frequency (MAF) < 0.01), as recalling genotypes simultaneously across all samples increases the ability to call rare variants with a more discrete distinction between allele calls and sensitivity for low-frequency (high-intensity) loci. All data was recalled using GenomeStudio® v2011.1 (Illumina Inc., San Diego, CA, USA). Nine of these predominantly European CREAM cohorts were combined in a single cohort (N = 11,505), henceforth referred to as the CREAM cohort, for analysis. Because of legal requirements the BDES data (N = 1740) made it impossible to analyse this cohort alongside the rest of cohorts participating in CREAM.

Statistical analyses

To minimise confounding arising from population structure, the study sample was restricted to 50,893 unrelated UK Biobank participants of European descent. The ancestry and relatedness information was calculated based on the genetic data made available from the UK Biobank [8]. Individuals with European descent were identified by projecting UK Biobank participants onto the coordination of 1000 Genome Project principal components. The genetic data was used to identify related individuals by calculating kinship coefficients for all pairs and third-degree or closer relatives were excluded.

Gene-based analyses were conducted in the optimised SNP-set kernel association test SKAT-O test [9] implemented in the rvtests package [10]. The spherical equivalent measurements were the dependent variable and the weighed allelic burden the independent variables. All analyses were adjusted for age and sex. Our analyses incorporated several variant annotations that previous works have shown to boost the power and accuracy of detecting causal associations in gene-based analyses [11]. Variants in protein-coding regions of genes including synonymous and non-synonymous, stop gain/loss, start gain/loss or splice-site mutations with minor allele frequency below 1% were selected for inclusion. The splicing sites were defined as 3 bases into exon and 8 bases into an intron. Mutations in these regions were annotated as "Normal_Splice_Site" unless they affected the functionally important "GU…AG" region of the intron which was annotated as "Essential_Splice_Site". We excluded UTR variants and polymorphisms with unknown or inconclusive molecular consequences such as intronic variants. We used the GRCh38/hg38 assembly of the human genome (https://www.ncbi.nlm.nih.gov/assembly/GCA_000001405.29) as a reference and variants were identified and annotated using the ANNO package (https://github.com/zhanxw/anno). Gene-based associations with probabilities below the selected Bonferroni multiple testing correction level derived by dividing 0.05 by the total number of 19,293 protein-coding genes that analysed were considered statistically significant (p < 2.59 × 10−6). We sought replication of significant genes using the results of the gene-based analysis performed in the predominantly European CREAM cohort, described elsewhere [12] and the BDES cohort described elsewhere [13, 14]. Replication was considered successful if the association probabilities were below the selected Bonferroni multiple testing correction level.

Sequential analyses evaluating the role of single variants in gene-based associations.

To elucidate which variants were driving observed associations with candidate genes, we performed sequential sensitivity analyses by progressively removing markers from the gene-based analyses. The associations with target genes were assessed using the SKAT-O test adjusted for age, sex and lead common variants within the same locus. The lead common variants were selected from previously published refractive error GWAS [3]. Minor allele frequencies observed were compared with those reported in the gnomAD database [15] and pairwise linkage disequilibrium between any two variants was calculated with reference to the entire European panel included in the “ldlink” online tool (https://ldlink.nci.nih.gov/).

Results

The study sample included 50893 unrelated UK Biobank participants of European descent; 54% were women with a median age of 57 years (±8 years). Detailed information about the study participants’ demographic and refractive characteristics can be found in Table 1 and characteristics of the spherical equivalent (SPHE) distribution in S1 Fig.

We ran SKAT-O tests using 2,923,839 rare (minor allele frequency, MAF < 0.01) variants in 19,293 genes (Fig 1). The statistically strongest association was observed between SPHE and SIX6 gene (p-value = 2.15 x10-10). The second Bonferroni-significant association was found with CRX (p = 6.65 x 10−08). This finding was novel and not described in prior refractive error GWAS. Suggestive statistical evidence of association was found for the RPSAP52 (p = 1.65 x 10−05), PCCA (p = 1.82 x 10−05), MIR4683 (p = 2.81 x 10−05), SELENOM (p = 3.52 x 10−05), NAPA (p = 4.55 x 10−05) and VWA8 (p = 5.68x10-05) genes, whose association however did not meet our criteria of statistical significance after multiple testing correction (Table 2).

thumbnail
Fig 1. Manhattan plot displaying SKAT-O association results.

Each point represents one of the 19,293 genes tested for the association with the spherical equivalent in the UK Biobank cohort (N = 50,893). The plot shows -log10 transformed p-values for each gene plotted against the chromosomal location. The red dashed line indicates the Bonferroni significance threshold (p < 2.59 × 10−6). Regions are named with symbols of genes that were most strongly associated with refractive error.

https://doi.org/10.1371/journal.pone.0272379.g001

thumbnail
Table 2. Top eight gene associations with refractive error.

https://doi.org/10.1371/journal.pone.0272379.t002

No other exome sequencing datasets of comparable size were available to us. However, two smaller cohorts were genotyped for a selection of rare exonic variants using SNP chip arrays, in both of which data from only four of our candidate genes, including SIX6, CRX and two genes associated with SPHE at suggestive levels in the UK Biobank analysis were available. In addition, only 2 and 3 variants were present in the exome chip data for the SIX6 and CRX genes respectively, that had demonstrated the statistically strongest relationship with spherical equivalent in the discovery cohort. Unsurprisingly given the low minor allele frequencies, none of these variants was in strong linkage disequilibrium with the rare variants that showed significant association with the phenotype in the discovery cohort (R2 < 0.002 for all of them), although they were likely located on the same haplotypes (D’ = 1). None of these genes was associated at statistically significant levels with SHPE in the pooled exome chip cohort from the Consortium for the Refractive Error and Myopia (CREAM, N = 11,505), but we observed a strong association for the NAPA gene (SKAT-O p = 3.73 x 10−05), in the smaller Beaver Dam Eye Study (BDES, N = 1,740) cohort (Table 3).

thumbnail
Table 3. Replication of four loci associated with refractive error using gene-based analyses performed in Beaver Dam (N = 1740) and CREAM Consortium dataset (N = 11,505).

https://doi.org/10.1371/journal.pone.0272379.t003

To identify independent variants driving gene-based associations at the SIX6 and CRX loci, sensitivity analyses were performed by progressively removing SNPs from gene-based analyses. The removal of rare variants from the gene-based SKAT-O analyses revealed a decrease in the statistical significance of the analyses. The results of these analyses suggested that association with the SIX6 gene was most strongly influenced by the rs146737847 variant, whose removal resulted in the loss of statistical significance in our samples (S2 Fig). Similarly, exonic marker rs61748438 was identified as a lead variant in the CRX locus (S3 Fig). The removal of other functionally important variants within this gene also resulted in a progressive decrease in statistical significance. This gradual decrease may suggest that although gene-based association at both loci is mostly due to the presence of a few lead variants, additional lower frequency variants within these genes may also contribute to associations with spherical equivalent, but our abilities to fully evaluate their role at a general population level may be constrained by sample size and statistical power limitations.

Discussion

Here we report significant associations between spherical equivalent and rare variants located within SIX6 and CRX, but also RPSAP52, PCCA, MIR4683, SELENOM, NAPA and VWA8 genes. In our study, the strongest association was observed with the SIX6 gene, located on 14q23.1 and which encodes a homeobox protein involved in ocular development [16], morphogenesis and visual perception [17]. The SIX Homeobox 6 (SIX6) is part of a group of evolutionarily conserved genes, which were known eye transcription factors [18], which regulate the proliferation of specific retinal cells during optic disc development [18] and retain their importance in the mature retina [18]. SIX6 is implicated both in the early stages of eye formation and subsequent differentiation of retinal progenitor cells (RPC). Interestingly, previous works have shown that the rs146737847 (Glu129Lys) adversely affects the SIX6 gene function [19] and is also associated with primary open-angle glaucoma potentially through its known effect over the vertical cup-disc ratio [20]. While observational correlation between glaucoma and myopia status is well known [21, 22] there is little evidence of large-scale shared allelic risk between spherical equivalent and vertical cup-to-disc ratio [3]. The associations observed with both spherical equivalent and glaucoma phenotypes for the rs146737847 suggest that genetic pleiotropy may explain a considerable proportion of the phenotypic correlation between these two common ocular conditions.

The Cone-Rod homeobox gene, or CRX, located on 19q13.33, encodes a photoreceptor-specific transcription factor [23]. Although a previous association with refractive error was detected for the broader chromosomal location (1), this is the first time that direct evidence links this gene with spherical equivalent or myopia. This gene is a master regulator of photoreceptor development [24] and differentiation [25]. Certain mutations in this gene cause several retinal disorders, including cone-rod dystrophy, retinitis pigmentosa, adult-onset macular dystrophy and Leber congenital amaurosis [23, 26]. Additionally, the Cone-Rod Homeobox (CRX) Transcription Factor regulates the expression of over 700 genes in the retina, including downstream effects over rhodopsin and cone arrestin [27]. CRX expression in the retinal cells was inhibited by light stimulation, a mechanism previously implicated in myopia development [28, 29].

We identified suggestive associations with rare variants located within other genes and SPHE. In particular, our analyses implicate the NAPA and PCCA genes. Common polymorphisms at genomic loci encompassing these genes are associated with refractive error [3] and the age of refractive correction [7], but this is the first time that, to our knowledge, rare variants within their coding regions are associated with SPHE. PCAA encoded the biotin-binding region of mitochondrial Propionyl-CoA carboxylase involved branched and odd chain fatty-acid and cholesterol catabolism [30]. The protein product of the NAPA gene is a member of the soluble NSF attachment protein family aiding the fusion and docking of vesicles to target membranes [31]. NAPA also participates in synaptic activity and plays a role in neurogenesis [31]. Notably, this particular gene was the only candidate that achieved replication in an independent dataset.

Several novel potential candidates for which we found suggestive evidence of association are implicated in cognitive development and learning difficulty disorders. In particular, Ribosomal Protein SA Pseudogene 52 (RPSAP52) is associated with brain structure variations in TWAS [32] and described in genetic investigations of cognitive impairment, neurodevelopmental and neurodegenerative disorders [33]. The polymorphisms within the RPSAP52 gene were associated with schizophrenia in founder populations [34] and associated with biomarkers of Alzheimer’s disease, such as cerebrospinal fluid beta-amyloid 1–42 levels [35]. Similarly, mutations in microRNA MIR4683 may be associated with epilepsy in children [36]. Another interesting candidate VWA8 encodes a mitochondrial ATPase, whose precise function is not fully understood [37]. Genome-wide association studies demonstrated that variation in VWA8 may influence susceptibility to autism [38] and bipolar disorder [39], and also educational attainment and mathematical ability [40]. SELENOM, another novel candidate suggested by our analyses, encoded a selenoprotein that is highly expressed in the brain and that is thought to be essential for normal neurocognitive development [41].

For this study, we used some of the largest ever sample sizes analyzed to date to assess the role of rare variants in refractive error. However, our study has a number of limitations. In our analyses, we assumed a simple dominant model of inheritance, while recessive or compound heterozygosity models of inheritance may also play a role in refractive disorders. Our analyses were restricted to the coding regions of the genome. However, non-coding areas of the genome also proved to be important for several other diseases [4244] and could potentially provide a new direction for additional myopia work. In our study, we used chip array information in two replication datasets and found relatively little evidence for replication. However, the arrays only include a small number of variants within the exome of the genes of interest and none of which was a particularly strong contributor to the overall gene-based association in the discovery data. This is a general limitation of array-based studies, and to fully validate our results, future work on large scale exome sequencing on independent cohorts will be needed. Additionally, despite the large sample size, statistical power for rare variants is often limited due to the very low allele frequencies. Power will benefit from additional sequencing data from the several national cohorts and biobanks whose data will become available in the future. Finally, the results obtained from an exclusively European population sample are unfortunately not representative of more diverse populations and may not be generalisable to other ancestral groups.

This study demonstrates that variants with significantly large effects on refractive error are extremely rare (Table 2). We identify associations between population spherical equivalent and rare variants located within the protein-coding regions of the SIX6 gene, which plays an important role in eye morphogenesis and is implicated in several ocular disorders, including myopia. We also identify the CRX gene, a transcription factor crucial for the development of photoreceptors, as the origin of an important association signal. Our investigation demonstrates high-quality whole-exome sequencing provides a superior alternative to array-based methods that have power limitations and are prone to bias arising from population admixture [45]. Beyond novel associations, the incorporation of rare variants in existing myopia risk prediction models that currently rely on common polymorphisms will improve their accuracy and augment our understanding of refractive disorders.

Supporting information

S1 Fig. Spherical equivalent distribution in UK Biobank cohort (N = 50,893).

The distribution of the spherical equivalent (x-axis) in the samples; the number of participants for each spherical equivalent bin is given in the y-axis.

https://doi.org/10.1371/journal.pone.0272379.s001

(PNG)

S2 Fig. Sensitivity analyses results for the SIX6 gene.

Y-axis shows the number of SIX6 variants included in gene-based analyses, testing associations with SPHE. The model was adjusted for age, sex and the best common variant within the same locus. The -log(p-values) from SKAT-O tests are displayed on X-axis.

https://doi.org/10.1371/journal.pone.0272379.s002

(PNG)

S3 Fig. Sensitivity analyses results for the CRX gene.

Y-axis shows the number of CRX variants included in gene-based analyses, testing associations with SPHE. The model was adjusted for age, sex and the best common variant within the same locus. The -log(p-values) from SKAT-O tests are displayed on X-axis.

https://doi.org/10.1371/journal.pone.0272379.s003

(PNG)

Acknowledgments

We would like to thank UK Biobank and CREAM consortium participants for their contribution to this research project.

UK Biobank Eye and Vision Consortium Membership

  • Prof Tariq ASLAM—Manchester University, Manchester, United Kingdom
  • Prof Sarah BARMAN—Kingston University, London, United Kingdom
  • Prof Jenny BARRETT—University of Leeds, Yorkshire, United Kingdom
  • Prof Paul BISHOP—Manchester University, Manchester, United Kingdom
  • Mr Peter BLOWS—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Dr Catey BUNCE—King’s College London, London, United Kingdom
  • Dr Roxana CARARE—University of Southampton, Southampton, United Kingdom
  • Prof Usha CHAKRAVARTHY—Queens University Belfast, Belfast, Ireland
  • Miss Michelle CHAN—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation
  • Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Dr Sharon CHUA—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Prof David CRABB–City, University of London, London, United Kingdom
  • Mrs Philippa CUMBERLAND—UCL Great Ormond Street Institute of Child Health, London, United
  • Kingdom
  • Dr Alexander DAY—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Miss Parul DESAI—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Prof Bal DHILLON—University of Edinburgh, Scotland, United Kingdom
  • Prof Andrew DICK—University of Bristol, Bristol, United Kingdom
  • Dr Cathy EGAN—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Prof Sarah ENNIS—University of Southampton, Southampton, United Kingdom
  • Prof Paul FOSTER—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Dr Marcus FRUTTIGER—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Dr John GALLACHER—University of Oxford, Oxford, United Kingdom
  • Prof David (Ted) GARWAY-HEATH—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS
  • Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Dr Jane GIBSON—University of Southampton, Southampton, United Kingdom
  • Mr Dan GORE—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Prof Jeremy GUGGENHEIM—Cardiff University, Wales, United Kingdom
  • Prof Chris HAMMOND—King’s College London, London, United Kingdom
  • Prof Alison HARDCASTLE—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS
  • Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Prof Simon HARDING—University of Liverpool, London, United Kingdom
  • Dr Ruth HOGG—Queens University Belfast, Belfast, Ireland
  • Dr Pirro HYSI—King’s College London, London, United Kingdom*
  • Mr Pearse A KEANE—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Prof Sir Peng Tee KHAW—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS
  • Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Mr Anthony KHAWAJA—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Mr Gerassimos LASCARATOS—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS
  • Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Prof Andrew LOTERY- University of Southampton, Southampton, United Kingdom
  • Prof Phil LUTHERT—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Dr Tom MACGILLIVRAY—University of Edinburgh, Scotland, United Kingdom
  • Dr Sarah MACKIE—University of Leeds, Yorkshire, United Kingdom
  • Prof Keith MARTIN—University of Cambridge, Cambridge, United Kingdom
  • Ms Michelle MCGAUGHEY—Queen’s University Belfast, Belfast, Ireland
  • Dr Bernadette MCGUINNESS—Queen’s University Belfast, Belfast, Ireland
  • Dr Gareth MCKAY—Queen’s University Belfast, Belfast, Ireland
  • Mr Martin MCKIBBIN—Leeds Teaching Hospitals NHS Trust, Yorkshire, United Kingdom
  • Dr Danny MITRY–NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom & Royal Free Hospital, London, United Kingdom
  • Prof Tony MOORE—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Prof James MORGAN—Cardiff University, Wales, United Kingdom
  • Ms Zaynah MUTHY–NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Mr Eoin O’SULLIVAN—King’s College Hospital NHS Foundation Trust, London, United Kingdom
  • Dr Chris OWEN—St George’s, University of London, London, United Kingdom
  • Mr Praveen PATEL—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Mr Euan PATERSON—Queens University Belfast, Belfast, Ireland
  • Dr Tunde PETO—Queen’s University Belfast, Belfast, Ireland
  • Dr Axel PETZOLD—UCL Institute of Neurology, London, United Kingdom
  • Prof Jugnoo RAHI—UCL Great Ormond Street Institute of Child Health, London, United Kingdom
  • Dr Alicja RUDNICKA—St George’s, University of London, London, United Kingdom
  • Mr Jay SELF—University of Southampton, Southampton, United Kingdom
  • Prof Sobha SIVAPRASAD—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS
  • Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Mr David STEEL—Newcastle University, Newcastle, United Kingdom
  • Mrs Irene STRATTON—Gloucestershire Hospitals NHS Foundation Trust 5
  • Mr Nicholas STROUTHIDIS—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS
  • Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Prof Cathie SUDLOW—University of Edinburgh, Scotland, United Kingdom
  • Dr Caroline THAUNG—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Miss Dhanes THOMAS—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Prof Emanuele TRUCCO—University of Dundee, Scotland, United Kingdom
  • Prof Adnan TUFAIL—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation
  • Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Dr Veronique VITART—University of Edinburgh, Scotland, United Kingdom
  • Prof Stephen VERNON–Nottingham University Hospitals NHS Trust, Nottingham, United Kingdom
  • Mr Ananth VISWANATHAN—NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS
  • Foundation Trust and UCL Institute of Ophthalmology, London, United Kingdom
  • Dr Cathy WILLIAMS—University of Bristol, Bristol, United Kingdom
  • Dr Katie WILLIAMS—King’s College London, London, United Kingdom
  • Prof Jayne WOODSIDE—Queen’s University Belfast, Belfast, Ireland
  • Dr Max YATES—University of East Anglia, Norwich, United Kingdom
  • Ms Jennifer YIP—University of Cambridge, Cambridge, United Kingdom
  • Dr Yalin ZHENG—University of Liverpool, London, United Kingdom

*—lead author of the group; pirro.hysi@kcl.ac.uk

The CREAM Consortium

Joan E. Bailey-Wilson1*, Paul Nigel Baird2, Amutha Barathi Veluchamy3-5, Ginevra Biino6, Kathryn P. Burdon7, Harry Campbell8, Li Jia Chen9, Ching-Yu Cheng10-12, Emily Y. Chew13, Jamie E. Craig14, Phillippa M. Cumberland15, Margaret M. Deangelis16, Cécile Delcourt17, Xiaohu Ding18, Priya Duggal19, Cornelia M. van Duijn20, David M. Evans21-23, Qiao Fan24, Maurizio Fossarello25, Paul J. Foster26, Puya Gharahkhani27, Adriana I. Iglesias20,28,29, Jeremy A. Guggenheim30, Xiaobo Guo18,31, Annechien E. G. Haarman20,29, Toomas Haller32, Christopher J. Hammond33, Xikun Han27, Caroline Hayward34, Mingguang He2,18, Alex W. Hewitt2,7,35, Quan V. Hoang3,36, Pirro G. Hysi33, Robert P. Igo Jr.37, Sudha K. Iyengar37-39, Jost B. Jonas40,41, Mika Kähönen42,43, Jaakko Kaprio44,45, Anthony P. Khawaja26,46, Caroline C. W. Klaver16,20,29,47, Alison P. Klein19, Barbara E. Klein48, Ronald Klein48, Jonathan H. Lass37,38, Kris Lee48, Terho Lehtimäki49,50, Deyana Lewis1, Qing Li51, Shi-Ming Li41, Leo-Pekka Lyytikäinen49,50, Stuart MacGregor27, David A. Mackey2,7,35, Nicholas G. Martin52, Akira Meguro53, Andres Metspalu32, Candace Middlebrooks1, Masahiro Miyake54, Nobuhisa Mizuki53, Anthony Musolf1, Stefan Nickels55, Konrad Oexle56, Chi Pui Pang9, Olavi Pärssinen57,58, Andrew D. Paterson59, Craig Pennell35, Norbert Pfeiffer55, Ozren Polasek60,61, Jugnoo S. Rahi15,26,62, Olli Raitakari63,64, Igor Rudan8, Srujana Sahebjada2, Seang-Mei Saw65,66, Dwight Stambolian67, Claire L. Simpson1,68, E-Shyong Tai66, Milly S. Tedja20,29, J. Willem L. Tideman20,29, Akitaka Tsujikawa54, Virginie J. M. Verhoeven20,28,29, Veronique Vitart34, Ningli Wang41, Juho Wedenoja44,69, Wen Bin Wei70, Cathy Williams23, Katie M. Williams33, James F. Wilson8,34, Robert Wojciechowski1,19,71, Ya Xing Wang41, Kenji Yamashiro72, Jason C. S. Yam9, Maurice K. H. Yap73, Seyhan Yazar35, Shea Ping Yip74, Terri L. Young48, Xiangtian Zhou75

Affiliations

1. Computational and Statistical Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, Maryland, USA.

2. Centre for Eye Research Australia, Ophthalmology, Department of Surgery, University of Melbourne, Royal Victorian Eye and Ear Hospital, Melbourne, Australia.

3. Singapore Eye Research Institute, Singapore National Eye Centre, Singapore.

4. Duke-NUS Medical School, Singapore, Singapore.

5. Department of Ophthalmology, National University Health Systems, National University of Singapore, Singapore.

6. Institute of Molecular Genetics, National Research Council of Italy, Pavia, Italy.

7. Department of Ophthalmology, Menzies Institute of Medical Research, University of Tasmania, Hobart, Australia.

8. Centre for Global Health Research, Usher Institute for Population Health Sciences and Informatics, University of Edinburgh, Edinburgh, UK.

9. Department of Ophthalmology and Visual Sciences, The Chinese University of Hong Kong, Hong Kong Eye Hospital, Kowloon, Hong Kong.

10. Department of Ophthalmology, Yong Loo Lin School of Medicine, National University of Singapore, Singapore.

11. Ocular Epidemiology Research Group, Singapore Eye Research Institute, Singapore National Eye Centre, Singapore.

12. Ophthalmology & Visual Sciences Academic Clinical Program (Eye ACP), Duke-NUS Medical School, Singapore.

13. Division of Epidemiology and Clinical Applications, National Eye Institute/National Institutes of Health, Bethesda, USA.

14. Department of Ophthalmology, Flinders University, Adelaide, Australia.

15. Great Ormond Street Institute of Child Health, University College London, London, UK.

16. Institute for Molecular and Clinical Ophthalmology, Basel, Switzerland.

17. Université de Bordeaux, Inserm, Bordeaux Population Health Research Center, team LEHA, UMR 1219, F-33000 Bordeaux, France.

18. State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou, China.

19. Department of Epidemiology and Medicine, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland, USA.

20. Department of Epidemiology, Erasmus Medical Center, Rotterdam, The Netherlands.

21. Translational Research Institute, University of Queensland Diamantina Institute, Brisbane, Queensland, Australia.

22. MRC Integrative Epidemiology Unit, University of Bristol, Bristol, UK.

23. Department of Population Health Sciences, Bristol Medical School, Bristol, UK.

24. Centre for Quantitative Medicine, DUKE-National University of Singapore, Singapore.

25. University Hospital ’San Giovanni di Dio’, Cagliari, Italy.

26. NIHR Biomedical Research Centre, Moorfields Eye Hospital NHS Foundation Trust and UCL Institute of Ophthalmology, London, UK.

27. Statistical Genetics, QIMR Berghofer Medical Research Institute, Brisbane, Australia.

28. Department of Clinical Genetics, Erasmus Medical Center, Rotterdam, The Netherlands.

29. Department of Ophthalmology, Erasmus Medical Center, Rotterdam, The Netherlands.

30. School of Optometry & Vision Sciences, Cardiff University, Cardiff, UK.

31. Department of Statistical Science, School of Mathematics, Sun Yat-Sen University, Guangzhou, China.

32. Institute of Genomics, University of Tartu, Tartu, Estonia.

33. Section of Academic Ophthalmology, School of Life Course Sciences, King’s College London, London, UK.

34. MRC Human Genetics Unit, MRC Institute of Genetics & Molecular Medicine, University of Edinburgh, Edinburgh, UK.

35. Centre for Ophthalmology and Visual Science, Lions Eye Institute, University of Western Australia, Perth, Australia.

36. Department of Ophthalmology, Columbia University, New York, USA.

37. Department of Population and Quantitative Health Sciences, Case Western Reserve University, Cleveland, Ohio, USA.

38. Department of Ophthalmology and Visual Sciences, Case Western Reserve University and University Hospitals Eye Institute, Cleveland, Ohio, USA.

39. Department of Genetics, Case Western Reserve University, Cleveland, Ohio, USA.

40. Department of Ophthalmology, Medical Faculty Mannheim of the Ruprecht-Karls-University of Heidelberg, Mannheim, Germany.

41. Beijing Tongren Eye Center, Beijing Tongren Hospital, Beijing Institute of Ophthalmology, Beijing Key Laboratory of Ophthalmology and Visual Sciences, Capital Medical University, Beijing, China.

42. Department of Clinical Physiology, Tampere University Hospital and School of Medicine, University of Tampere, Tampere, Finland.

43. Finnish Cardiovascular Research Center, Faculty of Medicine and Life Sciences, University of Tampere, Tampere, Finland.

44. Department of Public Health, University of Helsinki, Helsinki, Finland.

45. Institute for Molecular Medicine Finland FIMM, HiLIFE Unit, University of Helsinki, Helsinki, Finland.

46. Department of Public Health and Primary Care, University of Cambridge, Cambridge, UK.

47. Department of Ophthalmology, Radboud University Medical Center, Nijmegen, The Netherlands.

48. Department of Ophthalmology and Visual Sciences, University of Wisconsin–Madison, Madison, Wisconsin, USA.

49. Department of Clinical Chemistry, Finnish Cardiovascular Research Center-Tampere, Faculty of Medicine and Life Sciences, University of Tampere.

50. Department of Clinical Chemistry, Fimlab Laboratories, University of Tampere, Tampere, Finland.

51. National Human Genome Research Institute, National Institutes of Health, Baltimore, USA.

52. Genetic Epidemiology, QIMR Berghofer Medical Research Institute, Brisbane, Australia.

53. Department of Ophthalmology, Yokohama City University School of Medicine, Yokohama, Kanagawa, Japan.

54. Department of Ophthalmology and Visual Sciences, Kyoto University Graduate School of Medicine, Kyoto, Japan.

55. Department of Ophthalmology, University Medical Center of the Johannes Gutenberg University Mainz, Mainz, Germany.

56. Institute of Neurogenomics, Helmholtz Zentrum München, German Research Centre for Environmental Health, Neuherberg, Germany.

57. Department of Ophthalmology, Central Hospital of Central Finland, Jyväskylä, Finland.

58. Gerontology Research Center, Faculty of Sport and Health Sciences, University of Jyväskylä, Jyväskylä, Finland.

59. Program in Genetics and Genome Biology, Hospital for Sick Children and University of Toronto, Toronto, Ontario, Canada.

60. Gen-info Ltd, Zagreb, Croatia.

61. University of Split School of Medicine, Soltanska 2, Split, Croatia.

62. Ulverscroft Vision Research Group, University College London, London, UK.

63. Research Centre of Applied and Preventive Cardiovascular Medicine, University of Turku, Turku, Finland.

64. Department of Clinical Physiology and Nuclear Medicine, Turku University Hospital, Turku, Finland.

65. Myopia Research Group, Singapore Eye Research Institute, Singapore National Eye Centre, Singapore.

66. Saw Swee Hock School of Public Health, National University Health Systems, National University of Singapore, Singapore.

67. Department of Ophthalmology, University of Pennsylvania, Philadelphia, Pennsylvania, USA.

68. Department of Genetics, Genomics and Informatics, University of Tennessee Health Sciences Center, Memphis, Tenessee.

69. Department of Ophthalmology, University of Helsinki and Helsinki University Hospital, Helsinki, Finland.

70. Beijing Tongren Eye Center, Beijing Key Laboratory of Intraocular Tumor Diagnosis and Treatment, Beijing Ophthalmology & Visual Sciences Key Lab, Beijing Tongren Hospital, Capital Medical University, Beijing, China.

71. Wilmer Eye Institute, Johns Hopkins Medical Institutions, Baltimore, Maryland, USA.

72. Department of Ophthalmology, Otsu Red Cross Hospital, Nagara, Japan.

73. Centre for Myopia Research, School of Optometry, The Hong Kong Polytechnic University, Hong Kong, Hong Kong.

74. Department of Health Technology and Informatics, The Hong Kong Polytechnic University, Hong Kong, Hong Kong.

75. School of Ophthalmology and Optometry, Eye Hospital, Wenzhou Medical University, China.

*—lead author of the group; jebw@mail.nih.gov

References

  1. 1. Saw S-M, Gazzard G, Shih‐Yen EC, Chua W-H. Myopia and associated pathological complications. Ophthalmic Physiol Opt. 2005;25: 381–391. pmid:16101943
  2. 2. Wojciechowski R. Nature and nurture: the complex genetics of myopia and refractive error. Clin Genet. 2011;79: 301–320. pmid:21155761
  3. 3. Hysi PG, Choquet H, Khawaja AP, Wojciechowski R, Tedja MS, Yin J, et al. Meta-analysis of 542,934 subjects of European ancestry identifies new genes and mechanisms predisposing to refractive error and myopia. Nat Genet. 2020;52: 401–407. pmid:32231278
  4. 4. Sanfilippo PG, Hewitt AW, Hammond CJ, Mackey DA. The heritability of ocular traits. Surv Ophthalmol. 2010;55: 561–583. pmid:20851442
  5. 5. Cumberland PM, Bao Y, Hysi PG, Foster PJ, Hammond CJ, Rahi JS, et al. Frequency and Distribution of Refractive Error in Adult Life: Methodology and Findings of the UK Biobank Study. PLOS ONE. 2015;10: e0139780. pmid:26430771
  6. 6. Van Hout CV, Tachmazidou I, Backman JD, Hoffman JD, Liu D, Pandey AK, et al. Exome sequencing and characterization of 49,960 individuals in the UK Biobank. Nature. 2020;586: 749–756. pmid:33087929
  7. 7. Tedja MS, Wojciechowski R, Hysi PG, Eriksson N, Furlotte NA, Verhoeven VJM, et al. Genome-wide association meta-analysis highlights light-induced signaling as a driver for refractive error. Nat Genet. 2018;50: 834–848. pmid:29808027
  8. 8. Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018;562: 203–209. pmid:30305743
  9. 9. Lee S, Wu MC, Lin X. Optimal tests for rare variant effects in sequencing association studies. Biostatistics. 2012;13: 762–775. pmid:22699862
  10. 10. Zhan X, Hu Y, Li B, Abecasis GR, Liu DJ. RVTESTS: an efficient and comprehensive tool for rare variant association analysis using sequence data. Bioinformatics. 2016;32: 1423–1426. pmid:27153000
  11. 11. Quick C, Wen X, Abecasis G, Boehnke M, Kang HM. Integrating comprehensive functional annotations to boost power and accuracy in gene-based association analysis. PLOS Genet. 2020;16: e1009060. pmid:33320851
  12. 12. Verhoeven VJM, Hysi PG, Wojciechowski R, Fan Q, Guggenheim JA, Höhn R, et al. Genome-wide meta-analyses of multiancestry cohorts identify multiple new susceptibility loci for refractive error and myopia. Nat Genet. 2013;45: 314–318. pmid:23396134
  13. 13. Vergara C, Bomotti SM, Valencia C, Klein BEK, Lee KE, Klein R, et al. Association analysis of exome variants and refraction, axial length, and corneal curvature in a European-American population. Hum Mutat. 2018;39: 1973–1979. pmid:30157304
  14. 14. Bomotti S, Lau B, Klein BEK, Lee KE, Klein R, Duggal P, et al. Refraction and Change in Refraction Over a 20-Year Period in the Beaver Dam Eye Study. Invest Ophthalmol Vis Sci. 2018;59: 4518–4524. pmid:30208419
  15. 15. Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature. 2020;581: 434–443. pmid:32461654
  16. 16. Slavotinek AM. Eye development genes and known syndromes. Mol Genet Metab. 2011;104: 448–456. pmid:22005280
  17. 17. Gallardo ME, Lopez-Rios J, Fernaud-Espinosa I, Granadino B, Sanz R, Ramos C, et al. Genomic Cloning and Characterization of the Human Homeobox Gene SIX6 Reveals a Cluster of SIX Genes in Chromosome 14 and Associates SIX6 Hemizygosity with Bilateral Anophthalmia and Pituitary Anomalies. Genomics. 1999;61: 82–91. pmid:10512683
  18. 18. Ledford KL, Martinez-De Luna RI, Theisen MA, Rawlins KD, Viczian AS, Zuber ME. Distinct cis-acting regions control six6 expression during eye field and optic cup stages of eye formation. Dev Biol. 2017;426: 418–428. pmid:28438336
  19. 19. Carnes MU, Liu YP, Allingham RR, Whigham BT, Havens S, Garrett ME, et al. Discovery and functional annotation of SIX6 variants in primary open-angle glaucoma. PLoS Genet. 2014;10: e1004372. pmid:24875647
  20. 20. Iglesias AI, Springelkamp H, van der Linde H, Severijnen L-A, Amin N, Oostra B, et al. Exome sequencing and functional analyses suggest that SIX6 is a gene involved in an altered proliferation–differentiation balance early in life and optic nerve degeneration at old age. Hum Mol Genet. 2014;23: 1320–1332. pmid:24150847
  21. 21. Wong TY, Klein BEK, Klein R, Knudtson M, Lee KE. Refractive errors, intraocular pressure, and glaucoma in a white population. Ophthalmology. 2003;110: 211–217. pmid:12511368
  22. 22. Mastropasqua L, Lobefalo L, Mancini A, Ciancaglini M, Palma S. Prevalence of Myopia in Open Angle Glaucoma. Eur J Ophthalmol. 1992;2: 33–35. pmid:1638164
  23. 23. Griffith JF, DeBenedictis MJ, Traboulsi EI. A novel dominant CRX mutation causes adult-onset macular dystrophy. Ophthalmic Genet. 2018;39: 120–124. pmid:28945142
  24. 24. Furukawa T, Morrow EM, Cepko CL. Crx, a Novel otx-like Homeobox Gene, Shows Photoreceptor-Specific Expression and Regulates Photoreceptor Differentiation. Cell. 1997;91: 531–541. pmid:9390562
  25. 25. Andzelm MM, Cherry TJ, Harmin DA, Boeke AC, Lee C, Hemberg M, et al. MEF2D Drives Photoreceptor Development through a Genome-wide Competition for Tissue-Specific Enhancers. Neuron. 2015;86: 247–263. pmid:25801704
  26. 26. Sohocki MM, Sullivan LS, Mintz-Hittner HA, Birch D, Heckenlively JR, Freund CL, et al. A Range of Clinical Phenotypes Associated with Mutations in CRX, a Photoreceptor Transcription-Factor Gene. Am J Hum Genet. 1998;63: 1307–1315. pmid:9792858
  27. 27. Wu Y, Qiu J, Chen S, Chen X, Zhang J, Zhuang J, et al. Crx Is Posttranscriptionally Regulated by Light Stimulation in Postnatal Rat Retina. Front Cell Dev Biol. 2020;0. pmid:32318566
  28. 28. Dryja TP, Berson EL, Rao VR, Oprian DD. Heterozygous missense mutation in the rhodopsin gene as a cause of congenital stationary night blindness. Nat Genet. 1993;4: 280–283. pmid:8358437
  29. 29. Craft CM, Deming JD. Cone Arrestin: Deciphering the Structure and Functions of Arrestin 4 in Vision. In: Gurevich VV, editor. Arrestins—Pharmacology and Therapeutic Potential. Berlin, Heidelberg: Springer; 2014. pp. 117–131. https://doi.org/10.1007/978-3-642-41199-1_6 pmid:24292827
  30. 30. Ugarte M, Pérez-Cerdá C, Rodríguez-Pombo P, Desviat LR, Pérez B, Richard E, et al. Overview of mutations in the PCCA and PCCB genes causing propionic acidemia. Hum Mutat. 1999;14: 275–282. https://doi.org/10.1002/(SICI)1098-1004(199910)14:4<275::AID-HUMU1>3.0.CO;2-N pmid:10502773
  31. 31. Li X, Zhang J, Wang Y, Ji J, Yang F, Wan C, et al. Association study on the NAPG gene and bipolar disorder in the Chinese Han population. Neurosci Lett. 2009;457: 159–162. pmid:19429185
  32. 32. Zhao B, Shan Y, Yang Y, Yu Z, Li T, Wang X, et al. Transcriptome-wide association analysis of brain structures yields insights into pleiotropy with complex neuropsychiatric traits. Nat Commun. 2021;12: 2878. pmid:34001886
  33. 33. Bulayeva K, Lesch K-P, Bulayev O, Walsh C, Glatt S, Gurgenova F, et al. Genomic structural variants are linked with intellectual disability. J Neural Transm. 2015;122: 1289–1301. pmid:25626716
  34. 34. Goes FS, McGrath J, Avramopoulos D, Wolyniec P, Pirooznia M, Ruczinski I, et al. Genome-wide association study of schizophrenia in Ashkenazi Jews. Am J Med Genet B Neuropsychiatr Genet. 2015;168: 649–659. pmid:26198764
  35. 35. Han M-R, Schellenberg GD, Wang L-S, Alzheimer’s Disease Neuroimaging Initiative. Genome-wide association reveals genetic effects on human Aβ42 and τ protein levels in cerebrospinal fluids: a case control study. BMC Neurol. 2010;10: 90. pmid:20932310
  36. 36. Cui L, Tao H, Wang Y, Liu Z, Xu Z, Zhou H, et al. A functional polymorphism of the microRNA-146a gene is associated with susceptibility to drug-resistant epilepsy and seizures frequency. Seizure. 2015;27: 60–65. pmid:25891929
  37. 37. Luo M, Ma W, Sand Z, Finlayson J, Wang T, Brinton RD, et al. Von Willebrand factor A domain-containing protein 8 (VWA8) localizes to the matrix side of the inner mitochondrial membrane. Biochem Biophys Res Commun. 2020;521: 158–163. pmid:31630795
  38. 38. Anney R, Klei L, Pinto D, Regan R, Conroy J, Magalhaes TR, et al. A genome-wide scan for common alleles affecting risk for autism. Hum Mol Genet. 2010;19: 4072–4082. pmid:20663923
  39. 39. Greenwood TA, Akiskal HS, Akiskal KK, Bipolar Genome Study, Kelsoe JR. Genome-wide association study of temperament in bipolar disorder reveals significant associations with three novel Loci. Biol Psychiatry. 2012;72: 303–310. pmid:22365631
  40. 40. Lee JJ, Wedow R, Okbay A, Kong E, Maghzian O, Zacher M, et al. Gene discovery and polygenic prediction from a 1.1-million-person GWAS of educational attainment. Nat Genet. 2018;50: 1112–1121. pmid:30038396
  41. 41. Pitts MW, Byrns CN, Ogawa-Wong AN, Kremer P, Berry MJ. Selenoproteins in Nervous System Development and Function. Biol Trace Elem Res. 2014;161: 231–245. pmid:24974905
  42. 42. Dekker AM, Diekstra FP, Pulit SL, Tazelaar GHP, van der Spek RA, van Rheenen W, et al. Exome array analysis of rare and low frequency variants in amyotrophic lateral sclerosis. Sci Rep. 2019;9. pmid:30976013
  43. 43. Verheijen J, Van den Bossche T, van der Zee J, Engelborghs S, Sanchez-Valle R, Lladó A, et al. A comprehensive study of the genetic impact of rare variants in SORL1 in European early-onset Alzheimer’s disease. Acta Neuropathol (Berl). 2016;132: 213–224. pmid:27026413
  44. 44. Liskova P, Dudakova L, Evans CJ, Rojas Lopez KE, Pontikos N, Athanasiou D, et al. Ectopic GRHL2 Expression Due to Non-coding Mutations Promotes Cell State Transition and Causes Posterior Polymorphous Corneal Dystrophy 4. Am J Hum Genet. 2018;102: 447–459. pmid:29499165
  45. 45. Maróti Z, Boldogkői Z, Tombácz D, Snyder M, Kalmár T. Evaluation of whole exome sequencing as an alternative to BeadChip and whole genome sequencing in human population genetic analysis. BMC Genomics. 2018;19: 778. pmid:30373510