Genomics Publications

Wiggs JL, Pasquale LR. Genetics of Glaucoma. Hum Mol Genet 2017;Abstract
Genetic and genomic studies, including genome-wide association studies (GWAS) have accelerated the discovery of genes contributing to glaucoma, the leading cause of irreversible blindness world-wide. Glaucoma can occur at all ages, with Mendelian inheritance typical for rare early onset disease (before age 40) and complex inheritance evident in common adult-onset forms of disease. Recent studies have suggested possible therapeutic targets for some patients with early-onset glaucoma based on the molecular and cellular events caused by MYOC, OPTN and TBK1 mutations. Diagnostic genetic tests using early-onset glaucoma genes are also proving useful for pre-symptomatic disease detection and genetic counseling. Recent GWAS completed for three types of common adult-onset glaucoma have identified novel loci for POAG (primary-open-angle glaucoma) (ABCA1, AFAP1, GMDS, PMM2, TGFBR3, FNDC3B, ARHGEF12, GAS7, FOXC1, ATXN2, TXNRD2); PACG (primary angle-closure glaucoma (EPDR1, CHAT, GLIS3, FERMT2, DPM2-FAM102); and exfoliation syndrome (XFS) glaucoma (CACNA1A). In total sixteen genomic regions have been associated with POAG (including the normal tension glaucoma (NTG) subgroup), 8 with PACG and 2 with XFS. These studies are defining important biological pathways and processes that contribute to disease pathogenesis.
Wang T, Huang T, Kang JH, Zheng Y, Jensen MK, Wiggs JL, Pasquale LR, Fuchs CS, Campos H, Rimm EB, Willett WC, Hu FB, Qi L. Habitual coffee consumption and genetic predisposition to obesity: gene-diet interaction analyses in three US prospective studies. BMC Med 2017;15(1):97.Abstract
BACKGROUND: Whether habitual coffee consumption interacts with the genetic predisposition to obesity in relation to body mass index (BMI) and obesity is unknown. METHODS: We analyzed the interactions between genetic predisposition and habitual coffee consumption in relation to BMI and obesity risk in 5116 men from the Health Professionals Follow-up Study (HPFS), in 9841 women from the Nurses' Health Study (NHS), and in 5648 women from the Women's Health Initiative (WHI). The genetic risk score was calculated based on 77 BMI-associated loci. Coffee consumption was examined prospectively in relation to BMI. RESULTS: The genetic association with BMI was attenuated among participants with higher consumption of coffee than among those with lower consumption in the HPFS (P interaction  = 0.023) and NHS (P interaction  = 0.039); similar results were replicated in the WHI (P interaction  = 0.044). In the combined data of all cohorts, differences in BMI per increment of 10-risk allele were 1.38 (standard error (SE), 0.28), 1.02 (SE, 0.10), and 0.95 (SE, 0.12) kg/m(2) for coffee consumption of < 1, 1-3 and > 3 cup(s)/day, respectively (P interaction  < 0.001). Such interaction was partly due to slightly higher BMI with higher coffee consumption among participants at lower genetic risk and slightly lower BMI with higher coffee consumption among those at higher genetic risk. Each increment of 10-risk allele was associated with 78% (95% confidence interval (CI), 59-99%), 48% (95% CI, 36-62%), and 43% (95% CI, 28-59%) increased risk for obesity across these subgroups of coffee consumption (P interaction  = 0.008). From another perspective, differences in BMI per increment of 1 cup/day coffee consumption were 0.02 (SE, 0.09), -0.02 (SE, 0.04), and -0.14 (SE, 0.04) kg/m(2) across tertiles of the genetic risk score. CONCLUSIONS: Higher coffee consumption might attenuate the genetic associations with BMI and obesity risk, and individuals with greater genetic predisposition to obesity appeared to have lower BMI associated with higher coffee consumption.
Lonfat N, Cepko C. Epigenomics of Retinal Development in Mice and Humans. Neuron 2017;94(3):420-423.Abstract
In this issue of Neuron, Aldiri et al. (2017) present an analysis of epigenetic changes during retinal development, and use these data to probe reprogramming of retinal iPSC cells, as well as the origin of retinoblastoma cells.
Aung T, Ozaki M, Lee MC, Schlötzer-Schrehardt U, Thorleifsson G, Mizoguchi T, Igo RP, Haripriya A, Williams SE, Astakhov YS, Orr AC, Burdon KP, Nakano S, Mori K, Abu-Amero K, Hauser M, Li Z, Prakadeeswari G, Bailey JCN, Cherecheanu AP, Kang JH, Nelson S, Hayashi K, Manabe S-I, Kazama S, Zarnowski T, Inoue K, Irkec M, Coca-Prados M, Sugiyama K, Järvelä I, Schlottmann P, Lerner FS, Lamari H, Nilgün Y, Bikbov M, Park KH, Cha SC, Yamashiro K, Zenteno JC, Jonas JB, Kumar RS, Perera SA, Chan ASY, Kobakhidze N, George R, Vijaya L, Do T, Edward DP, de Juan Marcos L, Pakravan M, Moghimi S, Ideta R, Bach-Holm D, Kappelgaard P, Wirostko B, Thomas S, Gaston D, Bedard K, Greer WL, Yang Z, Chen X, Huang L, Sang J, Jia H, Jia L, Qiao C, Zhang H, Liu X, Zhao B, Wang Y-X, Xu L, Leruez S, Reynier P, Chichua G, Tabagari S, Uebe S, Zenkel M, Berner D, Mossböck G, Weisschuh N, Hoja U, Welge-Luessen U-C, Mardin C, Founti P, Chatzikyriakidou A, Pappas T, Anastasopoulos E, Lambropoulos A, Ghosh A, Shetty R, Porporato N, Saravanan V, Venkatesh R, Shivkumar C, Kalpana N, Sarangapani S, Kanavi MR, Beni AN, Yazdani S, Lashay A, Naderifar H, Khatibi N, Fea A, Lavia C, Dallorto L, Rolle T, Frezzotti P, Paoli D, Salvi E, Manunta P, Mori Y, Miyata K, Higashide T, Chihara E, Ishiko S, Yoshida A, Yanagi M, Kiuchi Y, Ohashi T, Sakurai T, Sugimoto T, Chuman H, Aihara M, Inatani M, Miyake M, Gotoh N, Matsuda F, Yoshimura N, Ikeda Y, Ueno M, Sotozono C, Jeoung JW, Sagong M, Park KH, Ahn J, Cruz-Aguilar M, Ezzouhairi SM, Rafei A, Chong YF, Ng XY, Goh SR, Chen Y, Yong VHK, Khan MI, Olawoye OO, Ashaye AO, Ugbede I, Onakoya A, Kizor-Akaraiwe N, Teekhasaenee C, Suwan Y, Supakontanasan W, Okeke S, Uche NJ, Asimadu I, Ayub H, Akhtar F, Kosior-Jarecka E, Lukasik U, Lischinsky I, Castro V, Grossmann RP, Megevand GS, Roy S, Dervan E, Silke E, Rao A, Sahay P, Fornero P, Cuello O, Sivori D, Zompa T, Mills RA, Souzeau E, Mitchell P, Wang JJ, Hewitt AW, Coote M, Crowston JG, Astakhov SY, Akopov EL, Emelyanov A, Vysochinskaya V, Kazakbaeva G, Fayzrakhmanov R, Al-Obeidan SA, Owaidhah O, Aljasim LA, Chowbay B, Foo JN, Soh RQ, Sim KS, Xie Z, Cheong AWO, Mok SQ, Soo HM, Chen XY, Peh SQ, Heng KK, Husain R, Ho S-L, Hillmer AM, Cheng C-Y, Escudero-Domínguez FA, González-Sarmiento R, Martinon-Torres F, Salas A, Pathanapitoon K, Hansapinyo L, Wanichwecharugruang B, Kitnarong N, Sakuntabhai A, Nguyn HX, Nguyn GTT, Nguyn TV, Zenz W, Binder A, Klobassa DS, Hibberd ML, Davila S, Herms S, Nöthen MM, Moebus S, Rautenbach RM, Ziskind A, Carmichael TR, Ramsay M, Álvarez L, García M, González-Iglesias H, Rodríguez-Calvo PP, Cueto LF-V, Oguz Ç, Tamcelik N, Atalay E, Batu B, Aktas D, Kasım B, Wilson RM, Coleman AL, Liu Y, Challa P, Herndon L, Kuchtey RW, Kuchtey J, Curtin K, Chaya CJ, Crandall A, Zangwill LM, Wong TY, Nakano M, Kinoshita S, den Hollander AI, Vesti E, Fingert JH, Lee RK, Sit AJ, Shingleton BJ, Wang N, Cusi D, Qamar R, Kraft P, Pericak-Vance MA, Raychaudhuri S, Heegaard S, Kivelä T, Reis A, Kruse FE, Weinreb RN, Pasquale LR, Haines JL, Thorsteinsdottir U, Jonasson F, Allingham RR, Milea D, Ritch R, Kubota T, Tashiro K, Vithana EN, Micheal S, Topouzis F, Craig JE, Dubina M, Sundaresan P, Stefansson K, Wiggs JL, Pasutto F, Khor CC. Genetic association study of exfoliation syndrome identifies a protective rare variant at LOXL1 and five new susceptibility loci. Nat Genet 2017;Abstract
Exfoliation syndrome (XFS) is the most common known risk factor for secondary glaucoma and a major cause of blindness worldwide. Variants in two genes, LOXL1 and CACNA1A, have previously been associated with XFS. To further elucidate the genetic basis of XFS, we collected a global sample of XFS cases to refine the association at LOXL1, which previously showed inconsistent results across populations, and to identify new variants associated with XFS. We identified a rare protective allele at LOXL1 (p.Phe407, odds ratio (OR) = 25, P = 2.9 × 10(-14)) through deep resequencing of XFS cases and controls from nine countries. A genome-wide association study (GWAS) of XFS cases and controls from 24 countries followed by replication in 18 countries identified seven genome-wide significant loci (P < 5 × 10(-8)). We identified association signals at 13q12 (POMP), 11q23.3 (TMEM136), 6p21 (AGPAT1), 3p24 (RBMS3) and 5q23 (near SEMA6A). These findings provide biological insights into the pathology of XFS and highlight a potential role for naturally occurring rare LOXL1 variants in disease biology.
Ramke M, Lee JY, Dyer DW, Seto D, Rajaiya J, Chodosh J. The 5'UTR in human adenoviruses: leader diversity in late gene expression. Sci Rep 2017;7(1):618.Abstract

Human adenoviruses (HAdVs) shut down host cellular cap-dependent mRNA translation while initiating the translation of viral late mRNAs in a cap-independent manner. HAdV 5' untranslated regions (5'UTRs) are crucial for cap-independent initiation, and influence mRNA localization and stability. However, HAdV translational regulation remains relatively uncharacterized. The HAdV tripartite leader (TPL), composed of three introns (TPL 1-3), is critical to the translation of HAdV late mRNA. Herein, we annotated and analyzed 72 HAdV genotypes for the HAdV TPL and another previously described leader, the i-leader. Using HAdV species D, type 37 (HAdV-D37), we show by reverse transcription PCR and Sanger sequencing that mRNAs of the HAdV-D37 E3 transcription unit are spliced to the TPL. We also identified a polycistronic mRNA for RID-α and RID-β. Analysis of the i-leader revealed a potential open reading frame within the leader sequence and the termination of this potential protein in TPL3. A potential new leader embedded within the E3 region was also detected and tentatively named the j-leader. These results suggest an underappreciated complexity of post-transcriptional regulation, and the importance of HAdV 5'UTRs for precisely coordinated viral protein expression along the path from genotype to phenotype.

