Telomere length and risk of cirrhosis, hepatocellular carcinoma, and cholangiocarcinoma in 63,272 individuals from the general population

BACKGROUND AND AIMS

Inherited short telomeres are associated with a risk of liver disease, whereas longer telomeres predispose to cancer. The association between telomere length and risk of HCC and cholangiocarcinoma remains unknown.

APPROACH AND RESULTS

We measured leukocyte telomere length using multiplex PCR in 63,272 individuals from the Danish general population. Telomere length and plasma ALT concentration were not associated (β = 4 ×10 -6 , p -value = 0.06) in a linear regression model, without any signs of a nonlinear relationship. We tested the association between telomere length and risk of cirrhosis, HCC, and cholangiocarcinoma using Cox regression. During a median follow-up of 11 years, 241, 76, and 112 individuals developed cirrhosis, HCC, and cholangiocarcinoma, respectively. Telomere length and risk of cirrhosis were inversely and linearly associated ( p -value = 0.004, p for nonlinearity = 0.27). Individuals with telomeres in the shortest vs. longest quartile had a 2.25-fold higher risk of cirrhosis. Telomere length and risk of HCC were nonlinearly associated ( p -value = 0.009, p -value for nonlinearity = 0.01). This relationship resembled an inverted J-shape, with the highest risk observed in individuals with short telomeres. Individuals with telomeres in the shortest versus longest quartile had a 2.29-fold higher risk of HCC. Telomere length was inversely and linearly associated with the risk of cholangiocarcinoma. Individuals with telomeres in the shortest versus longest quartile had a 1.86-fold higher risk of cholangiocarcinoma.

CONCLUSIONS

Shorter telomere length is associated with a higher risk of cirrhosis, HCC, and cholangiocarcinoma.

Worldwide trends in underweight and obesity from 1990 to 2022: a pooled analysis of 3663 population-representative studies with 222 million children, adolescents, and adults

BACKGROUND

Underweight and obesity are associated with adverse health outcomes throughout the life course. We estimated the individual and combined prevalence of underweight or thinness and obesity, and their changes, from 1990 to 2022 for adults and school-aged children and adolescents in 200 countries and territories.

METHODS

We used data from 3663 population-based studies with 222 million participants that measured height and weight in representative samples of the general population. We used a Bayesian hierarchical model to estimate trends in the prevalence of different BMI categories, separately for adults (age ≥20 years) and school-aged children and adolescents (age 5-19 years), from 1990 to 2022 for 200 countries and territories. For adults, we report the individual and combined prevalence of underweight (BMI <18·5 kg/m2) and obesity (BMI ≥30 kg/m2). For school-aged children and adolescents, we report thinness (BMI <2 SD below the median of the WHO growth reference) and obesity (BMI >2 SD above the median).

FINDINGS

From 1990 to 2022, the combined prevalence of underweight and obesity in adults decreased in 11 countries (6%) for women and 17 (9%) for men with a posterior probability of at least 0·80 that the observed changes were true decreases. The combined prevalence increased in 162 countries (81%) for women and 140 countries (70%) for men with a posterior probability of at least 0·80. In 2022, the combined prevalence of underweight and obesity was highest in island nations in the Caribbean and Polynesia and Micronesia, and countries in the Middle East and north Africa. Obesity prevalence was higher than underweight with posterior probability of at least 0·80 in 177 countries (89%) for women and 145 (73%) for men in 2022, whereas the converse was true in 16 countries (8%) for women, and 39 (20%) for men. From 1990 to 2022, the combined prevalence of thinness and obesity decreased among girls in five countries (3%) and among boys in 15 countries (8%) with a posterior probability of at least 0·80, and increased among girls in 140 countries (70%) and boys in 137 countries (69%) with a posterior probability of at least 0·80. The countries with highest combined prevalence of thinness and obesity in school-aged children and adolescents in 2022 were in Polynesia and Micronesia and the Caribbean for both sexes, and Chile and Qatar for boys. Combined prevalence was also high in some countries in south Asia, such as India and Pakistan, where thinness remained prevalent despite having declined. In 2022, obesity in school-aged children and adolescents was more prevalent than thinness with a posterior probability of at least 0·80 among girls in 133 countries (67%) and boys in 125 countries (63%), whereas the converse was true in 35 countries (18%) and 42 countries (21%), respectively. In almost all countries for both adults and school-aged children and adolescents, the increases in double burden were driven by increases in obesity, and decreases in double burden by declining underweight or thinness.

INTERPRETATION

The combined burden of underweight and obesity has increased in most countries, driven by an increase in obesity, while underweight and thinness remain prevalent in south Asia and parts of Africa. A healthy nutrition transition that enhances access to nutritious foods is needed to address the remaining burden of underweight while curbing and reversing the increase in obesity.

FUNDING

UK Medical Research Council, UK Research and Innovation (Research England), UK Research and Innovation (Innovate UK), and European Union.

Identification of genetically predicted DNA methylation markers associated with non-small cell lung cancer risk among 34,964 cases and 448,579 controls

BACKGROUND

Although the associations between genetic variations and lung cancer risk have been explored, the epigenetic consequences of DNA methylation in lung cancer development are largely unknown. Here, the genetically predicted DNA methylation markers associated with non-small cell lung cancer (NSCLC) risk by a two-stage case-control design were investigated.

METHODS

The genetic prediction models for methylation levels based on genetic and methylation data of 1595 subjects from the Framingham Heart Study were established. The prediction models were applied to a fixed-effect meta-analysis of screening data sets with 27,120 NSCLC cases and 27,355 controls to identify the methylation markers, which were then replicated in independent data sets with 7844 lung cancer cases and 421,224 controls. Also performed was a multi-omics functional annotation for the identified CpGs by integrating genomics, epigenomics, and transcriptomics and investigation of the potential regulation pathways.

RESULTS

Of the 29,894 CpG sites passing the quality control, 39 CpGs associated with NSCLC risk (Bonferroni-corrected p ≤ 1.67 × 10-6 ) were originally identified. Of these, 16 CpGs remained significant in the validation stage (Bonferroni-corrected p ≤ 1.28 × 10-3 ), including four novel CpGs. Multi-omics functional annotation showed nine of 16 CpGs were potentially functional biomarkers for NSCLC risk. Thirty-five genes within a 1-Mb window of 12 CpGs that might be involved in regulatory pathways of NSCLC risk were identified.

CONCLUSIONS

Sixteen promising DNA methylation markers associated with NSCLC were identified. Changes of the methylation level at these CpGs might influence the development of NSCLC by regulating the expression of genes nearby.

PLAIN LANGUAGE SUMMARY

The epigenetic consequences of DNA methylation in lung cancer development are still largely unknown. This study used summary data of large-scale genome-wide association studies to investigate the associations between genetically predicted levels of methylation biomarkers and non-small cell lung cancer risk at the first time. This study looked at how well larotrectinib worked in adult patients with sarcomas caused by TRK fusion proteins. These findings will provide a unique insight into the epigenetic susceptibility mechanisms of lung cancer.

Lung Cancer in Ever- and Never-Smokers: Findings from Multi-Population GWAS Studies

BACKGROUND

Clinical, molecular, and genetic epidemiology studies displayed remarkable differences between ever- and never-smoking lung cancer.

METHODS

We conducted a stratified multi-population (European, East Asian, and African descent) association study on 44,823 ever-smokers and 20,074 never-smokers to identify novel variants that were missed in the non-stratified analysis. Functional analysis including expression quantitative trait loci (eQTL) colocalization and DNA damage assays, and annotation studies were conducted to evaluate the functional roles of the variants. We further evaluated the impact of smoking quantity on lung cancer risk for the variants associated with ever-smoking lung cancer.

RESULTS

Five novel independent loci, GABRA4, intergenic region 12q24.33, LRRC4C, LINC01088, and LCNL1 were identified with the association at two or three populations (P < 5 × 10-8). Further functional analysis provided multiple lines of evidence suggesting the variants affect lung cancer risk through excessive DNA damage (GABRA4) or cis-regulation of gene expression (LCNL1). The risk of variants from 12 independent regions, including the well-known CHRNA5, associated with ever-smoking lung cancer was evaluated for never-smokers, light-smokers (packyear ≤ 20), and moderate-to-heavy-smokers (packyear > 20). Different risk patterns were observed for the variants among the different groups by smoking behavior.

CONCLUSIONS

We identified novel variants associated with lung cancer in only ever- or never-smoking groups that were missed by prior main-effect association studies.

IMPACT

Our study highlights the genetic heterogeneity between ever- and never-smoking lung cancer and provides etiologic insights into the complicated genetic architecture of this deadly cancer.

Elevated plasma triglycerides increase risk of psoriasis: A cohort and Mendelian randomization study

BACKGROUND

It is increasingly clear that triglyceride-rich lipoproteins are proinflammatory and cause low-grade systemic inflammation. However, it is currently unknown whether elevated plasma triglycerides are causally related to development of psoriasis, a skin disorder driven by chronic inflammation.

