Apolipoprotein L1 (APOL1) renal risk variant-mediated podocyte cytotoxicity depends on African haplotype and surface expression

Homozygous Apolipoprotein L1 (APOL1) variants G1 and G2 cause APOL1-mediated kidney disease, purportedly acting as surface cation channels in podocytes. APOL1-G0 exhibits various single nucleotide polymorphisms, most commonly haplotype E150K, M228I and R255K ("KIK"; the Reference Sequence is "EMR"), whereas variants G1 and G2 are mostly found in a single "African" haplotype background ("EIK"). Several labs reported cytotoxicity with risk variants G1 and G2 in KIK or EIK background haplotypes, but used HEK-293 cells and did not verify equal surface expression. To see if haplotype matters in a more relevant cell type, we induced APOL1-G0, G1 and G2 EIK, KIK and EMR at comparable surface levels in immortalized podocytes. G1 and G2 risk variants (but not G0) caused dose-dependent podocyte death within 48h only in their native African EIK haplotype and correlated with K+ conductance (thallium FLIPR). We ruled out differences in localization and trafficking, except for possibly greater surface clustering of cytotoxic haplotypes. APOL1 surface expression was required, since Brefeldin A rescued cytotoxicity; and cytoplasmic isoforms vB3 and vC were not cytotoxic. Thus, APOL1-EIK risk variants kill podocytes in a dose and haplotype-dependent manner (as in HEK-293 cells), whereas unlike in HEK-293 cells the KIK risk variants did not.

Genome-wide SNP-sex interaction analysis of susceptibility to idiopathic pulmonary fibrosis

Background

Idiopathic pulmonary fibrosis (IPF) is a chronic lung condition that is more prevalent in males than females. The reasons for this are not fully understood, with differing environmental exposures due to historically sex-biased occupations, or diagnostic bias, being possible explanations. To date, over 20 independent genetic variants have been identified to be associated with IPF susceptibility, but these have been discovered when combining males and females. Our aim was to test for the presence of sex-specific associations with IPF susceptibility and assess whether there is a need to consider sex-specific effects when evaluating genetic risk in clinical prediction models for IPF.

Methods

We performed genome-wide single nucleotide polymorphism (SNP)-by-sex interaction studies of IPF risk in six independent IPF case-control studies and combined them using inverse-variance weighted fixed effect meta-analysis. In total, 4,561 cases (1,280 females and 2,281 males) and 23,500 controls (8,360 females and 14,528 males) of European genetic ancestry were analysed. We used polygenic risk scores (PRS) to assess differences in genetic risk prediction between males and females.

Findings

Three independent genetic association signals were identified. All showed a consistent direction of effect across all individual IPF studies and an opposite direction of effect in IPF susceptibility between females and males. None had been previously identified in IPF susceptibility genome-wide association studies (GWAS). The predictive accuracy of the PRSs were similar between males and females, regardless of whether using combined or sex-specific GWAS results.

Interpretation

We prioritised three genetic variants whose effect on IPF risk may be modified by sex, however these require further study. We found no evidence that the predictive accuracy of common SNP-based PRSs varies significantly between males and females.

Association study of human leukocyte antigen variants and idiopathic pulmonary fibrosis

Introduction

Idiopathic pulmonary fibrosis (IPF) is a chronic interstitial pneumonia marked by progressive lung fibrosis and a poor prognosis. Recent studies have highlighted the potential role of infection in the pathogenesis of IPF, and a prior association of the HLA-DQB1 gene with idiopathic fibrotic interstitial pneumonia (including IPF) has been reported. Owing to the important role that the human leukocyte antigen (HLA) region plays in the immune response, here we evaluated if HLA genetic variation was associated specifically with IPF risk.

Methods

We performed a meta-analysis of associations of the HLA region with IPF risk in individuals of European ancestry from seven independent case-control studies of IPF (comprising 5159 cases and 27 459 controls, including a prior study of fibrotic interstitial pneumonia). Single nucleotide polymorphisms, classical HLA alleles and amino acids were analysed and signals meeting a region-wide association threshold of p<4.5×10-4 and a posterior probability of replication >90% were considered significant. We sought to replicate the previously reported HLA-DQB1 association in the subset of studies independent of the original report.

Results

The meta-analysis of all seven studies identified four significant independent single nucleotide polymorphisms associated with IPF risk. However, none met the posterior probability for replication criterion. The HLA-DQB1 association was not replicated in the independent IPF studies.

Conclusion

Variation in the HLA region was not consistently associated with risk in studies of IPF. However, this does not preclude the possibility that other genomic regions linked to the immune response may be involved in the aetiology of IPF.

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.

Association study of human leukocyte antigen (HLA) variants and idiopathic pulmonary fibrosis

Introduction

Idiopathic pulmonary fibrosis (IPF) is a chronic interstitial pneumonia marked by progressive lung fibrosis and a poor prognosis. Recent studies have highlighted the potential role of infection in the pathogenesis of IPF and a prior association of the HLA-DQB1 gene with idiopathic fibrotic interstitial pneumonia (including IPF) has been reported. Due to the important role that the Human Leukocyte Antigen (HLA) region plays in the immune response, here we evaluated if HLA genetic variation was associated specifically with IPF risk.

Methods

We performed a meta-analysis of associations of the HLA region with IPF risk in individuals of European ancestry from seven independent case-control studies of IPF (comprising a total of 5,159 cases and 27,459 controls, including the prior study of fibrotic interstitial pneumonia). Single nucleotide polymorphisms, classical HLA alleles and amino acids were analysed and signals meeting a region-wide association threshold p<4.5×10-4 and a posterior probability of replication >90% were considered significant. We sought to replicate the previously reported HLA-DQB1 association in the subset of studies independent of the original report.

Results

The meta-analysis of all seven studies identified four significant independent single nucleotide polymorphisms associated with IPF risk. However, none met the posterior probability for replication criterion. The HLA-DQB1 association was not replicated in the independent IPF studies.

