Sonic Hedgehog Intron Variant Associated With an Unusual Pediatric Cortical Cataract.

Purpose: To identify the genetic basis of an unusual pediatric cortical cataract demonstrating autosomal dominant inheritance in a large European-Australian pedigree. Methods: DNA from four affected individuals were exome sequenced utilizing a NimbleGen SeqCap EZ Exome V3 kit and HiSeq 2500. DNA from 12 affected and four unaffected individuals were genotyped using Human OmniExpress-24 BeadChips. Multipoint linkage and haplotyping were performed (Superlink-Online SNP). DNA from one affected individual and his unaffected father were whole-genome sequenced on a HiSeq X Ten system. Rare small insertions/deletions and single-nucleotide variants (SNVs) were identified in the disease-linked region (Golden Helix SVS). Combined Annotation Dependent Depletion (CADD) analysis predicted variant deleteriousness. Putative enhancer function and variant effects were determined using the Dual-Glo Luciferase Assay system. Results: Linkage mapping identified a 6.23-centimorgan support interval at chromosome 7q36. A co-segregating haplotype refined the critical region to 6.03 Mbp containing 21 protein-coding genes. Whole-genome sequencing uncovered 114 noncoding variants from which CADD predicted one was highly deleterious, a novel substitution within intron-1 of the sonic hedgehog signaling molecule (SHH) gene. ENCODE data suggested this site was a putative enhancer, subsequently confirmed by luciferase reporter assays with variant-associated gene overexpression. Conclusions: In a large pedigree, we have identified a SHH intron variant that co-segregates with an unusual pediatric cortical cataract phenotype. SHH is important for lens formation, and mutations in its receptor (PTCH1) cause syndromic cataract. Our data implicate increased function of an enhancer important for SHH expression primarily within developing eye tissues.

Cellular crosstalk regulates the aqueous humor outflow pathway and provides new targets for glaucoma therapies.

Primary congenital glaucoma (PCG) is a severe disease characterized by developmental defects in the trabecular meshwork (TM) and Schlemm's canal (SC), comprising the conventional aqueous humor outflow pathway of the eye. Recently, heterozygous loss of function variants in TEK and ANGPT1 or compound variants in TEK/SVEP1 were identified in children with PCG. Moreover, common variants in ANGPT1and SVEP1 have been identified as risk alleles for primary open angle glaucoma (POAG) in GWAS studies. Here, we show tissue-specific deletion of Angpt1 or Svep1 from the TM causes PCG in mice with severe defects in the adjacent SC. Single-cell transcriptomic analysis of normal and glaucomatous Angpt1 deficient eyes allowed us to identify distinct TM and SC cell populations and discover additional TM-SC signaling pathways. Furthermore, confirming the importance of angiopoietin signaling in SC, delivery of a recombinant ANGPT1-mimetic promotes developmental SC expansion in healthy and Angpt1 deficient eyes, blunts intraocular pressure (IOP) elevation and RGC loss in a mouse model of PCG and lowers IOP in healthy adult mice. Our data highlight the central role of ANGPT1-TEK signaling and TM-SC crosstalk in IOP homeostasis and provide new candidates for SC-targeted glaucoma therapy.

SVEP1 as a Genetic Modifier of TEK-Related Primary Congenital Glaucoma.

Purpose: Affecting children by age 3, primary congenital glaucoma (PCG) can cause debilitating vision loss by the developmental impairment of aqueous drainage resulting in high intraocular pressure (IOP), globe enlargement, and optic neuropathy. TEK haploinsufficiency accounts for 5% of PCG in diverse populations, with low penetrance explained by variable dysgenesis of Schlemm's canal (SC) in mice. We report eight families with TEK-related PCG, and provide evidence for SVEP1 as a disease modifier in family 8 with a higher penetrance and severity. Methods: Exome sequencing identified coding/splice site variants with an allele frequency less than 0.0001 (gnomAD). TEK variant effects were assayed in construct-transfected HEK293 cells via detection of autophosphorylated (active) TEK protein. An enucleated eye from an affected member of family 8 was examined via histology. SVEP1 expression in developing outflow tissues was detected by immunofluorescent staining of 7-day mouse anterior segments. SVEP1 stimulation of TEK expression in human umbilical vascular endothelial cells (HUVECs) was measured by TaqMan quantitative PCR. Results: Heterozygous TEK loss-of-function alleles were identified in eight PCG families, with parent-child disease transmission observed in two pedigrees. Family 8 exhibited greater disease penetrance and severity, histology revealed absence of SC in one eye, and SVEP1:p.R997C was identified in four of the five affected individuals. During SC development, SVEP1 is secreted by surrounding tissues. SVEP1:p.R997C abrogates stimulation of TEK expression by HUVECs. Conclusions: We provide further evidence for PCG caused by TEK haploinsufficiency, affirm autosomal dominant inheritance in two pedigrees, and propose SVEP1 as a modifier of TEK expression during SC development, affecting disease penetrance and severity.

