Natural variation in the parameters of innate immune cells is preferentially driven by genetic factors.

The quantification and characterization of circulating immune cells provide key indicators of human health and disease. To identify the relative effects of environmental and genetic factors on variation in the parameters of innate and adaptive immune cells in homeostatic conditions, we combined standardized flow cytometry of blood leukocytes and genome-wide DNA genotyping of 1,000 healthy, unrelated people of Western European ancestry. We found that smoking, together with age, sex and latent infection with cytomegalovirus, were the main non-genetic factors that affected variation in parameters of human immune cells. Genome-wide association studies of 166 immunophenotypes identified 15 loci that showed enrichment for disease-associated variants. Finally, we demonstrated that the parameters of innate cells were more strongly controlled by genetic variation than were those of adaptive cells, which were driven by mainly environmental exposure. Our data establish a resource that will generate new hypotheses in immunology and highlight the role of innate immunity in susceptibility to common autoimmune diseases.

Publisher Correction: Natural variation in the parameters of innate immune cells is preferentially driven by genetic factors.

In the version of this article initially published, the name of one author was incorrect (James P. Santo). The correct name is James P. Di Santo. The error has been corrected in the HTML and PDF versions of the article.

In silico tools for accurate HLA and KIR inference from clinical sequencing data empower immunogenetics on individual-patient and population scales.

Immunogenetic variation in humans is important in research, clinical diagnosis and increasingly a target for therapeutic intervention. Two highly polymorphic loci play critical roles, namely the human leukocyte antigen (HLA) system, which is the human version of the major histocompatibility complex (MHC), and the Killer-cell immunoglobulin-like receptors (KIR) that are relevant for responses of natural killer (NK) and some subsets of T cells. Their accurate classification has typically required the use of dedicated biological specimens and a combination of in vitro and in silico efforts. Increased availability of next generation sequencing data has led to the development of ancillary computational solutions. Here, we report an evaluation of recently published algorithms to computationally infer complex immunogenetic variation in the form of HLA alleles and KIR haplotypes from whole-genome or whole-exome sequencing data. For both HLA allele and KIR gene typing, we identified tools that yielded >97% overall accuracy for four-digit HLA types, and >99% overall accuracy for KIR gene presence, suggesting the readiness of in silico solutions for use in clinical and high-throughput research settings.

Polygenic risk for skin autoimmunity impacts immune checkpoint blockade in bladder cancer.

PD-1 and PD-L1 act to restrict T cell responses in cancer and contribute to self-tolerance. Consistent with this role, PD-1 checkpoint inhibitors have been associated with immune-related adverse events (irAEs), immune toxicities thought to be autoimmune in origin. Analyses of dermatological irAEs have identified an association with improved overall survival (OS) following anti-PD-(L)1 therapy, but the factors that contribute to this relationship are poorly understood. We collected germline whole-genome sequencing data from IMvigor211, a recent phase 3 randomized controlled trial comparing atezolizumab (anti-PD-L1) monotherapy to chemotherapy in bladder cancer. We found that high vitiligo, high psoriasis, and low atopic dermatitis polygenic risk scores (PRSs) were associated with longer OS under anti-PD-L1 monotherapy as compared to chemotherapy, reflecting the Th17 polarization of these diseases. PRSs were not correlated with tumor mutation burden, PD-L1 immunohistochemistry, nor T-effector gene signatures. Shared genetic factors impact risk for dermatological autoimmunity and anti-PD-L1 monotherapy in bladder cancer.

Genome-wide association study and functional validation implicates JADE1 in tauopathy.

