A prospective study of smoking-related white blood cell DNA methylation markers and risk of bladder cancer.

Bladder cancer, a common neoplasm, is primarily caused by tobacco smoking. Epigenetic alterations including DNA methylation have the potential to be used as prospective markers of increased risk, particularly in at-risk populations such as smokers. We aimed to investigate the potential of smoking-related white blood cell (WBC) methylation markers to contribute to an increase in bladder cancer risk prediction over classical questionnaire-based smoking metrics (i.e., duration, intensity, packyears) in a nested case-control study within the prospective prostate, lung, colorectal, and ovarian (PLCO) Cancer Screening Trial and the alpha-tocopherol, beta-carotene cancer (ATBC) Prevention Study (789 cases; 849 controls). We identified 200 differentially methylated sites associated with smoking status and 28 significantly associated (after correction for multiple testing) with bladder cancer risk among 2670 previously reported smoking-related cytosine-phosphate-guanines sites (CpGs). Similar patterns were observed across cohorts. Receiver operating characteristic (ROC) analyses indicated that cg05575921 (AHHR), the strongest smoking-related association we identified for bladder cancer risk, alone yielded similar predictive performance (AUC: 0.60) than classical smoking metrics (AUC: 0.59-0.62). Best prediction was achieved by including the first principal component (PC1) from the 200 smoking-related CpGs alongside smoking metrics (AUC: 0.63-0.65). Further, PC1 remained significantly associated with elevated bladder cancer risk after adjusting for smoking metrics. These findings suggest DNA methylation profiles reflect aspects of tobacco smoke exposure in addition to those captured by smoking duration, intensity and packyears, and/or individual susceptibility relevant to bladder cancer etiology, warranting further investigation.

Association of hair polychlorinated biphenyls and multiclass pesticides with obesity, diabetes, hypertension and dyslipidemia in NESCAV study.

Obesity, diabetes, hypertension and dyslipidemia are well-established risk factors for cardiovascular diseases (CVDs), and have been associated with exposure to persistent organic pollutants. However, studies have been lacking as regards effects of non-persistent pesticides on CVD risk factors. Here, we investigated whether background chronic exposure to polychlorinated biphenyls (PCBs) and multiclass pesticides were associated with the prevalence of these CVD risk factors in 502 Belgian and 487 Luxembourgish adults aged 18-69 years from the Nutrition, environment and cardiovascular health (NESCAV) study 2007-2013. We used hair analysis to evaluate the chronic internal exposure to three PCBs, seven organochlorine pesticides (OCs) and 18 non-persistent pesticides. We found positive associations of obesity with hexachlorobenzene (HCB), ß-hexachlorocyclohexane (ß-HCH) and chlorpyrifos, diabetes with pentachlorophenol (PCP), fipronil and fipronil sulfone, hypertension with PCB180 and chlorpyrifos, and dyslipidemia with diflufenican and oxadiazon, among others. However, we also found some inverse associations, such as obesity with PCP, diabetes with γ-HCH, hypertension with diflufenican, and dyslipidemia with chlorpyrifos. These results add to the existing evidence that OC exposure may contribute to the development of CVDs. Additionally, the present study revealed associations between CVD risk factors and chronic environmental exposure to currently used pesticides such as organophosphorus and pyrethroid pesticides.

Automated calibration for stability selection in penalised regression and graphical models.

Stability selection represents an attractive approach to identify sparse sets of features jointly associated with an outcome in high-dimensional contexts. We introduce an automated calibration procedure via maximisation of an in-house stability score and accommodating a priori-known block structure (e.g. multi-OMIC) data. It applies to [Least Absolute Shrinkage Selection Operator (LASSO)] penalised regression and graphical models. Simulations show our approach outperforms non-stability-based and stability selection approaches using the original calibration. Application to multi-block graphical LASSO on real (epigenetic and transcriptomic) data from the Norwegian Women and Cancer study reveals a central/credible and novel cross-OMIC role of LRRN3 in the biological response to smoking. Proposed approaches were implemented in the R package sharp.

