Search

Genome-wide association studies (GWAS) have identified genetic variations associated with adverse drug effects in pharmacogenomics (PGx) research. However, interpreting the biological implications of these associations remains a challenge. This review highlights 2 promising post-GWAS methods for PGx. First, we discuss the polygenic architecture of the PGx traits, especially for drug-induced liver injury. Experimental modelling using multiple donors' human primary hepatocytes and human liver organoids demonstrated the polygenic architecture of drug-induced liver injury susceptibility and found biological vulnerability in genetically high-risk tissue donors. Second, we discuss the challenges of interpreting the roles of variants in noncoding regions. Beyond methods involving expression quantitative trait locus analysis and massively parallel reporter assays, we suggest the use of in silico mutagenesis through machine learning methods to understand the roles of variants in transcriptional regulation. This review underscores the importance of these post-GWAS methods in providing critical insights into PGx, potentially facilitating drug development and personalized treatment., (© 2023 The Author. British Journal of Clinical Pharmacology published by John Wiley & Sons Ltd on behalf of British Pharmacological Society.)

The FANTOM web resource (https://fantom.gsc.riken.jp/) has been a unique resource for studying mammalian genomes, which is built on the research activities conducted in the international collaborative project FANTOM (Functional ANnoTation Of the Mammalian genome). In recent updates, we expanded annotations for long non-coding RNAs (lncRNAs) and transcribed cis-regulatory elements (CREs). The former was derived from the large-scale lncRNA perturbations in induced pluripotent stem cells (iPSCs) and integrative analysis of Hi-C data conducted in the sixth iteration of the project (FANTOM6). The resulting annotations of lncRNAs, according to the impact on cellular and molecular phenotypes and the potential RNA-chromatin interactions, are accessible via the interactive ZENBU-Reports framework. The latter involves a new platform, fanta.bio (https://fanta.bio/), which collects transcribed CREs identified via use of an extended dataset of CAGE profiles. The CREs, with their annotations including genetic and epigenetic information, are accessible via a dedicated interface as well as the UCSC Genome Browser Database. These updates offer enhanced opportunities to investigate the functions of non-coding regions within mammalian genomes., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

To improve imputation quality for genome-wide association studies (GWAS) conducted on the Japanese population, we developed and evaluated four Japanese population-specific reference panels. These panels were constructed through the augmentation of the 1000 Genomes Project (1KG) panel using Japanese whole genome sequencing (WGS) data, with sample sizes ranging from 1 K to 7 K individuals enrolled through the Biobank Japan (BBJ) project, and sequencing depths ranging from 3× to 30×. Among these panels, an augmented reference panel comprising 7472 WGS samples of mixed depth (1KG+7K) exhibit the greatest improvement in imputation quality relative to the Trans-Omics for Precision Medicine (TOPMed) reference panel. Notably, we observe these improvements primarily for rare variants with a minor allele frequency (MAF) <5%. To demonstrate the benefits of improved imputation quality in association analyses of complex traits, we conducted GWAS for serum uric acid and total cholesterol levels following imputation up to the 1KG+7K panel. The analysis reveals several loci reaching genome-wide significance (P < 5 × 10 -8 ) in the 1KG+7K imputation output yet remaining undetected when the same sample set is imputed up to the TOPMed reference panel. In summary, the 1KG+7K panel demonstrates significant advantages in the discovery of trait-associated loci, particularly those influenced by low-frequency association signals., Competing Interests: Competing interests: The authors declare no competing interests. Kaoru Ito is an Editorial Board Member for Communications Biology, but was not involved in the editorial review of, nor the decision to publish this article., (© 2024. The Author(s).)

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., (© 2024. The Author(s).)

Human genetic variants are associated with many traits through largely unknown mechanisms. Here, combining approximately 260,000 Japanese study participants, a Japanese-specific genotype reference panel and statistical fine-mapping, we identified 4,423 significant loci across 63 quantitative traits, among which 601 were new, and 9,406 putatively causal variants. New associations included Japanese-specific coding, splicing and noncoding variants, exemplified by a damaging missense variant rs730881101 in TNNT2 associated with lower heart function and increased risk for heart failure (P = 1.4 × 10 -15 and odds ratio = 4.5, 95% confidence interval = 3.1-6.5). Putative causal noncoding variants were supported by state-of-art in silico functional assays and had comparable effect sizes to coding variants. A plausible example of new mechanisms of causal variants is an enrichment of causal variants in 3' untranslated regions (UTRs), including the Japanese-specific rs13306436 in IL6 associated with pro-inflammatory traits and protection against tuberculosis. We experimentally showed that transcripts with rs13306436 are resistant to mRNA degradation by regnase-1, an RNA-binding protein. Our study provides a list of fine-mapped causal variants to be tested for functionality and underscores the importance of sequencing, genotyping and association efforts in diverse populations., (© 2024. The Author(s).)

Importance: Vasospastic angina (VSA) is vasospasm of the coronary artery and is particularly prevalent in East Asian populations. However, the specific genetic architecture for VSA at genome-wide levels is not fully understood., Objective: To identify genetic factors associated with VSA., Design, Setting, and Participants: This was a case-control genome-wide association study of VSA. Data from Biobank Japan (BBJ; enrolled patients from 2002-2008 and 2013-2018) were used, and controls without coronary artery disease (CAD) were enrolled. Patients from the BBJ were genotyped using arrays or a set of arrays. Patients recruited between 2002 and 2005 were classified within the first dataset, and those recruited between 2006 and 2008 were classified within the second dataset. To replicate the genome-wide association study in the first and second datasets, VSA cases and control samples from the latest patients in the BBJ recruited between 2013 and 2018 were analyzed in a third dataset., Exposures: Single-nucleotide variants associated with VSA., Main Outcomes and Measures: Cases with VSA and controls without CAD., Results: A total of 5720 cases (mean [SD] age, 67 [10] years; 3672 male [64.2%]) and 153 864 controls (mean [SD] age, 62 [15] years; 77 362 male [50.3%]) in 3 datasets were included in this study. The variants at the RNF213 locus showed the strongest association with VSA across the 3 datasets (odds ratio [OR], 2.34; 95% CI, 1.99-2.74; P = 4.4 × 10-25). Additionally, rs112735431, an Asian-specific rare deleterious variant (p.Arg4810Lys) experimentally shown to be associated with reduced angiogenesis and a well-known causal risk for Moyamoya disease was the most promising candidate for a causal variant explaining the association. The effect size of rs112735431 on VSA was distinct from that of other CADs. Furthermore, homozygous carriers of rs112735431 showed an association with VSA characterized by a large effect estimate (OR, 18.34; 95% CI, 5.15-65.22; P = 7.0 × 10-6), deviating from the additive model (OR, 4.35; 95% CI, 1.18-16.05; P = .03). Stratified analyses revealed that rs112735431 exhibited a stronger association in males (χ21 = 7.24; P = .007) and a younger age group (OR, 3.06; 95% CI, 2.24-4.19), corresponding to the epidemiologic features of VSA. In the registry, carriers without CAD of the risk allele rs112735431 had a strikingly high mortality rate due to acute myocardial infarction during the follow-up period (hazard ratio, 2.71; 95% CI, 1.57-4.65; P = 3.3 × 10-4). As previously reported, a possible overlap between VSA and Moyamoya disease was not found., Conclusions and Relevance: Results of this study suggest that vascular cell dysfunction mediated by variants in the RNF213 locus may promote coronary vasospasm, and the presence of the risk allele could serve as a predictive factor for the prognosis.

Elevated blood pressure (BP) is major risk factor for cardiovascular diseases (CVD). Genome-wide association studies (GWAS) conducted predominantly in populations of European ancestry have identified >2,000 BP-associated loci, but other ancestries have been less well-studied. We conducted GWAS of systolic, diastolic, pulse, and mean arterial BP in 100,453 Chinese adults. We identified 128 non-overlapping loci associated with one or more BP traits, including 74 newly-reported associations. Despite strong genetic correlations between populations, we identified appreciably higher heritability and larger variant effect sizes in Chinese compared with European or Japanese ancestry populations. Using instruments derived from these GWAS, multivariable Mendelian randomisation demonstrated that BP traits contribute differently to the causal associations of BP with CVD. In particular, only pulse pressure was independently causally associated with carotid plaque. These findings reinforce the need for studies in diverse populations to understand the genetic determinants of BP traits and their roles in disease risk., (© 2024. The Author(s).)

Transcribed enhancer maps can reveal nuclear interactions underpinning each cell type and connect specific cell types to diseases. Using a 5' single-cell RNA sequencing approach, we defined transcription start sites of enhancer RNAs and other classes of coding and noncoding RNAs in human CD4 + T cells, revealing cellular heterogeneity and differentiation trajectories. Integration of these datasets with single-cell chromatin profiles showed that active enhancers with bidirectional RNA transcription are highly cell type-specific and that disease heritability is strongly enriched in these enhancers. The resulting cell type-resolved multimodal atlas of bidirectionally transcribed enhancers, which we linked with promoters using fine-scale chromatin contact maps, enabled us to systematically interpret genetic variants associated with a range of immune-mediated diseases.

The induction of tissue-specific vessels in in vitro living tissue systems remains challenging. Here, we directly differentiated human pluripotent stem cells into CD32b + putative liver sinusoidal progenitors (iLSEP) by dictating developmental pathways. By devising an inverted multilayered air-liquid interface (IMALI) culture, hepatic endoderm, septum mesenchyme, arterial and sinusoidal quadruple progenitors self-organized to generate and sustain hepatocyte-like cells neighbored by divergent endothelial subsets composed of CD32b low CD31 high , LYVE1 + STAB1 + CD32b high CD31 low THBD - vWF - , and LYVE1 - THBD + vWF + cells. Wnt2 mediated sinusoidal-to-hepatic intercellular crosstalk potentiates hepatocyte differentiation and branched endothelial network formation. Intravital imaging revealed iLSEP developed fully patent human vessels with functional sinusoid-like features. Organoid-derived hepatocyte- and sinusoid-derived coagulation factors enabled correction of in vitro clotting time with Factor V, VIII, IX, and XI deficient patients' plasma and rescued the severe bleeding phenotype in hemophilia A mice upon transplantation. Advanced organoid vascularization technology allows for interrogating key insights governing organ-specific vessel development, paving the way for coagulation disorder therapeutics.

Cerebral small vessel disease (cSVD) is a leading cause of stroke and dementia with no specific mechanism-based treatment. We used Mendelian randomization to combine a unique cerebrospinal fluid (CSF) and plasma pQTL resource with the latest European-ancestry GWAS of MRI-markers of cSVD (white matter hyperintensities, perivascular spaces). We describe a new biological fingerprint of 49 protein-cSVD associations, predominantly in the CSF. We implemented a multipronged follow-up, across fluids, platforms, and ancestries (Europeans and East-Asian), including testing associations of direct plasma protein measurements with MRI-cSVD. We highlight 16 proteins robustly associated in both CSF and plasma, with 24/4 proteins identified in CSF/plasma only. cSVD-proteins were enriched in extracellular matrix and immune response pathways, and in genes enriched in microglia and specific microglial states (integration with single-nucleus RNA sequencing). Immune-related proteins were associated with MRI-cSVD already at age twenty. Half of cSVD-proteins were associated with stroke, dementia, or both, and seven cSVD-proteins are targets for known drugs (used for other indications in directions compatible with beneficial therapeutic effects. This first cSVD proteogenomic signature opens new avenues for biomarker and therapeutic developments., Competing Interests: Competing interests C.C. has received research support from GSK and EISAI and is a member of the advisory board of Circular Genomics and owns stocks in this company. CC is part of the scientific advisory board for ADmit. B.P. serves on the Steering Committee of the Yale Open Data Project funded by Johnson & Johnson. P.M.M. has received an honourarium as Chair of the UKRI Medical Research Council Neuroscience and Mental Health Board until March 2024. He acknowledges consultancy fees from, Biogen, Sudo. Nimbus and GSK. He has received speakers’ honoraria from Sanofi and Redburn, and has received research or educational funds from Biogen, Merck, Bristol Meyers Squibb and Nimbus. J.W declares no commercial COI; various academic research grants and is CI for LACunar Intervention Trials. The authors declared no potential conflicts of interest with respect to research, authorship, and/or publication of this article. The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.

Background: Recently, various polygenic risk score (PRS)-based methods were developed to improve stroke prediction. However, current PRSs (including cross-ancestry PRS) poorly predict recurrent stroke. Here, we aimed to determine whether the best PRS for Japanese individuals can also predict stroke recurrence in this population by extensively comparing the methods and maximizing the predictive performance for stroke onset., Methods: We used data from the BioBank Japan (BBJ) 1 st cohort (n=179,938) to derive and optimize the PRSs using a 10-fold cross-validation. We integrated the optimized PRSs for multiple traits, such as vascular risk factors and stroke subtypes to generate a single PRS using the meta-scoring approach (metaGRS). We used an independent BBJ 2 nd cohort (n=41,929) as a test sample to evaluate the association of the metaGRS with stroke and recurrent stroke., Results: We analyzed recurrent stroke cases (n=174) and non-recurrent stroke controls (n=1,153) among subjects within the BBJ 2 nd cohort. After adjusting for known risk factors, metaGRS was associated with stroke recurrence (adjusted OR per SD 1.18 [95% CI: 1.00-1.39, p=0.044]), although no significant correlation was observed with the published PRSs. We administered three distinct tests to consider the potential index event bias; however, the outcomes derived from these examinations did not provide any significant indication of the influence of index event bias. The high metaGRS group without a history of hypertension had a higher risk of stroke recurrence than that of the low metaGRS group (adjusted OR 2.24 [95% CI: 1.07-4.66, p=0.032]). However, this association was weak in the hypertension group (adjusted OR 1.21 [95% CI: 0.69-2.13, p=0.50])., Conclusions: The metaGRS developed in a Japanese cohort predicted stroke recurrence in an independent cohort of patients. In particular, it predicted an increased risk of recurrence among stroke patients without hypertension. These findings provide clues for additional genetic risk stratification and help in developing personalized strategies for stroke recurrence prevention.

Objective: Adequate conditioning results from various physical, environmental, and psychological factors in sports activities. In this study, we aimed to clarify the concurrence and relevance of injuries, psychological problems, and sleep disturbance in university football (soccer) players. Biomechanical characteristics and risk factors for those injuries were also investigated., Methods: Overall, 1,084 university football players participated in this cross-sectional study using a web-based health history questionnaire. Assessments were made via the presence of injuries for physical problems, the Athlete Psychological Strain Questionnaire (APSQ) for psychological problems, and the Athlete Sleep Screening Questionnaire (ASSQ) for sleep difficulties., Results: There were 200 (19%) athletes with injuries; 413 (38%) of these were APSQ-positive, 246 (23%) were ASSQ-positive, and 468 (43%) athletes had no physical, psychological, or sleep problems. APSQ-positive respondents were significantly more likely to have a current injury than those who were APSQ-negative. No significant difference was observed between the number of ASSQ-positive and -negative respondents and current injuries. Those with sleep problems were significantly more likely to have psychological problems., Conclusion: More than 60% of university football players had one or more physical, psychological, or sleep problems. Players may have co-occurring problems, and it is important to address them for these athletes to reach peak performance.

Genetic variants in HTRA1 are associated with stroke risk. However, the mechanisms mediating this remain largely unknown, as does the full spectrum of phenotypes associated with genetic variation in HTRA1. Here we show that rare HTRA1 variants are linked to ischemic stroke in the UK Biobank and BioBank Japan. Integrating data from biochemical experiments, we next show that variants causing loss of protease function associated with ischemic stroke, coronary artery disease and skeletal traits in the UK Biobank and MyCode cohorts. Moreover, a common variant modulating circulating HTRA1 mRNA and protein levels enhances the risk of ischemic stroke and coronary artery disease while lowering the risk of migraine and macular dystrophy in genome-wide association study, UK Biobank, MyCode and BioBank Japan data. We found no interaction between proxied HTRA1 activity and levels. Our findings demonstrate the role of HTRA1 for cardiovascular diseases and identify two mechanisms as potential targets for therapeutic interventions., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

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&#95;3D&#95;lncRNA., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Agrawal et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

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.

Peptic ulcer disease (PUD) refers to acid-induced injury of the digestive tract, occurring mainly in the stomach (gastric ulcer (GU)) or duodenum (duodenal ulcer (DU)). In the present study, we conducted a large-scale, cross-ancestry meta-analysis of PUD combining genome-wide association studies with Japanese and European studies (52,032 cases and 905,344 controls), and discovered 25 new loci highly concordant across ancestries. An examination of GU and DU genetic architecture demonstrated that GUs shared the same risk loci as DUs, although with smaller genetic effect sizes and higher polygenicity than DUs, indicating higher heterogeneity of GUs. Helicobacter pylori (HP)-stratified analysis found an HP-related host genetic locus. Integrative analyses using bulk and single-cell transcriptome profiles highlighted the genetic factors of PUD being enriched in the highly expressed genes in stomach tissues, especially in somatostatin-producing D cells. Our results provide genetic evidence that gastrointestinal cell differentiations and hormone regulations are critical in PUD etiology., (© 2023. The Author(s).)

Large-scale imputation reference panels are currently available and have contributed to efficient genome-wide association studies through genotype imputation. However, whether large-size multi-ancestry or small-size population-specific reference panels are the optimal choices for under-represented populations continues to be debated. We imputed genotypes of East Asian (180k Japanese) subjects using the Trans-Omics for Precision Medicine reference panel and found that the standard imputation quality metric (Rsq) overestimated dosage r2 (squared correlation between imputed dosage and true genotype) particularly in marginal-quality bins. Variance component analysis of Rsq revealed that the increased imputed-genotype certainty (dosages closer to 0, 1 or 2) caused upward bias, indicating some systemic bias in the imputation. Through systematic simulations using different template switching rates (θ value) in the hidden Markov model, we revealed that the lower θ value increased the imputed-genotype certainty and Rsq; however, dosage r2 was insensitive to the θ value, thereby causing a deviation. In simulated reference panels with different sizes and ancestral diversities, the θ value estimates from Minimac decreased with the size of a single ancestry and increased with the ancestral diversity. Thus, Rsq could be deviated from dosage r2 for a subpopulation in the multi-ancestry panel, and the deviation represents different imputed-dosage distributions. Finally, despite the impact of the θ value, distant ancestries in the reference panel contributed only a few additional variants passing a predefined Rsq threshold. We conclude that the θ value substantially impacts the imputed dosage and the imputation quality metric value., (© The Author(s) 2024. Published by Oxford University Press.)

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., Competing Interests: COMPETING INTERESTS C.D.A. has received sponsored research support from Bayer AG, and has consulted for ApoPharma, unrelated to the content of this manuscript. All other authors declare no competing interests.

The link between aerosol dynamics and viral exposure risk is not fully understood, particularly during movement and face-to-face interactions. To investigate this, we employed Particle Trace Velocimetry with a laser sheet and a high-speed camera to measure microparticles from a human mannequin's mouth. The average peak time in the non-ventilated condition (expiratory volume, 30 L; passing speed, 5 km/h) was 1.33 s (standard deviation = 0.32 s), while that in the ventilated condition was 1.38 s (standard deviation = 0.35 s). Our results showed that the peak of viral exposure risk was within 5 s during face-to-face encounters under both ventilated and non-ventilated conditions. Moreover, the risk of viral exposure greatly decreased in ventilated conditions compared to non-ventilated conditions. Based on these findings, considering a risk mitigation strategy for the duration of 5 s during face-to-face encounters is expected to significantly reduce the risk of virus exposure in airborne transmission., (© 2023. Springer Nature Limited.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Otomo, Khanshour, Koido, Takeda, Momozawa, Kubo, Kamatani, Herring, Ogura, Takahashi, Minami, Uno, Kawakami, Ito, Sato, Watanabe, Kaito, Yanagida, Taneichi, Harimaya, Taniguchi, Shigematsu, Iida, Demura, Sugawara, Fujita, Yagi, Okada, Hosogane, Kono, Nakamura, Chiba, Kotani, Sakuma, Akazawa, Suzuki, Nishida, Kakutani, Tsuji, Sudo, Iwata, Inami, Wise, Kochi, Matsumoto, Ikegawa, Watanabe and Terao.)

Gene transcription is regulated through complex mechanisms involving non-coding RNAs (ncRNAs). As the transcription of ncRNAs, especially of enhancer RNAs, is often low and cell type specific, how the levels of RNA transcription depend on genotype remains largely unexplored. Here we report the development and utility of a machine-learning model (MENTR) that reliably links genome sequence and ncRNA expression at the cell type level. Effects on ncRNA transcription predicted by the model were concordant with estimates from published studies in a cell-type-dependent manner, regardless of allele frequency and genetic linkage. Among 41,223 variants from genome-wide association studies, the model identified 7,775 enhancer RNAs and 3,548 long ncRNAs causally associated with complex traits across 348 major human primary cells and tissues, such as rare variants plausibly altering the transcription of enhancer RNAs to influence the risks of Crohn's disease and asthma. The model may aid the discovery of causal variants and the generation of testable hypotheses for biological mechanisms driving complex traits., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

Ossification of the posterior longitudinal ligament of the spine (OPLL) is a common intractable disease that causes spinal stenosis and myelopathy. We have previously conducted genome-wide association studies for OPLL and identified 14 significant loci, but their biological implications remain mostly unclear. Here, we examined the 12p11.22 locus and identified a variant in the 5' UTR of a novel isoform of CCDC91 that was associated with OPLL. Using machine learning prediction models, we determined that higher expression of the novel CCDC91 isoform was associated with the G allele of rs35098487. The risk allele of rs35098487 showed higher affinity in the binding of nuclear proteins and transcription activity. Knockdown and overexpression of the CCDC91 isoform in mesenchymal stem cells and MG-63 cells showed paralleled expression of osteogenic genes, including RUNX2, the master transcription factor of osteogenic differentiation. The CCDC91 isoform directly interacted with MIR890, which bound to RUNX2 and decreased RUNX2 expression. Our findings suggest that the CCDC91 isoform acts as a competitive endogenous RNA by sponging MIR890 to increase RUNX2 expression., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Background: Recently, common genetic risk factors for intracranial aneurysm (IA) and aneurysmal subarachnoid hemorrhage (ASAH) were found to explain a large amount of disease heritability and therefore have potential to be used for genetic risk prediction. We constructed a genetic risk score to (1) predict ASAH incidence and IA presence (combined set of unruptured IA and ASAH) and (2) assess its association with patient characteristics., Methods: A genetic risk score incorporating genetic association data for IA and 17 traits related to IA (so-called metaGRS) was created using 1161 IA cases and 407 392 controls from the UK Biobank population study. The metaGRS was validated in combination with risk factors blood pressure, sex, and smoking in 828 IA cases and 68 568 controls from the Nordic HUNT population study. Furthermore, we assessed association between the metaGRS and patient characteristics in a cohort of 5560 IA patients., Results: Per SD increase of metaGRS, the hazard ratio for ASAH incidence was 1.34 (95% CI, 1.20-1.51) and the odds ratio for IA presence 1.09 (95% CI, 1.01-1.18). Upon including the metaGRS on top of clinical risk factors, the concordance index to predict ASAH hazard increased from 0.63 (95% CI, 0.59-0.67) to 0.65 (95% CI, 0.62-0.69), while prediction of IA presence did not improve. The metaGRS was statistically significantly associated with age at ASAH (β=-4.82×10 -3 per year [95% CI, -6.49×10 -3 to -3.14×10 -3 ]; P =1.82×10 -8 ), and location of IA at the internal carotid artery (odds ratio=0.92 [95% CI, 0.86-0.98]; P =0.0041)., Conclusions: The metaGRS was predictive of ASAH incidence, although with limited added value over clinical risk factors. The metaGRS was not predictive of IA presence. Therefore, we do not recommend using this metaGRS in daily clinical care. Genetic risk does partly explain the clinical heterogeneity of IA warranting prioritization of clinical heterogeneity in future genetic prediction studies of IA and ASAH.

Rheumatoid arthritis (RA) is a highly heritable complex disease with unknown etiology. Multi-ancestry genetic research of RA promises to improve power to detect genetic signals, fine-mapping resolution and performances of polygenic risk scores (PRS). Here, we present a large-scale genome-wide association study (GWAS) of RA, which includes 276,020 samples from five ancestral groups. We conducted a multi-ancestry meta-analysis and identified 124 loci (P < 5 × 10 -8 ), of which 34 are novel. Candidate genes at the novel loci suggest essential roles of the immune system (for example, TNIP2 and TNFRSF11A) and joint tissues (for example, WISP1) in RA etiology. Multi-ancestry fine-mapping identified putatively causal variants with biological insights (for example, LEF1). Moreover, PRS based on multi-ancestry GWAS outperformed PRS based on single-ancestry GWAS and had comparable performance between populations of European and East Asian ancestries. Our study provides several insights into the etiology of RA and improves the genetic predictability of RA., (© 2022. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Previous genome-wide association studies (GWASs) of stroke - the second leading cause of death worldwide - were conducted predominantly in populations of European ancestry 1,2 . Here, in cross-ancestry GWAS meta-analyses of 110,182 patients who have had a stroke (five ancestries, 33% non-European) and 1,503,898 control individuals, we identify association signals for stroke and its subtypes at 89 (61 new) independent loci: 60 in primary inverse-variance-weighted analyses and 29 in secondary meta-regression and multitrait analyses. On the basis of internal cross-ancestry validation and an independent follow-up in 89,084 additional cases of stroke (30% non-European) and 1,013,843 control individuals, 87% of the primary stroke risk loci and 60% of the secondary stroke risk loci were replicated (P < 0.05). Effect sizes were highly correlated across ancestries. Cross-ancestry fine-mapping, in silico mutagenesis analysis 3 , and transcriptome-wide and proteome-wide association analyses revealed putative causal genes (such as SH3PXD2A and FURIN) and variants (such as at GRK5 and NOS3). Using a three-pronged approach 4 , we provide genetic evidence for putative drug effects, highlighting F11, KLKB1, PROC, GP1BA, LAMC2 and VCAM1 as possible targets, with drugs already under investigation for stroke for F11 and PROC. A polygenic score integrating cross-ancestry and ancestry-specific stroke GWASs with vascular-risk factor GWASs (integrative polygenic scores) strongly predicted ischaemic stroke in populations of European, East Asian and African ancestry 5 . Stroke genetic risk scores were predictive of ischaemic stroke independent of clinical risk factors in 52,600 clinical-trial participants with cardiometabolic disease. Our results provide insights to inform biology, reveal potential drug targets and derive genetic risk prediction tools across ancestries., (© 2022. The Author(s).)

Objective: Genome-wide association studies (GWAS) have identified >100 risk loci for systemic lupus erythematosus (SLE), but the disease genes at most loci remain unclear, hampering translation of these genetic discoveries. We aimed to prioritise genes underlying the 110 SLE loci that were identified in the latest East Asian GWAS meta-analysis., Methods: We built gene expression predictive models in blood B cells, CD4 + and CD8 + T cells, monocytes, natural killer cells and peripheral blood cells of 105 Japanese individuals. We performed a transcriptome-wide association study (TWAS) using data from the latest genome-wide association meta-analysis of 208 370 East Asians and searched for candidate genes using TWAS and three data-driven computational approaches., Results: TWAS identified 171 genes for SLE (p<1.0×10 -5 ); 114 (66.7%) showed significance only in a single cell type; 127 (74.3%) were in SLE GWAS loci. TWAS identified a strong association between CD83 and SLE (p<7.7×10 -8 ). Meta-analysis of genetic associations in the existing 208 370 East Asian and additional 1498 cases and 3330 controls found a novel single-variant association at rs72836542 (OR=1.11, p=4.5×10 -9 ) around CD83 . For the 110 SLE loci, we identified 276 gene candidates, including 104 genes at recently-identified SLE novel loci. We demonstrated in vitro that putative causal variant rs61759532 exhibited an allele-specific regulatory effect on ACAP1 , and that presence of the SLE risk allele decreased ACAP1 expression., Conclusions: Cell-level TWAS in six types of immune cells complemented SLE gene discovery and guided the identification of novel genetic associations. The gene findings shed biological insights into SLE genetic associations., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2022. Re-use permitted under CC BY. Published by BMJ.)

The unfolded protein response (UPR) is an adaptive stress response pathway that is essential for cancer cell survival under endoplasmic reticulum stress such as during glucose starvation. In this study, we identified spautin-1, an autophagy inhibitor that suppresses ubiquitin-specific peptidase 10 (USP10) and USP13, as a novel UPR inhibitor under glucose starvation conditions. Spautin-1 prevented the induction of UPR-associated proteins, including glucose-regulated protein 78, activating transcription factor 4, and a splicing variant of x-box-binding protein-1, and showed preferential cytotoxicity in glucose-starved cancer cells. However, USP10 and USP13 silencing and treatment with other autophagy inhibitors failed to result in UPR inhibition and preferential cytotoxicity during glucose starvation. Using transcriptome and chemosensitivity-based COMPARE analyses, we identified a similarity between spautin-1 and mitochondrial complex I inhibitors and found that spautin-1 suppressed the activity of complex I extracted from isolated mitochondria. Our results indicated that spautin-1 may represent an attractive mitochondria-targeted seed compound that inhibits the UPR and cancer cell survival during glucose starvation., (© 2022. The Author(s).)

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) ORF6 is an antagonist of interferon (IFN)-mediated antiviral signaling, achieved through the prevention of STAT1 nuclear localization. However, the exact mechanism through which ORF6 prevents STAT1 nuclear trafficking remains unclear. Herein, we demonstrate that ORF6 directly binds to STAT1 with or without IFN stimulation, resulting in the nuclear exclusion of STAT1. ORF6 also recognizes importin α subtypes with different modes, in particular, high affinity to importin α1 but a low affinity to importin α5. Although ORF6 potentially disrupts the importin α/importin β1-mediated nuclear transport, thereby suppressing the nuclear translocation of the other classical nuclear localization signal-containing cargo proteins, the inhibitory effect of ORF6 is modest when compared with that of STAT1. The results indicate that the drastic nuclear exclusion of STAT1 is attributed to the specific binding with ORF6, which is a distinct strategy for the importin α1-mediated pathway. Combined with the results from a newly-produced replicon system and a hamster model, we conclude that SARS-CoV-2 ORF6 acts as a virulence factor via regulation of nucleocytoplasmic trafficking to accelerate viral replication, resulting in disease progression., (© 2022. The Author(s).)

Osteonecrosis of the femoral head (ONFH) involves necrosis of bone and bone marrow of the femoral head caused by ischemia with unknown etiology. Previous genetic studies on ONFH failed to produce consistent results, presumably because ONFH has various causes with different genetic backgrounds and the underlying diseases confounded the associations. Steroid-associated ONFH (S-ONFH) accounts for one-half of all ONFH, and systemic lupus erythematosus (SLE) is a representative disease underlying S-ONFH. We performed a genome-wide association study (GWAS) to identify genetic risk factors for S-ONFH in patients with SLE. We conducted a two-staged GWAS on 636 SLE patients with S-ONFH and 95 588 non-SLE controls. Among the novel loci identified, we determined S-ONFH-specific loci by comparing allele frequencies between SLE patients without S-ONFH and non-SLE controls. We also used Korean datasets comprising 148 S-ONFH cases and 37 015 controls to assess overall significance. We evaluated the functional annotations of significant variants by in silico analyses. The Japanese GWAS identified 4 significant loci together with 12 known SLE susceptibility loci. The four significant variants showed comparable effect sizes on S-ONFH compared with SLE controls and non-SLE controls. Three of the four loci, MIR4293/MIR1265 [odds ratio (OR) = 1.99, P-value = 1.1 × 10-9)], TRIM49/NAALAD2 (OR = 1.65, P-value = 4.8 × 10-8) and MYO16 (OR = 3.91, P-value = 4.9 × 10-10), showed significant associations in the meta-analysis with Korean datasets. Bioinformatics analyses identified MIR4293, NAALAD2 and MYO16 as candidate causal genes. MIR4293 regulates a PPARG-related adipogenesis pathway relevant to S-ONFH. We identified three novel susceptibility loci for S-ONFH in SLE., (© The Author(s) 2021. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Spatial transcriptomics is an emerging technology requiring costly reagents and considerable skills, limiting the identification of transcriptional markers related to histology. Here, we show that predicted spatial gene-expression in unmeasured regions and tissues can enhance biologists' histological interpretations. We developed the Deep learning model for Spatial gene Clusters and Expression, DeepSpaCE, and confirmed its performance using the spatial-transcriptome profiles and immunohistochemistry images of consecutive human breast cancer tissue sections. For example, the predicted expression patterns of SPARC, an invasion marker, highlighted a small tumor-invasion region difficult to identify using raw spatial transcriptome data alone because of a lack of measurements. We further developed semi-supervised DeepSpaCE using unlabeled histology images and increased the imputation accuracy of consecutive sections, enhancing applicability for a small sample size. Our method enables users to derive hidden histological characters via spatial transcriptome and gene annotations, leading to accelerated biological discoveries without additional experiments., (© 2022. The Author(s).)

No genome-wide association studies (GWAS) were reported for colorectal polyps and the overlap in polygenic backgrounds conferring risk of colorectal cancer and polyps remains unclear. We performed GWAS on subjects with colorectal polyps using the BioBank Japan data with 4447 cases and 157,226 controls. We evaluated genetic correlations between colorectal polyps and cancer, and effects on colorectal polyps of single nucleotide polymorphisms (SNPs) known to be associated with colorectal cancer. We identified CUX2, a known genetic locus to colorectal cancer, as a susceptibility locus to colorectal polyps (p value = 1.1 × 10 -15 ). Subsequent fine-mapping analysis indicated that rs11065828 in CUX2 is the causal variant for colorectal polyps. We found that known colorectal cancer-susceptible SNPs were also associated with colorectal polyps. The genetic correlation between colorectal cancer and polyps is very high (r = 0.98 and p value = 0.0006). We additionally identified 14 significant loci of colorectal polyps and three significant loci of colorectal cancer by applying the multi-trait analysis of GWAS of colorectal cancer and colorectal polyps. We showed very similar germline polygenic features, which gives us the additional insight into potential cancers at polygenic levels for patients with polyps who are followed up at outpatients' clinic; thus, close observation and polypectomy is critical to prevent colorectal cancers., (© 2021. The Author(s), under exclusive licence to The Japan Society of Human Genetics.)

Background: Idiopathic pulmonary fibrosis (IPF) is a progressive, fatal fibrotic interstitial lung disease. Few circulating biomarkers have been identified to have causal effects on IPF., Methods: To identify candidate IPF-influencing circulating proteins, we undertook an efficient screen of circulating proteins by applying a two-sample Mendelian randomisation (MR) approach with existing publicly available data. For instruments, we used genetic determinants of circulating proteins which reside cis to the encoded gene ( cis -single nucleotide polymorphisms (SNPs)), identified by two genome-wide association studies (GWASs) in European individuals (3301 and 3200 subjects). We then applied MR methods to test if the levels of these circulating proteins influenced IPF susceptibility in the largest IPF GWAS (2668 cases and 8591 controls). We validated the MR results using colocalisation analyses to ensure that both the circulating proteins and IPF shared a common genetic signal., Results: MR analyses of 834 proteins found that a 1 sd increase in circulating galactoside 3(4)-l-fucosyltransferase (FUT3) and α-(1,3)-fucosyltransferase 5 (FUT5) was associated with a reduced risk of IPF (OR 0.81, 95% CI 0.74-0.88; p=6.3×10 -7 and OR 0.76, 95% CI 0.68-0.86; p=1.1×10 -5 , respectively). Sensitivity analyses including multiple cis -SNPs provided similar estimates both for FUT3 (inverse variance weighted (IVW) OR 0.84, 95% CI 0.78-0.91; p=9.8×10 -6 and MR-Egger OR 0.69, 95% CI 0.50-0.97; p=0.03) and FUT5 (IVW OR 0.84, 95% CI 0.77-0.92; p=1.4×10 -4 and MR-Egger OR 0.59, 95% CI 0.38-0.90; p=0.01). FUT3 and FUT5 signals colocalised with IPF signals, with posterior probabilities of a shared genetic signal of 99.9% and 97.7%, respectively. Further transcriptomic investigations supported the protective effects of FUT3 for IPF., Conclusions: An efficient MR scan of 834 circulating proteins provided evidence that genetically increased circulating FUT3 level is associated with reduced risk of IPF., Competing Interests: Conflict of interest: T. Nakanishi has nothing to disclose. Conflict of interest: A. Cerani has nothing to disclose. Conflict of interest: V. Forgetta has nothing to disclose. Conflict of interest: S. Zhou has nothing to disclose. Conflict of interest: R.J. Allen has nothing to disclose. Conflict of interest: O.C. Leavy has nothing to disclose. Conflict of interest: M. Koido has nothing to disclose. Conflict of interest: D. Assayag reports grants and personal fees for consultancy from Boehringer Ingelheim Canada, personal fees for consultancy from Hoffman-LaRoche Canada, personal fees for lectures from AstraZeneca Canada, outside the submitted work. Conflict of interest: R.G. Jenkins reports grants from AstraZeneca, Biogen, Galecto and GlaxoSmithKline, personal fees from Boehringer Ingelheim, Daewoong, Galapagos, Heptares, Promedior and Roche, nonfinancial support from NuMedii and Redx, grants and personal fees from Pliant, other (trustee) from Action for Pulmonary Fibrosis, outside the submitted work. Conflict of interest: L.V. Wain reports grants from GlaxoSmithKline, outside the submitted work. Conflict of interest: I.V. Yang reports grants from the NIH and personal fees from Eleven P15 related to research in pulmonary fibrosis, outside the submitted work; and has a patent “Circulating biomarkers of preclinical pulmonary fibrosis” pending. Conflict of interest: G.M. Lathrop has nothing to disclose. Conflict of interest: P.J. Wolters reports grants and personal fees from Boehringer Ingelheim and Roche/Genentech, personal fees from Gossamer Bio, Blade Therapeutics and Pliant, outside the submitted work. Conflict of interest: D.A. Schwartz reports grants from the NIH-NHLBI (R38-HL143511, T32-HL007085, UG3/UH3-HL151865, R01-HL149836, P01-HL0928701, UH2/3-HL123442, X01-HL134585 and R25-ES025476) and DOD Focused Program Grant W81XWH-17-1-0597, during the conduct of the study; personal fees from Eleven P15, outside the submitted work; and has a patent “Compositions and methods of treating or preventing fibrotic diseases” pending, a patent “Biomarkers for the diagnosis and treatment of fibrotic lung disease” pending and a patent “Methods and compositions for risk prediction, diagnosis, prognosis, and treatment of pulmonary disorders” issued. Conflict of interest: J.B. Richards has served as an advisor to GlaxoSmithKline and Deerfield Capital., (Copyright ©The authors 2022.)

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., (Copyright © 2021 American Academy of Allergy, Asthma & Immunology. Published by Elsevier Inc. All rights reserved.)

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)., (© 2021 American Society for Bone and Mineral Research (ASBMR).)