Conclusion

Variation in the HLA region was not consistently associated with risk in studies of IPF. However, this does not preclude the possibility that other genomic regions linked to the immune response may be involved in the aetiology of IPF.

A systems biology approach uncovers novel disease mechanisms in age-related macular degeneration

Age-related macular degeneration (AMD) is a leading cause of blindness, affecting 200 million people worldwide. To identify genes that could be targeted for treatment, we created a molecular atlas at different stages of AMD. Our resource is comprised of RNA sequencing (RNA-seq) and DNA methylation microarrays from bulk macular retinal pigment epithelium (RPE)/choroid of clinically phenotyped normal and AMD donor eyes (n = 85), single-nucleus RNA-seq (164,399 cells), and single-nucleus assay for transposase-accessible chromatin (ATAC)-seq (125,822 cells) from the retina, RPE, and choroid of 6 AMD and 7 control donors. We identified 23 genome-wide significant loci differentially methylated in AMD, over 1,000 differentially expressed genes across different disease stages, and an AMD Müller state distinct from normal or gliosis. Chromatin accessibility peaks in genome-wide association study (GWAS) loci revealed putative causal genes for AMD, including HTRA1 and C6orf223. Our systems biology approach uncovered molecular mechanisms underlying AMD, including regulators of WNT signaling, FRZB and TLE2, as mechanistic players in disease.

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.

A genome-wide association study of chronic spontaneous urticaria risk and heterogeneity

BACKGROUND

Chronic spontaneous urticaria (CSU) is a dermatologic condition characterized by spontaneous, pruritic hives and/or angioedema that persists for 6 weeks or longer with no identifiable trigger. Antihistamines and second-line therapies such as omalizumab are effective for some CSU patients, but others remain symptomatic, with significant impact on quality of life. This variable response to treatment and autoantibody levels across patients highlight clinically heterogeneous subgroups.

OBJECTIVE

We aimed to highlight pathways involved in CSU by investigating the genetics of CSU risk and subgroups.

METHODS

We performed a genome-wide association study (GWAS) of 679 CSU patients and 4446 controls and a GWAS of chronic urticaria (CU)-index, which measures IgG autoantibodies levels, by comparing 447 CU index-low to 183 CU index-high patients. We also tested whether polygenic scores for autoimmune-related disorders were associated with CSU risk and CU index.

RESULTS

We identified 2 loci significantly associated with disease risk. The strongest association mapped to position 56 of HLA-DQA1 (P = 1.69 × 10-9), where the arginine residue was associated with increased risk (odds ratio = 1.64). The second association signal colocalized with expression-quantitative trait loci for ITPKB in whole blood (Pcolocalization = .997). The arginine residue at position 56 of HLA-DQA1 was also associated with increased risk of CU index-high (P = 6.15 × 10-5, odds ratio = 1.86), while the ITKPB association was not (P = .64). Polygenic scores for 3 autoimmune-related disorders (hypothyroidism, type 1 diabetes, and vitiligo) were associated with CSU risk and CU index (P < 2.34 × 10-3, odds ratio > 1.72).

CONCLUSION

A GWAS of CSU identified 2 genome-wide significant loci, highlighting the shared genetics between CU index and autoimmune disorders.

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.

Dual antibody inhibition of KLK5 and KLK7 for Netherton syndrome and atopic dermatitis

The epidermis is a barrier that prevents water loss while keeping harmful substances from penetrating the host. The impermeable cornified layer of the stratum corneum is maintained by balancing continuous turnover driven by epidermal basal cell proliferation, suprabasal cell differentiation, and corneal shedding. The epidermal desquamation process is tightly regulated by balance of the activities of serine proteases of the Kallikrein-related peptidases (KLK) family and their cognate inhibitor lymphoepithelial Kazal type-related inhibitor (LEKTI), which is encoded by the serine peptidase inhibitor Kazal type 5 gene. Imbalance of proteolytic activity caused by a deficiency of LEKTI leads to excessive desquamation due to increased activities of KLK5, KLK7, and KLK14 and results in Netherton syndrome (NS), a debilitating condition with an unmet clinical need. Increased activity of KLKs may also be pathological in other dermatoses such as atopic dermatitis (AD). Here, we describe the discovery of inhibitory antibodies against murine KLK5 and KLK7 that could compensate for the deficiency of LEKTI in NS. These antibodies are protective in mouse models of NS and AD and, when combined, promote improved skin barrier integrity and reduced inflammation. To translate these findings, we engineered a humanized bispecific antibody capable of potent inhibition of human KLK5 and KLK7. A crystal structure of KLK5 bound to the inhibitory Fab revealed that the antibody binds distal to its active site and uses a relatively unappreciated allosteric inhibition mechanism. Treatment with the bispecific anti-KLK5/7 antibody represents a promising therapy for clinical development in NS and other inflammatory dermatoses.

Leveraging global multi-ancestry meta-analysis in the study of idiopathic pulmonary fibrosis genetics

The research of rare and devastating orphan diseases, such as idiopathic pulmonary fibrosis (IPF) has been limited by the rarity of the disease itself. The prognosis is poor-the prevalence of IPF is only approximately four times the incidence, limiting the recruitment of patients to trials and studies of the underlying biology. Global biobanking efforts can dramatically alter the future of IPF research. We describe a large-scale meta-analysis of IPF, with 8,492 patients and 1,355,819 population controls from 13 biobanks around the globe. Finally, we combine this meta-analysis with the largest available meta-analysis of IPF, reaching 11,160 patients and 1,364,410 population controls. We identify seven novel genome-wide significant loci, only one of which would have been identified if the analysis had been limited to European ancestry individuals. We observe notable pleiotropy across IPF susceptibility and severe COVID-19 infection and note an unexplained sex-heterogeneity effect at the strongest IPF locus MUC5B.

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.

