Search

Aims/hypothesis: ALR/Lt, a mouse strain with strong resistance to type 1 diabetes, is closely related to autoimmune type 1 diabetes-prone NOD/Lt mice. ALR pancreatic beta cells are resistant to the beta cell toxin alloxan, combinations of cytotoxic cytokines, and diabetogenic NOD T-cell lines. Reciprocal F1 hybrids between either ALR and NOD or ALR and NON/Lt, showed that alloxan resistance was transmitted to Fl progeny only when ALR was the maternal parent. Here we show that the mitochondrial genome (mtDNA) of ALR mice contributes resistance to diabetes. Methods: When Fl progeny from reciprocal outcrosses between ALR and NOD were backcrossed to NOD, a four-fold lower frequency of spontaneous type 1 diabetes development occurred when ALR contributed the mtDNA. Because of the apparent inter- action between nuclear and mtDNA, the mitochondrial genomes were sequenced. Results: An ALR-specific sequence variation in the mt-Nd2 gene producing a leucine to methionine substitution at amino acid residue 276 in the NADH dehydrogenase 2 was discovered. An isoleucine to valine mutation in the mt-Co3 gene encoding COX3 distinguished ALR and NOD from NON and ALS. All four strains were distinguished by variation in a mt-encoded arginyl tRNA polyadenine tract. Shared alleles of mt-Co3 and mt-Tr comparing NOD and ALR allowed for exclusion of these two genes as candidates, implicating the mt-Nd2 variation as a potential ALR-derived type 1 diabetes protective gene. Conclusions/interpretation: The unusual resistance of ALR mice to both ROS-mediated and autoimmune type 1 diabete stresses reflects an interaction between the nuclear and mt genomes. The latter contribution is most likely via a single nucleotide polymorphism in mt-Nd2. [ABSTRACT FROM AUTHOR]

Keratoconus is characterised by reduced rigidity of the cornea with distortion and focal thinning that causes blurred vision, however, the pathogenetic mechanisms are unknown. It can lead to severe visual morbidity in children and young adults and is a common indication for corneal transplantation worldwide. Here we report the first large scale genome-wide association study of keratoconus including 4,669 cases and 116,547 controls. We have identified significant association with 36 genomic loci that, for the first time, implicate both dysregulation of corneal collagen matrix integrity and cell differentiation pathways as primary disease-causing mechanisms. The results also suggest pleiotropy, with some disease mechanisms shared with other corneal diseases, such as Fuchs endothelial corneal dystrophy. The common variants associated with keratoconus explain 12.5% of the genetic variance, which shows potential for the future development of a diagnostic test to detect susceptibility to disease.

In the past few years we have seen a great acceleration of discoveries in the field of keratoconus including new treatments, diagnostic tools, genomic and molecular determinants of disease risk. Recent genome-wide association studies (GWAS) of keratoconus cases and population wide studies of variation in central corneal thickness and in corneal biomechanical properties confirmed already identified genes and found many new susceptibility variants and biological pathways. Recent findings in genetic determinants of familial keratoconus revealed functionally important variants and established first mouse model of keratoconus. Latest transcriptomic and expression studies started assessing novel non-coding RNA targets in addition to identifying tissue specific effects of coding genes. First genomic insights into better prediction of treatment outcomes are bringing the advent of genomic medicine into keratoconus clinical practice., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Importance: Keratoconus is a condition in which the cornea progressively thins and protrudes in a conical shape, severely affecting refraction and vision. It is a major indication for corneal transplant. To discover new genetic loci associated with keratoconus and better understand the causative mechanism of this disease, we performed a genome-wide association study on patients with keratoconus., Objective: To identify genetic susceptibility regions for keratoconus in the human genome., Design, Setting, and Participants: This study was conducted with data from eye clinics in Australia, the United States, and Northern Ireland. The discovery cohort of individuals with keratoconus and control participants from Australia was genotyped using the Illumina HumanCoreExome single-nucleotide polymorphism array. After quality control and data cleaning, genotypes were imputed against the 1000 Genomes Project reference panel (phase III; version 5), and association analyses were completed using PLINK version 1.90. Single-nucleotide polymorphisms with P < 1.00 × 10-6 were assessed for replication in 3 additional cohorts. Control participants were drawn from the cohorts of the Blue Mountains Eye Study and a previous study of glaucoma. Replication cohorts were from a previous keratoconus genome-wide association study data set from the United States, a cohort of affected and control participants from Australia and Northern Ireland, and a case-control cohort from Victoria, Australia. Data were collected from January 2006 to March 2019., Main Outcomes and Measures: Associations between keratoconus and 6 252 612 genetic variants were estimated using logistic regression after adjusting for ancestry using the first 3 principal components., Results: The discovery cohort included 522 affected individuals and 655 control participants, while the replication cohorts included 818 affected individuals (222 from the United States, 331 from Australia and Northern Ireland, and 265 from Victoria, Australia) and 3858 control participants (2927 from the United States, 229 from Australia and Northern Ireland, and 702 from Victoria, Australia). Two novel loci reached genome-wide significance (defined as P < 5.00 × 10-8), with a P value of 7.46 × 10-9 at rs61876744 in patatin-like phospholipase domain-containing 2 gene (PNPLA2) on chromosome 11 and a P value of 6.35 × 10-12 at rs138380, 2.2 kb upstream of casein kinase I isoform epsilon gene (CSNK1E) on chromosome 22. One additional locus was identified with a P value less than 1.00 × 10-6 in mastermind-like transcriptional coactivator 2 (MAML2) on chromosome 11 (P = 3.91 × 10-7). The novel locus in PNPLA2 reached genome-wide significance in an analysis of all 4 cohorts (P = 2.45 × 10-8)., Conclusions and Relevance: In this relatively large keratoconus genome-wide association study, we identified a genome-wide significant locus for keratoconus in the region of PNPLA2 on chromosome 11.

