Identification and single-base gene-editing functional validation of a cis-EPO variant as a genetic predictor for EPO-increasing therapies.

Hypoxia-inducible factor prolyl hydroxylase inhibitors (HIF-PHIs) are currently under clinical development for treating anemia in chronic kidney disease (CKD), but it is important to monitor their cardiovascular safety. Genetic variants can be used as predictors to help inform the potential risk of adverse effects associated with drug treatments. We therefore aimed to use human genetics to help assess the risk of adverse cardiovascular events associated with therapeutically altered EPO levels to help inform clinical trials studying the safety of HIF-PHIs. By performing a genome-wide association meta-analysis of EPO (n = 6,127), we identified a cis-EPO variant (rs1617640) lying in the EPO promoter region. We validated this variant as most likely causal in controlling EPO levels by using genetic and functional approaches, including single-base gene editing. Using this variant as a partial predictor for therapeutic modulation of EPO and large genome-wide association data in Mendelian randomization tests, we found no evidence (at p < 0.05) that genetically predicted long-term rises in endogenous EPO, equivalent to a 2.2-unit increase, increased risk of coronary artery disease (CAD, OR [95% CI] = 1.01 [0.93, 1.07]), myocardial infarction (MI, OR [95% CI] = 0.99 [0.87, 1.15]), or stroke (OR [95% CI] = 0.97 [0.87, 1.07]). We could exclude increased odds of 1.15 for cardiovascular disease for a 2.2-unit EPO increase. A combination of genetic and functional studies provides a powerful approach to investigate the potential therapeutic profile of EPO-increasing therapies for treating anemia in CKD.

Genetic variation in severe cystic fibrosis liver disease is associated with novel mechanisms for disease pathogenesis.

BACKGROUND AND AIMS: It is not known why severe cystic fibrosis (CF) liver disease (CFLD) with portal hypertension occurs in only ~7% of people with CF. We aimed to identify genetic modifiers for severe CFLD to improve understanding of disease mechanisms. APPROACH AND RESULTS: Whole-genome sequencing was available in 4082 people with CF with pancreatic insufficiency (n = 516 with severe CFLD; n = 3566 without CFLD). We tested ~15.9 million single nucleotide polymorphisms (SNPs) for association with severe CFLD versus no-CFLD, using pre-modulator clinical phenotypes including (1) genetic variant ( SERPINA1 ; Z allele) previously associated with severe CFLD; (2) candidate SNPs (n = 205) associated with non-CF liver diseases; (3) genome-wide association study of common/rare SNPs; (4) transcriptome-wide association; and (5) gene-level and pathway analyses. The Z allele was significantly associated with severe CFLD ( p = 1.1 × 10 -4 ). No significant candidate SNPs were identified. A genome-wide association study identified genome-wide significant SNPs in 2 loci and 2 suggestive loci. These 4 loci contained genes [significant, PKD1 ( p = 8.05 × 10 -10 ) and FNBP1 ( p = 4.74 × 10 -9 ); suggestive, DUSP6 ( p = 1.51 × 10 -7 ) and ANKUB1 ( p = 4.69 × 10 -7 )] relevant to severe CFLD pathophysiology. The transcriptome-wide association identified 3 genes [ CXCR1 ( p = 1.01 × 10 -6 ) , AAMP ( p = 1.07 × 10 -6 ), and TRBV24 ( p = 1.23 × 10 -5 )] involved in hepatic inflammation and innate immunity. Gene-ranked analyses identified pathways enriched in genes linked to multiple liver pathologies. CONCLUSION: These results identify loci/genes associated with severe CFLD that point to disease mechanisms involving hepatic fibrosis, inflammation, innate immune function, vascular pathology, intracellular signaling, actin cytoskeleton and tight junction integrity and mechanisms of hepatic steatosis and insulin resistance. These discoveries will facilitate mechanistic studies and the development of therapeutics for severe CFLD.

