Search

Patients with inherited retinal dystrophies (IRDs) were recruited from two understudied populations: Mexico and Pakistan as well as a third well-studied population of European Americans to define the genetic architecture of IRD by performing whole-genome sequencing (WGS). Whole-genome analysis was performed on 409 individuals from 108 unrelated pedigrees with IRDs. All patients underwent an ophthalmic evaluation to establish the retinal phenotype. Although the 108 pedigrees in this study had previously been examined for mutations in known IRD genes using a wide range of methodologies including targeted gene(s) or mutation(s) screening, linkage analysis and exome sequencing, the gene mutations responsible for IRD in these 108 pedigrees were not determined. WGS was performed on these pedigrees using Illumina X10 at a minimum of 30X depth. The sequence reads were mapped against hg19 followed by variant calling using GATK. The genome variants were annotated using SnpEff, PolyPhen2, and CADD score; the structural variants (SVs) were called using GenomeSTRiP and LUMPY. We identified potential causative sequence alterations in 61 pedigrees (57%), including 39 novel and 54 reported variants in IRD genes. For 57 of these pedigrees the observed genotype was consistent with the initial clinical diagnosis, the remaining 4 had the clinical diagnosis reclassified based on our findings. In seven pedigrees (12%) we observed atypical causal variants, i.e. unexpected genotype(s), including 4 pedigrees with causal variants in more than one IRD gene within all affected family members, one pedigree with intrafamilial genetic heterogeneity (different affected family members carrying causal variants in different IRD genes), one pedigree carrying a dominant causative variant present in pseudo-recessive form due to consanguinity and one pedigree with a de-novo variant in the affected family member. Combined atypical and large structural variants contributed to about 20% of cases. Among the novel mutations, 75% were detected in Mexican and 50% found in European American pedigrees and have not been reported in any other population while only 20% were detected in Pakistani pedigrees and were not previously reported. The remaining novel IRD causative variants were listed in gnomAD but were found to be very rare and population specific. Mutations in known IRD associated genes contributed to pathology in 63% Mexican, 60% Pakistani and 45% European American pedigrees analyzed. Overall, contribution of known IRD gene variants to disease pathology in these three populations was similar to that observed in other populations worldwide. This study revealed a spectrum of mutations contributing to IRD in three populations, identified a large proportion of novel potentially causative variants that are specific to the corresponding population or not reported in gnomAD and shed light on the genetic architecture of IRD in these diverse global populations.

We investigate mutations in trβ2, a splice variant of thrb, identifying changes in function, structure, and behavior in larval and adult zebrafish retinas. Two N-terminus CRISPR mutants were identified. The first is a 6BP+1 insertion deletion frameshift resulting in a truncated protein. The second is a 3BP in frame deletion with intact binding domains. ERG recordings of isolated cone signals showed that the 6BP+1 mutants did not respond to red wavelengths of light while the 3BP mutants did respond. 6BP+1 mutants lacked optomotor and optokinetic responses to red/black and green/black contrasts. Both larval and adult 6BP+1 mutants exhibit a loss of red-cone contribution to the ERG and an increase in UV-cone contribution. Transgenic reporters show loss of cone trβ2 activation in the 6BP+1 mutant but increase in the density of cones with active blue, green, and UV opsin genes. Antibody reactivity for red-cone LWS1 and LWS2 opsin was absent in the 6BP+1 mutant, as was reactivity for arrestin3a. Our results confirm a critical role for trβ2 in long-wavelength cone development.

We identified a homozygous missense alteration (c.75C>A, p.D25E) in CLCC1, encoding a presumptive intracellular chloride channel highly expressed in the retina, associated with autosomal recessive retinitis pigmentosa (arRP) in eight consanguineous families of Pakistani descent. The p.D25E alteration decreased CLCC1 channel function accompanied by accumulation of mutant protein in granules within the ER lumen, while siRNA knockdown of CLCC1 mRNA induced apoptosis in cultured ARPE-19 cells. TALEN KO in zebrafish was lethal 11 days post fertilization. The depressed electroretinogram (ERG) cone response and cone spectral sensitivity of 5 dpf KO zebrafish and reduced eye size, retinal thickness, and expression of rod and cone opsins could be rescued by injection of wild type CLCC1 mRNA. Clcc1+/- KO mice showed decreased ERGs and photoreceptor number. Together these results strongly suggest that intracellular chloride transport by CLCC1 is a critical process in maintaining retinal integrity, and CLCC1 is crucial for survival and function of retinal cells.

