Molecular Vision 2017; 23:318-333
Received 05 December 2016 | Accepted 12 June 2017 | Published 14 June 2017
1Department of Ophthalmology; NIBR Informatics, Novartis Institutes for Biomedical Research, Cambridge, MA; 2Scientific Data Analysis, NIBR Informatics, Novartis Institutes for Biomedical Research, Cambridge, MA
Correspondence to: Sha-Mei Liao, Department of Ophthalmology, Novartis Institutes for Biomedical Research, 22 Windsor Street, Cambridge, MA 02139, USA; Phone: (617) 871 4004; FAX: (617) 871 5784; email: email@example.com
Dr. Zheng now at GNS Healthcare, Cambridge, MA
Purpose: A region within chromosome 10q26 has a set of single nucleotide polymorphisms (SNPs) that define a haplotype that confers high risk for age-related macular degeneration (AMD). We used a bioinformatics approach to search for genes in this region that may be responsible for risk for AMD by assessing levels of gene expression in individuals carrying different haplotypes and by searching for open chromatin regions in the retinal pigment epithelium (RPE) that might include one or more of the SNPs.
Methods: We surveyed the PubMed and the 1000 Genomes databases to find all common (minor allele frequency > 0.01) SNPs in 10q26 strongly associated with AMD. We used the HaploReg and LDlink databases to find sets of SNPs with alleles in linkage disequilibrium and used the Genotype-Tissue Expression (GTEx) database to search for correlations between genotypes at individual SNPs and the relative level of expression of the genes. We also accessed Encyclopedia of DNA Elements (ENCODE) to find segments of open chromatin in the region with the AMD-associated SNPs. Predicted transcription factor binding motifs were identified using HOMER, PROMO, and RegulomeDB software programs.
Results: There are 34 polymorphisms within a 30-kb region that are in strong linkage disequilibrium (r2>0.8) with the reference SNP rs10490924 previously associated with risk for AMD. The expression of three genes in this region, PLEKHA1, ARMS2, and HTRA1 varies between people who have the low-AMD-risk haplotype compared with those with the high-AMD-risk haplotype. For PLEKHA1, 44 tissues have an expression pattern with the high-AMD-risk haplotype associated with low expression (rs10490924 effect size -0.43, p = 3.8 x 10-5 in ovary). With regard to ARMS2, the variation is most pronounced in testes: homozygotes with the high-AMD-risk haplotype express ARMS2 at lower levels than homozygotes with the low-AMD-risk haplotype; expression in heterozygotes falls in between (rs10490924 effect size -0.79, p = 7.5 x 10-24). For HTRA1, the expression pattern is the opposite; the high-AMD-risk haplotype has higher levels of expression in 27 tissues (rs10490924 effect size 0.40, p = 1.5 × 10−7 in testes). None of the other 22 genes within one megabase of rs10490924, or any gene in the entire genome, have mRNA expression levels that correlate with the high-AMD-risk haplotype. More than 100 other SNPs in the 10q26 region affect the expression of PLEKHA1 and ARMS2 but not that of HTRA1; none of these SNPs affects the risk for AMD according to published genome-wide association studies (GWASs). Two of the AMD-risk SNPs (rs36212732 and rs36212733) affect transcription factor binding sites in proximity to a DNase I hypersensitive region (i.e., a region of open chromatin) in RPE cells.
Conclusions: SNPs in chromosome 10q26 that influence the expression of only PLEKHA1 or ARMS2 are not associated with risk for AMD, while most SNPs that influence the expression of HTRA1 are associated with risk for AMD. Two of the AMD-risk SNPs affect transcription factor binding sites that may control expression of one of the linked genes in the RPE. These findings suggest that the variation in the risk for AMD associated with chromosome 10q26 is likely due to variation in HTRA1 expression. Modulating HTRA1 activity might be a potential therapy for AMD.
Age-related macular degeneration (AMD) is the leading cause of severe vision loss in older individuals. Single nucleotide polymorphisms (SNPs) in at least 34 loci influence the risk for AMD , and the mechanisms by which these variants influence AMD are still being elucidated. Although some AMD-risk SNPs change the coding region of genes, most reside outside coding regions, suggesting potential effects in regulating the expression of linked genes [2,3]. A set of closely linked SNPs on chromosome 10q26 is of special interest since this chromosome has more influence on the risk for AMD than any other AMD region [1,4-7]; however, which gene in the region confers the risk remains unclear [5,8-10]. Three genes are within 100 kilobase pairs (kb) of the SNPs associated with AMD. From the centromeric to telomeric ends of this region, the genes are pleckstrin homology domain-containing family A member 1 (PLEKHA1; gene ID: 59338, OMIM: 607772), age-related maculopathy-2 (ARMS2; gene ID: 387715, OMIM: 611313), and high temperature requirement protein A1 (HTRA1; gene ID: 5654, OMIM: 602194). It is even conceivable that a more distant gene actually confers the risk for AMD.
We searched for associations between AMD-risk SNPs in the 10q26 region and the expression of closely linked genes. We also looked for open chromatin in the region in RPE cells and transcription factor binding sites in the open chromatin since such regions often correspond with regions that regulate transcription.
Using the online tools HaploReg and LDlink, SNPs associated with AMD and other variants were placed in a haplotype block using Query SNP rs10490924 and the linkage disequilibrium inclusion criterion r2 ≥0.8 [11,12] in the 1000 Genome European population.
The GTEx database allows one to search for relationships between human genotypes and gene expression in specific tissues [13,14]. At the time of this analysis (January–October 2016), the GTEx project V6p (GtexPortal) contained genotypes from 449 human adult donors and whole genome transcription data from up to 53 tissues from an overlapping set of 544 donors. Transcript levels of PLEKHA1, ARMS2, and HTRA1 were included. The steady-state mRNA level is measured as reads per kilobase of transcript per million mapped reads (RPKM). RPKM are normalized according to the number of sequencing reads and the read lengths. RPKM below 1 indicate levels of mRNA expression that are difficult to distinguish from background noise.
We searched for expression quantitative trait loci (eQTLs) in the chromosome 10q26 region by looking for correlations between SNP alleles and the expression of genes using the “Test your own eQTLs” option in the GTEx portal Version 6. The mRNA levels are expressed as rank normalized gene expression . The effect size of an eQTL is the change in the value of the standardized gene expression level with each extra copy of the alternative (ALT) allele relative to the reference allele, conditional to all other adjustments (gender, probabilistic estimation of expression residuals factors, genotype principal components, and genotyping platforms). Because the gene expression levels of tissues have been transformed to a standard normal distribution (mean of 0 and standard deviation of 1), the effect size is also equivalent to the change in the population standard deviation; that is, an effect size of 0.2 means a change in 0.2 standard deviations from the baseline level with both reference alleles. The effect size provides the variation in the strength of expression with positive numbers indicating higher mRNA levels in samples from people with the minor allele compared to those with the major allele, and negative numbers indicating lower mRNA levels in samples with the minor allele. For the eQTL analysis, selected tissues were those for which high-quality mRNA results were available in at least 70 genotyped donors.
We used the Encyclopedia of DNA Elements (ENCODE) to search for DNase I sensitive sites; this database has results from more than 125 cell types, including primary cultures of RPE. We used three online programs that predict transcription factor binding sites: HOMER, PROMO, and RegulomeDB. The consensus binding sequences are from JASPAR.
The GTEx consortium database has precalculated nominal eQTL p values for every human gene and all the informative SNPs analyzed. The p values for each SNP-gene pair were calculated using a two-tailed t test as described at that site (GtexPortal). We considered associations statistically significant when the p value was less than 0.05. In some cases, as stated in the text, we made adjustments for the multiple analyses.
There are ten SNPs in a 17-kb region within 10q26 for which published data indicate strong associations of alleles with risk for AMD (Table 1 and Figure 1) [1,16-39]. One of the SNPs, rs10490929, changes the ARMS2 coding sequence (G versus T, Ala69Ser). The other nine SNPs are in strong linkage disequilibrium with rs10490929, and each has a correlation coefficient r2 ≥0.89 with rs10490924. Of these nine SNPs, rs3750846, rs3750847, and rs3750848 are in the only ARMS2 intron, rs3793917 and rs11200638 are located in the HTRA1 promoter region, rs1049331 and rs2293870 are synonymous variations in the HTRA1 coding region, and rs2284665 and rs932275 are in the first HTRA1 intron. We considered the set of alleles of these ten SNPs, all from published studies showing association with elevated risk for AMD, as a presumed haplotype that we refer to as the “high-AMD-risk” haplotype (Table 1). Although data from all ten SNPs are in the GTEx database, we used rs10490924 as the reference SNP because it is the most commonly evaluated SNP in studies of AMD [16-28]. The HaploReg and LDlink websites indicate that the haplotype defined by the ten SNPs includes an additional 24 polymorphisms (21 SNPs and three insertion-deletion polymorphisms) extending across a region of about 30 kb. Many of these 24 additional polymorphisms were not specifically evaluated in most previous genetic evaluations of AMD, but they likely have some association with risk for AMD because they have strong allelic correlations with rs10490924, with some having perfect correlations (i.e., r2 = 1; Table 1). Of these 24 additional polymorphisms, 15 are included in the GTEx database. The present analysis of eQTLs (SNP-gene expression correlations) concentrated on the 25 SNPs in the GTEx database, comprising the ten that have been reported to be associated with risk for AMD and the 15 in strong linkage disequilibrium with those ten. We focused especially on the reference SNP rs10490924 because the strong linkage disequilibrium of the other SNPs with the reference SNP meant that results from any of them would approximately predict the results from the others. The present studies of the other SNPs in the haplotype confirmed that prediction.
The three genes in 10q26 closest to the region of the AMD-associated SNPs are (centromeric to telomeric) PLEKHA1, ARMS2, and HTRA1. There is an intergenic region of about 22 kilobases between PLEKHA1 and ARMS2 and an intergenic region of about 5 kilobases between ARMS2 and HTRA1 (Figure 1). The 30-kb AMD-risk haplotype stretches from the intergenic region between the PLEKHA1 and ARMS2 genes to the middle of the HTRA1 gene.
The PLEKHA1 gene is about 60 kb in length and its termination codon, which is at the gene’s telomeric boundary, is about 12 kb from the centromeric boundary of the AMD-risk haplotype. PLEKHA1 is detectably expressed in all 53 tissues in the GTEx database (Table 2 and Appendix 1) with an average expression level of about 9.5 RPKM. In 32 tissues, the average PLEKHA1 expression was lower in people with the high-AMD-risk haplotype defined by rs10490924-T, and the associated p values were below 0.05 in ten of those tissues (Table 2 and Appendix 2). Other SNPs in the AMD haplotype showed similar results, which was expected because of their strong linkage disequilibrium (Appendix 2). The effect was most striking in ovarian tissue (effect size −0.43, p = 3.8 × 10−5 for rs100490924; Figure 2). The effect appears to be semidominant with the expression level in heterozygotes falling between the TT and GG homozygotes. In the nine tissues that showed the opposite pattern (higher expression levels in people with the high-AMD-risk haplotype), the variation in mRNA levels among the genotypes was small and of no statistical significance, with all but one of the p values greater than 0.4. Thus, although the PLEKHA1 gene is more than 22 kb away from rs10490924 and about 12 kb from SNP rs61871744 that defines the centromeric boundary of the AMD-risk haplotype, the expression of this gene in many tissues is influenced by sequences within the haplotype.
The ARMS2 transcriptional unit is about 2 kb and is completely contained within the 30-kb region containing the AMD-risk SNPs. In the GTEx transcriptome data set, ARMS2 mRNA is not detected in 18 tissues and is at a low level in most of the other tissues, averaging only about 0.2 RPKM and never above 1.0 RPKM except testes that have a level of 3.0 RPKM (Table 3 and Appendix 3). In the testes, homozygotes with the rs10490929-T allele (a marker of the high-AMD-risk haplotype) generally have lower ARMS2 mRNA levels than homozygotes with the low-AMD-risk G allele (effect size −0.79, p = 7.5 × 10−24; Figure 3). The effect appears to be semidominant with the expression level in heterozygotes falling between the TT and GG homozygotes. As expected because of strong linkage disequilibrium, another high-AMD-risk allele in the region (rs1120638-A), which is in the same haplotype as the rs10490924-T allele, is also associated with a low ARMS2 mRNA level in the testes (effect size −0.76, p = 8.3 × 10−20; Appendix 2). Eight other AMD-associated variants in this region that we evaluated are associated with a similarly low expression of ARMS2 in the testes (Appendix 2).
Among tissues other than the testes, 34 express detectable but low levels of ARMS2 and thus provide poor statistical power (Table 3). In 32 of the 34 tissues, the ARMS2 transcript levels were numerically lower in homozygotes with the rs10490929-T allele compared with heterozygotes or G-allele homozygotes, a trend that is concordant with the results from the testes (Table 3). The difference in expression between people with the high-AMD-risk T allele versus the low-AMD-risk G allele reached p values less than or equal to 0.05 in 14 of the tissues, such as sun-exposed skin (p = 7.0 × 10−7) and coronary artery (p = 5.5 × 10−6). Analysis of the other SNPs in the AMD-risk haplotype that we evaluated showed a similar pattern (Appendix 2).
The SNP rs2142308, which forms the telomeric boundary of the AMD-risk haplotype, is inside intron 1 of HTRA1. HTRA1 mRNA is ubiquitously expressed with RPKM values ranging from 1 (whole blood) to 269 (the aorta) and averages about 400-fold higher than ARMS2 and 8.5-fold higher than PLEKHA1 (81.2 RPKM versus 0.2 or 9.5 RPKM, respectively; Table 4 and Appendix 4). In 27 tissues, the average HTRA1 expression was higher in people with the high-AMD-risk haplotype defined by rs10490924-T, and the effect was associated with p values of less than 0.05 in six of those tissues (Table 4 and Appendix 2). The effect was most striking in the testes (effect size 0.4, p = 1.5 × 10−7 for rs10490924; Figure 4 and Table 4). In one tissue (prostate), the opposite association was observed with an effect size of −0.23 (p = 0.02). Analysis of the other SNPs in the AMD-risk haplotype that we evaluated showed a similar pattern in HTRA1 mRNA levels in human tissues (Appendix 2). In one tissue (prostate) the opposite association was observed with an effect size of -0.23 (p = 0.02) with rs10490924.
To test whether the high-AMD-risk haplotype could influence the expression of genes even more distant than PLEKHA1, ARMS2, or HTRA1, we examined the transcription of the 22 genes located within 1 megabase of the PLEKHA1-ARMS2-HTRA1 cluster in all 44 human tissues in the GTEx database. Of the 22 genes, 16 were detectably transcribed in at least one tissue. There was no statistically significant association between alleles at the reference SNP rs10490924 and the expression of any of these 16 genes in any of the 44 tissues. In fact, searching through the entire transcriptomes of all 44 tissues in the GTEx database with the criterion of a false discovery rate (FDR) of less than or equal to 0.5, none of the AMD-associated SNPs appeared to be eQTLs for any gene on any chromosome except PLEKHA1, ARMS2, and HTRA1 (Table 5). Thus, the 10q26 high-AMD-risk haplotype appears to affect only PLEKHA1, ARMS2, and HTRA1, and it is unlikely that the high-AMD-risk haplotype has a distant eQTL or inter-chromosome effect on gene expression.
Many more SNPs in 10q26 influence the expression of PLEKHA1 and ARMS2 than HTRA1 (Figure 5A). For example, the expression of PLEKHA1 across 12 tissues is significantly influenced by more than 254 SNPs, with “significance” defined by the GTEx database significance criterion of an FDR of 5% or less. Of the 254 SNPs, 28 also affect the expression of HTRA1, while the remaining 226 SNPs affect the expression of PLEKHA1 but not HTRA1. None of the 226 has been found to influence the risk for AMD. Only 25 of the 254 SNPs that influence the expression of PLEKHA1 also modulate the risk for AMD, and all 25 are among those that modulate HTRA1 expression. The expression of ARMS2 in 13 tissues is influenced by 192 SNPs extending across 855 kb (Figure 5A and Appendix 2). Of the 192 SNPs, 35 also affect the expression of HTRA1, while the remaining 157 affect the expression of ARMS2 but not HTRA1. None of the 157 has been found to influence the risk for AMD. Only 25 of the 192 SNPs that influence the expression of ARMS2 also modulate the risk for AMD, and all 25 are among those that modulate HTRA1 expression. Thus, more than 85% of the SNPs that influence the expression of ARMS2 and PLEKHA1 are not associated with risk for AMD. Ninety-nine of these SNPs have high minor allele frequencies (higher than 0.2), and 35 have a high effect size (absolute value higher than 1; Appendix 2). Some of these SNPs were definitely included in published genome-wide association studies (GWASs). For example, 20 SNPs that affect PLEKHA1 and ARMS2 but not HTRA1 are in the Illumina Human 610-Quad BeadChip used by Fritsch et al. who did not report an association with AMD with any of these SNPs . It is likely that among the numerous, large GWASs conducted on patients with AMD, some of these SNPs would have been identified as influencing risk for AMD if they actually had an effect on risk for AMD.
For HTRA1, only 41 SNPs influence its expression. Thirty-four of these 41 SNPs are concentrated in the 30-kb AMD haplotype region, including all 25 SNPs in the AMD-risk haplotype (Figure 5B). All AMD-risk SNPs that influence the expression of PLEKHA1 or ARMS2 also influence the expression of HTRA1.
The ENCODE database of open chromatin regions includes data from the RPE, a monolayer of cells that plays a key role in the pathogenesis of AMD. Within the 30-kb AMD-risk region, there are three DNase I hypersensitive sites in RPE cells (Table 6). We evaluated whether any of the 34 polymorphisms in the AMD-risk haplotype, including the 25 in the GTEx database and the nine additional polymorphisms in the LDlink database, were within the DNase I hypersensitive sites. The site with the most DNase I hypersensitivity is a 170-bp segment extending from chromosome 10 positions 124,215,021-124,215,190 (Figure 5B). This site is 5.8 kb upstream of the transcription start site of HTRA1 and within the intron of ARMS2. Using the transcription-factor-binding-motif-prediction software HOMER, PROMO, and RegulomeDB, we searched for consensus binding sequences for transcription factors in the open chromatin region that might include AMD-risk SNPs. The results implicated two AMD-risk SNPs, rs36212732 and rs36212733, which are 8 bp and 21 bp, respectively, away from the telomeric boundary of the hypersensitive site (Figure 6). In comparison with the low-AMD-risk haplotype, the high-AMD-risk haplotype loses sites for transcription factors YY1, LHX2, LHX3, NKX6–1, ALX1, and ALX3, and it creates a site for the transcription factor c-MYB (Figure 6). These transcription factors are expressed in human RPE cells [41,42].
The other two open chromatin sites in this region are less than one tenth as sensitive to DNase I (Table 6). One extends from 124,220,906-124,222,950 and contains two AMD-risk SNPs (rs1049331 and rs2293870), but these polymorphisms do not affect any potential transcription factor binding sites according to the HOMER, PROMO, and RegulomeDB programs. The other open chromatin site extends from 124,228,506-124,228,935 and has no AMD-associated SNPs within or nearby.
We used the GTEx genotype-tissue expression database to explore the effects of SNPs that influence the susceptibility for AMD on the expression of nearby genes in the 10q26 ARMS2-HTRA1 region. The AMD-risk alleles at ten SNPs in this region are in strong linkage disequilibrium, and there are 24 additional nearby SNPs or insertion-deletion polymorphisms that are similarly correlated, 15 of which are in the GTEx database. Based on results from 25 of these SNPs in the GTEx database, high-AMD-risk alleles are associated with lower levels of ARMS2 and PLEKHA1 mRNA and higher levels of HTRA1 mRNA than low-AMD-risk alleles in many human tissues. The AMD-risk SNPS do not influence the level of expression of any other gene in this region or anywhere in the human genome with statistical significance after adjustment for the multiple comparisons. Thus, if the risk for AMD conferred by this region is due to variation in gene expression (rather than a change in the transcribed protein), then the risk is due to variation in the expression of one of these three genes.
The GTEx data additionally provide evidence that variations in ARMS2 and PLEKHA1 expression are less likely than HTRA1 to influence risk for AMD. Hundreds of additional SNPs in this region influence the expression of ARMS2 and PLEKHA1 but not HTRA1. None of those additional SNPs has been associated with the risk for AMD in any published human genetics study although some of these SNPs were included in those studies. However, most of the SNPs (25/41) that influence HTRA1 expression are associated with risk for AMD either from direct evidence from GWASs or because the SNPs are highly correlated with directly implicated SNPs. Two other items provide further support for disregarding ARMS2 as an AMD-susceptibility gene: 1) ARMS2 mRNA and protein are expressed at extremely low levels in eye tissues [9,18]; and 2) human genetics studies of AMD indicate that a specific SNP that creates a nonsense mutation (Arg38End) in ARMS2 is associated with low risk for AMD [5,10,43-45]. In short, HTRA1 is the most likely candidate gene for risk for AMD in this region, and if so, elevated expression of HTRA1 likely increases risk for AMD.
The present analysis was based heavily on the idea that the level of expression of a gene on 10q26 determines risk for AMD. An alternative explanation is that there is a change in the primary structure of an encoded protein. Only one such polymorphism is known among the three candidates. It involves the reference SNP rs10490924, which is a missense polymorphism (Ala69Ser) that affects the ARMS2 coding sequence. The Ala69 ARMS2 allele corresponds to low risk for AMD and high ARMS2 expression while the Ser69 allele corresponds to high risk for AMD and low ARMS2 expression. Evidence against this polymorphism as the basis for risk for AMD is as follows. If expression of Ala69-ARMS2 protects against AMD, one would expect that loss of ARMS2 expression would always be associated with elevated risk for AMD. However, a separate ARMS2 variant, the nonsense change Arg38End, would be expected to produce no functional ARMS2, yet the variant has been found to confer low risk for AMD. The possibility remains that the Ser69 allele promotes the development of AMD because of some toxic effect of the Ser69-ARMS2 protein. This mechanism remains a possibility, but we feel it is unlikely because the Ser69 variant is expressed at low levels across all evaluated tissues. In particular, a recent report showed that the Ser69-ARMS2 protein could not be detected in monocytes from patients carrying the homozygous risk for AMD rs10490924-T variant .
Support for HTRA1 as the risk for AMD factor comes from reports of two- to threefold higher HTRA1 expression in eyes with AMD [10,17,47-51], although other groups report no effect [28,45,52-55]. Some support for low HTRA1 expression protecting against AMD comes from patients who lack HTRA1 due to recessive, null mutations. Such patients have cerebral arteries with small lumens and thick walls with reduplicated elastic laminas. No AMD has been observed in such patients .
There are weaknesses in the present analysis. Although the GTEx database includes six tissues with a strong correlation between high-AMD-risk SNP alleles and higher HTRA1 expression level and 21 others have a trend in the same direction, one tissue (prostate) showed a correlation in the opposite direction. It would be ideal to have ocular tissues for evaluation, but unfortunately, no data from ocular tissues are available in the GTEx database. It is still unclear to what extent systemic factors influence one’s risk for AMD . It is conceivable that variation in the systemic expression of HTRA1, not its local ocular expression, is responsible for increasing the risk for AMD.
A potential mechanism for the variation in the expression HTRA1 due to the AMD-risk haplotype involves an open chromatin (DNase I-sensitive) region in the RPE that we found in the ENCODE database. This region is 8–21 bp away from the AMD-risk SNPs rs36212732 and rs36212733. Specifically, the change from the low-AMD-risk allele to the high-AMD-risk allele switches transcription factor binding from YY-1, LHX2, LHX3, ALX1, ALX3, or NKX6–1 to c-MYB (Figure 6).
A recently published evaluation of the ARMS2-HTRA1 region provides additional evidence for the importance of the open chromatin region and the SNPs near it that potentially affect transcription factor binding sites . Based on 33,000 AMD cases and controls, the interval most likely responsible for risk for AMD was narrowed down to a 7136-bp segment bounded by SNPs rs11200630 and esv2663177. This interval is within the AMD haplotype we defined and contains 13 of the 25 SNPs. This interval includes the open chromatin region we uncovered, as well as the SNPs rs36212732 and rs36212733 that affect transcription factor binding sites. It is conceivable that variation in these two SNPs is the fundamental cause for risk for AMD in 10q26 and that the changes in transcription factor affinity mediated by the SNP alleles cause the variation in the expression of the HTRA1 gene located 5.8 kb away.
There are other possible mechanisms for the variation in the expression HTRA1 or other genes due to the AMD-risk haplotype, three of which are the following. 1) AMD-risk SNPs, such as rs10490924, are in strong genetic linkage disequilibrium with the insertion/deletion polymorphism del443ins54 in the 3′ untranslated region of ARMS mRNA . It is possible that the del443ins54 polymorphism introduces a conformational change in chromatin thus affecting the expression of genes in 10q26. However, the recently reported minimal region responsible for risk for AMD does not include this polymorphism, making it unlikely that it modulates risk for AMD . 2) The reported pattern of DNA methylation in the promoter region of ARMS2 correlates with the risk for AMD allele rs10490924-T . However, the variation in methylation would more likely affect the expression of ARMS2 rather than HTRA1 or PLEKHA1. 3) It is possible that transcriptional activity of the ARMS2 or PLEKHA1 genes might interfere with transcription of the nearby HTRA1. Chimeric transcripts starting from PLEKHA1 and ending in ARMS2 were recently reported . The reduction in ARMS2 and PLEKHA1 transcription by the high-AMD-risk variants might consequently allow more HTRA1 mRNA to be transcribed. This indirect effect on HTRA1 gene expression may explain why fewer tissues with increased HTRA1 gene expression reached statistical significance compared to ARMS2 in the GTEx database analysis.
Appendix 1. PLEKHA1 mRNA level is expressed in all tissues evaluated in the GTEx project.
Appendix 2. Maximum effect sizes for PLEKHA1, ARMS2, and HTRA1 cis-eQTLs.
Appendix 3. ARMS2 mRNA level is highest in testis compared to other tissues in the GTEx project.
Appendix 4. HTRA1 mRNA level is expressed in all tissues evaluated in the GTEx project.
The authors thank Eric Marshall for assisting with the GTEx data analysis. The GTEx website asks for the following acknowledgment: “The Genotype-Tissue Expression (GTEx) Project was supported by the Common Fund of the Office of the Director of the National Institutes of Health. Additional funds were provided by the NCI, NHGRI, NHLBI, NIDA, NIMH, and NINDS. Donors were enrolled at Biospecimen Source Sites funded by NCI\SAIC-Frederick, Inc. (SAIC-F) subcontracts to the National Disease Research Interchange (10XS170), Roswell Park Cancer Institute (10XS171), and Science Care, Inc. (X10S172). The Laboratory, Data Analysis, and Coordinating Center (LDACC) was funded through a contract (HHSN268201000029C) to The Broad Institute, Inc. Biorepository operations were funded through an SAIC-F subcontract to Van Andel Institute (10ST1035). Additional data repository and project management were provided by SAIC-F (HHSN261200800001E). The Brain Bank was supported by supplements to University of Miami grants DA006227 & DA033684 and to contract N01MH000028. Statistical Methods development grants were made to the University of Geneva (MH090941 & MH101814), the University of Chicago (MH090951, MH090937, MH101820, MH101825), the University of North Carolina - Chapel Hill (MH090936 & MH101819), Harvard University (MH090948), Stanford University (MH101782), Washington University St Louis (MH101810), and the University of Pennsylvania (MH101822). The data used for the analyses described in this manuscript were obtained from the GTEx Portal from 11/01/2015 to 12/04/2016. Parts of the data (Figure 2-4 and Table 2-4) have been presented at ARVO 2016. All work was performed at Novartis.