The "antidepressant efficacy" survey (AES) was deployed to > 50,000 23andMe, Inc. research participants to investigate the genetic basis of treatment-resistant depression (TRD) and non-treatment-resistant depression (NTRD). Genome-wide association studies (GWAS) were performed, including TRD vs. NTRD, selective serotonin reuptake inhibitor (SSRI) responders vs. non-responders, serotonin-norepinephrine reuptake inhibitor (SNRI) responders vs. non-responders, and norepinephrine-dopamine reuptake inhibitor responders vs. non-responders. Only the SSRI association reached the genome-wide significance threshold (p < 5 × 10-8): one genomic region in RNF219-AS1 (SNP rs4884091, p = 2.42 × 10-8, OR = 1.21); this association was also observed in the meta-analysis (13,130 responders vs. 6,610 non-responders) of AES and an earlier "antidepressant efficacy and side effects" survey (AESES) cohort. Meta-analysis for SNRI response phenotype derived from AES and AESES (4030 responders vs. 3049 non-responders) identified another genomic region (lead SNP rs4955665, p = 1.62 × 10-9, OR = 1.25) in an intronic region of MECOM passing the genome-wide significance threshold. Meta-analysis for the TRD phenotype (31,068 NTRD vs 5,714 TRD) identified one additional genomic region (lead SNP rs150245813, p = 8.07 × 10-9, OR = 0.80) in 10p11.1 passing the genome-wide significance threshold. A stronger association for rs150245813 was observed in current study (p = 7.35 × 10-7, OR = 0.79) than the previous study (p = 1.40 × 10-3, OR = 0.81), and for rs4955665, a stronger association in previous study (p = 1.21 × 10-6, OR = 1.27) than the current study (p = 2.64 × 10-4, OR = 1.21). In total, three novel loci associated with SSRI or SNRI (responders vs. non-responders), and NTRD vs TRD were identified; gene level association and gene set enrichment analyses implicate enrichment of genes involved in immune process.

A variety of eye-related symptoms due to the overuse of digital devices is collectively referred to as computer vision syndrome (CVS). In this study, a web-based survey about mind and body functions, including eye strain, was conducted on 1998 Japanese volunteers. To investigate the biological mechanisms behind CVS, a multi-trait genome-wide association study (GWAS), a multivariate analysis on individual-level multivariate data, was performed based on the structural equation modeling methodology assuming a causal pathway for a genetic variant to influence each symptom via a single common latent variable. Twelve loci containing lead variants with a suggestive level of significance were identified. Two loci showed relatively strong signals and were associated with TRABD2B relative to the Wnt signaling pathway and SDK1 having neuronal adhesion and immune functions, respectively. By utilizing publicly available eQTL data, colocalization between GWAS and eQTL signals for four loci was detected, and a locus on 2p25.3 showed a strong colocalization (PPH4 > 0.9) on retinal MYT1L, known to play an important role in neuronal differentiation. This study suggested that the use of multivariate questionnaire data and multi-trait GWAS can lead to biologically reasonable findings and enhance our genetic understanding of complex relationships among symptoms related to CVS.

BACKGROUND: Basal cell carcinoma (BCC) of the skin is the most common form of human cancer, with more than 90% of tumours presenting with clear genetic activation of the Hedgehog pathway. However, polygenic risk factors affecting mechanisms such as DNA repair and cell cycle checkpoints or which modulate the tumour microenvironment or host immune system play significant roles in determining whether genetic mutations culminate in BCC development. We set out to define background genetic factors that play a role in influencing BCC susceptibility via promoting or suppressing the effects of oncogenic drivers of BCC. METHODS: We performed genome-wide association studies (GWAS) on 17,416 cases and 375,455 controls. We subsequently performed statistical analysis by integrating data from population-based genetic studies of multi-omics data, including blood- and skin-specific expression quantitative trait loci and methylation quantitative trait loci, thereby defining a list of functionally relevant candidate BCC susceptibility genes from our GWAS loci. We also constructed a local GWAS functional interaction network (consisting of GWAS nearest genes) and another functional interaction network, consisting specifically of candidate BCC susceptibility genes. RESULTS: A total of 71 GWAS loci and 46 functional candidate BCC susceptibility genes were identified. Increased risk of BCC was associated with the decreased expression of 26 susceptibility genes and increased expression of 20 susceptibility genes. Pathway analysis of the functional candidate gene regulatory network revealed strong enrichment for cell cycle, cell death, and immune regulation processes, with a global enrichment of genes and proteins linked to TReg cell biology. CONCLUSIONS: Our genome-wide association analyses and functional interaction network analysis reveal an enrichment of risk variants that function in an immunosuppressive regulatory network, likely hindering cancer immune surveillance and effective antitumour immunity.

