While the recurrent 22q11.2 deletion is one of the strongest genetic risk factors for schizophrenia (SCZ), variability of its associated neuropsychiatric endophenotypes reflects its incomplete penetrance for psychosis development. To assess whether this phenotypic variability is linked to common variants associated with SCZ, we studied the association between SCZ polygenic risk score (PRS) and longitudinally acquired phenotypic information of the Swiss 22q11.2DS cohort (n = 97, 50% females, mean age 17.7 yr, mean visit interval 3.8 yr). The SCZ PRS with the best predictive performance was ascertained in the Estonian Biobank (n = 201,146) with LDpred. The infinitesimal SCZ PRS model showed the strongest capacity in discriminating SCZ cases from controls with one SD difference in SCZ PRS corresponding to an odds ratio (OR) of 1.73 (95% CI 1.57-1.90, P = 1.47 × 10-29). In 22q11.2 patients, random-effects ordinal regression modelling using longitudinal data showed SCZ PRS to have the strongest effect on social anhedonia (OR = 2.09, P = 0.0002), and occupational functioning (OR = 1.82, P = 0.0003) within the negative symptoms course, and dysphoric mood (OR = 2.00, P = 0.002) and stress intolerance (OR = 1.76, P = 0.0002) within the general symptoms course. Genetic liability for SCZ was additionally associated with full scale cognitive decline (β = -0.25, P = 0.02) and with longitudinal volumetric reduction of the right and left hippocampi (β = -0.28, P = 0.005; β = -0.23, P = 0.02, respectively). Our results indicate that the polygenic contribution to SCZ acts upon the threshold-lowering first hit (i.e., the deletion). It modifies the endophenotypes of 22q11.2DS and augments the derailment of developmental trajectories of negative and general symptoms, cognition, and hippocampal volume.

Hearing loss is one of the top contributors to years lived with disability and is a risk factor for dementia. Molecular evidence on the cellular origins of hearing loss in humans is growing. Here, we performed a genome-wide association meta-analysis of clinically diagnosed and self-reported hearing impairment on 723,266 individuals and identified 48 significant loci, 10 of which are novel. A large proportion of associations comprised missense variants, half of which lie within known familial hearing loss loci. We used single-cell RNA-sequencing data from mouse cochlea and brain and mapped common-variant genomic results to spindle, root, and basal cells from the stria vascularis, a structure in the cochlea necessary for normal hearing. Our findings indicate the importance of the stria vascularis in the mechanism of hearing impairment, providing future paths for developing targets for therapeutic intervention in hearing loss.

Elevated serum urate levels, a complex trait and major risk factor for incident gout, are correlated with cardiometabolic traits via incompletely understood mechanisms. DNA methylation in whole blood captures genetic and environmental influences and is assessed in transethnic meta-analysis of epigenome-wide association studies (EWAS) of serum urate (discovery, n = 12,474, replication, n = 5522). The 100 replicated, epigenome-wide significant (p < 1.1E-7) CpGs explain 11.6% of the serum urate variance. At SLC2A9, the serum urate locus with the largest effect in genome-wide association studies (GWAS), five CpGs are associated with SLC2A9 gene expression. Four CpGs at SLC2A9 have significant causal effects on serum urate levels and/or gout, and two of these partly mediate the effects of urate-associated GWAS variants. In other genes, including SLC7A11 and PHGDH, 17 urate-associated CpGs are associated with conditions defining metabolic syndrome, suggesting that these CpGs may represent a blood DNA methylation signature of cardiometabolic risk factors. This study demonstrates that EWAS can provide new insights into GWAS loci and the correlation of serum urate with other complex traits.

Atopic dermatitis (AD) is a common chronic inflammatory skin disease with high heritability. Previous genome-wide association studies have identified several loci predisposing to AD. These findings explain approximately 30% of the variance in AD susceptibility, suggesting that further work is required to fully understand the genetic underpinnings.