OBJECTIVE

To determine if elevated plasma triglycerides are associated with increased risk of psoriasis in observational and Mendelian randomization analysis.

METHODS

Consecutive individuals from the Copenhagen General Population Study (CGPS) were included. We used plasma triglycerides (n = 108,043) and a weighted triglyceride allele score (n = 92,579) on nine known triglyceride-altering genetic variants. Genetic results were replicated in 337,159 individuals from the UK biobank. Psoriasis was ICD10-code hospital contact in main analyses, and prescription of topical antipsoriatics for mild psoriasis in sensitivity analysis.

RESULTS

During a median 9.3 years (0.1-15.1) of follow-up (from 2003-2015 through 2018), 855 (1%) individuals were diagnosed with psoriasis by ICD-10 in observational analysis and 772 (1%) in Mendelian randomization analysis. In observational analysis, multivariable adjusted hazard ratio for psoriasis by ICD-10 were 1.26 (95% CI:1.15-1.39) per doubling in plasma triglycerides with a corresponding causal, genetic risk ratio of 2.10 (1.30-3.38). Causality was confirmed in the UK biobank. Results were similar but slightly attenuated when we used topical antipsoriatics prescription for mild psoriasis.

CONCLUSION

Elevated plasma triglycerides are associated with increased risk of psoriasis in observational and Mendelian randomization analysis.

Prognostic Value of Pre-treatment Plasma C-Reactive Protein in Patients with Early-Stage Breast Cancer

BACKGROUND

Breast cancer incidence is now the highest among all cancers and accountable for 6.6% of all cancer-related deaths worldwide. Studies of the prognostic utility of plasma CRP measurement in early-stage breast cancer have given discrepant results.

METHODS

We identified 6942 patients in the Danish Breast Cancer Cooperative Group database with early-stage breast cancer diagnosed between 2002 and 2016 who had a measure of pre-treatment plasma CRP. Outcomes were recurrence-free interval and survival for a period up to ten years. We analyzed associations with plasma CRP using Fine-Gray proportional sub-distribution hazards model with recurrence-free interval. Data on plasma CRP were analyzed per doubling of concentration and in relation to CRP levels of < 3 mg/L; 3-10 mg/L, and > 10 mg/L and stratified according to standard clinical parameters in sensitivity analyses.

RESULTS

A doubling of the plasma CRP concentration was associated with increased risk of recurrence (multivariate adjusted HR = 1.05, 95% CI 1.01;1.08) and shorter survival (HR = 1.13, 1.09;1.16) in multivariate analyses. Survival was shorter in patients with plasma CRP levels of 3-10 and > 10 vs. < 3 mg/L, with multivariate adjusted HRs of 1.30 (1.17;1.45) and 1.65 (1.39;1.95), respectively.

CONCLUSIONS

Elevated plasma CRP measured before treatment in patients with early-stage breast cancer is an independent biomarker of increased risk of recurrence and early death.

IMPACT

CRP measures before treatment might be used to individualize follow-up of patients with early-stage breast cancer.

CYP2A6 Activity and Cigarette Consumption Interact in Smoking-Related Lung Cancer Susceptibility

Cigarette smoke, containing both nicotine and carcinogens, causes lung cancer. However, not all smokers develop lung cancer, highlighting the importance of the interaction between host susceptibility and environmental exposure in tumorigenesis. Here, we aimed to delineate the interaction between metabolizing ability of tobacco carcinogens and smoking intensity in mediating genetic susceptibility to smoking-related lung tumorigenesis. Single-variant and gene-based associations of 43 tobacco carcinogen-metabolizing genes with lung cancer were analyzed using summary statistics and individual-level genetic data, followed by causal inference of Mendelian randomization, mediation analysis, and structural equation modeling. Cigarette smoke-exposed cell models were used to detect gene expression patterns in relation to specific alleles. Data from the International Lung Cancer Consortium (29,266 cases and 56,450 controls) and UK Biobank (2,155 cases and 376,329 controls) indicated that the genetic variant rs56113850 C>T located in intron 4 of CYP2A6 was significantly associated with decreased lung cancer risk among smokers (OR = 0.88, 95% confidence interval = 0.85-0.91, P = 2.18 × 10-16), which might interact (Pinteraction = 0.028) with and partially be mediated (ORindirect = 0.987) by smoking status. Smoking intensity accounted for 82.3% of the effect of CYP2A6 activity on lung cancer risk but entirely mediated the genetic effect of rs56113850. Mechanistically, the rs56113850 T allele rescued the downregulation of CYP2A6 caused by cigarette smoke exposure, potentially through preferential recruitment of transcription factor helicase-like transcription factor. Together, this study provides additional insights into the interplay between host susceptibility and carcinogen exposure in smoking-related lung tumorigenesis.

SIGNIFICANCE

The causal pathway connecting CYP2A6 genetic variability and activity, cigarette consumption, and lung cancer susceptibility in smokers highlights the need for behavior modification interventions based on host susceptibility for cancer prevention.

Differences in polygenic score distributions in European ancestry populations: implications for breast cancer risk prediction

The 313-variant polygenic risk score (PRS313) provides a promising tool for breast cancer risk prediction. However, evaluation of the PRS313 across different European populations which could influence risk estimation has not been performed. Here, we explored the distribution of PRS313 across European populations using genotype data from 94,072 females without breast cancer, of European-ancestry from 21 countries participating in the Breast Cancer Association Consortium (BCAC) and 225,105 female participants from the UK Biobank. The mean PRS313 differed markedly across European countries, being highest in south-eastern Europe and lowest in north-western Europe. Using the overall European PRS313 distribution to categorise individuals leads to overestimation and underestimation of risk in some individuals from south-eastern and north-western countries, respectively. Adjustment for principal components explained most of the observed heterogeneity in mean PRS. Country-specific PRS distributions may be used to calibrate risk categories in individuals from different countries.

Correlation between allostatic load index and cumulative mortality: a register-based study of Danish municipalities

OBJECTIVES

The aim of this study was to examine population-based allostatic load (AL) indices as an indicator of community health across 14 municipalities in Denmark.

DESIGN

Register-based study.

SETTING

Data derived from: the Lolland-Falster Health Study, the Copenhagen General Population Study and the Danish General Suburban Population Study. Nine biomarkers (systolic blood pressure, diastolic blood pressure, pulse rate, total serum cholesterol, high-density lipoprotein cholesterol, waist-to-hip ratio, triglycerides, C-reactive protein and serum albumin) were divided into high-risk and low-risk values based on clinically accepted criteria, and the AL index was defined as the average between the nine values. All-cause mortality data were obtained from Statistics Denmark.

PARTICIPANTS

We examined a total of 106 808 individuals aged 40-79 years.

PRIMARY OUTCOME MEASURE

Linear regression models were performed to investigate the association between mean AL index and cumulative mortality risk.

RESULTS

Mean AL index was higher in men (range 2.3-3.3) than in women (range 1.7-2.6). We found AL index to be strongly correlated with the cumulative mortality rate, correlation coefficient of 0.82. A unit increase in mean AL index corresponded to an increase in the cumulative mortality rate of 19% (95% CI 13% to 25%) for men, and 16% (95% CI 8% to 23%) for women but this difference was not statistically significant. The overall mean increase in cumulative mortality rate for both men and women was 17% (95% CI 14% to 20%).

CONCLUSIONS

Our findings indicate the population-based AL index to be a strong indicator of community health, and suggest identification of targets for reducing AL.

Impact of individual level uncertainty of lung cancer polygenic risk score (PRS) on risk stratification

BACKGROUND

Although polygenic risk score (PRS) has emerged as a promising tool for predicting cancer risk from genome-wide association studies (GWAS), the individual-level accuracy of lung cancer PRS and the extent to which its impact on subsequent clinical applications remains largely unexplored.

METHODS

Lung cancer PRSs and confidence/credible interval (CI) were constructed using two statistical approaches for each individual: (1) the weighted sum of 16 GWAS-derived significant SNP loci and the CI through the bootstrapping method (PRS-16-CV) and (2) LDpred2 and the CI through posteriors sampling (PRS-Bayes), among 17,166 lung cancer cases and 12,894 controls with European ancestry from the International Lung Cancer Consortium. Individuals were classified into different genetic risk subgroups based on the relationship between their own PRS mean/PRS CI and the population level threshold.

RESULTS

Considerable variances in PRS point estimates at the individual level were observed for both methods, with an average standard deviation (s.d.) of 0.12 for PRS-16-CV and a much larger s.d. of 0.88 for PRS-Bayes. Using PRS-16-CV, only 25.0% of individuals with PRS point estimates in the lowest decile of PRS and 16.8% in the highest decile have their entire 95% CI fully contained in the lowest and highest decile, respectively, while PRS-Bayes was unable to find any eligible individuals. Only 19% of the individuals were concordantly identified as having high genetic risk (> 90th percentile) using the two PRS estimators. An increased relative risk of lung cancer comparing the highest PRS percentile to the lowest was observed when taking the CI into account (OR = 2.73, 95% CI: 2.12-3.50, P-value = 4.13 × 10-15) compared to using PRS-16-CV mean (OR = 2.23, 95% CI: 1.99-2.49, P-value = 5.70 × 10-46). Improved risk prediction performance with higher AUC was consistently observed in individuals identified by PRS-16-CV CI, and the best performance was achieved by incorporating age, gender, and detailed smoking pack-years (AUC: 0.73, 95% CI = 0.72-0.74).