Correction: Genetic variants linked to myopic macular degeneration in persons with high myopia: CREAM Consortium.

[This corrects the article DOI: 10.1371/journal.pone.0220143.].

Genetic variants linked to myopic macular degeneration in persons with high myopia: CREAM Consortium.

PURPOSE: To evaluate the roles of known myopia-associated genetic variants for development of myopic macular degeneration (MMD) in individuals with high myopia (HM), using case-control studies from the Consortium of Refractive Error and Myopia (CREAM). METHODS: A candidate gene approach tested 50 myopia-associated loci for association with HM and MMD, using meta-analyses of case-control studies comprising subjects of European and Asian ancestry aged 30 to 80 years from 10 studies. Fifty loci with the strongest associations with myopia were chosen from a previous published GWAS study. Highly myopic (spherical equivalent [SE] ≤ -5.0 diopters [D]) cases with MMD (N = 348), and two sets of controls were enrolled: (1) the first set included 16,275 emmetropes (SE ≤ -0.5 D); and (2) second set included 898 highly myopic subjects (SE ≤ -5.0 D) without MMD. MMD was classified based on the International photographic classification for pathologic myopia (META-PM). RESULTS: In the first analysis, comprising highly myopic cases with MMD (N = 348) versus emmetropic controls without MMD (N = 16,275), two SNPs were significantly associated with high myopia in adults with HM and MMD: (1) rs10824518 (P = 6.20E-07) in KCNMA1, which is highly expressed in human retinal and scleral tissues; and (2) rs524952 (P = 2.32E-16) near GJD2. In the second analysis, comprising highly myopic cases with MMD (N = 348) versus highly myopic controls without MMD (N = 898), none of the SNPs studied reached Bonferroni-corrected significance. CONCLUSIONS: Of the 50 myopia-associated loci, we did not find any variant specifically associated with MMD, but the KCNMA1 and GJD2 loci were significantly associated with HM in highly myopic subjects with MMD, compared to emmetropes.

In-utero epigenetic factors are associated with early-onset myopia in young children.

OBJECTIVES: To assess whether epigenetic mechanisms affecting gene expression may be involved in the pathogenesis of early-onset myopia, we performed genome-wide DNA methylation analyses of umbilical cord tissues, and assessed any associations between CpG site-specific methylation and the development of the disorder when the children were 3 years old. METHODS: Genome-wide DNA methylation profiling of umbilical cord samples from 519 Singaporean infants involved in a prospective birth cohort 'Growing Up in Singapore Towards healthy Outcomes' (GUSTO) was performed using the Illumina Infinium HumanMethylation450K chip microarray. Multivariable logistic regression models were used to assess any associations between site-specific CpG methylation of umbilical cord tissue at birth and myopia risk in 3 year old children, adjusting for potential confounders. Gene expression of genes located near CpG sites that demonstrated statistically significant associations were measured in relevant ocular tissues using human and mouse fetal and adult eye samples. RESULTS: We identified statistically significant associations between DNA methylation levels at five CpG sites and early-onset myopia risk after correcting for multiple comparisons using a false discovery rate of 5%. Two statistically significant CpG sites were identified in intergenic regions: 8p23(p = 1.70×10-7) and 12q23.2(p = 2.53×10-7). The remaining 3 statistically significant CpG sites were identified within the following genes: FGB (4q28, p = 3.60×10-7), PQLC1 (18q23, p = 8.9×10-7) and KRT12 (17q21.2, p = 1.2×10-6). Both PQLC1 and KRT12 were found to be significantly expressed in fetal and adult cornea and sclera tissues in both human and mouse. CONCLUSIONS: We identified five CpG methylation sites that demonstrate a statistically significant association with increased risk of developing early-onset myopia. These findings suggest that variability in the neonatal cord epigenome may influence early-onset myopia risk in children. Further studies of the epigenetic influences on myopia risk in larger study populations, and the associations with adulthood myopia risk are warranted.

