Identification of four novel variants that influence central corneal thickness in multi-ethnic Asian populations

What intrinsic factors affect the central corneal thickness?

The cornea is one of the tissues responsible for covering and protecting the inner structures of the eye. Central corneal thickness (CCT) is defined as the distance between the anterior epithelial surface and the posterior surface of the endothelial layer. This parameter plays a very important role regarding intraocular pressure (IOP) measurement, evaluation of corneal uniformity, selection of a suitable technique for corneal refractive surgery and the planning of surgical procedures to overcome corneal disease. This comprehensive review elucidates the multifaceted factors influencing the central corneal thickness. Recognising the impact of these factors not only enhances our understanding of corneal dynamics but also contributes significantly to the refinement of diagnostic and therapeutic strategies in ophthalmology.

Relationship Between Retinal Layer Thickness and Genetic Susceptibility to Age-Related Macular Degeneration in Asian Populations

Purpose

For OCT retinal thickness measurements to be used as a prodromal age-related macular degeneration (AMD) risk marker, the 3-dimensional (3D) topographic variation of the relationship between genetic susceptibility to AMD and retinal thickness needs to be assessed. We aimed to evaluate individual retinal layer thickness changes and topography at the macula that are associated with AMD genetic susceptibility.

Design

Genetic association study.

Participants

A total of 1579 healthy participants (782 Chinese, 353 Malays, and 444 Indians) from the multiethnic Singapore Epidemiology of Eye Diseases study were included.

Methods

Spectral-domain OCT and automatic segmentation of individual retinal layers were performed to produce 10 retinal layer thickness measurements at each ETDRS subfield, producing 3D topographic information. Age-related macular degeneration genetic susceptibility was represented via single nucleotide polymorphisms (SNPs) and aggregated via whole genome (overall) and pathway-specific age-related macular degeneration polygenic risk score (PRSAMD).

Main Outcome Measures

Associations of individual SNPs, overall PRSAMD, and pathway-specific PRSAMD with retinal thickness were analyzed by individual retinal layer and ETDRS subfield.

Results

CFH rs10922109, ARMS2-HTRA1 rs3750846, and LIPC rs2043085 were the top AMD susceptibility SNPs associated with retinal thickness of individual layers (P < 1.67 × 10-3), all at the central subfield. The overall PRSAMD was most associated with thinner L9 (outer segment photoreceptor/retinal pigment epithelium complex) thickness at the central subfield (β = -0.63 μm; P = 5.45 × 10-9). Pathway-specific PRSAMD for the complement cascade (β = -0.53 μm; P = 9.42 × 10-7) and lipoprotein metabolism (β = -0.05 μm; P = 0.0061) were associated with thinner photoreceptor layers (L9 and L7 [photoreceptor inner/outer segments], respectively) at the central subfield. The mean PRSAMD score was larger among Indians compared with that of the Chinese and had the thinnest thickness at the L9 central subfield (β = -1.00 μm; P = 2.91 × 10-7; R2 = 5.5%). Associations at other retinal layers and ETDRS regions were more heterogeneous.

Conclusions

Overall genetic susceptibility to AMD and the aggregate effects of the complement cascade and lipoprotein metabolism pathway are associated most significantly with L7 and L9 photoreceptor thinning at the central macula in healthy individuals. Photoreceptor thinning has potential to be a prodromal AMD risk marker, and topographic variation should be considered.

Financial Disclosures

Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Genetic prescreening of a candidate for laser refractive surgery identifies risk for inadequate tissue response: a case report

Background

Inadequate response to corneal laser refractive surgery, e.g., ectatic corneal diseases, may not be identified by conventional examinations, hence creating therapeutic uncertainty. Herein we demonstrate the application of genetic prescreening to augment preassessment for corneal laser refractive surgery and highlight the ability to prevent the possibility of enrolling a subject at risk for developing ectatic corneal diseases.

Case presentation

Preoperative tests were performed alongside deoxyribonucleic acid (DNA) sequencing of 75 genes specific to the structure and health of the eye of a 44-year-old Caucasian male candidate for corneal laser refractive surgery. The patient had no medical, family, or psychosocial history, nor symptoms that could lead to suspect any corneal abnormalities, and conventional preoperative tests confirmed that no corneal abnormalities were present. The sequencing results uncovered rare DNA variants within the ADGRV1, PTK2, ZNF469, and KRT15 genes. These variants were considered potential risk factors for inadequate response in the patient post corneal laser refractive surgery. Subsequent reevaluation with three different last-generation corneal tomographers identified in the left eye a “warning” for a deformity of the posterior profile of the cornea.

Conclusions

Genetic prescreening identifies potential risk of inadequate response to corneal laser refractive surgery where current technologies in use may lead to a hazardous predictive diagnostic uncertainty.

Evaluating the association between MPDZ-NF1B rs1324183 and keratoconus in an independent northwestern Chinese population

Background

Keratoconus (KC) is a complex, non-inflammatory corneal degenerative disease. Although numerous studies have analyzed the correlation of SNP rs1324183, which located in MPDZ-NF1B gene, and KC in different populations, only few findings were repeated. In this study, to evaluate the association between rs1324183 and KC in a new independent Chinese population, we performed a replication study of the significantly associated rs1324183.

Methods

In total of 114 unrelated KC patients and 88 unrelated controls were recruited from Ningxia, China. We detected the genotypes and alleles of rs1324183 using PCR technology and Sanger sequencing and also analyzed the association between this locus and KC, its clinical parameters by statistical methods.

Results

The frequency of genotype AA (11, 9.6%) and genotypes containing allele A (47, 41.2%) of rs1324183 in KC were both higher than those of the control group. And genotype AA of rs1324183 conferred a higher risk of KC (OR > 1). Moreover, corneal parameter Belin/Ambrósio enhanced ectasia display final D value (BAD-D) had significant correlation (p = 0.002) with AA genotype of rs1324183 in KC.

Conclusions

Our replication study indicates that the results of rs1324183 associated with KC in our population is robust and further better illustrates the significance of BAD-D as a diagnostic indicator for KC. rs1324183 should be considered as the first genetic mark of KC risk in its future diagnosis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12886-022-02359-1.

A Precise Method to Evaluate 360 Degree Measures of Optic Cup and Disc Morphology in an African American Cohort and Its Genetic Applications

(1) Background: Vertical cup-to-disc ratio (CDR) is an important measure for evaluating damage to the optic nerve head (ONH) in glaucoma patients. However, this measure often does not fully capture the irregular cupping observed in glaucomatous nerves. We developed and evaluated a method to measure cup-to-disc ratio (CDR) at all 360 degrees of the ONH. (2) Methods: Non-physician graders from the Scheie Reading Center outlined the cup and disc on digital stereo color disc images from African American patients enrolled in the Primary Open-Angle African American Glaucoma Genetics (POAAGG) study. After converting the resultant coordinates into polar representation, the CDR at each 360-degree location of the ONH was obtained. We compared grader VCDR values with clinical VCDR values, using Spearman correlation analysis, and validated significant genetic associations with clinical VCDR, using grader VCDR values. (3) Results: Graders delineated outlines of the cup contour and disc boundaries twice in each of 1815 stereo disc images. For both cases and controls, the mean CDR was highest at the horizontal bisector, particularly in the temporal region, as compared to other degree locations. There was a good correlation between grader CDR at the vertical bisector and clinical VCDR (Spearman Correlation OD: r = 0.78 [95% CI: 0.76–0.79]). An SNP in the MPDZ gene, associated with clinical VCDR in a prior genome-wide association study, showed a significant association with grader VCDR (p = 0.01) and grader CDR area ratio (p = 0.02). (4) Conclusions: The CDR of both glaucomatous and non-glaucomatous eyes varies by degree location, with the highest measurements in the temporal region of the eye. This method can be useful for capturing innate eccentric ONH morphology, tracking disease progression, and identifying genetic associations.

Rare single nucleotide variants in COL5A1 promoter do not play a major role in keratoconus susceptibility associated with rs1536482

BACKGROUND

