Multi-ancestry GWAS analysis identifies two novel loci associated with diabetic eye disease and highlights APOL1 as a high risk locus in patients with diabetic macular edema.

Diabetic retinopathy (DR) is a common complication of diabetes. Approximately 20% of DR patients have diabetic macular edema (DME) characterized by fluid leakage into the retina. There is a genetic component to DR and DME risk, but few replicable loci. Because not all DR cases have DME, we focused on DME to increase power, and conducted a multi-ancestry GWAS to assess DME risk in a total of 1,502 DME patients and 5,603 non-DME controls in discovery and replication datasets. Two loci reached GWAS significance (p<5x10-8). The strongest association was rs2239785, (K150E) in APOL1. The second finding was rs10402468, which co-localized to PLVAP and ANKLE1 in vascular / endothelium tissues. We conducted multiple sensitivity analyses to establish that the associations were specific to DME status and did not reflect diabetes status or other diabetic complications. Here we report two novel loci for risk of DME which replicated in multiple clinical trial and biobank derived datasets. One of these loci, containing the gene APOL1, is a risk factor in African American DME and DKD patients, indicating that this locus plays a broader role in diabetic complications for multiple ancestries. Trial Registration: NCT00473330, NCT00473382, NCT03622580, NCT03622593, NCT04108156.

Rare CIDEC coding variants enriched in age-related macular degeneration patients with small low-luminance deficit cause lipid droplet and fat storage defects.

BACKGROUND: The basis of Age-related macular degeneration (AMD) genetic risk has been well documented; however, few studies have looked at genetic biomarkers of disease progression or treatment response within advanced AMD patients. Here we report the first genome-wide analysis of genetic determinants of low-luminance vision deficit (LLD), which is seen as predictive of visual acuity loss and anti-VEGF treatment response in neovascular AMD patients. METHODS: AMD patients were separated into small- and large-LLD groups for comparison and whole genome sequencing was performed. Genetic determinants of LLD were assessed by common and rare variant genetic analysis. Follow-up functional analysis of rare coding variants identified by the burden test was then performed in vitro. RESULTS: We identified four coding variants in the CIDEC gene. These rare variants were only present in patients with a small LLD, which has been previously shown to indicate better prognosis and better anti-VEGF treatment response. Our in vitro functional characterization of these CIDEC alleles revealed that all decrease the binding affinity between CIDEC and the lipid droplet fusion effectors PLIN1, RAB8A and AS160. The rare CIDEC alleles all cause a hypomorphic defect in lipid droplet fusion and enlargement, resulting in a decreased fat storage capability in adipocytes. CONCLUSIONS: As we did not detect CIDEC expression in the ocular tissue affected by AMD, our results suggest that the CIDEC variants do not play a direct role in the eye and influence low-luminance vision deficit via an indirect and systemic effect related to fat storage capacity.

PCSK6 and Survival in Idiopathic Pulmonary Fibrosis.

Rationale: Idiopathic pulmonary fibrosis (IPF) is a devastating disease characterized by limited treatment options and high mortality. A better understanding of the molecular drivers of IPF progression is needed. Objectives: To identify and validate molecular determinants of IPF survival. Methods: A staged genome-wide association study was performed using paired genomic and survival data. Stage I cases were drawn from centers across the United States and Europe and stage II cases from Vanderbilt University. Cox proportional hazards regression was used to identify gene variants associated with differential transplantation-free survival (TFS). Stage I variants with nominal significance (P < 5 × 10-5) were advanced for stage II testing and meta-analyzed to identify those reaching genome-wide significance (P < 5 × 10-8). Downstream analyses were performed for genes and proteins associated with variants reaching genome-wide significance. Measurements and Main Results: After quality controls, 1,481 stage I cases and 397 stage II cases were included in the analysis. After filtering, 9,075,629 variants were tested in stage I, with 158 meeting advancement criteria. Four variants associated with TFS with consistent effect direction were identified in stage II, including one in an intron of PCSK6 (proprotein convertase subtilisin/kexin type 6) reaching genome-wide significance (hazard ratio, 4.11 [95% confidence interval, 2.54-6.67]; P = 9.45 × 10-9). PCSK6 protein was highly expressed in IPF lung parenchyma. PCSK6 lung staining intensity, peripheral blood gene expression, and plasma concentration were associated with reduced TFS. Conclusions: We identified four novel variants associated with IPF survival, including one in PCSK6 that reached genome-wide significance. Downstream analyses suggested that PCSK6 protein plays a potentially important role in IPF progression.