Suzuki J, Hashimoto K, Xiao R, Vandenberghe LH, Liberman CM. Cochlear gene therapy with ancestral AAV in adult mice: complete transduction of inner hair cells without cochlear dysfunction. Sci Rep 2017;7:45524.Abstract

The use of viral vectors for inner ear gene therapy is receiving increased attention for treatment of genetic hearing disorders. Most animal studies to date have injected viral suspensions into neonatal ears, via the round window membrane. Achieving transduction of hair cells, or sensory neurons, throughout the cochlea has proven difficult, and no studies have been able to efficiently transduce sensory cells in adult ears while maintaining normal cochlear function. Here, we show, for the first time, successful transduction of all inner hair cells and the majority of outer hair cells in an adult cochlea via virus injection into the posterior semicircular canal. We used a "designer" AAV, AAV2/Anc80L65, in which the main capsid proteins approximate the ancestral sequence state of AAV1, 2, 8, and 9. Our injections also transduced ~10% of spiral ganglion cells and a much larger fraction of their satellite cells. In the vestibular sensory epithelia, the virus transduced large numbers of hair cells and virtually all the supporting cells, along with close to half of the vestibular ganglion cells. We conclude that this viral vector and this delivery route hold great promise for gene therapy applications in both cochlear and vestibular sense organs.

