Genomics Publications

Mychaleckyj J, Valo E, Ichimura T, Ahluwalia T, Dina C, Miller R, Shabalin I, Gyorgy B, Cao J, Onengut-Gumuscu S, Satake E, Smiles A, Haukka J, Tregouet D-A, Costacou T, O'Neil K, Paterson A, Forsblom C, Keenan H, Pezzolesi M, Pragnell M, Galecki A, Rich S, Sandholm N, Klein R, Klein B, Susztak K, Orchard T, Korstanje R, King G, Hadjadj S, Rossing P, Bonventre J, Groop P-H, Warram J, Krolewski A. Association of Coding Variants in Hydroxysteroid 17-beta Dehydrogenase 14 (HSD17B14) with Reduced Progression to End Stage Kidney Disease in Type 1 Diabetes. J Am Soc Nephrol 2021;Abstract
BACKGROUND: Rare variants in gene coding regions likely have a greater impact on disease-related phenotypes than common variants through disruption of their encoded protein. We searched for rare variants associated with onset of end stage kidney disease (ESKD) in individuals with type 1 diabetes at advanced kidney disease stage. METHODS: Gene-based exome array analysis of 15,449 genes in 5 large incidence cohorts of individuals with type 1 diabetes and proteinuria were analyzed for survival time-to-ESKD, testing the top gene in a 6th cohort (N=2,372/1,115 events all cohorts) and replicating in two retrospective case-control studies (N=1,072 cases, 752 controls). Deep resequencing of the top associated gene in 5 cohorts confirmed the findings. We performed immunohistochemistry and gene expression experiments in human control and diseased cells, and in mouse ischemia reperfusion and aristolochic acid nephropathy models. RESULTS: Protein coding variants in the hydroxysteroid 17-beta dehydrogenase 14 gene (HSD17B14), predicted to affect protein structure, had a net protective effect against development of ESKD at exome-wide significance (N=4,196; p-value=3.3x10-7). The HSD17B14 gene and encoded enzyme were robustly expressed in healthy human kidney, maximally in proximal tubular cells. Paradoxically, gene and protein expression were attenuated in human diabetic proximal tubules and in mouse kidney injury models. Expressed HSD17B14 gene and protein levels remained low without recovery after 21 days in a murine ischemic reperfusion injury model. Decreased gene expression was found in other chronic kidney disease-associated renal pathologies. CONCLUSIONS: HSD17B14 gene is mechanistically involved in diabetic kidney disease. The encoded sex steroid enzyme is a druggable target, potentially opening a new avenue for therapeutic development.
van der Heide C, Goar W, Meyer KJ, Alward WLM, Boese EA, Sears NC, Roos BR, Kwon YH, DeLuca AP, Siggs OM, Gonzaga-Jauregui C, Sheffield VC, Wang K, Stone EM, Mullins RF, Anderson MG, Fan BJ, Ritch R, Craig JE, Wiggs JL, Scheetz TE, Fingert JH. Exome-based investigation of the genetic basis of human pigmentary glaucoma. BMC Genomics 2021;22(1):477.Abstract
BACKGROUND: Glaucoma is a leading cause of visual disability and blindness. Release of iris pigment within the eye, pigment dispersion syndrome (PDS), can lead to one type of glaucoma known as pigmentary glaucoma. PDS has a genetic component, however, the genes involved with this condition are largely unknown. We sought to discover genes that cause PDS by testing cohorts of patients and controls for mutations using a tiered analysis of exome data. RESULTS: Our primary analysis evaluated melanosome-related genes that cause dispersion of iris pigment in mice (TYRP1, GPNMB, LYST, DCT, and MITF). We identified rare mutations, but they were not statistically enriched in PDS patients. Our secondary analyses examined PMEL (previously linked with PDS), MRAP, and 19 other genes. Four MRAP mutations were identified in PDS cases but not in controls (p = 0.016). Immunohistochemical analysis of human donor eyes revealed abundant MRAP protein in the iris, the source of pigment in PDS. However, analysis of MRAP in additional cohorts (415 cases and 1645 controls) did not support an association with PDS. We also did not confirm a link between PMEL and PDS in our cohorts due to lack of reported mutations and similar frequency of the variants in PDS patients as in control subjects. CONCLUSIONS: We did not detect a statistical enrichment of mutations in melanosome-related genes in human PDS patients and we found conflicting data about the likely pathogenicity of MRAP mutations. PDS may have a complex genetic basis that is not easily unraveled with exome analyses.
Sangermano R, Deitch I, Peter VG, Ba-Abbad R, Place EM, Zampaglione E, Wagner NE, Fulton AB, Coutinho-Santos L, Rosin B, Dunet V, AlTalbishi A'a, Banin E, Sousa AB, Neves M, Larson A, Quinodoz M, Michaelides M, Ben-Yosef T, Pierce EA, Rivolta C, Webster AR, Arno G, Sharon D, Huckfeldt RM, Bujakowska KM. Broadening INPP5E phenotypic spectrum: detection of rare variants in syndromic and non-syndromic IRD. NPJ Genom Med 2021;6(1):53.Abstract
Pathogenic variants in INPP5E cause Joubert syndrome (JBTS), a ciliopathy with retinal involvement. However, despite sporadic cases in large cohort sequencing studies, a clear association with non-syndromic inherited retinal degenerations (IRDs) has not been made. We validate this association by reporting 16 non-syndromic IRD patients from ten families with bi-allelic mutations in INPP5E. Additional two patients showed early onset IRD with limited JBTS features. Detailed phenotypic description for all probands is presented. We report 14 rare INPP5E variants, 12 of which have not been reported in previous studies. We present tertiary protein modeling and analyze all INPP5E variants for deleteriousness and phenotypic correlation. We observe that the combined impact of INPP5E variants in JBTS and non-syndromic IRD patients does not reveal a clear genotype-phenotype correlation, suggesting the involvement of genetic modifiers. Our study cements the wide phenotypic spectrum of INPP5E disease, adding proof that sequence defects in this gene can lead to early-onset non-syndromic IRD.
Rossin EJ, Sobrin L, Kim LA. Single-cell RNA sequencing: An overview for the ophthalmologist. Semin Ophthalmol 2021;36(4):191-197.Abstract
Understanding the molecular composition of pathogenic tissues is a critical step in understanding the pathophysiology of disease and designing therapeutics. First described in 2009, single cell RNA sequencing (scRNAseq) is a methodology whereby thousands of cells are simultaneously isolated into individual micro-environments that can be altered experimentally and the genome-wide RNA expression of each cell is captured. It has undergone significant technological improvement over the last decade and gained tremendous popularity. scRNAseq is an improvement over prior pooled RNA analyses which cannot identify the cellular composition and heterogeneity of a tissue of interest. This new approach offers new opportunity for new discovery, as tissue samples can now be sub-categorized into groups of cell types based on genome-wide gene expression in an unbiased fashion. As ophthalmologists, we are uniquely positioned to obtain pathologic samples from the eye for further study. ScRNAseq has already been applied in ophthalmology to characterize retinal tissue, and it may offer the key to understanding various pathological processes in the future.
Imamura M, Takahashi A, Matsunami M, Horikoshi M, Iwata M, Araki S-I, Toyoda M, Susarla G, Ahn J, Park KH, Kong J, Moon S, Sobrin L, and (iDRAGON) IDRGCON, Yamauchi T, Tobe K, Maegawa H, Kadowaki T, Maeda S. Genome-wide association studies identify two novel loci conferring susceptibility to diabetic retinopathy in Japanese patients with type 2 diabetes. Hum Mol Genet 2021;30(8):716-726.Abstract
Several reports have suggested that genetic susceptibility contributes to the development and progression of diabetic retinopathy. We aimed to identify genetic loci that confer susceptibility to diabetic retinopathy in Japanese patients with type 2 diabetes. We analysed 5 790 508 single nucleotide polymorphisms (SNPs) in 8880 Japanese patients with type 2 diabetes, 4839 retinopathy cases and 4041 controls, as well as 2217 independent Japanese patients with type 2 diabetes, 693 retinopathy cases and 1524 controls. The results of these two genome-wide association studies (GWAS) were combined with an inverse variance meta-analysis (Stage-1), followed by de novo genotyping for the candidate SNP loci (P < 1.0 × 10-4) in an independent case-control study (Stage-2, 2260 cases and 723 controls). After combining the association data (Stages 1 and 2) using meta-analysis, the associations of two loci reached a genome-wide significance level: rs12630354 near STT3B on chromosome 3, P = 1.62 × 10-9, odds ratio (OR) = 1.17, 95% confidence interval (CI) 1.11-1.23, and rs140508424 within PALM2 on chromosome 9, P = 4.19 × 10-8, OR = 1.61, 95% CI 1.36-1.91. However, the association of these two loci was not replicated in Korean, European or African American populations. Gene-based analysis using Stage-1 GWAS data identified a gene-level association of EHD3 with susceptibility to diabetic retinopathy (P = 2.17 × 10-6). In conclusion, we identified two novel SNP loci, STT3B and PALM2, and a novel gene, EHD3, that confers susceptibility to diabetic retinopathy; however, further replication studies are required to validate these associations.
Peng C, Wang Y, Ji L, Kuang L, Yu Z, Li H, Zhang J, Zhao J. LncRNA-MALAT1/miRNA-204-5p/Smad4 Axis Regulates Epithelial-Mesenchymal Transition, Proliferation and Migration of Lens Epithelial Cells. Curr Eye Res 2021;46(8):1137-1147.Abstract
MATERIALS AND METHODS: LECs were cultured and induced with TGF-β2 (10 ng/mL). SiRNA against MALAT1 (Si-MALAT1) was transfected into LECs to knockdown the expression of MALAT1. To overexpress or knockdown miR-204-5p, miR-204-5p mimics (miR-204-5p mimics) and anti-miR-204-5p (miR-204-5p inhibitor) were transfected into LECs. We used RNA FISH to identify the location of MALAT1. RNA levels of MALAT1 and miR-204-5p were analyzed by RT-qPCR. Additionally, target protein levels of Smad4, epithelial differentiation and mesenchymal markers were analyzed with Western blot. We employed EdU Labeling to measured cell proliferation and performed Transwell Assay to analyze the cell migration. Dual-luciferase reporter assays in LECs were conducted to verify whether miRNA-204-5p was negatively regulated by MALAT1 and Smad4 was a direct target of miR-204-5p. RESULTS: The expression of MALAT1 was upregulated in PCO specimens. MALAT1 was overexpressed in TGF-β2 induced LECs, and the knockdown of MALAT1 could attenuate TGF-β2 induced EMT. Besides, the upregulation of MALAT1 was correlated with the downregulation of miR-204-5p and upregulation of Smad4. Importantly, MALAT1 was revealed to be located in the cytoplasm of LECs. Furthermore, luciferase reporter assays confirmed that MALAT1 could negatively regulate the expression of miR-204-5p and then regulate its direct target Smad4. Finally, the knockdown of MALAT1 could inhibit the EMT, proliferation, and migration of LECs; however, those can be reversed by anti-miR-204-5p. CONCLUSIONS: Our findings reveal that MALAT1 may regulate EMT, proliferation, and migration of LECs as a ceRNA by "sponging" miR-204-5p and targeting Smad4, and serve as a promising therapeutic target in preventing PCO.
Cogné B, Latypova X, Senaratne LDS, Martin L, Koboldt DC, Kellaris G, Fievet L, Le Meur G, Caldari D, Debray D, Nizon M, Frengen E, Bowne SJ, Consortium L99, Cadena EL, Daiger SP, Bujakowska KM, Pierce EA, Gorin M, Katsanis N, Bézieau S, Petersen-Jones SM, Occelli LM, Lyons LA, Legeai-Mallet L, Sullivan LS, Davis EE, Isidor B. Mutations in the Kinesin-2 Motor KIF3B Cause an Autosomal-Dominant Ciliopathy. Am J Hum Genet 2020;106(6):893-904.Abstract
Kinesin-2 enables ciliary assembly and maintenance as an anterograde intraflagellar transport (IFT) motor. Molecular motor activity is driven by a heterotrimeric complex comprised of KIF3A and KIF3B or KIF3C plus one non-motor subunit, KIFAP3. Using exome sequencing, we identified heterozygous KIF3B variants in two unrelated families with hallmark ciliopathy phenotypes. In the first family, the proband presents with hepatic fibrosis, retinitis pigmentosa, and postaxial polydactyly; he harbors a de novo c.748G>C (p.Glu250Gln) variant affecting the kinesin motor domain encoded by KIF3B. The second family is a six-generation pedigree affected predominantly by retinitis pigmentosa. Affected individuals carry a heterozygous c.1568T>C (p.Leu523Pro) KIF3B variant segregating in an autosomal-dominant pattern. We observed a significant increase in primary cilia length in vitro in the context of either of the two mutations while variant KIF3B proteins retained stability indistinguishable from wild type. Furthermore, we tested the effects of KIF3B mutant mRNA expression in the developing zebrafish retina. In the presence of either missense variant, rhodopsin was sequestered to the photoreceptor rod inner segment layer with a concomitant increase in photoreceptor cilia length. Notably, impaired rhodopsin trafficking is also characteristic of recessive KIF3B models as exemplified by an early-onset, autosomal-recessive, progressive retinal degeneration in Bengal cats; we identified a c.1000G>A (p.Ala334Thr) KIF3B variant by genome-wide association study and whole-genome sequencing. Together, our genetic, cell-based, and in vivo modeling data delineate an autosomal-dominant syndromic retinal ciliopathy in humans and suggest that multiple KIF3B pathomechanisms can impair kinesin-driven ciliary transport in the photoreceptor.
Hysi PG, Choquet H, Khawaja AP, Wojciechowski R, Tedja MS, Yin J, Simcoe MJ, Patasova K, Mahroo OA, Thai KK, Cumberland PM, Melles RB, Verhoeven VJM, Vitart V, Segre A, Stone RA, Wareham N, Hewitt AW, Mackey DA, Klaver CCW, Macgregor S, for and Myopia CRE, Khaw PT, Foster PJ, and Consortium UKEV, Guggenheim JA, Guggenheim JA, Rahi JS, Jorgenson E, Hammond CJ. Meta-analysis of 542,934 subjects of European ancestry identifies new genes and mechanisms predisposing to refractive error and myopia. Nat Genet 2020;52(4):401-407.Abstract
Refractive errors, in particular myopia, are a leading cause of morbidity and disability worldwide. Genetic investigation can improve understanding of the molecular mechanisms that underlie abnormal eye development and impaired vision. We conducted a meta-analysis of genome-wide association studies (GWAS) that involved 542,934 European participants and identified 336 novel genetic loci associated with refractive error. Collectively, all associated genetic variants explain 18.4% of heritability and improve the accuracy of myopia prediction (area under the curve (AUC) = 0.75). Our results suggest that refractive error is genetically heterogeneous, driven by genes that participate in the development of every anatomical component of the eye. In addition, our analyses suggest that genetic factors controlling circadian rhythm and pigmentation are also involved in the development of myopia and refractive error. These results may enable the prediction of refractive error and the development of personalized myopia prevention strategies in the future.
Almeida LM, Lebreton F, Gaca A, Bispo PM, Saavedra JT, Calumby RN, Grillo LM, Nascimento TG, Filsner PH, Moreno AM, Gilmore MS. Transferable Resistance Gene in Enterococcus faecalis from Swine in Brazil. Antimicrob Agents Chemother 2020;64(6)Abstract
OptrA is an ATP-binding cassette (ABC)-F protein that confers resistance to oxazolidinones and phenicols and can be either plasmid-encoded or chromosomally encoded. Here, we isolated 13 strains possessing a linezolid MIC of ≥4 mg/liter from nursery pigs in swine herds located across Brazil. Genome sequence comparison showed that these strains possess in different genetic contexts occurring in 5 different sequence type backgrounds. The gene invariably occurred in association with an regulator and a gene encoding a hypothetical protein. In some contexts, this genetic island was able to excise and form a covalently closed circle within the cell; this circle appeared to occur in high abundance and to be transmissible by coresident plasmids.
Bronstein R, Capowski EE, Mehrotra S, Jansen AD, Navarro-Gomez D, Maher M, Place E, Sangermano R, Bujakowska KM, Gamm DM, Pierce EA. A combined RNA-seq and whole genome sequencing approach for identification of non-coding pathogenic variants in single families. Hum Mol Genet 2020;29(6):967-979.Abstract
Inherited retinal degenerations (IRDs) are at the focus of current genetic therapeutic advancements. For a genetic treatment such as gene therapy to be successful, an accurate genetic diagnostic is required. Genetic diagnostics relies on the assessment of the probability that a given DNA variant is pathogenic. Non-coding variants present a unique challenge for such assessments as compared to coding variants. For one, non-coding variants are present at much higher number in the genome than coding variants. In addition, our understanding of the rules that govern the non-coding regions of the genome is less complete than our understanding of the coding regions. Methods that allow for both the identification of candidate non-coding pathogenic variants and their functional validation may help overcome these caveats allowing for a greater number of patients to benefit from advancements in genetic therapeutics. We present here an unbiased approach combining whole genome sequencing (WGS) with patient-induced pluripotent stem cell (iPSC)-derived retinal organoids (ROs) transcriptome analysis. With this approach, we identified and functionally validated a novel pathogenic non-coding variant in a small family with a previously unresolved genetic diagnosis.
Levy JM, Yeh W-H, Pendse N, Davis JR, Hennessey E, Butcher R, Koblan LW, Comander J, Liu Q, Liu DR. Cytosine and adenine base editing of the brain, liver, retina, heart and skeletal muscle of mice via adeno-associated viruses. Nat Biomed Eng 2020;4(1):97-110.Abstract
The success of base editors for the study and treatment of genetic diseases depends on the ability to deliver them in vivo to the relevant cell types. Delivery via adeno-associated viruses (AAVs) is limited by AAV packaging capacity, which precludes the use of full-length base editors. Here, we report the application of dual AAVs for the delivery of split cytosine and adenine base editors that are then reconstituted by trans-splicing inteins. Optimized dual AAVs enable in vivo base editing at therapeutically relevant efficiencies and dosages in the mouse brain (up to 59% of unsorted cortical tissue), liver (38%), retina (38%), heart (20%) and skeletal muscle (9%). We also show that base editing corrects, in mouse brain tissue, a mutation that causes Niemann-Pick disease type C (a neurodegenerative ataxia), slowing down neurodegeneration and increasing lifespan. The optimized delivery vectors should facilitate the efficient introduction of targeted point mutations into multiple tissues of therapeutic interest.
Perez-Cervantes C, Smith LA, Nadadur RD, Hughes AEO, Wang S, Corbo JC, Cepko C, Lonfat N, Moskowitz IP. Enhancer transcription identifies -regulatory elements for photoreceptor cell types. Development 2020;147(3)Abstract
Identification of cell type-specific regulatory elements (CREs) is crucial for understanding development and disease, although identification of functional regulatory elements remains challenging. We hypothesized that context-specific CREs could be identified by context-specific non-coding RNA (ncRNA) profiling, based on the observation that active CREs produce ncRNAs. We applied ncRNA profiling to identify rod and cone photoreceptor CREs from wild-type and mutant mouse retinas, defined by presence or absence, respectively, of the rod-specific transcription factor (TF) -dependent ncRNA expression strongly correlated with epigenetic profiles of rod and cone photoreceptors, identified thousands of candidate rod- and cone-specific CREs, and identified motifs for rod- and cone-specific TFs. Colocalization of NRL and the retinal TF CRX correlated with rod-specific ncRNA expression, whereas CRX alone favored cone-specific ncRNA expression, providing quantitative evidence that heterotypic TF interactions distinguish cell type-specific CRE activity. We validated the activity of novel -dependent ncRNA-defined CREs in developing cones. This work supports differential ncRNA profiling as a platform for the identification of cell type-specific CREs and the discovery of molecular mechanisms underlying TF-dependent CRE activity.
Shakhmantsir I, Dooley SJ, Kishore S, Chen D, Pierce E, Bennett J, Sehgal A. RNA Splicing Factor Mutations That Cause Retinitis Pigmentosa Result in Circadian Dysregulation. J Biol Rhythms 2020;35(1):72-83.Abstract
Circadian clocks regulate multiple physiological processes in the eye, but their requirement for retinal health remains unclear. We previously showed that Drosophila homologs of spliceosome proteins implicated in human retinitis pigmentosa (RP), the most common genetically inherited cause of blindness, have a role in the brain circadian clock. In this study, we report circadian phenotypes in murine models of RP. We found that mice carrying a homozygous H2309P mutation in () display a lengthened period of the circadian wheel-running activity rhythm. We show also that the daily cycling of circadian gene expression is dampened in the retina of H2309P mice. Surprisingly, molecular rhythms are intact in the eye cup, which includes the retinal pigment epithelium (RPE), even though the RPE is thought to be the primary tissue affected in this form of RP. Downregulation of , another RNA splicing factor implicated in RP, leads to period lengthening in a human cell culture model. The period of circadian bioluminescence in primary fibroblasts of human RP patients is not significantly altered. Together, these studies link a prominent retinal disorder to circadian deficits, which could contribute to disease pathology.
Dudek AM, Zabaleta N, Zinn E, Pillay S, Zengel J, Porter C, Franceschini JS, Estelien R, Carette JE, Zhou GL, Vandenberghe LH. GPR108 Is a Highly Conserved AAV Entry Factor. Mol Ther 2020;28(2):367-381.Abstract
Adeno-associated virus (AAV) is a highly promising gene transfer vector, yet major cellular requirements for AAV entry are poorly understood. Using a genome-wide CRISPR screen for entry of evolutionarily divergent serotype AAVrh32.33, we identified GPR108, a member of the G protein-coupled receptor superfamily, as an AAV entry factor. Of greater than 20 divergent AAVs across all AAV clades tested in human cell lines, only AAV5 transduction was unaffected in the GPR108 knockout (KO). GPR108 dependency was further shown in murine and primary cells in vitro. These findings are further validated in vivo, as the Gpr108 KO mouse demonstrates 10- to 100-fold reduced expression for AAV8 and rh32.33 but not AAV5. Mechanistically, both GPR108 N- and C-terminal domains are required for transduction, and on the capsid, a VP1 unique domain that is not conserved on AAV5 can be transferred to confer GPR108 independence onto AAV2 chimeras. In vitro binding and fractionation studies indicate reduced nuclear import and cytosolic accumulation in the absence of GPR108. We thus have identified the second of two AAV entry factors that is conserved between mice and humans relevant both in vitro and in vivo, further providing a mechanistic understanding to the tropism of AAV gene therapy vectors.
Choquet H, Wiggs JL, Khawaja AP. Clinical implications of recent advances in primary open-angle glaucoma genetics. Eye (Lond) 2020;34(1):29-39.Abstract
Over the last decade, genetic studies, including genome-wide association studies (GWAS), have accelerated the discovery of genes and genomic regions contributing to primary open-angle glaucoma (POAG), a leading cause of irreversible vision loss. Here, we review the findings of genetic studies of POAG published in English prior to September 2019. In total, 74 genomic regions have been associated at a genome-wide level of significance with POAG susceptibility. Recent POAG GWAS provide not only insight into global and ethnic-specific genetic risk factors for POAG susceptibility across populations of diverse ancestry, but also important functional insights underlying biological mechanisms of glaucoma pathogenesis. In this review, we also summarize the genetic overlap between POAG, glaucoma endophenotypes, such as intraocular pressure and vertical cup-disc ratio (VCDR), and other eye disorders. We also discuss approaches recently developed to increase power for POAG locus discovery and to predict POAG risk. Finally, we discuss the recent development of POAG gene-based therapies and future strategies to treat glaucoma effectively. Understanding the genetic architecture of POAG is essential for an earlier diagnosis of this common eye disorder, predictive testing of at-risk patients, and design of gene-based targeted medical therapies none of which are currently available.