Genetic drivers and cellular selection of female mosaic X chromosome loss.

Mosaic loss of the X chromosome (mLOX) is the most common clonal somatic alteration in leukocytes of female individuals1,2, but little is known about its genetic determinants or phenotypic consequences. Here, to address this, we used data from 883,574 female participants across 8 biobanks; 12% of participants exhibited detectable mLOX in approximately 2% of leukocytes. Female participants with mLOX had an increased risk of myeloid and lymphoid leukaemias. Genetic analyses identified 56 common variants associated with mLOX, implicating genes with roles in chromosomal missegregation, cancer predisposition and autoimmune diseases. Exome-sequence analyses identified rare missense variants in FBXO10 that confer a twofold increased risk of mLOX. Only a small fraction of associations was shared with mosaic Y chromosome loss, suggesting that distinct biological processes drive formation and clonal expansion of sex chromosome missegregation. Allelic shift analyses identified X chromosome alleles that are preferentially retained in mLOX, demonstrating variation at many loci under cellular selection. A polygenic score including 44 allelic shift loci correctly inferred the retained X chromosomes in 80.7% of mLOX cases in the top decile. Our results support a model in which germline variants predispose female individuals to acquiring mLOX, with the allelic content of the X chromosome possibly shaping the magnitude of clonal expansion.

Large registry-based analysis of genetic predisposition to tuberculosis identifies genetic risk factors at HLA.

Tuberculosis is a significant public health concern resulting in the death of over 1 million individuals each year worldwide. While treatment options and vaccines exist, a substantial number of infections still remain untreated or are caused by treatment resistant strains. Therefore, it is important to identify mechanisms that contribute to risk and prognosis of tuberculosis as this may provide tools to understand disease mechanisms and provide novel treatment options for those with severe infection. Our goal was to identify genetic risk factors that contribute to the risk of tuberculosis and to understand biological mechanisms and causality behind the risk of tuberculosis. A total of 1895 individuals in the FinnGen study had International Classification of Diseases-based tuberculosis diagnosis. Genome-wide association study analysis identified genetic variants with statistically significant association with tuberculosis at the human leukocyte antigen (HLA) region (P < 5e-8). Fine mapping of the HLA association provided evidence for one protective haplotype tagged by HLA DQB1*05:01 (P = 1.82E-06, OR = 0.81 [CI 95% 0.74-0.88]), and predisposing alleles tagged by HLA DRB1*13:02 (P = 0.00011, OR = 1.35 [CI 95% 1.16-1.57]). Furthermore, genetic correlation analysis showed association with earlier reported risk factors including smoking (P < 0.05). Mendelian randomization supported smoking as a risk factor for tuberculosis (inverse-variance weighted P < 0.05, OR = 1.83 [CI 95% 1.15-2.93]) with no significant evidence of pleiotropy. Our findings indicate that specific HLA alleles associate with the risk of tuberculosis. In addition, lifestyle risk factors such as smoking contribute to the risk of developing tuberculosis.

The impact of Mendelian sleep and circadian genetic variants in a population setting.

Rare variants in ten genes have been reported to cause Mendelian sleep conditions characterised by extreme sleep duration or timing. These include familial natural short sleep (ADRB1, DEC2/BHLHE41, GRM1 and NPSR1), advanced sleep phase (PER2, PER3, CRY2, CSNK1D and TIMELESS) and delayed sleep phase (CRY1). The association of variants in these genes with extreme sleep conditions were usually based on clinically ascertained families, and their effects when identified in the population are unknown. We aimed to determine the effects of these variants on sleep traits in large population-based cohorts. We performed genetic association analysis of variants previously reported to be causal for Mendelian sleep and circadian conditions. Analyses were performed using 191,929 individuals with data on sleep and whole-exome or genome-sequence data from 4 population-based studies: UK Biobank, FINRISK, Health-2000-2001, and the Multi-Ethnic Study of Atherosclerosis (MESA). We identified sleep disorders from self-report, hospital and primary care data. We estimated sleep duration and timing measures from self-report and accelerometery data. We identified carriers for 10 out of 12 previously reported pathogenic variants for 8 of the 10 genes. They ranged in frequency from 1 individual with the variant in CSNK1D to 1,574 individuals with a reported variant in the PER3 gene in the UK Biobank. No carriers for variants reported in NPSR1 or PER2 were identified. We found no association between variants analyzed and extreme sleep or circadian phenotypes. Using sleep timing as a proxy measure for sleep phase, only PER3 and CRY1 variants demonstrated association with earlier and later sleep timing, respectively; however, the magnitude of effect was smaller than previously reported (sleep midpoint ~7 mins earlier and ~5 mins later, respectively). We also performed burden tests of protein truncating (PTVs) or rare missense variants for the 10 genes. Only PTVs in PER2 and PER3 were associated with a relevant trait (for example, 64 individuals with a PTV in PER2 had an odds ratio of 4.4 for being "definitely a morning person", P = 4x10-8; and had a 57-minute earlier midpoint sleep, P = 5x10-7). Our results indicate that previously reported variants for Mendelian sleep and circadian conditions are often not highly penetrant when ascertained incidentally from the general population.