Automated calibration of consensus weighted distance-based clustering approaches using sharp.

MOTIVATION: In consensus clustering, a clustering algorithm is used in combination with a subsampling procedure to detect stable clusters. Previous studies on both simulated and real data suggest that consensus clustering outperforms native algorithms. RESULTS: We extend here consensus clustering to allow for attribute weighting in the calculation of pairwise distances using existing regularized approaches. We propose a procedure for the calibration of the number of clusters (and regularization parameter) by maximizing the sharp score, a novel stability score calculated directly from consensus clustering outputs, making it extremely computationally competitive. Our simulation study shows better clustering performances of (i) approaches calibrated by maximizing the sharp score compared to existing calibration scores and (ii) weighted compared to unweighted approaches in the presence of features that do not contribute to cluster definition. Application on real gene expression data measured in lung tissue reveals clear clusters corresponding to different lung cancer subtypes. AVAILABILITY AND IMPLEMENTATION: The R package sharp (version ≥1.4.3) is available on CRAN at https://CRAN.R-project.org/package=sharp.

DNA methylation signature of chronic low-grade inflammation and its role in cardio-respiratory diseases.

We performed a multi-ethnic Epigenome Wide Association study on 22,774 individuals to describe the DNA methylation signature of chronic low-grade inflammation as measured by C-Reactive protein (CRP). We find 1,511 independent differentially methylated loci associated with CRP. These CpG sites show correlation structures across chromosomes, and are primarily situated in euchromatin, depleted in CpG islands. These genomic loci are predominantly situated in transcription factor binding sites and genomic enhancer regions. Mendelian randomization analysis suggests altered CpG methylation is a consequence of increased blood CRP levels. Mediation analysis reveals obesity and smoking as important underlying driving factors for changed CpG methylation. Finally, we find that an activated CpG signature significantly increases the risk for cardiometabolic diseases and COPD.

Epigenetic mechanisms of lung carcinogenesis involve differentially methylated CpG sites beyond those associated with smoking.

Smoking-related epigenetic changes have been linked to lung cancer, but the contribution of epigenetic alterations unrelated to smoking remains unclear. We sought for a sparse set of CpG sites predicting lung cancer and explored the role of smoking in these associations. We analysed CpGs in relation to lung cancer in participants from two nested case-control studies, using (LASSO)-penalised regression. We accounted for the effects of smoking using known smoking-related CpGs, and through conditional-independence network. We identified 29 CpGs (8 smoking-related, 21 smoking-unrelated) associated with lung cancer. Models additionally adjusted for Comprehensive Smoking Index-(CSI) selected 1 smoking-related and 49 smoking-unrelated CpGs. Selected CpGs yielded excellent discriminatory performances, outperforming information provided by CSI only. Of the 8 selected smoking-related CpGs, two captured lung cancer-relevant effects of smoking that were missed by CSI. Further, the 50 CpGs identified in the CSI-adjusted model complementarily explained lung cancer risk. These markers may provide further insight into lung cancer carcinogenesis and help improving early identification of high-risk patients.

Genetic analysis of over half a million people characterises C-reactive protein loci.

Chronic low-grade inflammation is linked to a multitude of chronic diseases. We report the largest genome-wide association study (GWAS) on C-reactive protein (CRP), a marker of systemic inflammation, in UK Biobank participants (N = 427,367, European descent) and the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium (total N = 575,531 European descent). We identify 266 independent loci, of which 211 are not previously reported. Gene-set analysis highlighted 42 gene sets associated with CRP levels (p ≤ 3.2 ×10-6) and tissue expression analysis indicated a strong association of CRP related genes with liver and whole blood gene expression. Phenome-wide association study identified 27 clinical outcomes associated with genetically determined CRP and subsequent Mendelian randomisation analyses supported a causal association with schizophrenia, chronic airway obstruction and prostate cancer. Our findings identified genetic loci and functional properties of chronic low-grade inflammation and provided evidence for causal associations with a range of diseases.