Background: The underlying pathology of inguinal hernia is still not fully known; thus, further investigations of genetic backgrounds is needed. Here, we aimed to identify genetic factors attributing to inguinal hernias and explore the polygenic architecture of which some components are population-specific, while others are more common among populations., Methods: We performed a genome-wide association study (GWAS) on subjects with inguinal hernias using BioBank Japan (BBJ) data with 1,983 cases and 172,507 controls, followed by a trans-ethnic meta-analysis with UK Biobank (UKBB) data. We performed downstream analyses in order to identify the mechanisms underlying inguinal hernias supported by genetic findings., Findings: We identified a locus closest to ELN, which encodes elastin, at the GWAS significant level. The trans-ethnic meta-analysis revealed 23 additional significant loci, including five loci newly identified not significant in BBJ or UKBB GWAS: TGFB2, RNA5SP214/VGLL2, LOC646588, HMCN2, and ATP5F1CP1/CDKN3. Downstream analyses revealed the overlap of GWAS significant signals in extracellular components, including elastin fiber formation. We also found a highly shared polygenic architecture across different populations (trans-ethnic genetic-effect correlation = 0•77, standard error = 0•26) and population-specific lead variants in ELN, indicating the critical role of elastin in inguinal hernias., Interpretation: We identified a significant locus of the ELN gene in the Japanese population and five additional loci across different populations. Downstream analyses revealed highly shared genetic architectures across populations and highlighted the important roles of extracellular components in the development of inguinal hernias. These findings deepen our understanding of the mechanisms underlying inguinal hernia., Funding: The Japan Agency for Medical Research and Development (AMED) (Grant Number: JP19km0605001)., Competing Interests: Declaration of Competing Interest The authors declare no conflict of interest., (Copyright © 2021 The Authors. Published by Elsevier B.V. All rights reserved.)

