Functional and dynamic profiling of transcript isoforms reveals essential roles of alternative splicing in interferon response.

Type I interferon (IFN-I) plays an important role in the innate immune response through inducing IFN-I-stimulated genes (ISGs). However, how alternative splicing (AS) events, especially over time, affect their function remains poorly understood. We generated an annotation (113,843 transcripts) for IFN-I-stimulated human B cells called isoISG using high-accuracy long-read sequencing data from PacBio Sequel II/IIe. Transcript isoform profiling using isoISG revealed that isoform switching occurred in the early response to IFN-I so that ISGs would gain functional domains (e.g., C4B) or higher protein production (e.g., IRF3). Conversely, isoforms lacking functional domains increased during the late phase of IFN-I response, mainly due to intron retention events. This suggests that isoform switching both triggers and terminates IFN-I responses at the translation and protein levels. Furthermore, genetic variants influencing the isoform ratio of ISGs were associated with immunological and infectious diseases. AS has essential roles in regulating innate immune response and associated diseases.

The muscle tissue transcriptome of idiopathic inflammatory myopathy reflects the muscle damage process by monocytes and presence of skin lesions.

OBJECTIVE: To investigate transcriptomic and immunophenotypic features of muscle specimens from patients with idiopathic inflammatory myopathy (IIM). METHODS: Bulk RNA-sequencing was performed on muscle biopsy samples from 16 patients with dermatomyositis (DM) and 9 patients with polymyositis (PM). Seven tested positive for anti-aminoacyl t-RNA synthetase antibodies in the DM patients (ARS-DM). We conducted weighted gene coexpression network analysis (WGCNA), differentially expressed gene (DEG) analysis, and gene set variation analysis (GSVA) to assess contributions of specific pathways. Cell proportions in muscle specimens were estimated using a deconvolution approach. RESULTS: WGCNA revealed significant positive correlations between serum creatine kinase (CK) levels and gene modules involved in cellular respiration, phagocytosis, and oxidative phosphorylation (OXPHOS). Significant positive correlations were also observed between CK levels and proportions of CD16-positive and -negative monocytes and myeloid dendritic cells. Notably, DM patients demonstrated enrichment of complement and interferon-α and -γ pathway genes compared to those with PM. Furthermore, ARS-DM demonstrated a higher proportion of Th1 cells and DEGs related to OXPHOS. Additionally, serum Krebs von den Lungen-6 levels correlated with gene modules associated with extracellular matrix and transforming growth factor-ß signaling pathway. CONCLUSION: Our study highlights a significant involvement of monocytes in muscle damage and delineates pathological differences among IIM subtypes. DM was characterized by complement, interferon-α and -γ signaling, whilst ARS-DM was associated with OXPHOS. Distinctive gene expression variations in muscle specimens suggest that different pathologic mechanisms underlie muscle damage in each IIM phenotype.

Humoral responses are enhanced by facilitating B cell viability by Fcrl5 overexpression in B cells.

B cell initial activity is regulated through a balance of activation and suppression mediated by regulatory molecules expressed in B cells; however, the molecular mechanisms underlying this process remain incompletely understood. In this study, we investigated the function of the Fc receptor-like (Fcrl) family molecule Fcrl5, which is constitutively expressed in naive B cells, in humoral immune responses. Our study demonstrated that B cell-specific overexpression of Fcrl5 enhanced antibody (Ab) production in both T cell-independent type 1 (TI1) and T cell-dependent (TD) responses. Additionally, it promoted effector B cell formation under competitive conditions in TD responses. Mechanistically, in vitro ligation of Fcrl5 by agonistic Abs reduced cell death and enhanced proliferation in lipopolysaccharide-stimulated B cells. In the presence of anti-CD40 Abs and IL-5, the Fcrl5 ligation not only suppressed cell death but also enhanced differentiation into plasma cells. These findings reveal a novel role of Fcrl5 in promoting humoral immune responses by enhancing B cell viability and plasma cell differentiation.

Long-read sequencing for 29 immune cell subsets reveals disease-linked isoforms.