Primary age-related tauopathy (PART) is a neurodegenerative pathology with features distinct from but also overlapping with Alzheimer disease (AD). While both exhibit Alzheimer-type temporal lobe neurofibrillary degeneration alongside amnestic cognitive impairment, PART develops independently of amyloid-ß (Aß) plaques. The pathogenesis of PART is not known, but evidence suggests an association with genes that promote tau pathology and others that protect from Aß toxicity. Here, we performed a genetic association study in an autopsy cohort of individuals with PART (n = 647) using Braak neurofibrillary tangle stage as a quantitative trait. We found some significant associations with candidate loci associated with AD (SLC24A4, MS4A6A, HS3ST1) and progressive supranuclear palsy (MAPT and EIF2AK3). Genome-wide association analysis revealed a novel significant association with a single nucleotide polymorphism on chromosome 4 (rs56405341) in a locus containing three genes, including JADE1 which was significantly upregulated in tangle-bearing neurons by single-soma RNA-seq. Immunohistochemical studies using antisera targeting JADE1 protein revealed localization within tau aggregates in autopsy brains with four microtubule-binding domain repeats (4R) isoforms and mixed 3R/4R, but not with 3R exclusively. Co-immunoprecipitation in post-mortem human PART brain tissue revealed a specific binding of JADE1 protein to four repeat tau lacking N-terminal inserts (0N4R). Finally, knockdown of the Drosophila JADE1 homolog rhinoceros (rno) enhanced tau-induced toxicity and apoptosis in vivo in a humanized 0N4R mutant tau knock-in model, as quantified by rough eye phenotype and terminal deoxynucleotidyl transferase dUTP nick end-labeling (TUNEL) in the fly brain. Together, these findings indicate that PART has a genetic architecture that partially overlaps with AD and other tauopathies and suggests a novel role for JADE1 as a modifier of neurofibrillary degeneration.

Previously reported placebo-response-associated variants do not predict patient outcomes in inflammatory disease Phase III trial placebo arms.

In clinical trials, a placebo response refers to improvement in disease symptoms arising from the psychological effect of receiving a treatment rather than the actual treatment under investigation. Previous research has reported genomic variation associated with the likelihood of observing a placebo response, but these studies have been limited in scope and have not been validated. Here, we analyzed whole-genome sequencing data from 784 patients undergoing placebo treatment in Phase III Asthma or Rheumatoid Arthritis trials to assess the impact of previously reported variation on patient outcomes in the placebo arms and to identify novel variants associated with the placebo response. Contrary to expectations based on previous reports, we did not observe any statistically significant associations between genomic variants and placebo treatment outcome. Our findings suggest that the biological origin of the placebo response is complex and likely to be variable between disease areas.

Identifying and mitigating batch effects in whole genome sequencing data.

BACKGROUND: Large sample sets of whole genome sequencing with deep coverage are being generated, however assembling datasets from different sources inevitably introduces batch effects. These batch effects are not well understood and can be due to changes in the sequencing protocol or bioinformatics tools used to process the data. No systematic algorithms or heuristics exist to detect and filter batch effects or remove associations impacted by batch effects in whole genome sequencing data. RESULTS: We describe key quality metrics, provide a freely available software package to compute them, and demonstrate that identification of batch effects is aided by principal components analysis of these metrics. To mitigate batch effects, we developed new site-specific filters that identified and removed variants that falsely associated with the phenotype due to batch effect. These include filtering based on: a haplotype based genotype correction, a differential genotype quality test, and removing sites with missing genotype rate greater than 30% after setting genotypes with quality scores less than 20 to missing. This method removed 96.1% of unconfirmed genome-wide significant SNP associations and 97.6% of unconfirmed genome-wide significant indel associations. We performed analyses to demonstrate that: 1) These filters impacted variants known to be disease associated as 2 out of 16 confirmed associations in an AMD candidate SNP analysis were filtered, representing a reduction in power of 12.5%, 2) In the absence of batch effects, these filters removed only a small proportion of variants across the genome (type I error rate of 3%), and 3) in an independent dataset, the method removed 90.2% of unconfirmed genome-wide SNP associations and 89.8% of unconfirmed genome-wide indel associations. CONCLUSIONS: Researchers currently do not have effective tools to identify and mitigate batch effects in whole genome sequencing data. We developed and validated methods and filters to address this deficiency.

Polycomb-like 3 promotes polycomb repressive complex 2 binding to CpG islands and embryonic stem cell self-renewal.

