Large-scale genome-wide association study in a Japanese population identifies novel susceptibility loci across different diseases

Associations between high-altitude adaptation and risk of cardiovascular diseases: a bidirectional Mendelian randomization study

High-altitude adaptation (HAA) was reported to be significantly associated with reduced risks for multiple cardiovascular diseases (CVDs). However, the causality and direction of the associations are largely uncharacterized. We aimed to examine the potential causal relationships between HAA and six types of CVD, including coronary artery disease (CAD), cerebral aneurysm, ischemic stroke, peripheral artery disease, arrhythmia and atrial fibrillation. We obtained the summary data from largest available genome-wide association study of HAA and six types of CVD. Two-sample bidirectional Mendelian randomization (MR) analyses were performed to infer the causality between them. In the sensitivity analyses, MR-Egger regression analyses and MR-Pleiotropy RESidual Sum and Outlier (MR-PRESSO) global analyses were used to assess the pleiotropic effects; Cochran's Q tests were used to test the heterogeneity by inverse variance-weighted (IVW) and MR-Egger methods; and the leave-one-out analyses were used to examine whether some single nucleotide polymorphisms (SNPs) could influence the results independently. The MR main analyses showed that the genetically instrumented HAA was significantly causally associated with the reduced risks of CAD (odds ratio [OR] = 0.029; 95% confidence interval [CI] = 0.004-0.234; P = 8.64 × 10-4). In contrast, there was no statistically significant relationship between CVDs and HAA. Our findings provide evidence for the causal effects of HAA on the reduced risks of CAD. However, there is no causality of CVDs on HAA. These findings might be helpful in developing the prevention and intervention strategies for CAD.

Multi-ancestry meta-analysis identifies 5 novel loci for ischemic stroke and reveals heterogeneity of effects between sexes and ancestries

Stroke is the second leading cause of death and disability worldwide. Stroke prevalence varies by sex and ancestry, possibly due to genetic heterogeneity between subgroups. We performed a genome-wide meta-analysis of 16 biobanks across multiple ancestries to study the genetics of ischemic stroke (60,176 cases, 1,310,725 controls) as part of the Global Biobank Meta-analysis Initiative (GBMI) and further combined the results with previously published MegaStroke. Five novel loci for ischemic stroke (LAMC1, CALCRL, PLSCR1, CDKN1A, and SWAP70) were identified after replication in four additional datasets. One previously reported locus showed significant ancestry heterogeneity (ABO), and one showed significant sex heterogeneity (ALDH2). The ALDH2 association was male specific (males p = 1.67e-24, females p = 0.126) and was additionally observed only in the East Asian ancestry (male) samples. These findings emphasize the need for more diverse datasets with large sample sizes to further understand the genetic predisposition of stroke in different ancestry and sex groups.

Association Between Telomere Length and Risk of Lung Cancer in an Asian Population: A Mendelian Randomization Study

BACKGROUND

Several traditional observational studies and Mendelian randomization (MR) studies have indicated an association between leukocyte telomere length (LTL) and the risk of lung cancer in the European population. However, the results in the Asian population are still unclear. The objective was to reveal the genetic causal association between LTL and the risk of lung cancer in the Asian population.

METHODS

We conducted a two-sample MR analysis using summary statistics. Instrumental variables (IVs) were obtained from the genome-wide association studies (GWAS) of LTL (n = 23,096) and lung cancer (n = 212,453) of Asian ancestry. We applied the random-effects inverse-variance weighted (IVW) model as the main method. As well, several other models were performed as complementary methods to assess the impact of potential MR assumption violations, including MR-Egger regression, weighted median, and weighted mode models.

RESULTS

We included eight single-nucleotide polymorphisms (SNPs) as IVs for LTL and found that LTL was significantly associated with the risk of lung cancer in the IVW model (odds ratio (OR): 1.60; 95% confidence interval (CI): 1.31 - 1.97; P = 5.96 × 10-6), which was in line with the results in the weighted median and weighted mode models. However, the relationship was not statistically significant in the MR-Egger regression model (OR: 1.44; 95% CI: 0.92 - 2.26; P = 0.160). Sensitivity analyses indicated the robustness of the results.

CONCLUSIONS

This two-sample MR study confirmed that longer telomere length significantly increased the risk of lung cancer in the Asian population, which was in accord with findings in the Western population.

A Large-Scale Exome-Wide Association Study Identifies Novel Germline Mutations in Lung Cancer

Rationale: Genome-wide association studies have identified common variants of lung cancer. However, the contribution of rare exome-wide variants, especially protein-coding variants, to cancers remains largely unexplored. Objectives: To evaluate the role of human exomes in genetic predisposition to lung cancer. Methods: We performed exome-wide association studies to detect the association of exomes with lung cancer in 30,312 patients and 652,902 control subjects. A scalable and accurate implementation of a generalized mixed model was used to detect the association signals for loss-of-function, missense, and synonymous variants and gene-level sets. Furthermore, we performed association and Bayesian colocalization analyses to evaluate their relationships with intermediate exposures. Measurements and Main Results: We systematically analyzed 216,739 single-nucleotide variants in the human exome. The loss-of-function variants exhibited the most notable effects on lung cancer risk. We identified four novel variants, including two missense variants (rs202197044TET3 [Pmeta (P values of meta-analysis) = 3.60 × 10-8] and rs202187871POT1 [Pmeta = 2.21 × 10-8]) and two synonymous variants (rs7447927TMEM173 [Pmeta = 1.32 × 10-9] and rs140624366ATRN [Pmeta = 2.97 × 10-9]). rs202197044TET3 was significantly associated with emphysema (odds ratio, 3.55; Pfdr = 0.015), whereas rs7447927POT1 was strongly associated with telomere length (β = 1.08; Pfdr (FDR corrected P value) = 3.76 × 10-53). Functional evidence of expression of quantitative trait loci, splicing quantitative trait loci, and isoform expression was found for the four novel genes. Gene-level association tests identified several novel genes, including POT1 (protection of telomeres 1), RTEL1, BSG, and ZNF232. Conclusions: Our findings provide insights into the genetic architecture of human exomes and their role in lung cancer predisposition.

Low-density lipoprotein cholesterol and risk of hepatocellular carcinoma: a Mendelian randomization and mediation analysis

BACKGROUND

A previous study demonstrated that low-density lipoprotein cholesterol (LDL-C) is associated with hepatocellular carcinoma (HCC); however, the causality between them has not been proven due to conflicting research results and the interference of confounders. This study utilized Mendelian randomization (MR) to investigate the causal relationship between LDL-C and HCC and identify the mediating factors.

METHODS

LDL-C, HCC, and coronary artery disease (CAD) genome-wide association study (GWAS) data were obtained from a public database. To investigate causality, inverse variance weighting (IVW) was the main analysis approach. MR‒Egger, simple mode, weighted median (WM), and weighted mode were employed as supplementary analytic methods. In addition, horizontal pleiotropy and heterogeneity were tested. To evaluate the stability of the MR results, a "leave-one-out" approach was used. Multivariate MR (MVMR) was utilized to correct the confounders that might affect causality, and mediation analysis was used to investigate the potential mediating effects. Finally, we used HCC risk to infer the reverse causality with LDL-C level.

RESULTS

Random effects IVW results were (LDL-C-HCC: odds ratio (OR) = 0.703, 95% confidence interval (CI) = [0.508, 0.973], P = 0.034; CAD-HCC: OR = 0.722, 95% CI = [0.645, 0.808], P = 1.50 × 10^-8; LDL-C-CAD: OR = 2.103, 95% CI = [1.862, 2.376], P = 5.65 × 10^-33), demonstrating a causal link between LDL-C levels and a lower risk of HCC. Through MVMR, after mutual correction, the causal effect of LDL-C and CAD on HCC remained significant (P < 0.05). Through mediation analysis, it was proven that CAD mediated the causative connection between LDL-C and HCC, and the proportion of mediating effect on HCC was 58.52%. Reverse MR showed that HCC could affect LDL-C levels with a negative correlation (OR_IVW = 0.979, 95% CI = [0.961, 0.997], P = 0.025).

CONCLUSION

This MR study confirmed the causal effect between LDL-C levels and HCC risk, with CAD playing a mediating role. It may provide a new view on HCC occurrence and development mechanisms, as well as new metabolic intervention targets for treatment.

Characterizing the polygenic architecture of complex traits in populations of East Asian and European descent

To investigate the polygenicity of complex traits in populations of East Asian (EAS) and European (EUR) descents, we leveraged genome-wide data from Biobank Japan, UK Biobank, and FinnGen cohorts. Specifically, we analyzed up to 215 outcomes related to 18 health domains, assessing their polygenic architecture via descriptive statistics, such as the proportion of susceptibility SNPs per trait (πc). While we did not observe EAS-EUR differences in the overall distribution of polygenicity parameters across the phenotypes investigated, there were ancestry-specific patterns in the polygenicity differences between health domains. In EAS, pairwise comparisons across health domains showed enrichment for πc differences related to hematological and metabolic traits (hematological fold-enrichment = 4.45, p = 2.15 × 10-7; metabolic fold-enrichment = 4.05, p = 4.01 × 10-6). For both categories, the proportion of susceptibility SNPs was lower than that observed for several other health domains (EAS-hematological median πc = 0.15%, EAS-metabolic median πc = 0.18%) with the strongest πc difference with respect to respiratory traits (EAS-respiratory median πc = 0.50%; hematological-p = 2.26 × 10-3; metabolic-p = 3.48 × 10-3). In EUR, pairwise comparisons showed multiple πc differences related to the endocrine category (fold-enrichment = 5.83, p = 4.76 × 10-6), where these traits showed a low proportion of susceptibility SNPs (EUR-endocrine median πc = 0.01%) with the strongest difference with respect to psychiatric phenotypes (EUR-psychiatric median πc = 0.50%; p = 1.19 × 10-4). Simulating sample sizes of 1,000,000 and 5,000,000 individuals, we also showed that ancestry-specific polygenicity patterns translate into differences across health domains in the genetic variance explained by susceptibility SNPs projected to be genome-wide significant (e.g., EAS hematological-neoplasm p = 2.18 × 10-4; EUR endocrine-gastrointestinal p = 6.80 × 10-4). These findings highlight that traits related to the same health domains may present ancestry-specific variability in their polygenicity.

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

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