Keratoconus (KC) is the most common corneal ectatic disorder affecting >300,000 people in the US. KC normally has its onset in adolescence, progressively worsening through the third to fourth decades of life. KC patients report significant impaired vision-related quality of life. Genetic factors play an important role in KC pathogenesis. To identify novel genes in familial KC patients, we performed whole exome and genome sequencing in a four-generation family. We identified potential variants in the PPIP5K2 and PCSK1 genes. Using in vitro cellular model and in vivo gene-trap mouse model, we found critical evidence to support the role of PPIP5K2 in normal corneal function and KC pathogenesis. The gene-trap mouse showed irregular corneal surfaces and pathological corneal thinning resembling KC. For the first time, we have integrated corneal tomography and pachymetry mapping into characterization of mouse corneal phenotypes which could be widely implemented in basic and translational research for KC diagnosis and therapy in the future.

Purpose: To report a case of bilateral and repetitive corneal perforations after corneal cross-linking (CXL) for keratoconus in a woman harboring potentially pathogenic variants in the ZNF469 gene and to characterize the keratoconus phenotype in this woman and her daughter who shared the same ZNF469 mutations., Methods: Clinical characterization of the proband and her daughter followed by sequencing of the genes associated with brittle cornea syndrome, ZNF469 and PRDM5, in both individuals., Results: An Ashkenazi Jewish woman in her sixth decade presented with diffuse corneal thinning and progressive steepening consistent with keratoconus. After CXL, epithelium-off in the first eye and epithelium-on in the second, she developed spontaneous corneal perforations in each eye. Her daughter in her fourth decade demonstrated a similar pattern of diffuse corneal thinning and progressive corneal steepening but did not undergo CXL and did not develop corneal perforation. Screening of the ZNF469 and PRDM5 genes revealed 3 missense ZNF469 variants (c.2035G>A, c.10244G>C, and c.11119A>G) in cis arrangement on 1 allele of ZNF469 in both proband and her daughter. Although the 3 variants share low (<0.01) global minor allele frequencies, each has significantly higher minor allele frequencies (0.01-0.03) in the Ashkenazi Jewish population, leading to uncertainty regarding a pathogenic role for the identified variants., Conclusions: CXL may be associated with the development of corneal perforation in particular at-risk individuals with keratoconus. Identifying clinical and genetic risk factors, including screening of ZNF469 and PRDM5, may be useful in the prevention of significant complications after CXL.

Purpose: Keratoconus (KC) is the most common corneal ectasia. We aimed to determine the differential expression of coding and long noncoding RNAs (lncRNAs) in human corneas affected with KC., Methods: From the corneas of 10 KC patients and 8 non-KC healthy controls, 200 ng total RNA was used to prepare sequencing libraries with the SMARTer Stranded RNA-Seq kit after ribosomal RNA depletion, followed by paired-end 50-bp sequencing with Illumina Sequencer. Differential analysis was done using TopHat/Cufflinks with a gene file from Ensembl and a lncRNA file from NONCODE. Pathway analysis was performed using WebGestalt. Using the expression level of differentially expressed coding and noncoding RNAs in each sample, we correlated their expression levels in KC and controls separately and identified significantly different correlations in KC against controls followed by visualization using Cytoscape., Results: Using |fold change| ≥ 2 and a false discovery rate ≤ 0.05, we identified 436 coding RNAs and 584 lncRNAs with differential expression in the KC-affected corneas. Pathway analysis indicated the enrichment of genes involved in extracellular matrix, protein binding, glycosaminoglycan binding, and cell migration. Our correlation analysis identified 296 pairs of significant KC-specific correlations containing 117 coding genes enriched in functions related to cell migration/motility, extracellular space, cytokine response, and cell adhesion. Our study highlighted the potential roles of several genes (CTGF, SFRP1, AQP5, lnc-WNT4-2:1, and lnc-ALDH3A2-2:1) and pathways (TGF-β, WNT signaling, and PI3K/AKT pathways) in KC pathogenesis., Conclusions: Our RNA-Seq-based differential expression and correlation analyses have identified many potential KC contributing coding and noncoding RNAs.

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.