Circulating proteins have important functions in inflammation and a broad range of diseases. To identify genetic influences on inflammation-related proteins, we conducted a genome-wide protein quantitative trait locus (pQTL) study of 91 plasma proteins measured using the Olink Target platform in 14,824 participants. We identified 180 pQTLs (59 cis, 121 trans). Integration of pQTL data with eQTL and disease genome-wide association studies provided insight into pathogenesis, implicating lymphotoxin-α in multiple sclerosis. Using Mendelian randomization (MR) to assess causality in disease etiology, we identified both shared and distinct effects of specific proteins across immune-mediated diseases, including directionally discordant effects of CD40 on risk of rheumatoid arthritis versus multiple sclerosis and inflammatory bowel disease. MR implicated CXCL5 in the etiology of ulcerative colitis (UC) and we show elevated gut CXCL5 transcript expression in patients with UC. These results identify targets of existing drugs and provide a powerful resource to facilitate future drug target prioritization.

Osteoporotic fracture is among the most common and costly of diseases. While reasonably heritable, its genetic determinants have remained elusive. Forearm fractures are the most common clinically recognized osteoporotic fractures with a relatively high heritability. To establish an atlas of the genetic determinants of forearm fractures, we performed genome-wide association analyses including 100,026 forearm fracture cases. We identified 43 loci, including 26 new fracture loci. Although most fracture loci associated with bone mineral density, we also identified loci that primarily regulate bone quality parameters. Functional studies of one such locus, at TAC4, revealed that Tac4-/- mice have reduced mechanical bone strength. The strongest forearm fracture signal, at WNT16, displayed remarkable bone-site-specificity with no association with hip fractures. Tall stature and low body mass index were identified as new causal risk factors for fractures. The insights from this atlas may improve fracture prediction and enable therapeutic development to prevent fractures.

Rare copy-number variants (rCNVs) include deletions and duplications that occur infrequently in the global human population and can confer substantial risk for disease. In this study, we aimed to quantify the properties of haploinsufficiency (i.e., deletion intolerance) and triplosensitivity (i.e., duplication intolerance) throughout the human genome. We harmonized and meta-analyzed rCNVs from nearly one million individuals to construct a genome-wide catalog of dosage sensitivity across 54 disorders, which defined 163 dosage sensitive segments associated with at least one disorder. These segments were typically gene dense and often harbored dominant dosage sensitive driver genes, which we were able to prioritize using statistical fine-mapping. Finally, we designed an ensemble machine-learning model to predict probabilities of dosage sensitivity (pHaplo & pTriplo) for all autosomal genes, which identified 2,987 haploinsufficient and 1,559 triplosensitive genes, including 648 that were uniquely triplosensitive. This dosage sensitivity resource will provide broad utility for human disease research and clinical genetics.

Glaucoma is a leading cause of blindness. Current glaucoma medications work by lowering intraocular pressure (IOP), a risk factor for glaucoma, but most treatments do not directly target the pathological changes leading to increased IOP, which can manifest as medication resistance as disease progresses. To identify physiological modulators of IOP, we performed genome- and exome-wide association analysis in >129,000 individuals with IOP measurements and extended these findings to an analysis of glaucoma risk. We report the identification and functional characterization of rare coding variants (including loss-of-function variants) in ANGPTL7 associated with reduction in IOP and glaucoma protection. We validated the human genetics findings in mice by establishing that Angptl7 knockout mice have lower (~2 mmHg) basal IOP compared to wild-type, with a trend towards lower IOP also in heterozygotes. Conversely, increasing murine Angptl7 levels via injection into mouse eyes increases the IOP. We also show that acute Angptl7 silencing in adult mice lowers the IOP (~2-4 mmHg), reproducing the observations in knockout mice. Collectively, our data suggest that ANGPTL7 is important for IOP homeostasis and is amenable to therapeutic modulation to help maintain a healthy IOP that can prevent onset or slow the progression of glaucoma.