Genome-wide association studies highlight novel risk loci for septal defects and left-sided congenital heart defects.

BACKGROUND: Congenital heart defects (CHD) are structural defects of the heart affecting approximately 1% of newborns. They exhibit low penetrance and non-Mendelian patterns of inheritance as varied and complex traits. While genetic factors are known to play an important role in the development of CHD, the specific genetics remain unknown for the majority of patients. To elucidate the underlying genetic risk, we performed a genome wide association study (GWAS) of CHDs in general and specific CHD subgroups using the FinnGen Release 10 (R10) (N > 393,000), followed by functional fine-mapping through eQTL and co-localization analyses using the GTEx database. RESULTS: We discovered three genome-wide significant loci associated with general CHD. Two of them were located in chromosome 17: 17q21.32 (rs2316327, intronic: LRRC37A2, Odds ratio (OR) [95% Confidence Interval (CI)] = 1.17[1.12-1.23], p = 1.5 × 10-9) and 17q25.3 (rs1293973611, nearest: BAHCC1, OR[95%CI] = 4.48[2.80-7.17], p = 7.0 × 10-10), respectively, and in addition to general CHD, the rs1293973611 locus was associated with the septal defect subtype. The third locus was in band 1p21.2 (rs35046143, nearest: PALMD, OR[95%CI] = 1.15[1.09-1.21], p = 7.1 × 10-9), and it was associated with general CHD and left-sided lesions. In the subgroup analysis, two additional loci were associated with septal defects (rs75230966 and rs6824295), and one with left-sided lesions (rs1305393195). In the eQTL analysis the variants rs2316327 (general CHD), and rs75230966 (septal defects) both located in 17q21.32 (with a LD r2 of 0.41) were both predicted to significantly associate with the expression of WNT9B in the atrial appendage tissue category. This effect was further confirmed by co-localization analysis, which also implicated WNT3 expression in the atrial appendage. A meta-analysis of general CHD together with the UK Biobank (combined N = 881,678) provided a different genome-wide significant locus in LRRC37A2; rs16941382 (OR[95%CI] = 1.15[1.11-1.20], p = 1.5 × 10-9) which is in significant LD with rs2316327. CONCLUSIONS: Our results of general CHD and different CHD subcategories identified a complex risk locus on chromosome 17 near BAHCC1 and LRRC37A2, interacting with the genes WNT9B, WNT3 and MYL4, may constitute potential novel CHD risk associated loci, warranting future experimental tests to determine their role.

Kleine-Levin syndrome is associated with birth difficulties and genetic variants in the TRANK1 gene loci.

Kleine-Levin syndrome (KLS) is a rare disorder characterized by severe episodic hypersomnia, with cognitive impairment accompanied by apathy or disinhibition. Pathophysiology is unknown, although imaging studies indicate decreased activity in hypothalamic/thalamic areas during episodes. Familial occurrence is increased, and risk is associated with reports of a difficult birth. We conducted a worldwide case-control genome-wide association study in 673 KLS cases collected over 14 y, and ethnically matched 15,341 control individuals. We found a strong genome-wide significant association (rs71947865, Odds Ratio [OR] = 1.48, P = 8.6 × 10-9) within the 3'region of TRANK1 gene locus, previously associated with bipolar disorder and schizophrenia. Strikingly, KLS cases with rs71947865 variant had significantly increased reports of a difficult birth. As perinatal outcomes have dramatically improved over the last 40 y, we further stratified our sample by birth years and found that recent cases had a significantly reduced rs71947865 association. While the rs71947865 association did not replicate in the entire follow-up sample of 171 KLS cases, rs71947865 was significantly associated with KLS in the subset follow-up sample of 59 KLS cases who reported birth difficulties (OR = 1.54, P = 0.01). Genetic liability of KLS as explained by polygenic risk scores was increased (pseudo R2 = 0.15; P < 2.0 × 10-22 at P = 0.5 threshold) in the follow-up sample. Pathway analysis of genetic associations identified enrichment of circadian regulation pathway genes in KLS cases. Our results suggest links between KLS, circadian regulation, and bipolar disorder, and indicate that the TRANK1 polymorphisms in conjunction with reported birth difficulties may predispose to KLS.

