Fundamentals for predicting transcriptional regulations from DNA sequence patterns.

Cell-type-specific regulatory elements, cataloged through extensive experiments and bioinformatics in large-scale consortiums, have enabled enrichment analyses of genetic associations that primarily utilize positional information of the regulatory elements. These analyses have identified cell types and pathways genetically associated with human complex traits. However, our understanding of detailed allelic effects on these elements' activities and on-off states remains incomplete, hampering the interpretation of human genetic study results. This review introduces machine learning methods to learn sequence-dependent transcriptional regulation mechanisms from DNA sequences for predicting such allelic effects (not associations). We provide a concise history of machine-learning-based approaches, the requirements, and the key computational processes, focusing on primers in machine learning. Convolution and self-attention, pivotal in modern deep-learning models, are explained through geometrical interpretations using dot products. This facilitates understanding of the concept and why these have been used for machine learning for DNA sequences. These will inspire further research in this genetics and genomics field.

Annotation of nuclear lncRNAs based on chromatin interactions.

The human genome is pervasively transcribed and produces a wide variety of long non-coding RNAs (lncRNAs), constituting the majority of transcripts across human cell types. Some specific nuclear lncRNAs have been shown to be important regulatory components acting locally. As RNA-chromatin interaction and Hi-C chromatin conformation data showed that chromatin interactions of nuclear lncRNAs are determined by the local chromatin 3D conformation, we used Hi-C data to identify potential target genes of lncRNAs. RNA-protein interaction data suggested that nuclear lncRNAs act as scaffolds to recruit regulatory proteins to target promoters and enhancers. Nuclear lncRNAs may therefore play a role in directing regulatory factors to locations spatially close to the lncRNA gene. We provide the analysis results through an interactive visualization web portal at https://fantom.gsc.riken.jp/zenbu/reports/#F6_3D_lncRNA.

Decoding triancestral origins, archaic introgression, and natural selection in the Japanese population by whole-genome sequencing.

We generated Japanese Encyclopedia of Whole-Genome/Exome Sequencing Library (JEWEL), a high-depth whole-genome sequencing dataset comprising 3256 individuals from across Japan. Analysis of JEWEL revealed genetic characteristics of the Japanese population that were not discernible using microarray data. First, rare variant-based analysis revealed an unprecedented fine-scale genetic structure. Together with population genetics analysis, the present-day Japanese can be decomposed into three ancestral components. Second, we identified unreported loss-of-function (LoF) variants and observed that for specific genes, LoF variants appeared to be restricted to a more limited set of transcripts than would be expected by chance, with PTPRD as a notable example. Third, we identified 44 archaic segments linked to complex traits, including a Denisovan-derived segment at NKX6-1 associated with type 2 diabetes. Most of these segments are specific to East Asians. Fourth, we identified candidate genetic loci under recent natural selection. Overall, our work provided insights into genetic characteristics of the Japanese population.

Comparative evaluation of SNVs, indels, and structural variations detected with short- and long-read sequencing data.

Short- and long-read sequencing technologies are routinely used to detect DNA variants, including SNVs, indels, and structural variations (SVs). However, the differences in the quality and quantity of variants detected between short- and long-read data are not fully understood. In this study, we comprehensively evaluated the variant calling performance of short- and long-read-based SNV, indel, and SV detection algorithms (6 for SNVs, 12 for indels, and 13 for SVs) using a novel evaluation framework incorporating manual visual inspection. The results showed that indel-insertion calls greater than 10 bp were poorly detected by short-read-based detection algorithms compared to long-read-based algorithms; however, the recall and precision of SNV and indel-deletion detection were similar between short- and long-read data. The recall of SV detection with short-read-based algorithms was significantly lower in repetitive regions, especially for small- to intermediate-sized SVs, than that detected with long-read-based algorithms. In contrast, the recall and precision of SV detection in nonrepetitive regions were similar between short- and long-read data. These findings suggest the need for refined strategies, such as incorporating multiple variant detection algorithms, to generate a more complete set of variants using short-read data.

Relationship of Major Depressive Disorder and Schizophrenia Polygenic Risk Scores to Suicide: A Comparison Between European and Asian Ancestry Populations.