Keratoconus is a chronic degenerative disorder of the cornea characterized by thinning and cone-shaped protrusions. Although genetic factors play a key role in keratoconus development, the etiology is still under investigation. The occurrence of single-nucleotide polymorphisms (SNPs) associated with keratoconus in Russian patients is poorly studied. The purpose of this study was to validate whether three reported keratoconus-associated SNPs (rs1536482 near the COL5A1 gene, rs2721051 near the FOXO1 gene, rs1324183 near the MPDZ gene) are also actual for a Russian cohort of patients. Additionally, we investigated the COL5A1 promoter sequence for single-nucleotide variants (SNVs) in a subgroup of keratoconus patients with at least one rs1536482 minor allele (rs1536482+) to assess the role of these SNVs in keratoconus susceptibility associated with rs1536482.

METHODS

This case-control study included 150 keratoconus patients and two control groups (main and additional, 205 and 474 participants, respectively). We performed PCR targeting regions flanking SNVs and the COL5A1 promoter, followed by Sanger sequencing of amplicons. The additional control group was genotyped using an SNP array.

RESULTS

The minor allele frequency was significantly different between the keratoconus and control cohorts (main and combined) for rs1536482, rs2721051, and rs1324183 (p-value < 0.05). The rare variants rs1043208782 and rs569248712 were found in the COL5A1 promoter in two out of 94 rs1536482+ keratoconus patients.

CONCLUSION

rs1536482, rs2721051, and rs1324183 were associated with keratoconus in a Russian cohort. SNVs in the COL5A1 promoter do not play a major role in keratoconus susceptibility associated with rs1536482.

Mendelian randomization analysis identified genes pleiotropically associated with central corneal thickness

OBJECTIVE

To prioritize genes that were pleiotropically or potentially causally associated with central corneal thickness (CCT).

METHODS

We applied the summary data-based Mendelian randomization (SMR) method integrating summarized data of genome-wide association study (GWAS) on CCT and expression quantitative trait loci (eQTL) data to identify genes that were pleiotropically associated with CCT. We performed separate SMR analysis using CAGE eQTL data and GTEx eQTL data. SMR analyses were done for participants of European and East Asian ancestries, separately.

RESULTS

We identified multiple genes showing pleiotropic association with CCT in the participants of European ancestry. CLIC3 (ILMN_1796423; PSMR = 4.15 × 10- 12), PTGDS (ILMN_1664464; PSMR = 6.88 × 10- 9) and C9orf142 (ILMN_1761138; PSMR = 8.09 × 10- 9) were the top three genes using the CAGE eQTL data, and RP11-458F8.4 (ENSG00000273142.1; PSMR = 5.89 × 10- 9), LCNL1 (ENSG00000214402.6; PSMR = 5.67 × 10- 8), and PTGDS (ENSG00000107317.7; PSMR = 1.92 × 10- 7) were the top three genes using the GTEx eQTL data. No genes showed significantly pleiotropic association with CCT in the participants of East Asian ancestry after correction for multiple testing.

CONCLUSION

We identified several genes pleiotropically associated with CCT, some of which represented novel genes influencing CCT. Our findings provided important leads to a better understanding of the genetic factors influencing CCT, and revealed potential therapeutic targets for the treatment of primary open-angle glaucoma and keratoconus.

Molecular Genetics and Functional Analysis Implicate CDKN2BAS1-CDKN2B Involvement in POAG Pathogenesis

The genes in the 9p21 locus (CDKN2B-AS1 & CDKN2B) are widely associated with Primary open-angle glaucoma (POAG). However, the functional importance of this locus in POAG pathogenesis is still unexplored. This study investigated the role of CDKN2BAS1-CDKN2B axis in POAG. We observed significant association of CDKN2B-AS1 SNP rs4977756 with POAG and its endophenotypic traits (vertical cup-disc ratio (p = 0.033) and central corneal thickness (p = 0.008)) by screening African American POAG cases (n = 1567) and controls (n = 1600). A luciferase reporter assay in Human embryonic kidney 293T (HEK293T) cells revealed that the region surrounding rs4977756 likely serves as a transcriptional repressor. siRNA-mediated knockdown of CDKN2B-AS1 in HEK293T cells and trabecular meshwork (TM) cells resulted in significantly increased expression of CDKN2B, which was also observed in human POAG ocular tissues. Pathway focused qRT-PCR gene expression analysis showed increased cellular senescence, TGFβ signaling and ECM deposition in TM cells after CDKN2B-AS1 suppression. In conclusion, we report that CDKN2B-AS1 may act as a regulator, and it could function by modulating the expression of CDKN2B. In addition, increase in CDKN2B levels due to CDKN2B-AS1 suppression may result in the senescence of trabecular meshwork cells leading to POAG pathogenesis.

Evaluation of Primary Angle-Closure Glaucoma Susceptibility Loci for Estimating Angle Closure Disease Severity

PURPOSE

To investigate whether recently identified genetic loci for primary angle-closure glaucoma (PACG) are associated with disease severity.

DESIGN

Case-control study.

PARTICIPANTS

Eight hundred four PACG patients and 943 control participants of Chinese ethnicity from Singapore.

METHODS

The 8 PACG-associated single nucleotide polymorphisms (SNPs; rs11024102 at PLEKHA7, rs3753841 at COL11A1, rs1015213 located between PCMTD1 and ST18 on chromosome 8q, rs3816415 at EPDR1, rs1258267 at CHAT, rs736893 at GLIS3, rs7494379 at FERMT2, and rs3739821 mapping in between DPM2 and FAM102A) identified from genome-wide association studies were tested for association with disease severity using logistic regression adjusted for age and gender. A P value of 0.006 was set as significant after Bonferroni correction for testing of 8 loci. We also calculated the weighted genetic risk score (GRS) weighted by the estimated individual SNP effect size on PACG calculated as logarithm of the odds ratio (OR). Disease severity was based on the visual field mean deviation (MD) and classified as early to moderate (MD, >-12 dB) and severe (MD, <-20 dB).

MAIN OUTCOME MEASURES

Association of PACG loci with severe disease.

RESULTS

Of the 804 PACG patients, genotyping data were available for 768 individuals and included 436 with mild-to-moderate PACG and 206 with severe PACG. The PACG patients were significantly older (mean age, 64.3 ± 9.1 years vs. 56.4 ± 8.9 years; P < 0.001) and there were proportionately more women compared with control participants (58.4% vs. 49.0%; P < 0.001). Of the 8 loci investigated, we observed significant evidence of association with severe PACG at 1 SNP, namely rs3816415 in EPDR1 (OR, 2.03; 95% confidence interval [CI], 1.49-2.78; P = 1 × 10^-5). A higher-weighted GRS was associated significantly with severe PACG, with an OR of 3.11 (95% CI, 1.95-4.96) comparing the lowest quartile with the highest quartile.

CONCLUSIONS

Our results show that EPDR1 is associated significantly with severe PACG, suggesting that it may predispose patients to more aggressive disease development. Individuals with PACG with a higher GRS were associated with a higher risk of severe PACG.

Genome-Wide Association of Genetic Variants With Refraction, Axial Length, and Corneal Curvature: A Longitudinal Study of Chinese Schoolchildren

Background

Myopia is a common eye disorder that is approaching epidemic proportions worldwide. A genome-wide association study identified AREG (rs12511037), GABRR1 (rs13215566), and PDE10A (rs12206610) as being associated with refractive error in Asian populations. The present study investigated the associations between these three genetic variants and the occurrence and development of myopia, spherical equivalent refraction (SER), axial length (AL), and corneal curvature (CC) in a cohort of southeastern Chinese schoolchildren.

Methods

We examined and followed 550 children in grade 1 enrolled in the Wenzhou Epidemiology of Refractive Error (WERE) project. During the 4-year follow-up, non-cycloplegic refraction was evaluated twice each year, and the AL and CC were measured once every year. Age, sex, and the amounts of time spent on near work and outdoors were documented with a questionnaire. Sanger DNA sequencing was used to genotype single nucleotide polymorphisms (SNPs). SNPtest software was used to identify potential genetic variants associated with myopia, SER, AL, and CC. Ten thousand permutations were used to correct for multiple testing.