Genetic Modifiers of Cystic Fibrosis Lung Disease Severity: Whole-Genome Analysis of 7,840 Patients.

Rationale: Lung disease is the major cause of morbidity and mortality in persons with cystic fibrosis (pwCF). Variability in CF lung disease has substantial non-CFTR (CF transmembrane conductance regulator) genetic influence. Identification of genetic modifiers has prognostic and therapeutic importance. Objectives: Identify genetic modifier loci and genes/pathways associated with pulmonary disease severity. Methods: Whole-genome sequencing data on 4,248 unique pwCF with pancreatic insufficiency and lung function measures were combined with imputed genotypes from an additional 3,592 patients with pancreatic insufficiency from the United States, Canada, and France. This report describes association of approximately 15.9 million SNPs using the quantitative Kulich normal residual mortality-adjusted (KNoRMA) lung disease phenotype in 7,840 pwCF using premodulator lung function data. Measurements and Main Results: Testing included common and rare SNPs, transcriptome-wide association, gene-level, and pathway analyses. Pathway analyses identified novel associations with genes that have key roles in organ development, and we hypothesize that these genes may relate to dysanapsis and/or variability in lung repair. Results confirmed and extended previous genome-wide association study findings. These whole-genome sequencing data provide finely mapped genetic information to support mechanistic studies. No novel primary associations with common single variants or rare variants were found. Multilocus effects at chr5p13 (SLC9A3/CEP72) and chr11p13 (EHF/APIP) were identified. Variant effect size estimates at associated loci were consistently ordered across the cohorts, indicating possible age or birth cohort effects. Conclusions: This premodulator genomic, transcriptomic, and pathway association study of 7,840 pwCF will facilitate mechanistic and postmodulator genetic studies and the development of novel therapeutics for CF lung disease.

Airway Mucosal Host Defense Is Key to Genomic Regulation of Cystic Fibrosis Lung Disease Severity.

RATIONALE: The severity of cystic fibrosis (CF) lung disease varies widely, even for Phe508del homozygotes. Heritability studies show that more than 50% of the variability reflects non-cystic fibrosis transmembrane conductance regulator (CFTR) genetic variation; however, the full extent of the pertinent genetic variation is not known. OBJECTIVES: We sought to identify novel CF disease-modifying mechanisms using an integrated approach based on analyzing "in vivo" CF airway epithelial gene expression complemented with genome-wide association study (GWAS) data. METHODS: Nasal mucosal RNA from 134 patients with CF was used for RNA sequencing. We tested for associations of transcriptomic (gene expression) data with a quantitative phenotype of CF lung disease severity. Pathway analysis of CF GWAS data (n = 5,659 patients) was performed to identify novel pathways and assess the concordance of genomic and transcriptomic data. Association of gene expression with previously identified CF GWAS risk alleles was also tested. MEASUREMENTS AND MAIN RESULTS: Significant evidence of heritable gene expression was identified. Gene expression pathways relevant to airway mucosal host defense were significantly associated with CF lung disease severity, including viral infection, inflammation/inflammatory signaling, lipid metabolism, apoptosis, ion transport, Phe508del CFTR processing, and innate immune responses, including HLA (human leukocyte antigen) genes. Ion transport and CFTR processing pathways, as well as HLA genes, were identified across differential gene expression and GWAS signals. CONCLUSIONS: Transcriptomic analyses of CF airway epithelia, coupled to genomic (GWAS) analyses, highlight the role of heritable host defense variation in determining the pathophysiology of CF lung disease. The identification of these pathways provides opportunities to pursue targeted interventions to improve CF lung health.

Genomic Characterization of Metformin Hepatic Response.