Psychiatric diagnosis rates in suicide decedents appear higher in European ancestry populations compared with East Asians. Shared genetic components exist between major depressive disorder (MDD)/schizophrenia (SCZ) and suicide, but no study has compared these shared polygenic architectures between Europeans and East Asians. We compared polygenic risk scores (PRSs) for MDD/SCZ determined from large data sets specific to each ancestry in European and East Asian suicide decedent samples. MDD/SCZ PRSs appeared more prominent in European suicides compared with Japanese suicides. A greater coexistence of psychiatric disorders in European suicide decedents than in East Asian suicide decedents may be partly explained by genetics. Our results are limited by the smaller sample size of our suicide decedents and sample size disparities between the European discovery data set and the East Asian data set for MDD/SCZ, resulting in less statistical power to detect robust difference between the two ancestries.

Psychiatric diagnosis in suicides appears more common in Europeans than in East Asians.This is the first comparison of suicide genome-wide association studies between Europeans and East Asians.Major depressive disorder/schizophrenia polygenic risk scores for suicide were more significant for Europeans than for East Asians.

Protein-altering variants at copy number-variable regions influence diverse human phenotypes.

Copy number variants (CNVs) are among the largest genetic variants, yet CNVs have not been effectively ascertained in most genetic association studies. Here we ascertained protein-altering CNVs from UK Biobank whole-exome sequencing data (n = 468,570) using haplotype-informed methods capable of detecting subexonic CNVs and variation within segmental duplications. Incorporating CNVs into analyses of rare variants predicted to cause gene loss of function (LOF) identified 100 associations of predicted LOF variants with 41 quantitative traits. A low-frequency partial deletion of RGL3 exon 6 conferred one of the strongest protective effects of gene LOF on hypertension risk (odds ratio = 0.86 (0.82-0.90)). Protein-coding variation in rapidly evolving gene families within segmental duplications-previously invisible to most analysis methods-generated some of the human genome's largest contributions to variation in type 2 diabetes risk, chronotype and blood cell traits. These results illustrate the potential for new genetic insights from genomic variation that has escaped large-scale analysis to date.

Disease-specific variant interpretation highlighted the genetic findings in 2325 Japanese patients with retinitis pigmentosa and allied diseases.

BACKGROUND: As gene-specific therapy for inherited retinal dystrophy (IRD) advances, unified variant interpretation across institutes is becoming increasingly important. This study aims to update the genetic findings of 86 retinitis pigmentosa (RP)-related genes in a large number of Japanese patients with RP by applying the standardised variant interpretation guidelines for Japanese patients with IRD (J-IRD-VI guidelines) built upon the American College of Medical Genetics and Genomics and the Association for Molecular Pathology rules, and assess the contribution of these genes in RP-allied diseases. METHODS: We assessed 2325 probands with RP (n=2155, including n=1204 sequenced previously with the same sequencing panel) and allied diseases (n=170, newly analysed), including Usher syndrome, Leber congenital amaurosis and cone-rod dystrophy (CRD). Target sequencing using a panel of 86 genes was performed. The variants were interpreted according to the J-IRD-VI guidelines. RESULTS: A total of 3564 variants were detected, of which 524 variants were interpreted as pathogenic or likely pathogenic. Among these 524 variants, 280 (53.4%) had been either undetected or interpreted as variants of unknown significance or benign variants in our earlier study of 1204 patients with RP. This led to a genetic diagnostic rate in 38.6% of patients with RP, with EYS accounting for 46.7% of the genetically solved patients, showing a 9% increase in diagnostic rate from our earlier study. The genetic diagnostic rate for patients with CRD was 28.2%, with RP-related genes significantly contributing over other allied diseases. CONCLUSION: A large-scale genetic analysis using the J-IRD-VI guidelines highlighted the population-specific genetic findings for Japanese patients with IRD; these findings serve as a foundation for the clinical application of gene-specific therapies.

Attenuation of HOIL-1L ligase activity promotes systemic autoimmune disorders by augmenting linear ubiquitin signaling.