Polycomb repressive complex 2 (PRC2) trimethylates lysine 27 of histone H3 (H3K27me3) to regulate gene expression during diverse biological transitions in development, embryonic stem cell (ESC) differentiation, and cancer. Here, we show that Polycomb-like 3 (Pcl3) is a component of PRC2 that promotes ESC self-renewal. Using mass spectrometry, we identified Pcl3 as a Suz12 binding partner and confirmed Pcl3 interactions with core PRC2 components by co-immunoprecipitation. Knockdown of Pcl3 in ESCs increases spontaneous differentiation, yet does not affect early differentiation decisions as assessed in teratomas and embryoid bodies, indicating that Pcl3 has a specific role in regulating ESC self-renewal. Consistent with Pcl3 promoting PRC2 function, decreasing Pcl3 levels reduces H3K27me3 levels while overexpressing Pcl3 increases H3K27me3 levels. Furthermore, chromatin immunoprecipitation and sequencing (ChIP-seq) reveal that Pcl3 co-localizes with PRC2 core component, Suz12, and depletion of Pcl3 decreases Suz12 binding at over 60% of PRC2 targets. Mutation of conserved residues within the Pcl3 Tudor domain, a domain implicated in recognizing methylated histones, compromises H3K27me3 formation, suggesting that the Tudor domain of Pcl3 is essential for function. We also show that Pcl3 and its paralog, Pcl2, exist in different PRC2 complexes but bind many of the same PRC2 targets, particularly CpG islands regulated by Pcl3. Thus, Pcl3 is a component of PRC2 critical for ESC self-renewal, histone methylation, and recruitment of PRC2 to a subset of its genomic sites.

Floxin, a resource for genetically engineering mouse ESCs.

We describe a method for the highly efficient and precise targeted modification of gene trap loci in mouse embryonic stem cells (ESCs). Through the Floxin method, gene trap mutations were reverted and new DNA sequences inserted using Cre recombinase and a shuttle vector, pFloxin. Floxin technology is applicable to the existing collection of 24,149 compatible gene trap cell lines, which should enable high-throughput modification of many genes in mouse ESCs.

Allelic variation in HLA-DRB1 is associated with development of antidrug antibodies in cancer patients treated with atezolizumab that are neutralizing in vitro.

The treatment of diseases with biologic agents can result in the formation of antidrug antibodies (ADA). Although drivers for ADA formation are unknown, a role for antigen presentation is likely, and variation in human leukocyte antigen (HLA) genes has been shown to be associated with occurrence of ADA for several biologics. Here, we performed an HLA-wide association study in 1982 patients treated with the anti-PD-L1 antibody atezolizumab across eight clinical trials. On average, 29.8% of patients were ADA-positive (N = 591, range of 13.5%-38.4% per study) and 14.6% of patients were positive for ADA that were neutralizing in vitro (neutralizing antibodies [NAb], N = 278, range of 6.4%-21.9% per study). In a meta-analysis of logistic regression coefficients, we found statistically significant associations between HLA class II alleles and ADA status. The top-associated alleles were HLA-DRB1*01:01 in a comparison of ADA-positive versus ADA-negative patients (p = 3.4 × 10-5 , odds ratio [OR] 1.96, 95% confidence interval [95% CI] 1.64-2.28) and HLA-DQA1*01:01 when comparing NAb-positive with ADA-negative patients (p = 2.8 × 10-7 , OR 2.31, 95% CI 1.98-2.66). Both alleles occur together on a common HLA haplotype, and analyses considering only NAb-negative, ADA-positive patients did not yield significant results, suggesting that the genetic association is mainly driven by NAb status. In conclusion, our study showed that HLA class II genotype is associated with the risk of developing ADA, and specifically NAb, in patients treated with atezolizumab, but the effect estimates suggest that immunogenetic factors are not sufficient as clinically meaningful predictors.

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.

Genetic variation associated with thyroid autoimmunity shapes the systemic immune response to PD-1 checkpoint blockade.