Results

In total, 469 children, including 249 (53.1%) boys and 220 (46.9%) girls, were included in analyses. The mean age of all the children was 6.33 ± 0.48 years. After adjusting for age, sex, time spent on near work and time spent outdoors, neither the genotypes nor the allele frequencies of the three SNPs were significantly associated with myopic shift, incident myopia or the change in SER. After adjusting for age, sex, near-work time and outdoor time with 10,000 permutations, the genotype AREG (rs12511037) was associated with an increase in AL (P′-values for the dominant, recessive, additive and general models were 0.0032, 0.0275, 0.0045, and 0.0099, respectively); the genotype PDE10A (rs12206610) was associated with a change in CC in the additive (P′ = 0.0096), dominant (P′ = 0.0096), and heterozygous models (P′ = 0.0096).

Conclusion

These findings preliminarily indicate that AREG SNP rs12511037 and PDE10A SNP rs12206610 are etiologically relevant for ocular traits, providing a basis for further exploration of the development of myopia and its molecular mechanism. However, elucidating the role of AREG and PDE10A in the pathogenesis of myopia requires further animal model and human genetic epidemiology studies. This trial is registered as ChiCTR1900020584 at www.Chictr.org.cn.

Integration of Genetic and Biometric Risk Factors for Detection of Primary Angle Closure Glaucoma

PURPOSE

The purpose of this study was to investigate whether the addition of primary angle closure glaucoma (PACG)-associated genetic loci allows improved detection of PACG, compared to anterior segment parameters measured by imaging.

DESIGN

Case-control study.

METHODS

Genotype data of the 8 PACG single-nucleotide polymorphisms (SNPs) (rs11024102 at PLEKHA7, rs3753841 at COL11A1, rs1015213 located between PCMTD1 and ST18 on Chromosome 8q, rs3816415 at EPDR1, rs1258267 at CHAT, rs736893 at GLIS3, rs7494379 at FERMT2, and rs3739821 mapping in between DPM2 and FAM102A) were available. Customized software was used to measure anterior segment optical coherence tomography (ASOCT) parameters, namely, anterior chamber depth, width, and area (ACD, ACW, and ACA) and lens vault (LV). Statistical analysis for positive predictive values was modeled using the receiver operating characteristic curve (AUC). Statistical significance comparing predictive power of the different parameters was calculated using permutation.

RESULTS

A total of 388 PACG subjects and 751 controls with both ASOCT and genetic data were available for analysis. Anterior segment parameters including ACD, ACA, and LV had excellent predictive value (AUCs >0.94), except ACW (AUC=0.65), for identifying PACG. The inclusion of genetic risk alleles (either singly or as a composite genetic risk score for 8 genomewide association study SNPs) to ACD only provided a +0.50% improvement in reclassifying PACG cases and controls over and above the discriminatory value of ACD. This +0.50% improvement was not statistically significant (P > .05).

CONCLUSIONS

Although significant on their own, the incorporation of genetic data in the context of anterior segment imaging parameters like ACD provided only a marginal improvement of PACG detection.

Genetic Aspects of Keratoconus: A Literature Review Exploring Potential Genetic Contributions and Possible Genetic Relationships with Comorbidities

INTRODUCTION

Keratoconus (KC) is a complex, genetically heterogeneous, multifactorial degenerative disorder that is accompanied by corneal ectasia which usually progresses asymmetrically. With an incidence of approximately 1 per 2000 and 2 cases per 100,000 population presenting annually, KC follows an autosomal recessive or dominant pattern of inheritance and is, apparently, associated with genes that interact with environmental, genetic, and/or other factors. This is an important consideration in refractive surgery in the case of familial KC, given the association of KC with other genetic disorders and the imbalance between dizygotic twins. The present review attempts to identify the genetic loci contributing to the different KC clinical presentations and relate them to the common genetically determined comorbidities associated with KC.

METHODS

The PubMed, MEDLINE, Google Scholar, and GeneCards databases were screened for KC-related articles published in English between January 2006 and November 2017. Keyword combinations of "keratoconus," "risk factor(s)," "genetics," "genes," "genetic association(s)," and "cornea" were used. In total, 217 articles were retrieved and analyzed, with greater weight placed on the more recent literature. Further bibliographic research based on the 217 articles revealed another 124 relevant articles that were included in this review. Using the reviewed literature, an attempt was made to correlate genes and genetic risk factors with KC characteristics and genetically related comorbidities associated with KC based on genome-wide association studies, family-based linkage analysis, and candidate-gene approaches.

RESULTS

An association matrix between known KC-related genes and KC symptoms and/or clinical signs together with an association matrix between identified KC genes and genetically related KC comorbidities/syndromes were constructed.

CONCLUSION

Twenty-four genes were identified as potential contributors to KC and 49 KC-related comorbidities/syndromes were found. More than 85% of the known KC-related genes are involved in glaucoma, Down syndrome, connective tissue disorders, endothelial dystrophy, posterior polymorphous corneal dystrophy, and cataract.

Correlation of Color Fundus Photograph Grading with Risks of Early Age-related Macular Degeneration by using Automated OCT-derived Drusen Measurements

We evaluated automated OCT-derived drusen volume measures in a population-based study (n = 4,512) aged ≥40 years, and its correlation with conventional color fundus photographs (CFP)-derived early AMD features. Participants had protocol-based assessment to capture medical and ocular history, genotyping for SNPs in CFH, ARMS2, and CETP, CFP-based AMD grading and automated drusen volume based on SD-OCT using built-in software (Cirrus OCT advanced RPE analysis software). Significantly fewer eyes with early AMD features (drusen, hyperpigmentation, soft or reticular drusen) had drusen volume = 0 mm³ (p < 0.001). In eyes with drusen volume > 0 mm³, increasing AMD severity was associated with increase in drusen volume (correlation coefficient 0.17, p < 0.001). However 220 (59.14%) of 372 participants with AMD based on CFP grading had drusen volume = 0 mm³. Factors associated with drusen volume included age (OR 1.42 per 5 years, 95% confidence interval [CI] 2.76, 4.48), systolic blood pressure (OR1.00, 95% CI 1.00, 1.01), ethnic Malay (OR 1.54, 95% CI 1.29, 1.83) and Chinese (OR 1.66, 95% CI 1.37, 2.01) compared to Indian. The ARMS2 rs10490924 T allele was associated with increased drusen volume in subjects with AMD (multivariable adjusted OR1.54, 95% CI 1.08, 2.19). Automated OCT-derived drusen volume is correlated with CFP-based AMD grading in many, but not all subjects. However the agreement is not good. These two modalities provide complementary information and should be incorporated into future studies.

Functional relevance for central cornea thickness-associated genetic variants by using integrative analyses

Background

The genetic architecture underlying central cornea thickness (CCT) is far from understood. Most of the CCT-associated variants are located in the non-coding regions, raising the difficulty of following functional characterizations. Thus, integrative functional analyses on CCT-associated loci might benefit in overcoming these issues by prioritizing the hub genes that are located in the center of CCT genetic network.

Methods

Integrative analyses including functional annotations, enrichment analysis, and protein-protein interaction analyses were performed on all reported CCT GWAS lead SNPs, together with their proxy variants. Functional annotations were conducted by CADD, GWAVA, and Eigen. Enrichment analyses for CCT-associated genes were performed using ToppGene suite. Protein-protein interaction network and gene co-expression analyses were performed by GeneMANIA.

Results

Functional annotations prioritized eight genes (ADAMSTS6, ARID5B, FOXO1, AKAP13, COL4A3, COL8A2, TBL1XR1, and KCMB2) harboring SNPs with strong evidence of regulatory potential. It was also shown that CCT-associated genes were significantly enriched in collagen-related pathways and the phenotype of keratoconus, and some of them were found to be involved in one interaction network.

Conclusion

This study revealed the hub genes that were located in the center of CCT genetic network and provided a new insight into the genetic regulation underlying CCT GWAS findings.

Analysis of multiple genetic loci reveals MPDZ-NF1B rs1324183 as a putative genetic marker for keratoconus

OBJECTIVE