Saavedra JT, Schwartzman JA, Gilmore MS. Mapping Transposon Insertions in Bacterial Genomes by Arbitrarily Primed PCR. Curr Protoc Mol Biol 2017;118:15.15.1-15.15.15.Abstract

Transposons can be used to easily generate and label the location of mutations throughout bacterial and other genomes. Transposon insertion mutants may be screened for a phenotype as individual isolates, or by selection applied to a pool of thousands of mutants. Identifying the location of a transposon insertion is critical for connecting phenotype to the genetic lesion. In this unit, we present an easy and detailed approach for mapping transposon insertion sites using arbitrarily-primed PCR (AP-PCR). Two rounds of PCR are used to (1) amplify DNA spanning the transposon insertion junction, and (2) increase the specific yield of transposon insertion junction fragments for sequence analysis. The resulting sequence is mapped to a bacterial genome to identify the site of transposon insertion. In this protocol, AP-PCR as it is routinely used to map sites of transposon insertion within Staphylococcus aureus, is used to illustrate the principle. Guidelines are provided for adapting this protocol for mapping insertions in other bacterial genomes. Mapping transposon insertions using this method is typically achieved in 2 to 3 days if starting from a culture of the transposon insertion mutant. © 2017 by John Wiley & Sons, Inc.

Ung C, Sanchez AV, Shen L, Davoudi S, Ahmadi T, Navarro-Gomez D, Chen CJ, Hancock H, Penman A, Hoadley S, Consugar M, Restrepo C, Shah VA, Arboleda-Velasquez JF, Sobrin L, Gai X, Kim LA. Whole exome sequencing identification of novel candidate genes in patients with proliferative diabetic retinopathy. Vision Res 2017;Abstract
Rare or novel gene variants in patients with proliferative diabetic retinopathy may contribute to disease development. We performed whole exome sequencing (WES) on patients at the phenotypic extremes of diabetic retinal complications: 57 patients diagnosed with proliferative diabetic retinopathy (PDR) as cases and 13 patients with no diabetic retinopathy despite at least 10years of type 2 diabetes as controls. Thirty-one out of the 57 cases and all 13 controls were from the African American Proliferative Diabetic Retinopathy Study (AA). The rest of the cases were of mixed ethnicities (ME). WES identified 721 candidate genes with rare or novel non-synonymous variants found in at least one case with PDR and not present in any controls. After filtering for genes with null alleles in greater than two cases, 28 candidate genes were identified in our ME cases and 16 genes were identified in our AA cases. Our analysis showed rare and novel variants within these genes that could contribute to the development of PDR, including rare non-synonymous variants in FAM132A, SLC5A9, ZNF600, and TMEM217. We also found previously unidentified variants in VEGFB and APOB. We found that VEGFB, VPS13B, PHF21A, NAT1, ZNF600, PKHD1L1 expression was reduced in human retinal endothelial cells (HRECs) cultured under high glucose conditions. In an exome sequence analysis of patients with PDR, we identified variants in genes that could contribute to pathogenesis. Six of these genes were further validated and found to have reduced expression in HRECs under high glucose conditions, suggestive of an important role in the development of PDR.
Pacouret S, Bouzelha M, Shelke R, Andres-Mateos E, Xiao R, Maurer A, Mevel M, Turunen H, Barungi T, Penaud-Budloo M, Broucque F, Blouin V, Moullier P, Ayuso E, Vandenberghe LH. AAV-ID: A Rapid and Robust Assay for Batch-to-Batch Consistency Evaluation of AAV Preparations. Mol Ther 2017;Abstract