CONCLUSIONS

Lung cancer PRS estimates using different methods have modest correlations at the individual level, highlighting the importance of considering individual-level uncertainty when evaluating the practical utility of PRS.

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.

Evaluation of European-based polygenic risk score for breast cancer in Ashkenazi Jewish women in Israel

BACKGROUND

Polygenic risk score (PRS), calculated based on genome-wide association studies (GWASs), can improve breast cancer (BC) risk assessment. To date, most BC GWASs have been performed in individuals of European (EUR) ancestry, and the generalisation of EUR-based PRS to other populations is a major challenge. In this study, we examined the performance of EUR-based BC PRS models in Ashkenazi Jewish (AJ) women.

METHODS

We generated PRSs based on data on EUR women from the Breast Cancer Association Consortium (BCAC). We tested the performance of the PRSs in a cohort of 2161 AJ women from Israel (1437 cases and 724 controls) from BCAC (BCAC cohort from Israel (BCAC-IL)). In addition, we tested the performance of these EUR-based BC PRSs, as well as the established 313-SNP EUR BC PRS, in an independent cohort of 181 AJ women from Hadassah Medical Center (HMC) in Israel.

RESULTS

In the BCAC-IL cohort, the highest OR per 1 SD was 1.56 (±0.09). The OR for AJ women at the top 10% of the PRS distribution compared with the middle quintile was 2.10 (±0.24). In the HMC cohort, the OR per 1 SD of the EUR-based PRS that performed best in the BCAC-IL cohort was 1.58±0.27. The OR per 1 SD of the commonly used 313-SNP BC PRS was 1.64 (±0.28).

CONCLUSIONS

Extant EUR GWAS data can be used for generating PRSs that identify AJ women with markedly elevated risk of BC and therefore hold promise for improving BC risk assessment in AJ women.

Cell-Free DNA Fragmentomes in the Diagnostic Evaluation of Patients With Symptoms Suggestive of Lung Cancer

BACKGROUND

The diagnostic workup of individuals suspected of having lung cancer can be complex and protracted because conventional symptoms of lung cancer have low specificity and sensitivity.

RESEARCH QUESTION

Among individuals with symptoms of lung cancer, can a blood-based approach to analyze cell-free DNA (cfDNA) fragmentation (the DNA evaluation of fragments for early interception [DELFI] score) enhance evaluation for the possible presence of lung cancer?

STUDY DESIGN AND METHODS

Adults were referred to Bispebjerg Hospital (Copenhagen, Denmark) for diagnostic evaluation of initial imaging anomalies and symptoms consistent with lung cancer. Numbers and types of symptoms were extracted from medical records. cfDNA from plasma samples obtained at the prediagnostic visit was isolated, sequenced, and analyzed for genome-wide cfDNA fragmentation patterns. The relationships among clinical presentation, cancer status, and DELFI score were examined.

RESULTS

A total of 296 individuals were analyzed. Median DELFI scores were higher for those with lung cancer (n = 98) than those without cancer (n = 198; 0.94 vs 0.19; P < .001). In a multivariate model adjusted for age, smoking history, and presenting symptoms, the addition of the DELFI score improved the prediction of lung cancer for those who demonstrated symptoms (area under the receiver operating characteristic curve, 0.74-0.94).

INTERPRETATION

The DELFI score distinguishes individuals with lung cancer from those without cancer better than suspicious symptoms do. These results represent proof-of-concept support that fragmentation-based biomarker approaches may facilitate diagnostic resolution for patients with concerning symptoms of lung cancer.

A genome-wide gene-environment interaction study of breast cancer risk for women of European ancestry

BACKGROUND

Genome-wide studies of gene-environment interactions (G×E) may identify variants associated with disease risk in conjunction with lifestyle/environmental exposures. We conducted a genome-wide G×E analysis of ~ 7.6 million common variants and seven lifestyle/environmental risk factors for breast cancer risk overall and for estrogen receptor positive (ER +) breast cancer.

METHODS

Analyses were conducted using 72,285 breast cancer cases and 80,354 controls of European ancestry from the Breast Cancer Association Consortium. Gene-environment interactions were evaluated using standard unconditional logistic regression models and likelihood ratio tests for breast cancer risk overall and for ER + breast cancer. Bayesian False Discovery Probability was employed to assess the noteworthiness of each SNP-risk factor pairs.

RESULTS

Assuming a 1 × 10-5 prior probability of a true association for each SNP-risk factor pairs and a Bayesian False Discovery Probability < 15%, we identified two independent SNP-risk factor pairs: rs80018847(9p13)-LINGO2 and adult height in association with overall breast cancer risk (ORint = 0.94, 95% CI 0.92-0.96), and rs4770552(13q12)-SPATA13 and age at menarche for ER + breast cancer risk (ORint = 0.91, 95% CI 0.88-0.94).

CONCLUSIONS

Overall, the contribution of G×E interactions to the heritability of breast cancer is very small. At the population level, multiplicative G×E interactions do not make an important contribution to risk prediction in breast cancer.

Association of the CHEK2 c.1100delC variant, radiotherapy, and systemic treatment with contralateral breast cancer risk and breast cancer-specific survival

BACKGROUND

Breast cancer (BC) patients with a germline CHEK2 c.1100delC variant have an increased risk of contralateral BC (CBC) and worse BC-specific survival (BCSS) compared to non-carriers.

AIM

To assessed the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS.

METHODS

Analyses were based on 82,701 women diagnosed with a first primary invasive BC including 963 CHEK2 c.1100delC carriers; median follow-up was 9.1 years. Differential associations with treatment by CHEK2 c.1100delC status were tested by including interaction terms in a multivariable Cox regression model. A multi-state model was used for further insight into the relation between CHEK2 c.1100delC status, treatment, CBC risk and death.

RESULTS

There was no evidence for differential associations of therapy with CBC risk by CHEK2 c.1100delC status. The strongest association with reduced CBC risk was observed for the combination of chemotherapy and endocrine therapy [HR (95% CI): 0.66 (0.55-0.78)]. No association was observed with radiotherapy. Results from the multi-state model showed shorter BCSS for CHEK2 c.1100delC carriers versus non-carriers also after accounting for CBC occurrence [HR (95% CI): 1.30 (1.09-1.56)].

CONCLUSION

Systemic therapy was associated with reduced CBC risk irrespective of CHEK2 c.1100delC status. Moreover, CHEK2 c.1100delC carriers had shorter BCSS, which appears not to be fully explained by their CBC risk.

Mosaic Chromosomal Alterations Are Associated With Increased Lung Cancer Risk: Insight From the INTEGRAL-ILCCO Cohort Analysis

INTRODUCTION

Mosaic chromosomal alterations (mCAs) detected in white blood cells represent a type of clonal hematopoiesis (CH) that is understudied compared with CH-related somatic mutations. A few recent studies indicated their potential link with nonhematological cancers, especially lung cancer.

METHODS

In this study, we investigated the association between mCAs and lung cancer using the high-density genotyping data from the OncoArray study of INTEGRAL-ILCCO, the largest single genetic study of lung cancer with 18,221 lung cancer cases and 14,825 cancer-free controls.

RESULTS

We identified a comprehensive list of autosomal mCAs, ChrX mCAs, and mosaic ChrY (mChrY) losses from these samples. Autosomal mCAs were detected in 4.3% of subjects, in addition to ChrX mCAs in 3.6% of females and mChrY losses in 9.6% of males. Multivariable logistic regression analysis indicated that the presence of autosomal mCAs in white blood cells was associated with an increased lung cancer risk after adjusting for key confounding factors, including age, sex, smoking status, and race. This association was mainly driven by a specific type of mCAs: copy-neutral loss of heterozygosity on autosomal chromosomes. The association between autosome copy-neutral loss of heterozygosity and increased risk of lung cancer was further confirmed in two major histologic subtypes, lung adenocarcinoma and squamous cell carcinoma. In addition, we observed a significant increase of ChrX mCAs and mChrY losses in smokers compared with nonsmokers and racial differences in certain types of mCA events.

CONCLUSIONS

Our study established a link between mCAs in white blood cells and increased risk of lung cancer.

Spectrum and Frequency of Germline FANCM Protein-Truncating Variants in 44,803 European Female Breast Cancer Cases