Genome-wide association study identifies kallikrein 5 in type 2 inflammation-low asthma

BACKGROUND

Clinical studies of type 2 (T2) cytokine-related neutralizing antibodies in asthma have identified a substantial subset of patients with low levels of T2 inflammation who do not benefit from T2 cytokine neutralizing antibody treatment. Non type 2 mechanisms are poorly understood in asthma but represent the redefined unmet medical need.

OBJECTIVE

To gain understanding of the genetic contribution to T2-low asthma.

METHODS

We utilized an unbiased genome-wide association study (GWAS) of moderate-severe asthma patients stratified by T2 serum biomarker periostin. We also performed additional expression and biological analysis for the top genetic hit.

RESULTS

This analysis identified a novel protective SNP at chr19q13.41 which is selectively associated with T2-low asthma and establishes KLK5 as the causal gene mediating this association. Heterozygous carriers of the SNP have reduced KLK5 expression. KLK5 is secreted by human bronchial epithelial cells and elevated in asthma bronchial alveolar lavage. T2 cytokines IL-4 and IL-13 downregulate KLK5 in human bronchial epithelial cells. KLK5, dependent on its catalytic function, induces epithelial chemokine/cytokine expression. Lastly, overexpression of KLK5 in airway, or lack of an endogenous KLK5 inhibitor, SPINK5, leads to spontaneous airway neutrophilic inflammation.

CONCLUSION

Our data identifies KLK5 as the causal gene at a novel locus at chr19q13.41 associated with T2-low asthma.

A whole genome sequencing study of moderate to severe asthma identifies a lung function locus associated with asthma risk

Genome-wide association studies (GWAS) have identified many common variant loci associated with asthma susceptibility, but few studies investigate the genetics underlying moderate-to-severe asthma risk. Here, we present a whole-genome sequencing study comparing 3181 moderate-to-severe asthma patients to 3590 non-asthma controls. We demonstrate that asthma risk is genetically correlated with lung function measures and that this component of asthma risk is orthogonal to the eosinophil genetics that also contribute to disease susceptibility. We find that polygenic scores for reduced lung function are associated with younger asthma age of onset. Genome-wide, seven previously reported common asthma variant loci and one previously reported lung function locus, near THSD4, reach significance. We replicate association of the lung function locus in a recently published GWAS of moderate-to-severe asthma patients. We additionally replicate the association of a previously reported rare (minor allele frequency < 1%) coding variant in IL33 and show significant enrichment of rare variant burden in genes from common variant allergic disease loci. Our findings highlight the contribution of lung function genetics to moderate-to-severe asthma risk, and provide initial rare variant support for associations with moderate-to-severe asthma risk at several candidate genes from common variant loci.

An Allosteric Anti-tryptase Antibody for the Treatment of Mast Cell-Mediated Severe Asthma

Severe asthma patients with low type 2 inflammation derive less clinical benefit from therapies targeting type 2 cytokines and represent an unmet need. We show that mast cell tryptase is elevated in severe asthma patients independent of type 2 biomarker status. Active β-tryptase allele count correlates with blood tryptase levels, and asthma patients carrying more active alleles benefit less from anti-IgE treatment. We generated a noncompetitive inhibitory antibody against human β-tryptase, which dissociates active tetramers into inactive monomers. A 2.15 Å crystal structure of a β-tryptase/antibody complex coupled with biochemical studies reveal the molecular basis for allosteric destabilization of small and large interfaces required for tetramerization. This anti-tryptase antibody potently blocks tryptase enzymatic activity in a humanized mouse model, reducing IgE-mediated systemic anaphylaxis, and inhibits airway tryptase in Ascaris-sensitized cynomolgus monkeys with favorable pharmacokinetics. These data provide a foundation for developing anti-tryptase as a clinical therapy for severe asthma.

Genome-Wide Association Study of Susceptibility to Idiopathic Pulmonary Fibrosis

Rationale: Idiopathic pulmonary fibrosis (IPF) is a complex lung disease characterized by scarring of the lung that is believed to result from an atypical response to injury of the epithelium. Genome-wide association studies have reported signals of association implicating multiple pathways including host defense, telomere maintenance, signaling, and cell-cell adhesion.Objectives: To improve our understanding of factors that increase IPF susceptibility by identifying previously unreported genetic associations.Methods: We conducted genome-wide analyses across three independent studies and meta-analyzed these results to generate the largest genome-wide association study of IPF to date (2,668 IPF cases and 8,591 controls). We performed replication in two independent studies (1,456 IPF cases and 11,874 controls) and functional analyses (including statistical fine-mapping, investigations into gene expression, and testing for enrichment of IPF susceptibility signals in regulatory regions) to determine putatively causal genes. Polygenic risk scores were used to assess the collective effect of variants not reported as associated with IPF.Measurements and Main Results: We identified and replicated three new genome-wide significant (P < 5 × 10-8) signals of association with IPF susceptibility (associated with altered gene expression of KIF15, MAD1L1, and DEPTOR) and confirmed associations at 11 previously reported loci. Polygenic risk score analyses showed that the combined effect of many thousands of as yet unreported IPF susceptibility variants contribute to IPF susceptibility.Conclusions: The observation that decreased DEPTOR expression associates with increased susceptibility to IPF supports recent studies demonstrating the importance of mTOR signaling in lung fibrosis. New signals of association implicating KIF15 and MAD1L1 suggest a possible role of mitotic spindle-assembly genes in IPF susceptibility.

Rare Protein-Altering Telomere-related Gene Variants in Patients with Chronic Hypersensitivity Pneumonitis