Adeno-associated virus (AAV) vectors are promising clinical candidates for therapeutic gene transfer, and a number of AAV-based drugs may emerge on the market over the coming years. To insure the consistency in efficacy and safety of any drug vial that reaches the patient, regulatory agencies require extensive characterization of the final product. Identity is a key characteristic of a therapeutic product, as it ensures its proper labeling and batch-to-batch consistency. Currently, there is no facile, fast, and robust characterization assay enabling to probe the identity of AAV products at the protein level. Here, we investigated whether the thermostability of AAV particles could inform us on the composition of vector preparations. AAV-ID, an assay based on differential scanning fluorimetry (DSF), was evaluated in two AAV research laboratories for specificity, sensitivity, and reproducibility, for six different serotypes (AAV1, 2, 5, 6.2, 8, and 9), using 67 randomly selected AAV preparations. In addition to enabling discrimination of AAV serotypes based on their melting temperatures, the obtained fluorescent fingerprints also provided information on sample homogeneity, particle concentration, and buffer composition. Our data support the use of AAV-ID as a reproducible, fast, and low-cost method to ensure batch-to-batch consistency in manufacturing facilities and academic laboratories.

Lindström S, Loomis S, Turman C, Huang H, Huang J, Aschard H, Chan AT, Choi H, Cornelis M, Curhan G, De Vivo I, Eliassen HA, Fuchs C, Gaziano M, Hankinson SE, Hu F, Jensen M, Kang JH, Kabrhel C, Liang L, Pasquale LR, Rimm E, Stampfer MJ, Tamimi RM, Tworoger SS, Wiggs JL, Hunter DJ, Kraft P. A comprehensive survey of genetic variation in 20,691 subjects from four large cohorts. PLoS One 2017;12(3):e0173997.Abstract