FANCM germline protein truncating variants (PTVs) are moderate-risk factors for ER-negative breast cancer. We previously described the spectrum of FANCM PTVs in 114 European breast cancer cases. In the present, larger cohort, we report the spectrum and frequency of four common and 62 rare FANCM PTVs found in 274 carriers detected among 44,803 breast cancer cases. We confirmed that p.Gln1701* was the most common PTV in Northern Europe with lower frequencies in Southern Europe. In contrast, p.Gly1906Alafs*12 was the most common PTV in Southern Europe with decreasing frequencies in Central and Northern Europe. We verified that p.Arg658* was prevalent in Central Europe and had highest frequencies in Eastern Europe. We also confirmed that the fourth most common PTV, p.Gln498Thrfs*7, might be a founder variant from Lithuania. Based on the frequency distribution of the carriers of rare PTVs, we showed that the FANCM PTVs spectra in Southwestern and Central Europe were much more heterogeneous than those from Northeastern Europe. These findings will inform the development of more efficient FANCM genetic testing strategies for breast cancer cases from specific European populations.

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.

The PANcreatic Disease ReseArch (PANDoRA) consortium: Ten years' experience of association studies to understand the genetic architecture of pancreatic cancer

Pancreatic cancer has an incidence that almost matches its mortality. Only a small number of risk factors and 33 susceptibility loci have been identified. so Moreover, the relative rarity of pancreatic cancer poses significant hurdles for research aimed at increasing our knowledge of the genetic mechanisms contributing to the disease. Additionally, the inability to adequately power research questions prevents small monocentric studies from being successful. Several consortia have been established to pursue a better understanding of the genetic architecture of pancreatic cancers. The Pancreatic disease research (PANDoRA) consortium is the largest in Europe. PANDoRA is spread across 12 European countries, Brazil and Japan, bringing together 29 basic and clinical research groups. In the last ten years, PANDoRA has contributed to the discovery of 25 susceptibility loci, a feat that will be instrumental in stratifying the population by risk and optimizing preventive strategies.

Genome-wide association study of lung adenocarcinoma in East Asia and comparison with a European population

Lung adenocarcinoma is the most common type of lung cancer. Known risk variants explain only a small fraction of lung adenocarcinoma heritability. Here, we conducted a two-stage genome-wide association study of lung adenocarcinoma of East Asian ancestry (21,658 cases and 150,676 controls; 54.5% never-smokers) and identified 12 novel susceptibility variants, bringing the total number to 28 at 25 independent loci. Transcriptome-wide association analyses together with colocalization studies using a Taiwanese lung expression quantitative trait loci dataset (n = 115) identified novel candidate genes, including FADS1 at 11q12 and ELF5 at 11p13. In a multi-ancestry meta-analysis of East Asian and European studies, four loci were identified at 2p11, 4q32, 16q23, and 18q12. At the same time, most of our findings in East Asian populations showed no evidence of association in European populations. In our studies drawn from East Asian populations, a polygenic risk score based on the 25 loci had a stronger association in never-smokers vs. individuals with a history of smoking (Pinteraction = 0.0058). These findings provide new insights into the etiology of lung adenocarcinoma in individuals from East Asian populations, which could be important in developing translational applications.

PREDICT validity for prognosis of breast cancer patients with pathogenic BRCA1/2 variants

We assessed the PREDICT v 2.2 for prognosis of breast cancer patients with pathogenic germline BRCA1 and BRCA2 variants, using follow-up data from 5453 BRCA1/2 carriers from the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) and the Breast Cancer Association Consortium (BCAC). PREDICT for estrogen receptor (ER)-negative breast cancer had modest discrimination for BRCA1 carrier patients overall (Gönen & Heller unbiased concordance 0.65 in CIMBA, 0.64 in BCAC), but it distinguished clearly the high-mortality group from lower risk categories. In an analysis of low to high risk categories by PREDICT score percentiles, the observed mortality was consistently lower than the expected mortality, but the confidence intervals always included the calibration slope. Altogether, our results encourage the use of the PREDICT ER-negative model in management of breast cancer patients with germline BRCA1 variants. For the PREDICT ER-positive model, the discrimination was slightly lower in BRCA2 variant carriers (concordance 0.60 in CIMBA, 0.65 in BCAC). Especially, inclusion of the tumor grade distorted the prognostic estimates. The breast cancer mortality of BRCA2 carriers was underestimated at the low end of the PREDICT score distribution, whereas at the high end, the mortality was overestimated. These data suggest that BRCA2 status should also be taken into consideration with tumor characteristics, when estimating the prognosis of ER-positive breast cancer patients.

Associations of a Breast Cancer Polygenic Risk Score With Tumor Characteristics and Survival

PURPOSE

A polygenic risk score (PRS) consisting of 313 common genetic variants (PRS313) is associated with risk of breast cancer and contralateral breast cancer. This study aimed to evaluate the association of the PRS313 with clinicopathologic characteristics of, and survival following, breast cancer.

METHODS

Women with invasive breast cancer were included, 98,397 of European ancestry and 12,920 of Asian ancestry, from the Breast Cancer Association Consortium (BCAC), and 683 women from the European MINDACT trial. Associations between PRS313 and clinicopathologic characteristics, including the 70-gene signature for MINDACT, were evaluated using logistic regression analyses. Associations of PRS313 (continuous, per standard deviation) with overall survival (OS) and breast cancer-specific survival (BCSS) were evaluated with Cox regression, adjusted for clinicopathologic characteristics and treatment.

RESULTS

The PRS313 was associated with more favorable tumor characteristics. In BCAC, increasing PRS313 was associated with lower grade, hormone receptor-positive status, and smaller tumor size. In MINDACT, PRS313 was associated with a low risk 70-gene signature. In European women from BCAC, higher PRS313 was associated with better OS and BCSS: hazard ratio (HR) 0.96 (95% CI, 0.94 to 0.97) and 0.96 (95% CI, 0.94 to 0.98), but the association disappeared after adjustment for clinicopathologic characteristics (and treatment): OS HR, 1.01 (95% CI, 0.98 to 1.05) and BCSS HR, 1.02 (95% CI, 0.98 to 1.07). The results in MINDACT and Asian women from BCAC were consistent.

CONCLUSION

An increased PRS313 is associated with favorable tumor characteristics, but is not independently associated with prognosis. Thus, PRS313 has no role in the clinical management of primary breast cancer at the time of diagnosis. Nevertheless, breast cancer mortality rates will be higher for women with higher PRS313 as increasing PRS313 is associated with an increased risk of disease. This information is crucial for modeling effective stratified screening programs.

Incorporating Alternative Polygenic Risk Scores into the BOADICEA Breast Cancer Risk Prediction Model

BACKGROUND

The multifactorial risk prediction model BOADICEA enables identification of women at higher or lower risk of developing breast cancer. BOADICEA models genetic susceptibility in terms of the effects of rare variants in breast cancer susceptibility genes and a polygenic component, decomposed into an unmeasured and a measured component - the polygenic risk score (PRS). The current version was developed using a 313 SNP PRS. Here, we evaluated approaches to incorporating this PRS and alternative PRS in BOADICEA.

METHODS

The mean, SD, and proportion of the overall polygenic component explained by the PRS (α2) need to be estimated. $\alpha $ was estimated using logistic regression, where the age-specific log-OR is constrained to be a function of the age-dependent polygenic relative risk in BOADICEA; and using a retrospective likelihood (RL) approach that models, in addition, the unmeasured polygenic component.

RESULTS

Parameters were computed for 11 PRS, including 6 variations of the 313 SNP PRS used in clinical trials and implementation studies. The logistic regression approach underestimates $\alpha $, as compared with the RL estimates. The RL $\alpha $ estimates were very close to those obtained by assuming proportionality to the OR per 1 SD, with the constant of proportionality estimated using the 313 SNP PRS. Small variations in the SNPs included in the PRS can lead to large differences in the mean.

CONCLUSIONS

BOADICEA can be readily adapted to different PRS in a manner that maintains consistency of the model.

IMPACT

: The methods described facilitate comprehensive breast cancer risk assessment.

The impact of coding germline variants on contralateral breast cancer risk and survival

Evidence linking coding germline variants in breast cancer (BC)-susceptibility genes other than BRCA1, BRCA2, and CHEK2 with contralateral breast cancer (CBC) risk and breast cancer-specific survival (BCSS) is scarce. The aim of this study was to assess the association of protein-truncating variants (PTVs) and rare missense variants (MSVs) in nine known (ATM, BARD1, BRCA1, BRCA2, CHEK2, PALB2, RAD51C, RAD51D, and TP53) and 25 suspected BC-susceptibility genes with CBC risk and BCSS. Hazard ratios (HRs) and 95% confidence intervals (CIs) were estimated with Cox regression models. Analyses included 34,401 women of European ancestry diagnosed with BC, including 676 CBCs and 3,449 BC deaths; the median follow-up was 10.9 years. Subtype analyses were based on estrogen receptor (ER) status of the first BC. Combined PTVs and pathogenic/likely pathogenic MSVs in BRCA1, BRCA2, and TP53 and PTVs in CHEK2 and PALB2 were associated with increased CBC risk [HRs (95% CIs): 2.88 (1.70-4.87), 2.31 (1.39-3.85), 8.29 (2.53-27.21), 2.25 (1.55-3.27), and 2.67 (1.33-5.35), respectively]. The strongest evidence of association with BCSS was for PTVs and pathogenic/likely pathogenic MSVs in BRCA2 (ER-positive BC) and TP53 and PTVs in CHEK2 [HRs (95% CIs): 1.53 (1.13-2.07), 2.08 (0.95-4.57), and 1.39 (1.13-1.72), respectively, after adjusting for tumor characteristics and treatment]. HRs were essentially unchanged when censoring for CBC, suggesting that these associations are not completely explained by increased CBC risk, tumor characteristics, or treatment. There was limited evidence of associations of PTVs and/or rare MSVs with CBC risk or BCSS for the 25 suspected BC genes. The CBC findings are relevant to treatment decisions, follow-up, and screening after BC diagnosis.