Objective: Systemic lupus erythematosus (SLE), an autoimmune disorder, has been associated with nearly 100 susceptibility loci. Nevertheless, these loci only partially explain SLE heritability and their putative causal variants are rarely prioritised, which make challenging to elucidate disease biology. To detect new SLE loci and causal variants, we performed the largest genome-wide meta-analysis for SLE in East Asian populations., Methods: We newly genotyped 10 029 SLE cases and 180 167 controls and subsequently meta-analysed them jointly with 3348 SLE cases and 14 826 controls from published studies in East Asians. We further applied a Bayesian statistical approach to localise the putative causal variants for SLE associations., Results: We identified 113 genetic regions including 46 novel loci at genome-wide significance (p<5×10 -8 ). Conditional analysis detected 233 association signals within these loci, which suggest widespread allelic heterogeneity. We detected genome-wide associations at six new missense variants. Bayesian statistical fine-mapping analysis prioritised the putative causal variants to a small set of variants (95% credible set size ≤10) for 28 association signals. We identified 110 putative causal variants with posterior probabilities ≥0.1 for 57 SLE loci, among which we prioritised 10 most likely putative causal variants (posterior probability ≥0.8). Linkage disequilibrium score regression detected genetic correlations for SLE with albumin/globulin ratio (r g =-0.242) and non-albumin protein (r g =0.238)., Conclusion: This study reiterates the power of large-scale genome-wide meta-analysis for novel genetic discovery. These findings shed light on genetic and biological understandings of SLE., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2021. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.)