Rare germ-line mutations in the coding regions of the human EPHA2 gene (EPHA2) have been associated with inherited forms of pediatric cataract, whereas, frequent, non-coding, single nucleotide variants (SNVs) have been associated with age-related cataract. Here we sought to determine if germ-line EPHA2 coding SNVs were associated with age-related cataract in a case-control DNA panel (> 50 years) and if somatic EPHA2 coding SNVs were associated with lens aging and/or cataract in a post-mortem lens DNA panel (> 48 years). Micro-fluidic PCR amplification followed by targeted amplicon (exon) next-generation (deep) sequencing of EPHA2 (17-exons) afforded high read-depth coverage (1000x) for > 82% of reads in the cataract case-control panel (161 cases, 64 controls) and > 70% of reads in the post-mortem lens panel (35 clear lens pairs, 22 cataract lens pairs). Novel and reference (known) missense SNVs in EPHA2 that were predicted in silico to be functionally damaging were found in both cases and controls from the age-related cataract panel at variant allele frequencies (VAFs) consistent with germ-line transmission (VAF > 20%). Similarly, both novel and reference missense SNVs in EPHA2 were found in the post-mortem lens panel at VAFs consistent with a somatic origin (VAF > 3%). The majority of SNVs found in the cataract case-control panel and post-mortem lens panel were transitions and many occurred at di-pyrimidine sites that are susceptible to ultraviolet (UV) radiation induced mutation. These data suggest that novel germ-line (blood) and somatic (lens) coding SNVs in EPHA2 that are predicted to be functionally deleterious occur in adults over 50 years of age. However, both types of EPHA2 coding variants were present at comparable levels in individuals with or without age-related cataract making simple genotype-phenotype correlations inconclusive.

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

This study was initiated to identify causal mutations responsible for autosomal recessive congenital cataracts in consanguineous familial cases.Affected individuals underwent a detailed ophthalmological and clinical examination, and slit-lamp photographs were ascertained for affected individuals who have not yet been operated for the removal of the cataractous lens. Blood samples were obtained, and genomic DNA was extracted from white blood cells. A genome-wide scan was completed with short tandem repeat (STR) markers, and the logarithm of odds (LOD) scores were calculated. Protein coding exons of CRYAB were sequenced, bi-directionally. Evolutionary conservation was investigated by aligning CRYAB orthologues, and the expression of Cryab in embryonic and postnatal mice lens was investigated with TaqMan probe.The clinical and ophthalmological examinations suggested that all affected individuals had nuclear cataracts. Genome-wide linkage analysis suggested a potential region on chromosome 11q23 harboring CRYAB. DNA sequencing identified a missense variation: c.34C>T (p.R12C) in CRYAB that segregated with the disease phenotype in the family. Subsequent interrogation of our entire cohort of familial cases identified a second familial case localized to chromosome 11q23 harboring a c.31C>T (p.R11C) mutation. In silico analyses suggested that the mutations identified in familial cases, p.R11C and p.R12C will not be tolerated by the three-dimensional structure of CRYAB. Real-time PCR analysis identified the expression of Cryab in mouse lens as early as embryonic day 15 (E15) that increased significantly until postnatal day 6 (P6) with steady level of expression thereafter.Here, we report two novel missense mutations, p.R11C and p.R12C, in CRYAB associated with autosomal recessive congenital nuclear cataracts.