Diminishing benefits of urban living for children and adolescents' growth and development

Optimal growth and development in childhood and adolescence is crucial for lifelong health and well-being1-6. Here we used data from 2,325 population-based studies, with measurements of height and weight from 71 million participants, to report the height and body-mass index (BMI) of children and adolescents aged 5-19 years on the basis of rural and urban place of residence in 200 countries and territories from 1990 to 2020. In 1990, children and adolescents residing in cities were taller than their rural counterparts in all but a few high-income countries. By 2020, the urban height advantage became smaller in most countries, and in many high-income western countries it reversed into a small urban-based disadvantage. The exception was for boys in most countries in sub-Saharan Africa and in some countries in Oceania, south Asia and the region of central Asia, Middle East and north Africa. In these countries, successive cohorts of boys from rural places either did not gain height or possibly became shorter, and hence fell further behind their urban peers. The difference between the age-standardized mean BMI of children in urban and rural areas was <1.1 kg m-2 in the vast majority of countries. Within this small range, BMI increased slightly more in cities than in rural areas, except in south Asia, sub-Saharan Africa and some countries in central and eastern Europe. Our results show that in much of the world, the growth and developmental advantages of living in cities have diminished in the twenty-first century, whereas in much of sub-Saharan Africa they have amplified.

Association of the CHEK2 c.1100delC variant, radiotherapy, and systemic treatment with contralateral breast cancer risk and breast cancer-specific survival

Breast cancer (BC) patients with a germline CHEK2 c.1100delC variant have an increased risk of contralateral BC (CBC) and worse BC-specific survival (BCSS) compared to non-carriers. We aimed to assess the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS. Analyses were based on 82,701 women diagnosed with invasive BC including 963 CHEK2 c.1100delC carriers; median follow-up was 9.1 years. Differential associations of treatment by CHEK2 c.1100delC status were tested by including interaction terms in a multivariable Cox regression model. A multi-state model was used for further insight into the relation between CHEK2 c.1100delC status, treatment, CBC risk and death. There was no evidence for differential associations of therapy with CBC risk by CHEK2 c.1100delC status The strongest association with reduced CBC risk was observed for the combination of chemotherapy and endocrine therapy [HR(95%CI): 0.66 (0.55-0.78)]. No association was observed with radiotherapy. Results from the multi-state model showed shorter BCSS for CHEK2 c.1100delC carriers versus non-carriers also after accounting for CBC occurrence [HR(95%CI) :1.30 (1.09-1.56)]. In conclusion, systemic therapy was associated with reduced CBC risk irrespective of CHEK2 c.1100delC status. Moreover, CHEK2 c.1100delC carriers had shorter BCSS, which appears not to be fully explained by their CBC risk. (Main MS: 3201 words).

Prognostic value of circulating proteins in patients undergoing surgery for pancreatic cancer

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) is the fourth leading cause of cancer death. Less than 20% of patients are diagnosed with resectable disease. Identifying truly resectable disease is challenging because 20%-40% of the patients subjected to resection are found to have advanced disease during surgery. The aim of our study was to identify panels of circulating proteins that could be used to distinguish patients with unresectable PDAC from patients with resectable PDAC and to identify prognostic signatures for both groups.

METHODS

We measured 92 circulating immuno-oncology-related proteins using the proximity extension assay from Olink Proteomics in 273 patients eligible for surgery for PDAC. Two bioinformaticians worked independently of one another on the same data. LASSO and Ridge regression were used in the statistical analyses.

RESULTS

One protein index for determining resectability had an AUC value of 0.66. Several indices for prognosis had AUC values between 0.50 and 0.75 and were therefore not better than existing prognostic markers.

DISCUSSION

Our study did not reveal any new high-performing protein panels that could be used to identify patients with inoperable PDAC before surgery. The panel of 92 proteins investigated has previously been found to be applicable for diagnostic use in patients with PDAC, but it does not seem to warrant further investigation regarding resectability in the subgroup of patients with PDAC referred to surgery.

Aggregation tests identify new gene associations with breast cancer in populations with diverse ancestry

BACKGROUND

Low-frequency variants play an important role in breast cancer (BC) susceptibility. Gene-based methods can increase power by combining multiple variants in the same gene and help identify target genes.

METHODS

We evaluated the potential of gene-based aggregation in the Breast Cancer Association Consortium cohorts including 83,471 cases and 59,199 controls. Low-frequency variants were aggregated for individual genes' coding and regulatory regions. Association results in European ancestry samples were compared to single-marker association results in the same cohort. Gene-based associations were also combined in meta-analysis across individuals with European, Asian, African, and Latin American and Hispanic ancestry.

RESULTS

In European ancestry samples, 14 genes were significantly associated (q < 0.05) with BC. Of those, two genes, FMNL3 (P = 6.11 × 10-6) and AC058822.1 (P = 1.47 × 10-4), represent new associations. High FMNL3 expression has previously been linked to poor prognosis in several other cancers. Meta-analysis of samples with diverse ancestry discovered further associations including established candidate genes ESR1 and CBLB. Furthermore, literature review and database query found further support for a biologically plausible link with cancer for genes CBLB, FMNL3, FGFR2, LSP1, MAP3K1, and SRGAP2C.

CONCLUSIONS

Using extended gene-based aggregation tests including coding and regulatory variation, we report identification of plausible target genes for previously identified single-marker associations with BC as well as the discovery of novel genes implicated in BC development. Including multi ancestral cohorts in this study enabled the identification of otherwise missed disease associations as ESR1 (P = 1.31 × 10-5), demonstrating the importance of diversifying study cohorts.

Preoperative plasma vitamin D in patients with localized colorectal cancer: Age-dependent association with inflammation, postoperative complications, and survival

INTRODUCTION

Aging is often associated with low-grade chronic inflammation and a senescent immune system. Vitamin D is a regulator of immune function, and low plasma vitamin D is associated with poor health. The association between plasma vitamin D and inflammatory biomarkers and risk of postoperative complications and survival in patients with colorectal cancer (CRC) is unknown. Our aim was to investigate these associations and how they are influenced by age.

MATERIALS AND METHODS

Circulating vitamin D and the inflammatory biomarkers C-reactive protein (CRP), interleukin (IL)-6, and YKL-40 were measured in 398 patients with stage I-III CRC preoperatively. Older patients (≥70 years, n = 208) were compared to younger patients (<70 years, n = 190). The relation between vitamin D and complications and high inflammatory biomarker levels was presented by odds ratios ([OR], 95% confidence interval [CI]). Associations with survival were presented with hazard ratios ([HR], 95% CI).

RESULTS

Plasma vitamin D was higher in older patients than in younger patients (75 vs. 67 nmol/L, P = 0.001). High vitamin D was associated with low plasma CRP in younger patients (OR = 0.35, 95% CI 0.17-0.76), but not in older patients (OR = 0.93, 0.49-1.76). High vitamin D in older patients with CRC was associated with reduced risk of major complications (OR = 0.52, 0.28-0.95). This was not found in younger patients (OR = 1.47, 0.70-3.11). Deficient vitamin D (<25 nmol/L) was associated with short overall survival compared to sufficient (>50 nmol/L) irrespective of age (HR = 3.39, 1.27-9.37, P = 0.02).

CONCLUSION

For patients with localized CRC, high vitamin D levels before resection were associated with reduced risk of high inflammatory biomarkers for younger patients and reduced risk of major postoperative complications for older patients. Vitamin D deficiency was associated with reduced survival regardless of age.

Kidney function and risk of dementia: Observational study, meta-analysis, and two-sample mendelian randomization study