Purpose: To test candidate genes TSC1 and TSC2 in a family affected by tuberous sclerosis complex (TSC) where proband was also diagnosed with bilateral keratoconus (KC) and to test the hypothesis that defects in the same gene may lead to a nonsyndromic KC., Methods: Next-generation sequencing of TSC1 and TSC2 genes was performed in a proband affected by TSC and KC. Identified mutation was confirmed by Sanger DNA sequencing. Whole exome sequencing (WES) was performed in patients with nonsyndromic KC. Sanger DNA sequencing was used to confirm WES results and to screen additional patients. RT-PCR was used to investigate TSC1 expression in seven normal human corneas and eight corneas from patients with KC. Various in silico tools were employed to model functional consequences of identified mutations., Results: A heterozygous nonsense TSC1 mutation g.132902703C>T (c.2293C>T, p.Gln765Ter) was identified in a patient with TSC and KC. Two heterozygous missense TSC1 variants g.132896322A>T (c.3408A>T, p.Asp1136Glu) and g.132896452G>A (c.3278G>A, p.Arg1093Gln) were identified in three patients with nonsyndromic KC. Two mutations were not present in The Genome Aggregation (GnomAD), The Exome Aggregation (ExAC), and 1000 Genomes (1000G) databases, while the third one was present in GnomAD and 1000G with minor allele frequencies (MAF) of 0.00001 and 0.0002, respectively. We found TSC1 expressed in normal corneas and KC corneas, albeit with various levels., Conclusions: Here for the first time we found TSC1 gene to be involved in bilateral KC and TSC as well as with nonsyndromic KC, supporting the hypothesis that diverse germline mutations of the same gene can cause genetic disorders with overlapping clinical features.

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.

Mitochondrial myopathy with lactic acidosis and sideroblastic anemia (MLASA) is an oxidative phosphorylation disorder, with primary clinical manifestations of myopathic exercise intolerance and a macrocytic sideroblastic anemia. One cause of MLASA is recessive mutations in PUS1, which encodes pseudouridine (Ψ) synthase 1 (Pus1p). Here we describe a mouse model of MLASA due to mutations in PUS1. As expected, certain Ψ modifications were missing in cytoplasmic and mitochondrial tRNAs from Pus1(-/-) animals. Pus1(-/-) mice were born at the expected Mendelian frequency and were non-dysmorphic. At 14 weeks the mutants displayed reduced exercise capacity. Examination of tibialis anterior (TA) muscle morphology and histochemistry demonstrated an increase in the cross sectional area and proportion of myosin heavy chain (MHC) IIB and low succinate dehydrogenase (SDH) expressing myofibers, without a change in the size of MHC IIA positive or high SDH myofibers. Cytochrome c oxidase activity was significantly reduced in extracts from red gastrocnemius muscle from Pus1(-/-) mice. Transmission electron microscopy on red gastrocnemius muscle demonstrated that Pus1(-/-) mice also had lower intermyofibrillar mitochondrial density and smaller mitochondria. Collectively, these results suggest that alterations in muscle metabolism related to mitochondrial content and oxidative capacity may account for the reduced exercise capacity in Pus1(-/-) mice.

Aim: To identify changes in the expression of genes coding for extracellular matrix (ECM) proteins in patients with non-inflammatory corneal disorder keratoconus (KC), patients with corneal scarring, and normal controls., Materials and Methods: Total RNA extracted from corneal tissue of 13 KC patients, 2 patients with corneal scaring and 4 normal controls was analyzed using Human Extracellular Matrix & Adhesion Molecules Profiler PCR Array. Statistically significant changes in gene expression were identified using the Data Analysis software., Results: Comparison of KC and control corneas with thresholds of 1.5 or greater fold change and a p-value of 0.05 or lower, revealed 21 differentially expressed genes, 16 genes were downregulated and 5 were upregulated. Among transcripts downregulated in KC patients we identified THBS1, ADAMTS1, SPP1, several collagens and integrins. We found TGFBI (BIGH3) gene was the most significantly upregulated transcript., Conclusion: Development of keratoconus results in deregulation of gene expression of extracellular matrix and adhesion molecules., Clinical Significance: Downregulation of collagens and upregulation of TGFBI repeatedly identified in KC patients may be used as clinical markers of the disease., Competing Interests: Conflict of interest statement All authors certify that they have no conflicts of interest.

A c.57 C > T mutation in the seed region of MIR184 located at the 15q25.1 chromosomal region has been independently associated with autosomal dominant keratoconus with early-onset anterior polar cataract in the Northern Irish family and with autosomal dominant EDICT (Endothelial Dystrophy, Iris hypoplasia, Congenital cataracts, and stromal Thinning) syndrome. In this study we report a five-generation family originating in Galicia, Spain with early onset cataracts and variable corneal abnormalities which include non-ectatic corneal thinning and severe early-onset keratoconus. We identified a heterozygous c.57 C > T mutation in miR-184 in the proband and two additional affected relatives on the maternal side. This finding represents a third independent occurrence of this mutation in familiar ocular disease thus strengthening the link between miR-184 abnormalities and inherited eye defects.

