Integration of genomics and transcriptomics predicts diabetic retinopathy susceptibility genes
Abstract
We determined differential gene expression in response to high glucose in lymphoblastoid cell lines derived from matched individuals with type 1 diabetes with and without retinopathy. Those genes exhibiting the largest difference in glucose response were assessed for association to diabetic retinopathy in a genome-wide association study meta-analysis. Expression Quantitative Trait Loci (eQTLs) of the glucose response genes were tested for association with diabetic retinopathy. We detected an enrichment of the eQTLs from the glucose response genes among small association p-values and identified FLCN as a susceptibility gene for diabetic retinopathy. Expression of FLCN in response to glucose was greater in individuals with diabetic retinopathy. Independent cohorts of individuals with diabetes revealed an association of FLCN eQTLs to diabetic retinopathy. Mendelian randomization confirmed a direct positive effect of increased FLCN expression on retinopathy. Integrating genetic association with gene expression implicated FLCN as a disease gene for diabetic retinopathy.
Data availability
Source files and code for all the figures and tables have been provided, except for drawings, flowcharts and histopathology findings. We have also included links and references where appropriate.Figure 3 source data 5 and 6 are available on Dryad at https://doi.org/10.5061/dryad.zkh18938jAdditional data files can be found here: microarray expression data at Gene Expression Omnibus (GEO) under accession code GSE146615 and diabetic retinopathy GWAS data at UKBB archive (https://oxfile.ox.ac.uk/oxfile/work/extBox?id=825146B4380F72048D).
-
Figure 3 additional source data filesDryad Digital Repository, doi:10.5061/dryad.zkh18938j.
Article and author information
Author details
Funding
National Eye Institute (R01EY023644)
- Michael A Grassi
National Eye Institute (ZIAEY000546)
- Anand Swaroop
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Reviewing Editor
- David E James, The University of Sydney, Australia
Publication history
- Received: June 13, 2020
- Accepted: November 6, 2020
- Accepted Manuscript published: November 9, 2020 (version 1)
- Version of Record published: December 10, 2020 (version 2)
Copyright
© 2020, Skol et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,004
- Page views
-
- 203
- Downloads
-
- 2
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Genetics and Genomics
Splicing is highly regulated and is modulated by numerous factors. Quantitative predictions for how a mutation will affect precursor mRNA (pre-mRNA) structure and downstream function are particularly challenging. Here, we use a novel chemical probing strategy to visualize endogenous precursor and mature MAPT mRNA structures in cells. We used these data to estimate Boltzmann suboptimal structural ensembles, which were then analyzed to predict consequences of mutations on pre-mRNA structure. Further analysis of recent cryo-EM structures of the spliceosome at different stages of the splicing cycle revealed that the footprint of the Bact complex with pre-mRNA best predicted alternative splicing outcomes for exon 10 inclusion of the alternatively spliced MAPT gene, achieving 74% accuracy. We further developed a β-regression weighting framework that incorporates splice site strength, RNA structure, and exonic/intronic splicing regulatory elements capable of predicting, with 90% accuracy, the effects of 47 known and 6 newly discovered mutations on inclusion of exon 10 of MAPT. This combined experimental and computational framework represents a path forward for accurate prediction of splicing-related disease-causing variants.
-
- Cell Biology
- Genetics and Genomics
In the first meiotic cell division, proper segregation of chromosomes in most organisms depends on chiasmata, exchanges of continuity between homologous chromosomes that originate from the repair of programmed double-strand breaks (DSBs) catalyzed by the Spo11 endonuclease. Since DSBs can lead to irreparable damage in germ cells, while chromosomes lacking DSBs also lack chiasmata, the number of DSBs must be carefully regulated to be neither too high nor too low. Here, we show that in Caenorhabditis elegans, meiotic DSB levels are controlled by the phosphoregulation of DSB-1, a homolog of the yeast Spo11 cofactor Rec114, by the opposing activities of PP4PPH-4.1 phosphatase and ATRATL-1 kinase. Increased DSB-1 phosphorylation in pph-4.1 mutants correlates with reduction in DSB formation, while prevention of DSB-1 phosphorylation drastically increases the number of meiotic DSBs both in pph-4.1 mutants as well as in the wild type background. C. elegans and its close relatives also possess a diverged paralog of DSB-1, called DSB-2, and loss of dsb-2 is known to reduce DSB formation in oocytes with increasing age. We show that the proportion of the phosphorylated, and thus inactivated, form of DSB-1 increases with age and upon loss of DSB-2, while non-phosphorylatable DSB-1 rescues the age-dependent decrease in DSBs in dsb-2 mutants. These results suggest that DSB-2 evolved in part to compensate for the inactivation of DSB-1 through phosphorylation, to maintain levels of DSBs in older animals. Our work shows that PP4PPH-4.1, ATRATL-1, and DSB-2 act in concert with DSB-1 to promote optimal DSB levels throughout the reproductive lifespan.