Construction and evaluation of the functional polygenic risk score for gastric cancer in a prospective cohort of the European population

BACKGROUND

A polygenic risk score (PRS) derived from 112 single-nucleotide polymorphisms (SNPs) for gastric cancer has been reported in Chinese populations (PRS-112). However, its performance in other populations is unknown. A functional PRS (fPRS) using functional SNPs (fSNPs) may improve the generalizability of the PRS across populations with distinct ethnicities.

METHODS

We performed functional annotations on SNPs in strong linkage disequilibrium (LD) with the 112 previously reported SNPs to identify fSNPs that affect protein-coding or transcriptional regulation. Subsequently, we constructed an fPRS based on the fSNPs by using the LDpred2-infinitesimal model and then analyzed the performance of the PRS-112 and fPRS in the risk prediction of gastric cancer in 457,521 European participants of the UK Biobank cohort. Finally, the performance of the fPRS in combination with lifestyle factors were evaluated in predicting the risk of gastric cancer.

RESULTS

During 4,582,045 person-years of follow-up with a total of 623 incident gastric cancer cases, we found no significant association between the PRS-112 and gastric cancer risk in the European population (hazard ratio [HR] = 1.00 [95% confidence interval (CI) 0.93-1.09], P = 0.846). We identified 125 fSNPs, including seven deleterious protein-coding SNPs and 118 regulatory non-coding SNPs, and used them to constructed the fPRS-125. Our result showed that the fPRS-125 was significantly associated with gastric cancer risk (HR = 1.11 [95% CI, 1.03-1.20], P = 0.009). Compared to participants with a low fPRS-125 (bottom quintile), those with a high fPRS-125 (top quintile) had a higher risk of incident gastric cancer (HR = 1.43 [95% CI, 1.12-1.84], P = 0.005). Moreover, we observed that participants with both an unfavorable lifestyle and a high genetic risk had the highest risk of incident gastric cancer (HR = 4.99 [95% CI, 1.55-16.10], P = 0.007) compared to those with both a favorable lifestyle and a low genetic risk.

CONCLUSION

These results indicate that the fPRS-125 derived from fSNPs may act as an indicator to measure the genetic risk of gastric cancer in the European population.

Sex differences of the shared genetic landscapes between type 2 diabetes and peripheral artery disease in East Asians and Europeans

Type 2 diabetes (T2D) is a critical risk factor for peripheral artery disease (PAD). However, the sex differences in genetic basis, causality, and underlying mechanisms of the two diseases are still unclear. Using sex-stratified and ethnic-based GWAS summary, we explored the genetic correlation and causal relationship between T2D and PAD in both ethnicities and sexes by linkage disequilibrium score regression, LAVA and six Mendelian Randomization approaches. We observed stronger genetic correlations between T2D and PAD in females than males in East Asians and Europeans. East Asian females exhibit higher causal effects of T2D on PAD than males. The gene-level analysis found KCNJ11 and ANK1 genes associated with the cross-trait of T2D and PAD in both sexes. Our study provides genetic evidence for the sex difference of genetic correlations and causal relationships between PAD and T2D, indicating the importance of using sex-specific strategies for monitoring PAD in T2D patients.

Long-term effect of styrene and ethylbenzene exposure on fasting plasma glucose: A gene-environment interaction study

Styrene and ethylbenzene (S/EB) are hazardous pollutants that have attracted worldwide concern. In this prospective cohort study, S/EB exposure biomarker (the sum of mandelic acid and phenylglyoxylic acid [MA+PGA]) and fasting plasma glucose (FPG) were repeatedly measured three times. The polygenic risk score (PRS) based on 137 single nucleotide polymorphisms for type 2 diabetes mellitus (T2DM) was calculated to evaluate cumulative genetic effect. In repeated-measures cross-sectional analyses, MA+PGA (β [95% confidence interval]: 0.106 [0.022, 0.189]) and PRS (0.111 [0.047, 0.176]) were significantly related to FPG. For long-term effect assessment, participants with sustained high MA+PGA or with high PRS had 0.021 (95% CI: -0.398, 0.441) or 0.465 (0.064, 0.866) mmol/L increase in FPG, respectively, over 3 years follow-up, and had 0.256 (0.017, 0.494) or 0.265 (0.004, 0.527) mmol/L increase in FPG, respectively, over 6 years follow-up. We further detected a significant interaction effect between MA+PGA and PRS on FPG change, compared with participants with sustained low MA+PGA and low PRS, those with sustained high MA+PGA and high PRS had 0.778 (0.319, 1.258) mmol/L increase in FPG (P for interaction=0.028) over 6 years follow-up. Our study provides the first evidence that long-term exposure to S/EB potentially increases FPG, which might be aggravated by genetic susceptibility.

Detection of trait-associated structural variations using short-read sequencing

Genomic structural variation (SV) affects genetic and phenotypic characteristics in diverse organisms, but the lack of reliable methods to detect SV has hindered genetic analysis. We developed a computational algorithm (MOPline) that includes missing call recovery combined with high-confidence SV call selection and genotyping using short-read whole-genome sequencing (WGS) data. Using 3,672 high-coverage WGS datasets, MOPline stably detected ∼16,000 SVs per individual, which is over ∼1.7-3.3-fold higher than previous large-scale projects while exhibiting a comparable level of statistical quality metrics. We imputed SVs from 181,622 Japanese individuals for 42 diseases and 60 quantitative traits. A genome-wide association study with the imputed SVs revealed 41 top-ranked or nearly top-ranked genome-wide significant SVs, including 8 exonic SVs with 5 novel associations and enriched mobile element insertions. This study demonstrates that short-read WGS data can be used to identify rare and common SVs associated with a variety of traits.

MRBEE: A novel bias-corrected multivariable Mendelian Randomization method

Mendelian randomization (MR) is an instrumental variable approach used to infer causal relationships between exposures and outcomes and can apply to summary data from genome-wide association studies (GWAS). Since GWAS summary statistics are subject to estimation errors, most existing MR approaches suffer from measurement error bias, whose scale and direction are influenced by weak instrumental variables and GWAS sample overlap, respectively. We introduce MRBEE (MR using Bias-corrected Estimating Equation), a novel multivariable MR method capable of simultaneously removing measurement error bias and identifying horizontal pleiotropy. In simulations, we showed that MRBEE is capable of effectively removing measurement error bias in the presence of weak instrumental variables and sample overlap. In two independent real data analyses, we discovered that the causal effect of BMI on coronary artery disease risk is entirely mediated by blood pressure, and that existing MR methods may underestimate the causal effect of cannabis use disorder on schizophrenia risk compared to MRBEE. MRBEE possesses significant potential for advancing genetic research by providing a valuable tool to study causality between multiple risk factors and disease outcomes, particularly as a large number of GWAS summary statistics become publicly available.

Polygenic risk scores for prediction of breast cancer in Korean women

BACKGROUND

Polygenic risk scores (PRSs) for breast cancer, developed using European and Asian genome-wide association studies (GWAS), have been shown to have good discrimination in Asian women. However, prospective calibration of absolute risk prediction models, based on a PRS or PRS combined with lifestyle, clinical and environmental factors, in Asian women is limited.

METHODS

We consider several PRSs trained using European and/or Asian GWAS. For each PRS, we evaluate the discrimination and calibration of three absolute risk models among 41 031 women from the Korean Cancer Prevention Study (KCPS)-II Biobank: (i) a model using incidence, mortality and risk factor distributions (reference inputs) among US women and European relative risks; (ii) a recalibrated model, using Korean reference but European relative risks; and (iii) a fully Korean-based model using Korean reference and relative risk estimates from KCPS.

RESULTS

All Asian and European PRS improved discrimination over lifestyle, clinical and environmental (Qx) factors in Korean women. US-based absolute risk models overestimated the risks for women aged ≥50 years, and this overestimation was larger for models that only included PRS (expected-to-observed ratio E/O = 1.2 for women <50, E/O = 2.7 for women ≥50). Recalibrated and Korean-based risk models had better calibration in the large, although the risk in the highest decile was consistently overestimated. Absolute risk projections suggest that risk-reducing lifestyle changes would lead to larger absolute risk reductions among women at higher PRS.

CONCLUSIONS

Absolute risk models incorporating PRS trained in European and Asian GWAS and population-appropriate average age-specific incidences may be useful for risk-stratified interventions in Korean women.

A Multi-Ancestry Polygenic Risk Score for Coronary Heart Disease Based on an Ancestrally Diverse Genome-Wide Association Study and Population-Specific Optimization