Poor trans-ancestry portability of polygenic risk scores is a consequence of Eurocentric genetic studies and limited knowledge of shared causal variants. Leveraging regulatory annotations may improve portability by prioritizing functional over tagging variants. We constructed a resource of 707 cell-type-specific IMPACT regulatory annotations by aggregating 5,345 epigenetic datasets to predict binding patterns of 142 transcription factors across 245 cell types. We then partitioned the common SNP heritability of 111 genome-wide association study summary statistics of European (average n ≈ 189,000) and East Asian (average n ≈ 157,000) origin. IMPACT annotations captured consistent SNP heritability between populations, suggesting prioritization of shared functional variants. Variant prioritization using IMPACT resulted in increased trans-ancestry portability of polygenic risk scores from Europeans to East Asians across all 21 phenotypes analyzed (49.9% mean relative increase in R 2 ). Our study identifies a crucial role for functional annotations such as IMPACT to improve the trans-ancestry portability of genetic data.

Rupture of an intracranial aneurysm leads to subarachnoid hemorrhage, a severe type of stroke. To discover new risk loci and the genetic architecture of intracranial aneurysms, we performed a cross-ancestry, genome-wide association study in 10,754 cases and 306,882 controls of European and East Asian ancestry. We discovered 17 risk loci, 11 of which are new. We reveal a polygenic architecture and explain over half of the disease heritability. We show a high genetic correlation between ruptured and unruptured intracranial aneurysms. We also find a suggestive role for endothelial cells by using gene mapping and heritability enrichment. Drug-target enrichment shows pleiotropy between intracranial aneurysms and antiepileptic and sex hormone drugs, providing insights into intracranial aneurysm pathophysiology. Finally, genetic risks for smoking and high blood pressure, the two main clinical risk factors, play important roles in intracranial aneurysm risk, and drive most of the genetic correlation between intracranial aneurysms and other cerebrovascular traits.