Purpose: To develop parameters using a combination of optical coherence tomography (OCT) and videokeratography to detect early keratoconus., Methods: Videokeratography, wavefront analysis, and measured OCT indices were performed on 180 normal eyes, 46 eyes with moderate keratoconus, 54 eyes with early keratoconus, 7 eyes with forme fruste keratoconus, and 16 eyes with keratoconus "suspect" to determine the most sensitive parameters for separating these groups., Results: A combination of videokeratography and OCT indices (inferior-superior [I-S] value and minimum pachymetry) was statistically the most significant in separating the keratoconus groups from normal eyes (P < .001). Using a newly derived index, the minimum pachymetry divided by the I-S value (pachymetry/asymmetry [PA]/I-S index) with a cut-off of 100, 100% of early and forme fruste keratoconus could be identified as being abnormal with 7 normals misclassified (misclassification rate 2.7%). By adding keratoconus "suspect" to the analysis and an I-S value of 1.2 as a cut-off point, 5 "suspects" were classified as normal and 11 normals as abnormal (misclassification rate 7.8%). The PA/I-S index, with a cut-off point of 100, reduced this misclassification rate to 4.4%., Conclusions: These results suggest that OCT combined with videokeratography may be more useful for differentiating mild forms of keratoconus than videokeratography alone., (Copyright 2014, SLACK Incorporated.)

Purpose: Keratoconus (KC) is a genetically heterogeneous corneal dystrophy. Previously, we performed 2 genome-wide linkage scans in a 4-generation autosomal dominant pedigree and repeatedly mapped a KC locus to a genomic region located on chromosome 5q overlapping the gene encoding the inhibitor of calpains, calpastatin (CAST). To test whether variants in CAST gene are involved in genetic susceptibility to KC, we performed genetic testing of polymorphic markers in CAST gene in family and case-control panels of patients with KC., Methods: We genotyped single-nucleotide polymorphisms (SNPs) located in CAST gene in 262 patients in 40 white KC families and in a white case-control panel with 304 cases and 518 controls. Generalized estimating equation models accounting for familial correlations implemented in GWAF program were used for association testing in families. Logistic regression models implemented in PLINK were performed to test the associations in case-control samples., Results: Genetic testing of the first set of 7 SNPs in familial samples revealed 2 tentative nominally significant markers (rs4869307, P = 0.03; rs27654, P = 0.07). Additional genotyping of 12 tightly spaced SNPs identified CAST SNP rs4434401 to be associated with KC in both familial and case-control panels with P values of 0.005 and 0.05, respectively, and with combined meta P value of familial and case-control cohorts of 0.002 or after Bonferroni correction of 0.04., Conclusions: Linkage analysis and genetic association support involvement of CAST gene in the genetic susceptibility to KC. In silico analysis of CAST expression suggests differential regulation of calpain/calpastatin system in cornea as a potential mechanism of functional defect.

Purpose: Single nucleotide polymorphisms (SNPs) located near or within the COL5A1 gene, at 9q34.2-q34.3 chromosomal region have been reported in association with central corneal thickness (CCT). Using family linkage analysis, we identified a keratoconus susceptibility locus at 9q34. These findings led us to perform an association study between COL5A1 variation and keratoconus susceptibility., Methods: A Caucasian case-control cohort of 222 keratoconus patients and 3324 controls was selected as the discovery panel. An independent case-control panel of 304 cases and 518 controls and a family panel of 186 subjects were replicated for genotyping and association. Forty-four SNPs (21 for discovery and 23 for fine-mapping) spanning 300 kilobases in and around COL5A1 were genotyped and tested for genetic association. Logistic regression models implemented in PLINK were used to test for association in case controls. Generalized estimating equation models accounting for familial correlations implemented in genome-wide interaction analyses with family data were used for association testing in families., Results: Two CCT associated SNPs (rs1536482 and rs7044529 near and within COL5A1) were identified in the keratoconus discovery cohort (P values of 6.5 × 10(-3) and 7.4 × 10(-3)). SNP rs1536482 was replicated in the second case-control sample (P = 0.02), and SNP rs7044529 was replicated in a keratoconus family panel (P = 0.03). Meta P values of rs1536482 and rs7044529 in the keratoconus cohorts were 1.5 × 10(-4) (odds ratio [OR] = 1.30) and 2.9 × 10(-3) (OR = 1.39). After Bonferroni correction, the association of SNP rs1536482 remained significant (P = 6.5 × 10(-3))., Conclusions: SNPs in the COL5A1 region, which regulate normal variation in CCT, may play a role in the thinning associated with 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.