Whether impaired kidney function is associated with increased risk of developing dementia is unclear. We investigated the association between estimated glomerular filtration rate (eGFR) and dementia. Using a triangulation approach, we performed (1) a prospective study in 90,369 Danes from the Copenhagen General Population Study (CGPS), (2) a meta-analysis in 468,699 Scandinavians (including CGPS) and (3) a two-sample Mendelian randomization study in 218,792-1,004,040 Europeans using summary data from largest publicly available genome wide association studies (GWASs). During up to 15 years of follow-up (CGPS), 2,468 individuals developed dementia. Age and sex standardized percentile of eGFR below versus above the median conferred a multifactorially adjusted hazard ratio of 1.09 (95% confidence interval: 1.01-1.18). In meta-analysis, random-effects risk of dementia was 1.14 (1.06-1.22) for mildly decreased eGFR (60-90 mL/min/1.73 m²), 1.31 (0.92-1.87) for moderately decreased eGFR (30-59 mL/min/1.73 m²) and 1.91 (1.21-3.01) for severely decreased eGFR (< 30 mL/min/1.73 m²), compared to reference eGFR (> 90 mL/min/1.73 m²). Using directly comparable eGFR measures (log[eGFR] scaled to one standard deviation, as well as eGFR below versus above 60 mL/min/1.73 m²), we found no association with risk of dementia in observational CGPS or in Mendelian randomization analyses. In conclusion, impaired kidney function was associated with modestly increased risk of developing dementia. This was not supported by causal, genetic analyses using a Mendelian randomization approach. However, future stronger genetic instruments for kidney function and larger GWASs with more dementia cases, particularly for the vascular dementia subtype, warrant a re-evaluation of the causal hypothesis.

Combined plasma C‐reactive protein, interleukin 6 and YKL ‐40 for detection of cancer and prognosis in patients with serious nonspecific symptoms and signs of cancer

BACKGROUND AND METHODS

Inflammation is a hallmark of cancer and its progression. Plasma levels of C-reactive protein (CRP), interleukin-6 (IL-6) and YKL-40 reflect inflammation, and are elevated in patients with cancer. This study investigated whether plasma CRP, IL-6 and YKL-40 had diagnostic value in 753 patients referred with nonspecific signs and symptoms of cancer to a diagnostic outpatient clinic.

RESULTS

In total, 111 patients were diagnosed with cancer within 3 months and 30 after 3 months. CRP, IL-6 and YKL-40 were elevated in 44%, 60% and 45% of the cancer patients, and in 15%, 33% and 25% of the patients without cancer. Elevated levels of all three markers were associated with risk of cancer within 3 months: CRP (odds ratio (OR) 4.41, 95% confidence interval (CI) 2.86-6.81), IL-6 (OR = 2.89, 1.91-4.37) and YKL-40 (OR = 2.42, 1.59-3.66). Multivariate explorative analyses showed that increasing values were associated with the risk of getting a cancer diagnosis (continuous scale: CRP (OR = 1.28, 1.12-1.47), carcinoembryonic antigen (CEA) (OR = 1.61, 1.41-1.98), CA19-9 (OR = 1.15, 1.03-1.29), age (OR = 1.29, 1.02-1.63); dichotomized values: CRP (OR = 2.54, 1.39-4.66), CEA (OR = 4.22, 2.13-8.34), age (OR = 1.42, 1.13-1.80)). CRP had the highest diagnostic value (area under the curve = 0.69). Combined high CRP, IL-6 and YKL-40 was associated with short overall survival (HR = 3.8, 95% CI 2.5-5.9, p < 0.001).

CONCLUSION

In conclusion, plasma CRP, IL-6 and YKL-40 alone or combined cannot be used to identify patients with cancer, but high levels were associated with poor prognosis. CRP may be useful to indicate whether further diagnostic evaluation is needed when patients present with nonspecific signs and symptoms of cancer.

Clonal haematopoiesis of indeterminate potential and impaired kidney function-A Danish general population study with 11 years follow-up

The myeloproliferative neoplasms are associated with chronic kidney disease but whether clonal haematopoiesis of indeterminate potential (CHIP) is associated with impaired kidney function is unknown. In the Danish General Suburban Population Study (N = 19 958) from 2010 to 2013, 645 individuals were positive for JAK2V617F (N = 613) or CALR (N = 32) mutations. Mutation-positive individuals without haematological malignancy were defined as having CHIP (N = 629). We used multiple and inverse probability weighted (IPW)-adjusted linear regression analysis to estimate adjusted mean (95% confidence interval) differences in estimated glomerular filtration rate (eGFR; ml/min/1.73 m² ) by mutation status, variant allele frequency (VAF%), blood cell counts, and neutrophil-to-lymphocyte ratio (NLR). We performed 11-year longitudinal follow-up of eGFR in all individuals. Compared to CHIP-negative individuals, the mean differences in eGFR were -5.6 (-10.3, -0.8, p = .02) for CALR, -11.9 (-21.4, -2.4, p = 0.01) for CALR type 2, and -10.1 (-18.1, -2.2, p = .01) for CALR with VAF ≥ 1%. The IPW-adjusted linear regression analyses showed similar results. NLR was negatively associated with eGFR. Individuals with CALR type 2 had a worse 11-year longitudinal follow-up on eGFR compared to CHIP-negative individuals (p = .004). In conclusion, individuals with CALR mutations, especially CALR type 2, had impaired kidney function compared to CHIP-negative individuals as measured by a lower eGFR at baseline and during 11-year follow-up.

PredictCBC-2.0: a contralateral breast cancer risk prediction model developed and validated in ~ 200,000 patients

BACKGROUND

Prediction of contralateral breast cancer (CBC) risk is challenging due to moderate performances of the known risk factors. We aimed to improve our previous risk prediction model (PredictCBC) by updated follow-up and including additional risk factors.

METHODS

We included data from 207,510 invasive breast cancer patients participating in 23 studies. In total, 8225 CBC events occurred over a median follow-up of 10.2 years. In addition to the previously included risk factors, PredictCBC-2.0 included CHEK2 c.1100delC, a 313 variant polygenic risk score (PRS-313), body mass index (BMI), and parity. Fine and Gray regression was used to fit the model. Calibration and a time-dependent area under the curve (AUC) at 5 and 10 years were assessed to determine the performance of the models. Decision curve analysis was performed to evaluate the net benefit of PredictCBC-2.0 and previous PredictCBC models.

RESULTS

The discrimination of PredictCBC-2.0 at 10 years was higher than PredictCBC with an AUC of 0.65 (95% prediction intervals (PI) 0.56-0.74) versus 0.63 (95%PI 0.54-0.71). PredictCBC-2.0 was well calibrated with an observed/expected ratio at 10 years of 0.92 (95%PI 0.34-2.54). Decision curve analysis for contralateral preventive mastectomy (CPM) showed the potential clinical utility of PredictCBC-2.0 between thresholds of 4 and 12% 10-year CBC risk for BRCA1/2 mutation carriers and non-carriers.

CONCLUSIONS

Additional genetic information beyond BRCA1/2 germline mutations improved CBC risk prediction and might help tailor clinical decision-making toward CPM or alternative preventive strategies. Identifying patients who benefit from CPM, especially in the general breast cancer population, remains challenging.

Physical activity, sedentary time and breast cancer risk: a Mendelian randomisation study

OBJECTIVES

Physical inactivity and sedentary behaviour are associated with higher breast cancer risk in observational studies, but ascribing causality is difficult. Mendelian randomisation (MR) assesses causality by simulating randomised trial groups using genotype. We assessed whether lifelong physical activity or sedentary time, assessed using genotype, may be causally associated with breast cancer risk overall, pre/post-menopause, and by case-groups defined by tumour characteristics.

METHODS

We performed two-sample inverse-variance-weighted MR using individual-level Breast Cancer Association Consortium case-control data from 130 957 European-ancestry women (69 838 invasive cases), and published UK Biobank data (n=91 105-377 234). Genetic instruments were single nucleotide polymorphisms (SNPs) associated in UK Biobank with wrist-worn accelerometer-measured overall physical activity (nsnps=5) or sedentary time (nsnps=6), or accelerometer-measured (nsnps=1) or self-reported (nsnps=5) vigorous physical activity.

RESULTS

Greater genetically-predicted overall activity was associated with lower breast cancer overall risk (OR=0.59; 95% confidence interval (CI) 0.42 to 0.83 per-standard deviation (SD;~8 milligravities acceleration)) and for most case-groups. Genetically-predicted vigorous activity was associated with lower risk of pre/perimenopausal breast cancer (OR=0.62; 95% CI 0.45 to 0.87,≥3 vs. 0 self-reported days/week), with consistent estimates for most case-groups. Greater genetically-predicted sedentary time was associated with higher hormone-receptor-negative tumour risk (OR=1.77; 95% CI 1.07 to 2.92 per-SD (~7% time spent sedentary)), with elevated estimates for most case-groups. Results were robust to sensitivity analyses examining pleiotropy (including weighted-median-MR, MR-Egger).

CONCLUSION

Our study provides strong evidence that greater overall physical activity, greater vigorous activity, and lower sedentary time are likely to reduce breast cancer risk. More widespread adoption of active lifestyles may reduce the burden from the most common cancer in women.

Incorporating progesterone receptor expression into the PREDICT breast prognostic model

BACKGROUND

Predict Breast (www.predict.nhs.uk) is an online prognostication and treatment benefit tool for early invasive breast cancer. The aim of this study was to incorporate the prognostic effect of progesterone receptor (PR) status into a new version of PREDICT and to compare its performance to the current version (2.2).

METHOD

