Detailed knowledge of how diversity in the sequence of the human genome affects phenotypic diversity depends on a comprehensive and reliable characterization of both sequences and phenotypic variation. Over the past decade, insights into this relationship have been obtained from whole-exome sequencing or whole-genome sequencing of large cohorts with rich phenotypic data1,2. Here we describe the analysis of whole-genome sequencing of 150,119 individuals from the UK Biobank3. This constitutes a set of high-quality variants, including 585,040,410 single-nucleotide polymorphisms, representing 7.0% of all possible human single-nucleotide polymorphisms, and 58,707,036 indels. This large set of variants allows us to characterize selection based on sequence variation within a population through a depletion rank score of windows along the genome. Depletion rank analysis shows that coding exons represent a small fraction of regions in the genome subject to strong sequence conservation. We define three cohorts within the UK Biobank: a large British Irish cohort, a smaller African cohort and a South Asian cohort. A haplotype reference panel is provided that allows reliable imputation of most variants carried by three or more sequenced individuals. We identified 895,055 structural variants and 2,536,688 microsatellites, groups of variants typically excluded from large-scale whole-genome sequencing studies. Using this formidable new resource, we provide several examples of trait associations for rare variants with large effects not found previously through studies based on whole-exome sequencing and/or imputation.

SARS-CoV-2 vaccines are crucial in controlling COVID-19, but knowledge of which factors determine waning immunity is limited. We examined antibody levels and T-cell gamma-interferon release after two doses of BNT162b2 vaccine or a combination of ChAdOx1-nCoV19 and BNT162b2 vaccines for up to 230 days after the first dose. Generalized mixed models with and without natural cubic splines were used to determine immunity over time. Antibody responses were influenced by natural infection, sex, and age. IgA only became significant in naturally infected. A one-year IgG projection suggested an initial two-phase response in those given the second dose delayed (ChAdOx1/BNT162b2) followed by a more rapid decrease of antibody levels. T-cell responses correlated significantly with IgG antibody responses. Our results indicate that IgG levels will drop at different rates depending on prior infection, age, sex, T-cell response, and the interval between vaccine injections. Only natural infection mounted a significant and lasting IgA response.

Significance There is growing evidence that preexisting autoantibodies neutralizing type I interferons (IFNs) are strong determinants of life-threatening COVID-19 pneumonia. It is important to estimate their quantitative impact on COVID-19 mortality upon SARS-CoV-2 infection, by age and sex, as both the prevalence of these autoantibodies and the risk of COVID-19 death increase with age and are higher in men. Using an unvaccinated sample of 1,261 deceased patients and 34,159 individuals from the general population, we found that autoantibodies against type I IFNs strongly increased the SARS-CoV-2 infection fatality rate at all ages, in both men and women. Autoantibodies against type I IFNs are strong and common predictors of life-threatening COVID-19. Testing for these autoantibodies should be considered in the general population.

Back pain is a common and debilitating disorder with largely unknown underlying biology. Here we report a genome-wide association study of back pain using diagnoses assigned in clinical practice; dorsalgia (119,100 cases, 909,847 controls) and intervertebral disc disorder (IDD) (58,854 cases, 922,958 controls). We identify 41 variants at 33 loci. The most significant association (ORIDD = 0.92, P = 1.6 × 10−39; ORdorsalgia = 0.92, P = 7.2 × 10−15) is with a 3’UTR variant (rs1871452-T) in CHST3, encoding a sulfotransferase enzyme expressed in intervertebral discs. The largest effects on IDD are conferred by rare (MAF = 0.07 − 0.32%) loss-of-function (LoF) variants in SLC13A1, encoding a sodium-sulfate co-transporter (LoF burden OR = 1.44, P = 3.1 × 10−11); variants that also associate with reduced serum sulfate. Genes implicated by this study are involved in cartilage and bone biology, as well as neurological and inflammatory processes.