Alternative splicing events are a major causal mechanism for complex traits, but they have been understudied due to the limitation of short-read sequencing. Here, we generate a full-length isoform annotation of human immune cells from an individual by long-read sequencing for 29 cell subsets. This contains a number of unannotated transcripts and isoforms such as a read-through transcript of TOMM40-APOE in the Alzheimer's disease locus. We profile characteristics of isoforms and show that repetitive elements significantly explain the diversity of unannotated isoforms, providing insight into the human genome evolution. In addition, some of the isoforms are expressed in a cell-type specific manner, whose alternative 3'-UTRs usage contributes to their specificity. Further, we identify disease-associated isoforms by isoform switch analysis and by integration of several quantitative trait loci analyses with genome-wide association study data. Our findings will promote the elucidation of the mechanism of complex diseases via alternative splicing.

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.

Upregulated Fcrl5 disrupts B cell anergy and causes autoimmune disease.

B cell anergy plays a critical role in maintaining self-tolerance by inhibiting autoreactive B cell activation to prevent autoimmune diseases. Here, we demonstrated that Fc receptor-like 5 (Fcrl5) upregulation contributes to autoimmune disease pathogenesis by disrupting B cell anergy. Fcrl5-a gene whose homologs are associated with human autoimmune diseases-is highly expressed in age/autoimmunity-associated B cells (ABCs), an autoreactive B cell subset. By generating B cell-specific Fcrl5 transgenic mice, we demonstrated that Fcrl5 overexpression in B cells caused systemic autoimmunity with age. Additionally, Fcrl5 upregulation in B cells exacerbated the systemic lupus erythematosus-like disease model. Furthermore, an increase in Fcrl5 expression broke B cell anergy and facilitated toll-like receptor signaling. Thus, Fcrl5 is a potential regulator of B cell-mediated autoimmunity by regulating B cell anergy. This study provides important insights into the role of Fcrl5 in breaking B cell anergy and its effect on the pathogenesis of autoimmune diseases.

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.

Expression of the readthrough transcript CiDRE in alveolar macrophages boosts SARS-CoV-2 susceptibility and promotes COVID-19 severity.

Lung infection during severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) via the angiotensin-I-converting enzyme 2 (ACE2) receptor induces a cytokine storm. However, the precise mechanisms involved in severe COVID-19 pneumonia are unknown. Here, we showed that interleukin-10 (IL-10) induced the expression of ACE2 in normal alveolar macrophages, causing them to become vectors for SARS-CoV-2. The inhibition of this system in hamster models attenuated SARS-CoV-2 pathogenicity. Genome-wide association and quantitative trait locus analyses identified a IFNAR2-IL10RB readthrough transcript, COVID-19 infectivity-enhancing dual receptor (CiDRE), which was highly expressed in patients harboring COVID-19 risk variants at the IFNAR2 locus. We showed that CiDRE exerted synergistic effects via the IL-10-ACE2 axis in alveolar macrophages and functioned as a decoy receptor for type I interferons. Collectively, our data show that high IL-10 and CiDRE expression are potential risk factors for severe COVID-19. Thus, IL-10R and CiDRE inhibitors might be useful COVID-19 therapies.

Rare germline variants in pancreatic cancer and multiple primary cancers: an autopsy study.