Predictive symptoms for COVID-19 in the community: REACT-1 study of over 1 million people.

BACKGROUND: Rapid detection, isolation, and contact tracing of community COVID-19 cases are essential measures to limit the community spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We aimed to identify a parsimonious set of symptoms that jointly predict COVID-19 and investigated whether predictive symptoms differ between the B.1.1.7 (Alpha) lineage (predominating as of April 2021 in the US, UK, and elsewhere) and wild type. METHODS AND FINDINGS: We obtained throat and nose swabs with valid SARS-CoV-2 PCR test results from 1,147,370 volunteers aged 5 years and above (6,450 positive cases) in the REal-time Assessment of Community Transmission-1 (REACT-1) study. This study involved repeated community-based random surveys of prevalence in England (study rounds 2 to 8, June 2020 to January 2021, response rates 22%-27%). Participants were asked about symptoms occurring in the week prior to testing. Viral genome sequencing was carried out for PCR-positive samples with N-gene cycle threshold value < 34 (N = 1,079) in round 8 (January 2021). In univariate analysis, all 26 surveyed symptoms were associated with PCR positivity compared with non-symptomatic people. Stability selection (1,000 penalized logistic regression models with 50% subsampling) among people reporting at least 1 symptom identified 7 symptoms as jointly and positively predictive of PCR positivity in rounds 2-7 (June to December 2020): loss or change of sense of smell, loss or change of sense of taste, fever, new persistent cough, chills, appetite loss, and muscle aches. The resulting model (rounds 2-7) predicted PCR positivity in round 8 with area under the curve (AUC) of 0.77. The same 7 symptoms were selected as jointly predictive of B.1.1.7 infection in round 8, although when comparing B.1.1.7 with wild type, new persistent cough and sore throat were more predictive of B.1.1.7 infection while loss or change of sense of smell was more predictive of the wild type. The main limitations of our study are (i) potential participation bias despite random sampling of named individuals from the National Health Service register and weighting designed to achieve a representative sample of the population of England and (ii) the necessary reliance on self-reported symptoms, which may be prone to recall bias and may therefore lead to biased estimates of symptom prevalence in England. CONCLUSIONS: Where testing capacity is limited, it is important to use tests in the most efficient way possible. We identified a set of 7 symptoms that, when considered together, maximize detection of COVID-19 in the community, including infection with the B.1.1.7 lineage.

Genetic variation in cervical preinvasive and invasive disease: a genome-wide association study.