The prognostic effect of PR status was based on the analysis of data from 45,088 European patients with breast cancer from 49 studies in the Breast Cancer Association Consortium. Cox proportional hazard models were used to estimate the hazard ratio for PR status. Data from a New Zealand study of 11,365 patients with early invasive breast cancer were used for external validation. Model calibration and discrimination were used to test the model performance.

RESULTS

Having a PR-positive tumour was associated with a 23% and 28% lower risk of dying from breast cancer for women with oestrogen receptor (ER)-negative and ER-positive breast cancer, respectively. The area under the ROC curve increased with the addition of PR status from 0.807 to 0.809 for patients with ER-negative tumours (p = 0.023) and from 0.898 to 0.902 for patients with ER-positive tumours (p = 2.3 × 10-6) in the New Zealand cohort. Model calibration was modest with 940 observed deaths compared to 1151 predicted.

CONCLUSION

The inclusion of the prognostic effect of PR status to PREDICT Breast has led to an improvement of model performance and more accurate absolute treatment benefit predictions for individual patients. Further studies should determine whether the baseline hazard function requires recalibration.

Genome-wide interaction analysis identified low-frequency variants with sex disparity in lung cancer risk

Differences by sex in lung cancer incidence and mortality have been reported which cannot be fully explained by sex differences in smoking behavior, implying existence of genetic and molecular basis for sex disparity in lung cancer development. However, the information about sex dimorphism in lung cancer risk is quite limited despite the great success in lung cancer association studies. By adopting a stringent two-stage analysis strategy, we performed a genome-wide gene-sex interaction analysis using genotypes from a lung cancer cohort including ~ 47 000 individuals with European ancestry. Three low-frequency variants (minor allele frequency < 0.05), rs17662871 [odds ratio (OR) = 0.71, P = 4.29×10-8); rs79942605 (OR = 2.17, P = 2.81×10-8) and rs208908 (OR = 0.70, P = 4.54×10-8) were identified with different risk effect of lung cancer between men and women. Further expression quantitative trait loci and functional annotation analysis suggested rs208908 affects lung cancer risk through differential regulation of Coxsackie virus and adenovirus receptor gene expression in lung tissues between men and women. Our study is one of the first studies to provide novel insights about the genetic and molecular basis for sex disparity in lung cancer development.

A Large-Scale Genome-Wide Gene-Gene Interaction Study of Lung Cancer Susceptibility in Europeans With a Trans-Ethnic Validation in Asians

INTRODUCTION

Although genome-wide association studies have been conducted to investigate genetic variation of lung tumorigenesis, little is known about gene-gene (G × G) interactions that may influence the risk of non-small cell lung cancer (NSCLC).

METHODS

Leveraging a total of 445,221 European-descent participants from the International Lung Cancer Consortium OncoArray project, Transdisciplinary Research in Cancer of the Lung and UK Biobank, we performed a large-scale genome-wide G × G interaction study on European NSCLC risk by a series of analyses. First, we used BiForce to evaluate and rank more than 58 billion G × G interactions from 340,958 single-nucleotide polymorphisms (SNPs). Then, the top interactions were further tested by demographically adjusted logistic regression models. Finally, we used the selected interactions to build lung cancer screening models of NSCLC, separately, for never and ever smokers.

RESULTS

With the Bonferroni correction, we identified eight statistically significant pairs of SNPs, which predominantly appeared in the 6p21.32 and 5p15.33 regions (e.g., rs521828C6orf10 and rs204999PRRT1, ORinteraction = 1.17, p = 6.57 × 10-13; rs3135369BTNL2 and rs2858859HLA-DQA1, ORinteraction = 1.17, p = 2.43 × 10-13; rs2858859HLA-DQA1 and rs9275572HLA-DQA2, ORinteraction = 1.15, p = 2.84 × 10-13; rs2853668TERT and rs62329694CLPTM1L, ORinteraction = 0.73, p = 2.70 × 10-13). Notably, even with much genetic heterogeneity across ethnicities, three pairs of SNPs in the 6p21.32 region identified from the European-ancestry population remained significant among an Asian population from the Nanjing Medical University Global Screening Array project (rs521828C6orf10 and rs204999PRRT1, ORinteraction = 1.13, p = 0.008; rs3135369BTNL2 and rs2858859HLA-DQA1, ORinteraction = 1.11, p = 5.23 × 10-4; rs3135369BTNL2 and rs9271300HLA-DQA1, ORinteraction = 0.89, p = 0.006). The interaction-empowered polygenetic risk score that integrated classical polygenetic risk score and G × G information score was remarkable in lung cancer risk stratification.

CONCLUSIONS

Important G × G interactions were identified and enriched in the 5p15.33 and 6p21.32 regions, which may enhance lung cancer screening models.

Breast cancer risks associated with missense variants in breast cancer susceptibility genes

BACKGROUND

Protein truncating variants in ATM, BRCA1, BRCA2, CHEK2, and PALB2 are associated with increased breast cancer risk, but risks associated with missense variants in these genes are uncertain.

METHODS

We analyzed data on 59,639 breast cancer cases and 53,165 controls from studies participating in the Breast Cancer Association Consortium BRIDGES project. We sampled training (80%) and validation (20%) sets to analyze rare missense variants in ATM (1146 training variants), BRCA1 (644), BRCA2 (1425), CHEK2 (325), and PALB2 (472). We evaluated breast cancer risks according to five in silico prediction-of-deleteriousness algorithms, functional protein domain, and frequency, using logistic regression models and also mixture models in which a subset of variants was assumed to be risk-associated.

RESULTS

The most predictive in silico algorithms were Helix (BRCA1, BRCA2 and CHEK2) and CADD (ATM). Increased risks appeared restricted to functional protein domains for ATM (FAT and PIK domains) and BRCA1 (RING and BRCT domains). For ATM, BRCA1, and BRCA2, data were compatible with small subsets (approximately 7%, 2%, and 0.6%, respectively) of rare missense variants giving similar risk to those of protein truncating variants in the same gene. For CHEK2, data were more consistent with a large fraction (approximately 60%) of rare missense variants giving a lower risk (OR 1.75, 95% CI (1.47-2.08)) than CHEK2 protein truncating variants. There was little evidence for an association with risk for missense variants in PALB2. The best fitting models were well calibrated in the validation set.

CONCLUSIONS

These results will inform risk prediction models and the selection of candidate variants for functional assays and could contribute to the clinical reporting of gene panel testing for breast cancer susceptibility.

Pathology of Tumors Associated With Pathogenic Germline Variants in 9 Breast Cancer Susceptibility Genes

IMPORTANCE

Rare germline genetic variants in several genes are associated with increased breast cancer (BC) risk, but their precise contributions to different disease subtypes are unclear. This information is relevant to guidelines for gene panel testing and risk prediction.

OBJECTIVE

To characterize tumors associated with BC susceptibility genes in large-scale population- or hospital-based studies.

DESIGN, SETTING, AND PARTICIPANTS

The multicenter, international case-control analysis of the BRIDGES study included 42 680 patients and 46 387 control participants, comprising women aged 18 to 79 years who were sampled independently of family history from 38 studies. Studies were conducted between 1991 and 2016. Sequencing and analysis took place between 2016 and 2021.

EXPOSURES

Protein-truncating variants and likely pathogenic missense variants in ATM, BARD1, BRCA1, BRCA2, CHEK2, PALB2, RAD51C, RAD51D, and TP53.

MAIN OUTCOMES AND MEASURES

The intrinsic-like BC subtypes as defined by estrogen receptor, progesterone receptor, and ERBB2 (formerly known as HER2) status, and tumor grade; morphology; size; stage; lymph node involvement; subtype-specific odds ratios (ORs) for carrying protein-truncating variants and pathogenic missense variants in the 9 BC susceptibility genes.

RESULTS

The mean (SD) ages at interview (control participants) and diagnosis (cases) were 55.1 (11.9) and 55.8 (10.6) years, respectively; all participants were of European or East Asian ethnicity. There was substantial heterogeneity in the distribution of intrinsic subtypes by gene. RAD51C, RAD51D, and BARD1 variants were associated mainly with triple-negative disease (OR, 6.19 [95% CI, 3.17-12.12]; OR, 6.19 [95% CI, 2.99-12.79]; and OR, 10.05 [95% CI, 5.27-19.19], respectively). CHEK2 variants were associated with all subtypes (with ORs ranging from 2.21-3.17) except for triple-negative disease. For ATM variants, the association was strongest for the hormone receptor (HR)+ERBB2- high-grade subtype (OR, 4.99; 95% CI, 3.68-6.76). BRCA1 was associated with increased risk of all subtypes, but the ORs varied widely, being highest for triple-negative disease (OR, 55.32; 95% CI, 40.51-75.55). BRCA2 and PALB2 variants were also associated with triple-negative disease. TP53 variants were most strongly associated with HR+ERBB2+ and HR-ERBB2+ subtypes. Tumors occurring in pathogenic variant carriers were of higher grade. For most genes and subtypes, a decline in ORs was observed with increasing age. Together, the 9 genes were associated with 27.3% of all triple-negative tumors in women 40 years or younger.