The pathogenesis of type 2 diabetes (T2D) might change with increasing age. Here, we used a stratification based on age of diagnosis to gain insight into the genetics and causal risk factors of T2D across different age-groups. We performed genome-wide association studies (GWAS) on T2D and T2D subgroups based on age of diagnosis (<50, 50-60, 60-70, and >70 years) (total of 24,986 cases). As control subjects, participants were at least 70 years of age at the end of follow-up without developing T2D (N =187,130). GWAS identified 208 independent lead single nucleotide polymorphism (SNPs) mapping to 69 loci associated with T2D (P < 1.0e-8). Among others, SNPs mapped to CDKN2B-AS1 and multiple independent SNPs mapped to TCF7L2 were more strongly associated with cases diagnosed after age 70 years than with cases diagnosed before age 50 years. Based on the different case groups, we performed two-sample Mendelian randomization. Most notably, we observed that of the investigated risk factors, the association between BMI and T2D attenuated with increasing age of diagnosis. Collectively, our results indicate that stratification of T2D based on age of diag-nosis reveals subgroup-specific genetics and causal determinants, supporting the hypothesis that the pathogenesis of T2D changes with increasing age.

Understanding factors associated with COVID-19 vaccination can highlight issues in public health systems. Using machine learning, we considered the effects of 2,890 health, socio-economic and demographic factors in the entire Finnish population aged 30-80 and genome-wide information from 273,765 individuals. The strongest predictors of vaccination status were labour income and medication purchase history. Mental health conditions and having unvaccinated first-degree relatives were associated with reduced vaccination. A prediction model combining all predictors achieved good discrimination (area under the receiver operating characteristic curve, 0.801; 95% confidence interval, 0.799-0.803). The 1% of individuals with the highest predicted risk of not vaccinating had an observed vaccination rate of 18.8%, compared with 90.3% in the study population. We identified eight genetic loci associated with vaccination uptake and derived a polygenic score, which was a weak predictor in an independent subset. Our results suggest that individuals at higher risk of suffering the worst consequences of COVID-19 are also less likely to vaccinate.

Diagnostic criteria for major depressive disorder allow for heterogeneous symptom profiles but genetic analysis of major depressive symptoms has the potential to identify clinical and aetiological subtypes. There are several challenges to integrating symptom data from genetically-informative cohorts, such as sample size differences between clinical and community cohorts and various patterns of missing data. We conducted genome-wide association studies of major depressive symptoms in three clinical cohorts that were enriched for affected participants (Psychiatric Genomics Consortium, Australian Genetics of Depression Study, Generation Scotland) and three community cohorts (Avon Longitudinal Study of Parents and Children, Estonian Biobank, and UK Biobank). We fit a series of confirmatory factor models with factors that accounted for how symptom data was sampled and then compared alternative models with different symptom factors. The best fitting model had a distinct factor for Appetite/Weight symptoms and an additional measurement factor that accounted for missing data patterns in the community cohorts (use of Depression and Anhedonia as gating symptoms). The results show the importance of assessing the directionality of symptoms (such as hypersomnia versus insomnia) and of accounting for study and measurement design when meta-analysing genetic association data.

Atopic dermatitis (AD) is a common inflammatory skin condition and prior genome-wide association studies (GWAS) have identified 71 associated loci. In the current study we conducted the largest AD GWAS to date (discovery N = 1,086,394, replication N = 3,604,027), combining previously reported cohorts with additional available data. We identified 81 loci (29 novel) in the European-only analysis (which all replicated in a separate European analysis) and 10 additional loci in the multi-ancestry analysis (3 novel). Eight variants from the multi-ancestry analysis replicated in at least one of the populations tested (European, Latino or African), while two may be specific to individuals of Japanese ancestry. AD loci showed enrichment for DNAse I hypersensitivity and eQTL associations in blood. At each locus we prioritised candidate genes by integrating multi-omic data. The implicated genes are predominantly in immune pathways of relevance to atopic inflammation and some offer drug repurposing opportunities.