Association of accelerometer-derived sleep measures with lifetime psychiatric diagnoses: A cross-sectional study of 89,205 participants from the UK Biobank.

BACKGROUND: Sleep problems are both symptoms of and modifiable risk factors for many psychiatric disorders. Wrist-worn accelerometers enable objective measurement of sleep at scale. Here, we aimed to examine the association of accelerometer-derived sleep measures with psychiatric diagnoses and polygenic risk scores in a large community-based cohort. METHODS AND FINDINGS: In this post hoc cross-sectional analysis of the UK Biobank cohort, 10 interpretable sleep measures-bedtime, wake-up time, sleep duration, wake after sleep onset, sleep efficiency, number of awakenings, duration of longest sleep bout, number of naps, and variability in bedtime and sleep duration-were derived from 7-day accelerometry recordings across 89,205 participants (aged 43 to 79, 56% female, 97% self-reported white) taken between 2013 and 2015. These measures were examined for association with lifetime inpatient diagnoses of major depressive disorder, anxiety disorders, bipolar disorder/mania, and schizophrenia spectrum disorders from any time before the date of accelerometry, as well as polygenic risk scores for major depression, bipolar disorder, and schizophrenia. Covariates consisted of age and season at the time of the accelerometry recording, sex, Townsend deprivation index (an indicator of socioeconomic status), and the top 10 genotype principal components. We found that sleep pattern differences were ubiquitous across diagnoses: each diagnosis was associated with a median of 8.5 of the 10 accelerometer-derived sleep measures, with measures of sleep quality (for instance, sleep efficiency) generally more affected than mere sleep duration. Effect sizes were generally small: for instance, the largest magnitude effect size across the 4 diagnoses was ß = -0.11 (95% confidence interval -0.13 to -0.10, p = 3 × 10-56, FDR = 6 × 10-55) for the association between lifetime inpatient major depressive disorder diagnosis and sleep efficiency. Associations largely replicated across ancestries and sexes, and accelerometry-derived measures were concordant with self-reported sleep properties. Limitations include the use of accelerometer-based sleep measurement and the time lag between psychiatric diagnoses and accelerometry. CONCLUSIONS: In this study, we observed that sleep pattern differences are a transdiagnostic feature of individuals with lifetime mental illness, suggesting that they should be considered regardless of diagnosis. Accelerometry provides a scalable way to objectively measure sleep properties in psychiatric clinical research and practice, even across tens of thousands of individuals.

Cross-disorder analysis of schizophrenia and 19 immune-mediated diseases identifies shared genetic risk.

Many immune diseases occur at different rates among people with schizophrenia compared to the general population. Here, we evaluated whether this phenomenon might be explained by shared genetic risk factors. We used data from large genome-wide association studies to compare the genetic architecture of schizophrenia to 19 immune diseases. First, we evaluated the association with schizophrenia of 581 variants previously reported to be associated with immune diseases at genome-wide significance. We identified five variants with potentially pleiotropic effects. While colocalization analyses were inconclusive, functional characterization of these variants provided the strongest evidence for a model in which genetic variation at rs1734907 modulates risk of schizophrenia and Crohn's disease via altered methylation and expression of EPHB4-a gene whose protein product guides the migration of neuronal axons in the brain and the migration of lymphocytes towards infected cells in the immune system. Next, we investigated genome-wide sharing of common variants between schizophrenia and immune diseases using cross-trait LD score regression. Of the 11 immune diseases with available genome-wide summary statistics, we observed genetic correlation between six immune diseases and schizophrenia: inflammatory bowel disease (rg = 0.12 ± 0.03, P = 2.49 × 10-4), Crohn's disease (rg = 0.097 ± 0.06, P = 3.27 × 10-3), ulcerative colitis (rg = 0.11 ± 0.04, P = 4.05 × 10-3), primary biliary cirrhosis (rg = 0.13 ± 0.05, P = 3.98 × 10-3), psoriasis (rg = 0.18 ± 0.07, P = 7.78 × 10-3) and systemic lupus erythematosus (rg = 0.13 ± 0.05, P = 3.76 × 10-3). With the exception of ulcerative colitis, the degree and direction of these genetic correlations were consistent with the expected phenotypic correlation based on epidemiological data. Our findings suggest shared genetic risk factors contribute to the epidemiological association of certain immune diseases and schizophrenia.