Background

Predictive performance of polygenic risk scores (PRS) varies across populations. To facilitate equitable clinical use, we developed PRS for coronary heart disease (PRSCHD) for 5 genetic ancestry groups.

Methods

We derived ancestry-specific and multi-ancestry PRSCHD based on pruning and thresholding (PRSP+T) and continuous shrinkage priors (PRSCSx) applied on summary statistics from the largest multi-ancestry genome-wide meta-analysis for CHD to date, including 1.1 million participants from 5 continental populations. Following training and optimization of PRSCHD in the Million Veteran Program, we evaluated predictive performance of the best performing PRSCHD in 176,988 individuals across 9 cohorts of diverse genetic ancestry.

Results

Multi-ancestry PRSP+T outperformed ancestry specific PRSP+T across a range of tuning values. In training stage, for all ancestry groups, PRSCSx performed better than PRSP+T and multi-ancestry PRS outperformed ancestry-specific PRS. In independent validation cohorts, the selected multi-ancestry PRSP+T demonstrated the strongest association with CHD in individuals of South Asian (SAS) and European (EUR) ancestry (OR per 1SD[95% CI]; 2.75[2.41-3.14], 1.65[1.59-1.72]), followed by East Asian (EAS) (1.56[1.50-1.61]), Hispanic/Latino (HIS) (1.38[1.24-1.54]), and weakest in African (AFR) ancestry (1.16[1.11-1.21]). The selected multi-ancestry PRSCSx showed stronger associacion with CHD in comparison within each ancestry group where the association was strongest in SAS (2.67[2.38-3.00]) and EUR (1.65[1.59-1.71]), progressively decreasing in EAS (1.59[1.54-1.64]), HIS (1.51[1.35-1.69]), and lowest in AFR (1.20[1.15-1.26]).

Conclusions

Utilizing diverse summary statistics from a large multi-ancestry genome-wide meta-analysis led to improved performance of PRSCHD in most ancestry groups compared to single-ancestry methods. Improvement of predictive performance was limited, specifically in AFR and HIS, despite use of one of the largest and most diverse set of training and validation cohorts to date. This highlights the need for larger GWAS datasets of AFR and HIS individuals to enhance performance of PRSCHD.

Evaluating the role of non-alcoholic fatty liver disease in cardiovascular diseases and type 2 diabetes: a Mendelian randomization study in Europeans and East Asians

BACKGROUND

Whether non-alcoholic fatty liver disease (NAFLD) causes cardiovascular disease (CVD) and type 2 diabetes (T2D) is unclear and possible differences between ethnicities have not been thoroughly explored. We used Mendelian randomization (MR) to assess the role of NAFLD in CVD and T2D risk in Europeans and East Asians.

METHODS

We conducted a MR study using genetic predictors of alanine aminotransferase (ALT), liability to NAFLD, aspartate transaminase (AST), liver magnetic resonance imaging corrected T1 and proton density fat fraction and combined them with genome-wide association studies (GWAS) summary statistics of CVD, T2D and glycaemic traits (sample size ranging from 14 400 to 977 320). Inverse-variance weighted analysis was used to assess the effect of NAFLD in these outcomes, with sensitivity analyses and replication in FinnGen. We conducted analyses in East Asians using ethnicity-specific genetic predictors of ALT and AST, and the respective outcome GWAS summary statistics.

RESULTS

In Europeans, higher ALT was associated with higher T2D risk (odds ratio: 1.77 per standard deviation, 95% CI 1.5 to 2.08), with similar results for other exposures, across sensitivity analyses and in FinnGen. Although NAFLD proxies were related to higher coronary artery disease (CAD) and stroke risk, sensitivity analyses suggested possible bias by horizontal pleiotropy. In East Asians, higher ALT was possibly associated with higher T2D risk, and ALT and AST were inversely associated with CAD.

CONCLUSIONS

NAFLD likely increases the risk of T2D in Europeans and East Asians. Potential differential effects on CAD between Europeans and East Asians require further investigation.

Characterizing the polygenic architecture of complex traits in populations of East Asian and European descent

To investigate the polygenicity of complex traits in populations of East Asian (EAS) and European (EUR) descents, we leveraged genome-wide data from Biobank Japan, UK Biobank, and FinnGen cohorts. Specifically, we analyzed up to 215 outcomes related to 18 health domains, assessing their polygenic architecture via descriptive statistics, such as the proportion of susceptibility SNPs per trait (π _c ). While we did not observe EAS-EUR differences in the overall distribution of polygenicity parameters across the phenotypes investigated, there were ancestry-specific patterns in the polygenicity differences between health domains. In EAS, pairwise comparisons across health domains showed enrichment for π _c differences related to hematological and metabolic traits (hematological fold-enrichment=4.45, p=2.15×10 ^-7 ; metabolic fold-enrichment=4.05, p=4.01×10 ^-6 ). For both categories, the proportion of susceptibility SNPs was lower than that observed for several other health domains (EAS-hematological median π _c =0.15%, EAS-metabolic median π _c =0.18%) with the strongest π _c difference with respect to respiratory traits (EAS-respiratory median π _c =0.50%; Hematological-p=2.26×10 ^-3 ; Metabolic-p=3.48×10 ^-3 ). In EUR, pairwise comparisons showed multiple π _c differences related to the endocrine category (fold-enrichment=5.83, p=4.76×10 ^-6 ), where these traits showed a low proportion of susceptibility SNPs (EUR-endocrine median π _c =0.01%) with the strongest difference with respect to psychiatric phenotypes (EUR-psychiatric median π _c =0.50%; p=1.19×10 ^-4 ). Simulating sample sizes of 1,000,000 and 5,000,000 individuals, we also showed that ancestry-specific polygenicity patterns translate into differences across health domains in the genetic variance explained by susceptibility SNPs projected to be genome-wide significant (e.g., EAS hematological-neoplasm p=2.18×10 ^-4 ; EUR endocrine-gastrointestinal p=6.80×10 ^-4 ). These findings highlight that traits related to the same health domains may present ancestry-specific variability in their polygenicity.

Identification of a Common Variant for Coronary Heart Disease at PDE1A Contributes to Individualized Treatment Goals and Risk Stratification of Cardiovascular Complications in Chinese Patients With Type 2 Diabetes

OBJECTIVE

In this study we aim to unravel genetic determinants of coronary heart disease (CHD) in type 2 diabetes (T2D) and explore their applications.

RESEARCH DESIGN AND METHODS

We performed a two-stage genome-wide association study for CHD in Chinese patients with T2D (3,596 case and 8,898 control subjects), followed by replications in European patients with T2D (764 case and 4,276 control subjects) and general populations (n = 51,442-547,261). Each identified variant was examined for its association with a wide range of phenotypes and its interactions with glycemic, blood pressure (BP), and lipid controls in incident cardiovascular diseases.

RESULTS

We identified a novel variant (rs10171703) for CHD (odds ratio 1.21 [95% CI 1.13-1.30]; P = 2.4 × 10-8) and BP (β ± SE 0.130 ± 0.017; P = 4.1 × 10-14) at PDE1A in Chinese T2D patients but found only a modest association with CHD in general populations. This variant modulated the effects of BP goal attainment (130/80 mmHg) on CHD (Pinteraction = 0.0155) and myocardial infarction (MI) (Pinteraction = 5.1 × 10-4). Patients with CC genotype of rs10171703 had >40% reduction in either cardiovascular events in response to BP control (2.9 × 10-8 < P < 3.6 × 10-5), those with CT genotype had no difference (0.0726 < P < 0.2614), and those with TT genotype had a threefold increase in MI risk (P = 6.7 × 10-3).

CONCLUSIONS

We discovered a novel CHD- and BP-related variant at PDE1A that interacted with BP goal attainment with divergent effects on CHD risk in Chinese patients with T2D. Incorporating this information may facilitate individualized treatment strategies for precision care in diabetes, only when our findings are validated.

Causal associations between chronic hepatitis B and COVID-19 in East Asian populations

BACKGROUND

The relationship between chronic hepatitis B (CHB) and Coronavirus disease 2019 (COVID-19) has been inconsistent in traditional observational studies.

METHODS

We explored the total causal and direct causal associations between CHB and the three COVID-19 outcomes using univariate and multivariate Mendelian randomization (MR) analyses, respectively. Genome-wide association study datasets for CHB and COVID-19 were obtained from the Japan Biobank and the COVID-19 Host Genetics Initiative, respectively.

RESULTS

Univariate MR analysis showed that CHB increased the risk of SARS-CoV-2 infection (OR = 1.04, 95% CI 1.01-1.07, P = 3.39E-03), hospitalized COVID-19 (OR = 1.10, 95% CI 1.06-1.13, P = 7.31E-08), and severe COVID-19 (OR = 1.16, 95%CI 1.08-1.26, P = 1.43E-04). A series of subsequent sensitivity analyses ensured the stability and reliability of these results. In multivariable MR analyses adjusting for type 2 diabetes, body mass index, basophil count, and smoking, genetically related CHB is still positively associated with increased risk of SARS-CoV-2 infection (OR = 1.06, 95% CI 1.02-1.11, P = 1.44E-03) and hospitalized COVID-19 (OR = 1.12, 95% CI 1.07-1.16, P = 5.13E-07). However, the causal link between CHB and severe COVID-19 was attenuated after adjustment for the above variables. In addition, the MR analysis did not support the causal effect of COVID-19 on CHB.