We aimed to identify attributing factors to the interindividual variabilities of the infusion rates in unfractionated heparin therapy. We included patients who required unfractionated heparin therapy to achieve the target APTT after cardiac surgery between May 2014 and February 2018. Fifty-nine patients were included, of whom 8 underwent Blalock-Taussig shunt; 27, Glenn procedure; 19, Fontan procedure; 3, mechanical valve replacement; and 2, Rastelli procedure. Previously reported variables that influenced the response to unfractionated heparin treatment were initially compared, which included age; weight; sex; type of surgery; platelet count; fibrinogen, antithrombin III, total protein, albumin, alanine transaminase, and creatinine levels; and use of fresh frozen plasma. The type of surgical procedure was found to be significantly associated with the differences in heparin infusion rate (P = 0.00073). Subsequently, the variance explained by these factors was estimated through a selection based on the minimum Akaike information criterion value; models constructed by various combinations of the surgery types were compared. The model including the Blalock-Taussig shunt, Glenn procedure, and mechanical valve replacement showed the highest summed variance explained (29.1%). More than 70% of the interindividual variability in initial heparin maintenance dosing was unexplained.

Drug-induced liver injury (DILI) is a leading cause of termination in drug development programs and removal of drugs from the market; this is partially due to the inability to identify patients who are at risk 1 . In this study, we developed a polygenic risk score (PRS) for DILI by aggregating effects of numerous genome-wide loci identified from previous large-scale genome-wide association studies 2 . The PRS predicted the susceptibility to DILI in patients treated with fasiglifam, amoxicillin-clavulanate or flucloxacillin and in primary hepatocytes and stem cell-derived organoids from multiple donors treated with over ten different drugs. Pathway analysis highlighted processes previously implicated in DILI, including unfolded protein responses and oxidative stress. In silico screening identified compounds that elicit transcriptomic signatures present in hepatocytes from individuals with elevated PRS, supporting mechanistic links and suggesting a novel screen for safety of new drug candidates. This genetic-, cellular-, organoid- and human-scale evidence underscored the polygenic architecture underlying DILI vulnerability at the level of hepatocytes, thus facilitating future mechanistic studies. Moreover, the proposed 'polygenicity-in-a-dish' strategy might potentially inform designs of safer, more efficient and robust clinical trials.

CYP2C9*3 and HLA-B alleles are reportedly associated with phenytoin-induced eruption in some East Asian populations; however, this finding is not readily applicable to the Japanese population. Thus, we aimed to investigate the risk alleles using samples and data from BioBank Japan. A total of 747 patients (24 cases and 723 tolerant controls) were selected for analysis. Case-control association studies were conducted, using CYP2C9*3, CYP2C9*27, CYP2C19*2, CYP2C19*3, and HLA-B allele genotype data. CYP2C9*3 carrier status was significantly associated with phenytoin-induced eruption (P = 0.0022, odds ratio 7.05, 95% confidence interval, 2.44-20.4). HLA-B*51:01 showed the most prominent association (P = 0.010, odds ratio 3.19, 95% confidence interval, 1.37-7.48). Including both of these features improved predictive performance, measured as area under the receiver operating characteristic curve, by 10%. CYP2C9*3 and HLA-B*51:01 allele carrier statuses are significantly associated with phenytoin-induced eruption; thus, checking this carrier status before prescription would decrease the incidence of phenytoin-induced eruption in clinical practice., (© 2019 The Authors Clinical Pharmacology & Therapeutics © 2019 American Society for Clinical Pharmacology and Therapeutics.)