Genetic analysis of obstructive sleep apnoea discovers a strong association with cardiometabolic health.

There is currently limited understanding of the genetic aetiology of obstructive sleep apnoea (OSA). We aimed to identify genetic loci associated with OSA risk, and to test if OSA and its comorbidities share a common genetic background.We conducted the first large-scale genome-wide association study of OSA using the FinnGen study (217 955 individuals) with 16 761 OSA patients identified using nationwide health registries.We estimated 0.08 (95% CI 0.06-0.11) heritability and identified five loci associated with OSA (p<5.0×10-8): rs4837016 near GAPVD1 (GTPase activating protein and VPS9 domains 1), rs10928560 near CXCR4 (C-X-C motif chemokine receptor type 4), rs185932673 near CAMK1D (calcium/calmodulin-dependent protein kinase ID) and rs9937053 near FTO (fat mass and obesity-associated protein; a variant previously associated with body mass index (BMI)). In a BMI-adjusted analysis, an association was observed for rs10507084 near RMST/NEDD1 (rhabdomyosarcoma 2 associated transcript/NEDD1 γ-tubulin ring complex targeting factor). We found high genetic correlations between OSA and BMI (rg=0.72 (95% CI 0.62-0.83)), and with comorbidities including hypertension, type 2 diabetes, coronary heart disease, stroke, depression, hypothyroidism, asthma and inflammatory rheumatic disease (rg>0.30). The polygenic risk score for BMI showed 1.98-fold increased OSA risk between the highest and the lowest quintile, and Mendelian randomisation supported a causal relationship between BMI and OSA.Our findings support the causal link between obesity and OSA, and the joint genetic basis between OSA and comorbidities.

HLA-DPB1 and HLA class I confer risk of and protection from narcolepsy.

Type 1 narcolepsy, a disorder caused by a lack of hypocretin (orexin), is so strongly associated with human leukocyte antigen (HLA) class II HLA-DQA1(∗)01:02-DQB1(∗)06:02 (DQ0602) that very few non-DQ0602 cases have been reported. A known triggering factor for narcolepsy is pandemic 2009 influenza H1N1, suggesting autoimmunity triggered by upper-airway infections. Additional effects of other HLA-DQ alleles have been reported consistently across multiple ethnic groups. Using over 3,000 case and 10,000 control individuals of European and Chinese background, we examined the effects of other HLA loci. After careful matching of HLA-DR and HLA-DQ in case and control individuals, we found strong protective effects of HLA-DPA1(∗)01:03-DPB1(∗)04:02 (DP0402; odds ratio [OR] = 0.51 [0.38-0.67], p = 1.01 × 10(-6)) and HLA-DPA1(∗)01:03-DPB1(∗)04:01 (DP0401; OR = 0.61 [0.47-0.80], p = 2.07 × 10(-4)) and predisposing effects of HLA-DPB1(∗)05:01 in Asians (OR = 1.76 [1.34-2.31], p = 4.71 × 10(-05)). Similar effects were found by conditional analysis controlling for HLA-DR and HLA-DQ with DP0402 (OR = 0.45 [0.38-0.55] p = 8.99 × 10(-17)) and DP0501 (OR = 1.38 [1.18-1.61], p = 7.11 × 10(-5)). HLA-class-II-independent associations with HLA-A(∗)11:01 (OR = 1.32 [1.13-1.54], p = 4.92 × 10(-4)), HLA-B(∗)35:03 (OR = 1.96 [1.41-2.70], p = 5.14 × 10(-5)), and HLA-B(∗)51:01 (OR = 1.49 [1.25-1.78], p = 1.09 × 10(-5)) were also seen across ethnic groups in the HLA class I region. These effects might reflect modulation of autoimmunity or indirect effects of HLA class I and HLA-DP alleles on response to viral infections such as that of influenza.