BACKGROUND: Age at onset (AAO) is an essential feature of Parkinson's disease (PD) and can help predict disease progression and mortality. Identification of genetic variants influencing AAO of PD could lead to a better understanding of the disease's biological mechanism and provide clinical guidance. However, genetic determinants for AAO of PD remain mostly unknown, especially in the Asian population. OBJECTIVES: To identify genetic determinants for AAO of PD in the Asian population. METHODS: We performed a genome-wide association meta-analysis on AAO of PD in 5166 Chinese patients with PD (Ndiscovery  = 3628, Nreplication  = 1538). We then conducted a further cross-ethnic meta-analysis using our results and summary statistics for the AAO of PD from the European population. RESULTS: The total heritability of AAO of PD was around 0.10 ~ 0.14, similar to that (~0.11) estimated in populations of European ancestry. One novel significant intergenic locus rs9783733 (NDN; PWRN4) was identified (P = 3.14E-09, beta = 2.30, SE = 0.39). Remarkably, this variant could delay AAO of PD by ~2.43 years, with a more considerable effect on males (~3.18 years) than females (~1.45 years). The variant was suggestively significant in the cross-ethnic meta-analysis and suggested a positive selection in the East Asian population. Additionally, cross-ethnic meta-analysis identified a significant locus rs356203 in SNCA (P = 2.35E-11, beta = -0.71, SE = 0.01). CONCLUSIONS: These findings improve the current understanding of the genetic etiology of AAO of PD in different ethnic groups, and provide a new target for further research on PD pathogenesis and potential therapeutic options. © 2021 International Parkinson and Movement Disorder Society.

[Figure: see text].

OBJECTIVE: The aim of this study was to investigate if there is a causal relationship between circulating levels of TIM-1 (T-cell immunoglobulin and mucin domain 1) and incidence of stroke. Approach and Results: Plasma TIM-1 was analyzed in 4591 subjects (40% men; mean age, 57.5 years) attending the Malmö Diet and Cancer Study. Incidence of stroke was studied in relation to TIM-1 levels during a mean of 19.5 years follow-up. Genetic variants associated with TIM-1 (pQTLs [protein quantitative trait loci]) were examined, and a 2-sample Mendelian randomization analysis was performed to explore the role of TIM-1 in stroke using summary statistics from our pQTLs and the MEGASTROKE consortium. A total of 416 stroke events occurred during follow-up, of which 338 were ischemic strokes. After risk factor adjustment, TIM-1 was associated with increased incidence of all-cause stroke (hazards ratio for third versus first tertile, 1.44 [95% CI, 1.10-1.87]; P for trend, 0.004), and ischemic stroke (hazards ratio, 1.42 [95% CI, 1.06-1.90]; P for trend, 0.011). Nineteen independent lead SNPs, located in three genomic risk loci showed significant associations with TIM-1 (P<5×10-8). A 2-sample Mendelian Randomization analysis suggested a causal effect of TIM-1 on stroke (β=0.083, P=0.0004) and ischemic stroke (β=0.102, P=7.7×10-5). CONCLUSIONS: Plasma level of TIM-1 is associated with incidence of stroke. The genetic analyses suggest that this could be a causal relationship.