BACKGROUND: There is a lack of information on rare germline variants of pancreatic cancer-predisposing genes. Risk genes for multiple primary cancers may overlap with those for pancreatic cancer. METHODS: A retrospective study of autopsy cases with a negative family history in the Japanese single nucleotide polymorphism for geriatric research database examined rare germline variants in the protein-coding regions of 61 genes. Targeted sequencing of these genes was performed and classified for pathogenicity using the American College of Medical Genetics and Genomics guidelines. Polyphen-2, SIFT and LoFtool algorithms were used to predict damage to protein function. RESULTS: Of the 189 subjects used (90 cancer and 99 non-cancer controls), 72 patients had pancreatic cancer (23 had multiple primary cancers) and 18 had no pancreatic cancer in multiple primary cancers. APC, BRCA2, BUB1B, ENG and MSH6 were associated with cancer predisposition, and pathogenic/likely pathogenic (P/LP) variants occurred in 6% [pancreatic cancer (4/72); all-cancer (5/90)] and 54% (49/90) carried only variants of uncertain significance (VUS) among cancer patients. Of these VUS, in pancreatic cancer patients, four DNA mismatch repair (MMR) genes ( MLH1, MSH2, MSH6 and PMS2 ), and POLQ in men were significantly associated (odds ratio = 3.83; P = 0.025; P = 0.027, respectively). The most abundant predictor of functionally damaging variants was POLQ . CONCLUSIONS: The frequency of P/LP variants in patients with sporadic pancreatic cancer suggests the need for genetic evaluation of individuals with no family history. VUS of MMR genes ( MLH1, MSH2, MSH6 and PMS2 ) and POLQ may be useful in predicting genetic trends in the potential risk of pancreatic cancer, especially in individuals lacking P/LP.

Combined plasma metabolomic and transcriptomic analysis identify histidine as a biomarker and potential contributor in SLE pathogenesis.

OBJECTIVES: To investigate metabolite alterations in the plasma of SLE patients to identify novel biomarkers and provide insight into SLE pathogenesis. METHODS: Patients with SLE (n = 41, discovery cohort and n = 37, replication cohort), healthy controls (n = 30 and n = 29) and patients with RA (n = 19, disease control) were recruited. Metabolic profiles of the plasma samples were analysed using liquid chromatography-time-of-flight mass spectrometry and capillary electrophoresis-time-of-flight mass spectrometry. Transcriptome data was analysed using RNA-sequencing for 18 immune cell subsets. The importance of histidine (His) in plasmablast differentiation was investigated by using mouse splenic B cells. RESULTS: We demonstrate that a specific amino acid combination including His can effectively distinguish between SLE patients and healthy controls. Random forest and partial least squares-discriminant analysis identified His as an effective classifier for SLE patients. A decrease in His plasma levels correlated with damage accrual independent of prednisolone dosage and type I IFN signature. The oxidative phosphorylation signature in plasmablasts negatively correlated with His levels. We also showed that plasmablast differentiation induced by innate immune signals was dependent on His. CONCLUSIONS: Plasma His levels are a potential biomarker for SLE patients and are associated with damage accrual. Our data suggest the importance of His as a pathogenic metabolite in SLE pathogenesis.

Uncovering the Localization and Function of a Novel Read-Through Transcript 'TOMM40-APOE'.

Recent advancements in genome analysis technology have revealed the presence of read-through transcripts in which transcription continues by skipping the polyA signal. We here identified and characterized a new read-through transcript, TOMM40-APOE. With cDNA amplification from THP-1 cells, the TOMM40-APOE3 product was successfully generated. We also generated TOMM40-APOE4, another isoform, by introducing point mutations. Notably, while APOE3 and APOE4 exhibited extracellular secretion, both TOMM40-APOE3 and TOMM40-APOE4 were localized exclusively to the mitochondria. But functionally, they did not affect mitochondrial membrane potential. Cell death induction studies illustrated increased cell death with TOMM40-APOE3 and TOMM40-APOE4, and we did not find any difference in cellular function between the two isoforms. These findings indicated that the new mitochondrial protein TOMM40-APOE has cell toxic ability.

Utility of tissue-specific gene expression scores for gene prioritization in Mendelian diseases.

In genetic testing of Mendelian diseases, it is a bioinformatics challenge to effectively prioritize disease-causing candidate genes listed from massively parallel sequencing. Tissue specificity of the gene expression levels may give a clue because it may reflect tissue-specific disease manifestation. However, considering poor correlations between mRNA and protein expression in some genes, it is not clear whether transcriptomics- or proteomics-based tissue specificity should be used to prioritize candidate genes. Therefore, we compared the efficiency of tissue-specific scores (TS scores) obtained from transcriptome and proteome data in prioritizing candidate genes for whole exome sequencing (WES) analysis of Mendelian disease patients. We show that both Protein and RNA TS scores are useful in prioritizing candidate genes in WES analysis, although diseases like coagulopathies get more benefit from Protein TS score. This study may provide useful evidence in developing new methods to effectively identify novel disease-causing genes.