BACKGROUND: Most uterine cervical high-risk human papillomavirus (HPV) infections are transient, with only a small fraction developing into cervical cancer. Family aggregation studies and heritability estimates suggest a significant inherited genetic component. Candidate gene studies and previous genome-wide association studies (GWASs) report associations between the HLA region and cervical cancer. Adopting a genome-wide approach, we aimed to compare genetic variation in women with invasive cervical cancer and cervical intraepithelial neoplasia (CIN) grade 3 with that in healthy controls. METHODS: We did a GWAS in a cohort of unrelated European individuals using data from UK Biobank, a population-based cohort including 273 377 women aged 40-69 years at recruitment between March 13, 2006, and Oct 1, 2010. We used an additive univariate logistic regression model to analyse genetic variants associated with invasive cervical cancer or CIN3. We sought replication of candidate associations in FinnGen, a large independent dataset of 128 123 individuals. We also did a two-sample mendelian randomisation approach to explore the role of risk factors in the genetic risk of cervical cancer. FINDINGS: We included 4769 CIN3 and invasive cervical cancer case samples and 145 545 control samples in the GWAS. Of 9 600 464 assayed and imputed single-nucleotide polymorphisms (SNPs), six independent variants were associated with CIN3 and invasive cervical cancer. These included novel loci rs10175462 (PAX8; odds ratio [OR] 0·87, 95% CI 0·84-0·91; p=1·07 × 10-9) and rs27069 (CLPTM1L; 0·88, 0·84-0·92; p=2·51 × 10-9), and previously reported signals at rs9272050 (HLA-DQA1; 1·27, 1·21-1·32; p=2·51 × 10-28), rs6938453 (MICA; 0·79, 0·75-0·83; p=1·97 × 10-17), rs55986091 (HLA-DQB1; 0·66, 0·60-0·72; p=6·42 × 10-28), and rs9266183 (HLA-B; 0·73, 0·64-0·83; p=1·53 × 10-6). Three SNPs were replicated in the independent Finnish dataset of 1648 invasive cervical cancer cases: PAX8 (rs10175462; p=0·015), CLPTM1L (rs27069; p=2·54 × 10-7), and HLA-DQA1 (rs9272050; p=7·90 × 10-8). Mendelian randomisation further supported the complementary role of smoking (OR 2·46, 95% CI 1·64-3·69), older age at first pregnancy (0·80, 0·68-0·95), and number of sexual partners (1·95, 1·44-2·63) in the risk of developing cervical cancer. INTERPRETATION: Our results provide new evidence for the genetic susceptibility to cervical cancer, specifically the PAX8, CLPTM1L, and HLA genes, suggesting disruption in apoptotic and immune function pathways. Future studies integrating host and viral, genetic, and epigenetic variation, could further elucidate complex host-viral interactions. FUNDING: NIHR Imperial BRC Wellcome 4i Clinician Scientist Training Programme.

COVID-19 mortality in the UK Biobank cohort: revisiting and evaluating risk factors.

Most studies of severe/fatal COVID-19 risk have used routine/hospitalisation data without detailed pre-morbid characterisation. Using the community-based UK Biobank cohort, we investigate risk factors for COVID-19 mortality in comparison with non-COVID-19 mortality. We investigated demographic, social (education, income, housing, employment), lifestyle (smoking, drinking, body mass index), biological (lipids, cystatin C, vitamin D), medical (comorbidities, medications) and environmental (air pollution) data from UK Biobank (N = 473,550) in relation to 459 COVID-19 and 2626 non-COVID-19 deaths to 21 September 2020. We used univariate, multivariable and penalised regression models. Age (OR = 2.76 [2.18-3.49] per S.D. [8.1 years], p = 2.6 × 10-17), male sex (OR = 1.47 [1.26-1.73], p = 1.3 × 10-6) and Black versus White ethnicity (OR = 1.21 [1.12-1.29], p = 3.0 × 10-7) were independently associated with and jointly explanatory of (area under receiver operating characteristic curve, AUC = 0.79) increased risk of COVID-19 mortality. In multivariable regression, alongside demographic covariates, being a healthcare worker, current smoker, having cardiovascular disease, hypertension, diabetes, autoimmune disease, and oral steroid use at enrolment were independently associated with COVID-19 mortality. Penalised regression models selected income, cardiovascular disease, hypertension, diabetes, cystatin C, and oral steroid use as jointly contributing to COVID-19 mortality risk; Black ethnicity, hypertension and oral steroid use contributed to COVID-19 but not non-COVID-19 mortality. Age, male sex and Black ethnicity, as well as comorbidities and oral steroid use at enrolment were associated with increased risk of COVID-19 death. Our results suggest that previously reported associations of COVID-19 mortality with body mass index, low vitamin D, air pollutants, renin-angiotensin-aldosterone system inhibitors may be explained by the aforementioned factors.

Prospective Identification of Elevated Circulating CDCP1 in Patients Years before Onset of Lung Cancer.