Comparative Intraoperative Anterior Segment OCT Findings in Pediatric Patients with and without Glaucoma.

PURPOSE: To compare the angle structures during surgery using anterior segment (AS) OCT in children with and without glaucoma. DESIGN: Prospective comparative study, with single masked grader evaluating AS OCT images. PARTICIPANTS: This study included a total of 44 eyes of 8 children with glaucoma and 14 children without glaucoma undergoing ophthalmic surgery with general anesthesia. The median age for the glaucoma group was 6 months (range, 0.5-108 months), and that for the nonglaucoma group was 42 months (range, 14-95 months). METHODS: An Envisu C2300 handheld OCT with a high-resolution AS lens was used to capture videos of temporal and nasal quadrant anterior chamber angles from each eye. Still images of the temporal and nasal angles were selected and then graded by a masked ophthalmologist examiner to determine whether Schlemm's canal (SC) was visible, if abnormal tissue was present in the angle, and whether intrascleral lumens were present. Clinical correlation with severity of glaucoma was performed after AS OCT images were graded. MAIN OUTCOME MEASURES: Presence or absence of SC on AS OCT in glaucomatous and nonglaucomatous eyes. RESULTS: Schlemm's canal was present in 7 of 13 glaucomatous eyes compared with 27 of 27 nonglaucomatous eyes. In patients with unilateral glaucoma, absence of SC was noted only in the affected eye. Eight eyes with glaucoma showed abnormal tissue over the angle and SC. Intrascleral lumen and collector channels were found in both groups. Two patients with severe stenosis of SC during glaucoma angle surgery showed absent SC on AS OCT. CONCLUSIONS: Intraoperative AS OCT identified differences in the presence and characteristics of SC and the anterior chamber angle in pediatric patients with and without glaucoma. Further studies are needed to determine clinical correlation.

Genome-wide association meta-analysis highlights light-induced signaling as a driver for refractive error.

Refractive errors, including myopia, are the most frequent eye disorders worldwide and an increasingly common cause of blindness. This genome-wide association meta-analysis in 160,420 participants and replication in 95,505 participants increased the number of established independent signals from 37 to 161 and showed high genetic correlation between Europeans and Asians (>0.78). Expression experiments and comprehensive in silico analyses identified retinal cell physiology and light processing as prominent mechanisms, and also identified functional contributions to refractive-error development in all cell types of the neurosensory retina, retinal pigment epithelium, vascular endothelium and extracellular matrix. Newly identified genes implicate novel mechanisms such as rod-and-cone bipolar synaptic neurotransmission, anterior-segment morphology and angiogenesis. Thirty-one loci resided in or near regions transcribing small RNAs, thus suggesting a role for post-transcriptional regulation. Our results support the notion that refractive errors are caused by a light-dependent retina-to-sclera signaling cascade and delineate potential pathobiological molecular drivers.

Angiopoietin-1 is required for Schlemm's canal development in mice and humans.

Primary congenital glaucoma (PCG) is a leading cause of blindness in children worldwide and is caused by developmental defects in 2 aqueous humor outflow structures, Schlemm's canal (SC) and the trabecular meshwork. We previously identified loss-of-function mutations in the angiopoietin (ANGPT) receptor TEK in families with PCG and showed that ANGPT/TEK signaling is essential for SC development. Here, we describe roles for the major ANGPT ligands in the development of the aqueous outflow pathway. We determined that ANGPT1 is essential for SC development, and that Angpt1-knockout mice form a severely hypomorphic canal with elevated intraocular pressure. By contrast, ANGPT2 was dispensable, although mice deficient in both Angpt1 and Angpt2 completely lacked SC, indicating that ANGPT2 compensates for the loss of ANGPT1. In addition, we identified 3 human subjects with rare ANGPT1 variants within an international cohort of 284 PCG patients. Loss of function in 2 of the 3 patient alleles was observed by functional analysis of ANGPT1 variants in a combined in silico, in vitro, and in vivo approach, supporting a causative role for ANGPT1 in disease. By linking ANGPT1 with PCG, these results highlight the importance of ANGPT/TEK signaling in glaucoma pathogenesis and identify a candidate target for therapeutic development.