CONCLUSIONS

This study provides evidence that CHB increases COVID-19 susceptibility and severity among individuals of East Asian ancestry.

Gene-environment interaction in long-term effects of polychlorinated biphenyls exposure on glucose homeostasis and type 2 diabetes: The modifying effects of genetic risk and lifestyle

The longitudinal relationships of polychlorinated biphenyls (PCBs) exposure with glucose homeostasis and type 2 diabetes (T2D) risk among Chinese population have not been assessed, and interactions of PCB exposure with genetic susceptibility and lifestyle are unclear. In this prospective cohort study, fasting plasma glucose (FPG) and insulin (FPI) and seven serum indicator-PCBs were measured for each participant. We constructed polygenic risk score (PRS) of T2D and healthy lifestyle score. Each 1-unit increment of ln-transformed PCB-118 was related with a 0.141 mmol/L, 11.410 pmol/L, 0.661, and 74.5% increase in FPG, FPI, homeostasis model assessment of insulin resistance, and incident T2D risk over 6 years, respectively. Each 1-unit increment in T2D-PRS was related with a 0.169 mmol/L elevation of FPG and 65.5% elevation of incident T2D risk during 6 years. Compared with participants who had low T2D-PRS and low PCB-118, participants with high T2D-PRS and high PCB-118 showed a significant increase in FPG (0.162 mmol/L; P for interaction <0.001) and incident T2D risk [hazard ratio (HR)= 2.222]. Participants with low PCB-118, low PRS, and healthy lifestyle had the lowest incident T2D risk (HR=0.232). Our findings highlighted the significance of reducing PCB exposure and improvement in lifestyle for T2D prevention and management, especially for individuals with higher genetic risk of T2D.

The relationship between the serum lipid profile and hepatocellular carcinoma in east Asian population: A mendelian randomization study

Background

Although several studies have found that the serum lipid profile may be correlated with hepatocellular carcinoma (HCC), the causal relationships between the serum lipid profile and HCC have not been determined due to potential confounder. Here, Mendelian randomization (MR) analysis was performed to identify the relationship between the serum lipid profile and HCC in the East Asian population.

Method

Our study made a MR analysis with the validation of two data sets. We obtained genome-wide association study (GWAS) data related to the serum lipid profile from Asian Genetic Epidemiology Network (AGEN). Then, the data from a recent large GWAS of the East Asian ancestry in Japan (BioBank Japan, BBJ) were extracted. Summary-level statistical data for HCC were obtained from a large GWAS of the East Asian ancestry in Japan. Univariable MR analysis were performed to identify whether the genetic evidence of serum lipid profile was significantly associated with HCC risk. Multivariable MR analysis was conducted to estimate the independent effects of exposures on HCC.

Results

Univariable and multivariable MR analyses indicated that the serum lipid profile was not a risk factor for HCC incidence in either data set based on the East Asian population. Multivariable MR analysis revealed that the hazard ratios of the probability of HCC in AGEN were 1.134 (95% confidence interval (CI), 0.903-1.424) for TG, 1.010 (95% CI: 0.824-1.237) for HDL-C, 0.974 (95% CI: 0.746-1.271) for TC, 0.918 (95% CI: 0.734-1.147) for LDL-C, while the results in BBJ were also non-significant: 1.111 (95% CI: 0.869-1.419) for TG, 0.957 (95% CI: 0.790-1.158) for HDL-C, 0.917 (95% CI: 0.643-1.308) for TC, 0.932 (95% CI: 0.699-1.243) for LDL-C.

Conclusion

Our MR study with the validation of two data sets found no strong evidence to support causal associations between the serum lipid profile and HCC risk.

Aberrant inflammasome activation as a driving force of human autoimmune skin disease

Autoimmune skin diseases are understood as conditions in which the adaptive immune system with autoantigen-specific T cells and autoantibody-producing B cells reacting against self-tissues plays a crucial pathogenic role. However, there is increasing evidence that inflammasomes, which are large multiprotein complexes that were first described 20 years ago, contribute to autoimmune disease progression. The inflammasome and its contribution to the bioactivation of interleukins IL-1β and IL-18 play an essential role in combating foreign pathogens or tissue damage, but may also act as a pathogenic driver of myriad chronic inflammatory diseases when dysfunctionally regulated. Inflammasomes containing the NOD-like receptor family members NLRP1 and NLRP3 as well as the AIM2-like receptor family member AIM2 have been increasingly investigated in inflammatory skin conditions. In addition to autoinflammatory diseases, which are often associated with skin involvement, the aberrant activation of the inflammasome has also been implied in autoimmune diseases that can either affect the skin besides other organs such as systemic lupus erythematosus and systemic sclerosis or are isolated to the skin in humans. The latter include, among others, the T-cell mediated disorders vitiligo, alopecia areata, lichen planus and cutaneous lupus erythematosus as well as the autoantibody-driven blistering skin disease bullous pemphigoid. Some diseases are characterized by both autoinflammatory and autoimmune responses such as the chronic inflammatory skin disease psoriasis. Further insights into inflammasome dysregulation and associated pathways as well as their role in forming adaptive immune responses in human autoimmune skin pathology could potentially offer a new field of therapeutic options in the future.

Dissecting the genetic heterogeneity of gastric cancer

BACKGROUND

Gastric cancer (GC) is clinically heterogenous according to location (cardia/non-cardia) and histopathology (diffuse/intestinal). We aimed to characterize the genetic risk architecture of GC according to its subtypes. Another aim was to examine whether cardia GC and oesophageal adenocarcinoma (OAC) and its precursor lesion Barrett's oesophagus (BO), which are all located at the gastro-oesophageal junction (GOJ), share polygenic risk architecture.

METHODS

We did a meta-analysis of ten European genome-wide association studies (GWAS) of GC and its subtypes. All patients had a histopathologically confirmed diagnosis of gastric adenocarcinoma. For the identification of risk genes among GWAS loci we did a transcriptome-wide association study (TWAS) and expression quantitative trait locus (eQTL) study from gastric corpus and antrum mucosa. To test whether cardia GC and OAC/BO share genetic aetiology we also used a European GWAS sample with OAC/BO.

FINDINGS

Our GWAS consisting of 5816 patients and 10,999 controls highlights the genetic heterogeneity of GC according to its subtypes. We newly identified two and replicated five GC risk loci, all of them with subtype-specific association. The gastric transcriptome data consisting of 361 corpus and 342 antrum mucosa samples revealed that an upregulated expression of MUC1, ANKRD50, PTGER4, and PSCA are plausible GC-pathomechanisms at four GWAS loci. At another risk locus, we found that the blood-group 0 exerts protective effects for non-cardia and diffuse GC, while blood-group A increases risk for both GC subtypes. Furthermore, our GWAS on cardia GC and OAC/BO (10,279 patients, 16,527 controls) showed that both cancer entities share genetic aetiology at the polygenic level and identified two new risk loci on the single-marker level.

INTERPRETATION

Our findings show that the pathophysiology of GC is genetically heterogenous according to location and histopathology. Moreover, our findings point to common molecular mechanisms underlying cardia GC and OAC/BO.

FUNDING

German Research Foundation (DFG).

Mobile element variation contributes to population-specific genome diversification, gene regulation and disease risk

Mobile genetic elements (MEs) are heritable mutagens that recursively generate structural variants (SVs). ME variants (MEVs) are difficult to genotype and integrate in statistical genetics, obscuring their impact on genome diversification and traits. We developed a tool that accurately genotypes MEVs using short-read whole-genome sequencing (WGS) and applied it to global human populations. We find unexpected population-specific MEV differences, including an Alu insertion distribution distinguishing Japanese from other populations. Integrating MEVs with expression quantitative trait loci (eQTL) maps shows that MEV classes regulate tissue-specific gene expression by shared mechanisms, including creating or attenuating enhancers and recruiting post-transcriptional regulators, supporting class-wide interpretability. MEVs more often associate with gene expression changes than SNVs, thus plausibly impacting traits. Performing genome-wide association study (GWAS) with MEVs pinpoints potential causes of disease risk, including a LINE-1 insertion associated with keloid and fasciitis. This work implicates MEVs as drivers of human divergence and disease risk.

Polygenic Risk Score, Lifestyles, and Type 2 Diabetes Risk: A Prospective Chinese Cohort Study

The aim of this study was to generate a polygenic risk score (PRS) for type 2 diabetes (T2D) and test whether it could be used in identifying high-risk individuals for lifestyle intervention in a Chinese cohort. We genotyped 80 genetic variants among 5024 participants without non-communicable diseases at baseline in the Wuxi Non-Communicable Diseases cohort (Wuxi NCDs cohort). During the follow-up period of 14 years, 440 cases of T2D were newly diagnosed. Using Cox regression, we found that the PRS of 46 SNPs identified by the East Asians was relevant to the future T2D. Participants with a high PRS (top quintile) had a two-fold higher risk of T2D than the bottom quintile (hazard ratio: 2.06, 95% confidence interval: 1.42-2.97). Lifestyle factors were considered, including cigarette smoking, alcohol consumption, physical exercise, diet, body mass index (BMI), and waist circumference (WC). Among high-PRS individuals, the 10-year incidence of T2D slumped from 6.77% to 3.28% for participants having ideal lifestyles (4-6 healthy lifestyle factors) compared with poor lifestyles (0-2 healthy lifestyle factors). When integrating the high PRS, the 10-year T2D risk of low-clinical-risk individuals exceeded that of high-clinical-risk individuals with a low PRS (3.34% vs. 2.91%). These findings suggest that the PRS of 46 SNPs could be used in identifying high-risk individuals and improve the risk stratification defined by traditional clinical risk factors for T2D. Healthy lifestyles can reduce the risk of a high PRS, which indicates the potential utility in early screening and precise prevention.