Genetic variants in the HLA class II region associated with risk of cutaneous squamous cell carcinoma.

BACKGROUND: The immune system has been implicated in the pathophysiology of cutaneous squamous cell carcinoma (cSCC) as evidenced by the substantially increased risk of cSCC in immunosuppressed individuals. Associations between cSCC risk and single nucleotide polymorphisms (SNPs) in the HLA region have been identified by genome-wide association studies (GWAS). The translation of the associated HLA SNPs to structural amino acids changes in HLA molecules has not been previously elucidated. METHODS: Using data from a GWAS that included 7238 cSCC cases and 56,961 controls of non-Hispanic white ancestry, we imputed classical alleles and corresponding amino acid changes in HLA genes. Logistic regression models were used to examine associations between cSCC risk and genotyped or imputed SNPs, classical HLA alleles, and amino acid changes. RESULTS: Among the genotyped SNPs, cSCC risk was associated with rs28535317 (OR = 1.20, p = 9.88 × 10- 11) corresponding to an amino-acid change from phenylalanine to leucine at codon 26 of HLA-DRB1 (OR = 1.17, p = 2.48 × 10- 10). An additional independent association was observed for a threonine to isoleucine change at codon 107 of HLA-DQA1 (OR = 1.14, p = 2.34 × 10- 9). Among the classical HLA alleles, cSCC was associated with DRB1*01 (OR = 1.18, p = 5.86 × 10- 10). Conditional analyses revealed additional independent cSCC associations with DQA1*05:01 and DQA1*05:05. Extended haplotype analysis was used to complement the imputed haplotypes, which identified three extended haplotypes in the HLA-DR and HLA-DQ regions. CONCLUSIONS: Associations with specific HLA-DR and -DQ alleles are likely to explain previously observed GWAS signals in the HLA region associated with cSCC risk.

Autoimmunity in narcolepsy.

PURPOSE OF REVIEW: Summarize the recent findings in narcolepsy focusing on the environmental and genetic risk factors in disease development. RECENT FINDINGS: Both genetic and epidemiological evidence point towards an autoimmune mechanism in the destruction of orexin/hypocretin neurons. Recent studies suggest both humoral and cellular immune responses in the disease development. SUMMARY: Narcolepsy is a severe sleep disorder, in which neurons producing orexin/hypocretin in the hypothalamus are destroyed. The core symptoms of narcolepsy are debilitating, extreme sleepiness, cataplexy, and abnormalities in the structure of sleep. Both genetic and epidemiological evidence point towards an autoimmune mechanism in the destruction of orexin/hypocretin neurons. Importantly, the highest environmental risk is seen with influenza-A infection and immunization. However, how the cells are destroyed is currently unknown. In this review we summarize the disease symptoms, and focus on the immunological findings in narcolepsy. We also discuss the environmental and genetic risk factors as well as propose a model for disease development.

Genome wide analysis of narcolepsy in China implicates novel immune loci and reveals changes in association prior to versus after the 2009 H1N1 influenza pandemic.

Previous studies in narcolepsy, an autoimmune disorder affecting hypocretin (orexin) neurons and recently associated with H1N1 influenza, have demonstrated significant associations with five loci. Using a well-characterized Chinese cohort, we refined known associations in TRA@ and P2RY11-DNMT1 and identified new associations in the TCR beta (TRB@; rs9648789 max P = 3.7 × 10(-9) OR 0.77), ZNF365 (rs10995245 max P = 1.2 × 10(-11) OR 1.23), and IL10RB-IFNAR1 loci (rs2252931 max P = 2.2 × 10(-9) OR 0.75). Variants in the Human Leukocyte Antigen (HLA)- DQ region were associated with age of onset (rs7744020 P = 7.9×10(-9) beta -1.9 years) and varied significantly among cases with onset after the 2009 H1N1 influenza pandemic compared to previous years (rs9271117 P = 7.8 × 10(-10) OR 0.57). These reflected an association of DQB1*03:01 with earlier onset and decreased DQB1*06:02 homozygosity following 2009. Our results illustrate how genetic association can change in the presence of new environmental challenges and suggest that the monitoring of genetic architecture over time may help reveal the appearance of novel triggers for autoimmune diseases.