Purpose: Keratoconus is a bilateral noninflammatory progressive corneal disorder with complex genetic inheritance and a common cause for cornea transplantation in young adults. A genomewide linkage scan in keratoconus families identified a locus at 5q23.2, overlapping the gene coding for the lysyl oxidase (LOX). LOX encodes an enzyme responsible for collagen cross-linking in a variety of tissues including the cornea. Corneal collagen cross-linking with long-wave ultraviolet light and riboflavin is a promising new treatment for keratoconus. To determine whether LOX is a genetic determinant of the pathogenesis of keratoconus, we analyzed association results of LOX polymorphisms in two independent case-control samples and in keratoconus families., Methods: Association results were analyzed of single-nucleotide polymorphisms (SNPs) in the LOX gene from a Genome-Wide Association Study (GWAS) investigation in two independent panels of patients with keratoconus and controls and in keratoconus families., Results: Evidence of association was found at SNPs rs10519694 and rs2956540 located in intron 4 of LOX in the GWAS discovery case-control panel with P values of 2.3×10(-3) and 7×10(-3), respectively. The same two SNPs were found to be associated with keratoconus by family-based association testing with P values of 2.7×10(-3) and 7.7×10(-4), respectively. Meta P values of 4.0×10(-5) and 4.0×10(-7) were calculated for SNPs rs10519694 and rs2956540 by analyzing case-control and family samples simultaneously. Sequencing of LOX exons in a subset of keratoconus patients identified two polymorphisms, rs1800449 and rs2288393, located in LOX transcripts I and II, associated with keratoconus in case-control and family samples with a meta P value of 0.02., Conclusions: Results provided strong genetic evidence that LOX variants lead to increased susceptibility to developing of keratoconus.

Keratoconus is a condition in which the cornea progressively thins over time, and is a major cause for cornea transplantation. To identify keratoconus susceptibility regions, we performed a comprehensive genome-wide association study (GWAS) using a discovery and replication design. A discovery panel of 222 keratoconus Caucasian patients and 3324 Caucasian controls was genotyped using Illumina 370K beadchips. Further associated and fine-mapping single nucleotide polymorphisms (SNPs) (n= 4905) were genotyped in an independent replication case-control panel of 304 cases and 518 controls and a family panel of 307 subjects in 70 families. Logistic regression models implemented in PLINK were performed to test associations in case-control samples with and without principal component (PC) adjustments. Generalized estimation equation models accounting for familial correlations implemented in GWAF were used for association testing in families. No genome-wide associations were identified in the discovery GWAS panel. From the initial testing without adjustments for PCs, the top three SNPs located at 3p26 (rs6442925), 2q21.3 (rs4954218) and 19q13.3 (rs1428642) were identified with unadjusted P-values of 6.5 × 10(-8), 2.4 × 10(-7) and 3.1 × 10(-7), respectively. After adjustments for PCs, rs1428642 became the most significant through the genome with a P-value of 1.4 × 10(-6), while rs6442925 and rs4954218 were less significant (P= 1.9 × 10(-5) and 2.6 × 10(-4)). SNP rs4954218 was confirmed in two independent replication panels with P-values of 0.004 and 0.009, respectively. Meta-analysis revealed a highest association at rs4954218 with adjusted P= 1.6 × 10(-7) (unadjusted P= 1.2 × 10(-9)). These findings suggest SNP rs4954218, located near the RAB3GAP1 gene, previously reported to be associated with corneal malformation, is a potential susceptibility locus for keratoconus., (© The Author 2011. Published by Oxford University Press. All rights reserved.)

Purpose: Keratoconus is a progressive disorder of the cornea that can lead to severe visual impairment or blindness. Although several genomic regions have been linked to rare familial forms of keratoconus, no genes have yet been definitively identified for common forms of the disease., Methods: Two genome-wide association scans were undertaken in parallel. The first used pooled DNA from an Australian cohort, followed by typing of top-ranked single-nucleotide polymorphisms (SNPs) in individual DNA samples. The second was conducted in individually genotyped patients, and controls from the USA. Tag SNPs around the hepatocyte growth factor (HGF) gene were typed in three additional replication cohorts. Serum levels of HGF protein in normal individuals were assessed with ELISA and correlated with genotype., Results: The only SNP observed to be associated in both the pooled discovery and primary replication cohort was rs1014091, located upstream of the HGF gene. The nearby SNP rs3735520 was found to be associated in the individually typed discovery cohort (P = 6.1 × 10(-7)). Genotyping of tag SNPs around HGF revealed association at rs3735520 and rs17501108/rs1014091 in four of the five cohorts. Meta-analysis of all five datasets together yielded suggestive P values for rs3735520 (P = 9.9 × 10(-7)) and rs17501108 (P = 9.9 × 10(-5)). In addition, SNP rs3735520 was found to be associated with serum HGF level in normal individuals (P = 0.036)., Conclusions: Taken together, these results implicate genetic variation at the HGF locus with keratoconus susceptibility.