Assaying the Effects of Splice Site Variants by Exon Trapping in a Mammalian Cell Line.

There are several in silico programs that endeavor to predict the functional impact of an individual's sequence variation at splice donor/acceptor sites, but experimental confirmation is problematic without a source of RNA from the individual that carries the variant. With the aid of an exon trapping vector, such as pSPL3, an investigator can test whether a splice site sequence change leads to altered RNA splicing, through expression of reference and variant mini-genes in mammalian cells and analysis of the resultant RNA products.

Exome Sequence Analysis of 14 Families With High Myopia.

Purpose: To identify causal gene mutations in 14 families with autosomal dominant (AD) high myopia using exome sequencing. Methods: Select individuals from 14 large Caucasian families with high myopia were exome sequenced. Gene variants were filtered to identify potential pathogenic changes. Sanger sequencing was used to confirm variants in original DNA, and to test for disease cosegregation in additional family members. Candidate genes and chromosomal loci previously associated with myopic refractive error and its endophenotypes were comprehensively screened. Results: In 14 high myopia families, we identified 73 rare and 31 novel gene variants as candidates for pathogenicity. In seven of these families, two of the novel and eight of the rare variants were within known myopia loci. A total of 104 heterozygous nonsynonymous rare variants in 104 genes were identified in 10 out of 14 probands. Each variant cosegregated with affection status. No rare variants were identified in genes known to cause myopia or in genes closest to published genome-wide association study association signals for refractive error or its endophenotypes. Conclusions: Whole exome sequencing was performed to determine gene variants implicated in the pathogenesis of AD high myopia. This study provides new genes for consideration in the pathogenesis of high myopia, and may aid in the development of genetic profiling of those at greatest risk for attendant ocular morbidities of this disorder.

Reduced penetrance in a large Caucasian pedigree with Stickler syndrome.

BACKGROUND: In a four-generation Caucasian family variably diagnosed with autosomal dominant (AD) Stickler or Wagner disease, commercial gene screening failed to identify a mutation in COL2A1 or VCAN. We utilized linkage mapping and exome sequencing to identify the causal variant. MATERIALS AND METHODS: Genomic DNA samples collected from 40 family members were analyzed. A whole-genome linkage scan was performed using Illumina HumanLinkage-24 BeadChip followed by two-point and multipoint linkage analyses using FASTLINK and MERLIN. Exome sequencing was performed on two affected individuals, followed by co-segregation analysis. RESULTS: Parametric multipoint linkage analysis using an AD inheritance model demonstrated HLOD scores > 2.00 at chromosomes 1p36.13-1p36.11 and 12q12-12q14.1. SIMWALK multipoint analysis replicated the peak in chromosome 12q (peak LOD = 1.975). FASTLINK two-point analysis highlighted several clustered chromosome 12q SNPs with HLOD > 1.0. Exome sequencing revealed a novel nonsense mutation (c.115C>T, p.Gln39*) in exon 2 of COL2A1 that is expected to result in nonsense-mediated decay of the RNA transcript. This mutation co-segregated with all clinically affected individuals and seven individuals who were clinically unaffected. CONCLUSIONS: The utility of combining traditional linkage mapping and exome sequencing is highlighted to identify gene mutations in large families displaying a Mendelian inheritance of disease. Historically, nonsense mutations in exon 2 of COL2A1 have been reported to cause a fully penetrant ocular-only Stickler phenotype with few or no systemic manifestations. We report a novel nonsense mutation in exon 2 of COL2A1 that displays incomplete penetrance and/or variable age of onset with extraocular manifestations.

Angiopoietin receptor TEK mutations underlie primary congenital glaucoma with variable expressivity.