FHL5 Controls Vascular Disease-Associated Gene Programs in Smooth Muscle Cells

BACKGROUND

Genome-wide association studies have identified hundreds of loci associated with common vascular diseases, such as coronary artery disease, myocardial infarction, and hypertension. However, the lack of mechanistic insights for many GWAS loci limits their translation into the clinic. Among these loci with unknown functions is UFL1-four-and-a-half LIM (LIN-11, Isl-1, MEC-3) domain 5 (FHL5; chr6q16.1), which reached genome-wide significance in a recent coronary artery disease/ myocardial infarction GWAS meta-analysis. UFL1-FHL5 is also associated with several vascular diseases, consistent with the widespread pleiotropy observed for GWAS loci.

METHODS

We apply a multimodal approach leveraging statistical fine-mapping, epigenomic profiling, and ex vivo analysis of human coronary artery tissues to implicate FHL5 as the top candidate causal gene. We unravel the molecular mechanisms of the cross-phenotype genetic associations through in vitro functional analyses and epigenomic profiling experiments in coronary artery smooth muscle cells.

RESULTS

We prioritized FHL5 as the top candidate causal gene at the UFL1-FHL5 locus through expression quantitative trait locus colocalization methods. FHL5 gene expression was enriched in the smooth muscle cells and pericyte population in human artery tissues with coexpression network analyses supporting a functional role in regulating smooth muscle cell contraction. Unexpectedly, under procalcifying conditions, FHL5 overexpression promoted vascular calcification and dysregulated processes related to extracellular matrix organization and calcium handling. Lastly, by mapping FHL5 binding sites and inferring FHL5 target gene function using artery tissue gene regulatory network analyses, we highlight regulatory interactions between FHL5 and downstream coronary artery disease/myocardial infarction loci, such as FOXL1 and FN1 that have roles in vascular remodeling.

CONCLUSIONS

Taken together, these studies provide mechanistic insights into the pleiotropic genetic associations of UFL1-FHL5. We show that FHL5 mediates vascular disease risk through transcriptional regulation of downstream vascular remodeling gene programs. These transacting mechanisms may explain a portion of the heritable risk for complex vascular diseases.

Diseases of the musculoskeletal system and connective tissue and risk of breast cancer: Mendelian randomization study in European and East Asian populations

Objective

Associations between diseases of the musculoskeletal system and connective tissue (MSCTD) and breast cancer (BC) have not been elucidated completely. The purpose of this study was to investigate the associations of MSCTD, rheumatoid arthritis (RA), Sjogren syndrome (SS), systemic lupus erythematosus (SLE), systemic sclerosis (SSc), dermatomyositis (DM), polymyositis (PM), osteoarthritis (OA) of hip or knee, and ankylosing spondylitis (AS) with BC in European populations and East Asian populations using Mendelian randomized (MR) analysis.

Methods

The genetic instruments linked to MSCTD, RA, SS, SLE, SSc, DM, PM, OA, and AS were chosen from the EBI database of complete genome-wide association studies (GWAS) summary data and the FinnGen consortium. The associations of genetic variants with BC were extracted from the Breast Cancer Association Consortium (BCAC). Two Sample MR was performed using summary data from GWAS, principally using the inverse variant weighted (IVW) method. Heterogeneity, pleiotropy, and sensitivity analyses were performed to evaluate the robustness of the results by weighted median, MR Egger, simple mode, weighted mode, and leave-one-out analysis.

Results

In the European population, causal relationships between RA and BC (OR=1.04, 95%CI: 1.01-1.07, P=0.023), AS and BC (OR=1.21, 95%CI: 1.06-1.36, P=0.013) were confirmed. IVW analysis showed DM (OR=0.98, 95%CI: 0.96-0.99, P=0.026) and PM (OR=0.98, 95%CI: 0.97-0.99, P=0.002) were associated with slightly decreased risks of estrogen receptor (ER)+ BC, and MSCTD was associated with an increased risk of ER- BC (OR=1.85, 95%CI: 1.27-2.44, P=0.039). There was no causal relationship between SLE, SS, SSc, OA, and BC, neither ER+ BC nor ER- BC. However, in the East Asian population, IVW analysis showed that RA (OR=0.94, 95%CI: 0.89-0.99, P=0.0096) and SLE (OR=0.95, 95%CI: 0.92-0.99, P=0.0058) was associated with decreased risks of BC.

Conclusions

This study suggests that causal relationships between patients with MSCTD and BC in the European population are different from those in the East Asian population, patients with RA and AS in the European population have an increased risk of BC, patients with MSCTD have increased risk of ER- BC in the European population, while patients with RA and SLE in the East Asian population have decreased risk of BC.

Leveraging IgG N-glycosylation to infer the causality between T2D and hypertension

BACKGROUND

Observational studies demonstrated a bidirectional association between type 2 diabetes (T2D) and hypertension, whereas Mendelian randomization (MR) analyses supported the causality from T2D to hypertension but not causal from hypertension to T2D. We previously found that IgG N-glycosylation is associated with both T2D and hypertension, and thus IgG N-glycosylation might link the causality between them.

METHODS

We carried out a genome-wide association study (GWAS) to identify IgG N-glycosylation-quantitative-trait loci (QTLs) integrating GWAS for T2D and hypertension and then performed bidirectional univariable and multivariable MR analyses to infer the causal association among them. The inverse-variance-weighted (IVW) analysis was performed as the primary analysis, followed by some sensitivity analyses to explore the stability of the results.

RESULTS

Six putatively causal IgG N-glycans for T2D and four for hypertension were identified in the IVW method. Genetically predicted T2D increased the risk of hypertension (odds ratio [OR] = 1.177, 95% confidence interval (95% CI) = 1.037-1.338, P = 0.012) and vice versa (OR = 1.391, 95% CI = 1.081-1.790, P = 0.010). Multivariable MR showed that T2D remained at risk effect with hypertension ([OR] = 1.229, 95% CI = 1.140-1.325, P = 7.817 × 10^-8) after conditioning on T2D-related IgG-glycans. Conversely, hypertension was associated with higher T2D risk (OR = 1.287, 95% CI = 1.107-1.497, P = 0.001) after adjusting for related IgG-glycans. No evidence of horizontal pleiotropy was observed, as MR‒Egger regression provided P values for intercept > 0.05.

CONCLUSION

Our study validated the mutual causality between T2D and hypertension from the perspective of IgG N-glycosylation, further validating the "common soil" hypothesis underlying the pathogenesis of T2D and hypertension.

An attempt to establish real-world databases of poly(ADP-ribose) polymerase inhibitors for advanced or recurrent epithelial ovarian cancer: the Japanese Gynecologic Oncology Group

The development of new treatments for gynecological malignancies has been conducted mainly through collaborative international phase III trials led by the United States and Europe. The survival outcomes of many gynecological malignancies have greatly improved as a result. Recent large-scale genome-wide association studies have revealed that drug efficacy and adverse event profiles are not always uniform. Thus, it is important to validate new treatment options in each country to safely and efficiently provide newly developed treatment options to patients with gynecological malignancies. The Japanese Gynecologic Oncology Group (JGOG) is conducting 5 cohort studies (JGOG 3026, 3027, 3028, 3030, and 3031) to establish real-world data (RWD) of poly(ADP-ribose) polymerase (PARP) inhibitor use in patients with advanced or recurrent epithelial ovarian cancer. The RWD constructed will be used to provide newly developed PARP inhibitors for women with advanced or recurrent ovarian cancer in a safer and more efficient manner as well as to develop further treatment options. In 2022, The JGOG, Korean Gynecologic Oncology Group, Chinese Gynecologic Cancer Society, and Taiwanese Gynecologic Oncology Group established the East Asian Gynecologic Oncology Trial Group to collaborate with East Asian countries in clinical research on gynecologic malignancies and disseminate new knowledge on gynecologic malignancies from Asia. The JGOG will conduct a collaborative integrated analysis of the RWD generated from Asian countries and disseminate real-world clinical knowledge regarding new treatment options that have been clinically implemented.

Genetic Risk Score for Prediction of Coronary Heart Disease in the Korean Genome and Epidemiology Study

Background

Using a genetic risk score (GRS) to predict coronary heart disease (CHD) may detect disease earlier. The current study aims to assess whether GRS is associated with CHD incidence and whether it is clinically useful for improving prediction using traditional risk factors (TRFs) as well as family history.

Methods

Data from a total of 48,941 participants in the Korean Genome and Epidemiology Study were analyzed in the current study. The weighted GRS was constructed using 55 single-nucleotide polymorphisms based on published genome-wide association studies. The association of GRS with incident CHD was analyzed using Cox proportional hazard model. Discrimination and reclassification were assessed to demonstrate the clinical utility of GRS. The analyses were performed separately by sex.

Results