Thyroid hormones play a critical role in regulation of multiple physiological functions and thyroid dysfunction is associated with substantial morbidity. Here, we use electronic health records to undertake a genome-wide association study of thyroid-stimulating hormone (TSH) levels, with a total sample size of 247,107. We identify 158 novel genetic associations, more than doubling the number of known associations with TSH, and implicate 112 putative causal genes, of which 76 are not previously implicated. A polygenic score for TSH is associated with TSH levels in African, South Asian, East Asian, Middle Eastern and admixed American ancestries, and associated with hypothyroidism and other thyroid disease in South Asians. In Europeans, the TSH polygenic score is associated with thyroid disease, including thyroid cancer and age-of-onset of hypothyroidism and hyperthyroidism. We develop pathway-specific genetic risk scores for TSH levels and use these in phenome-wide association studies to identify potential consequences of pathway perturbation. Together, these findings demonstrate the potential utility of genetic associations to inform future therapeutics and risk prediction for thyroid diseases.

Copy-number variations (CNVs) have been associated with rare and debilitating genomic disorders (GDs) but their impact on health later in life in the general population remains poorly described.

Natural killer (NK) cells contribute to the development of obesity-associated insulin resistance. We demonstrate that in mice obesity promotes expansion of a distinct, interleukin-6 receptor (IL6R)a-expressing NK subpopulation, which also expresses a number of other myeloid lineage genes such as the colony-stimulating factor 1 receptor (Csf1r). Selective ablation of this Csf1r-expressing NK cell population prevents obesity and insulin resistance. Moreover, conditional inactivation of IL6Ra or Stat3 in NK cells limits obesity-associated formation of these myeloid signature NK cells, protecting from obesity, insulin resistance, and obesity-associated inflammation. Also in humans IL6Ra+ NK cells increase in obesity and correlate with markers of systemic low-grade inflammation, and their gene expression profile overlaps with characteristic gene sets of NK cells in obese mice. Collectively, we demonstrate that obesity-associated inflammation and metabolic disturbances depend on interleukin-6/Stat3-dependent formation of a distinct NK population, which may provide a target for the treatment of obesity, metaflammation-associated pathologies, and diabetes.

Acute kidney injury (AKI) is a frequent and critical complication in the clinical setting. In rodents, AKI can be effectively prevented through caloric restriction (CR), which has also been shown to increase lifespan in many species. In Caenorhabditis elegans (C. elegans), longevity studies revealed that a marked CR-induced reduction of endocannabinoids may be a key mechanism. Thus, we hypothesized that regulation of endocannabinoids, particularly arachidonoyl ethanolamide (AEA), might also play a role in CR-mediated protection from renal ischemia-reperfusion injury (IRI) in mammals including humans. In male C57Bl6J mice, CR significantly reduced renal IRI and led to a significant decrease of AEA. Supplementation of AEA to near-normal serum concentrations by repetitive intraperitoneal administration in CR mice, however, did not abrogate the protective effect of CR. We also analyzed serum samples taken before and after CR from patients of three different pilot trials of dietary interventions. In contrast to mice and C. elegans, we detected an increase of AEA. We conclude that endocannabinoid levels in mice are modulated by CR, but CR-mediated renal protection does not depend on this effect. Moreover, our results indicate that modulation of endocannabinoids by CR in humans may differ fundamentally from the effects in animal models.

To regulate food intake, our brain constantly integrates external cues, such as the incentive value of a potential food reward, with internal state signals, such as hunger feelings. Incentive motivation refers to the processes that translate an expected reward into the effort spent to obtain the reward; the magnitude and probability of a reward involved in prompting motivated behaviour are encoded by the dopaminergic (DA) midbrain and its mesoaccumbens DA projections. This type of reward circuity is particularly sensitive to the metabolic state signalled by peripheral mediators, such as insulin or glucagon-like peptide 1 (GLP-1). While in rodents the modulatory effect of metabolic state signals on motivated behaviour is well documented, evidence of state-dependent modulation and the role of incentive motivation underlying overeating in humans is lacking.