Splicing QTL analysis focusing on coding sequences reveals mechanisms for disease susceptibility loci.

Splicing quantitative trait loci (sQTLs) are one of the major causal mechanisms in genome-wide association study (GWAS) loci, but their role in disease pathogenesis is poorly understood. One reason is the complexity of alternative splicing events producing many unknown isoforms. Here, we propose two approaches, namely integration and selection, for this complexity by focusing on protein-structure of isoforms. First, we integrate isoforms with the same coding sequence (CDS) and identify 369-601 integrated-isoform ratio QTLs (i2-rQTLs), which altered protein-structure, in six immune subsets. Second, we select CDS incomplete isoforms annotated in GENCODE and identify 175-337 isoform-ratio QTL (i-rQTL). By comprehensive long-read capture RNA-sequencing among these incomplete isoforms, we reveal 29 full-length isoforms with unannotated CDSs associated with GWAS traits. Furthermore, we show that disease-causal sQTL genes can be identified by evaluating their trans-eQTL effects. Our approaches highlight the understudied role of protein-altering sQTLs and are broadly applicable to other tissues and diseases.

47XXY and 47XXX in Scleroderma and Myositis.

OBJECTIVE: We undertook this study to examine the X chromosome complement in participants with systemic sclerosis (SSc) as well as idiopathic inflammatory myopathies. METHODS: The participants met classification criteria for the diseases. All participants underwent single-nucleotide polymorphism typing. We examined X and Y single-nucleotide polymorphism heterogeneity to determine the number of X chromosomes. For statistical comparisons, we used χ2 analyses with calculation of 95% confidence intervals. RESULTS: Three of seventy men with SSc had 47,XXY (P = 0.0001 compared with control men). Among the 435 women with SSc, none had 47,XXX. Among 709 men with polymyositis or dermatomyositis (PM/DM), seven had 47,XXY (P = 0.0016), whereas among the 1783 women with PM/DM, two had 47,XXX. Of 147 men with inclusion body myositis (IBM), six had 47,XXY, and 1 of the 114 women with IBM had 47,XXX. For each of these myositis disease groups, the excess 47,XXY and/or 47,XXX was significantly higher compared with in controls as well as the known birth rate of Klinefelter syndrome or 47,XXX. CONCLUSION: Klinefelter syndrome (47,XXY) is associated with SSc and idiopathic inflammatory myopathies, similar to other autoimmune diseases with type 1 interferon pathogenesis, namely, systemic lupus erythematosus and Sjögren syndrome.

Immune cell multiomics analysis reveals contribution of oxidative phosphorylation to B-cell functions and organ damage of lupus.

OBJECTIVE: Systemic lupus erythematosus (SLE) is the prototypical systemic autoimmune disease. While the long-term prognosis has greatly improved, better long-term survival is still necessary. The type I interferon (IFN) signature, a prominent feature of SLE, is not an ideal therapeutic target or outcome predictor. To explore immunological pathways in SLE more precisely, we performed transcriptomic, epigenomic and genomic analyses using 19 immune cell subsets from peripheral blood. METHODS: We sorted 19 immune cell subsets and identified the mRNA expression profiles and genetic polymorphisms in 107 patients with SLE and 92 healthy controls. Combined differentially expressed genes and expression quantitative trait loci analysis was conducted to find key driver genes in SLE pathogenesis. RESULTS: We found transcriptomic, epigenetic and genetic importance of oxidative phosphorylation (OXPHOS)/mitochondrial dysfunction in SLE memory B cells. Particularly, we identified an OXPHOS-regulating gene, PRDX6 (peroxiredoxin 6), as a key driver in SLE B cells. Prdx6-deficient B cells showed upregulated mitochondrial respiration as well as antibody production. We revealed OXPHOS signature was associated with type I IFN signalling-related genes (ISRGs) signature in SLE memory B cells. Furthermore, the gene sets related to innate immune signalling among ISRGs presented correlation with OXPHOS and these two signatures showed associations with SLE organ damage as well as specific clinical phenotypes. CONCLUSION: This work elucidated the potential prognostic marker for SLE. Since OXPHOS consists of the electron transport chain, a functional unit in mitochondria, these findings suggest the importance of mitochondrial dysfunction as a key immunological pathway involved in SLE.