To investigate the associations between 16 single-nucleotide polymorphisms (SNPs) in 14 genetic loci and keratoconus in an independent Chinese cohort.

METHODS

This cross-sectional, case-control association study included a Chinese cohort of 133 patients with keratoconus and 371 control subjects. In a recent meta-analysis study, we identified association of 16 SNPs in 14 gene loci with keratoconus. In this study, we genotyped these 16 SNPs in all the patients and controls and analysed their association with keratoconus, its clinical severities and progression profiles. We also analysed the genotype-phenotype correlation between individual SNPs and steep keratometry, flat keratometry (Kf), average keratometry (Avg K) and best-fit sphere diameter (BFS) of the anterior and posterior corneal surface.

RESULTS

Among the 16 selected SNPs, rs1324183 in the MPDZ-NF1B locus showed a significant association with keratoconus (OR=2.22; 95% CI 1.42 to 3.45, p=4.30×10^-4), especially severe keratoconus (OR=5.10, 95% CI 1.63 to 15.93, p=0.005). The rs1324183 A allele was positively associated with anterior Kf (p=0.008), anterior Avg K (p=0.017), posterior Kf (p=0.01) and negatively associated with apex pachymetry (p=0.007) and anterior BFS (p=0.023) in keratoconus. The other 15 SNPs had no significant association with keratoconus or genotype-phenotype correlations.

CONCLUSIONS

This study confirmed the association of SNP rs1324183 in MPDZ-NF1B with keratoconus and revealed the association of this SNP with keratoconus severity and corneal parameters. It is thus a putative genetic marker for monitoring the progression of keratoconus to a severe form and facilitating early intervention.

The effects of collagen cross-linking on corneal density: a comparison between accelerated and conventional methods

PURPOSE

To compare the corneal density before and after accelerated versus conventional corneal cross-linking (CXL).

METHODS

Pentacam densitometry and corneal mapping were performed before and 6 months following the randomly assigned intervention. Corneal density was measured over different zones and layers of the cornea by two independent ophthalmologists. Total corneal density was calculated. The patients were followed up for the occurrence of procedure-related and unrelated complication 1 week, 1 month and a final visit around 6 months (180 ± 10 days) after the procedure, during which the corneal density was remeasured and recorded.

RESULTS

Twenty-three patients were enrolled. There was a significant increase in corneal density over the anterior layer (P < 0.001) and middle corneal layer (P = 0.035). Additionally, the corneal density significantly increased in 0-2 and 2-6 mm zones (P < 0.001), as well as total density (P = 0.002) following corneal CXL. Although there was no significant difference in the corneal density between the two CXL groups, the increase in corneal density was only significant in group A for the 6-10 mm zone and the posterior corneal layer.

CONCLUSION

Both the conventional and accelerated CXL methods increase the overall corneal density with exception of posterior central corneal layer and the 6-10 mm zone, which is only increased in the accelerated method.

Cross-ancestry genome-wide association analysis of corneal thickness strengthens link between complex and Mendelian eye diseases

Central corneal thickness (CCT) is a highly heritable trait associated with complex eye diseases such as keratoconus and glaucoma. We perform a genome-wide association meta-analysis of CCT and identify 19 novel regions. In addition to adding support for known connective tissue-related pathways, pathway analyses uncover previously unreported gene sets. Remarkably, >20% of the CCT-loci are near or within Mendelian disorder genes. These included FBN1, ADAMTS2 and TGFB2 which associate with connective tissue disorders (Marfan, Ehlers-Danlos and Loeys-Dietz syndromes), and the LUM-DCN-KERA gene complex involved in myopia, corneal dystrophies and cornea plana. Using index CCT-increasing variants, we find a significant inverse correlation in effect sizes between CCT and keratoconus (r = -0.62, P = 5.30 × 10-5) but not between CCT and primary open-angle glaucoma (r = -0.17, P = 0.2). Our findings provide evidence for shared genetic influences between CCT and keratoconus, and implicate candidate genes acting in collagen and extracellular matrix regulation.

The thick and thin of the central corneal thickness in glaucoma

Central corneal thickness (CCT) is an important parameter in the assessment of any potential glaucoma patient. While it affects prognosis in ocular hypertension, its value in patients diagnosed with glaucoma is less certain. There are several biological factors and genetic components that may influence glaucoma progression, which have been associated with thinner CCT. The CCT itself can be affected by several factors including ethnicity, age, sex, glaucoma medications, genetics, and the subtype of glaucoma. Besides, there is variability in the measurement of CCT between difference types of devices. These factors need to be considered in the evaluation of glaucoma patients' CCT and its effect on interpretation of intraocular pressure levels and risk stratification.

Genome-wide association study identifies WNT7B as a novel locus for central corneal thickness in Latinos

The cornea is the outermost layer of the eye and is a vital component of focusing incoming light on the retina. Central corneal thickness (CCT) is now recognized to have a significant role in ocular health and is a risk factor for various ocular diseases, such as keratoconus and primary open angle glaucoma. Most previous genetic studies utilized European and Asian subjects to identify genetic loci associated with CCT. Minority populations, such as Latinos, may aid in identifying additional loci and improve our understanding of the genetic architecture of CCT. In this study, we conducted a genome-wide association study (GWAS) in Latinos, a traditionally understudied population in genetic research, to further identify loci contributing to CCT. Study participants were genotyped using either the Illumina OmniExpress BeadChip (∼730K markers) or the Illumina Hispanic/SOL BeadChip (∼2.5 million markers). All study participants were 40 years of age and older. We assessed the association between individual single nucleotide polymorphisms (SNPs) and CCT using linear regression, adjusting for age, gender and principal components of genetic ancestry. To expand genomic coverage and to interrogate additional SNPs, we imputed SNPs from the 1000 Genomes Project reference panels. We identified a novel SNP, rs10453441 (P = 6.01E-09), in an intron of WNT7B that is associated with CCT. Furthermore, WNT7B is expressed in the human cornea. We also replicated 11 previously reported loci, including IBTK, RXRA-COL5A1, COL5A1, FOXO1, LRRK1 and ZNF469 (P < 1.25E-3). These findings provide further insight into the genetic architecture of CCT and illustrate that the use of minority groups in GWAS will help identify additional loci.

Genomic locus modulating corneal thickness in the mouse identifies POU6F2 as a potential risk of developing glaucoma

Central corneal thickness (CCT) is one of the most heritable ocular traits and it is also a phenotypic risk factor for primary open angle glaucoma (POAG). The present study uses the BXD Recombinant Inbred (RI) strains to identify novel quantitative trait loci (QTLs) modulating CCT in the mouse with the potential of identifying a molecular link between CCT and risk of developing POAG. The BXD RI strain set was used to define mammalian genomic loci modulating CCT, with a total of 818 corneas measured from 61 BXD RI strains (between 60–100 days of age). The mice were anesthetized and the eyes were positioned in front of the lens of the Phoenix Micron IV Image-Guided OCT system or the Bioptigen OCT system. CCT data for each strain was averaged and used to QTLs modulating this phenotype using the bioinformatics tools on GeneNetwork (www.genenetwork.org). The candidate genes and genomic loci identified in the mouse were then directly compared with the summary data from a human POAG genome wide association study (NEIGHBORHOOD) to determine if any genomic elements modulating mouse CCT are also risk factors for POAG.This analysis revealed one significant QTL on Chr 13 and a suggestive QTL on Chr 7. The significant locus on Chr 13 (13 to 19 Mb) was examined further to define candidate genes modulating this eye phenotype. For the Chr 13 QTL in the mouse, only one gene in the region (Pou6f2) contained nonsynonymous SNPs. Of these five nonsynonymous SNPs in Pou6f2, two resulted in changes in the amino acid proline which could result in altered secondary structure affecting protein function. The 7 Mb region under the mouse Chr 13 peak distributes over 2 chromosomes in the human: Chr 1 and Chr 7. These genomic loci were examined in the NEIGHBORHOOD database to determine if they are potential risk factors for human glaucoma identified using meta-data from human GWAS. The top 50 hits all resided within one gene (POU6F2), with the highest significance level of p = 10⁻⁶ for SNP rs76319873. POU6F2 is found in retinal ganglion cells and in corneal limbal stem cells. To test the effect of POU6F2 on CCT we examined the corneas of a Pou6f2-null mice and the corneas were thinner than those of wild-type littermates. In addition, these POU6F2 RGCs die early in the DBA/2J model of glaucoma than most RGCs. Using a mouse genetic reference panel, we identified a transcription factor, Pou6f2, that modulates CCT in the mouse. POU6F2 is also found in a subset of retinal ganglion cells and these RGCs are sensitive to injury.