The homoplasmic mitochondrial A1555G mutation in the 12S rRNA gene leads to a mitochondrial translation disorder associated with deafness. The absence of disease in non-cochlear tissues in all patients, and in the cochlea in some patients, is not well understood. We used a system-based approach, including whole genome expression and biological function analysis, to elucidate the pathways underlying tissue specificity and clinical severity of this condition. Levels of over 48K RNA transcripts from EBV-transformed lymphoblasts of deaf and hearing individuals with the A1555G mutation and controls were obtained. Differentially expressed transcripts were functionally grouped using gene set enrichment analysis. Over 50 RNA binding proteins were differentially expressed between deaf and hearing individuals with the A1555G mutation (P-value of 2.56E-7), confirming previous genetic data implicating this pathway in the determination of the severity of hearing loss. Unexpectedly, the majority of cytoplasmic ribosomal genes were up-regulated in a coordinated fashion in individuals with the A1555G mutation versus controls (P-value of 3.91E-135). This finding was verified through real time RT-PCR, and through measuring of protein levels by flow cytometry. Analysis of expression levels of other differentially expressed genes suggests that this coordinated over-expression of cytoplasmic ribosomal proteins might occur through the Myc/Max pathway. We propose that expression levels of RNA binding proteins help determine the severity of the cochlear phenotype, and that coordinated up-regulation of the cytoplasmic translation apparatus operates as a compensation mechanism in unaffected tissues of patients with maternal deafness associated with the A1555G mutation.

The tissue specificity of mitochondrial diseases is poorly understood. Recently, tissue-specific quantitative differences of the components of the mitochondrial translation system have been found to correlate with disease presentation in fatal hepatopathy caused by mutations in mitochondrial translation factor EFG1. MLASA is an autosomal recessive inherited progressive oxidative phosphorylation disorder that affects muscle and erythroid cells. The disease is caused by the homozygous point mutation C656T (R116W) in the catalytic domain of the pseudouridylate synthase 1 (PUS1) gene, which leads to a complete lack of pseudouridylation at the expected sites in mitochondrial and cytoplasmic tRNAs. Despite the presence of these altered tRNAs, most tissues are unaffected, and even in muscle and erythroid cells the disease phenotype only slowly emerges over the course of years. In order to elucidate intracellular pathways through which the homozygous mutation leads to tissue-restricted phenotype, we performed microarray expression analysis of EBV-transformed lymphoblasts from MLASA patients, heterozygous parents, and controls using human Beadchip microarray with 47,296 transcripts. Genes coding for proteins involved in DNA transcription and its regulation, and metal binding proteins, demonstrated major differences in expression between patients and all other individuals with normal phenotype. Genes coding for ribosomal proteins differed significantly between individual with at least one copy of the mutated PUS1 gene and controls. These findings indicate that the lack of tRNA pseudouridylation can be overcome by compensatory changes in levels of ribosomal proteins, and that the disease phenotype in affected tissues is likely due to pleiotropic effects of PUS1p on non-tRNA molecules involved in DNA transcription and iron metabolism. Similar combinations of mechanisms may play a role in the tissue specificity of other mitochondrial disorders.

The human mitochondrial 12S ribosomal RNA (rRNA) A1555G mutation has been associated with aminoglycoside-induced and nonsyndromic deafness in many families worldwide. Our previous investigation revealed that the A1555G mutation is a primary factor underlying the development of deafness but is not sufficient to produce a deafness phenotype. However, it has been proposed that nuclear-modifier genes modulate the phenotypic manifestation of the A1555G mutation. Here, we identified the nuclear-modifier gene TRMU, which encodes a highly conserved mitochondrial protein related to transfer RNA (tRNA) modification. Genotyping analysis of TRMU in 613 subjects from 1 Arab-Israeli kindred, 210 European (Italian pedigrees and Spanish pedigrees) families, and 31 Chinese pedigrees carrying the A1555G or the C1494T mutation revealed a missense mutation (G28T) altering an invariant amino acid residue (A10S) in the evolutionarily conserved N-terminal region of the TRMU protein. Interestingly, all 18 Arab-Israeli/Italian-Spanish matrilineal relatives carrying both the TRMU A10S and 12S rRNA A1555G mutations exhibited prelingual profound deafness. Functional analysis showed that this mutation did not affect importation of TRMU precursors into mitochondria. However, the homozygous A10S mutation leads to a marked failure in mitochondrial tRNA metabolisms, specifically reducing the steady-state levels of mitochondrial tRNA. As a consequence, these defects contribute to the impairment of mitochondrial-protein synthesis. Resultant biochemical defects aggravate the mitochondrial dysfunction associated with the A1555G mutation, exceeding the threshold for expressing the deafness phenotype. These findings indicate that the mutated TRMU, acting as a modifier factor, modulates the phenotypic manifestation of the deafness-associated 12S rRNA mutations.