Association of Polygenic Risk Scores With Radiographic Progression in Patients With Rheumatoid Arthritis.

OBJECTIVE: To investigate whether polygenic risk scores obtained using data from a genome-wide association study (GWAS) of rheumatoid arthritis (RA) susceptibility can be predictors of radiographic progression. METHODS: We constructed polygenic risk scores using GWAS summary data on associations of single-nucleotide polymorphisms with RA susceptibility. The polygenic risk scores were stratified into quintiles based on levels of significance (ranging from top quintile of polygenic risk scores to bottom quintile). In addition, change in the Sharp/van der Heijde score (SHS) of radiographic progression over the first 5 years after onset of RA was assessed. The change in SHS over 5 years was stratified according to quartiles, with the top quartile of change in SHS defined as severe radiographic progression (score change of >35 points) and the remaining 3 quartiles defined as nonsevere radiographic progression. Polygenic risk scores were assessed for their ability to predict the SHS status over 5 years in a training set (n = 500 RA patients) for selection of the best model, and in a testing set (n = 740 RA patients) for validation of the data. We evaluated the performance of the polygenic risk score as a predictor of severe radiographic progression in univariable and multivariable analyses with inclusion of other factors. RESULTS: Polygenic risk scores constructed from 43,784 single-nucleotide polymorphisms significantly differed between patients who experienced severe radiographic progression and those with nonsevere radiographic progression in both the training set (P = 0.0064) and the testing set (P = 0.017). Patients with polygenic risk scores in the top quintile had a higher risk of severe progression compared to those with polygenic risk scores in the bottom quintile (odds ratio [OR] 1.90, P = 0.0022), and the risk of severe radiographic progression was even higher when restricted to patients who were younger at disease onset (OR 5.06, P = 0.00038). The group with polygenic risk scores in the top quintile and the anti-citrullinated protein antibody (ACPA)-positive group had significantly higher proportions of patients with severe radiographic progression (P = 0.00052 and P = 0.0022, respectively) compared to the remaining groups. Multivariable analysis showed that polygenic risk score (P = 0.00019) as well as female sex (P = 0.0033), ACPA positivity (P = 0.0023), and body mass index (P = 0.024) were independent risk factors for severe radiographic progression. CONCLUSION: A polygenic risk score that is derived from GWAS data on RA susceptibility is associated with the level of severity of radiographic progression in patients with RA.

Anti-Double-Stranded DNA Antibodies Recognize DNA Presented on HLA Class II Molecules of Systemic Lupus Erythematosus Risk Alleles.

OBJECTIVE: Specific HLA class II alleles are associated with susceptibility to systemic lupus erythematosus (SLE). The role of HLA class II molecules in SLE pathogenesis remains unclear, although anti-DNA antibodies are specific to SLE and correlate with disease activity. We previously demonstrated that misfolded proteins bound to HLA class II molecules are specific targets for the autoantibodies produced in autoimmune diseases. This study was undertaken to validate our hypothesis that DNA binds to HLA class II molecules in a manner similar to that of misfolded proteins and that DNA bound to HLA class II molecules is involved in SLE pathogenesis. METHODS: We analyzed the binding of DNA to HLA class II molecules, as well as the response of cells expressing anti-DNA B cell receptors (BCRs) to cells expressing the DNA/HLA class II complex. RESULTS: Efficient binding of DNA to HLA class II molecules was observed in risk alleles of SLE, such as HLA-DRB1*15:01. The efficiency of DNA binding to each HLA-DR allele was positively associated with the risk of SLE conferred by the HLA-DR allele. In addition, reporter cells carrying anti-DNA BCRs were activated by cells expressing DNA/HLA class II complexes. CONCLUSION: These results provide evidence that DNA bound to HLA class II molecules is involved in SLE pathogenesis.