Glaucoma is a complex group of diseases with several known causal mutations and many known risk factors. One well-known risk factor for developing primary open angle glaucoma is the thickness of the central cornea. The present study leverages a unique blend of systems biology methods using BXD recombinant inbred mice and genome-wide association studies from humans to define a putative molecular link between a phenotypic risk factor (central corneal thickness) and glaucoma. We identified a transcription factor, POU6F2, that is found in the developing retinal ganglion cells and cornea. POU6F2 is also present in a subpopulation of retinal ganglion cells and in stem cells of the cornea. Functional studies reveal that POU6F2 is associated with the central corneal thickness and susceptibility of retinal ganglion cells to injury.

Evaluation of Primary Angle-Closure Glaucoma Susceptibility Loci in Patients with Early Stages of Angle-Closure Disease

PURPOSE

To investigate whether newly identified genetic loci for primary angle-closure glaucoma (PACG) are associated with early stage angle-closure disease defined as primary angle closure suspect (PACS).

DESIGN

Case-control study.

PARTICIPANTS

A total of 1397 PACS patients and 943 controls of Chinese ethnicity from Singapore and 604 PACS patients and 287 controls of Indian ethnicity.

METHODS

The 8 PACG single nucleotide polymorphisms (SNPs; rs11024102 at PLEKHA7, rs3753841 at COL11A1, rs1015213 located between PCMTD1 and ST18 son chromosome 8q, rs3816415 at EPDR1, rs1258267 at CHAT, rs736893 at GLIS3, rs7494379 at FERMT2, and rs3739821 mapping in between DPM2 and FAM102A) were genotyped by Taqman assays. The association between SNP genotypes and PACS status was measured using logistic regression. A P value of 0.006 was set to account for the testing of 8 genetic loci using a Bonferroni correction. A meta-analysis was conducted to calculate the overall P value and accompanying per-allele odds ratios for each SNP analyzed.

MAIN OUTCOME MEASURES

Association of PACG loci with PACS status.

RESULTS

The PACS patients were significantly older in both cohorts (Chinese, P < 0.001; Indian, P = 0.002), and there were also more women (P < 0.001, both Chinese and Indian cohorts). In the Chinese cohort, significant evidence of association was noted at 3 SNPs: rs1015213 [A] in PCMTD1-ST18 (odds ratio [OR], 2.36; 95% confidence interval [CI], 1.36-4.11; P = 0.002), rs3816415 [A] in EPDR1 (OR, 1.49; 95% CI, 1.19-1.85; P < 0.001), and rs3739821 [G] in DPM2-FAM102A (OR, 1.40; 95% CI, 1.18-1.65; P < 0.001). Only PCMTD1-ST-18 was replicated modestly in the Indian population (P = 0.056). Meta-analysis showed significant evidence of association for PCMTD1-ST-18 (OR, 1.55; 95% CI, 1.18-2.04; P = 0.002) and DPM2-FAM102A (OR, 1.27; 95% CI, 1.12-1.45; P = 0.0002).

CONCLUSIONS

In this study, 2 of 8 PACG-associated loci were associated significantly with PACS status, the earliest stage in the angle-closure glaucoma disease course. The association of these PACG loci with PACS status suggests that these loci may confer susceptibility to a narrow angle configuration.

A large multi-ethnic genome-wide association study identifies novel genetic loci for intraocular pressure

Elevated intraocular pressure (IOP) is a major risk factor for glaucoma, a leading cause of blindness. IOP heritability has been estimated to up to 67%, and to date only 11 IOP loci have been reported, accounting for 1.5% of IOP variability. Here, we conduct a genome-wide association study of IOP in 69,756 untreated individuals of European, Latino, Asian, and African ancestry. Multiple longitudinal IOP measurements were collected through electronic health records and, in total, 356,987 measurements were included. We identify 47 genome-wide significant IOP-associated loci (P < 5 × 10⁻⁸); of the 40 novel loci, 14 replicate at Bonferroni significance in an external genome-wide association study analysis of 37,930 individuals of European and Asian descent. We further examine their effect on the risk of glaucoma within our discovery sample. Using longitudinal IOP measurements from electronic health records improves our power to identify new variants, which together explain 3.7% of IOP variation.

Intraocular pressure (IOP) is a major risk factor for glaucoma. Here, Choquet and co-authors perform a multi-ethnic genome-wide association study of repeat IOP measurements in 69,756 individuals and identify 40 novel loci, 36 of which show directionally consistent effects in glaucoma.

Novel Zinc Finger Protein Gene 469 (ZNF469) Variants in Advanced Keratoconus

PURPOSE

Common polymorphic variants upstream of Zinc finger protein gene 469 (ZNF469) have been associated with central corneal thickness. Rare ZNF469 variants have been shown in keratoconus patients. The aim of the current study was to investigate the frequency of ZNF 469 gene variants in rapidly progressive advance keratoconus patients who underwent corneal transplant surgery by the age of 30, compared to their frequency in the normal Turkish population.

METHODS

A search in a patient database was performed to identify patients with a rapidly progressive keratoconus requiring corneal transplant surgery by the age of 30 in at least one eye. Twenty-six advance keratoconus patients (study group) and 109 health subjects (control group) were included in the study. Blood samples were donated, and genomic DNA was extracted. The entire coding sequence of the ZNF469 gene including the 84 bp of the putative intron was amplified using PCR primers and analyzed using next generation sequencing (NGS).

RESULTS

Fifteen single nucleotide polymorphisms previously reported and registered to the dbSNP database were detected in the study group. The allele frequencies of these polymorphisms were higher in the keratoconus group compared to the control group and to the ExAC genome database. Three new missense heterozygote variants and one new synonym variant were detected in keratoconus group. According to prediction software, the P873T and Q2188H variants were shown to be non-tolerated, whereas G3424S could be tolerated. The synonymous variant R1060R is not predicted to lead to abnormal splicing by Human Splicing Finder in silico analysis.

CONCLUSION

New detected ZNF 469 P873T and Q2188H heterozygote coding variants in isolated advance keratoconus patients may be associated with the disease pathogenesis.

Characterizing the "POAGome": A bioinformatics-driven approach to primary open-angle glaucoma

Primary open-angle glaucoma (POAG) is a genetically, physiologically, and phenotypically complex neurodegenerative disorder. This study addressed the expanding collection of genes associated with POAG, referred to as the "POAGome." We used bioinformatics tools to perform an extensive, systematic literature search and compiled 542 genes with confirmed associations with POAG and its related phenotypes (normal tension glaucoma, ocular hypertension, juvenile open-angle glaucoma, and primary congenital glaucoma). The genes were classified according to their associated ocular tissues and phenotypes, and functional annotation and pathway analyses were subsequently performed. Our study reveals that no single molecular pathway can encompass the pathophysiology of POAG. The analyses suggested that inflammation and senescence may play pivotal roles in both the development and perpetuation of the retinal ganglion cell degeneration seen in POAG. The TGF-β signaling pathway was repeatedly implicated in our analyses, suggesting that it may be an important contributor to the manifestation of POAG in the anterior and posterior segments of the globe. We propose a molecular model of POAG revolving around TGF-β signaling, which incorporates the roles of inflammation and senescence in this disease. Finally, we highlight emerging molecular therapies that show promise for treating POAG.

Central corneal thickness in glaucoma

PURPOSE OF REVIEW

The purpose of this review is to summarize the role of central corneal thickness (CCT) in the clinical management of a glaucoma patient.

RECENT FINDINGS