PURPOSE:To define the molecular basis of retinal degeneration in consanguineous Pakistani pedigrees with early onset retinal degeneration. METHODS:A cohort of 277 individuals representing 26 pedigrees from the Punjab province of Pakistan was analyzed. Exomes were captured with commercial kits and sequenced on an Illumina HiSeq 2500. Candidate variants were identified using standard tools and analyzed using exomeSuite to detect all potentially pathogenic changes in genes implicated in retinal degeneration. Segregation analysis was performed by dideoxy sequencing and novel variants were additionally investigated for their presence in ethnicity-matched controls. RESULTS:We identified a total of nine causal mutations, including six novel variants in RPE65, LCA5, USH2A, CNGB1, FAM161A, CERKL and GUCY2D as the underlying cause of inherited retinal degenerations in 13 of 26 pedigrees. In addition to the causal variants, a total of 200 variants each observed in five or more unrelated pedigrees investigated in this study that were absent from the dbSNP, HapMap, 1000 Genomes, NHLBI ESP6500, and ExAC databases were identified, suggesting that they are common in, and unique to the Pakistani population. CONCLUSIONS:We identified causal mutations associated with retinal degeneration in nearly half of the pedigrees investigated in this study through next generation whole exome sequencing. All novel variants detected in this study through exome sequencing have been cataloged providing a reference database of variants common in, and unique to the Pakistani population.

Congenital cataracts are one of the leading causes of visual impairment and blindness in children, and genetic factors play an important role in their development. This study aimed to identify the genetic defects associated with autosomal dominant congenital progressive punctate cataracts in a Chinese family and to explore the potential pathogenesis. Detailed family history and clinical data were recorded, and all the family members' blood samples were collected for DNA extraction. Linkage analysis was performed by microsatellite markers that are associated with punctate cataracts, and logarithm (base 10) of odds (LOD) scores were calculated using the LINKAGE program. Positive two-point LOD scores were obtained at markers D12S1622 (Zmax = 2.71 at θ = 0.0), D12S1724 (Zmax = 2.71 at θ = 0.0), and D12S90 (Zmax = 2.71 at θ = 0.0), which flank the major intrinsic protein of lens fiber (MIP) gene on chromosomal region 12q13. Direct sequencing of the encoding region of the MIP gene revealed a novel mutation (G>D) in exon 4 at nucleotide 644, which caused a substitution of glycine to aspartic acid at codon 215 (p.G215D) for the MIP protein. The mutation cosegregated with all patients with congenital progressive punctate cataracts, but it was absent in the healthy members. Bioinformatics analysis predicted that the mutation affects the function of the MIP protein. The wild type (WT) and G215D mutant of MIP were transfected with green fluorescent protein (GFP) into Hela cells separately, and it was found that the G215D mutant was aberrantly located in the cytoplasm instead of in the plasma membrane. In summary, our study presented genetic and functional evidence linking the new MIP mutation of G215D to autosomal dominant congenital cataracts, which adds to the list of MIP mutations linked to congenital progressive punctate cataracts.

Small heat shock proteins (sHsps) maintain cellular homeostasis by preventing stress and disease-induced protein aggregation. While it is known that hydrophobicity impacts the ability of sHsps to bind aggregation-prone denaturing proteins, the complex quaternary structure of globular sHsps has made understanding the significance of specific changes in hydrophobicity difficult. Here we used recombinant protein of the lenticular sHsp α A-crystallin from six teleost fishes environmentally adapted to temperatures ranging from -2°C to 40°C to identify correlations between physiological temperature, protein stability and chaperone-like activity. Using sequence and structural modeling analysis we identified specific amino acid differences between the warm adapted zebrafish and cold adapted Antarctic toothfish that could contribute to these correlations and validated the functional consequences of three specific hydrophobicity-altering amino acid substitutions in αA-crystallin. Site directed mutagenesis of three residues in the zebrafish (V62T, C143S, T147V) confirmed that each impacts either protein stability or chaperone-like activity or both, with the V62T substitution having the greatest impact. Our results indicate a role for changing hydrophobicity in the thermal adaptation of α A-crystallin and suggest ways to produce sHsp variants with altered chaperone-like activity. These data also demonstrate that a comparative approach can provide new information about sHsp function and evolution.