Genome-wide association study across five cohorts identifies five novel loci associated with idiopathic pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF) is a chronic lung condition with poor survival times. We previously published a genome-wide meta-analysis of IPF risk across three studies with independent replication of associated variants in two additional studies. To maximise power and to generate more accurate effect size estimates, we performed a genome-wide meta-analysis across all five studies included in the previous IPF risk genome-wide association studies. We used the distribution of effect sizes across the five studies to assess the replicability of the results and identified five robust novel genetic association signals implicating mTOR (mammalian target of rapamycin) signalling, telomere maintenance and spindle assembly genes in IPF risk.

Analysis of protein-altering variants in telomerase genes and their association with MUC5B common variant status in patients with idiopathic pulmonary fibrosis: a candidate gene sequencing study.

BACKGROUND: Idiopathic pulmonary fibrosis (IPF) risk has a strong genetic component. Studies have implicated variations at several loci, including TERT, surfactant genes, and a single nucleotide polymorphism at chr11p15 (rs35705950) in the intergenic region between TOLLIP and MUC5B. Patients with IPF who have risk alleles at rs35705950 have longer survival from the time of IPF diagnosis than do patients homozygous for the non-risk allele, whereas patients with shorter telomeres have shorter survival times. We aimed to assess whether rare protein-altering variants in genes regulating telomere length are enriched in patients with IPF homozygous for the non-risk alleles at rs35705950. METHODS: Between Nov 1, 2014, and Nov 1, 2016, we assessed blood samples from patients aged 40 years or older and of European ancestry with sporadic IPF from three international phase 3 clinical trials (INSPIRE, CAPACITY, ASCEND), one phase 2 study (RIFF), and US-based observational studies (Vanderbilt Clinical Interstitial Lung Disease Registry and the UCSF Interstitial Lung Disease Clinic registry cohorts) at the Broad Institute (Cambridge, MA, USA) and Human Longevity (San Diego, CA, USA). We also assessed blood samples from non-IPF controls in several clinical trials. We did whole-genome sequencing to assess telomere length and identify rare protein-altering variants, stratified by rs35705950 genotype. We also assessed rare functional variation in TERT exons and compared telomere length and disease progression across genotypes. FINDINGS: We assessed samples from 1510 patients with IPF and 1874 non-IPF controls. 30 (3%) of 1046 patients with an rs35705950 risk allele had a rare protein-altering variant in TERT compared with 34 (7%) of 464 non-risk allele carriers (odds ratio 0·40 [95% CI 0·24-0·66], p=0·00039). Subsequent analyses identified enrichment of rare protein-altering variants in PARN and RTEL1, and rare variation in TERC in patients with IPF compared with controls. We expanded our study population to provide a more accurate estimation of rare variant frequency in these four loci, and to calculate telomere length. The proportion of patients with at least one rare variant in TERT, PARN, TERC, or RTEL1 was higher in patients with IPF than in controls (149 [9%] of 1739 patients vs 205 [2%] of 8645 controls, p=2·44 × 10-8). Patients with IPF who had a variant in any of the four identified telomerase component genes had telomeres that were 3·69-16·10% shorter than patients without a variant in any of the four genes and had an earlier mean age of disease onset than patients without one or more variants (65·1 years [SD 7·8] vs 67·1 years [7·9], p=0·004). In the placebo arms of clinical trials, shorter telomeres were significantly associated with faster disease progression (1·7% predicted forced vital capacity per kb per year, p=0·002). Pirfenidone had treatment benefit regardless of telomere length (p=4·24 × 10-8 for telomere length lower than the median, p=0·0044 for telomere length greater than the median). INTERPRETATION: Rare protein-altering variants in TERT, PARN, TERC, and RTEL1 are enriched in patients with IPF compared with controls, and, in the case of TERT, particularly in individuals without a risk allele at the rs35705950 locus. This suggests that multiple genetic factors contribute to sporadic IPF, which might implicate distinct mechanisms of pathogenesis and disease progression. FUNDING: Genentech, National Institutes of Health, Francis Family Foundation, Pulmonary Fibrosis Foundation, Nina Ireland Program for Lung Health, US Department of Veterans Affairs.