Linear ubiquitin chains, which are generated specifically by the linear ubiquitin assembly complex (LUBAC) ubiquitin ligase, play crucial roles in immune signaling, including NF-κB activation. LUBAC comprises catalytic large isoform of heme-oxidized iron regulatory protein 2 ubiquitin ligase 1 (HOIL-1L) interacting protein (HOIP), accessory HOIL-1L, and SHANK-associated RH domain-interacting protein (SHARPIN). Deletion of the ubiquitin ligase activity of HOIL-1L, an accessory ligase of LUBAC, augments LUBAC functions by enhancing LUBAC-mediated linear ubiquitination, which is catalyzed by HOIP. Here, we show that HOIL-1L ΔRING1 mice, which exhibit augmented LUBAC functions upon loss of the HOIL-1L ligase, developed systemic lupus erythematosus (SLE) and Sjögren's syndrome in a female-dominant fashion. Augmented LUBAC activity led to hyperactivation of both lymphoid and myeloid cells. In line with the findings in mice, we sought to identify missense single nucleotide polymorphisms/variations of the RBCK1/HOIL-1L gene in humans that attenuate HOIL-1L ligase activity. We found that the R464H variant, which is encoded by rs774507518 within the RBCK1/HOIL-1L gene, attenuated HOIL-1L ligase activity and augmented LUBAC-mediated immune signaling, including that mediated by Toll-like receptors. We also found that rs774507518 was enriched significantly in patients with SLE, strongly suggesting that RBCK1/HOIL-1L is an SLE susceptibility gene and that augmented linear ubiquitin signaling generated specifically by LUBAC underlies the pathogenesis of this prototype systemic autoimmune disease.

Fine-mapping genomic loci refines bipolar disorder risk genes.