Metformin is used as a first-line therapy for type 2 diabetes (T2D) and prescribed for numerous other diseases. However, its mechanism of action in the liver has yet to be characterized in a systematic manner. To comprehensively identify genes and regulatory elements associated with metformin treatment, we carried out RNA-seq and ChIP-seq (H3K27ac, H3K27me3) on primary human hepatocytes from the same donor treated with vehicle control, metformin or metformin and compound C, an AMP-activated protein kinase (AMPK) inhibitor (allowing to identify AMPK-independent pathways). We identified thousands of metformin responsive AMPK-dependent and AMPK-independent differentially expressed genes and regulatory elements. We functionally validated several elements for metformin-induced promoter and enhancer activity. These include an enhancer in an ataxia telangiectasia mutated (ATM) intron that has SNPs in linkage disequilibrium with a metformin treatment response GWAS lead SNP (rs11212617) that showed increased enhancer activity for the associated haplotype. Expression quantitative trait locus (eQTL) liver analysis and CRISPR activation suggest that this enhancer could be regulating ATM, which has a known role in AMPK activation, and potentially also EXPH5 and DDX10, its neighboring genes. Using ChIP-seq and siRNA knockdown, we further show that activating transcription factor 3 (ATF3), our top metformin upregulated AMPK-dependent gene, could have an important role in gluconeogenesis repression. Our findings provide a genome-wide representation of metformin hepatic response, highlight important sequences that could be associated with interindividual variability in glycemic response to metformin and identify novel T2D treatment candidates.

Gene expression in transformed lymphocytes reveals variation in endomembrane and HLA pathways modifying cystic fibrosis pulmonary phenotypes.

Variation in cystic fibrosis (CF) phenotypes, including lung disease severity, age of onset of persistent Pseudomonas aeruginosa (P. aeruginosa) lung infection, and presence of meconium ileus (MI), has been partially explained by genome-wide association studies (GWASs). It is not expected that GWASs alone are sufficiently powered to uncover all heritable traits associated with CF phenotypic diversity. Therefore, we utilized gene expression association from lymphoblastoid cells lines from 754 p.Phe508del CF-affected homozygous individuals to identify genes and pathways. LPAR6, a G protein coupled receptor, associated with lung disease severity (false discovery rate q value = 0.0006). Additional pathway analyses, utilizing a stringent permutation-based approach, identified unique signals for all three phenotypes. Pathways associated with lung disease severity were annotated in three broad categories: (1) endomembrane function, containing p.Phe508del processing genes, providing evidence of the importance of p.Phe508del processing to explain lung phenotype variation; (2) HLA class I genes, extending previous GWAS findings in the HLA region; and (3) endoplasmic reticulum stress response genes. Expression pathways associated with lung disease were concordant for some endosome and HLA pathways, with pathways identified using GWAS associations from 1,978 CF-affected individuals. Pathways associated with age of onset of persistent P. aeruginosa infection were enriched for HLA class II genes, and those associated with MI were related to oxidative phosphorylation. Formal testing demonstrated that genes showing differential expression associated with lung disease severity were enriched for heritable genetic variation and expression quantitative traits. Gene expression provided a powerful tool to identify unrecognized heritable variation, complementing ongoing GWASs in this rare disease.

Polymorphisms associated with expression of BPIFA1/BPIFB1 and lung disease severity in cystic fibrosis.

BPI fold containing family A, member 1 (BPIFA1) and BPIFB1 are putative innate immune molecules expressed in the upper airways. Because of their hypothesized roles in airway defense, these molecules may contribute to lung disease severity in cystic fibrosis (CF). We interrogated BPIFA1/BPIFB1 single-nucleotide polymorphisms in data from an association study of CF modifier genes and found an association of the G allele of rs1078761 with increased lung disease severity (P = 2.71 × 10(-4)). We hypothesized that the G allele of rs1078761 is associated with decreased expression of BPIFA1 and/or BPIFB1. Genome-wide lung gene expression and genotyping data from 1,111 individuals with lung disease, including 51 patients with CF, were tested for associations between genotype and BPIFA1 and BPIFB1 gene expression levels. Findings were validated by quantitative PCR in a subset of 77 individuals. Western blotting was used to measure BPIFA1 and BPIFB1 protein levels in 93 lung and 101 saliva samples. The G allele of rs1078761 was significantly associated with decreased mRNA levels of BPIFA1 (P = 4.08 × 10(-15)) and BPIFB1 (P = 0.0314). These findings were confirmed with quantitative PCR and Western blotting. We conclude that the G allele of rs1078761 may be detrimental to lung function in CF owing to decreased levels of BPIFA1 and BPIFB1.