Rationale: Rare genetic variants in telomere-related genes have been identified in familial, idiopathic, and rheumatoid arthritis-associated pulmonary fibrosis. Short peripheral blood leukocyte (PBL) telomere length predicts poor outcomes in chronic hypersensitivity pneumonitis (CHP).Objectives: Determine the prevalence and clinical relevance of rare protein-altering variants in telomere-related genes in patients with CHP.Methods: Next-generation sequences from two CHP cohorts were analyzed to identify variants in TERT (telomerase reverse transcriptase), TERC (telomerase RNA component), DKC1 (dyskerin pseudouridine synthase 1), RTEL1 (regulator of telomere elongation helicase 1), PARN (poly[A]-specific RNase), and TINF2 (TERF1-interacting nuclear factor 2). To qualify, variants were required to have a minor allele frequency less than 0.005 and be predicted to be damaging to protein function. Variant status (binary variable) was used in statistical association tests, including Cox proportional hazard models for transplant-free survival. PBL telomere length was measured using quantitative PCR.Measurements and Main Results: Qualifying variants were identified in 16 of 144 patients (11.1%; 95% confidence interval [CI], 6.5-17.4) in the discovery cohort and 17 of 209 patients (8.1%; 95% CI, 4.8-12.7) in the replication cohort. Age- and ancestry-adjusted PBL telomere length was significantly shorter in the presence of a variant in both cohorts (discovery: -561 bp; 95% CI, -933 to -190; P = 0.003; replication: -612 bp; 95% CI, -870 to -354; P = 5.30 × 10-6). Variant status was significantly associated with transplant-free survival in both cohorts (discovery: age-, sex-, and ancestry-adjusted hazard ratio, 3.73; 95% CI, 1.92-7.28; P = 0.0001; replication: hazard ratio, 2.72; 95% CI, 1.26-5.88; P = 0.011).Conclusions: A substantial proportion of patients diagnosed with CHP have rare, protein-altering variants in telomere-related genes, which are associated with short peripheral blood telomere length and significantly reduced transplant-free survival.

Multiethnic Genome-Wide Association Study of Diabetic Retinopathy Using Liability Threshold Modeling of Duration of Diabetes and Glycemic Control

To identify genetic variants associated with diabetic retinopathy (DR), we performed a large multiethnic genome-wide association study. Discovery included eight European cohorts (n = 3,246) and seven African American cohorts (n = 2,611). We meta-analyzed across cohorts using inverse-variance weighting, with and without liability threshold modeling of glycemic control and duration of diabetes. Variants with a P value <1 × 10-5 were investigated in replication cohorts that included 18,545 European, 16,453 Asian, and 2,710 Hispanic subjects. After correction for multiple testing, the C allele of rs142293996 in an intron of nuclear VCP-like (NVL) was associated with DR in European discovery cohorts (P = 2.1 × 10-9), but did not reach genome-wide significance after meta-analysis with replication cohorts. We applied the Disease Association Protein-Protein Link Evaluator (DAPPLE) to our discovery results to test for evidence of risk being spread across underlying molecular pathways. One protein-protein interaction network built from genes in regions associated with proliferative DR was found to have significant connectivity (P = 0.0009) and corroborated with gene set enrichment analyses. These findings suggest that genetic variation in NVL, as well as variation within a protein-protein interaction network that includes genes implicated in inflammation, may influence risk for DR.

Strategies for Pathway Analysis Using GWAS and WGS Data

Single-allele study designs, commonly used in genome-wide association studies (GWAS) as well as the more recently developed whole genome sequencing (WGS) studies, are a standard approach for investigating the relationship of common variation within the human genome to a given phenotype of interest. However, single-allele association results published for many GWAS studies represent only the tip of the iceberg for the information that can be extracted from these datasets. The primary analysis strategy for GWAS entails association analysis in which only the single nucleotide polymorphisms (SNPs) with the strongest p-values are declared statistically significant due to issues arising from multiple testing and type I errors. Factors such as locus heterogeneity, epistasis, and multiple genes conferring small effects contribute to the complexity of the genetic models underlying phenotype expression. Thus, many biologically meaningful associations having lower effect sizes at individual genes are overlooked, making it difficult to separate true associations from a sea of false-positive associations. Organizing these individual SNPs into biologically meaningful groups to look at the overall effects of minor perturbations to genes and pathways is desirable. This pathway-based approach provides researchers with insight into the functional foundations of the phenotype being studied and allows testing of various genetic scenarios. © 2018 by John Wiley & Sons, Inc.

Regulation of T-cell receptor downmodulation by DENND1B in TH2 cells and allergic disease

The DENN-domain is an evolutionary conserved protein module found in all eukaryotes and serves as an exchange factor for Rab-GTPases to regulate diverse cellular functions. Variants in DENND1B are associated with development of childhood asthma and other immune disorders. To understand how DENND1B may contribute to diseases, Dennd1b−/− mice were generated and exhibit hyper-allergic responses following antigen challenge. Dennd1b−/− TH2, but not other TH cells, exhibit delayed receptor-induced TCR downmodulation, enhanced TCR signaling, and increased production of effector cytokines. As DENND1B interacts with AP-2 and Rab35, TH2 cells deficient in AP-2 or Rab35 also exhibit enhanced TCR-mediated effector functions. Moreover, human TH2 cells carrying asthma-associated DENND1B variants express less DENND1B and phenocopy Dennd1b−/− TH2 cells. These results provide a molecular basis for how DENND1B, a previously unrecognized regulator of TCR downmodulation in TH2 cells, contributes to asthma pathogenesis and how DENN-domain containing proteins may contribute to other human disorders.

Pathway analysis by randomization incorporating structure-PARIS: an update

MOTIVATION

We present an update to the pathway enrichment analysis tool 'Pathway Analysis by Randomization Incorporating Structure (PARIS)' that determines aggregated association signals generated from genome-wide association study results. Pathway-based analyses highlight biological pathways associated with phenotypes. PARIS uses a unique permutation strategy to evaluate the genomic structure of interrogated pathways, through permutation testing of genomic features, thus eliminating many of the over-testing concerns arising with other pathway analysis approaches.

RESULTS

We have updated PARIS to incorporate expanded pathway definitions through the incorporation of new expert knowledge from multiple database sources, through customized user provided pathways, and other improvements in user flexibility and functionality.