After adjusting for family history and TRFs, GRS was significantly associated with CHD incidence in men; compared to the low GRS group, men in the high GRS group had a 2.07-fold increased risk of CHD (95% confidence interval [CI]: 1.51-2.85). In men, the combination of TRFs, family history, and GRS had better performance than TRFs alone (C statistics for TRF-only model, 0.66, 95% CI, 0.64-0.69; C statistics for combination model, 0.68, 95% CI, 0.65-0.71; category-free reclassification index, 15%). In women, however, there was no significant association between GRS and CHD and no improvement between models.

Conclusions

GRS was associated with CHD incidence and contributed to a small improvement of CHD prediction in men. The potential clinical use of GRS may not outweigh the value of family history.

Genome-wide association study of tuberculosis in the western Chinese Han and Tibetan population

Tuberculosis (TB) remains a serious global public health threat. Accumulated evidence has demonstrated that human susceptibility to TB has a strong genetic basis. And different susceptibility single nucleotide polymorphisms (SNP) have been reported in different studies. To gain greater insight into the host susceptibility to TB, we perform a two-stage genome-wide association study to identify the susceptible loci of TB. In the discovery stage, 3116 (1532 TB patients and 1584 healthy controls) and 439 (211 TB patients and 228 healthy controls) individuals were genome-wide genotyped from a western Chinese Han and Tibetan population, respectively. Based on the additive genetic model, we discovered 14 and three independent loci that had potential associations with TB susceptibility in the Chinese Han and Tibetan populations, respectively (p < 1 × 10-5). Furthermore, we conducted an imputation-based meta-analysis on another two East Asia cohorts to replicate our findings. We identified one independent locus harbored by the human leukocyte antigen (HLA) class II genes that was genome-wide significantly associated with TB (lead SNP rs111875628 with a p-value of 2.20 × 10-9). Our findings suggest a novel mechanism of the interaction with the HLA class II genes and reinforce the importance of the HLA class II alleles in response to TB.

New insights from genetic studies of eczema

Genome-wide association studies (GWAS) provided fundamental insight into the genetic determinants of complex allergic diseases. For eczema, 58 susceptibility loci were reported. Protein-changing variants were associated with eczema at genome-wide significance at 12 loci. The majority of risk variants were, however, located in non-coding, regulatory regions of the genome. Prioritized target genes were enriched in pathways of the immune response and of epithelial barrier function. Interestingly, a large overlap in the genetic architecture underlying different allergic diseases was identified pointing to common pathomechanisms for eczema, asthma, hay fever, and food allergy. Here, we review the most recent findings from GWAS for eczema including the role of rare variants and genetic heterogeneity in ethnically diverse populations. In addition, we provide an overview of genes underlying Mendelian disorders featuring eczematous skin inflammation.

Association between serum uric acid and prostate cancer risk in East Asian populations: a Mendelian randomization study

PURPOSE

Previous observational studies showed that serum uric acid (SUA) was associated with prostate cancer, but the causal relationship is unclear. This study aimed to explore the potential causal association between SUA and prostate cancer risk using Mendelian randomization (MR) analyses in the East Asian populations.

METHODS

Publicly available summary-level genome-wide association studies (GWAS) data on SUA were obtained from a genome-wide meta-analysis of three Japanese cohorts (121,745 subjects). The GWAS data on prostate cancer were derived from Biobank Japan (109,347 subjects with 5,408 cases and 103,939 controls). A total of 34 SUA-related single-nucleotide polymorphisms (SNPs) (P value < 5 × 10^-8) were identified as instrumental variables. The inverse variance weighted method was used as the primary method to compute the odds ratios (ORs) and 95% confidence intervals (95% CIs) for per standard deviation increase in SUA. MR Egger, weighted median, and weighted mode were also applied to test the robustness of the results.

RESULTS

Genetically predicted SUA was positively associated with prostate cancer risk using inverse variance weighted (OR = 1.12; 95% CI 1.00-1.26; P = 0.043). The positive association was robust when MR Egger (OR = 1.16; 95% CI 1.01-1.34; P = 0.048), weighted median (OR = 1.18; 95% CI 1.03-1.36; P = 0.018), and weighted mode (OR = 1.14; 95% CI 1.01-1.29; P = 0.041) were used.

CONCLUSION

There were potential causal associations between higher genetically predicted SUA levels and increased prostate cancer risk. Further, MR studies with more valid SNPs and more cancer cases are needed. Validation of the findings is also recommended.

Mendelian randomization analyses reveal causal relationships between the human microbiome and longevity

Although recent studies have revealed the association between the human microbiome especially gut microbiota and longevity, their causality remains unclear. Here, we assess the causal relationships between the human microbiome (gut and oral microbiota) and longevity, by leveraging bidirectional two-sample Mendelian randomization (MR) analyses based on genome-wide association studies (GWAS) summary statistics of the gut and oral microbiome from the 4D-SZ cohort and longevity from the CLHLS cohort. We found that some disease-protected gut microbiota such as Coriobacteriaceae and Oxalobacter as well as the probiotic Lactobacillus amylovorus were related to increased odds of longevity, whereas the other gut microbiota such as colorectal cancer pathogen Fusobacterium nucleatum, Coprococcus, Streptococcus, Lactobacillus, and Neisseria were negatively associated with longevity. The reverse MR analysis further revealed genetically longevous individuals tended to have higher abundances of Prevotella and Paraprevotella but lower abundances of Bacteroides and Fusobacterium species. Few overlaps of gut microbiota-longevity interactions were identified across different populations. We also identified abundant links between the oral microbiome and longevity. The additional analysis suggested that centenarians genetically had a lower gut microbial diversity, but no difference in oral microbiota. Our findings strongly implicate these bacteria to play a role in human longevity and underscore the relocation of commensal microbes among different body sites that would need to be monitored for long and healthy life.

Associations of genetic markers of diabetes mellitus with carotid atherosclerosis: a community-based case-control study

BACKGROUND

Diabetes mellitus (DM) is a well-established determinant of atherosclerosis and cardiovascular diseases (CVD). Recently, genome-wide association studies (GWAS) identified several single nucleotide polymorphism (SNP) significantly correlated with DM. The study aimed to explore the relationships of the top significant DM SNPs with carotid atherosclerosis (CA).

METHODS

We used a case-control design and randomly selected 309 cases and 439 controls with and without, respectively, carotid plaque (CP) from a community-based cohort. Eight recent GWAS on DM in East Asians reported hundreds of SNPs with genome-wide significance. The study used the top significant DM SNPs, with a p-value < 10^-16, as the candidate genetic markers of CA. The independent effects of these DM SNPs on CA were assessed by multivariable logistic regression analyses to control the effects of conventional cardio-metabolic risk factors.

RESULTS

Multivariable analyses showed that, 9 SNPs, including rs4712524, rs1150777, rs10842993, rs2858980, rs9583907, rs1077476, rs7180016, rs4383154, and rs9937354, showed promising associations with the presence of carotid plaque (CP). Among them, rs9937354, rs10842993, rs7180016, and rs4383154 showed significantly independent effects. The means (SD) of the 9-locus genetic risk score (9-GRS) of CP-positive and -negative subjects were 9.19 (1.53) and 8.62 (1.63), respectively (p < 0.001). The corresponding values of 4-locus GRS (4-GRS) were 4.02 (0.81) and. 3.78 (0.92), respectively (p < 0.001). The multivariable-adjusted odds ratio of having CP for per 1.0 increase in 9-GRS and 4-GRS were 1.30 (95% CI 1.18-1.44; p = 4.7 × 10^-7) and 1.47 (95% CI 1.74-9.40; p = 6.1 × 10^-5), respectively. The means of multi-locus GRSs of DM patients were similar to those of CP-positive subjects and higher than those of CP-negative or DM-negative subjects.

CONCLUSIONS

We identified 9 DM SNPs showing promising associations with CP. The multi-locus GRSs may be used as biomarkers for the identification and prediction of high-risks subjects for atherosclerosis and atherosclerotic diseases. Future studies on these specific SNPs and their associated genes may provide valuable information for the preventions of DM and atherosclerosis.

Associations of alcohol and coffee with colorectal cancer risk in East Asian populations: a Mendelian randomization study

PURPOSE

Previous observational studies have shown that alcohol and coffee were associated with colorectal cancer (CRC) risk, but the causal relationships have not been adequately explored. This study aimed to assess the potential causal associations of alcohol and coffee with CRC risk using Mendelian randomization (MR) analyses in an East Asian population.

METHODS

Publicly available summary-level genome-wide association studies data on ever/never alcohol drinker (n = 165,084), alcohol consumption (n = 58,610), coffee consumption (n = 152,634), and CRC (7062 cases and 195,745 controls) were obtained from the BioBank Japan (BBJ). Single-nucleotide polymorphisms (SNPs) that were significantly related to the exposures were identified as instrumental variables. Five, two, and six SNPs were used for ever/never alcohol drinkers, alcohol consumption, and coffee consumption, respectively. The inverse variance weighted method was used as the main MR method to calculate the odds ratios (ORs) and 95% confidence intervals (95% CIs) of CRC risk per one-unit change in exposures.

RESULTS

Genetically predicted ever/never alcohol drinkers (OR: 1.08; 95% CI 1.06, 1.11; P < 0.001) and alcohol consumption (OR: 1.39; 95% CI 1.21, 1.60; P < 0.001) were positively associated with CRC risk. Conversely, genetically predicted coffee consumption was inversely related to CRC risk, with an OR (95% CI) of 0.80 (0.64, 0.99) (P = 0.037).

CONCLUSION