Dementia revealed: novel chromosome 6 locus for late-onset Alzheimer disease provides genetic evidence for folate-pathway abnormalities.

Genome-wide association studies (GWAS) of late-onset Alzheimer disease (LOAD) have consistently observed strong evidence of association with polymorphisms in APOE. However, until recently, variants at few other loci with statistically significant associations have replicated across studies. The present study combines data on 483,399 single nucleotide polymorphisms (SNPs) from a previously reported GWAS of 492 LOAD cases and 496 controls and from an independent set of 439 LOAD cases and 608 controls to strengthen power to identify novel genetic association signals. Associations exceeding the experiment-wide significance threshold (alpha=1.03x10(-7)) were replicated in an additional 1,338 cases and 2,003 controls. As expected, these analyses unequivocally confirmed APOE's risk effect (rs2075650, P=1.9x10(-36)). Additionally, the SNP rs11754661 at 151.2 Mb of chromosome 6q25.1 in the gene MTHFD1L (which encodes the methylenetetrahydrofolate dehydrogenase (NADP+ dependent) 1-like protein) was significantly associated with LOAD (P=4.70x10(-8); Bonferroni-corrected P=0.022). Subsequent genotyping of SNPs in high linkage disequilibrium (r2>0.8) with rs11754661 identified statistically significant associations in multiple SNPs (rs803424, P=0.016; rs2073067, P=0.03; rs2072064, P=0.035), reducing the likelihood of association due to genotyping error. In the replication case-control set, we observed an association of rs11754661 in the same direction as the previous association at P=0.002 (P=1.90x10(-10) in combined analysis of discovery and replication sets), with associations of similar statistical significance at several adjacent SNPs (rs17349743, P=0.005; rs803422, P=0.004). In summary, we observed and replicated a novel statistically significant association in MTHFD1L, a gene involved in the tetrahydrofolate synthesis pathway. This finding is noteworthy, as MTHFD1L may play a role in the generation of methionine from homocysteine and influence homocysteine-related pathways and as levels of homocysteine are a significant risk factor for LOAD development.

Mitochondrial haplogroup X is associated with successful aging in the Amish.

Avoiding disease, maintaining physical and cognitive function, and continued social engagement in long-lived individuals describe successful aging (SA). Mitochondrial lineages described by patterns of common genetic variants ("haplogroups") have been associated with increased longevity in different populations. We investigated the influence of mitochondrial haplogroups on SA in an Amish community sample. Cognitively intact volunteers aged ≥80 years (n = 261) were enrolled in a door-to-door survey of Amish communities in Indiana and Ohio. Individuals scoring in the top third for lower extremity function, needing little assistance with self-care tasks, having no depression symptoms, and expressing high life satisfaction were considered SA (n = 74). The remainder (n = 187) were retained as controls. These individuals descend from 51 matrilines in a single 13-generation pedigree. Mitochondrial haplogroups were assigned using the ten mitochondrial single nucleotide polymorphisms (mtSNPs) defining the nine most common European haplogroups. An additional 17 mtSNPs from a genome-wide association panel were also investigated. Associations between haplogroups, mtSNPs, and SA were determined by logistic regression models accounting for sex, age, body mass index, and matriline via generalized estimating equations. SA cases were more likely to carry Haplogroup X (OR = 7.56, p = 0.0015), and less likely to carry Haplogroup J (OR = 0.40, p = 0.0003). Our results represent a novel association of Haplogroup X with SA and suggest that variants in the mitochondrial genome may promote maintenance of both physical and cognitive function in older adults.