Hypothesis-independent pathway analysis implicates GABA and acetyl-CoA metabolism in primary open-angle glaucoma and normal-pressure glaucoma.

Primary open-angle glaucoma (POAG) is a leading cause of blindness worldwide. Using genome-wide association single-nucleotide polymorphism data from the Glaucoma Genes and Environment study and National Eye Institute Glaucoma Human Genetics Collaboration comprising 3,108 cases and 3,430 controls, we assessed biologic pathways as annotated in the KEGG database for association with risk of POAG. After correction for genic overlap among pathways, we found 4 pathways, butanoate metabolism (hsa00650), hematopoietic cell lineage (hsa04640), lysine degradation (hsa00310) and basal transcription factors (hsa03022) related to POAG with permuted p < 0.001. In addition, the human leukocyte antigen (HLA) gene family was significantly associated with POAG (p < 0.001). In the POAG subset with normal-pressure glaucoma (NPG), the butanoate metabolism pathway was also significantly associated (p < 0.001) as well as the MAPK and Hedgehog signaling pathways (hsa04010 and hsa04340), glycosaminoglycan biosynthesis-heparan sulfate pathway (hsa00534) and the phenylalanine, tyrosine and tryptophan biosynthesis pathway (hsa0400). The butanoate metabolism pathway overall, and specifically the aspects of the pathway that contribute to GABA and acetyl-CoA metabolism, was the only pathway significantly associated with both POAG and NPG. Collectively these results implicate GABA and acetyl-CoA metabolism in glaucoma pathogenesis, and suggest new potential therapeutic targets.

Rare coding variants of the adenosine A3 receptor are increased in autism: on the trail of the serotonin transporter regulome.

BACKGROUND: Rare genetic variation is an important class of autism spectrum disorder (ASD) risk factors and can implicate biological networks for investigation. Altered serotonin (5-HT) signaling has been implicated in ASD, and we and others have discovered multiple, rare, ASD-associated variants in the 5-HT transporter (SERT) gene leading to elevated 5-HT re-uptake and perturbed regulation. We hypothesized that loci encoding SERT regulators harbor variants that impact SERT function and/or regulation and therefore could contribute to ASD risk. The adenosine A3 receptor (A3AR) regulates SERT via protein kinase G (PKG) and other signaling pathways leading to enhanced SERT surface expression and catalytic activity. METHODS: To test our hypothesis, we asked whether rare variants in the A3AR gene (ADORA3) were increased in ASD cases vs. controls. Discovery sequencing in a case-control sample and subsequent analysis of comparison exome sequence data were conducted. We evaluated the functional impact of two variants from the discovery sample on A3AR signaling and SERT activity. RESULTS: Sequencing discovery showed an increase of rare coding variants in cases vs. controls (P=0.013). While comparison exome sequence data did not show a significant enrichment (P=0.071), combined analysis strengthened evidence for association (P=0.0025). Two variants discovered in ASD cases (Leu90Val and Val171Ile) lie in or near the ligand-binding pocket, and Leu90Val was enriched individually in cases (P=0.040). In vitro analysis of cells expressing Val90-A3AR revealed elevated basal cGMP levels compared with the wildtype receptor. Additionally, a specific A3AR agonist increased cGMP levels across the full time course studied in Val90-A3AR cells, compared to wildtype receptor. In Val90-A3AR/SERT co-transfections, agonist stimulation elevated SERT activity over the wildtype receptor with delayed 5-HT uptake activity recovery. In contrast, Ile171-A3AR was unable to support agonist stimulation of SERT. Although both Val90 and Ile171 were present in greater numbers in these ASD cases, segregation analysis in families showed incomplete penetrance, consistent with other rare ASD risk alleles. CONCLUSIONS: Our results validate the hypothesis that the SERT regulatory network harbors rare, functional variants that impact SERT activity and regulation in ASD, and encourages further investigation of this network for other variation that may impact ASD risk.