Genome-wide scan of job-related exhaustion with three replication studies implicate a susceptibility variant at the UST gene locus.

Job-related exhaustion is the core dimension of burnout, a work-related stress syndrome that has several negative health consequences. In this study, we explored the molecular genetic background of job-related exhaustion. A genome-wide analysis of job-related exhaustion was performed in the GENMETS subcohort (n = 1256) of the Finnish population-based Health 2000 study. Replication analyses included an analysis of the strongest associations in the rest of the Health 2000 sample (n = 1660 workers) and in three independent populations (the FINRISK population cohort, n = 10 753; two occupational cohorts, total n = 1451). Job-related exhaustion was ascertained using a standard self-administered questionnaire (the Maslach Burnout Inventory (MBI)-GS exhaustion scale in the Health 2000 sample and the occupational cohorts) or a single question (FINRISK). A variant located in an intron of UST, uronyl-2-sulfotransferase (rs13219957), gave the strongest statistical evidence in the initial genome-wide study (P = 1.55 × 10(-7)), and was associated with job-related exhaustion in all the replication sets (P < 0.05; P = 6.75 × 10(-7) from the meta-analysis). Consistent with studies of mood disorders, individual common genetic variants did not have any strong effect on job-related exhaustion. However, the nominally significant signals from the allelic variant of UST in four separate samples suggest that this variant might be a weak risk factor for job-related exhaustion. Together with the previously reported associations of other dermatan/chondroitin sulfate genes with mood disorders, these results indicate a potential molecular pathway for stress-related traits and mark a candidate region for further studies of job-related and general exhaustion.

Comparison of Pandemrix and Arepanrix, two pH1N1 AS03-adjuvanted vaccines differentially associated with narcolepsy development.

Narcolepsy onset in children has been associated with the 2009 influenza A H1N1 pandemic and vaccination with Pandemrix. However it was not clearly observed with other adjuvanted pH1N1 vaccines such as Arepanrix or Focetria. Our aim was to characterize the differences between Pandemrix and Arepanrix that might explain the risk for narcolepsy after Pandemrix vaccination using 2D-DIGE and mass spectrometry (MS). We found that Pandemrix (2009 batch) and Arepanrix (2010 batch) showed 5 main viral proteins: hemagglutinin HA1 and HA2 subunits, neuraminidase NA, nucleoprotein NP, and matrix protein MA1 and non-viral proteins from the Gallus gallus growth matrix used in the manufacturing of the vaccines. Latticed patterns of HA1, HA2 and NA indicated charge and molecular weight heterogeneity, a phenomenon likely caused by glycosylation and sulfation. Overall, Pandemrix contained more NP and NA, while Arepanrix displayed a larger diversity of viral and chicken proteins, with the exception of five chicken proteins (PDCD6IP, TSPAN8, H-FABP, HSP and TUB proteins) that were relatively more abundant in Pandemrix. Glycosylation patterns were similar in both vaccines. A higher degree of deamidation and dioxidation was found in Pandemrix, probably reflecting differential degradation across batches. Interestingly, HA1 146N (residue 129N in the mature protein) displayed a 10-fold higher deamidation in Arepanrix versus Pandemrix. In recent vaccine strains and Focetria, 146N is mutated to D which is associated with increased production yields suggesting that 146N deamidation may have also occurred during the manufacturing of Arepanrix. The presence of 146N in large relative amounts in Pandemrix and the wild type virus and in lower relative quantities in Arepanrix or other H1N1 vaccines may have affected predisposition to narcolepsy.

Genome-wide association study of sleep duration in the Finnish population.