The ARMS2 A69S variant and bilateral advanced age-related macular degeneration.

PURPOSE: To identify genetic associations between specific risk genes and bilateral advanced age-related macular degeneration (AMD) in a retrospective, observational case series of 1,003 patients: 173 patients with geographic atrophy in at least 1 eye and 830 patients with choroidal neovascularization in at least 1 eye. METHODS: Patients underwent clinical examination and fundus photography. The images were subsequently graded using a modified grading system adapted from the Age-Related Eye Disease Study. Genetic analysis was performed to identify genotypes at 4 AMD-associated variants (ARMS2 A69S, CFH Y402H, C3 R102G, and CFB R32Q) in these patients. RESULTS: There were no statistically significant relationships between clinical findings and genotypes at CFH, C3, and CFB. The genotype at ARMS2 correlated with bilateral advanced AMD using a variety of comparisons: unilateral geographic atrophy versus bilateral geographic atrophy (P = 0.08), unilateral choroidal neovascularization versus bilateral choroidal neovascularization (P = 9.0 × 10(-8)), and unilateral late AMD versus bilateral late AMD (P = 5.9 × 10(-8)). CONCLUSION: In this series, in patients with geographic atrophy or choroidal neovascularization in at least 1 eye, the ARMS2 A69S substitution strongly associated with geographic atrophy or choroidal neovascularization in the fellow eye. The ARMS2 A69S substitution may serve as a marker for bilateral advanced AMD.

A Population-Based Human In Vitro Approach to Quantify Inter-Individual Variability in Responses to Chemical Mixtures.

Human cell-based population-wide in vitro models have been proposed as a strategy to derive chemical-specific estimates of inter-individual variability; however, the utility of this approach has not yet been tested for cumulative exposures in mixtures. This study aimed to test defined mixtures and their individual components and determine whether adverse effects of the mixtures were likely to be more variable in a population than those of the individual chemicals. The in vitro model comprised 146 human lymphoblastoid cell lines from four diverse subpopulations of European and African descent. Cells were exposed, in concentration−response, to 42 chemicals from diverse classes of environmental pollutants; in addition, eight defined mixtures were prepared from these chemicals using several exposure- or hazard-based scenarios. Points of departure for cytotoxicity were derived using Bayesian concentration−response modeling and population variability was quantified in the form of a toxicodynamic variability factor (TDVF). We found that 28 chemicals and all mixtures exhibited concentration−response cytotoxicity, enabling calculation of the TDVF. The median TDVF across test substances, for both individual chemicals or defined mixtures, ranged from a default assumption (101/2) of toxicodynamic variability in human population to >10. The data also provide a proof of principle for single-variant genome-wide association mapping for toxicity of the chemicals and mixtures, although replication would be necessary due to statistical power limitations with the current sample size. This study demonstrates the feasibility of using a set of human lymphoblastoid cell lines as an in vitro model to quantify the extent of inter-individual variability in hazardous properties of both individual chemicals and mixtures. The data show that population variability of the mixtures is unlikely to exceed that of the most variable component, and that similarity in genome-wide associations among components may be used to accrue additional evidence for grouping of constituents in a mixture for cumulative assessments.

A resource for integrated genomic analysis of the human liver.

In this study, we generated whole-transcriptome RNA-Seq from n = 192 genotyped liver samples and used these data with existing data from the GTEx Project (RNA-Seq) and previous liver eQTL (microarray) studies to create an enhanced transcriptomic sequence resource in the human liver. Analyses of genotype-expression associations show pronounced enrichment of associations with genes of drug response. The associations are primarily consistent across the two RNA-Seq datasets, with some modest variation, indicating the importance of obtaining multiple datasets to produce a robust resource. We further used an empirical Bayesian model to compare eQTL patterns in liver and an additional 20 GTEx tissues, finding that MHC genes, and especially class II genes, are enriched for liver-specific eQTL patterns. To illustrate the utility of the resource to augment GWAS analysis with small sample sizes, we developed a novel meta-analysis technique to combine several liver eQTL data sources. We also illustrate its application using a transcriptome-enhanced re-analysis of a study of neutropenia in pancreatic cancer patients. The associations of genotype with liver expression, including splice variation and its genetic associations, are made available in a searchable genome browser.