Activation of systemic immune responses using PD-1 checkpoint inhibitors is an essential approach to cancer therapy. Yet, the extent of benefit relative to risk of immune related adverse events (irAE) varies widely among patients. Here, we study endocrine irAE from 7 clinical trials across 6 cancers where atezolizumab (anti-PD-L1) was combined with chemotherapies and compared to standard of care. We show that atezolizumab-induced thyroid dysfunction is associated with longer survival. We construct a polygenic risk score (PRS) for lifetime risk of hypothyroidism using a GWAS from the UK Biobank and apply this PRS to genetic data collected from 2,616 patients of European ancestry from these trials. Patients with high PRS are at increased risk of atezolizumab-induced thyroid dysfunction and lower risk of death in triple negative breast cancer. Our results indicate that genetic variation associated with thyroid autoimmunity interacts with biological pathways driving the systemic immune response to PD-1 blockade.

Enabling genome-wide association testing with multiple diseases and no healthy controls.

MOTIVATION: While large-scale whole genome sequencing is feasible the high costs compel investigators to focus on disease subjects. As a result large sequencing datasets of samples with different diseases are often readily available, but not healthy controls to contrast them with. While it is possible to perform an association study using only diseases, the associations could be driven by a disease acting as a control and not the focal disease. METHODS: We developed a genotype-on-phenotype reverse regression with a Bayesian spike and slab prior to enable association testing in datasets with multiple diseases. This method, referred to as revreg, flagged associations (both common and rare) that were driven by diseases that were not of primary interest. RESULTS: Based on simulations, revreg had 80% power to detect an odds ratio of 1.74 for common variants (3500 samples total) and 3.73 for rare variants (14,000 samples total), with minimal type I error. For common variants, we tested this method on 3657 whole genome sequenced samples aimed at discovering variants associated with disease risk of Chronic Obstructive Pulmonary Disease using three other diseases as controls. We demonstrated detection of six highly significant associations likely due to Age-Related Macular Degeneration. In an exome dataset of 8836 samples aimed at characterizing rare variants associated with disease risk of Asthma, using five other diseases as controls, we detected and removed genic regions due to AMD (C3, CFH, CFHR5, CFI, and DNMT3A) and RA (KRTAP13-4).

Influence of genetic copy number variants of the human GLUT3 glucose transporter gene SLC2A3 on protein expression, glycolysis and rheumatoid arthritis risk: A genetic replication study.

OBJECTIVES: The gene encoding glucose transporter 3 (GLUT3, SLC2A3) is present in the human population at variable copy number. An overt disease phenotype of SLC2A3 copy number variants has not been reported; however, deletion of SLC2A3 has been previously reported to protect carriers from rheumatoid arthritis, implicating GLUT3 as a therapeutic target in rheumatoid arthritis. Here we aim to perform functional analysis of GLUT3 copy number variants in immune cells, and test the reported protective association of the GLUT3 copy number variants for rheumatoid arthritis in a genetic replication study. METHODS: Cells from genotyped healthy controls were analyzed for SLC2A3/GLUT3 expression and glycolysis capacity. We genotyped the SLC2A3 copy number variant in four independent cohorts of rheumatoid arthritis and controls and one cohort of multiple sclerosis and controls. RESULTS: Heterozygous deletion of SLC2A3 correlates directly with expression levels of GLUT3 and influences glycolysis rates in the human immune system. The frequency of the SLC2A3 copy number variant is not different between rheumatoid arthritis, multiple sclerosis and control groups. CONCLUSIONS: Despite a robust SLC2A3 gene copy number dependent phenotype, our study of large groups of rheumatoid arthritis cases and controls provides no evidence for rheumatoid arthritis disease protection in deletion carriers. These data emphasize the importance of well powered replication studies to confirm or refute genetic associations, particularly for relatively rare variants.

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.

A meta-analysis of genome-wide association studies identifies 17 new Parkinson's disease risk loci.