CONCLUSIONS AND RELEVANCE

The results of this case-control study suggest that variants in the 9 BC risk genes differ substantially in their associated pathology but are generally associated with triple-negative and/or high-grade disease. Knowing the age and tumor subtype distributions associated with individual BC genes can potentially aid guidelines for gene panel testing, risk prediction, and variant classification and guide targeted screening strategies.

Epigenetic Regulation of F2RL3 Associates With Myocardial Infarction and Platelet Function

BACKGROUND

DNA hypomethylation at the F2RL3 (F2R like thrombin or trypsin receptor 3) locus has been associated with both smoking and atherosclerotic cardiovascular disease; whether these smoking-related associations form a pathway to disease is unknown. F2RL3 encodes protease-activated receptor 4, a potent thrombin receptor expressed on platelets. Given the role of thrombin in platelet activation and the role of thrombus formation in myocardial infarction, alterations to this biological pathway could be important for ischemic cardiovascular disease.

METHODS

We conducted multiple independent experiments to assess whether DNA hypomethylation at F2RL3 in response to smoking is associated with risk of myocardial infarction via changes to platelet reactivity. Using cohort data (N=3205), we explored the relationship between smoking, DNA hypomethylation at F2RL3, and myocardial infarction. We compared platelet reactivity in individuals with low versus high DNA methylation at F2RL3 (N=41). We used an in vitro model to explore the biological response of F2RL3 to cigarette smoke extract. Finally, a series of reporter constructs were used to investigate how differential methylation could impact F2RL3 gene expression.

RESULTS

Observationally, DNA methylation at F2RL3 mediated an estimated 34% of the smoking effect on increased risk of myocardial infarction. An association between methylation group (low/high) and platelet reactivity was observed in response to PAR4 (protease-activated receptor 4) stimulation. In cells, cigarette smoke extract exposure was associated with a 4.9% to 9.3% reduction in DNA methylation at F2RL3 and a corresponding 1.7-(95% CI, 1.2-2.4, P=0.04) fold increase in F2RL3 mRNA. Results from reporter assays suggest the exon 2 region of F2RL3 may help control gene expression.

CONCLUSIONS

Smoking-induced epigenetic DNA hypomethylation at F2RL3 appears to increase PAR4 expression with potential downstream consequences for platelet reactivity. Combined evidence here not only identifies F2RL3 DNA methylation as a possible contributory pathway from smoking to cardiovascular disease risk but from any feature potentially influencing F2RL3 regulation in a similar manner.

Gene-gene interaction of AhRwith and within the Wntcascade affects susceptibility to lung cancer

BACKGROUND

Aberrant Wnt signalling, regulating cell development and stemness, influences the development of many cancer types. The Aryl hydrocarbon receptor (AhR) mediates tumorigenesis of environmental pollutants. Complex interaction patterns of genes assigned to AhR/Wnt-signalling were recently associated with lung cancer susceptibility.

AIM

To assess the association and predictive ability of AhR/Wnt-genes with lung cancer in cases and controls of European descent.

METHODS

Odds ratios (OR) were estimated for genomic variants assigned to the Wnt agonist and the antagonistic genes DKK2, DKK3, DKK4, FRZB, SFRP4 and Axin2. Logistic regression models with variable selection were trained, validated and tested to predict lung cancer, at which other previously identified SNPs that have been robustly associated with lung cancer risk could also enter the model. Furthermore, decision trees were created to investigate variant × variant interaction. All analyses were performed for overall lung cancer and for subgroups.

RESULTS

No genome-wide significant association of AhR/Wnt-genes with overall lung cancer was observed, but within the subgroups of ever smokers (e.g., maker rs2722278 SFRP4; OR  = 1.20; 95% CI 1.13-1.27; p  = 5.6 × 10-10) and never smokers (e.g., maker rs1133683 Axin2; OR  = 1.27; 95% CI 1.19-1.35; p  = 1.0 × 10-12). Although predictability is poor, AhR/Wnt-variants are unexpectedly overrepresented in optimized prediction scores for overall lung cancer and for small cell lung cancer. Remarkably, the score for never-smokers contained solely two AhR/Wnt-variants. The optimal decision tree for never smokers consists of 7 AhR/Wnt-variants and only two lung cancer variants.

CONCLUSIONS

The role of variants belonging to Wnt/AhR-pathways in lung cancer susceptibility may be underrated in main-effects association analysis. Complex interaction patterns in individuals of European descent have moderate predictive capacity for lung cancer or subgroups thereof, especially in never smokers.

Iam hiQ-a novel pair of accuracy indices for imputed genotypes

BACKGROUND

Imputation of untyped markers is a standard tool in genome-wide association studies to close the gap between directly genotyped and other known DNA variants. However, high accuracy with which genotypes are imputed is fundamental. Several accuracy measures have been proposed and some are implemented in imputation software, unfortunately diversely across platforms. In the present paper, we introduce Iam hiQ, an independent pair of accuracy measures that can be applied to dosage files, the output of all imputation software. Iam (imputation accuracy measure) quantifies the average amount of individual-specific versus population-specific genotype information in a linear manner. hiQ (heterogeneity in quantities of dosages) addresses the inter-individual heterogeneity between dosages of a marker across the sample at hand.

RESULTS

Applying both measures to a large case-control sample of the International Lung Cancer Consortium (ILCCO), comprising 27,065 individuals, we found meaningful thresholds for Iam and hiQ suitable to classify markers of poor accuracy. We demonstrate how Manhattan-like plots and moving averages of Iam and hiQ can be useful to identify regions enriched with less accurate imputed markers, whereas these regions would by missed when applying the accuracy measure info (implemented in IMPUTE2).

CONCLUSION

We recommend using Iam hiQ additional to other accuracy scores for variant filtering before stepping into the analysis of imputed GWAS data.

Rare germline copy number variants (CNVs) and breast cancer risk

Germline copy number variants (CNVs) are pervasive in the human genome but potential disease associations with rare CNVs have not been comprehensively assessed in large datasets. We analysed rare CNVs in genes and non-coding regions for 86,788 breast cancer cases and 76,122 controls of European ancestry with genome-wide array data. Gene burden tests detected the strongest association for deletions in BRCA1 (P = 3.7E-18). Nine other genes were associated with a p-value < 0.01 including known susceptibility genes CHEK2 (P = 0.0008), ATM (P = 0.002) and BRCA2 (P = 0.008). Outside the known genes we detected associations with p-values < 0.001 for either overall or subtype-specific breast cancer at nine deletion regions and four duplication regions. Three of the deletion regions were in established common susceptibility loci. To the best of our knowledge, this is the first genome-wide analysis of rare CNVs in a large breast cancer case-control dataset. We detected associations with exonic deletions in established breast cancer susceptibility genes. We also detected suggestive associations with non-coding CNVs in known and novel loci with large effects sizes. Larger sample sizes will be required to reach robust levels of statistical significance.

Common variants in breast cancer risk loci predispose to distinct tumor subtypes

BACKGROUND

Genome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER) status, but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear.

METHODS

Among 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes.

RESULTS

Eighty-five of 173 variants were associated with at least one tumor feature (false discovery rate < 5%), most commonly ER and grade, followed by PR and HER2. Models for intrinsic-like subtypes found nearly all of these variants (83 of 85) associated at p < 0.05 with risk for at least one luminal-like subtype, and approximately half (41 of 85) of the variants were associated with risk of at least one non-luminal subtype, including 32 variants associated with triple-negative (TN) disease. Ten variants were associated with risk of all subtypes in different magnitude. Five variants were associated with risk of luminal A-like and TN subtypes in opposite directions.

CONCLUSION

This report demonstrates a high level of complexity in the etiology heterogeneity of breast cancer susceptibility variants and can inform investigations of subtype-specific risk prediction.

Germline variants and breast cancer survival in patients with distant metastases at primary breast cancer diagnosis

Breast cancer metastasis accounts for most of the deaths from breast cancer. Identification of germline variants associated with survival in aggressive types of breast cancer may inform understanding of breast cancer progression and assist treatment. In this analysis, we studied the associations between germline variants and breast cancer survival for patients with distant metastases at primary breast cancer diagnosis. We used data from the Breast Cancer Association Consortium (BCAC) including 1062 women of European ancestry with metastatic breast cancer, 606 of whom died of breast cancer. We identified two germline variants on chromosome 1, rs138569520 and rs146023652, significantly associated with breast cancer-specific survival (P = 3.19 × 10-8 and 4.42 × 10-8). In silico analysis suggested a potential regulatory effect of the variants on the nearby target genes SDE2 and H3F3A. However, the variants showed no evidence of association in a smaller replication dataset. The validation dataset was obtained from the SNPs to Risk of Metastasis (StoRM) study and included 293 patients with metastatic primary breast cancer at diagnosis. Ultimately, larger replication studies are needed to confirm the identified associations.