The Nurses' Health Study (NHS), Nurses' Health Study II (NHSII), Health Professionals Follow Up Study (HPFS) and the Physicians Health Study (PHS) have collected detailed longitudinal data on multiple exposures and traits for approximately 310,000 study participants over the last 35 years. Over 160,000 study participants across the cohorts have donated a DNA sample and to date, 20,691 subjects have been genotyped as part of genome-wide association studies (GWAS) of twelve primary outcomes. However, these studies utilized six different GWAS arrays making it difficult to conduct analyses of secondary phenotypes or share controls across studies. To allow for secondary analyses of these data, we have created three new datasets merged by platform family and performed imputation using a common reference panel, the 1,000 Genomes Phase I release. Here, we describe the methodology behind the data merging and imputation and present imputation quality statistics and association results from two GWAS of secondary phenotypes (body mass index (BMI) and venous thromboembolism (VTE)). We observed the strongest BMI association for the FTO SNP rs55872725 (β = 0.45, p = 3.48x10-22), and using a significance level of p = 0.05, we replicated 19 out of 32 known BMI SNPs. For VTE, we observed the strongest association for the rs2040445 SNP (OR = 2.17, 95% CI: 1.79-2.63, p = 2.70x10-15), located downstream of F5 and also observed significant associations for the known ABO and F11 regions. This pooled resource can be used to maximize power in GWAS of phenotypes collected across the cohorts and for studying gene-environment interactions as well as rare phenotypes and genotypes.

Zeng Z, Zhang J, Jing S, Cheng Z, Bofill-Mas S, Maluquer de Motes C, Hundesa A, Girones R, Seto D, Zhang Q. Genome Sequence of a Cynomolgus Macaque Adenovirus (CynAdV-1) Isolate from a Primate Colony in the United Kingdom. Genome Announc 2016;4(6)Abstract

The genome sequence of a simian adenovirus from a cynomolgus macaque, denoted CynAdV-1, is presented here. Phylogenetic analysis supports CynAdV-1 in an independent clade, comprising a new simian adenovirus (SAdV) species. These genome data are critical for understanding the evolution and relationships of primate adenoviruses, including zoonosis and emergent human pathogens.

Bujakowska KM, Fernandez-Godino R, Place E, Consugar M, Navarro-Gomez D, White J, Bedoukian EC, Zhu X, Xie HM, Gai X, Leroy BP, Pierce EA. Copy-number variation is an important contributor to the genetic causality of inherited retinal degenerations. Genet Med 2016;Abstract

PURPOSE: Despite substantial progress in sequencing, current strategies can genetically solve only approximately 55-60% of inherited retinal degeneration (IRD) cases. This can be partially attributed to elusive mutations in the known IRD genes, which are not easily identified by the targeted next-generation sequencing (NGS) or Sanger sequencing approaches. We hypothesized that copy-number variations (CNVs) are a major contributor to the elusive genetic causality of IRDs. METHODS: Twenty-eight cases previously unsolved with a targeted NGS were investigated with whole-genome single-nucleotide polymorphism (SNP) and comparative genomic hybridization (CGH) arrays. RESULTS: Deletions in the IRD genes were detected in 5 of 28 families, including a de novo deletion. We suggest that the de novo deletion occurred through nonallelic homologous recombination (NAHR) and we constructed a genomic map of NAHR-prone regions with overlapping IRD genes. In this article, we also report an unusual case of recessive retinitis pigmentosa due to compound heterozygous mutations in SNRNP200, a gene that is typically associated with the dominant form of this disease. CONCLUSIONS: CNV mapping substantially increased the genetic diagnostic rate of IRDs, detecting genetic causality in 18% of previously unsolved cases. Extending the search to other structural variations will probably demonstrate an even higher contribution to genetic causality of IRDs.Genet Med advance online publication 13 October 2016Genetics in Medicine (2016); doi:10.1038/gim.2016.158.