The prognostic value of CCT is well recognized in patients with ocular hypertension. However, its predictive value in other glaucoma suspects and patients with established glaucoma is less certain. Tonometry artefacts can result from variations in CCT. However, an adequately validated correction algorithm for Goldmann applanation tonometry measurements does not exist. Newer methods of tonometry are potentially less influenced by CCT but are limited in their clinical use. There may also be biological and genetic associations between corneal thickness and glaucoma. Demographics, environmental factors, glaucoma treatment and the measurement device used have a significant influence on CCT, and should be considered when interpreting the effect of cornea thickness in patients with glaucoma. New measurements of the biomechanical properties of the cornea are likely to be better approximations of the globe biomechanics than CCT, but these require further evaluation.

SUMMARY

The clinical significance of CCT is well recognized in the context of glaucoma diagnosis and management, though the extent of its importance remains debatable. Corneal biomechanical properties may be more significantly associated with glaucoma than CCT.

Genetics in Keratoconus: where are we?

Keratoconus (KC) is a non-inflammatory thinning and protrusion of the cornea in which the cornea assumes a conical shape. Complex etiology of this condition at present remains an enigma. Although environmental factors have been involved in KC pathogenesis, strong underlining genetic susceptibility has been proven. The lack of consistent findings among early genetic studies suggested a heterogeneity and complex nature of the genetic contribution to the development of KC. Recently, genome-wide linkage studies (GWLS) and genome-wide association studies (GWAS) were undertaken. Next-generation sequencing (NGS)-based genomic screens are also currently being carried out. Application of these recently developed comprehensive genetic tools led to a much greater success and increased reproducibility of genetic findings in KC. Involvement of the LOX gene identified through GWLS has been confirmed in multiple cohorts of KC patients around the world. KC susceptibility region located at the 2q21.3 chromosomal region near the RAB3GAP1 gene identified through GWAS was independently replicated. Rare variants in the ZNF469 gene (mutated in corneal dystrophy Brittle Cornea Syndrome) and in the TGFBI gene (mutated in multiple corneal epithelial–stromal TGFBI dystrophies) have been repeatedly identified in familial and sporadic KC patients of different ethnicities. Additional comprehensive strategies using quantitative endophenotypes have been successfully employed to bring further understanding to the genetics of KC. Additional genetic determinants including the COL5A1 gene have been identified in the GWAS of KC-related trait central corneal thickness. These recent discoveries confirmed the importance of the endophenotype approach for studying complex genetic diseases such as KC and showed that different connective tissue disorders may have the same genetic determinants.

Stereo Photo Measured ONH Shape Predicts Development of POAG in Subjects With Ocular Hypertension

PURPOSE

To identify objective, quantitative optic nerve head (ONH) structural features and model the contributions of glaucoma.

METHODS

Baseline stereoscopic optic disc images of 1635 glaucoma-free participants at risk for developing primary open-angle glaucoma (POAG) were collected as part of the Ocular Hypertension Treatment Study. A stereo correspondence algorithm designed for fundus images was applied to extract the three-dimensional (3D) information about the ONH. Principal component analysis was used to identify ONH 3D structural features and the contributions of demographic features, clinical variables, and disease were modeled using linear regression and linear component analysis. The computationally identified features were evaluated based on associations with glaucoma and ability to predict which participants would develop POAG.

RESULTS

The computationally identified features were significantly associated with future POAG, POAG-related demographics (age, ethnicity), and clinical measurements (horizontal and vertical cup-to-disc ratio, central corneal thickness, and refraction). Models predicting future POAG development using the OHTS baseline data and STEP features achieved an AUC of 0.722 in cross-validation testing. This was a significant improvement over using only demographics (age, sex, and ethnicity), which had an AUC of 0.599.

CONCLUSIONS

Methods for identifying objective, quantitative measurements of 3D ONH structure were developed using a large dataset. The identified features were significantly associated with POAG and POAG-related variables. Further, these features increased predictive model accuracy in predicting future POAG. The results indicate that the computationally identified features might be useful in POAG early screening programs or as endophenotypes to investigate POAG genetics.

Genetic Evidence for Differential Regulation of Corneal Epithelial and Stromal Thickness

PURPOSE

Central corneal thickness (CCT) is a quantitative trait associated with keratoconus and primary open-angle glaucoma. Although CCT is highly heritable, known genetic variations explain only a fraction of the phenotypic variability. The purpose of this study was to identify additional CCT-influencing loci using inbred strains of mice.

METHODS

Cohorts of 82 backcrossed (N2) and 99 intercrossed (F2) mice were generated from crosses between recombinant inbred BXD24/TyJ and wild-derived CAST/EiJ mice. Using anterior chamber optical coherence tomography, mice were phenotyped at 10 to 12 weeks of age, genotyped based on 96 genome-wide single nucleotide polymorphisms (SNPs), and subjected to quantitative trait locus (QTL) analysis.

RESULTS

In an analysis of total CCT among all mice, two loci passed the significance threshold of P = 0.05. These were on Chr 3 and Chr 11 (Cctq4 and Cctq5, respectively). A third locus of interest was identified in a two-dimensional pairwise analysis; this locus on Chr 14 (Cctq6) exhibited a significant additive effect with Cctq5. Independent analyses of the dataset for epithelial and stromal thickness revealed that Cctq4 is specific to the epithelial layer and that Cctq5 and Cctq6 are specific to the stromal layer.

CONCLUSIONS

Our findings demonstrate a quantitative multigenic pattern of CCT inheritance in mice and identify three previously unrecognized CCT-influencing loci: Cctq4, Cctq5, and Cctq6. This is the first demonstration that distinct layers of the cornea are under differential genetic control and highlights the need to refine the design of future genome-wide association studies of CCT.

Genetic association of SNPs near ATOH7, CARD10, CDKN2B, CDC7 and SIX1/SIX6 with the endophenotypes of primary open angle glaucoma in Indian population

Primary open angle glaucoma (POAG) belonging to a group of optic neuropathies, result from interaction between genetic and environmental factors. Study of associations with quantitative traits (QTs) is one of the successful strategies to understand the complex genetics of POAG. The current study attempts to explore the association of variations near/in genes like ATOH7, SIX1/SIX6 complex, CDKN2B, CARD10, and CDC7 with POAG and its QTs including vertical cup to disc ratio (VCDR), central corneal thickness (CCT), intra ocular pressure (IOP), and axial length (AL). Case-control study design was carried out in a sample size of 97 POAG cases and 371 controls from South India. Model-based (additive, recessive, dominant) association of the genotypes and their interaction was carried out between cases and controls using chi-square, linear and logistic regression methods. Nominal significance (P<0.05) was observed for QTs like i) VCDR with SNPs rs1900004 (ATOH7); rs1192415 (CDC7); rs10483727 (SIX1/SIX6), rs9607469 (CARD10); ii) CCT with rs1192415; iii) IOP with rs1900004 and iv) AL with rs1900004 and rs1063192 (CDKN2B). We were able to replicate previously known interactions between ATOH7-SIX6 and SIX6-CDKN2B along with few novel interactions between ATOH7—CDC7 and SIX6 with genes including CARD10 and CDC7. In summary, our results suggest that a probable interaction among the candidate genes for QTs, play a major role in determining the individual’s susceptibility to POAG.

Evaluating the Association between Keratoconus and Reported Genetic Loci in a Han Chinese Population

BACKGROUND

Keratoconus (KC) is a complex degenerative disorder of the cornea. Genetic, environmental, and lifestyle factors may all contribute to the pathogenesis of KC. Most of the reported KC-associated SNPs have been detected in Caucasians and Australians. To investigate whether the reported associated SNPs can be found in a Chinese population, we performed a replication study of the significantly associated SNPs.

MATERIALS AND METHODS