The admixture between modern humans and Neandertals has resulted in ∼2% of the genomes of present-day non-Africans being composed of Neandertal DNA. Introgressed Neandertal DNA has been demonstrated to significantly affect the transcriptomic landscape in people today and via this molecular mechanism influence phenotype variation as well. However, little is known about how much of that regulatory impact is mediated through long-range regulatory effects that have been shown to explain ∼20% of expression variation. Here we identified 60 transcription factors (TFs) with their top cis-eQTL SNP in GTEx being of Neandertal ancestry and predicted long-range Neandertal DNA-induced regulatory effects by screening for the predicted target genes of those TFs. We show that the TFs form a significantly connected protein-protein interaction network. Among them are JUN and PRDM5, two brain-expressed TFs that have their predicted target genes enriched in regions devoid of Neandertal DNA. Archaic cis-eQTLs for the 60 TFs include multiple candidates for local adaptation, some of which show significant allele frequency increases over the last ∼10,000 years. A large proportion of the cis-eQTL-associated archaic SNPs have additional associations with various immune traits, schizophrenia, blood cell type composition and anthropometric measures. Finally, we demonstrate that our results are consistent with those of Neandertal DNA-associated empirical trans-eQTLs. Our results suggest that Neandertal DNA significantly influences regulatory networks, that its regulatory reach goes beyond the 40% of genomic sequence it still covers in present-day non-Africans and that via the investigated mechanism Neandertal DNA influences the phenotypic variation in people today.

Survival under selective pressure is driven by the ability of our brain to use sensory information to our advantage to control physiological needs. To that end, neural circuits receive and integrate external environmental cues and internal metabolic signals to form learned sensory associations, consequently motivating and adapting our behaviour. The dopaminergic midbrain plays a crucial role in learning adaptive behaviour and is particularly sensitive to peripheral metabolic signals, including intestinal peptides, such as glucagon-like peptide 1 (GLP-1). In a single-blinded, randomized, controlled, crossover basic human functional magnetic resonance imaging study relying on a computational model of the adaptive learning process underlying behavioural responses, we show that adaptive learning is reduced when metabolic sensing is impaired in obesity, as indexed by reduced insulin sensitivity (participants: N = 30 with normal insulin sensitivity; N = 24 with impaired insulin sensitivity). Treatment with the GLP-1 receptor agonist liraglutide normalizes impaired learning of sensory associations in men and women with obesity. Collectively, our findings reveal that GLP-1 receptor activation modulates associative learning in people with obesity via its central effects within the mesoaccumbens pathway. These findings provide evidence for how metabolic signals can act as neuromodulators to adapt our behaviour to our body's internal state and how GLP-1 receptor agonists work in clinics.

Pernicious anemia is a rare condition characterized by vitamin B12 deficiency anemia due to lack of intrinsic factor, often caused by autoimmune gastritis. Patients with pernicious anemia have a higher incidence of other autoimmune disorders, such as type 1 diabetes, vitiligo, and autoimmune thyroid issues. Therefore, the disease has a clear autoimmune basis, although the genetic susceptibility factors have thus far remained poorly studied. We conduct a genome-wide association study meta-analysis in 2166 cases and 659,516 European controls from population-based biobanks and identify genome-wide significant signals in or near the PTPN22 (rs6679677, p = 1.91 × 10-24, OR = 1.63), PNPT1 (rs12616502, p = 3.14 × 10-8, OR = 1.70), HLA-DQB1 (rs28414666, p = 1.40 × 10-16, OR = 1.38), IL2RA (rs2476491, p = 1.90 × 10-8, OR = 1.22) and AIRE (rs74203920, p = 2.33 × 10-9, OR = 1.83) genes, thus providing robust associations between pernicious anemia and genetic risk factors.