Lamont RE, Tan W-H, Innes MA, Parboosingh JS, Schneidman-Duhovny D, Rajkovic A, Pappas J, Altschwager P, DeWard S, Fulton A, Gray KJ, Krall M, Mehta L, Rodan LH, Saller DN, Steele D, Stein D, Yatsenko SA, Bernier FP, Slavotinek AM. Expansion of phenotype and genotypic data in CRB2-related syndrome. Eur J Hum Genet 2016;24(10):1436-44.Abstract

Sequence variants in CRB2 cause a syndrome with greatly elevated maternal serum alpha-fetoprotein and amniotic fluid alpha-fetoprotein levels, cerebral ventriculomegaly and renal findings similar to Finnish congenital nephrosis. All reported patients have been homozygotes or compound heterozygotes for sequence variants in the Crumbs, Drosophila, Homolog of, 2 (CRB2) genes. Variants affecting CRB2 function have also been identified in four families with steroid resistant nephrotic syndrome, but without any other known systemic findings. We ascertained five, previously unreported individuals with biallelic variants in CRB2 that were predicted to affect function. We compiled the clinical features of reported cases and reviewed available literature for cases with features suggestive of CRB2-related syndrome in order to better understand the phenotypic and genotypic manifestations. Phenotypic analyses showed that ventriculomegaly was a common clinical manifestation (9/11 confirmed cases), in contrast to the original reports, in which patients were ascertained due to renal disease. Two children had minor eye findings and one was diagnosed with a B-cell lymphoma. Further genetic analysis identified one family with two affected siblings who were both heterozygous for a variant in NPHS2 predicted to affect function and separate families with sequence variants in NPHS4 and BBS7 in addition to the CRB2 variants. Our report expands the clinical phenotype of CRB2-related syndrome and establishes ventriculomegaly and hydrocephalus as frequent manifestations. We found additional sequence variants in genes involved in kidney development and ciliopathies in patients with CRB2-related syndrome, suggesting that these variants may modify the phenotype.

Gao X, Nannini DR, Corrao K, Torres M, Chen Y-DI, Fan BJ, Wiggs JL, Wiggs JL, Taylor KD, James Gauderman W, Rotter JI, Varma R. Genome-Wide Association Study Identifies WNT7B as a Novel Locus for Central Corneal Thickness in Latinos. Hum Mol Genet 2016;Abstract

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

Arno G, Hull S, Carss K, Dev-Borman A, Chakarova C, Bujakowska K, van den Born I, Robson AG, Holder GE, Michaelides M, Cremers FPM, Pierce E, Raymond LF, Moore AT, Webster AR. Reevaluation of the Retinal Dystrophy Due to Recessive Alleles of RGR With the Discovery of a Cis-Acting Mutation in CDHR1. Invest Ophthalmol Vis Sci 2016;57(11):4806-13.Abstract

PURPOSE: Mutation of RGR, encoding retinal G-protein coupled receptor was originally reported in association with retinal dystrophy in 1999. A single convincing recessive variant segregated perfectly in one family of five affected and two unaffected siblings. At least one further individual, homozygous for the same variant has since been reported. The aim of this report was to reevaluate the findings in consideration of data from a whole genome sequencing (WGS) study of a large cohort of retinal dystrophy families. METHODS: Whole genome sequencing was performed on 599 unrelated probands with inherited retinal disease. Detailed phenotyping was performed, including clinical evaluation, electroretinography, fundus photography, fundus autofluorescence imaging (FAF) and spectral-domain optical coherence tomography (OCT). RESULTS: Overall we confirmed that affected individuals from six unrelated families were homozygous for both the reported RGR p.Ser66Arg variant and a nearby frameshifting deletion in CDHR1 (p.Ile841Serfs119*). All had generalized rod and cone dysfunction with severe macular involvement. An additional proband was heterozygous for the same CDHR1/RGR haplotype but also carried a second null CDHR1 mutation on a different haplotype. A comparison of the clinical presentation of the probands reported here with other CDHR1-related retinopathy patients shows the phenotypes to be similar in presentation, severity, and rod/cone involvement. CONCLUSIONS: These data suggest that the recessive retinal disorder previously reported to be due to homozygous mutation in RGR is, at least in part, due to variants in CDHR1 and that the true consequences of RGR knock-out on human retinal structure and function are yet to be determined.