Genetically predicted alcohol use and consumption were risk factors for CRC while genetically predicted coffee consumption was a protective factor. Our findings highlight the effectiveness of keeping healthy dietary habits to prevent CRC. Further studies with more valid SNPs and CRC cases are needed. Validation of our findings is also recommended.

Whole-Exome Sequencing Analyses Support a Role of Vitamin D Metabolism in Ischemic Stroke

BACKGROUND

Ischemic stroke (IS) is a highly heritable trait, and genome-wide association studies have identified several commonly occurring susceptibility risk loci for this condition. However, there are limited data on the contribution of rare genetic variation to IS.

METHODS

We conducted an exome-wide study using whole-exome sequencing data from 152 058 UK Biobank participants, including 1777 IS cases. We performed single-variant analyses for rare variants and gene-based analyses for loss-of-function and deleterious missense rare variants. We validated these results through (1) gene-based testing using summary statistics from MEGASTROKE-a genome-wide association study of IS that included 67 162 IS cases and 454 450 controls, (2) gene-based testing using individual-level data from 1706 IS survivors, including 142 recurrent IS cases, enrolled in the VISP trial (Vitamin Intervention for Stroke Prevention); and (3) gene-based testing against neuroimaging phenotypes related to cerebrovascular disease using summary-level data from 42 310 UK Biobank participants with available magnetic resonance imaging data.

RESULTS

In single-variant association analyses, none of the evaluated variants were associated with IS at genome-wide significance levels (P<5×10-8). In the gene-based analysis focused on loss-of-function and deleterious missense variants, rare genetic variation at CYP2R1 was significantly associated with IS risk (P=2.6×10-6), exceeding the Bonferroni-corrected threshold for 16 074 tests (P<3.1×10-6). Validations analyses indicated that CYP2R1 was associated with IS risk in MEGASTROKE (gene-based test, P=0.003), with IS recurrence in the VISP trial (gene-based test, P=0.001) and with neuroimaging traits (white matter hyperintensity, mean diffusivity, and fractional anisotropy) in the UK Biobank neuroimaging study (all gene-based tests, P<0.05).

CONCLUSIONS

Because CYP2R1 plays an important role in vitamin D metabolism and existing observational evidence suggests an association between vitamin D levels and cerebrovascular disease, our results support a role of this pathway in the occurrence of IS.

Multitrait genome-wide analyses identify new susceptibility loci and candidate drugs to primary sclerosing cholangitis

Primary sclerosing cholangitis (PSC) is a rare autoimmune bile duct disease that is strongly associated with immune-mediated disorders. In this study, we implemented multitrait joint analyses to genome-wide association summary statistics of PSC and numerous clinical and epidemiological traits to estimate the genetic contribution of each trait and genetic correlations between traits and to identify new lead PSC risk-associated loci. We identified seven new loci that have not been previously reported and one new independent lead variant in the previously reported locus. Functional annotation and fine-mapping nominated several potential susceptibility genes such as MANBA and IRF5. Network-based in silico drug efficacy screening provided candidate agents for further study of pharmacological effect in PSC.

Eye-brain connections revealed by multimodal retinal and brain imaging genetics in the UK Biobank

As an anatomical extension of the brain, the retina of the eye is synaptically connected to the visual cortex, establishing physiological connections between the eye and the brain. Despite the unique opportunity retinal structures offer for assessing brain disorders, less is known about their relationship to brain structure and function. Here we present a systematic cross-organ genetic architecture analysis of eye-brain connections using retina and brain imaging endophenotypes. Novel phenotypic and genetic links were identified between retinal imaging biomarkers and brain structure and function measures derived from multimodal magnetic resonance imaging (MRI), many of which were involved in the visual pathways, including the primary visual cortex. In 65 genomic regions, retinal imaging biomarkers shared genetic influences with brain diseases and complex traits, 18 showing more genetic overlaps with brain MRI traits. Mendelian randomization suggests that retinal structures have bidirectional genetic causal links with neurological and neuropsychiatric disorders, such as Alzheimer's disease. Overall, cross-organ imaging genetics reveals a genetic basis for eye-brain connections, suggesting that the retinal images can elucidate genetic risk factors for brain disorders and disease-related changes in intracranial structure and function.

Genetic Susceptibility to Hepatocellular Carcinoma in Patients with Chronic Hepatitis Virus Infection

Hepatocellular carcinoma (HCC) is one of the leading causes of cancer-related deaths globally. The risk factors for HCC include chronic hepatitis B and C virus infections, excessive alcohol consumption, obesity, metabolic disease, and aflatoxin exposure. In addition to these viral and environmental risk factors, individual genetic predisposition is a major determinant of HCC risk. Familial clustering of HCC has been observed, and a hereditary factor likely contributes to the risk of HCC development. The familial aggregation may depend on a shared environment and genetic background as well as the interactions of environmental and genetic factors. Genome-wide association studies (GWASs) are one of the most practical tools for mapping the patterns of inheritance for the most common form of genomic variation, single nucleotide polymorphisms. This approach is practical for investigating genetic variants across the human genome, which is affected by thousands of common genetic variants that do not follow Mendelian inheritance. This review article summarizes the academic knowledge of GWAS-identified genetic loci and their association with HCC. We summarize the GWASs in accordance with various chronic hepatitis virus infection statuses. This genetic profiling could be used to identify candidate biomarkers to refine HCC screening and management by enabling individual risk-based personalization and stratification. A more comprehensive understanding of the genetic mechanisms underlying individual predisposition to HCC may lead to improvements in the prevention and early diagnosis of HCC and the development of effective treatment strategies.

Genome-wide genotype-serum proteome mapping provides insights into the cross-ancestry differences in cardiometabolic disease susceptibility

Identification of protein quantitative trait loci (pQTL) helps understand the underlying mechanisms of diseases and discover promising targets for pharmacological intervention. For most important class of drug targets, genetic evidence needs to be generalizable to diverse populations. Given that the majority of the previous studies were conducted in European ancestry populations, little is known about the protein-associated genetic variants in East Asians. Based on data-independent acquisition mass spectrometry technique, we conduct genome-wide association analyses for 304 unique proteins in 2,958 Han Chinese participants. We identify 195 genetic variant-protein associations. Colocalization and Mendelian randomization analyses highlight 60 gene-protein-phenotype associations, 45 of which (75%) have not been prioritized in Europeans previously. Further cross-ancestry analyses uncover key proteins that contributed to the differences in the obesity-induced diabetes and coronary artery disease susceptibility. These findings provide novel druggable proteins as well as a unique resource for the trans-ancestry evaluation of protein-targeted drug discovery.

Epigenetic and integrative cross-omics analyses of cerebral white matter hyperintensities on MRI

Cerebral white matter hyperintensities on MRI are markers of cerebral small vessel disease, a major risk factor for dementia and stroke. Despite the successful identification of multiple genetic variants associated with this highly heritable condition, its genetic architecture remains incompletely understood. More specifically, the role of DNA methylation has received little attention. We investigated the association between white matter hyperintensity burden and DNA methylation in blood at ∼450 000 cytosine-phosphate-guanine (CpG) sites in 9732 middle-aged to older adults from 14 community-based studies. Single CpG and region-based association analyses were carried out. Functional annotation and integrative cross-omics analyses were performed to identify novel genes underlying the relationship between DNA methylation and white matter hyperintensities. We identified 12 single CpG and 46 region-based DNA methylation associations with white matter hyperintensity burden. Our top discovery single CpG, cg24202936 (P = 7.6 × 10-8), was associated with F2 expression in blood (P = 6.4 × 10-5) and co-localized with FOLH1 expression in brain (posterior probability = 0.75). Our top differentially methylated regions were in PRMT1 and in CCDC144NL-AS1, which were also represented in single CpG associations (cg17417856 and cg06809326, respectively). Through Mendelian randomization analyses cg06809326 was putatively associated with white matter hyperintensity burden (P = 0.03) and expression of CCDC144NL-AS1 possibly mediated this association. Differentially methylated region analysis, joint epigenetic association analysis and multi-omics co-localization analysis consistently identified a role of DNA methylation near SH3PXD2A, a locus previously identified in genome-wide association studies of white matter hyperintensities. Gene set enrichment analyses revealed functions of the identified DNA methylation loci in the blood-brain barrier and in the immune response. Integrative cross-omics analysis identified 19 key regulatory genes in two networks related to extracellular matrix organization, and lipid and lipoprotein metabolism. A drug-repositioning analysis indicated antihyperlipidaemic agents, more specifically peroxisome proliferator-activated receptor-alpha, as possible target drugs for white matter hyperintensities. Our epigenome-wide association study and integrative cross-omics analyses implicate novel genes influencing white matter hyperintensity burden, which converged on pathways related to the immune response and to a compromised blood-brain barrier possibly due to disrupted cell-cell and cell-extracellular matrix interactions. The results also suggest that antihyperlipidaemic therapy may contribute to lowering risk for white matter hyperintensities possibly through protection against blood-brain barrier disruption.

Genetic Regulation of SMC Gene Expression and Splicing Predict Causal CAD Genes

BACKGROUND

Coronary artery disease (CAD) is the leading cause of death worldwide. Recent meta-analyses of genome-wide association studies have identified over 175 loci associated with CAD. The majority of these loci are in noncoding regions and are predicted to regulate gene expression. Given that vascular smooth muscle cells (SMCs) play critical roles in the development and progression of CAD, we aimed to identify the subset of the CAD loci associated with the regulation of transcription in distinct SMC phenotypes.