Primary congenital glaucoma (PCG) is a devastating eye disease and an important cause of childhood blindness worldwide. In PCG, defects in the anterior chamber aqueous humor outflow structures of the eye result in elevated intraocular pressure (IOP); however, the genes and molecular mechanisms involved in the etiology of these defects have not been fully characterized. Previously, we observed PCG-like phenotypes in transgenic mice that lack functional angiopoietin-TEK signaling. Herein, we identified rare TEK variants in 10 of 189 unrelated PCG families and demonstrated that each mutation results in haploinsufficiency due to protein loss of function. Multiple cellular mechanisms were responsible for the loss of protein function resulting from individual TEK variants, including an absence of normal protein production, protein aggregate formation, enhanced proteasomal degradation, altered subcellular localization, and reduced responsiveness to ligand stimulation. Further, in mice, hemizygosity for Tek led to the formation of severely hypomorphic Schlemm's canal and trabecular meshwork, as well as elevated IOP, demonstrating that anterior chamber vascular development is sensitive to Tek gene dosage and the resulting decrease in angiopoietin-TEK signaling. Collectively, these results identify TEK mutations in patients with PCG that likely underlie disease and are transmitted in an autosomal dominant pattern with variable expressivity.

Mice expressing mutant Trpv4 recapitulate the human TRPV4 disorders.

Activating mutations in transient receptor potential vanilloid family member 4 (Trpv4) are known to cause a spectrum of skeletal dysplasias ranging from autosomal dominant brachyolmia to lethal metatropic dysplasia. To develop an animal model of these disorders, we created transgenic mice expressing either wild-type or mutant TRPV4. Mice transgenic for wild-type Trpv4 showed no morphological changes at embryonic day 16.5 but did have a delay in bone mineralization. Overexpression of a mutant TRPV4 caused a lethal skeletal dysplasia that phenocopied many abnormalities associated with metatropic dysplasia in humans, including dumbbell-shaped long bones, a small ribcage, abnormalities in the autopod, and abnormal ossification in the vertebrae. The difference in phenotype between embryos transgenic for wild-type or mutant Trpv4 demonstrates that an increased amount of wild-type protein can be tolerated and that an activating mutation of this protein is required to produce a skeletal dysplasia phenotype.

Exome sequencing identifies PDE4D mutations in acrodysostosis.

Acrodysostosis is a dominantly-inherited, multisystem disorder characterized by skeletal, endocrine, and neurological abnormalities. To identify the molecular basis of acrodysostosis, we performed exome sequencing on five genetically independent cases. Three different missense mutations in PDE4D, which encodes cyclic AMP (cAMP)-specific phosphodiesterase 4D, were found to be heterozygous in three of the cases. Two of the mutations were demonstrated to have occurred de novo, providing strong genetic evidence of causation. Two additional cases were heterozygous for de novo missense mutations in PRKAR1A, which encodes the cAMP-dependent regulatory subunit of protein kinase A and which has been recently reported to be the cause of a form of acrodysostosis resistant to multiple hormones. These findings demonstrate that acrodysostosis is genetically heterogeneous and underscore the exquisite sensitivity of many tissues to alterations in cAMP homeostasis.

Dominant and recessive forms of fibrochondrogenesis resulting from mutations at a second locus, COL11A2.

Fibrochondrogenesis is a severe, recessively inherited skeletal dysplasia shown to result from mutations in the gene encoding the proα1(XI) chain of type XI collagen, COL11A1. The first of two cases reported here was the affected offspring of first cousins and sequence analysis excluded mutations in COL11A1. Consequently, whole-genome SNP genotyping was performed to identify blocks of homozygosity, identical-by-descent, wherein the disease locus would reside. COL11A1 was not within a region of homozygosity, further excluding it as the disease locus, but the gene encoding the proα2(XI) chain of type XI collagen, COL11A2, was located within a large region of homozygosity. Sequence analysis identified homozygosity for a splice donor mutation in intron 18. Exon trapping demonstrated that the mutation resulted in skipping of exon 18 and predicted deletion of 18 amino acids from the triple helical domain of the protein. In the second case, heterozygosity for a de novo 9 bp deletion in exon 40 of COL11A2 was identified, indicating that there are autosomal dominant forms of fibrochondrogenesis. These findings thus demonstrate that fibrochondrogenesis can result from either recessively or dominantly inherited mutations in COL11A2.

Fibrochondrogenesis results from mutations in the COL11A1 type XI collagen gene.