Sleep duration is genetically regulated, but the genetic variants are largely unknown. We aimed to identify such genes using a genome-wide association study (GWAS) combined with RNA expression at the population level, and with experimental verification. A GWAS was performed in a Finnish sample (n = 1941), and variants with suggestive association (P < 5 × 10(-5) ) were tested in a follow-up sample from the same population with sleep duration (n = 6834) and time in bed (n = 1720). Variants with pointwise association of P < 0.05 in the follow-up sample were analysed further. First, we correlated genotypes with transcript expression levels with sleep duration (n = 207). The expression levels of significant transcripts were further studied in experimental sleep restriction. Of the 31 variants with P < 5 × 10(-5) in the discovery sample, three variants showed nominal allelic association (P < 0.05) in the follow-up sample: rs10914351, near PTPRU (P = 0.049), rs1037079 in PCDH7-CENTD1 (P = 0.011) and rs2031573 near KLF6 (P = 0.044). The risk alleles for shorter sleep (rs2031573 and rs1037079) were also associated with higher KLF6 and PCDH7 expression levels (P < 0.05). Experimental sleep restriction increased the expression of KLF6 (P < 0.01). These data suggest that rs2031573 near KLF6 or related loci and rs1037079 between PCDH7-CENTD1 or related loci may contribute to the regulation of sleep duration via gene expression. These results illustrate the utility of combining different analytical approaches to identify genetic determinants for traits related to sleep physiology. However, additional studies are needed in order to understand the roles of KLF6 and PCDH7 in sleep regulation.

Genetics of sleep medication purchases suggests causality from sleep problems to psychiatric traits.

STUDY OBJECTIVES: Over 10% of the population in Europe and in the United States use sleep medication to manage sleep problems. Our objective was to elucidate genetic risk factors and clinical correlates that contribute to sleep medication purchase and estimate the comorbid impact of sleep problems. METHODS: We performed epidemiological analysis for psychiatric diagnoses, and genetic association studies of sleep medication purchase in 797 714 individuals from FinnGen Release 7 (N = 311 892) and from the UK Biobank (N = 485 822). Post-association analyses included genetic correlation, co-localization, Mendelian randomization (MR), and polygenic risk estimation. RESULTS: In a GWAS we identified 27 genetic loci significantly associated with sleep medication, located in genes associated with sleep; AUTS2, CACNA1C, MEIS1, KIRREL3, PAX8, GABRA2, psychiatric traits; CACNA1C, HIST1H2BD, NUDT12. TOPAZ1 and TSNARE1. Co-localization and expression analysis emphasized effects on the KPNA2, GABRA2, and CACNA1C expression in the brain. Sleep medications use was epidemiologically related to psychiatric traits in FinnGen (OR [95% (CI)] = 3.86 [3.78 to 3.94], p < 2 × 10-16), and the association was accentuated by genetic correlation and MR; depression (rg = 0.55 (0.027), p = 2.86 × 10-89, p MR = 4.5 × 10-5), schizophrenia (rg = 0.25 (0.026), p = 2.52 × 10-21, p MR = 2 × 10-4), and anxiety (rg = 0.44 (0.047), p = 2.88 × 10-27, p MR = 8.6 × 10-12). CONCLUSIONS: These results demonstrate the genetics behind sleep problems and the association between sleep problems and psychiatric traits. Our results highlight the scientific basis for sleep management in treating the impact of psychiatric diseases.

Nightmares share genetic risk factors with sleep and psychiatric traits.

Nightmares are vivid, extended, and emotionally negative or negative dreams that awaken the dreamer. While sporadic nightmares and bad dreams are common and generally harmless, frequent nightmares often reflect underlying pathologies of emotional regulation. Indeed, insomnia, depression, anxiety, or alcohol use have been associated with nightmares in epidemiological and clinical studies. However, the connection between nightmares and their comorbidities are poorly understood. Our goal was to examine the genetic risk factors for nightmares and estimate correlation or causality between nightmares and comorbidities. We performed a genome-wide association study (GWAS) in 45,255 individuals using a questionnaire-based assessment on the frequency of nightmares during the past month and genome-wide genotyping data. While the GWAS did not reveal individual risk variants, heritability was estimated at 5%. In addition, the genetic correlation analysis showed a robust correlation (rg > 0.4) of nightmares with anxiety (rg = 0.671, p = 7.507e-06), depressive (rg = 0.562, p = 1.282e-07) and posttraumatic stress disorders (rg = 0.4083, p = 0.0152), and personality trait neuroticism (rg = 0.667, p = 4.516e-07). Furthermore, Mendelian randomization suggested causality from insomnia to nightmares (beta = 0.027, p = 0.0002). Our findings suggest that nightmares share genetic background with psychiatric traits and that insomnia may increase an individual's liability to experience frequent nightmares. Given the significant correlations with psychiatric and psychological traits, it is essential to grow awareness of how nightmares affect health and disease and systematically collect information about nightmares, especially from clinical samples and larger cohorts.