Increasing evidence points to a role for inflammation in lung carcinogenesis. A small number of circulating inflammatory proteins have been identified as showing elevated levels prior to lung cancer diagnosis, indicating the potential for prospective circulating protein concentration as a marker of early carcinogenesis. To identify novel markers of lung cancer risk, we measured a panel of 92 circulating inflammatory proteins in 648 prediagnostic blood samples from two prospective cohorts in Italy and Norway (women only). To preserve the comparability of results and protect against confounding factors, the main statistical analyses were conducted in women from both studies, with replication sought in men (Italian participants). Univariate and penalized regression models revealed for the first time higher blood levels of CDCP1 protein in cases that went on to develop lung cancer compared with controls, irrespective of time to diagnosis, smoking habits, and gender. This association was validated in an additional 450 samples. Associations were stronger for future cases of adenocarcinoma where CDCP1 showed better explanatory performance. Integrative analyses combining gene expression and protein levels of CDCP1 measured in the same individuals suggested a link between CDCP1 and the expression of transcripts of LRRN3 and SEM1. Enrichment analyses indicated a potential role for CDCP1 in pathways related to cell adhesion and mobility, such as the WNT/ß-catenin pathway. Overall, this study identifies lung cancer-related dysregulation of CDCP1 expression years before diagnosis. SIGNIFICANCE: Prospective proteomics analyses reveal an association between increased levels of circulating CDCP1 and lung carcinogenesis irrespective of smoking and years before diagnosis, and integrating gene expression indicates potential underlying mechanisms.See related commentary by Itzstein et al., p. 3441.

Cosmopolitan Distribution of Endozoicomonas-Like Organisms and Other Intracellular Microcolonies of Bacteria Causing Infection in Marine Mollusks.

Intracellular microcolonies of bacteria (IMC), in some cases developing large extracellular cysts (bacterial aggregates), infecting primarily gill and digestive gland, have been historically reported in a wide diversity of economically important mollusk species worldwide, sometimes associated with severe lesions and mass mortality events. As an effort to characterize those organisms, traditionally named as Rickettsia or Chlamydia-like organisms, 1950 specimens comprising 22 mollusk species were collected over 10 countries and after histology examination, a selection of 99 samples involving 20 species were subjected to 16S rRNA gene amplicon sequencing. Phylogenetic analysis showed Endozoicomonadaceae sequences in all the mollusk species analyzed. Geographical differences in the distribution of Operational Taxonomic Units (OTUs) and a particular OTU associated with pathology in king scallop (OTU_2) were observed. The presence of Endozoicomonadaceae sequences in the IMC was visually confirmed by in situ hybridization (ISH) in eight selected samples. Sequencing data also indicated other symbiotic bacteria. Subsequent phylogenetic analysis of those OTUs revealed a novel microbial diversity associated with molluskan IMC infection distributed among different taxa, including the phylum Spirochetes, the families Anaplasmataceae and Simkaniaceae, the genera Mycoplasma and Francisella, and sulfur-oxidizing endosymbionts. Sequences like Francisella halioticida/philomiragia and Candidatus Brownia rhizoecola were also obtained, however, in the absence of ISH studies, the association between those organisms and the IMCs were not confirmed. The sequences identified in this study will allow for further molecular characterization of the microbial community associated with IMC infection in marine mollusks and their correlation with severity of the lesions to clarify their role as endosymbionts, commensals or true pathogens.

Education, biological ageing, all-cause and cause-specific mortality and morbidity: UK biobank cohort study.