Fibrochondrogenesis is a severe, autosomal-recessive, short-limbed skeletal dysplasia. In a single case of fibrochondrogenesis, whole-genome SNP genotyping identified unknown ancestral consanguinity by detecting three autozygous regions. Because of the predominantly skeletal nature of the phenotype, the 389 genes localized to the autozygous intervals were prioritized for mutation analysis by correlation of their expression with known cartilage-selective genes via the UCLA Gene Expression Tool, UGET. The gene encoding the α1 chain of type XI collagen (COL11A1) was the only cartilage-selective gene among the three candidate intervals. Sequence analysis of COL11A1 in two genetically independent fibrochondrogenesis cases demonstrated that each was a compound heterozygote for a loss-of-function mutation on one allele and a mutation predicting substitution for a conserved triple-helical glycine residue on the other. The parents who were carriers of missense mutations had myopia. Early-onset hearing loss was noted in both parents who carried a loss-of-function allele, suggesting COL11A1 as a locus for mild, dominantly inherited hearing loss. These findings identify COL11A1 as a locus for fibrochondrogenesis and indicate that there might be phenotypic manifestations among carriers.

A recessive skeletal dysplasia, SEMD aggrecan type, results from a missense mutation affecting the C-type lectin domain of aggrecan.

Analysis of a nuclear family with three affected offspring identified an autosomal-recessive form of spondyloepimetaphyseal dysplasia characterized by severe short stature and a unique constellation of radiographic findings. Homozygosity for a haplotype that was identical by descent between two of the affected individuals identified a locus for the disease gene within a 17.4 Mb interval on chromosome 15, a region containing 296 genes. These genes were assessed and ranked by cartilage selectivity with whole-genome microarray data, revealing only two genes, encoding aggrecan and chondroitin sulfate proteoglycan 4, that were selectively expressed in cartilage. Sequence analysis of aggrecan complementary DNA from an affected individual revealed homozygosity for a missense mutation (c.6799G --> A) that predicts a p.D2267N amino acid substitution in the C-type lectin domain within the G3 domain of aggrecan. The D2267 residue is predicted to coordinate binding of a calcium ion, which influences the conformational binding loops of the C-type lectin domain that mediate interactions with tenascins and other extracellular-matrix proteins. Expression of the normal and mutant G3 domains in mammalian cells showed that the mutation created a functional N-glycosylation site but did not adversely affect protein trafficking and secretion. Surface-plasmon-resonance studies showed that the mutation influenced the binding and kinetics of the interactions between the aggrecan G3 domain and tenascin-C. These findings identify an autosomal-recessive skeletal dysplasia and a significant role for the aggrecan C-type lectin domain in regulating endochondral ossification and, thereby, height.

Sequencing EVC and EVC2 identifies mutations in two-thirds of Ellis-van Creveld syndrome patients.

Ellis-van Creveld syndrome (EvC) is caused by mutations in EVC and EVC2, genes in a divergent orientation separated by only 2.6 kb. We systematically sought mutations in both genes in a panel of 65 affected individuals to assess the proportion of cases resulting from mutations in each gene. We PCR amplified and sequenced the coding exons of both genes. We investigated mutations that could affect splicing by in vitro splicing assays and cDNA analysis. We have identified EVC mutations in 20 cases (31%); in all of these we have detected the mutation on each allele. We have identified EVC2 mutations in 25 cases (38%); in 22 of these we have isolated a mutation on each allele. The majority of the mutations introduce a premature termination codon. We sequenced the region between the two genes in 10 of the 20 cases in which we had not identified a mutation in either gene, revealing only one SNP that was not a common polymorphism. As we have not identified mutations in either gene in 20 cases (31%) it is possible that there is further genetic heterogeneity.

Mutations in two nonhomologous genes in a head-to-head configuration cause Ellis-van Creveld syndrome.

Ellis-van Creveld syndrome (EvC) is an autosomal recessive skeletal dysplasia. Elsewhere, we described mutations in EVC in patients with this condition (Ruiz-Perez et al. 2000). We now report that mutations in EVC2 also cause EvC. These two genes lie in a head-to-head configuration that is conserved from fish to man. Affected individuals with mutations in EVC and EVC2 have the typical spectrum of features and are phenotypically indistinguishable.