Common variant genome-wide association studies (GWASs) have, to date, identified >24 risk loci for Parkinson's disease (PD). To discover additional loci, we carried out a GWAS comparing 6,476 PD cases with 302,042 controls, followed by a meta-analysis with a recent study of over 13,000 PD cases and 95,000 controls at 9,830 overlapping variants. We then tested 35 loci (P < 1 × 10-6) in a replication cohort of 5,851 cases and 5,866 controls. We identified 17 novel risk loci (P < 5 × 10-8) in a joint analysis of 26,035 cases and 403,190 controls. We used a neurocentric strategy to assign candidate risk genes to the loci. We identified protein-altering or cis-expression quantitative trait locus (cis-eQTL) variants in linkage disequilibrium with the index variant in 29 of the 41 PD loci. These results indicate a key role for autophagy and lysosomal biology in PD risk, and suggest potential new drug targets for PD.

A rare mutation in UNC5C predisposes to late-onset Alzheimer's disease and increases neuronal cell death.

We have identified a rare coding mutation, T835M (rs137875858), in the UNC5C netrin receptor gene that segregated with disease in an autosomal dominant pattern in two families enriched for late-onset Alzheimer's disease and that was associated with disease across four large case-control cohorts (odds ratio = 2.15, Pmeta = 0.0095). T835M alters a conserved residue in the hinge region of UNC5C, and in vitro studies demonstrate that this mutation leads to increased cell death in human HEK293T cells and in rodent neurons. Furthermore, neurons expressing T835M UNC5C are more susceptible to cell death from multiple neurotoxic stimuli, including ß-amyloid (Aß), glutamate and staurosporine. On the basis of these data and the enriched hippocampal expression of UNC5C in the adult nervous system, we propose that one possible mechanism in which T835M UNC5C contributes to the risk of Alzheimer's disease is by increasing susceptibility to neuronal cell death, particularly in vulnerable regions of the Alzheimer's disease brain.

The ciliogenic protein Oral-Facial-Digital 1 regulates the neuronal differentiation of embryonic stem cells.

Oral-Facial-Digital 1 (OFD1) Syndrome is an X-linked developmental disorder caused by mutations in the gene Ofd1. OFD1 syndrome involves malformation of the face, oral cavity, and digits and may be characterized by cystic kidneys and mental retardation. Deletion or missense mutations in Ofd1 also result in loss of primary cilia, a microtubule-based cellular projection that mediates multiple signaling pathways. Ofd1 mutant mice display pleiotropic developmental phenotypes, including neural, skeletal, and cardiac defects. To address how loss of Ofd1 and loss of primary cilia affect early differentiation decisions, we analyzed embryoid bodies (EBs) derived from Ofd1 mutant embryonic stem (ES) cells. Ofd1 mutant EBs do not form primary cilia and display defects in Hedgehog and Wnt signaling. Additionally, we show that ES cells lacking Ofd1 display an increased capacity to differentiate into neurons. Nevertheless, neurons derived from Ofd1 mutant ES cells fail to differentiate into V3 interneurons, a cell type dependent on ciliary function and Hedgehog signaling. Thus, loss of Ofd1 affects ES cell interpretation of developmental cues and reveals that EBs model some aspects of ciliopathies, providing insights into the developmental origins of OFD1 syndrome and functions of cilia.

Polycomb-like 2 associates with PRC2 and regulates transcriptional networks during mouse embryonic stem cell self-renewal and differentiation.

Polycomb group (PcG) proteins are conserved epigenetic transcriptional repressors that control numerous developmental gene expression programs and have recently been implicated in modulating embryonic stem cell (ESC) fate. We identified the PcG protein PCL2 (polycomb-like 2) in a genome-wide screen for regulators of self-renewal and pluripotency and predicted that it would play an important role in mouse ESC-fate determination. Using multiple biochemical strategies, we provide evidence that PCL2 is a Polycomb Repressive Complex 2 (PRC2)-associated protein in mouse ESCs. Knockdown of Pcl2 in ESCs resulted in heightened self-renewal characteristics, defects in differentiation, and altered patterns of histone methylation. Integration of global gene expression and promoter occupancy analyses allowed us to identify PCL2 and PRC2 transcriptional targets and draft regulatory networks. We describe the role of PCL2 in both modulating transcription of ESC self-renewal genes in undifferentiated ESCs as well as developmental regulators during early commitment and differentiation.