Leveraging TOPMed imputation server and constructing a cohort-specific imputation reference panel to enhance genotype imputation among cystic fibrosis patients.

Cystic fibrosis (CF) is a severe genetic disorder that can cause multiple comorbidities affecting the lungs, the pancreas, the luminal digestive system and beyond. In our previous genome-wide association studies (GWAS), we genotyped approximately 8,000 CF samples using a mixture of different genotyping platforms. More recently, the Cystic Fibrosis Genome Project (CFGP) performed deep (approximately 30×) whole genome sequencing (WGS) of 5,095 samples to better understand the genetic mechanisms underlying clinical heterogeneity among patients with CF. For mixtures of GWAS array and WGS data, genotype imputation has proven effective in increasing effective sample size. Therefore, we first performed imputation for the approximately 8,000 CF samples with GWAS array genotype using the Trans-Omics for Precision Medicine (TOPMed) freeze 8 reference panel. Our results demonstrate that TOPMed can provide high-quality imputation for patients with CF, boosting genomic coverage from approximately 0.3-4.2 million genotyped markers to approximately 11-43 million well-imputed markers, and significantly improving polygenic risk score (PRS) prediction accuracy. Furthermore, we built a CF-specific CFGP reference panel based on WGS data of patients with CF. We demonstrate that despite having approximately 3% the sample size of TOPMed, our CFGP reference panel can still outperform TOPMed when imputing some CF disease-causing variants, likely owing to allele and haplotype differences between patients with CF and general populations. We anticipate our imputed data for 4,656 samples without WGS data will benefit our subsequent genetic association studies, and the CFGP reference panel built from CF WGS samples will benefit other investigators studying CF.

Successful aging shows linkage to chromosomes 6, 7, and 14 in the Amish.

Successful aging (SA) is a multidimensional phenotype involving preservation of cognitive ability, physical function, and social engagement throughout life. Multiple components of SA are heritable, supporting a genetic component. The Amish are genetically and socially isolated with homogeneous lifestyles, making them a suitable population for studying the genetics of SA. DNA and measures of SA were collected on 214 cognitively intact Amish individuals over age 80. Individuals were grouped into a 13-generation pedigree using the Anabaptist Genealogy Database. A linkage screen of 5944 single nucleotide polymorphisms (SNPs) was performed using 12 informative subpedigrees with an affected-only 2-point and multipoint linkage analysis. Eleven SNPs produced 2-point LOD scores >2, suggestive of linkage. Multipoint linkage analyses, allowing for heterogeneity, detected significant LOD scores on chromosomes 6 (HLOD = 4.50), 7 (LOD*= 3.11), and 14 (HLOD = 4.17), suggesting multiple new loci underlying SA.

A genome-wide linkage screen in the Amish with Parkinson disease points to chromosome 6.

Parkinson disease (PD) is a common complex neurodegenerative disorder with an underlying genetic etiology that has been difficult to dissect. Although some PD risk genes have been discovered, most of the underlying genetic etiology remains unknown. To further elucidate the genetic component, we have undertaken a genome-wide linkage screen in an isolated founder population of Amish living in the Midwestern United States. We performed tests for linkage and for association using a marker set of nearly 6000 single-nucleotide polymorphisms. Parametric multipoint linkage analysis generated a logarithm of the odds of linkage (LOD) score of 2.44 on chromosome 6 in the SYNE1 gene, approximately 8 Mbp from the PARK2 gene. In a different region on chromosome 6 (∼67 Mbp from PARK2) an association was found for rs4302647 (p < 4.0 × 10(-6) ), which is not within 300 kb of any gene. While the association exceeds Bonferroni correction, it may yet represent a false positive due to the small sample size and the low minor allele frequency. The minor allele frequency in affecteds is 0.07 compared to 0.01 in unaffecteds. Taken together, these results support involvement of loci on chromosome 6 in the genetic etiology of PD.