The purpose of this study was to confirm previously reported associations of common variants in or near CDC7/TGFBR3, ZP4, SRBD1, ELOVL5, CAV1/CAV2, TLR4, CDKN2B, CDKN2B-AS1, ATOH7, PLXDC2, TMTC2, SIX1, and CARD10, with primary open angle glaucoma (POAG) in the Afro-Caribbean population of Barbados, West Indies. A total of 437 unrelated subjects from the Barbados Family Study of Open Angle Glaucoma (BFSG), including 272 with POAG and 165 unaffected individuals were included in this study. Eighteen SNPs were genotyped by using the multiplex SNaPshot method. Allelic, genotypic and model-based (dominant, recessive, and additive) associations of the SNPs with POAG were analyzed using Chi-squared tests and logistic regression. SNP rs1063192 (near CDKN2B) was found to be significantly associated with POAG (allelic P = 0.0008, genotypic P = 0.0029), and the minor allele C of rs1063192 was protective against POAG (OR = 0.39; 95%CI = 0.22-0.69). Suggestive association was also noted for rs7916697 (near ATHO7, allelic P = 0.0096, genotypic P = 0.01) with the minor allele being protective (OR = 0.67; 95% CI = 0.50-0.91), although this finding did not withstand correction for multiple testing. However, a significant interactive effect on POAG risk was identified between rs1063192 and rs7916697 (P-interaction = 2.80 × 10(-5)). Individuals with the rs1063192 protective genotype CC or CT and also rs7916697 genotypes GG or GA show a significantly decreased risk of POAG (OR = 0.17, 95%CI: 0.07-0.41). Our study confirms the significant association between SNP rs1063192 (CDKN2B, previously shown to influence vertical cup-to-disc ratio and POAG at 9p21) and POAG in the Afro-Caribbean population of Barbados. The minor allele of rs1063192 interacts with that of rs7916697 (ATOH7)) to reduce POAG risk. Our results also suggest that rs1063912 is a common protective variant for POAG in populations of African as well as European descent.

BACKGROUND: Leber congenital amaurosis (LCA) is the earliest onset and most severe form of hereditary retinal dystrophy. So far, full spectrum of variations in the 15 genes known to cause LCA has not been systemically evaluated in East Asians. Therefore, we performed comprehensive detection of variants in these 15 genes in 87 unrelated Han Chinese patients with LCA. METHODOLOGY/PRINCIPAL FINDINGS: The 51 most frequently mutated exons and introns in the 15 genes were selected for an initial scan using cycle sequencing. All the remaining exons in 11 of the 15 genes were subsequently sequenced. Fifty-three different variants were identified in 44 of the 87 patients (50.6%), involving 78 of the 88 alleles (11 homozygous and 56 heterozygous variants). Of the 53 variants, 35 (66%) were novel pathogenic mutations. In these Chinese patients, variants in GUCY2D are the most common cause of LCA (16.1% cases), followed by CRB1 (11.5%), RPGRIP1 (8%), RPE65 (5.7%), SPATA7 (4.6%), CEP290 (4.6%), CRX (3.4%), LCA5 (2.3%), MERTK (2.3%), AIPL1 (1.1%), and RDH12 (1.1%). This differs from the variation spectrum described in other populations. An initial scan of 55 of 215 PCR amplicons, including 214 exons and 1 intron, detected 83.3% (65/78) of the mutant alleles ultimately found in these 87 patients. In addition, sequencing only 9 exons would detect over 50% of the identified variants and require less than 5% of the labor and cost of comprehensive sequencing for all exons. CONCLUSIONS/SIGNIFICANCE: Our results suggest that specific difference in the variation spectrum found in LCA patients from the Han Chinese and other populations are related by ethnicity. Sequencing exons in order of decreasing risk is a cost-effective way to identify causative mutations responsible for LCA, especially in the context of genetic counseling for individual patients in a clinical setting.