METHODS

We measured gene expression in SMCs isolated from the ascending aortas of 151 heart transplant donors of various genetic ancestries in quiescent or proliferative conditions and calculated the association of their expression and splicing with ~6.3 million imputed single-nucleotide polymorphism markers across the genome.

RESULTS

We identified 4910 expression and 4412 splicing quantitative trait loci (sQTLs) representing regions of the genome associated with transcript abundance and splicing. A total of 3660 expression quantitative trait loci (eQTLs) had not been observed in the publicly available Genotype-Tissue Expression dataset. Further, 29 and 880 eQTLs were SMC-specific and sex-biased, respectively. We made these results available for public query on a user-friendly website. To identify the effector transcript(s) regulated by CAD loci, we used 4 distinct colocalization approaches. We identified 84 eQTL and 164 sQTL that colocalized with CAD loci, highlighting the importance of genetic regulation of mRNA splicing as a molecular mechanism for CAD genetic risk. Notably, 20% and 35% of the eQTLs were unique to quiescent or proliferative SMCs, respectively. One CAD locus colocalized with a sex-specific eQTL (TERF2IP), and another locus colocalized with SMC-specific eQTL (ALKBH8). The most significantly associated CAD locus, 9p21, was an sQTL for the long noncoding RNA CDKN2B-AS1, also known as ANRIL, in proliferative SMCs.

CONCLUSIONS

Collectively, our results provide evidence for the molecular mechanisms of genetic susceptibility to CAD in distinct SMC phenotypes.

Development and validation of genome-wide polygenic risk scores for predicting breast cancer incidence in Japanese females: a population-based case-cohort study

PURPOSE

This study aimed to develop an ancestry-specific polygenic risk scores (PRSs) for the prediction of breast cancer events in Japanese females and validate it in a longitudinal cohort study.

METHODS

Using publicly available summary statistics of female breast cancer genome-wide association study (GWAS) of Japanese and European ancestries, we, respectively, developed 31 candidate genome-wide PRSs using pruning and thresholding (P + T) and LDpred methods with varying parameters. Among the candidate PRS models, the best model was selected using a case-cohort dataset (63 breast cancer cases and 2213 sub-cohorts of Japanese females during a median follow-up of 11.9 years) according to the maximal predictive ability by Harrell's C-statistics. The best-performing PRS for each derivation GWAS was evaluated in another independent case-cohort dataset (260 breast cancer cases and 7845 sub-cohorts of Japanese females during a median follow-up of 16.9 years).

RESULTS

For the best PRS model involving 46,861 single nucleotide polymorphisms (SNPs; P + T method with P_T = 0.05 and R² = 0.2) derived from Japanese-ancestry GWAS, the Harrell's C-statistic was 0.598 ± 0.018 in the evaluation dataset. The age-adjusted hazard ratio for breast cancer in females with the highest PRS quintile compared with those in the lowest PRS quintile was 2.47 (95% confidence intervals, 1.64-3.70). The PRS constructed using Japanese-ancestry GWAS demonstrated better predictive performance for breast cancer in Japanese females than that using European-ancestry GWAS (Harrell's C-statistics 0.598 versus 0.586).

CONCLUSION

This study developed a breast cancer PRS for Japanese females and demonstrated the usefulness of the PRS for breast cancer risk stratification.

Genetic variation in TERT modifies the risk of hepatocellular carcinoma in alcohol-related cirrhosis: results from a genome-wide case-control study

OBJECTIVE

Hepatocellular carcinoma (HCC) often develops in patients with alcohol-related cirrhosis at an annual risk of up to 2.5%. Some host genetic risk factors have been identified but do not account for the majority of the variance in occurrence. This study aimed to identify novel susceptibility loci for the development of HCC in people with alcohol related cirrhosis.

DESIGN

Patients with alcohol-related cirrhosis and HCC (cases: n=1214) and controls without HCC (n=1866), recruited from Germany, Austria, Switzerland, Italy and the UK, were included in a two-stage genome-wide association study using a case-control design. A validation cohort of 1520 people misusing alcohol but with no evidence of liver disease was included to control for possible association effects with alcohol misuse. Genotyping was performed using the InfiniumGlobal Screening Array (V.24v2, Illumina) and the OmniExpress Array (V.24v1-0a, Illumina).

RESULTS

Associations with variants rs738409 in PNPLA3 and rs58542926 in TM6SF2 previously associated with an increased risk of HCC in patients with alcohol-related cirrhosis were confirmed at genome-wide significance. A novel locus rs2242652(A) in TERT (telomerase reverse transcriptase) was also associated with a decreased risk of HCC, in the combined meta-analysis, at genome-wide significance (p=6.41×10-9, OR=0.61 (95% CI 0.52 to 0.70). This protective association remained significant after correction for sex, age, body mass index and type 2 diabetes (p=7.94×10-5, OR=0.63 (95% CI 0.50 to 0.79). Carriage of rs2242652(A) in TERT was associated with an increased leucocyte telomere length (p=2.12×10-44).

CONCLUSION

This study identifies rs2242652 in TERT as a novel protective factor for HCC in patients with alcohol-related cirrhosis.

MRBEE: A novel bias-corrected multivariable Mendelian Randomization method

Mendelian Randomization (MR) has been widely applied to infer causality of exposures on outcomes in the genome wide association (GWAS) era. Existing approaches are often subject to biases from multiple sources including weak instruments, sample overlap, and measurement error. We introduce MRBEE, a computationally efficient multivariable MR method that can correct for all known biases simultaneously, which is demonstrated in theory, simulations, and real data analysis. In comparison, all existing MR methods are biased. In two independent real data analyses, we observed that the causal effect of BMI on coronary artery disease risk is completely mediated by blood pressure, and that existing MR methods drastically underestimate the causal effect of cannabis use disorder on schizophrenia risk compared to MRBEE. We demonstrate that MRBEE can be a useful tool in studying causality between multiple risk factors and a disease outcome, especially as more GWAS summary statistics are being made publicly available.

Recent advances of machine learning applications in human gut microbiota study: from observational analysis toward causal inference and clinical intervention

Statistical methods, especially machine learning, learning(ML), are pivotal for the analyses of large data generated by multiomics human gut microbiota study. These analyses lead to the discovery of microbe-disease associations. Furthermore, recent efforts for more data transparency and accessible analytical tools improved data availability and study reproducibility. Our recent accumulated knowledge on microbe-disease associations brings light to the next questions: what is the role of microbes in disease progression and how can we apply our knowledge of microbiome in clinical settings? Here, we introduce recent studies that implemented ML to answer the questions of causal inference and clinical translation.

Evidence of a causal effect of genetic tendency to gain muscle mass on uterine leiomyomata

Uterine leiomyomata (UL) are the most common tumours of the female genital tract and the primary cause of surgical removal of the uterus. Genetic factors contribute to UL susceptibility. To add understanding to the heritable genetic risk factors, we conduct a genome-wide association study (GWAS) of UL in up to 426,558 European women from FinnGen and a previous UL meta-GWAS. In addition to the 50 known UL loci, we identify 22 loci that have not been associated with UL in prior studies. UL-associated loci harbour genes enriched for development, growth, and cellular senescence. Of particular interest are the smooth muscle cell differentiation and proliferation-regulating genes functioning on the myocardin-cyclin dependent kinase inhibitor 1 A pathway. Our results further suggest that genetic predisposition to increased fat-free mass may be causally related to higher UL risk, underscoring the involvement of altered muscle tissue biology in UL pathophysiology. Overall, our findings add to the understanding of the genetic pathways underlying UL, which may aid in developing novel therapeutics.

Genome-wide association study of age at menarche in the Taiwan Biobank suggests NOL4 as a novel associated gene

Sexual maturation is a complex physiological process that involves multiple variables, such as genetic and environmental factors. Among females, age at menarche (AM) is a critical milestone for sexual maturation. This study aimed to identify genetic markers of AM using nationwide population cohort data in Taiwan. Females with self-reported AM between 10 and 16 years (N = 39,827) were eligible for the final analysis. To identify genetic signals related to AM, we conducted a genome-wide association study using a linear regression model and split-half meta-analysis method to verify our findings. The Functional Mapping and Annotation web-based platform was used for positional mapping and gene-based and gene-set analyses. The meta-analysis identified four significant loci, i.e., LIN28B (pooled P = 1.39 × 10^-21), NOL4 (pooled P = 8.94 × 10^-9), GPR45 (pooled P = 4.19 × 10^-11), and LOC105373831 (pooled P = 4.37 × 10^-8), that were associated with AM. MAGMA gene-based analysis revealed that LIN28B (P = 1.13 × 10^-8), NOL4 (P = 2.27 × 10^-7), RXRG (P = 4.34 × 10-⁷), ETV5 (P = 1.75 × 10^-6), and HACE1 (P = 1.82 × 10^-6) were significantly associated with AM, while the gene-set analysis identified a significantly enriched pathway involving mTOR signaling complex (FDR corrected P = 1.28 × 10^-2). The results replicated evidence for several genetic markers associated with AM in the Taiwanese female population. Our analysis identified a novel locus (rs7239368) in NOL4 associated with AM (β = 0.051 ± 0.009 years, pooled P = 8.94 × 10^-9), whereas additional research is needed to validate its molecular role in sexual maturation.