Estrogen pathway polymorphisms in relation to primary open angle glaucoma: an analysis accounting for gender from the United States.

PURPOSE: Circulating estrogen levels are relevant in glaucoma phenotypic traits. We assessed the association between an estrogen metabolism single nucleotide polymorphism (SNP) panel in relation to primary open angle glaucoma (POAG), accounting for gender. METHODS: We included 3,108 POAG cases and 3,430 controls of both genders from the Glaucoma Genes and Environment (GLAUGEN) study and the National Eye Institute Glaucoma Human Genetics Collaboration (NEIGHBOR) consortium genotyped on the Illumina 660W-Quad platform. We assessed the relation between the SNP panels representative of estrogen metabolism and POAG using pathway- and gene-based approaches with the Pathway Analysis by Randomization Incorporating Structure (PARIS) software. PARIS executes a permutation algorithm to assess statistical significance relative to the pathways and genes of comparable genetic architecture. These analyses were performed using the meta-analyzed results from the GLAUGEN and NEIGHBOR data sets. We evaluated POAG overall as well as two subtypes of POAG defined as intraocular pressure (IOP) ≥22 mmHg (high-pressure glaucoma [HPG]) or IOP <22 mmHg (normal pressure glaucoma [NPG]) at diagnosis. We conducted these analyses for each gender separately and then jointly in men and women. RESULTS: Among women, the estrogen SNP pathway was associated with POAG overall (permuted p=0.006) and HPG (permuted p<0.001) but not NPG (permuted p=0.09). Interestingly, there was no relation between the estrogen SNP pathway and POAG when men were considered alone (permuted p>0.99). Among women, gene-based analyses revealed that the catechol-O-methyltransferase gene showed strong associations with HTG (permuted gene p≤0.001) and NPG (permuted gene p=0.01). CONCLUSIONS: The estrogen SNP pathway was associated with POAG among women.

Seven new loci associated with age-related macular degeneration.

Age-related macular degeneration (AMD) is a common cause of blindness in older individuals. To accelerate the understanding of AMD biology and help design new therapies, we executed a collaborative genome-wide association study, including >17,100 advanced AMD cases and >60,000 controls of European and Asian ancestry. We identified 19 loci associated at P < 5 × 10(-8). These loci show enrichment for genes involved in the regulation of complement activity, lipid metabolism, extracellular matrix remodeling and angiogenesis. Our results include seven loci with associations reaching P < 5 × 10(-8) for the first time, near the genes COL8A1-FILIP1L, IER3-DDR1, SLC16A8, TGFBR1, RAD51B, ADAMTS9 and B3GALTL. A genetic risk score combining SNP genotypes from all loci showed similar ability to distinguish cases and controls in all samples examined. Our findings provide new directions for biological, genetic and therapeutic studies of AMD.

Genetic variants and prostate cancer risk: candidate replication and exploration of viral restriction genes.

The genetic variants underlying the strong heritable component of prostate cancer remain largely unknown. Genome-wide association studies of prostate cancer have yielded several variants that have significantly replicated across studies, predominantly in cases unselected for family history of prostate cancer. Additional candidate gene variants have also been proposed, many evaluated within familial prostate cancer study populations. Such variants hold great potential value for risk stratification, particularly for early-onset or aggressive prostate cancer, given the comorbidities associated with current therapies. Here, we investigate a Caucasian study population of 523 independent familial prostate cancer cases and 523 age-matched controls without a personal or family history of prostate cancer. We replicate identified associations at genome-wide association study loci 8q24, 11q13, and 2p15 (P = 2.9 x 10(-4) to P = 4.7 x 10(-5)), showing study population power. We also find evidence to support reported associations at candidate genes RNASEL, EZH2, and NKX3-1 (P = 0.031 to P = 0.0085). We further explore a set of candidate genes related to RNASEL and to its role in retroviral restriction, identifying nominal associations at XPR1 and RBM9. The effects at 8q24 seem more pronounced for those diagnosed at an early age, whereas at 2p15 and RNASEL the effects were more pronounced at a later age. However, these trends did not reach statistical significance. The effects at 2p15 were statistically significantly more pronounced for those diagnosed with aggressive disease.

