Abstract
Age-related macular degeneration (AMD) is the most common cause of irreversible vision loss in the developed world which affects the quality of life for millions of elderly individuals worldwide. Genome-wide association studies (GWAS) have identified genetic variants at 34 loci contributing to AMD. To better understand the disease pathogenesis and identify causal genes for AMD, we applied random walk (RW) and support vector machine (SVM) to identify AMD-related genes based on gene interaction relationship and significance of genes. Our model achieved 0.927 of area under the curve (AUC), and 65 novel genes have been identified as AMD-related genes. To verify our results, a statistics method called summary data-based Mendelian randomization (SMR) has been implemented to integrate GWAS data and transcriptome data to verify AMD susceptibility-related genes. We found 45 genes are related to AMD by SMR. Among these genes, 37 genes overlap with those found by SVM-RW. Finally, we revealed the biological process of genetic mutations leading to changes in gene expression leading to AMD. Our results reveal the genetic pathogenic factors and related mechanisms of AMD.
Introduction
Age-related macular degeneration (AMD) is the most common cause of irreversible blindness with limited therapeutic options in the elderly in many countries (). AMD causes decreased photoreceptor function in the macular area of the retina (). Researchers have found many factors which are related to the development and severity of AMD.
Genetic factors are significantly related to AMD. In 2005, Klein et al. found that CFH gene was related to AMD, which was the first discovered AMD-related gene (). This gene is significantly expressed in retinal pigment epithelial cells. Y402H mutation of CFH impairs the complement pathway regulation function of CFH gene (). Subsequently, the ARMS2 gene cluster was also found to be related to AMD. Multiple studies have shown that there is a strong correlation between multiple genetic variants in this gene cluster and AMD (). Recently, it has been discovered that the apolipoprotein E (APOE) gene has a strong correlation with AMD (). The APOE gene plays a role in transporting lipids and cholesterol in the central nervous system, and multiple studies have shown that this gene is associated with neurodegenerative diseases such as Alzheimer’s disease and stroke (; , ). The gene is expressed on photoreceptor cells, retinal ganglion cells, retinal pigment epithelial cells, Bruch’s membrane, and the choroid. Most studies have proved APOE can prevent AMD (). The genetic risk of advanced AMD is increased (). Researchers have found that the heritability estimate for twin studies is 0.45 for early AMD () but 0.71 for late AMD ().
Computational methods have been widely used to discover functions of biological molecules (,, ). AMD-related genome-wide association studies (GWAS) analyses have identified a strong association of 52 independent single-nucleotide polymorphisms (SNPs) at 34 genetic loci accounting for over 50% of the genetic heritability (). Machine learning methods can help researchers find disease-related information on a large scale. However, these methods cannot explain the genetic mechanism of the results. GWAS studies are a valuable resource for understanding disease pathologies, but they may not precisely point out the causal genes responsible for the disease of interest. Besides, there have been studies that reported that causal genes are distinct from the nearest genes discovered by GWAS (; ). However, The gene expression is related to the genetic variant so the gene expression levels are different in different genotypes (). Expression quantitative trait locus (eQTL) mapping offers a powerful approach to elucidate the genetic component underlying altered gene expression. Gene expression is vital for complex diseases () and is also differentially regulated across tissues, such as the brain, heart, and pancreas. have found potential causal genes in six AMD GWAS loci from human retinal samples. However, that analysis only considered retinal samples and was not comprehensive since it is difficult to obtain multiple living tissues and most eQTL studies so far have been performed with RNA isolated from immortalized lymphoblasts or lymphocytes. In this study, we fused random walk (RW) with support vector machine (SVM) to identify AMD-related genes. Since many GWAS and eQTL studies have been made public, to verify our results, AMD GWAS data and blood eQTL studies are integrated to further find expression of the genes related to AMD. In this method, we referred to the concept of Mendelian randomization (MR) analysis (; ), where a genetic variant (such as a SNP) is considered as an instrumental variable (such as gene expression) to validate for the causative effect of an exposure on an outcome (such as a phenotype). Based on this assumption, we can obtain AMD-related genes based on MR. We collected eQTL data from the GTEx database and collected GWAS datasets including 12,711 advanced AMD cases and 14,590 controls of European descent from a study by ; 707 Caucasian AMD patients and 2,014 controls from a study by ; and 14,034 cases, 91,214 controls, and 11 sources of data including the International AMD Genomics Consortium, IAMDGC, and United Kingdom Biobank (UKBB) from a study by . Based on these GWAS studies and eQTL dataset, we can not only identify genes related to AMD but also speculate on their biological processes.
Materials and Methods
Encoding Gene Interaction Network by Random Walk
The RW algorithm is a method that is simple to operate but not easy to fall into a local minimum. We constructed a gene interaction network by known AMD-related genes and a string database. Then, we implemented RW on the gene interaction network.
f(x) is a multivariate function with n variables; x = (x1,x2,…,xn) is an n dimension vector.
Step 1: Given the initial iteration point x, λ is the first walking step length, and ϵ is the control accuracy (ϵ is a very small positive number, used to control the end of the algorithm).
Step 2: Given the number of iterations control N, k is the current iteration number; set k = 1.
Step 3: When k < N, randomly generate an n-dimensional vector between (−1, 1). u = (u1,u2,…,un), (−1 < ui < 1, i = 1, 2,…, n), and standardize it to get u′.
Let x1 = x + λu′ to complete the first step of walking.
Step 4: Calculate the value of the function, if f(x1) < f(x), which is a better point than the initial value, then reset k to 1, change x1 to x, and go back to step 2; otherwise, k = k + 1. Go back to step 3.
Step 5: If no better value can be found for N consecutive times, it is considered that the optimal solution is within the N-dimensional sphere with the current optimal solution as the center and the current step as the radius (if it is three-dimensional, it just happens to be in the space sphere). At this point, if λ < ϵ, the algorithm ends; otherwise, let λ = λ2, go back to step 1, and start a new round of walking.
Finally, we can get the gene feature after encoding the gene network.
Classification by Support Vector Machine
We obtained the gene feature in the last section. Then, we can input the gene feature and label into SVM to get the relationship between the gene and AMD. The workflow of SVM is shown in Figure 1.
FIGURE 1
First, we used Z-score normalization to process the gene feature. Then, we constructed a Lagrangian function to obtain the values and dualized the original problem. Sequential minimal optimization (SMO) algorithm was used to solve the dualization problem. Finally, we can obtain the classification model and output the prediction results.
Results
AMD-Related Genes Identification by SVM-RW
We obtained 34 known AMD-related genes from GWAS data. We constructed a gene network which has 239 nodes (genes). We did 10-cross validation by SVM-RW and tested the performance of SVM-RW. The area under the curve (AUC) of SVM-RW is shown in Figure 2.
FIGURE 2
SVM-RW achieved AUC of 0.927 in identifying AMD-related genes. We compared the results of SVM-RW with several other methods. The results are shown in Table 1.
TABLE 1
| Algorithm | AUC | AUPR |
| SVM-RW | 0.927 | 0.781 |
| Random forest-RW | 0.852 | 0.645 |
| Naive Bayes-RW | 0.711 | 0.586 |
| Backpropagation-artificial neural network-RW | 0.823 | 0.692 |
| Logistic regression-RW | 0.691 | 0.531 |
Comparison results.
AUC, area under the curve; AUPR, area under the precision-recall curve; RW, random walk.
Bold values highlight the result of SVM-RW.
After verifying the effectiveness of SVM-RW, we randomly selected 34 genes as negative samples and built a final SVM model. SVM-RW predicted 65 novel genes as AMD-related genes.
Verify SVM-RW Results by Summary Data Level-Mendelian Randomization Analysis
If we use g to denote a genetic variant (such as a SNP), x as the expression level of a gene, and y as the trait, then the two-step least-squares (2SLS) estimate of the effect of x on y from an MR analysis can be denoted as:
where and indicate the least-squares estimates of y and x on z, respectively, and Exy indicates the effect size of x on y free of confounding from non-genetic factors. Then the sampling variance of the 2SLS estimate of Exy can be denoted as:
where n denotes the sample size, indicates the proportion of variance in the explanation of y by x, and is the proportion of variance in the explanation of y by z. Therefore, we use the statistic TMR to test the significance of Exy; TMR can be denoted as:
where TMR=. Based on the suggestion that the power of detecting Exy can be significantly increased using a two-sample MR analysis (; ), if GWAS and eQTL datasets share the same population, we can use unbiased estimates to replace Ezx. We therefore have
where is the estimate of a SNP effect from a GWAS for a trait, and is the estimate of a SNP effect on the expression level of a gene from an eQTL study. The sampling variance of can be approximately computed by the Delta method () as:
where is 0. Based on experience, we can replace them by their estimates in practice, indicated as an approximate χ2 test statistic of:
where zxy is the z statistic of the GWAS and zzx is the z statistic of the eQTL study.
In an MR analysis, Exy is interpreted as the effect of a phenotype on the gene expression without considering non-genetic confounders. We first collected GWAS summary data and blood eQTL data from available online studies. We first collected a GWAS summary dataset composed of 12,711 advanced AMD cases and 14,590 controls of European descent from the study by ; 707 Caucasian AMD patients and 2,014 controls from the study by ; and 14,034 cases, 91,214 controls, and 11 sourced from datasets including the International AMD Genomics Consortium, IAMDGC, and UKBB from the study by . The distribution of the above datasets is shown in Figures 3A,B.
FIGURE 3
Then summary data-based Mendelian randomization (SMR) analysis is implemented on the blood eQTL data and GWAS data; in this paper, we identified 48 SNPs regulating 45 genes (including 41 coding genes and four non-coding genes) resulting in AMD susceptibility. The workflow is shown in Figure 4.
FIGURE 4
For the first GWAS datasets consisting of 12,711 AMD cases and 14,590 controls from European cohorts, in total we found 3,872 SNPs coexist in both GWAS data and eQTL data; 43 of 3,872 SNPs are significant and regulate 44 genes in gene expression level. In the second GWAS dataset, we found 714 SNPs coexist in both GWAS dataset and eQTL dataset, with none significant. In the third GWAS dataset, we found 1,149 SNPs coexist both in GWAS dataset and eQTL dataset, with one significant regulating one gene in gene expression level. The distribution of the p-value of SNPs regulating genes tested by SMR is shown in Figures 5A–C. A Supplementary Table 1 indicates the p-values of significant SNPs regulating genes tested by SMR; the last line resulted from GWAS dataset 3, and the rest resulted from GWAS dataset 1.
FIGURE 5
Case Study
Age-related macular degeneration has been described as a partly genetic disease (; ). Recently, a unifying hypothesis is that immune response gene polymorphisms modulate susceptibility to AMD. Human leukocyte antigen (HLA) polymorphisms, encoded within the major histocompatibility complex (MHC), are the most polymorphic within the human genome. In AMD, researchers detected intense HLA-DR immunoreactivity in not only soft but also hard drusen (). In the study of , considering the effect of smoking, age, and body mass index (BMI), HLA alleles B∗4001, DRB1∗1301, and Cw∗0701 were found to be related to AMD, which is consistent with our results displayed in Table 1.
In a study by , they researched P2RX7 and P2RX4 genes in 744 AMD patients and 557 Caucasian controls and reached a conclusion that a rare functional haplotype of the P2RX4 leads to loss of innate phagocytosis and confers increased risk of AMD. P2RX7 and P2RX4 damage the normal scavenger function of macrophages and microglia through interaction, making individuals susceptible to AMD.
Gene Interaction Network Based on AMD
Figure 6 shows the gene interaction network produced from the results of SMR on AMD. Based on the interaction network, the HLA class intensively interacted and is significantly associated with AMD.
FIGURE 6
The cluster consisting of DDX39B (aka BAT1), PRRC2A (aka BAT2), and SKIV2L are genes found in the class III region of the MHC (MHC Class III). These genes encode RNA-binding proteins with clear roles in post-transcriptional gene regulation and RNA surveillance. They are likely to have important functions in immunity and are associated with autoimmune diseases (). Early work by immunologists have shown that DDX39B promoted gene expression of anti-inflammatory pathways (). Therefore, understanding the genes interactions may help speculate on the proposed AMD mechanisms and immunotherapy.
Conclusion
We applied the SMR method on AMD to test the gene–AMD associations based on GWAS summary data and blood eQTL data. From a total of 27,452 AMD cases and 107,818 controls, we obtained 44 SNPs regulating 45 genes significantly associated with AMD. Among the results, HLA class genes have been proved to be associated with immunologically mediated diseases because of the critical role of HLA in mediating the immune response, and genes from MHC Class III are also associated with autoimmune diseases. These genes may play important roles in causing AMD susceptibility and need to be further verified with experiments. Since AMD has been considered as a genetic disease, from this perspective, it is helpful in understanding the disease from gene-expression level to speculate about the AMD mechanisms and pathology and propose future treatment options for AMD.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author contributions
YD, NK, and JZ participated in its design, analyzed the data, and wrote the manuscript. All authors read and approved the published version of the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.726599/full#supplementary-material
References
1
AllcockR. J.WilliamsJ. H.PriceP. (2001). The central MHC gene, BAT1, may encode a protein that down-regulates cytokine production.Genes Cells6487–494. 10.1046/j.1365-2443.2001.00435.x
2
ClaussnitzerM.DankelS. N.KimK.-H.QuonG.MeulemanW.HaugenC.et al (2015). FTO obesity variant circuitry and adipocyte browning in humans.N. Engl. J. Med.373895–907. 10.1056/nejmoa1502214
3
Davey SmithG.EbrahimS. (2003). ‘Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease?Int. J. Epidemiol.321–22. 10.1093/ije/dyg070
4
FeherJ.KovacsI.ArticoM.CavallottiC.PapaleA.GabrieliC. B. (2006). Mitochondrial alterations of retinal pigment epithelium in age-related macular degeneration.Neurobiol. Aging27983–993. 10.1016/j.neurobiolaging.2005.05.012
5
Fernández-VegaB.GarcíaM.OlivaresL.ÁlvarezL.González-FernándezA.ArtimeE.et al (2020). The association study of lipid metabolism gene polymorphisms with AMD identifies a protective role for APOE-E2 allele in the wet form in a Northern Spanish population.Acta Ophthalmol.98e282–e291. 10.1111/aos.14280
6
FritscheL. G.FarissR. N.StambolianD.AbecasisG. R.CurcioC. A.SwaroopA. (2014). Age-related macular degeneration: genetics and biology coming together.Annu. Rev. Genomics Hum. Genet.15151–171. 10.1146/annurev-genom-090413-025610
7
FritscheL. G.IglW.BaileyJ. N. C.GrassmannF.SenguptaS.Bragg-GreshamJ. L.et al (2016). A large genome-wide association study of age-related macular degeneration highlights contributions of rare and common variants.Nat. Genet.48134–143. 10.1038/ng.3448
8
GoverdhanS. V.HowellM. W.MullinsR. F.OsmondC.HodgkinsP. R.SelfJ.et al (2005). Association of HLA class I and class II polymorphisms with age-related macular degeneration.Invest. Ophthalmol. Vis. Sci.461726–1734. 10.1167/iovs.04-0928
9
GuB. J.BairdP. N.VesseyK. A.SkarrattK. K.FletcherE. L.FullerS. J.et al (2013). A rare functional haplotype of the P2RX4 and P2RX7 genes leads to loss of innate phagocytosis and confers increased risk of age-related macular degeneration.FASEB J.271479–1487. 10.1096/fj.12-215368
10
HainesJ. L.HauserM. A.SchmidtS.ScottW. K.OlsonL. M.GallinsP.et al (2005). Complement factor H variant increases the risk of age-related macular degeneration.Science308419–421. 10.1126/science.1110359
11
HammondC. J.WebsterA. R.SniederH.BirdA. C.GilbertC. E.SpectorT. D. (2002). Genetic influence on early age-related maculopathy: a twin study.Ophthalmology109730–736. 10.1016/s0161-6420(01)01049-1
12
HanX.OngJ.-S.AnJ.HewittA. W.GharahkhaniP.MacGregorS. (2020). Using mendelian randomization to evaluate the causal relationship between serum C-reactive protein levels and age-related macular degeneration.Eur. J. Epidemiol.35139–146. 10.1007/s10654-019-00598-z
13
HeibaI. M.ElstonR. C.KleinB. E.KleinR. (1994). Sibling correlations and segregation analysis of age-related maculopathy: the beaver dam eye study.Genet. Epidemiol.1151–67. 10.1002/gepi.1370110106
14
InoueA.SolonG. (2010). Two-sample instrumental variables estimators.Rev. Econ. Stat.92557–561. 10.1162/rest_a_00011
15
JohnsonL. V.LeitnerW. P.StaplesM. K.AndersonD. H. (2001). Complement activation and inflammatory processes in drusen formation and age related macular degeneration.Exp. Eye Res.73887–896. 10.1006/exer.2001.1094
16
KatanM. B. (2004). Apolipoprotein E isoforms, serum cholesterol, and cancer.Int. J. Epidemiol.33:9. 10.1093/ije/dyh312
17
LandowskiM.KellyU.KlingebornM.GroelleM.DingJ.-D.GrigsbyD.et al (2019). Human complement factor H Y402H polymorphism causes an age-related macular degeneration phenotype and lipoprotein dysregulation in mice.Proc. Natl. Acad. Sci. U.S.A.1163703–3711. 10.1073/pnas.1814014116
18
LimL. S.MitchellP.SeddonJ. M.HolzF. G.WongT. Y. (2012). Age-related macular degeneration.Lancet3791728–1738. 10.1016/S0140-6736(12)60282-7
19
LynchM.WalshB. (1998). Genetics and Analysis of Quantitative Traits.Sunderland, MA: Sinauer Associates Inc.
20
MullinsR. F.RussellS. R.AndersonD. H.HagemanG. S. (2000). Drusen associated with aging and age-related macular degeneration contain proteins common to extracellular deposits associated with atherosclerosis, elastosis, amyloidosis, and dense deposit disease.FASEB J.14835–846. 10.1096/fasebj.14.7.835
21
PangC.BaumL.ChanW.LauT.PoonP.LamD. (2000). The apolipoprotein E ε4 allele is unlikely to be a major risk factor of age-related macular degeneration in Chinese.Ophthalmologica214289–291. 10.1159/000027506
22
PierceB. L.BurgessS. (2013). Efficient design for mendelian randomization studies: subsample and 2-sample instrumental variable estimators.Am. J. Epidemiol.1781177–1184. 10.1093/aje/kwt084
23
RatnapriyaR.SosinaO. A.StarostikM. R.KwicklisM.KapphahnR. J.FritscheL. G.et al (2019). Retinal transcriptome and eQTL analyses identify genes associated with age-related macular degeneration.Nat. Genet.51606–610. 10.1038/s41588-019-0351-9
24
SchottG.Garcia-BlancoM. A. (2020). MHC class III RNA binding proteins and immunity.RNA Biol.18640–646. 10.1080/15476286.2020.1860388
25
SeddonJ. M.GeorgeS.RosnerB.RifaiN. (2005). Progression of age-related macular degeneration: prospective assessment of C-reactive protein, interleukin 6, and other cardiovascular biomarkers.Arch. Ophthalmol.123774–782. 10.1001/archopht.123.6.774
26
SmemoS.TenaJ. J.KimK.-H.GamazonE. R.SakabeN. J.Gómez-MarínC.et al (2014). Obesity-associated variants within FTO form long-range functional connections with IRX3.Nature507371–375. 10.1038/nature13138
27
StoneE. M.BraunT. A.RussellS. R.KuehnM. H.LoteryA. J.MooreP. A.et al (2004). Missense variations in the fibulin 5 gene and age-related macular degeneration.N. Engl. J. Med.351346–353. 10.1056/nejmoa040833
28
WinklerT. W.GrassmannF.BrandlC.KielC.GüntherF.StrunzT.et al (2020). Genome-wide association meta-analysis for early age-related macular degeneration highlights novel loci and insights for advanced disease.BMC Med. Genomics13:120. 10.1186/s12920-020-00760-7
29
YanQ.DingY.LiuY.SunT.FritscheL. G.ClemonsT.et al (2018). Genome-wide analysis of disease progression in age-related macular degeneration.Hum. Mol. Genet.27929–940. 10.1093/hmg/ddy002
30
ZhaoT.HuY.ChengL. (2020a). Deep-DRM: a computational method for identifying disease-related metabolites based on graph deep learning approaches.Brief. Bioinform.9:bbaa212. 10.1093/bib/bbaa212
31
ZhaoT.HuY.PengJ.ChengL. (2020b). DeepLGP: a novel deep learning method for prioritizing lncRNA target genes.Bioinformatics364466–4472. 10.1093/bioinformatics/btaa428
32
ZhaoT.HuY.ValsdottirL. R.ZangT.PengJ. (2021a). Identifying drug–target interactions based on graph convolutional network and deep neural network.Brief. Bioinform.222141–2150. 10.1093/bib/bbaa044
33
ZhaoT.HuY.ZangT.ChengL. (2019). Identifying alzheimer’s disease-related proteins by LRRGD.BMC Bioinformatics20:570. 10.1186/s12859-019-3124-7
34
ZhaoT.HuY.ZangT.ChengL. (2020c). MRTFB regulates the expression of NOMO1 in colon.Proc. Natl. Acad. Sci. U.S.A.1177568–7569. 10.1073/pnas.2000499117
35
ZhaoT.HuY.ZangT.WangY. (2020d). Identifying protein biomarkers in blood for alzheimer’s disease.Front. Cell Dev. Biol.8:472. 10.3389/fcell.2020.00472
36
ZhaoT.LyuS.LuG.JuanL.ZengX.WeiZ.et al (2021b). SC2disease: a manually curated database of single-cell transcriptome for human diseases.Nucleic Acids Res.49D1413–D1419. 10.1093/nar/gkaa838
Summary
Keywords
AMD, GWAS, eQTL, SNPs, disease susceptibility
Citation
Du Y, Kong N and Zhang J (2021) Genetic Mechanism Revealed of Age-Related Macular Degeneration Based on Fusion of Statistics and Machine Learning Method. Front. Genet. 12:726599. doi: 10.3389/fgene.2021.726599
Received
17 June 2021
Accepted
13 July 2021
Published
05 August 2021
Volume
12 - 2021
Edited by
Lei Deng, Central South University, China
Reviewed by
Hong Ju, Heilongjiang Vocational College of Biology Science and Technology, China; Hui Ding, University of Electronic Science and Technology of China, China
Updates
Copyright
© 2021 Du, Kong and Zhang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jibin Zhang, kn13@tom.com
This article was submitted to Statistical Genetics and Methodology, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.