To identify candidate causal genes of asthma, we performed a genome-wide association study (GWAS) in UK Biobank on a broad asthma definition (n = 56,167 asthma cases and 352,255 controls). We then carried out functional mapping through transcriptome-wide association studies (TWAS) and Mendelian randomization in lung (n = 1,038) and blood (n = 31,684) tissues. The GWAS reveals 72 asthma-associated loci from 116 independent significant variants (PGWAS < 5.0E-8). The most significant lung TWAS gene on 17q12-q21 is GSDMB (PTWAS = 1.42E-54). Other TWAS genes include TSLP on 5q22, RERE on 1p36, CLEC16A on 16p13, and IL4R on 16p12, which all replicated in GTEx lung (n = 515). We demonstrate that the largest fold enrichment of regulatory and functional annotations among asthma-associated variants is in the blood. We map 485 blood eQTL-regulated genes associated with asthma and 50 of them are causal by Mendelian randomization. Prioritization of druggable genes reveals known (IL4R, TSLP, IL6, TNFSF4) and potentially new therapeutic targets for asthma.

So far, 23 germline susceptibility loci have been associated with multiple myeloma (MM) risk. It is unclear whether the genetic variation associated with MM susceptibility also predisposes to its precursor, monoclonal gammopathy of undetermined significance (MGUS). Leveraging 2434 MM cases, 754 MGUS cases, and 2 independent sets of controls (2567/879), we investigated potential shared genetic susceptibility of MM and MGUS by (1) performing MM and MGUS genome-wide association studies (GWAS); (2) validating the association of a polygenic risk score (PRS) based on 23 established MM loci (MM-PRS) with risk of MM, and for the first time with MGUS; and (3) examining genetic correlation of MM and MGUS. Heritability and genetic estimates yielded 17% (standard error [SE] ±0.04) and 15% (SE ±0.11) for MM and MGUS risk, respectively, and a 55% (SE ±0.30) genetic correlation. The MM-PRS was associated with risk of MM when assessed continuously (odds ratio [OR], 1.17 per SD; 95% confidence interval [CI], 1.13-1.21) or categorically (OR, 1.70; 95% CI, 1.38-2.09 for highest; OR, 0.71; 95% CI, 0.55-0.90 for lowest compared with middle quintile). The MM-PRS was similarly associated with MGUS (OR, 1.19 per SD; 95% CI, 1.14-1.26 as a continuous measure, OR, 1.77, 95%CI: 1.29-2.43 for highest and OR, 0.70, 95%CI: 0.50-0.98 for lowest compared with middle quintile). MM and MGUS associations did not differ by age, sex, or MM immunoglobulin isotype. We validated a 23-SNP MM-PRS in an independent series of MM cases and provide evidence for its association with MGUS. Our results suggest shared common genetic susceptibility to MM and MGUS.

BACKGROUND: While evidence indicates that familial predisposition influences the risk of developing degenerative rotator cuff disease (RCD), knowledge of specific genetic markers is limited. We conducted a genome-wide association study of RCD surgery using the UK Biobank, a prospective cohort of 500,000 people (40 to 69 years of age at enrollment) with genotype data. METHODS: Cases with surgery for degenerative RCD were identified using linked hospital records. The cases were defined as an International Classification of Diseases, Tenth Revision (ICD-10) code of M75.1 determined by a trauma/orthopaedic specialist and surgery consistent with RCD treatment. Cases were excluded if a diagnosis of traumatic injury had been made during the same hospital visit. For each case, up to 5 controls matched by age, sex, and follow-up time were chosen from the UK Biobank. Analyses were limited to European-ancestry individuals who were not third-degree or closer relations. We used logistic regression to test for genetic association of 674,405 typed and >10 million imputed markers, after adjusting for age, sex, population principal components, and follow-up. RESULTS: We identified 2,917 RCD surgery cases and 14,158 matched controls. We observed 1 genome-wide significant signal (p < 5 × 10-8) for a novel locus tagged by rs2237352 in the CREB5 gene on chromosome 7 (odds ratio [OR] = 1.17, 95% confidence interval [CI] = 1.11 to 1.24). The single-nucleotide polymorphism (SNP) rs2237352 was imputed with a high degree of confidence (info score = 0.9847) and is common, with a minor allele frequency of 47%. After expanding the control sample to include additional unmatched non-cases, rs2237352 and another SNP in the CREB5 gene, rs12700903, were genome-wide significant. We did not detect genome-wide significant signals at loci associated with RCD in previous studies. CONCLUSIONS: We identified a novel association between a variant in the CREB5 gene and RCD surgery. Validation of this finding in studies with imaging data to confirm diagnoses will be an important next step. CLINICAL RELEVANCE: Identification of genetic RCD susceptibility markers can guide understanding of biological processes in rotator cuff degeneration and help inform disease risk in the clinical setting. LEVEL OF EVIDENCE: Prognostic Level III. See Instructions for Authors for a complete description of levels of evidence.