A total of 210 unrelated Chinese KC patients and 191 unrelated controls were included in the present study. SNPs rs4954218 (Near RAB3GAP1 (5')), rs4894535 (FNDC3B), rs2956540 (LOX), rs3735520 (Near HGF (5')), rs1324183 (MPDZ-NF1B), rs1536482 (RXRA-COL5A1), rs7044529 (COL5A1), rs2721051 (Near FOXO1 (3')), rs9938149 (BANP-ZNF469) and rs6050307 (VSX1) were assessed for their association with KC. The genotype of each SNP was detected using the Sequenom MassARRAY-Assay.

RESULTS

SNP rs1324183 located in MPDZ-NF1B was associated with an increased risk of KC (OR=3.108, 95% CI=1.366-7.072, p=0.005), and SNP rs2956540 in the LOX gene may confer a reduced risk of KC with a borderline p value in our population (OR=0.664, 95% CI=0.447-0.986, p=0.042). No significant difference was observed between patients and controls in the other eight SNP genotypes and allele frequencies.

CONCLUSIONS

The replication association of rs1324183 (MPDZ-NF1B) with KC in our population and the results, which are identical to those in different populations, suggest that rs1324183 (MPDZ-NF1B) is a common genetic risk for KC and should be further investigated.

Genome-wide association study for refractive astigmatism reveals genetic co-determination with spherical equivalent refractive error: the CREAM consortium

To identify genetic variants associated with refractive astigmatism in the general population, meta-analyses of genome-wide association studies were performed for: White Europeans aged at least 25 years (20 cohorts, N = 31,968); Asian subjects aged at least 25 years (7 cohorts, N = 9,295); White Europeans aged <25 years (4 cohorts, N = 5,640); and all independent individuals from the above three samples combined with a sample of Chinese subjects aged <25 years (N = 45,931). Participants were classified as cases with refractive astigmatism if the average cylinder power in their two eyes was at least 1.00 diopter and as controls otherwise. Genome-wide association analysis was carried out for each cohort separately using logistic regression. Meta-analysis was conducted using a fixed effects model. In the older European group the most strongly associated marker was downstream of the neurexin-1 (NRXN1) gene (rs1401327, P = 3.92E-8). No other region reached genome-wide significance, and association signals were lower for the younger European group and Asian group. In the meta-analysis of all cohorts, no marker reached genome-wide significance: The most strongly associated regions were, NRXN1 (rs1401327, P = 2.93E-07), TOX (rs7823467, P = 3.47E-07) and LINC00340 (rs12212674, P = 1.49E-06). For 34 markers identified in prior GWAS for spherical equivalent refractive error, the beta coefficients for genotype versus spherical equivalent, and genotype versus refractive astigmatism, were highly correlated (r = -0.59, P = 2.10E-04). This work revealed no consistent or strong genetic signals for refractive astigmatism; however, the TOX gene region previously identified in GWAS for spherical equivalent refractive error was the second most strongly associated region. Analysis of additional markers provided evidence supporting widespread genetic co-susceptibility for spherical and astigmatic refractive errors.

Linkage disequilibrium and signatures of positive selection around LINE-1 retrotransposons in the human genome

Insertions of the human-specific subfamily of LINE-1 (L1) retrotransposon are highly polymorphic across individuals and can critically influence the human transcriptome. We hypothesized that L1 insertions could represent genetic variants determining important human phenotypic traits, and performed an integrated analysis of L1 elements and single nucleotide polymorphisms (SNPs) in several human populations. We found that a large fraction of L1s were in high linkage disequilibrium with their surrounding genomic regions and that they were well tagged by SNPs. However, L1 variants were only partially captured by SNPs on standard SNP arrays, so that their potential phenotypic impact would be frequently missed by SNP array-based genome-wide association studies. We next identified potential phenotypic effects of L1s by looking for signatures of natural selection linked to L1 insertions; significant extended haplotype homozygosity was detected around several L1 insertions. This finding suggests that some of these L1 insertions may have been the target of recent positive selection.

ABCC5, a gene that influences the anterior chamber depth, is associated with primary angle closure glaucoma

Anterior chamber depth (ACD) is a key anatomical risk factor for primary angle closure glaucoma (PACG). We conducted a genome-wide association study (GWAS) on ACD to discover novel genes for PACG on a total of 5,308 population-based individuals of Asian descent. Genome-wide significant association was observed at a sequence variant within ABCC5 (rs1401999; per-allele effect size = −0.045 mm, P = 8.17×10⁻⁹). This locus was associated with an increase in risk of PACG in a separate case-control study of 4,276 PACG cases and 18,801 controls (per-allele OR = 1.13 [95% CI: 1.06–1.22], P = 0.00046). The association was strengthened when a sub-group of controls with open angles were included in the analysis (per-allele OR = 1.30, P = 7.45×10⁻⁹; 3,458 cases vs. 3,831 controls). Our findings suggest that the increase in PACG risk could in part be mediated by genetic sequence variants influencing anterior chamber dimensions.

The anterior chamber is the space within the eye which is bound by the cornea, and the anterior surfaces of the iris and lens. Anterior chamber depth (ACD) is the distance measured along the eye's optical axis, from the cornea to the lens surface. ACD is an important risk factor for primary angle closure glaucoma (PACG), a major cause of irreversible blindness worldwide, and in particular, individuals of Asian ethnicity. In order to identify the genes that underlie PACG susceptibility, we conducted a two-staged study. We first conducted a large scale genetic study on a total of 5,308 population-based individuals of Asian descent to identify the genetic variants that influence ACD. This was followed by testing for associations between the identified genetic variant and PACG in another independent collection of 4,276 PACG cases and 18,801 controls. We found that a genetic variant within ABCC5 was associated with an increased risk of having PACG. Our findings suggest that the increase in PACG risk could in part be mediated by genetic sequence variants that influence the anterior chamber dimensions of the eye.

Combined genotype and haplotype tests for region-based association studies

Background

Although single-SNP analysis has proven to be useful in identifying many disease-associated loci, region-based analysis has several advantages. Empirically, it has been shown that region-based genotype and haplotype approaches may possess much higher power than single-SNP statistical tests. Both high quality haplotypes and genotypes may be available for analysis given the development of next generation sequencing technologies and haplotype assembly algorithms.

Results

As generally it is unknown whether genotypes or haplotypes are more relevant for identifying an association, we propose to use both of them with the purpose of preserving high power under both genotype and haplotype disease scenarios. We suggest two approaches for a combined association test and investigate the performance of these two approaches based on a theoretical model, population genetics simulations and analysis of a real data set.

Conclusions

Based on a theoretical model, population genetics simulations and analysis of a central corneal thickness (CCT) Genome Wide Association Study (GWAS) data set we have shown that combined genotype and haplotype approach has a high potential utility for applications in association studies.

Large-scale genotyping identifies a new locus at 22q13.2 associated with female breast size

BACKGROUND

Individual differences in breast size are a conspicuous feature of variation in human females and have been associated with fecundity and advantage in selection of mates. To identify common variants that are associated with breast size, we conducted a large-scale genotyping association meta-analysis in 7169 women of European descent across three independent sample collections with digital or screen film mammograms.

METHODS

The samples consisted of the Swedish KARMA, LIBRO-1 and SASBAC studies genotyped on iCOGS, a custom illumina iSelect genotyping array comprising of 211 155 single nucleotide polymorphisms (SNPs) designed for replication and fine mapping of common and rare variants with relevance to breast, ovary and prostate cancer. Breast size of each subject was ascertained by measuring total breast area (mm(2)) on a mammogram.

RESULTS

We confirm genome-wide significant associations at 8p11.23 (rs10086016, p=1.3×10(-14)) and report a new locus at 22q13 (rs5995871, p=3.2×10(-8)). The latter region contains the MKL1 gene, which has been shown to impact endogenous oestrogen receptor α transcriptional activity and is recruited on oestradiol sensitive genes. We also replicated previous genome-wide association study findings for breast size at four other loci.

CONCLUSIONS

A new locus at 22q13 may be associated with female breast size.

Developments in Ocular Genetics: Annual Review

PURPOSE

The purpose of this study was to summarize major developments in ocular genetics over the past year.

DESIGN

A literature review was performed for articles relating to the genetics of eye diseases and morphology. The search focused on articles published between September 15, 2011, and September 15, 2012.

METHODS

PubMed and Google Scholar search tools were used to search for ocular genetics articles in the desired date range.

RESULTS

Major advances have been reported in numerous areas including glaucoma, age-related macular degeneration, and keratoconus. Numerous novel associations have been identified through large genome-wide association studies. In addition, numerous disease genes have been identified through next-generation sequencing technologies.

CONCLUSIONS

Ocular genetics continues to advance at a rapid pace and benefit from new technologies. Numerous discoveries in the past year point toward areas for continued research.

A genome-wide association study of central corneal thickness in Latinos

PURPOSE

Central corneal thickness (CCT) is a clinically important risk factor for primary open-angle glaucoma and keratoconus. Genetic factors controlling CCT in Latinos, the most populous minority population in the United States, are unclear. Here we describe the first genome-wide association study (GWAS) report of CCT in Latinos.

METHODS

We performed a GWAS for CCT on 1768 Latinos recruited in the Los Angeles Latino Eye Study (LALES) using Illumina's HumanOmniExpress BeadChip (∼730K markers). To discover additional associated single-nucleotide polymorphisms (SNPs), we imputed SNPs based on the 1000 Genomes Project reference panels. All subjects were 40 years of age and older. We used linear regression with adjustment for age, sex, and principal components of genetic ancestry.

RESULTS

we replicated the involvement of several previously reported loci, SUCH AS RXRA-COL5A1, FOXO1, and ZNF469, for CCT in Latinos (P 0.002). moreover, we discovered novel SNPS, RS3118515, RS943423, RS3118594, AND RS3132307, THAT REACHED GWAS SIGNIFICANCE (P 5 10(8)) in the uncharacterized LOC100506532 (GENE TYPE: miscRNA) for CCT in Latinos. By conditional analysis, we demonstrate that rs3118515 in this gene is responsible for the GWAS signal in the chromosome 9 RXRA-COL5A1 region in Latinos. Moreover, multiple sources of ENCODE evidence suggest that rs3118515 is in a regulatory region. Reverse-transcription PCR products indicated that transcripts of LOC100506532 surrounding rs3118515 were expressed in human corneas.

CONCLUSIONS

We discovered novel SNPs for CCT in Latinos and provided the first reported evidence of the corneal expression of LOC100506532. These results help to further increase our understanding of the genetic architecture of CCT.

Genome-wide association analyses identify multiple loci associated with central corneal thickness and keratoconus

Central corneal thickness (CCT) is associated with eye conditions including keratoconus and glaucoma. We performed a meta-analysis on >20,000 individuals in European and Asian populations that identified 16 new loci associated with CCT at genome-wide significance (P < 5 × 10(-8)). We further showed that 2 CCT-associated loci, FOXO1 and FNDC3B, conferred relatively large risks for keratoconus in 2 cohorts with 874 cases and 6,085 controls (rs2721051 near FOXO1 had odds ratio (OR) = 1.62, 95% confidence interval (CI) = 1.4-1.88, P = 2.7 × 10(-10), and rs4894535 in FNDC3B had OR = 1.47, 95% CI = 1.29-1.68, P = 4.9 × 10(-9)). FNDC3B was also associated with primary open-angle glaucoma (P = 5.6 × 10(-4); tested in 3 cohorts with 2,979 cases and 7,399 controls). Further analyses implicate the collagen and extracellular matrix pathways in the regulation of CCT.

Population-based meta-analysis in Caucasians confirms association with COL5A1 and ZNF469 but not COL8A2 with central corneal thickness

Central corneal thickness (CCT) has become an endophenotype of major interest for the genetically complex disorder glaucoma. CCT has a high heritability, and thin CCT is an independent risk factor for the diagnosis and progression of open-angle glaucoma. Genome-wide association studies thus provide genetic loci associated with CCT and potentially related to open-angle glaucoma. The distribution of CCT and prevalence of glaucoma in population-based studies have demonstrated ethnic differences suggesting ethnic-dependent variations in the genetic determinants of CCT. We conducted a genome-wide association study in Caucasians (n = 3,931) from the Gutenberg Health Study (Germany) followed by replication of 30 genome-wide significant SNPs or SNPs of interest (P < 10(-5)) in the Rotterdam Study (The Netherlands, n = 1,418). In a combined analysis, we confirmed quantitative trait loci on chromosomes 9q34 and 16q24 for association with CCT. On chromosome 16q24, the locus is located in an intergenic region near the ZNF469 gene (top SNP: rs9938149, P = 1.45 × 10(-12)). ZNF469 missense mutation is involved in a syndrome with very thin cornea (brittle cornea syndrome). The second locus on chromosome 9q34 represents the intergenic region between the RXRA and COL5A1 gene (top SNP: rs3132306, P = 2.71 × 10(-10)). Collagen type 5 determines the diameter of the corneal collagen fibrils. In our Caucasian population-based GWA study, we reinforce the involvement of collagen-related genes influencing CCT in Caucasians. We could not confirm the collagen type 8 locus on chromosome 1 as reported in Asian studies.

Genome-wide analysis of central corneal thickness in primary open-angle glaucoma cases in the NEIGHBOR and GLAUGEN consortia

PURPOSE

To investigate the effects of central corneal thickness (CCT)-associated variants on primary open-angle glaucoma (POAG) risk using single nucleotide polymorphisms (SNP) data from the Glaucoma Genes and Environment (GLAUGEN) and National Eye Institute (NEI) Glaucoma Human Genetics Collaboration (NEIGHBOR) consortia.

METHODS

A replication analysis of previously reported CCT SNPs was performed in a CCT dataset (n = 1117) and these SNPs were then tested for association with POAG using a larger POAG dataset (n = 6470). Then a CCT genome-wide association study (GWAS) was performed. Top SNPs from this analysis were selected and tested for association with POAG. cDNA libraries from fetal and adult brain and ocular tissue samples were generated and used for candidate gene expression analysis.

RESULTS

Association with one of 20 previously published CCT SNPs was replicated: rs12447690, near the ZNF469 gene (P = 0.001; β = -5.08 μm/allele). None of these SNPs were significantly associated with POAG. In the CCT GWAS, no SNPs reached genome-wide significance. After testing 50 candidate SNPs for association with POAG, one SNP was identified, rs7481514 within the neurotrimin (NTM) gene, that was significantly associated with POAG in a low-tension subset (P = 0.00099; Odds Ratio [OR] = 1.28). Additionally, SNPs in the CNTNAP4 gene showed suggestive association with POAG (top SNP = rs1428758; P = 0.018; OR = 0.84). NTM and CNTNAP4 were shown to be expressed in ocular tissues.

CONCLUSIONS

The results suggest previously reported CCT loci are not significantly associated with POAG susceptibility. By performing a quantitative analysis of CCT and a subsequent analysis of POAG, SNPs in two cell adhesion molecules, NTM and CNTNAP4, were identified and may increase POAG susceptibility in a subset of cases.

Genome-wide association studies in the hunt for genes causing primary open-angle glaucoma: a review

Genome-wide association studies are a powerful tool for the identification of genetic risk factors for complex disease. This methodology has been successfully applied to primary open-angle glaucoma through the analysis of primary open-angle glaucoma (POAG) as well as specific subgroups of patients including those with normal tension glaucoma and advanced glaucoma. In addition, the analysis of quantitative traits important in POAG, including optic disc area and vertical cup-to-disc ratio has also identified genes important in POAG development. This review explores findings of genome-wide association studies for POAG and related traits.

The Genetics of Keratoconus: A Review

Keratoconus is the most common ectatic disorder of the corneal. Genetic and environmental factors may contribute to its pathogenesis. The focus of this article is to summarize current research into the complex genetics of keratoconus. We discuss the evidence of genetic etiology including family-based linkage studies, twin studies, genetic mutations, and genome-wide association studies. The genes implicated potentially include VSX1, miR-184, DOCK9, SOD1, RAB3GAP1, and HGF. Besides the coding mutations, we also highlight the potential contribution of DNA copy number variants in the pathogenesis of keratoconus. Finally, we present future directions for genetic research in the understanding of the complex genetics of keratoconus and its clinical significance. As new functional, candidate genes for keratoconus are being discovered at a rapid pace, the molecular genetic mechanisms underlying keratoconus pathogenesis will advance our understanding of keratoconus and promote the development of a novel therapy.