Protocol for genome-wide association study of human blood metabolites.

Genetic variations influence the levels of blood metabolites. We present analytical pipelines for assessing genetic influences on human blood metabolites. We describe steps for the normalization of metabolome data, genome-wide association studies, and the identification of metabolite quantitative trait loci (mQTLs). We then detail procedures for functional enrichment analysis of mQTLs. This protocol could be applicable to other quantitative traits, such as clinical measurements or proteome data. For complete details on the use and execution of this protocol, please refer to Iwasaki et al.1.

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

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

Poor accuracy and sustainability of the first-step FIB4 EASL pathway for stratifying steatotic liver disease risk in the general population.

BACKGROUND AND AIMS: The European Association for the Study of the Liver introduced a clinical pathway (EASL CP) for screening significant/advanced fibrosis in people at risk of steatotic liver disease (SLD). We assessed the performance of the first-step FIB4 EASL CP in the general population across different SLD risk groups (MASLD, Met-ALD and ALD) and various age classes. METHODS: We analysed a total of 3372 individuals at risk of SLD from the 2017-2018 National Health and Nutrition Examination Survey (NHANES17-18), projected to 152.3 million U.S. adults, 300,329 from the UK Biobank (UKBB) and 57,644 from the Biobank Japan (BBJ). We assessed liver stiffness measurement (LSM) ≥8 kPa and liver-related events occurring within 3 and 10 years (3/10 year-LREs) as outcomes. We defined MASLD, MetALD, and ALD according to recent international recommendations. RESULTS: FIB4 sensitivity for LSM ≥ 8 kPa was low (27.7%), but it ranged approximately 80%-90% for 3-year LREs. Using FIB4, 22%-57% of subjects across the three cohorts were identified as candidates for vibration-controlled transient elastography (VCTE), which was mostly avoidable (positive predictive value of FIB4 ≥ 1.3 for LSM ≥ 8 kPa ranging 9.5%-13% across different SLD categories). Sensitivity for LSM ≥ 8 kPa and LREs increased with increasing alcohol intake (ALD>MetALD>MASLD) and age classes. For individuals aged ≥65 years, using the recommended age-adjusted FIB4 cut-off (≥2) substantially reduced sensitivity for LSM ≥ 8 kPa and LREs. CONCLUSIONS: The first-step FIB4 EASL CP is poorly accurate and feasible for individuals at risk of SLD in the general population. It is crucial to enhance the screening strategy with a first-step approach able to reduce unnecessary VCTEs and optimise their yield.

A large-scale microRNA transcriptome-wide association study identifies two susceptibility microRNAs, miR-1307-5p and miR-192-3p, for colorectal cancer risk.

Transcriptome-wide association studies (TWAS) have identified many putative susceptibility genes for colorectal cancer (CRC) risk. However, susceptibility miRNAs, critical dysregulators of gene expression, remain unexplored. We genotyped DNA samples from 313 CRC East Asian patients and performed small RNA sequencing in their normal colon tissues distant from tumors to build genetic models for predicting miRNA expression. We applied these models and data from genome-wide association studies (GWAS) including 23 942 cases and 217 267 controls of East Asian ancestry to investigate associations of predicted miRNA expression with CRC risk. Perturbation experiments separately by promoting and inhibiting miRNAs expressions and further in vitro assays in both SW480 and HCT116 cells were conducted. At a Bonferroni-corrected threshold of P < 4.5 × 10-4, we identified two putative susceptibility miRNAs, miR-1307-5p and miR-192-3p, located in regions more than 500 kb away from any GWAS-identified risk variants in CRC. We observed that a high predicted expression of miR-1307-5p was associated with increased CRC risk, while a low predicted expression of miR-192-3p was associated with increased CRC risk. Our experimental results further provide strong evidence of their susceptible roles by showing that miR-1307-5p and miR-192-3p play a regulatory role, respectively, in promoting and inhibiting CRC cell proliferation, migration, and invasion, which was consistently observed in both SW480 and HCT116 cells. Our study provides additional insights into the biological mechanisms underlying CRC development.

Characterizing prostate cancer risk through multi-ancestry genome-wide discovery of 187 novel risk variants.

The transferability and clinical value of genetic risk scores (GRSs) across populations remain limited due to an imbalance in genetic studies across ancestrally diverse populations. Here we conducted a multi-ancestry genome-wide association study of 156,319 prostate cancer cases and 788,443 controls of European, African, Asian and Hispanic men, reflecting a 57% increase in the number of non-European cases over previous prostate cancer genome-wide association studies. We identified 187 novel risk variants for prostate cancer, increasing the total number of risk variants to 451. An externally replicated multi-ancestry GRS was associated with risk that ranged from 1.8 (per standard deviation) in African ancestry men to 2.2 in European ancestry men. The GRS was associated with a greater risk of aggressive versus non-aggressive disease in men of African ancestry (P = 0.03). Our study presents novel prostate cancer susceptibility loci and a GRS with effective risk stratification across ancestry groups.