Eight novel susceptibility loci and putative causal variants in atopic dermatitis.

BACKGROUND: Atopic dermatitis (AD) is the most common allergic disease in the world. While genetic components play critical roles in its pathophysiology, a large proportion of its genetic background is still unexplored. OBJECTIVES: This study sought to illuminate the genetic associations with AD using genome-wide association study (GWAS) and its downstream analyses. METHODS: This study conducted a GWAS for AD comprising 2,639 cases and 115,648 controls in the Japanese population, followed by a trans-ethnic meta-analysis with UK Biobank data and downstream analyses including partitioning heritability analysis by linkage disequilibrium score regression. RESULTS: This study identified 17 significant susceptibility loci, among which 4 loci-AFF1, ITGB8, EHMT1, and EGR2-were novel in the Japanese GWAS. The trans-ethnic meta-analysis revealed 4 additional novel loci, namely-ZBTB38,LOC105755953/LOC101928272, TRAF3, andIQGAP1. This study found a missense variant (R243W) with a deleterious functional effect in NLRP10 and a variant altering expression of CCDC80 via enhancer expression as highly likely causal variants. These 2 regions were Asian-specific, and these population-specific associations could be explained by the frequency of causal variants. The gene-based test showed SMAD4 as an additional novel significant locus. Downstream analyses revealed substantial overlap of GWAS significant signals in enhancers of skin cells and immune cells, especially CD4 T cells. A highly shared polygenic architecture of AD between Europeans and Asians was also found. CONCLUSIONS: This study identified Japanese-specific loci and novel significant loci shared by different populations. Two putative causal variants were illuminated in Japanese-specific loci. Trans-ethnic analyses revealed strong heritability enrichment in immune-related pathways, and relevant cell types shared among populations.

Polygenic Risk Score of Adolescent Idiopathic Scoliosis for Potential Clinical Use.

Adolescent idiopathic scoliosis (AIS) is a common disease causing three-dimensional spinal deformity in as many as 3% of adolescents. Development of a method that can accurately predict the onset and progression of AIS is an immediate need for clinical practice. Because the heritability of AIS is estimated as high as 87.5% in twin studies, prediction of its onset and progression based on genetic data is a promising option. We show the usefulness of polygenic risk score (PRS) for the prediction of onset and progression of AIS. We used AIS genomewide association study (GWAS) data comprising 79,211 subjects in three cohorts and constructed a PRS based on association statistics in a discovery set including 31,999 female subjects. After calibration using a validation data set, we applied the PRS to a test data set. By integrating functional annotations showing heritability enrichment in the selection of variants, the PRS demonstrated an association with AIS susceptibility (p = 3.5 × 10-40 with area under the receiver-operating characteristic [AUROC] = 0.674, sensitivity = 0.644, and specificity = 0.622). The decile with the highest PRS showed an odds ratio of as high as 3.36 (p = 1.4 × 10-10 ) to develop AIS compared with the fifth in decile. The addition of a predictive model with only a single clinical parameter (body mass index) improved predictive ability for development of AIS (AUROC = 0.722, net reclassification improvement [NRI] 0.505 ± 0.054, p = 1.6 × 10-8 ), potentiating clinical use of the prediction model. Furthermore, we found the Cobb angle (CA), the severity measurement of AIS, to be a polygenic trait that showed a significant genetic correlation with AIS susceptibility (rg = 0.6, p = 3.0 × 10-4 ). The AIS PRS demonstrated a significant association with CA. These results indicate a shared polygenic architecture between onset and progression of AIS and the potential usefulness of PRS in clinical settings as a predictor to promote early intervention of AIS and avoid invasive surgery. © 2021 American Society for Bone and Mineral Research (ASBMR).