Deletion of CFHR3 and CFHR1 genes in age-related macular degeneration.

Age-related macular degeneration (AMD) impairs vision for approximately 7.5 million Americans. Both susceptibility variants and protective haplotypes in the complement factor H (CFH) gene modulate risk for AMD. Recently, deletion of the 'CFH-related' genes CFHR1 and CFHR3 was found to be segregating with a particular CFH haplotype, which reduced the risk of AMD. We tested the deletion for association in a Caucasian population of 780 cases and 265 controls and examined its effect in the context of known AMD risk factors. The deletion did not segregate perfectly with any one SNP, as previously suggested. CFH haplotype P2 was the most frequent haplotype in deletion homozygotes (47%), and the majority (14/16) of these individuals were homozygous for the non-risk allele of Y402H. Overall, deletion homozygosity was significantly more frequent in controls than cases (2.6% controls, 0.8% cases, P = 0.025, OR = 0.29, 95% CI = 0.10-0.86). After controlling for age, Y402H, smoking and LOC387715 A69S, the protective effect of the deletion was no longer statistically significant (P = 0.27). However, using a CFH haplotype that all deletion homozygotes share as a surrogate for the deletion, this marker remained modestly associated with AMD after adjustment for known risk factors (OR = 0.63, 95% CI 0.39-1.04, P = 0.07). Therefore, deletion of CFHR1 and CFHR3 may account for a small portion of the protection from AMD associated with particular haplotypes in CFH. The presence of protective haplotypes in CFH that do not carry the deletion, suggests that other protective variants in this region have yet to be discovered.

C3 R102G polymorphism increases risk of age-related macular degeneration.

Inflammation has long been suspected to play a role in the pathogenesis of age-related macular degeneration (AMD). Association of variants in the complement factor H (CFH) and complement factor B (CFB) genes has targeted the search for additional loci to the alternative complement cascade, of which C3 is a major component. Two non-synonymous coding polymorphisms within C3, R102G and L314P, have previously been strongly associated with increased risk. These variants are in strong linkage disequilibrium (LD), making the contribution of this locus to AMD even more difficult to ascertain. We sought to determine whether the C3 association resulted primarily from only one of these two variants or from a combined effect of both in 223 families and an independent dataset of 701 cases and 286 unrelated controls. The C3 polymorphisms were in strong LD (r(2) = 0.85), and both were associated in the family-based and case-control datasets (R102G genoPDT P = 0.02, case-control genotypic P = 0.004; L314P genoPDT P = 0.001, case-control genotypic P = 0.04). In conditional analyses in the case-control dataset, R102G remained associated with disease in the L314P risk allele carriers (P = 0.01), but there was no effect of L314P in the R102G risk allele carriers (P = 0.2). After adjusting for age, smoking, CFH Y402H, LOC387715 A69S, and CFB R32Q, the effect of R102G remained strong [P = 0.015, odds ratio = 1.55, 95% confidence interval 1.09 to 2.21, adjusted PAR(population attributable risk) = 0.17]. Therefore, while the strong LD between R102G and L314P makes it difficult to disentangle their individual effects on disease risk, the R102G polymorphism acting alone provides the best model for disease in our data.

Analysis of single nucleotide polymorphisms in the NOS2A gene and interaction with smoking in age-related macular degeneration.