East Asian-specific and cross-ancestry genome-wide meta-analyses provide mechanistic insights into peptic ulcer disease.

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.

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

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

Boosting the power of genome-wide association studies within and across ancestries by using polygenic scores.

Genome-wide association studies (GWASs) have been mostly conducted in populations of European ancestry, which currently limits the transferability of their findings to other populations. Here, we show, through theory, simulations and applications to real data, that adjustment of GWAS analyses for polygenic scores (PGSs) increases the statistical power for discovery across all ancestries. We applied this method to analyze seven traits available in three large biobanks with participants of East Asian ancestry (n = 340,000 in total) and report 139 additional associations across traits. We also present a two-stage meta-analysis strategy whereby, in contributing cohorts, a PGS-adjusted GWAS is rerun using PGSs derived from a first round of a standard meta-analysis. On average, across traits, this approach yields a 1.26-fold increase in the number of detected associations (range 1.07- to 1.76-fold increase). Altogether, our study demonstrates the value of using PGSs to increase the power of GWASs in underrepresented populations and promotes such an analytical strategy for future GWAS meta-analyses.

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

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

Hyperfructosemia in sleep disordered breathing: metabolome analysis of Nagahama study.

Sleep disordered breathing (SDB), mainly obstructive sleep apnea (OSA), constitutes a major health problem due to the large number of patients. Intermittent hypoxia caused by SDB induces alterations in metabolic function. Nevertheless, metabolites characteristic for SDB are largely unknown. In this study, we performed gas chromatography-mass spectrometry-based targeted metabolome analysis using data from The Nagahama Study (n = 6373). SDB-related metabolites were defined based on their variable importance score in orthogonal partial least squares discriminant analysis and fold changes in normalized peak-intensity levels between moderate-severe SDB patients and participants without SDB. We identified 20 metabolites as SDB-related, and interestingly, these metabolites were frequently included in pathways related to fructose. Multivariate analysis revealed that moderate-severe SDB was a significant factor for increased plasma fructose levels (ß = 0.210, P = 0.006, generalized linear model) even after the adjustment of confounding factors. We further investigated changes in plasma fructose levels after continuous positive airway pressure (CPAP) treatment using samples from patients with OSA (n = 60) diagnosed by polysomnography at Kyoto University Hospital, and found that patients with marked hypoxemia exhibited prominent hyperfructosemia and their plasma fructose levels lowered after CPAP treatment. These data suggest that hyperfructosemia is the abnormality characteristic to SDB, which can be reduced by CPAP treatment.

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

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

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

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

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

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

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

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

Improved genetic prediction of the risk of knee osteoarthritis using the risk factor-based polygenic score.

BACKGROUND: Polygenic risk score (PRS) analysis is used to predict disease risk. Although PRS has been shown to have great potential in improving clinical care, PRS accuracy assessment has been mainly focused on European ancestry. This study aimed to develop an accurate genetic risk score for knee osteoarthritis (OA) using a multi-population PRS and leveraging a multi-trait PRS in the Japanese population. METHODS: We calculated PRS using PRS-CS-auto, derived from genome-wide association study (GWAS) summary statistics for knee OA in the Japanese population (same ancestry) and multi-population. We further identified risk factor traits for which PRS could predict knee OA and subsequently developed an integrated PRS based on multi-trait analysis of GWAS (MTAG), including genetically correlated risk traits. PRS performance was evaluated in participants of the Nagahama cohort study who underwent radiographic evaluation of the knees (n = 3,279). PRSs were incorporated into knee OA integrated risk models along with clinical risk factors. RESULTS: A total of 2,852 genotyped individuals were included in the PRS analysis. The PRS based on Japanese knee OA GWAS was not associated with knee OA (p = 0.228). In contrast, PRS based on multi-population knee OA GWAS showed a significant association with knee OA (p = 6.7 × 10-5, odds ratio (OR) per standard deviation = 1.19), whereas PRS based on MTAG of multi-population knee OA, along with risk factor traits such as body mass index GWAS, displayed an even stronger association with knee OA (p = 5.4 × 10-7, OR = 1.24). Incorporating this PRS into traditional risk factors improved the predictive ability of knee OA (area under the curve, 74.4% to 74.7%; p = 0.029). CONCLUSIONS: This study showed that multi-trait PRS based on MTAG, combined with traditional risk factors, and using large sample size multi-population GWAS, significantly improved predictive accuracy for knee OA in the Japanese population, even when the sample size of GWAS of the same ancestry was small. To the best of our knowledge, this is the first study to show a statistically significant association between the PRS and knee OA in a non-European population. TRIAL REGISTRATION: No. C278.

Increased PDGFRB and NF-κB signaling caused by highly prevalent somatic mutations in intracranial aneurysms.