Nuclear modifier genes have been proposed to modulate the phenotypic manifestation of human mitochondrial 12S rRNA A1491G mutation associated with deafness in many families world-wide. Here we identified and characterized the putative nuclear modifier gene TRMU encoding a highly conserved mitochondrial protein related to tRNA modification. A 1937bp TRMU cDNA has been isolated and the genomic organization of TRMU has been elucidated. The human TRMU gene containing 11 exons encodes a 421 residue protein with a strong homology to the TRMU-like proteins of bacteria and other homologs. TRMU is ubiquitously expressed in various tissues, but abundantly in tissues with high metabolic rates including heart, liver, kidney, and brain. Immunofluorescence analysis of human 143B cells expressing TRMU-GFP fusion protein demonstrated that the human Trmu localizes and functions in mitochondrion. Furthermore, we show that in families with the deafness-associated 12S rRNA A1491G mutation there is highly suggestive linkage and linkage disequilibrium between microsatellite markers adjacent to TRMU and the presence of deafness. These observations suggest that human TRMU may modulate the phenotypic manifestation of the deafness-associated mitochondrial 12S rRNA mutations.

A missense mutation in the PUS1 gene affecting a highly conserved amino acid has been associated with mitochondrial myopathy and sideroblastic anemia (MLASA), a rare autosomal recessive oxidative phosphorylation disorder. The PUS1 gene encodes the enzyme pseudouridine synthase 1 (Pus1p) that is known to pseudouridylate tRNAs in other species. Total RNA was isolated from lymphoblastoid cell lines established from patients, parents, unaffected siblings, and unrelated controls, and the tRNAs were assayed for the presence of pseudouridine (Psi) at the expected positions. Mitochondrial and cytoplasmic tRNAs from MLASA patients are lacking modification at sites normally modified by Pus1p, whereas tRNAs from controls, unaffected siblings, or parents all have Psi at these positions. In addition, there was no Pus1p activity in an extract made from a cell line derived from a patient with MLASA. Immunohistochemical staining of Pus1p in cell lines showed nuclear, cytoplasmic, and mitochondrial distribution of the protein, and there is no difference in staining between patients and unaffected family members. MLASA is thus associated with absent or greatly reduced tRNA pseudouridylation at specific sites, implicating this pathway in its molecular pathogenesis.

Phenotypic expression of the deafness-associated mitochondrial A1555G mutation in the 12S rRNA gene is influenced by aminoglycosides and complex inheritance of nuclear-encoded modifier genes. The position of a major nuclear modifier gene has been localized to chromosome 8p23.1, but the identification of this gene has remained elusive. Recently, we identified a second modifier gene, mitochondrial transcription factor B1 (TFB1M), involved in mitochondrial rRNA modification. In the present study, we tested three genes involved in mitochondrial tRNA or rRNA modification, and two genes associated with non-syndromic deafness, for linkage and linkage disequilibrium (LD) in 214 DNA samples from Spanish, Italian, and Arab-Israeli families with maternally inherited non-syndromic hearing loss. The multipoint non-parametric linkage analysis and transmission disequilibrium test testing were done using all families combined as well as divided based on linkage to the chromosome 8 locus and ethnicity. Two genes, MTO1 and GTPBP3, showed strongly suggestive linkage and significant LD results. Since both genes, as well as TFB1M, are involved in the process of mitochondrial RNA modification, it appears that the modification of mitochondrial RNA is an important regulatory pathway in the phenotypic expression of the deafness-associated mitochondrial A1555G mutation. This conclusion was supported by comparing linkage results of simulated genotypes with actual results for the four genes involved in mitochondrial RNA modification.

Mitochondrial myopathy and sideroblastic anemia (MLASA) is a rare, autosomal recessive oxidative phosphorylation disorder specific to skeletal muscle and bone marrow. Linkage analysis and homozygosity testing of two families with MLASA localized the candidate region to 1.2 Mb on 12q24.33. Sequence analysis of each of the six known genes in this region, as well as four putative genes with expression in bone marrow or muscle, identified a homozygous missense mutation in the pseudouridine synthase 1 gene (PUS1) in all patients with MLASA from these families. The mutation is the only amino acid coding change in these 10 genes that is not a known polymorphism, and it is not found in 934 controls. The amino acid change affects a highly conserved amino acid, and appears to be in the catalytic center of the protein, PUS1p. PUS1 is widely expressed, and quantitative expression analysis of RNAs from liver, brain, heart, bone marrow, and skeletal muscle showed elevated levels of expression in skeletal muscle and brain. We propose deficient pseudouridylation of mitochondrial tRNAs as an etiology of MLASA. Identification of the pathophysiologic pathways of the mutation in these families may shed light on the tissue specificity of oxidative phosphorylation disorders.

