A deep catalogue of protein-coding variation in 983,578 individuals.

Rare coding variants that significantly impact function provide insights into the biology of a gene1-3. However, ascertaining their frequency requires large sample sizes4-8. Here, we present a catalogue of human protein-coding variation, derived from exome sequencing of 983,578 individuals across diverse populations. 23% of the Regeneron Genetics Center Million Exome data (RGC-ME) comes from non-European individuals of African, East Asian, Indigenous American, Middle Eastern, and South Asian ancestry. This catalogue includes over 10.4 million missense and 1.1 million predicted loss-of-function (pLOF) variants. We identify individuals with rare biallelic pLOF variants in 4,848 genes, 1,751 of which have not been previously reported. From precise quantitative estimates of selection against heterozygous loss-of-function, we identify 3,988 loss-of-function intolerant genes, including 86 that were previously assessed as tolerant and 1,153 lacking established disease annotation. We also define regions of missense depletion at high resolution. Notably, 1,482 genes have regions depleted of missense variants despite being tolerant to pLOF variants. Finally, we estimate that 3% of individuals have a clinically actionable genetic variant, and that 11,773 variants reported in ClinVar with unknown significance are likely to be deleterious cryptic splice sites. To facilitate variant interpretation and genetics-informed precision medicine, we make this important resource of coding variation from the RGC-ME accessible via a public variant allele frequency browser.

A deep catalog of protein-coding variation in 985,830 individuals.

Coding variants that have significant impact on function can provide insights into the biology of a gene but are typically rare in the population. Identifying and ascertaining the frequency of such rare variants requires very large sample sizes. Here, we present the largest catalog of human protein-coding variation to date, derived from exome sequencing of 985,830 individuals of diverse ancestry to serve as a rich resource for studying rare coding variants. Individuals of African, Admixed American, East Asian, Middle Eastern, and South Asian ancestry account for 20% of this Exome dataset. Our catalog of variants includes approximately 10.5 million missense (54% novel) and 1.1 million predicted loss-of-function (pLOF) variants (65% novel, 53% observed only once). We identified individuals with rare homozygous pLOF variants in 4,874 genes, and for 1,838 of these this work is the first to document at least one pLOF homozygote. Additional insights from the RGC-ME dataset include 1) improved estimates of selection against heterozygous loss-of-function and identification of 3,459 genes intolerant to loss-of-function, 83 of which were previously assessed as tolerant to loss-of-function and 1,241 that lack disease annotations; 2) identification of regions depleted of missense variation in 457 genes that are tolerant to loss-of-function; 3) functional interpretation for 10,708 variants of unknown or conflicting significance reported in ClinVar as cryptic splice sites using splicing score thresholds based on empirical variant deleteriousness scores derived from RGC-ME; and 4) an observation that approximately 3% of sequenced individuals carry a clinically actionable genetic variant in the ACMG SF 3.1 list of genes. We make this important resource of coding variation available to the public through a variant allele frequency browser. We anticipate that this report and the RGC-ME dataset will serve as a valuable reference for understanding rare coding variation and help advance precision medicine efforts.

Genome-wide analysis provides genetic evidence that ACE2 influences COVID-19 risk and yields risk scores associated with severe disease.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) enters human host cells via angiotensin-converting enzyme 2 (ACE2) and causes coronavirus disease 2019 (COVID-19). Here, through a genome-wide association study, we identify a variant (rs190509934, minor allele frequency 0.2-2%) that downregulates ACE2 expression by 37% (P = 2.7 × 10-8) and reduces the risk of SARS-CoV-2 infection by 40% (odds ratio = 0.60, P = 4.5 × 10-13), providing human genetic evidence that ACE2 expression levels influence COVID-19 risk. We also replicate the associations of six previously reported risk variants, of which four were further associated with worse outcomes in individuals infected with the virus (in/near LZTFL1, MHC, DPP9 and IFNAR2). Lastly, we show that common variants define a risk score that is strongly associated with severe disease among cases and modestly improves the prediction of disease severity relative to demographic and clinical factors alone.