AVAILABILITY AND IMPLEMENTATION

PARIS is freely available to all users at https://ritchielab.psu.edu/software/paris-download

CONTACT

jnc43@case.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-wide association analysis identifies TXNRD2, ATXN2 and FOXC1 as susceptibility loci for primary open-angle glaucoma

Primary open angle glaucoma (POAG) is a leading cause of blindness world-wide. To identify new susceptibility loci, we meta-analyzed GWAS results from 8 independent studies from the United States (3,853 cases and 33,480 controls) and investigated the most significant SNPs in two Australian studies (1,252 cases and 2,592 controls), 3 European studies (875 cases and 4,107 controls) and a Singaporean Chinese study (1,037 cases and 2,543 controls). A meta-analysis of top SNPs identified three novel loci: rs35934224[T] within TXNRD2 (odds ratio (OR) = 0.78, P = 4.05×10⁻¹¹ encoding a mitochondrial protein required for redox homeostasis; rs7137828[T] within ATXN2 (OR = 1.17, P = 8.73×10⁻¹⁰), and rs2745572[A] upstream of FOXC1 (OR = 1.17, P = 1.76×10⁻¹⁰). Using RT-PCR and immunohistochemistry, we show TXNRD2 and ATXN2 expression in retinal ganglion cells and the optic nerve head. These results identify new pathways underlying POAG susceptibility and suggest novel targets for preventative therapies.

Regulation of T Cell Receptor Signaling by DENND1B in TH2 Cells and Allergic Disease

The DENN domain is an evolutionary conserved protein module found in all eukaryotes and serves as an exchange factor for Rab-GTPases to regulate diverse cellular functions. Variants in DENND1B are associated with development of childhood asthma and other immune disorders. To understand how DENND1B may contribute to human disease, Dennd1b(-/-) mice were generated and exhibit hyper-allergic responses following antigen challenge. Dennd1b(-/-) TH2, but not other TH cells, exhibit delayed receptor-induced T cell receptor (TCR) downmodulation, enhanced TCR signaling, and increased production of effector cytokines. As DENND1B interacts with AP-2 and Rab35, TH2 cells deficient in AP-2 or Rab35 also exhibit enhanced TCR-mediated effector functions. Moreover, human TH2 cells carrying asthma-associated DENND1B variants express less DENND1B and phenocopy Dennd1b(-/-) TH2 cells. These results provide a molecular basis for how DENND1B, a previously unrecognized regulator of TCR downmodulation in TH2 cells, contributes to asthma pathogenesis and how DENN-domain-containing proteins may contribute to other human disorders.

Age-related macular degeneration: Complement in action

The complement system plays a key role in host-defense against common pathogens but must be tightly controlled to avoid inflammation and tissue damage. Polymorphisms in genes encoding two important negative regulators of the alternative complement pathway, complement factor H (CFH) and complement factor I (CFI), are associated with the risk for Age-Related Macular Degeneration (AMD), a leading cause of vision impairment in the ageing population. In this review, we will discuss the genetic basis of AMD and the potential impact of complement de-regulation on disease pathogenesis. Finally, we will highlight recent therapeutic approaches aimed at controlling complement activation in patients with AMD.

Stratified medicine in inflammatory disorders: From theory to practice

Chronic inflammatory disorders are complex and characterized by significant heterogeneity in molecular, pathological, and clinical features. This heterogeneity poses challenges for the development of targeted molecular interventions for these disorders, as not all patients with a given clinical diagnosis have disease driven by a single dominant molecular pathway, hence not all patients will benefit equally from a given intervention. Biomarkers related to molecular manifestations of disease are increasingly being applied to enable stratified approaches to drug development. Biomarkers may be used to identify which patients are most likely to benefit from an intervention (predictive), identify patients at increased risk of disease progression (prognostic), and monitor biological responsiveness to an intervention (pharmacodynamic). Here we consider how biomarker-guided stratification of patients may increase benefit from targeted therapies for asthma, rheumatoid arthritis and inflammatory bowel diseases.

DNA copy number variants of known glaucoma genes in relation to primary open-angle glaucoma

PURPOSE

We examined the role of DNA copy number variants (CNVs) of known glaucoma genes in relation to primary open angle glaucoma (POAG).

METHODS

Our study included DNA samples from two studies (NEIGHBOR and GLAUGEN). All the samples were genotyped with the Illumina Human660W_Quad_v1 BeadChip. After removing non-blood-derived and amplified DNA samples, we applied quality control steps based on the mean Log R Ratio and the mean B allele frequency. Subsequently, data from 3057 DNA samples (1599 cases and 1458 controls) were analyzed with PennCNV software. We defined CNVs as those ≥5 kilobases (kb) in size and interrogated by ≥5 consecutive probes. We further limited our investigation to CNVs in known POAG-related genes, including CDKN2B-AS1, TMCO1, SIX1/SIX6, CAV1/CAV2, the LRP12-ZFPM2 region, GAS7, ATOH7, FNDC3B, CYP1B1, MYOC, OPTN, WDR36, SRBD1, TBK1, and GALC.

RESULTS

Genomic duplications of CDKN2B-AS1 and TMCO1 were each found in a single case. Two cases carried duplications in the GAS7 region. Genomic deletions of SIX6 and ATOH7 were each identified in one case. One case carried a TBK1 deletion and another case carried a TBK1 duplication. No controls had duplications or deletions in these six genes. A single control had a duplication in the MYOC region. Deletions of GALC were observed in five cases and two controls.

CONCLUSIONS

The CNV analysis of a large set of cases and controls revealed the presence of rare CNVs in known POAG susceptibility genes. Our data suggest that these rare CNVs may contribute to POAG pathogenesis and merit functional evaluation.

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.

Association of CAV1/CAV2 genomic variants with primary open-angle glaucoma overall and by gender and pattern of visual field loss

PURPOSE

The CAV1/CAV2 (caveolin 1 and caveolin 2) genomic region previously was associated with primary open-angle glaucoma (POAG), although replication among independent studies has been variable. The aim of this study was to assess the association between CAV1/CAV2 single nucleotide polymorphisms (SNPs) and POAG in a large case-control dataset and to explore associations by gender and pattern of visual field (VF) loss further.

DESIGN

Case-control study.

PARTICIPANTS

We analyzed 2 large POAG data sets: the Glaucoma Genes and Environment (GLAUGEN) study (976 cases, 1140 controls) and the National Eye Institute Glaucoma Human Genetics Collaboration (NEIGHBOR) consortium (2132 cases, 2290 controls).

METHODS

We studied the association between 70 SNPs located within the CAV1/CAV2 genomic region in the GLAUGEN and NEIGHBOR studies, both genotyped on the Illumina Human 660WQuadv1C BeadChip array and imputed with the Markov Chain Haplotyping algorithm using the HapMap 3 reference panel. We used logistic regression models of POAG in the overall population and separated by gender, as well as by POAG subtypes defined by type of VF defect (peripheral or paracentral). Results from GLAUGEN and NEIGHBOR were meta-analyzed, and a Bonferroni-corrected significance level of 7.7 × 10(-4) was used to account for multiple comparisons.

MAIN OUTCOME MEASURES

Overall POAG, overall POAG by gender, and POAG subtypes defined by pattern of early VF loss.

RESULTS

We found significant associations between 10 CAV1/CAV2 SNPs and POAG (top SNP, rs4236601; pooled P = 2.61 × 10(-7)). Of these, 9 were significant only in women (top SNP, rs4236601; pooled P = 1.59 × 10(-5)). Five of the 10 CAV1/CAV2 SNPs were associated with POAG with early paracentral VF (top SNP, rs17588172; pooled P = 1.07 × 10(-4)), and none of the 10 were associated with POAG with peripheral VF loss only or POAG among men.

CONCLUSIONS

CAV1/CAV2 SNPs were associated significantly with POAG overall, particularly among women. Furthermore, we found an association between CAV1/CAV2 SNPs and POAG with paracentral VF defects. These data support a role for caveolin 1, caveolin 2, or both in POAG and suggest that the caveolins particularly may affect POAG pathogenesis in women and in patients with early paracentral VF defects.

Genome-wide association study and meta-analysis of intraocular pressure

Elevated intraocular pressure (IOP) is a major risk factor for glaucoma and is influenced by genetic and environmental factors. Recent genome-wide association studies (GWAS) reported associations with IOP at TMCO1 and GAS7, and with primary open-angle glaucoma (POAG) at CDKN2B-AS1, CAV1/CAV2, and SIX1/SIX6. To identify novel genetic variants and replicate the published findings, we performed GWAS and meta-analysis of IOP in >6,000 subjects of European ancestry collected in three datasets: the NEI Glaucoma Human genetics collaBORation, GLAUcoma Genes and ENvironment study, and a subset of the Age-related Macular Degeneration-Michigan, Mayo, AREDS and Pennsylvania study. While no signal achieved genome-wide significance in individual datasets, a meta-analysis identified significant associations with IOP at TMCO1 (rs7518099-G, p = 8.0 × 10(-8)). Focused analyses of five loci previously reported for IOP and/or POAG, i.e., TMCO1, CDKN2B-AS1, GAS7, CAV1/CAV2, and SIX1/SIX6, revealed associations with IOP that were largely consistent across our three datasets, and replicated the previously reported associations in both effect size and direction. These results confirm the involvement of common variants in multiple genomic regions in regulating IOP and/or glaucoma risk.

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.

CDKN2B-AS1 genotype-glaucoma feature correlations in primary open-angle glaucoma patients from the United States

PURPOSE

To assess the association between single nucleotide polymorphisms (SNPs) of the gene region containing cyclin-dependent kinase inhibitor 2B antisense noncoding RNA (CDKN2B-AS1) and glaucoma features among primary open-angle glaucoma (POAG) patients.

DESIGN

Retrospective observational case series.

METHODS

We studied associations between 10 CDKN2B-AS1 SNPs and glaucoma features among 976 POAG cases from the Glaucoma Genes and Environment (GLAUGEN) study and 1971 cases from the National Eye Institute Glaucoma Human Genetics Collaboration (NEIGHBOR) consortium. For each patient, we chose the feature from the eye with the higher value. We created cohort-specific multivariable models for glaucoma features and then meta-analyzed the results.

RESULTS

For 9 of the 10 protective CDKN2B-AS1 SNPs with minor alleles associated with reduced disease risk (eg, the G allele at rs2157719), POAG patients carrying these minor alleles had smaller cup-to-disc ratio (0.05 units smaller per G allele at diagnosis; 95% CI: -0.08, -0.03; P = 6.23E-05) despite having higher intraocular pressure (IOP) (0.70 mm Hg higher per G allele at DNA collection; 95% CI: 0.40, 1.00; P = 5.45E-06). For the 1 adverse rs3217992 SNP with minor allele A associated with increased disease risk, POAG patients with A alleles had larger cup-to-disc ratio (0.05 units larger per A allele at diagnosis; 95% CI: 0.02, 0.07; P = 4.74E-04) despite having lower IOP (-0.57 mm Hg per A allele at DNA collection; 95% CI: -0.84, -0.29; P = 6.55E-05).

CONCLUSION

Alleles of CDKN2B-AS1 SNPs, which influence risk of developing POAG, also modulate optic nerve degeneration among POAG patients, underscoring the role of CDKN2B-AS1 in POAG.

Genome-wide association analyses identify multiple loci associated with central corneal thickness and keratoconus

Central corneal thickness (CCT) is associated with eye conditions including keratoconus and glaucoma. We performed a meta-analysis on >20,000 individuals in European and Asian populations that identified 16 new loci associated with CCT at genome-wide significance (P < 5 × 10(-8)). We further showed that 2 CCT-associated loci, FOXO1 and FNDC3B, conferred relatively large risks for keratoconus in 2 cohorts with 874 cases and 6,085 controls (rs2721051 near FOXO1 had odds ratio (OR) = 1.62, 95% confidence interval (CI) = 1.4-1.88, P = 2.7 × 10(-10), and rs4894535 in FNDC3B had OR = 1.47, 95% CI = 1.29-1.68, P = 4.9 × 10(-9)). FNDC3B was also associated with primary open-angle glaucoma (P = 5.6 × 10(-4); tested in 3 cohorts with 2,979 cases and 7,399 controls). Further analyses implicate the collagen and extracellular matrix pathways in the regulation of CCT.

Genome-wide analysis of central corneal thickness in primary open-angle glaucoma cases in the NEIGHBOR and GLAUGEN consortia

PURPOSE

To investigate the effects of central corneal thickness (CCT)-associated variants on primary open-angle glaucoma (POAG) risk using single nucleotide polymorphisms (SNP) data from the Glaucoma Genes and Environment (GLAUGEN) and National Eye Institute (NEI) Glaucoma Human Genetics Collaboration (NEIGHBOR) consortia.

METHODS

A replication analysis of previously reported CCT SNPs was performed in a CCT dataset (n = 1117) and these SNPs were then tested for association with POAG using a larger POAG dataset (n = 6470). Then a CCT genome-wide association study (GWAS) was performed. Top SNPs from this analysis were selected and tested for association with POAG. cDNA libraries from fetal and adult brain and ocular tissue samples were generated and used for candidate gene expression analysis.

RESULTS

Association with one of 20 previously published CCT SNPs was replicated: rs12447690, near the ZNF469 gene (P = 0.001; β = -5.08 μm/allele). None of these SNPs were significantly associated with POAG. In the CCT GWAS, no SNPs reached genome-wide significance. After testing 50 candidate SNPs for association with POAG, one SNP was identified, rs7481514 within the neurotrimin (NTM) gene, that was significantly associated with POAG in a low-tension subset (P = 0.00099; Odds Ratio [OR] = 1.28). Additionally, SNPs in the CNTNAP4 gene showed suggestive association with POAG (top SNP = rs1428758; P = 0.018; OR = 0.84). NTM and CNTNAP4 were shown to be expressed in ocular tissues.

CONCLUSIONS

The results suggest previously reported CCT loci are not significantly associated with POAG susceptibility. By performing a quantitative analysis of CCT and a subsequent analysis of POAG, SNPs in two cell adhesion molecules, NTM and CNTNAP4, were identified and may increase POAG susceptibility in a subset of cases.

Common variants at 9p21 and 8q22 are associated with increased susceptibility to optic nerve degeneration in glaucoma

Optic nerve degeneration caused by glaucoma is a leading cause of blindness worldwide. Patients affected by the normal-pressure form of glaucoma are more likely to harbor risk alleles for glaucoma-related optic nerve disease. We have performed a meta-analysis of two independent genome-wide association studies for primary open angle glaucoma (POAG) followed by a normal-pressure glaucoma (NPG, defined by intraocular pressure (IOP) less than 22 mmHg) subgroup analysis. The single-nucleotide polymorphisms that showed the most significant associations were tested for association with a second form of glaucoma, exfoliation-syndrome glaucoma. The overall meta-analysis of the GLAUGEN and NEIGHBOR dataset results (3,146 cases and 3,487 controls) identified significant associations between two loci and POAG: the CDKN2BAS region on 9p21 (rs2157719 [G], OR = 0.69 [95%CI 0.63-0.75], p = 1.86×10⁻¹⁸), and the SIX1/SIX6 region on chromosome 14q23 (rs10483727 [A], OR = 1.32 [95%CI 1.21-1.43], p = 3.87×10⁻¹¹). In sub-group analysis two loci were significantly associated with NPG: 9p21 containing the CDKN2BAS gene (rs2157719 [G], OR = 0.58 [95% CI 0.50-0.67], p = 1.17×10⁻¹²) and a probable regulatory region on 8q22 (rs284489 [G], OR = 0.62 [95% CI 0.53-0.72], p = 8.88×10⁻¹⁰). Both NPG loci were also nominally associated with a second type of glaucoma, exfoliation syndrome glaucoma (rs2157719 [G], OR = 0.59 [95% CI 0.41-0.87], p = 0.004 and rs284489 [G], OR = 0.76 [95% CI 0.54-1.06], p = 0.021), suggesting that these loci might contribute more generally to optic nerve degeneration in glaucoma. Because both loci influence transforming growth factor beta (TGF-beta) signaling, we performed a genomic pathway analysis that showed an association between the TGF-beta pathway and NPG (permuted p = 0.009). These results suggest that neuro-protective therapies targeting TGF-beta signaling could be effective for multiple forms of glaucoma.

A novel approach of homozygous haplotype sharing identifies candidate genes in autism spectrum disorder

Autism spectrum disorder (ASD) is a highly heritable disorder of complex and heterogeneous aetiology. It is primarily characterized by altered cognitive ability including impaired language and communication skills and fundamental deficits in social reciprocity. Despite some notable successes in neuropsychiatric genetics, overall, the high heritability of ASD (~90%) remains poorly explained by common genetic risk variants. However, recent studies suggest that rare genomic variation, in particular copy number variation, may account for a significant proportion of the genetic basis of ASD. We present a large scale analysis to identify candidate genes which may contain low-frequency recessive variation contributing to ASD while taking into account the potential contribution of population differences to the genetic heterogeneity of ASD. Our strategy, homozygous haplotype (HH) mapping, aims to detect homozygous segments of identical haplotype structure that are shared at a higher frequency amongst ASD patients compared to parental controls. The analysis was performed on 1,402 Autism Genome Project trios genotyped for 1 million single nucleotide polymorphisms (SNPs). We identified 25 known and 1,218 novel ASD candidate genes in the discovery analysis including CADM2, ABHD14A, CHRFAM7A, GRIK2, GRM3, EPHA3, FGF10, KCND2, PDZK1, IMMP2L and FOXP2. Furthermore, 10 of the previously reported ASD genes and 300 of the novel candidates identified in the discovery analysis were replicated in an independent sample of 1,182 trios. Our results demonstrate that regions of HH are significantly enriched for previously reported ASD candidate genes and the observed association is independent of gene size (odds ratio 2.10). Our findings highlight the applicability of HH mapping in complex disorders such as ASD and offer an alternative approach to the analysis of genome-wide association data.

Strategies for pathway analysis from GWAS data

Genome-wide association studies (GWAS) are a standard approach for investigating the relationship of common variation within the human genome to a given phenotype of interest. However, single-allele association results published for many GWAS studies represent only the tip of the iceberg for the information that can be extracted from these datasets. The primary analysis strategy for GWAS entails association analysis in which only the single nucleotide polymorphisms (SNPs) with the strongest p values are declared statistically significant due to issues arising from multiple testing and type I error concerns. Factors such as locus heterogeneity, epistasis, and multiple genes conferring small effects contribute to the complexity of the genetic models underlying phenotype expression. Thus, many biologically meaningful associations having lower effect sizes at individual genes are overlooked, as they are difficult to separate from a sea of false-positive associations. Organizing these individual SNPs into biologically meaningful groups to look at overall effects of minor perturbations to genes and pathways is desirable. This pathway-based approach provides researchers with insight into the functional foundations of the phenotype being studied and allows testing of various genetic scenarios.

Common variants near CAV1 and CAV2 are associated with primary open-angle glaucoma in Caucasians from the USA

Primary open-angle glaucoma (POAG) is a genetically complex common disease characterized by progressive optic nerve degeneration that results in irreversible blindness. Recently, a genome-wide association study (GWAS) for POAG in an Icelandic population identified significant associations with single nucleotide polymorphisms (SNPs) between the CAV1 and CAV2 genes on chromosome 7q31. In this study, we confirm that the identified SNPs are associated with POAG in our Caucasian US population and that specific haplotypes located in the CAV1/CAV2 intergenic region are associated with the disease. We also present data suggesting that associations with several CAV1/CAV2 SNPs are significant mostly in women.

Genetic analysis of biological pathway data through genomic randomization

Genome Wide Association Studies (GWAS) are a standard approach for large-scale common variation characterization and for identification of single loci predisposing to disease. However, due to issues of moderate sample sizes and particularly multiple testing correction, many variants of smaller effect size are not detected within a single allele analysis framework. Thus, small main effects and potential epistatic effects are not consistently observed in GWAS using standard analytical approaches that consider only single SNP alleles. Here, we propose unique methodology that aggregates variants of interest (for example, genes in a biological pathway) using GWAS results. Multiple testing and type I error concerns are minimized using empirical genomic randomization to estimate significance. Randomization corrects for common pathway-based analysis biases, such as SNP coverage and density, linkage disequilibrium, gene size and pathway size. Pathway Analysis by Randomization Incorporating Structure (PARIS) applies this randomization and in doing so directly accounts for linkage disequilibrium effects. PARIS is independent of association analysis method and is thus applicable to GWAS datasets of all study designs. Using the KEGG database as an example, we apply PARIS to the publicly available Autism Genetic Resource Exchange GWAS dataset, revealing pathways with a significant enrichment of positive association results.

A genome-wide scan for common alleles affecting risk for autism

Although autism spectrum disorders (ASDs) have a substantial genetic basis, most of the known genetic risk has been traced to rare variants, principally copy number variants (CNVs). To identify common risk variation, the Autism Genome Project (AGP) Consortium genotyped 1558 rigorously defined ASD families for 1 million single-nucleotide polymorphisms (SNPs) and analyzed these SNP genotypes for association with ASD. In one of four primary association analyses, the association signal for marker rs4141463, located within MACROD2, crossed the genome-wide association significance threshold of P < 5 × 10⁻⁸. When a smaller replication sample was analyzed, the risk allele at rs4141463 was again over-transmitted; yet, consistent with the winner's curse, its effect size in the replication sample was much smaller; and, for the combined samples, the association signal barely fell below the P < 5 × 10⁻⁸ threshold. Exploratory analyses of phenotypic subtypes yielded no significant associations after correction for multiple testing. They did, however, yield strong signals within several genes, KIAA0564, PLD5, POU6F2, ST8SIA2 and TAF1C.

Functional impact of global rare copy number variation in autism spectrum disorders

The autism spectrum disorders (ASDs) are a group of conditions characterized by impairments in reciprocal social interaction and communication, and the presence of restricted and repetitive behaviours. Individuals with an ASD vary greatly in cognitive development, which can range from above average to intellectual disability. Although ASDs are known to be highly heritable ( approximately 90%), the underlying genetic determinants are still largely unknown. Here we analysed the genome-wide characteristics of rare (<1% frequency) copy number variation in ASD using dense genotyping arrays. When comparing 996 ASD individuals of European ancestry to 1,287 matched controls, cases were found to carry a higher global burden of rare, genic copy number variants (CNVs) (1.19 fold, P = 0.012), especially so for loci previously implicated in either ASD and/or intellectual disability (1.69 fold, P = 3.4 x 10(-4)). Among the CNVs there were numerous de novo and inherited events, sometimes in combination in a given family, implicating many novel ASD genes such as SHANK2, SYNGAP1, DLGAP2 and the X-linked DDX53-PTCHD1 locus. We also discovered an enrichment of CNVs disrupting functional gene sets involved in cellular proliferation, projection and motility, and GTPase/Ras signalling. Our results reveal many new genetic and functional targets in ASD that may lead to final connected pathways.

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.

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.