Bipolar disorder (BD) is a heritable mental illness with complex etiology. While the largest published genome-wide association study identified 64 BD risk loci, the causal SNPs and genes within these loci remain unknown. We applied a suite of statistical and functional fine-mapping methods to these loci, and prioritized 22 likely causal SNPs for BD. We mapped these SNPs to genes, and investigated their likely functional consequences by integrating variant annotations, brain cell-type epigenomic annotations, brain quantitative trait loci, and results from rare variant exome sequencing in BD. Convergent lines of evidence supported the roles of SCN2A, TRANK1, DCLK3, INSYN2B, SYNE1, THSD7A, CACNA1B, TUBBP5, PLCB3, PRDX5, KCNK4, AP001453.3, TRPT1, FKBP2, DNAJC4, RASGRP1, FURIN, FES, YWHAE, DPH1, GSDMB, MED24, THRA, EEF1A2, and KCNQ2 in BD. These represent promising candidates for functional experiments to understand biological mechanisms and therapeutic potential. Additionally, we demonstrated that fine-mapping effect sizes can improve performance and transferability of BD polygenic risk scores across ancestrally diverse populations, and present a high-throughput fine-mapping pipeline (https://github.com/mkoromina/SAFFARI).

Association of genetic variation in COL11A1 with adolescent idiopathic scoliosis.

Adolescent idiopathic scoliosis (AIS) is a common and progressive spinal deformity in children that exhibits striking sexual dimorphism, with girls at more than fivefold greater risk of severe disease compared to boys. Despite its medical impact, the molecular mechanisms that drive AIS are largely unknown. We previously defined a female-specific AIS genetic risk locus in an enhancer near the PAX1 gene. Here, we sought to define the roles of PAX1 and newly identified AIS-associated genes in the developmental mechanism of AIS. In a genetic study of 10,519 individuals with AIS and 93,238 unaffected controls, significant association was identified with a variant in COL11A1 encoding collagen (α1) XI (rs3753841; NM_080629.2_c.4004C>T; p.(Pro1335Leu); p=7.07E-11, OR = 1.118). Using CRISPR mutagenesis we generated Pax1 knockout mice (Pax1-/-). In postnatal spines we found that PAX1 and collagen (α1) XI protein both localize within the intervertebral disc-vertebral junction region encompassing the growth plate, with less collagen (α1) XI detected in Pax1-/- spines compared to wild-type. By genetic targeting we found that wild-type Col11a1 expression in costal chondrocytes suppresses expression of Pax1 and of Mmp3, encoding the matrix metalloproteinase 3 enzyme implicated in matrix remodeling. However, the latter suppression was abrogated in the presence of the AIS-associated COL11A1P1335L mutant. Further, we found that either knockdown of the estrogen receptor gene Esr2 or tamoxifen treatment significantly altered Col11a1 and Mmp3 expression in chondrocytes. We propose a new molecular model of AIS pathogenesis wherein genetic variation and estrogen signaling increase disease susceptibility by altering a PAX1-COL11a1-MMP3 signaling axis in spinal chondrocytes.

Adolescent idiopathic scoliosis (AIS) is a twisting deformity of the spine that occurs during periods of rapid growth in children worldwide. Children with severe cases of AIS require surgery to stop it from getting worse, presenting a significant financial burden to health systems and families. Although AIS is known to cluster in families, its genetic causes and its inheritance pattern have remained elusive. Additionally, AIS is known to be more prevalent in females, a bias that has not been explained. Advances in techniques to study the genetics underlying diseases have revealed that certain variations that increase the risk of AIS affect cartilage and connective tissue. In humans, one such variation is near a gene called Pax1, and it is female-specific. The extracellular matrix is a network of proteins and other molecules in the space between cells that help connect tissues together, and it is particularly important in cartilage and other connective tissues. One of the main components of the extracellular matrix is collagen. Yu, Kanshour, Ushiki et al. hypothesized that changes in the extracellular matrix could affect the cartilage and connective tissues of the spine, leading to AIS. To show this, the scientists screened over 100,000 individuals and found that AIS is associated with variants in two genes coding for extracellular matrix proteins. One of these variants was found in a gene called Col11a1, which codes for one of the proteins that makes up collagen. To understand the relationship between Pax1 and Col11a1, Yu, Kanshour, Ushiki et al. genetically modified mice so that they would lack the Pax1 gene. In these mice, the activation of Col11a1 was reduced in the mouse spine. They also found that the form of Col11a1 associated with AIS could not suppress the activation of a gene called Mmp3 in mouse cartilage cells as effectively as unmutated Col11a1. Going one step further, the researchers found that lowering the levels of an estrogen receptor altered the activation patterns of Pax1, Col11a1, and Mmp3 in mouse cartilage cells. These findings suggest a possible mechanism for AIS, particularly in females. The findings of Yu, Kanshour, Ushiki et al. highlight that cartilage cells in the spine are particularly relevant in AIS. The results also point to specific molecules within the extracellular matrix as important for maintaining proper alignment in the spine when children are growing rapidly. This information may guide future therapies aimed at maintaining healthy spinal cells in adolescent children, particularly girls.

The genome-wide association study of serum IgE levels demonstrated a shared genetic background in allergic diseases.

Immunoglobulin E (IgE) synthessis is highly related to a variety of atopic diseases, and several genome-wide association studies (GWASs) have demonstrated the association between genes and IgE level. In this study, we conducted the largest genome-wide association study of IgE involving a Taiwanese Han population. Eight independent variants exhibited genome-wide significance. Among them, an intronic SNP of CD28, rs1181388, and an intergenic SNP, rs1002957030, on 11q23.2 were identified as novel signals for IgE. Seven of the loci were replicated successfully in a meta-analysis using data on Japanese population. Among all the human leukocyte antigen (HLA) regions, HLA-DQA1*03:02 - HLA-DQB1*03:03 was the most significant haplotype (OR = 1.25, SE = 0.02, FDR = 1.6 × 10-14), corresponding to HLA-DQA1 Asp160 and HLA-DQB1 Leu87 amino acid residues. The genetic correlation showed significance between IgE and allergic diseases including asthma, atopic dermatitis, and pollinosis. IgE PRS was significantly correlated with total IgE levels. Furthermore, the top decile IgE polygenic risk score (PRS) group had the highest risk of asthma for the Taiwan Biobank and Biobank Japan cohorts. IgE PRS may be used to aid in predicting the occurrence of allergic reactions before symptoms occur and biomarkers are detectable. Our study provided a more comprehensive understanding of the impact of genomic variants, including complex HLA alleles, on serum IgE levels.

Family History of Developmental Dysplasia of the Hip is a Risk Factor for the Progression of Hip Osteoarthritis.

BACKGROUND: Developmental dysplasia of the hip (DDH) is considered to have genetic predisposition and presents many intrafamilial occurrences. However, there is no report that evaluates the effect of DDH family history on the progression after the onset of hip osteoarthritis (OA). METHODS: Medical interviews about detailed clinical information including family history were conducted on 298 consecutive patients who had undergone surgery for OA due to DDH. Clinical or radiographic items that are associated with the severity of DDH (total hip arthroplasty [THA], involvement of bilateral DDH, onset age of hip pain, and three radiological indices of DDH: center-edge angle, sharp angle, and acetabular roof obliquity) were collected and evaluated in multivariate analyses for their associations with DDH family history in a qualitative or quantitative manner. Survival time analyses for THA as the endpoint was also performed to evaluate the effects of DDH family history on the progression of OA. RESULTS: The DDH family history showed significant associations with bilateral involvement of DDH (odds ratio = 2.09 [95% confidence interval {CI} 1.05 to 4.16]; P = .037), early onset of hip pain (P = .0065), and radiological severity of DDH (P = .016). The DDH family history showed a significant association with undergoing THA (odds ratio = 2.25 [95% CI 1.09 to 4.66]; P = .029), further supported by the Cox regression analyses (hazards ratio = 1.56 [95% CI 1.15 to 2.11]; P = .0044). CONCLUSION: A DDH family history is a risk factor for the progression of hip OA. Stronger genetic predisposition to DDH leads to faster onset and progression of hip OA.

Interrogating Causal Effects of Body Composition and Puberty-Related Risk Factors on Adolescent Idiopathic Scoliosis: A Two-Sample Mendelian Randomization Study.

Adolescent idiopathic scoliosis (AIS) is the most common form of pediatric musculoskeletal disorder. Observational studies have pointed to several risk factors for AIS, but almost no evidence exists to support their causal association with AIS. Here, we applied Mendelian randomization (MR), known to limit bias from confounding and reverse causation, to investigate causal associations between body composition and puberty-related exposures and AIS risk in Europeans and Asians. For our two-sample MR studies, we used single nucleotide polymorphisms (SNPs) associated with body mass index (BMI), waist-hip ratio, lean mass, childhood obesity, bone mineral density (BMD), 25-hydroxyvitamin D (25OHD), age at menarche, and pubertal growth in large European genome-wide association studies (GWAS), and with adult osteoporosis risk and age of menarche in Biobank Japan. We extracted estimates of the aforementioned SNPs on AIS risk from the European or Asian subsets of the largest multiancestry AIS GWAS (N = 7956 cases/88,459 controls). The results of our inverse variance-weighted (IVW) MR estimates suggest no causal association between the aforementioned risk factors and risk of AIS. Pleiotropy-sensitive MR methods yielded similar results. However, restricting our analysis to European females with AIS, we observed a causal association between estimated BMD and the risk of AIS (IVW odds ratio for AIS = 0.1, 95% confidence interval 0.01 to 0.7, p = 0.02 per SD increase in estimated BMD), but this association was no longer significant after adjusting for BMI, body fat mass, and 25OHD and remained significant after adjusting for age at menarche in multivariable MR. In conclusion, we demonstrated a protective causal effect of BMD on AIS risk in females of European ancestry, but this effect was modified by BMI, body fat mass, and 25OHD levels. Future MR studies using larger AIS GWAS are needed to investigate small effects of the aforementioned exposures on AIS. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Genetically proxied HTRA1 protease activity and circulating levels independently predict risk of ischemic stroke and coronary artery disease.

HTRA1 has emerged as a major risk gene for stroke and cerebral small vessel disease with both rare and common variants contributing to disease risk. However, the precise mechanisms mediating this risk remain largely unknown as does the full spectrum of phenotypes associated with genetic variation in HTRA1 in the general population. Using a family-history informed approach, we first show that rare variants in HTRA1 are linked to ischemic stroke in 425,338 European individuals from the UK Biobank with replication in 143,149 individuals from the Biobank Japan. Integrating data from biochemical experiments on 76 mutations occurring in the UK Biobank, we next show that rare variants causing loss of protease function in vitro associate with ischemic stroke, coronary artery disease, and skeletal traits. In addition, a common causal variant (rs2672592) modulating circulating HTRA1 mRNA and protein levels enhances the risk of ischemic stroke, small vessel stroke, and coronary artery disease while lowering the risk of migraine and age-related macular dystrophy in GWAS and UK Biobank data from > 2,000,000 individuals. There was no evidence of an interaction between genetically proxied HTRA1 activity and levels. Our findings demonstrate a central role of HTRA1 for human disease including stroke and coronary artery disease and identify two independent mechanisms that might qualify as targets for future therapeutic interventions.

Natural Selection Signatures in the Hondo and Ryukyu Japanese Subpopulations.

Natural selection signatures across Japanese subpopulations are under-explored. Here we conducted genome-wide selection scans with 622,926 single nucleotide polymorphisms for 20,366 Japanese individuals, who were recruited from the main-islands of Japanese Archipelago (Hondo) and the Ryukyu Archipelago (Ryukyu), representing two major Japanese subpopulations. The integrated haplotype score (iHS) analysis identified several signals in one or both subpopulations. We found a novel candidate locus at IKZF2, especially in Ryukyu. Significant signals were observed in the major histocompatibility complex region in both subpopulations. The lead variants differed and demonstrated substantial allele frequency differences between Hondo and Ryukyu. The lead variant in Hondo tags HLA-A*33:03-C*14:03-B*44:03-DRB1*13:02-DQB1*06:04-DPB1*04:01, a haplotype specific to Japanese and Korean. While in Ryukyu, the lead variant tags DRB1*15:01-DQB1*06:02, which had been recognized as a genetic risk factor for narcolepsy. In contrast, it is reported to confer protective effects against type 1 diabetes and human T lymphotropic virus type 1-associated myelopathy/tropical spastic paraparesis. The FastSMC analysis identified 8 loci potentially affected by selection within the past 20-150 generations, including 2 novel candidate loci. The analysis also showed differences in selection patterns of ALDH2 between Hondo and Ryukyu, a gene recognized to be specifically targeted by selection in East Asian. In summary, our study provided insights into the selection signatures within the Japanese and nominated potential sources of selection pressure.

Androgen receptor binding sites enabling genetic prediction of mortality due to prostate cancer in cancer-free subjects.

Prostate cancer (PrCa) is the second most common cancer worldwide in males. While strongly warranted, the prediction of mortality risk due to PrCa, especially before its development, is challenging. Here, we address this issue by maximizing the statistical power of genetic data with multi-ancestry meta-analysis and focusing on binding sites of the androgen receptor (AR), which has a critical role in PrCa. Taking advantage of large Japanese samples ever, a multi-ancestry meta-analysis comprising more than 300,000 subjects in total identifies 9 unreported loci including ZFHX3, a tumor suppressor gene, and successfully narrows down the statistically finemapped variants compared to European-only studies, and these variants strongly enrich in AR binding sites. A polygenic risk scores (PRS) analysis restricting to statistically finemapped variants in AR binding sites shows among cancer-free subjects, individuals with a PRS in the top 10% have a strongly higher risk of the future death of PrCa (HR: 5.57, P = 4.2 × 10-10). Our findings demonstrate the potential utility of leveraging large-scale genetic data and advanced analytical methods in predicting the mortality of PrCa.

Genetic risk score of cerebral infarction in atrial fibrillation genome-wide association study.

INTRODUCTION: Stroke is a leading cause of death and the primary cause of adult-acquired disability. Patients with cardiogenic embolic stroke also have higher mortality and recurrence rates than patients with other stroke subtypes. Atrial fibrillation (AF) is a major risk factor for cerebral infarction (CI). The large-scale study identified 32 loci in the MEGASTROKE study. However, few studies have attempted to identify novel stroke risk variants in patients with a history of AF. Our overall aim was to identify novel CI risk variants in AF cases and explore whether their associations with the CI risk were affected by the CHADS2 and CHA2DS2-VASc scores. METHODS: We performed association study with CI using 8181 AF cases in previous genome-wide association study (GWAS) and imputation data without controls. We classified AF cases into those with or without past history of CI, and the genetic associations with the CI risk were examined. RESULTS: GWAS identified eight associated loci. The generated genetic risk score (GRS) for the eight loci was significantly associated with CI in patients with AF (1.46 × 10-8 ). We estimated bivariate logistic regression model which contained GRS and CHADS2 score (GRS: p-Value = 7.41 × 10-9 , CHADS2 score: p-Value <2.0 × 10-16 ) or CHA2DS2-VASc scores (GRS: p-Value = 2.52 × 10-10 , CHA2DS2-VASc score: p-Value <2.0 × 10-16 ). CONCLUSION: We identified eight genetic variants that were potentially associated with the risk of CI of AF cases and the significant GRS, whose associations were independent of the CHADS2 or CHA2DS2-VASc score.

Understanding the genetic complexity of puberty timing across the allele frequency spectrum.

Pubertal timing varies considerably and has been associated with a range of health outcomes in later life. To elucidate the underlying biological mechanisms, we performed multi-ancestry genetic analyses in ~800,000 women, identifying 1,080 independent signals associated with age at menarche. Collectively these loci explained 11% of the trait variance in an independent sample, with women at the top and bottom 1% of polygenic risk exhibiting a ~11 and ~14-fold higher risk of delayed and precocious pubertal development, respectively. These common variant analyses were supported by exome sequence analysis of ~220,000 women, identifying several genes, including rare loss of function variants in ZNF483 which abolished the impact of polygenic risk. Next, we implicated 660 genes in pubertal development using a combination of in silico variant-to-gene mapping approaches and integration with dynamic gene expression data from mouse embryonic GnRH neurons. This included an uncharacterized G-protein coupled receptor GPR83, which we demonstrate amplifies signaling of MC3R, a key sensor of nutritional status. Finally, we identified several genes, including ovary-expressed genes involved in DNA damage response that co-localize with signals associated with menopause timing, leading us to hypothesize that the ovarian reserve might signal centrally to trigger puberty. Collectively these findings extend our understanding of the biological complexity of puberty timing and highlight body size dependent and independent mechanisms that potentially link reproductive timing to later life disease.

Evidence of causality of low body mass index on risk of adolescent idiopathic scoliosis: a Mendelian randomization study.

Introduction: Adolescent idiopathic scoliosis (AIS) is a disorder with a three-dimensional spinal deformity and is a common disease affecting 1-5% of adolescents. AIS is also known as a complex disease involved in environmental and genetic factors. A relation between AIS and body mass index (BMI) has been epidemiologically and genetically suggested. However, the causal relationship between AIS and BMI remains to be elucidated. Material and methods: Mendelian randomization (MR) analysis was performed using summary statistics from genome-wide association studies (GWASs) of AIS (Japanese cohort, 5,327 cases, 73,884 controls; US cohort: 1,468 cases, 20,158 controls) and BMI (Biobank Japan: 173430 individual; meta-analysis of genetic investigation of anthropometric traits and UK Biobank: 806334 individuals; European Children cohort: 39620 individuals; Population Architecture using Genomics and Epidemiology: 49335 individuals). In MR analyses evaluating the effect of BMI on AIS, the association between BMI and AIS summary statistics was evaluated using the inverse-variance weighted (IVW) method, weighted median method, and Egger regression (MR-Egger) methods in Japanese. Results: Significant causality of genetically decreased BMI on risk of AIS was estimated: IVW method (Estimate (beta) [SE] = -0.56 [0.16], p = 1.8 × 10-3), weighted median method (beta = -0.56 [0.18], p = 8.5 × 10-3) and MR-Egger method (beta = -1.50 [0.43], p = 4.7 × 10-3), respectively. Consistent results were also observed when using the US AIS summary statistic in three MR methods; however, no significant causality was observed when evaluating the effect of AIS on BMI. Conclusions: Our Mendelian randomization analysis using large studies of AIS and GWAS for BMI summary statistics revealed that genetic variants contributing to low BMI have a causal effect on the onset of AIS. This result was consistent with those of epidemiological studies and would contribute to the early detection of AIS.

Genetic insights into ossification of the posterior longitudinal ligament of the spine.

Ossification of the posterior longitudinal ligament of the spine (OPLL) is an intractable disease leading to severe neurological deficits. Its etiology and pathogenesis are primarily unknown. The relationship between OPLL and comorbidities, especially type 2 diabetes (T2D) and high body mass index (BMI), has been the focus of attention; however, no trait has been proven to have a causal relationship. We conducted a meta-analysis of genome-wide association studies (GWASs) using 22,016 Japanese individuals and identified 14 significant loci, 8 of which were previously unreported. We then conducted a gene-based association analysis and a transcriptome-wide Mendelian randomization approach and identified three candidate genes for each. Partitioning heritability enrichment analyses observed significant enrichment of the polygenic signals in the active enhancers of the connective/bone cell group, especially H3K27ac in chondrogenic differentiation cells, as well as the immune/hematopoietic cell group. Single-cell RNA sequencing of Achilles tendon cells from a mouse Achilles tendon ossification model confirmed the expression of genes in GWAS and post-GWAS analyses in mesenchymal and immune cells. Genetic correlations with 96 complex traits showed positive correlations with T2D and BMI and a negative correlation with cerebral aneurysm. Mendelian randomization analysis demonstrated a significant causal effect of increased BMI and high bone mineral density on OPLL. We evaluated the clinical images in detail and classified OPLL into cervical, thoracic, and the other types. GWAS subanalyses identified subtype-specific signals. A polygenic risk score for BMI demonstrated that the effect of BMI was particularly strong in thoracic OPLL. Our study provides genetic insight into the etiology and pathogenesis of OPLL and is expected to serve as a basis for future treatment development.