Intracranial aneurysms (IAs) are a high-risk factor for life-threatening subarachnoid hemorrhage. Their etiology, however, remains mostly unknown at present. We conducted screening for sporadic somatic mutations in 65 IA tissues (54 saccular and 11 fusiform aneurysms) and paired blood samples by whole-exome and targeted deep sequencing. We identified sporadic mutations in multiple signaling genes and examined their impact on downstream signaling pathways and gene expression in vitro and an arterial dilatation model in mice in vivo. We identified 16 genes that were mutated in at least one IA case and found that these mutations were highly prevalent (92%: 60 of 65 IAs) among all IA cases examined. In particular, mutations in six genes (PDGFRB, AHNAK, OBSCN, RBM10, CACNA1E, and OR5P3), many of which are linked to NF-κB signaling, were found in both fusiform and saccular IAs at a high prevalence (43% of all IA cases examined). We found that mutant PDGFRBs constitutively activated ERK and NF-κB signaling, enhanced cell motility, and induced inflammation-related gene expression in vitro. Spatial transcriptomics also detected similar changes in vessels from patients with IA. Furthermore, virus-mediated overexpression of a mutant PDGFRB induced a fusiform-like dilatation of the basilar artery in mice, which was blocked by systemic administration of the tyrosine kinase inhibitor sunitinib. Collectively, this study reveals a high prevalence of somatic mutations in NF-κB signaling pathway-related genes in both fusiform and saccular IAs and opens a new avenue of research for developing pharmacological interventions.

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.

Genome-wide Association Studies Categorized by Class of Antihypertensive Drugs Reveal Complex Pathogenesis of Hypertension with Drug Resistance.

Resistant hypertension is defined as uncontrolled blood pressure (BP) despite the use of three or more antihypertensive drugs of different classes. Although genetic factors may greatly contribute to hypertension with resistance to multiple drug classes, more than for general hypertension, its pathogenesis remains unknown. To reveal the genetic background of resistant hypertension, we categorized 32,239 patients whose data were obtained from the BioBank Japan Project, by prescription of 7 classes of antihypertensive drugs and performed genome-wide association studies (GWAS). Our GWAS identified four loci with significant association (P < 5 × 10-8 ): rs6445583 in CACNA1D and rs12308051 in the intergenic region on chromosome 12 for angiotensin II receptor blockers, rs35497065 in FOXA3 for calcium channel blockers, and rs11066280 in HECTD4 for αß-blockers. Because these loci are known to be susceptibility loci for hypertension and/or BP, our results indicate that resistant hypertension is caused by a combination of excessive BP and drug resistance to each antihypertensive pharmacological class. Furthermore, to investigate the genetic difference between BP traits and the treatment effectiveness of antihypertensive drugs, we performed gene-set analysis and calculated the genetic correlation continuously. Most of the genetic factors were in common between BP traits and antihypertensive effectiveness, but it seems that the genetic architecture of the drug response to antihypertensive treatment is more complicated than BP traits. This corresponds to the well-known mosaic theory of hypertension. Our findings reveal the complex pathogenesis of hypertension with resistance to multiple classes of antihypertensive drugs.

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.

Helicobacter pylori, Homologous-Recombination Genes, and Gastric Cancer.

BACKGROUND: Helicobacter pylori infection is a well-known risk factor for gastric cancer. However, the contribution of germline pathogenic variants in cancer-predisposing genes and their effect, when combined with H. pylori infection, on the risk of gastric cancer has not been widely evaluated. METHODS: We evaluated the association between germline pathogenic variants in 27 cancer-predisposing genes and the risk of gastric cancer in a sample of 10,426 patients with gastric cancer and 38,153 controls from BioBank Japan. We also assessed the combined effect of pathogenic variants and H. pylori infection status on the risk of gastric cancer and calculated the cumulative risk in 1433 patients with gastric cancer and 5997 controls from the Hospital-based Epidemiologic Research Program at Aichi Cancer Center (HERPACC). RESULTS: Germline pathogenic variants in nine genes (APC, ATM, BRCA1, BRCA2, CDH1, MLH1, MSH2, MSH6, and PALB2) were associated with the risk of gastric cancer. We found an interaction between H. pylori infection and pathogenic variants in homologous-recombination genes with respect to the risk of gastric cancer in the sample from HERPACC (relative excess risk due to the interaction, 16.01; 95% confidence interval [CI], 2.22 to 29.81; P = 0.02). At 85 years of age, persons with H. pylori infection and a pathogenic variant had a higher cumulative risk of gastric cancer than noncarriers infected with H. pylori (45.5% [95% CI, 20.7 to 62.6] vs. 14.4% [95% CI, 12.2 to 16.6]). CONCLUSIONS: H. pylori infection modified the risk of gastric cancer associated with germline pathogenic variants in homologous-recombination genes. (Funded by the Japan Agency for Medical Research and Development and others.).