BACKGROUND: Nocturnal enuresis (bedwetting) is a common disorder affecting 10-16% of 7-year-old children globally. Nocturnal enuresis is highly heritable, but its genetic determinants remain unknown. We aimed to identify genetic variants associated with nocturnal enuresis and explore its genetic architecture and underlying biology. METHODS: We did a genome-wide association study (GWAS) of nocturnal enuresis. Nocturnal enuresis cases were identified in iPSYCH2012, a large Danish population-based case cohort established to investigate mental disorders, on the basis of 10th revision of the International Statistical Classification of Diseases (ICD-10) diagnoses and redeemed desmopressin prescriptions in Danish registers. The GWAS was done in a genetically homogeneous sample of unrelated individuals using logistic regression with relevant covariates. All genome-wide significant variants were analysed for their association with nocturnal enuresis in an independent Icelandic sample from deCODE genetics. Standardised polygenic risk scores for attention-deficit hyperactivity disorder (ADHD) and autism spectrum disorder were constructed from summary statistics of large GWASs and analysed for association with nocturnal enuresis. FINDINGS: The GWAS included 3882 nocturnal enuresis cases and 31 073 controls. We found two loci at chromosome 6 and chromosome 13 significantly associated with nocturnal enuresis. Six genetic variants at the two loci (five variants at chromosome 6q16.2 and one variant at chromosome 13q22.3) surpassed the threshold for genome-wide significance (p<5 × 10-8). There were two lead variants: rs9376454 (chromosome 6q16.2), with an odds ratio (OR) of 1·199 (95% CI 1·135-1·267; p=9·91 × 10-11), and rs60721117 (chromosome 13q22.3), with an OR of 1·149 (1·095-1·205; p=1·21 × 10-8). All associated variants in the chromosome 6 locus were replicated (p<8 × 10-3) in the independent Icelandic cohort of 5475 nocturnal enuresis cases and 303 996 controls, whereas the associated variant in the chromosome 13 locus showed nominal significant association (p=0·031). The percentage of nocturnal enuresis phenotypic variance explained by the common genetic variants was 23·9-30·4%. Polygenic risk for ADHD was associated with nocturnal enuresis (OR 1·06, 95% CI, 1·01-1·10; p=0·011). Among the potential nocturnal enuresis risk genes mapped, PRDM13 and EDNRB have biological functions associated with known pathophysiological mechanisms in nocturnal enuresis, and SIM1 regulates the formation of the hypothalamic neuroendocrine lineage that produces arginine vasopressin, a well known nocturnal enuresis drug target. INTERPRETATION: This study shows that common genetic variants contribute considerably to nocturnal enuresis, and it identifies potential nocturnal enuresis risk genes with roles in sleep, urine production, and bladder function. Given that available treatments target these mechanisms, any of the identified genes and their functional gene networks are potential drug targets. FUNDING: The Lundbeck Foundation Initiative for Integrative Psychiatric Research (iPSYCH), Stanley Foundation.