Genetic Analysis of Obstructive Sleep Apnea and Its Relationship with Severe COVID-19.

Rationale: Although patients with obstructive sleep apnea (OSA) have a higher risk for coronavirus disease (COVID-19) hospitalization, the causal relationship has remained unexplored. Objectives: To understand the causal relationship between OSA and COVID-19 by leveraging data from vaccination and electronic health records, genetic risk factors from genome-wide association studies, and Mendelian randomization. Methods: We elucidated genetic risk factors for OSA using FinnGen (total N = 377,277), performing genome-wide association. We used the associated variants as instruments for univariate and multivariate Mendelian randomization (MR) analyses and computed absolute risk reduction against COVID-19 hospitalization with or without vaccination. Results: We identified nine novel loci for OSA and replicated our findings in the Million Veteran Program. Furthermore, MR analysis showed that OSA was a causal risk factor for severe COVID-19 (P = 9.41 × 10-4). Probabilistic modeling showed that the strongest genetic risk factor for OSA at the FTO locus reflected a signal of higher body mass index (BMI), whereas BMI-independent association was seen with the earlier reported SLC9A4 locus and a MECOM locus, which is a transcriptional regulator with 210-fold enrichment in the Finnish population. Similarly, multivariate MR analysis showed that the causality for severe COVID-19 was driven by BMI (multivariate MR P = 5.97 × 10-6, ß = 0.47). Finally, vaccination reduced the risk for COVID-19 hospitalization more in the patients with OSA than in the non-OSA controls, with respective absolute risk reductions of 13.3% versus 6.3%. Conclusions: Our analysis identified novel genetic risk factors for OSA and showed that OSA is a causal risk factor for severe COVID-19. The effect is predominantly explained by higher BMI and suggests BMI-dependent effects at the level of individual variants and at the level of comorbid causality.

SCGB1D2 inhibits growth of Borrelia burgdorferi and affects susceptibility to Lyme disease.

Lyme disease is a tick-borne disease caused by bacteria of the genus Borrelia. The host factors that modulate susceptibility for Lyme disease have remained mostly unknown. Using epidemiological and genetic data from FinnGen and Estonian Biobank, we identify two previously known variants and an unknown common missense variant at the gene encoding for Secretoglobin family 1D member 2 (SCGB1D2) protein that increases the susceptibility for Lyme disease. Using live Borrelia burgdorferi (Bb) we find that recombinant reference SCGB1D2 protein inhibits the growth of Bb in vitro more efficiently than the recombinant protein with SCGB1D2 P53L deleterious missense variant. Finally, using an in vivo murine infection model we show that recombinant SCGB1D2 prevents infection by Borrelia in vivo. Together, these data suggest that SCGB1D2 is a host defense factor present in the skin, sweat, and other secretions which protects against Bb infection and opens an exciting therapeutic avenue for Lyme disease.

Proteogenomic analyses identify coagulation factor XI as a thromboinflammatory mediator of long COVID.

SARS-CoV-2 infection can result in long COVID, characterized by post-acute symptoms from multiple organ systems. Current hypotheses on mechanisms underlying long COVID include persistent inflammation and dysregulated coagulation; however, precise mechanisms and causal mediators remain unclear. Here, we tested the associations of genetic instruments for 49 complement and coagulation factors from the UK Biobank ( N =34,557) with long COVID in the Long COVID Host Genetics Initiative ( N =997,600). Primary analyses revealed that genetically predicted higher factor XI increased long COVID risk (odds ratio, 1.17 [95% confidence interval, 1.08-1.27] per standard deviation; P =1.7×10 -4 ). This association was robust to sensitivity analyses using pleiotropy-robust methods and different genetic instruments and was replicated using proteogenomic data from an Icelandic cohort. Genetically predicted factor XI was also associated with venous thromboembolism, but not with acute COVID-19 or long COVID-resembling conditions. Collectively, these findings provide genetic evidence implicating factor XI in the biology of long COVID.