Mitochondrial myopathy and sideroblastic anemia (MSA) is a rare autosomal recessive disorder of oxidative phosphorylation and iron metabolism. Individuals with MSA present with weakness and anemia in late childhood and may become dependent on blood transfusions. Recently, we reported affected sibling pairs from a Jewish-Iranian kindred living in the US [Casas and Fischel-Ghodsian, 2003]. A genome scan and fine mapping of DNA from this family revealed homozygous alleles in the affected individuals, and a multipoint logarithm of the odds (lod) score of 3.3, within 2.3 mb of chromosome 12q24.33. Previously, Inbal et al. [1995: Am J Med Genet 55:372-378] described siblings with a similar clinical phenotype who lived in Israel but originated from the same Iranian town as the US family. Focused analysis of DNA from the Israeli family confirmed the presence of identical, homozygous alleles in the affected of the US and Israeli families within 1.2 mb of chromosome 12q24.33. Combined multipoint linkage analysis revealed a maximum lod score of 5.41 at the 132 cM position of chromosome 12. Therefore, in these two families of Jewish-Iranian descent, a disease gene for MSA maps to a 1.2 mb region of chromosome 12q24.33. This region contains 6 well described genes (SFRS8, MMP17, ULK1, PUS1, EP400, and GALNT9) and at least 15 additional putative transcripts. The known genes are expressed in multiple tissues and lack a function specific to mitochondria, making none an obvious candidate. The eventual identification of the disease gene in MSA is expected to provide insight into the tissue specificity and phenotypic variability of mitochondrial disease., (Copyright 2004 Wiley-Liss, Inc.)

Purpose: To examine the role of the nuclear genome in affecting the phenotypic expression of the simplest model of a mitochondrial DNA disease, maternally transmitted deafness., Methods: Linkage analysis in families with maternally inherited deafness associated with the homoplasmic A1555G mutation., Results: Significant linkage and linkage disequilibrium on chromosome 8 was identified., Conclusions: This finding represents the first identification of a modifier locus for a human mitochondrial DNA disease and supports the concept of mitochondrial DNA diseases having complex genetic inheritance. The eventual identification of this modifier gene will provide insights into the pathophysiological pathways determining the clinical expression of mitochondrial DNA diseases, an important step toward diagnostic and therapeutic interventions.

The pathophysiologic pathways and clinical expression of mitochondrial DNA (mtDNA) mutations are not well understood. This is mainly the result of the heteroplasmic nature of most pathogenic mtDNA mutations and of the absence of clinically relevant animal models with mtDNA mutations. mtDNA mutations predisposing to hearing impairment in humans are generally homoplasmic, yet some individuals with these mutations have severe hearing loss, whereas their maternal relatives with the identical mtDNA mutation have normal hearing. Epidemiologic, biochemical and genetic data indicate that nuclear genes are often the main determinants of these differences in phenotype. To identify a mouse model for maternally inherited hearing loss, we screened reciprocal backcrosses of three inbred mouse strains, A/J, NOD/LtJ and SKH2/J, with age-related hearing loss (AHL). In the (A/J x CAST/Ei) x A/J backcross, mtDNA derived from the A/J strain exerted a significant detrimental effect on hearing when compared with mtDNA from the CAST/Ei strain. This effect was not seen in the (NOD/LtJ x CAST/Ei) x NOD/LtJ and (SKH2/J x CAST/Ei) x SKH2/J backcrosses. Genotyping revealed that this effect was seen only in mice homozygous for the A/J allele at the Ahl locus on mouse chromosome 10. Sequencing of the mitochondrial genome in the three inbred strains revealed a single nucleotide insertion in the tRNA-Arg gene (mt-Tr) as the probable mediator of the mitochondrial effect. This is the first mouse model with a naturally occurring mtDNA mutation affecting a clinical phenotype, and it provides an experimental model to dissect the pathophysiologic processes connecting mtDNA mutations to hearing loss.

Maternally inherited deafness associated with the A1555G mutation in the mitochondrial 12S ribosomal RNA (rRNA) gene appears to require additional environmental or genetic changes for phenotypic expression. Aminoglycosides have been identified as one such environmental factor. In one large Arab-Israeli pedigree with congenital hearing loss in some of the family members with the A1555G mutation and with no exposure to aminoglycosides, biochemical evidence has suggested the role of nuclear modifier gene(s), but a genomewide search has indicated the absence of a single major locus having such an effect. Thus it has been concluded that the penetrance of the mitochondrial mutation appears to depend on additive effects of several nuclear genes. We have now investigated 10 multiplex Spanish and Italian families with 35 members with the A1555G mutation and sensorineural deafness. Parametric analysis of a genomewide screen again failed to identify significant evidence for linkage to a single autosomal locus. However, nonparametric analysis supported the role of the chromosomal region around marker D8S277. The combined maximized allele-sharing LOD score of 3.1 in Arab-Israeli/Spanish/Italian families represents a highly suggestive linkage result. We suggest that this region should be considered a candidate for containing the first human nuclear modifier gene for a mitochondrial DNA disorder. The locus operates in Arab-Israeli, Spanish, and Italian families, resulting in the deafness phenotype on a background of the mitochondrial A1555G mutation. No obvious candidate genes are located in this region.