Heritable variation of ERBB2 and breast cancer risk.

Amplification of the epithelial growth factor receptor gene ERBB2 (HER2, NEU) in breast cancer is associated with a poor clinical prognosis. In mammary gland development, this receptor plays a role in ductal and lobuloalveolar differentiation. We conducted a systematic investigation of the role of genetic variation of the ERBB2 gene in breast cancer risk in a study of 842 histologically confirmed invasive breast cancer cases and 1,108 controls from the Shanghai Breast Cancer Study. We observed that the ERBB2 gene resides within a locus of high linkage disequilibrium, composed of three major ancestral haplotypes in the study population. These haplotypes are marked by simple tandem repeat and single nucleotide polymorphisms, including the missense variants I655V and P1170A. We observed a risk-modifying effect of a highly polymorphic simple tandem repeat within an evolutionarily conserved region, 4.4 kb upstream from the ERBB2 transcription start site. Under a dominant genetic model, the age-adjusted odds ratio was 1.74 (95% confidence interval, 1.27-2.37). Its association with breast cancer, and with breast cancer stratified by histology, by histologic grade, and by stage, remained significant after correction for multiple comparisons. In contrast, we observed no association of ERBB2 single nucleotide polymorphism haplotypes with breast cancer predisposition.

A haplotype at chromosome Xq27.2 confers susceptibility to prostate cancer.

We conducted an association study to identify risk variants for familial prostate cancer within the HPCX locus at Xq27 among Americans of Northern European descent. We investigated a total of 507 familial prostate cancer probands and 507 age-matched controls without a personal or family history of prostate cancer. The study population was subdivided into a set of training subjects to explore genetic variation of the locus potentially impacting risk of prostate cancer, and an independent set of test subjects to confirm the association and to assign significance, addressing multiple comparisons. We identified a 22.9 kb haplotype nominally associated with prostate cancer among training subjects (292 cases, 292 controls; chi(2) = 5.08, P = 0.020), that was confirmed among test subjects (215 cases, 215 controls; chi(2) = 3.73, P = 0.040). The haplotype predisposed to prostate cancer with an odds ratio of 3.41 (95% CI 1.04-11.17, P = 0.034) among test subjects. The haplotype extending from rs5907859 to rs1493189 is concordant with a prior study of the region within the Finnish founder population, and warrants further independent investigation.

Common sequence variants on 2p15 and Xp11.22 confer susceptibility to prostate cancer.

We conducted a genome-wide SNP association study on prostate cancer on over 23,000 Icelanders, followed by a replication study including over 15,500 individuals from Europe and the United States. Two newly identified variants were shown to be associated with prostate cancer: rs5945572 on Xp11.22 and rs721048 on 2p15 (odds ratios (OR) = 1.23 and 1.15; P = 3.9 x 10(-13) and 7.7 x 10(-9), respectively). The 2p15 variant shows a significantly stronger association with more aggressive, rather than less aggressive, forms of the disease.

Haplotype analysis of CYP11A1 identifies promoter variants associated with breast cancer risk.

The CYP11A1 gene encodes the cholesterol side chain cleavage enzyme that catalyzes the initial and rate-limiting step of steroidogenesis. A large number of epidemiologic studies have implicated the duration and degree of endogenous estrogen exposure in the development of breast cancer in women. Here, we conduct a systematic investigation of the role of genetic variation of the CYP11A1 gene in breast cancer risk in a study of 1193 breast cancer cases and 1310 matched controls from the Shanghai Breast Cancer Study. We characterize the genetic architecture of the CYP11A1 gene in a Chinese study population. We then genotype tagging polymorphisms to capture common variation at the locus for tests of association. Variants designating a haplotype encompassing the gene promoter are significantly associated with both increased expression (P = 1.6e-6) and increased breast cancer risk: heterozygote age-adjusted odds ratio (OR), 1.51 [95% confidence interval (95% CI), 1.19-1.91]; homozygote age-adjusted OR, 2.94 (95% CI, 1.22-7.12), test for trend, P = 5.0e-5. Among genes controlling endogenous estrogen metabolism, CYP11A1 harbors common variants that may influence expression to significantly modify risk of breast cancer.