Age-related macular degeneration (AMD) is a complex degenerative retinal disease influenced by both genetic and environmental risk factors. We assessed whether single nucleotide polymorphisms (SNPs) in the NOS2A gene increase risk and modulate the effect of smoking in AMD. 998 Caucasian subjects (712 AMD cases and 286 controls) were genotyped for 17 SNPs in NOS2A. Multivariable logistic regression models containing SNP genotypes, age, sex, smoking status and genotype/smoking interaction were constructed. SNP rs8072199 was significantly associated with AMD (OR = 1.3; 95% CI : 1.02, 1.65; P = 0.035). A significant interaction with smoking was detected at rs2248814 (P = 0.037). Stratified data by genotypes demonstrated that the association between AMD and smoking was stronger in carriers of AA genotypes (OR = 35.98; 95% CI: 3.19, 405.98) than in carriers of the AG genotype (OR = 3.05; 95% CI: 1.36, 6.74) or GG genotype (OR = 2.1; 95% CI: 0.91, 4.84). The results suggest a possible synergistic interaction of AA genotype with smoking, although the result bears replication in larger samples. Our data suggests that SNPs in the NOS2A gene are associated with increased risk for AMD and might modulate the effect of smoking on AMD.

Analysis of the indel at the ARMS2 3'UTR in age-related macular degeneration.

Controversy remains as to which gene at the chromosome 10q26 locus confers risk for age-related macular degeneration (AMD) and statistical genetic analysis is confounded by the strong linkage disequilibrium (LD) across the region. Functional analysis of related genetic variations could solve this puzzle. Recently, Fritsche et al. reported that AMD is associated with unstable ARMS2 transcripts possibly caused by a complex insertion/deletion (indel; consisting of a 443 bp deletion and an adjacent 54 bp insertion) in its 3'UTR (untranslated region). To validate this indel, we sequenced our samples. We found that this indel is even more complex and is composed of two side-by-side indels separated by 17 bp: (1) 9 bp deletion with 10 bp insertion; (2) 417 bp deletion with 27 bp insertion. The indel is significantly associated with the risk of AMD, but is also in strong LD with the non-synonymous single nucleotide polymorphism rs10490924 (A69S). We also found that ARMS2 is expressed not only in placenta and retina but also in multiple human tissues. Using quantitative PCR, we found no correlation between the indel and ARMS2 mRNA level in human retina and blood samples. The lack of functional effects of the 3'UTR indel, the amino acid substitution of rs10490924 (A69S), and strong LD between them suggest that A69S, not the indel, is the variant that confers risk of AMD. To our knowledge, it is the first time it has been shown that ARMS2 is widely expressed in human tissues. Conclusively, the indel at 3'UTR of ARMS2 actually contains two side-by-side indels. The indels are associated with risk of AMD, but not correlated with ARMS2 mRNA level.

Exploring the Limits of Combined Image/'omics Analysis for Non-cancer Histological Phenotypes.

The last several years have witnessed an explosion of methods and applications for combining image data with 'omics data, and for prediction of clinical phenotypes. Much of this research has focused on cancer histology, for which genetic perturbations are large, and the signal to noise ratio is high. Related research on chronic, complex diseases is limited by tissue sample availability, lower genomic signal strength, and the less extreme and tissue-specific nature of intermediate histological phenotypes. Data from the GTEx Consortium provides a unique opportunity to investigate the connections among phenotypic histological variation, imaging data, and 'omics profiling, from multiple tissue-specific phenotypes at the sub-clinical level. Investigating histological designations in multiple tissues, we survey the evidence for genomic association and prediction of histology, and use the results to test the limits of prediction accuracy using machine learning methods applied to the imaging data, genomics data, and their combination. We find that expression data has similar or superior accuracy for pathology prediction as our use of imaging data, despite the fact that pathological determination is made from the images themselves. A variety of machine learning methods have similar performance, while network embedding methods offer at best limited improvements. These observations hold across a range of tissues and predictor types. The results are supportive of the use of genomic measurements for prediction, and in using the same target tissue in which pathological phenotyping has been performed. Although this last finding is sensible, to our knowledge our study is the first to demonstrate this fact empirically. Even while prediction accuracy remains a challenge, the results show clear evidence of pathway and tissue-specific biology.