BACKGROUND: Socioeconomic position as measured by education may be embodied and affect the functioning of key physiological systems. Links between social disadvantage, its biological imprint, and cause-specific mortality and morbidity have not been investigated in large populations, and yet may point towards areas for public health interventions beyond targeting individual behaviours. METHODS: Using data from 366,748 UK Biobank participants with 13 biomarker measurements, we calculated a Biological Health Score (BHS, ranging from 0 to 1) capturing the level of functioning of five physiological systems. Associations between BHS and incidence of cardiovascular disease (CVD) and cancer, and mortality from all, CVD, cancer, and external causes were examined. We explored the role of education in these associations. Mendelian randomisation using genetic evidence was used to triangulate these findings. FINDINGS: An increase in BHS of 0.1 was associated with all-cause (HR = 1.14 [1.12-1.16] and 1.09 [1.07-1.12] in men and women respectively), cancer (HR = 1.11 [1.09-1.14] and 1.07 [1.04-1.10]) and CVD (HR = 1.25 [1.20-1.31] and 1.21 [1.11-1.31]) mortality, CVD incidence (HR = 1.15 [1.13-1.16] and 1.17 [1.15-1.19]). These associations survived adjustment for education, lifestyle-behaviours, body mass index (BMI), co-morbidities and medical treatments. Mendelian randomisation further supported the link between the BHS and CVD incidence (HR = 1.31 [1.21-1.42]). The BHS contributed to CVD incidence prediction (age-adjusted C-statistic = 0.58), other than through education and health behaviours. INTERPRETATION: The BHS captures features of the embodiment of education, health behaviours, and more proximal unknown factors which all complementarily contribute to all-cause, cancer and CVD morbidity and premature death.

Agnostic Cys34-albumin adductomics and DNA methylation: Implication of N-acetylcysteine in lung carcinogenesis years before diagnosis.

Although smoking and oxidative stress are known contributors to lung carcinogenesis, their mechanisms of action remain poorly understood. To shed light into these mechanisms, we applied a novel approach using Cys34-adductomics in a lung cancer nested case-control study (n = 212). Adductomics profiles were integrated with DNA-methylation data at established smoking-related CpG sites measured in the same individuals. Our analysis identified 42 Cys34-albumin adducts, of which 2 were significantly differentially abundant in cases and controls: adduct of N-acetylcysteine (NAC, p = 4.15 × 10-3 ) and of cysteinyl-glycine (p = 7.89 × 10-3 ). Blood levels of the former were found associated to the methylation levels at 11 smoking-related CpG sites. We detect, for the first time in prospective blood samples, and irrespective of time to diagnosis, decreased levels of NAC adduct in lung cancer cases. Altogether, our results highlight the potential role of these adducts in the oxidative stress response contributing to lung carcinogenesis years before diagnosis.

Pre-diagnostic blood immune markers, incidence and progression of B-cell lymphoma and multiple myeloma: Univariate and functionally informed multivariate analyses.

Recent prospective studies have shown that dysregulation of the immune system may precede the development of B-cell lymphomas (BCL) in immunocompetent individuals. However, to date, the studies were restricted to a few immune markers, which were considered separately. Using a nested case-control study within two European prospective cohorts, we measured plasma levels of 28 immune markers in samples collected a median of 6 years before diagnosis (range 2.01-15.97) in 268 incident cases of BCL (including multiple myeloma [MM]) and matched controls. Linear mixed models and partial least square analyses were used to analyze the association between levels of immune marker and the incidence of BCL and its main histological subtypes and to investigate potential biomarkers predictive of the time to diagnosis. Linear mixed model analyses identified associations linking lower levels of fibroblast growth factor-2 (FGF-2 p = 7.2 × 10-4 ) and transforming growth factor alpha (TGF-α, p = 6.5 × 10-5 ) and BCL incidence. Analyses stratified by histological subtypes identified inverse associations for MM subtype including FGF-2 (p = 7.8 × 10-7 ), TGF-α (p = 4.08 × 10-5 ), fractalkine (p = 1.12 × 10-3 ), monocyte chemotactic protein-3 (p = 1.36 × 10-4 ), macrophage inflammatory protein 1-alpha (p = 4.6 × 10-4 ) and vascular endothelial growth factor (p = 4.23 × 10-5 ). Our results also provided marginal support for already reported associations between chemokines and diffuse large BCL (DLBCL) and cytokines and chronic lymphocytic leukemia (CLL). Case-only analyses showed that Granulocyte-macrophage colony stimulating factor levels were consistently higher closer to diagnosis, which provides further evidence of its role in tumor progression. In conclusion, our study suggests a role of growth-factors in the incidence of MM and of chemokine and cytokine regulation in DLBCL and CLL.