xAtlas: scalable small variant calling across heterogeneous next-generation sequencing experiments.

BACKGROUND: The growing volume and heterogeneity of next-generation sequencing (NGS) data complicate the further optimization of identifying DNA variation, especially considering that curated high-confidence variant call sets frequently used to validate these methods are generally developed from the analysis of comparatively small and homogeneous sample sets. FINDINGS: We have developed xAtlas, a single-sample variant caller for single-nucleotide variants (SNVs) and small insertions and deletions (indels) in NGS data. xAtlas features rapid runtimes, support for CRAM and gVCF file formats, and retraining capabilities. xAtlas reports SNVs with 99.11% recall and 98.43% precision across a reference HG002 sample at 60× whole-genome coverage in less than 2 CPU hours. Applying xAtlas to 3,202 samples at 30× whole-genome coverage from the 1000 Genomes Project achieves an average runtime of 1.7 hours per sample and a clear separation of the individual populations in principal component analysis across called SNVs. CONCLUSIONS: xAtlas is a fast, lightweight, and accurate SNV and small indel calling method. Source code for xAtlas is available under a BSD 3-clause license at https://github.com/jfarek/xatlas.

Pan-ancestry exome-wide association analyses of COVID-19 outcomes in 586,157 individuals.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19), a respiratory illness that can result in hospitalization or death. We used exome sequence data to investigate associations between rare genetic variants and seven COVID-19 outcomes in 586,157 individuals, including 20,952 with COVID-19. After accounting for multiple testing, we did not identify any clear associations with rare variants either exome wide or when specifically focusing on (1) 13 interferon pathway genes in which rare deleterious variants have been reported in individuals with severe COVID-19, (2) 281 genes located in susceptibility loci identified by the COVID-19 Host Genetics Initiative, or (3) 32 additional genes of immunologic relevance and/or therapeutic potential. Our analyses indicate there are no significant associations with rare protein-coding variants with detectable effect sizes at our current sample sizes. Analyses will be updated as additional data become available, and results are publicly available through the Regeneron Genetics Center COVID-19 Results Browser.

PARP Targeted Alpha-Particle Therapy Enhances Response to PD-1 Immune-Checkpoint Blockade in a Syngeneic Mouse Model of Glioblastoma.

We have previously demonstrated potent antitumor effects of PARP targeted alpha-therapy with astatine-211-MM4 ([211At]MM4) in neuroblastoma preclinical models, although differential sensitivity suggests it is unlikely to be curative as a single-agent in all tumor types. Alpha-particle induced DNA damage can elicit an immune response that results in T-cell activation against tumor cells; however, tumor cells can evade immune surveillance through expression of programmed death ligand 1 (PD-L1). Therefore, we investigated the effects of α particle therapy in combination with immune-checkpoint blockade using astatine-211-MM4 and anti-programmed death receptor 1 (anti-PD-1) immunotherapy in a syngeneic mouse model of glioblastoma. We characterized the sensitivity of four human glioblastoma cell lines to [211At]MM4 in vitro. To evaluate [211At]MM4 treatment effects on hematological tissues, complete blood counts were performed after a single dose at 12, 24, or 36 MBq/kg. In vivo efficacy was evaluated in a syngeneic mouse model of glioblastoma using GL26 glioblastoma cells in CB57BL/6J mice treated with either 36 MBq/kg [211At]MM4, anti-PD-1 antibody, or a combination of the two. Following a single dose of [211At]MM4, lymphocytes are significantly decreased compared to control at both 72 h and 1 week following treatment followed by recovery of counts by 2 weeks. However, neutrophils showed an increase with all dose levels of [211At]MM4 exhibiting higher levels than control. The average best tumor responses for combination, anti-PD-1, and [211At]MM4 were 100%, 83.6%, and 58.2% decrease in tumor volume, respectively. Average progression free intervals for combination, anti-PD-1, [211At]MM4, and control groups was 65, 36.4, 23.2, and 3 days, respectively. The percentages of disease-free mice at the end of the study for combination and anti-PD-1 were 100% and 60%, while [211At]MM4 and control groups were both 0%. In summary, combination therapy was more effective than either single agent in all response categories analyzed, highlighting the potential for PARP targeted alpha-therapy to enhance PD-1 immune-checkpoint blockade.

PARP-1-Targeted Auger Emitters Display High-LET Cytotoxic Properties In Vitro but Show Limited Therapeutic Utility in Solid Tumor Models of Human Neuroblastoma.

The currently available therapeutic radiopharmaceutical for high-risk neuroblastoma, 131I-metaiodobenzylguanidine, is ineffective at targeting micrometastases because of the low-linear-energy-transfer (LET) properties of high-energy ß-particles. In contrast, Auger radiation has high-LET properties with nanometer ranges in tissue, efficiently causing DNA damage when emitted near DNA. The aim of this study was to evaluate the cytotoxicity of targeted Auger therapy in preclinical models of high-risk neuroblastoma. Methods: We used a radiolabled poly(adenosine diphosphate ribose) polymerase (PARP) inhibitor called 125I-KX1 to deliver Auger radiation to PARP-1, a chromatin-binding enzyme overexpressed in neuroblastoma. The in vitro cytotoxicity of 125I-KX1 was assessed in 19 neuroblastoma cell lines, followed by in-depth pharmacologic analysis in a sensitive and resistant pair of cell lines. Immunofluorescence microscopy was used to characterize 125I-KX1-induced DNA damage. Finally, in vitro and in vivo microdosimetry was modeled from experimentally derived pharmacologic variables. Results:125I-KX1 was highly cytotoxic in vitro across a panel of neuroblastoma cell lines, directly causing double-strand DNA breaks. On the basis of subcellular dosimetry, 125I-KX1 was approximately twice as effective as 131I-KX1, whereas cytoplasmic 125I-metaiodobenzylguanidine demonstrated low biological effectiveness. Despite the ability to deliver a focused radiation dose to the cell nuclei, 125I-KX1 remained less effective than its α-emitting analog 211At-MM4 and required significantly higher activity for equivalent in vivo efficacy based on tumor microdosimetry. Conclusion: Chromatin-targeted Auger therapy is lethal to high-risk neuroblastoma cells and has the potential to be used in micrometastatic disease. This study provides the first evidence for cellular lethality from a PARP-1-targeted Auger emitter, calling for further investigation into targeted Auger therapy.

Differential Selective Pressures Experienced by the Aurora Kinase Gene Family.

Aurora kinases (AKs) are serine/threonine kinases that are essential for cell division. Humans have three AK genes: AKA, AKB, and AKC. AKA is required for centrosome assembly, centrosome separation, and bipolar spindle assembly, and its mutation leads to abnormal spindle morphology. AKB is required for the spindle checkpoint and proper cytokinesis, and mutations cause chromosome misalignment and cytokinesis failure. AKC is expressed in germ cells, and has a role in meiosis analogous to that of AKB in mitosis. Mutation of any of the three isoforms can lead to cancer. AK proteins possess divergent N- and C-termini and a conserved central catalytic domain. We examined the evolution of the AK gene family using an identity matrix and by building a phylogenetic tree. The data suggest that AKA is the vertebrate ancestral gene, and that AKB and AKC resulted from gene duplication in placental mammals. In a nonsynonymous/synonymous rate substitution analysis, we found that AKB experienced the strongest, and AKC the weakest, purifying selection. Both the N- and C-termini and regions within the kinase domain experienced differential selection among the AK isoforms. These differentially selected sequences may be important for species specificity and isoform specificity, and are therefore potential therapeutic targets.