Validation of the BOADICEA model in a prospective cohort of BRCA1/2 pathogenic variant carriers.

BACKGROUND: No validation has been conducted for the BOADICEA multifactorial breast cancer risk prediction model specifically in BRCA1/2 pathogenic variant (PV) carriers to date. Here, we evaluated the performance of BOADICEA in predicting 5-year breast cancer risks in a prospective cohort of BRCA1/2 PV carriers ascertained through clinical genetic centres. METHODS: We evaluated the model calibration and discriminatory ability in the prospective TRANsIBCCS cohort study comprising 1614 BRCA1 and 1365 BRCA2 PV carriers (209 incident cases). Study participants had lifestyle, reproductive, hormonal, anthropometric risk factor information, a polygenic risk score based on 313 SNPs and family history information. RESULTS: The full multifactorial model considering family history together with all other risk factors was well calibrated overall (E/O=1.07, 95% CI: 0.92 to 1.24) and in quintiles of predicted risk. Discrimination was maximised when all risk factors were considered (Harrell's C-index=0.70, 95% CI: 0.67 to 0.74; area under the curve=0.79, 95% CI: 0.76 to 0.82). The model performance was similar when evaluated separately in BRCA1 or BRCA2 PV carriers. The full model identified 5.8%, 12.9% and 24.0% of BRCA1/2 PV carriers with 5-year breast cancer risks of <1.65%, <3% and <5%, respectively, risk thresholds commonly used for different management and risk-reduction options. CONCLUSION: BOADICEA may be used to aid personalised cancer risk management and decision-making for BRCA1 and BRCA2 PV carriers. It is implemented in the free-access CanRisk tool (https://www.canrisk.org/).

Safety of the Breast Cancer Adjuvant Radiotherapy in Ataxia-Telangiectasia Mutated Variant Carriers.

The Ataxia-Telangiectasia Mutated (ATM) gene is implicated in DNA double-strand break repair. Controversies in clinical radiosensitivity remain known for monoallelic carriers of the ATM pathogenic variant (PV). An evaluation of the single-nucleotide polymorphism (SNP) rs1801516 (G-A) showed different results regarding late subcutaneous fibrosis after breast radiation therapy (RT). The main objective of this study was to evaluate acute and late toxicities in carriers of a rare ATM PV or predicted PV and in carriers of minor allele A of rs1801516 facing breast RT. Fifty women with localized breast cancer treated with adjuvant RT between 2000 and 2014 at Institut Curie were selected. Acute and late toxicities in carriers of a rare PV or predicted PV (n= 9), in noncarriers (n = 41) and in carriers of SNP rs1801516 (G-A) (n = 8), were examined. The median age at diagnosis was 53 years old and 82% of patients had an invasive ductal carcinoma and 84% were at clinical stage I-IIB. With a median follow-up of 13 years, no significant difference between carriers and noncarriers was found for acute toxicities (p > 0.05). The same results were observed for late toxicities without an effect from the rs1801516 genotype on toxicities. No significant difference in acute or late toxicities was observed between rare ATM variant carriers and noncarriers after breast RT for localized breast cancer.

Characteristics and correlates of seclusion and mechanical restraint measures in a Parisian psychiatric hospital group.

Introduction: Seclusion or restraint (S/R) are last-resort measures used in psychiatry to ensure the safety of the patient and the staff. However, they have harmful physical and psychological effects on patients, and efforts to limit their use are needed. We describe the characteristics and correlates of S/R events in four Parisian psychiatric centers. Methods: Within a 3-month period, November 5, 2018 to February 3, 2019, we recorded data for patients experiencing an S/R measure as well as characteristics of the measures. We studied the mean duration of a S/R event, the time between hospital admission and the occurrence of the event, as well as correlates of these durations. We also examined factors associated with use of a restraint versus a seclusion measure. Results: For the 233 patients included, we recorded 217 seclusion measures and 64 mechanical restraints. Seclusion measures mostly occurred after the patient's transfer from the emergency department. The duration of a seclusion measure was about 10 days. Patients considered resistant to psychotropic treatments more frequently had a longer seclusion duration than others. The mean duration of a mechanical restraint measure was 4 days. Male sex and younger age were associated with experiencing mechanical restraint. Discussion: S/R measures mostly occur among patients perceived as resistant to psychotropic drugs who are arriving from the emergency department. Developing specific emergency department protocols might be useful in limiting the use of coercive measures.

Do Determinants of Smoking Cessation and Relapse Differ between Men and Women? Data from a French National Study.

Introduction: There is limited understanding of different predictors of smoking cessation success (SCS) among women and men, despite well-documented differences in smoking behavior.Methods: Using data from DePICT (Description des Perceptions, Images, et Comportements liés au Tabagisme), a national survey of French adults which recruited 2377 regular and former smokers we investigated whether major determinants of SCS differed by sex. Factors associated with unsuccessful vs. No successful quit attempt; vs. SCS were studied using multivariate multinomial logistic regression analyses stratified by sex.Results: Women and men share some determinants of SCS including no cannabis use, living in a nonsmoker household and importance giving to being a nonsmoker. However, no e-cigarette use, low-to-moderate alcohol consumption, early smoking initiation, and higher education were associated with SCS only among women. No use of nicotine replacement, having family members who smoke, family opinion on smoking and current employment, were associated with SCS only among men. Neutral or negative friends' opinion on smoking or living with a smoker were associated with unsuccessful smoking attempts among men.Conclusions: Our results show differences between determinants of SCS according to sex, which highlights the importance of developing tailored interventions that account for sex/gender differences in smoking cessation.

Investigation of common genetic risk factors between thyroid traits and breast cancer.

Breast cancer (BC) risk is suspected to be linked to thyroid disorders, however observational studies exploring the association between BC and thyroid disorders gave conflicting results. We proposed an alternative approach by investigating the shared genetic risk factors between BC and several thyroid traits. We report a positive genetic correlation between BC and thyroxine (FT4) levels (corr = 0.13, p-value = 2.0 × 10-4) and a negative genetic correlation between BC and thyroid-stimulating hormone (TSH) levels (corr = -0.09, p-value = 0.03). These associations are more striking when restricting the analysis to estrogen receptor-positive BC. Moreover, the polygenic risk scores (PRS) for FT4 and hyperthyroidism are positively associated to BC risk (OR = 1.07, 95%CI: 1.00-1.13, p-value = 2.8 × 10-2 and OR = 1.04, 95%CI: 1.00-1.08, p-value = 3.8 × 10-2, respectively), while the PRS for TSH is inversely associated to BC risk (OR = 0.93, 95%CI: 0.89-0.97, p-value = 2.0 × 10-3). Using the PLACO method, we detected 49 loci associated to both BC and thyroid traits (p-value < 5 × 10-8), in the vicinity of 130 genes. An additional colocalization and gene-set enrichment analyses showed a convincing causal role for a known pleiotropic locus at 2q35 and revealed an additional one at 8q22.1 associated to both BC and thyroid cancer. We also found two new pleiotropic loci at 14q32.33 and 17q21.31 that were associated to both TSH levels and BC risk. Enrichment analyses and evidence of regulatory signals also highlighted brain tissues and immune system as candidates for obtaining associations between BC and TSH levels. Overall, our study sheds light on the complex interplay between BC and thyroid traits and provides evidence of shared genetic risk between those conditions.

Evolution of smoking rates among immigrants in France in the context of comprehensive tobacco control measures, and a decrease in the overall prevalence.

BACKGROUND: The evolution of smoking rates according to migrant status has not been examined in France, despite a recent reduction in overall smoking rates. METHODS: DePICT is a two waves (2016: n = 4356; 2017: n = 4114) nationwide telephone survey, representative of the French adult population. We compared smoking-related behaviors before and after implementation of tobacco-control measures (2017), according to the geographical region of birth. RESULTS: Compared to 2016, individuals originating from Africa or the Middle East had a slightly higher smoking prevalence in 2017 (34.7% vs 31.3%), despite a higher intention to quit or attempt in the preceding year (adjusted OR(ORa) = 2.72[1.90; 3.90]) compared to non-immigrants. They were also less likely to experience an unsuccessful quit attempt (ORa = 1.76[1.18; 2.62]). CONCLUSION: Tobacco-control measures could have widened smoking inequalities related to migrant status. The evolution of smoking-related behaviors among immigrants should be examined when studying the long-term effects of such policies.

Investigation of Shared Genetic Risk Factors Between Parkinson's Disease and Cancers.

BACKGROUND: Epidemiological studies that examined the association between Parkinson's disease (PD) and cancers led to inconsistent results, but they face a number of methodological difficulties. OBJECTIVE: We used results from genome-wide association studies (GWASs) to study the genetic correlation between PD and different cancers to identify common genetic risk factors. METHODS: We used individual data for participants of European ancestry from the Courage-PD (Comprehensive Unbiased Risk Factor Assessment for Genetics and Environment in Parkinson's Disease; PD, N = 16,519) and EPITHYR (differentiated thyroid cancer, N = 3527) consortia and summary statistics of GWASs from iPDGC (International Parkinson Disease Genomics Consortium; PD, N = 482,730), Melanoma Meta-Analysis Consortium (MMAC), Breast Cancer Association Consortium (breast cancer), the Prostate Cancer Association Group to Investigate Cancer Associated Alterations in the Genome (prostate cancer), International Lung Cancer Consortium (lung cancer), and Ovarian Cancer Association Consortium (ovarian cancer) (N comprised between 36,017 and 228,951 for cancer GWASs). We estimated the genetic correlation between PD and cancers using linkage disequilibrium score regression. We studied the association between PD and polymorphisms associated with cancers, and vice versa, using cross-phenotypes polygenic risk score (PRS) analyses. RESULTS: We confirmed a previously reported positive genetic correlation of PD with melanoma (Gcorr = 0.16 [0.04; 0.28]) and reported an additional significant positive correlation of PD with prostate cancer (Gcorr = 0.11 [0.03; 0.19]). There was a significant inverse association between the PRS for ovarian cancer and PD (odds ratio [OR] = 0.89 [0.84; 0.94]). Conversely, the PRS of PD was positively associated with breast cancer (OR = 1.08 [1.06; 1.10]) and inversely associated with ovarian cancer (OR = 0.95 [0.91; 0.99]). The association between PD and ovarian cancer was mostly driven by rs183211 located in an intron of the NSF gene (17q21.31). CONCLUSIONS: We show evidence in favor of a contribution of pleiotropic genes to the association between PD and specific cancers. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

A pleiotropic variant in DNAJB4 is associated with multiple myeloma risk.

Pleiotropy, which consists of a single gene or allelic variant affecting multiple unrelated traits, is common across cancers, with evidence for genome-wide significant loci shared across cancer and noncancer traits. This feature is particularly relevant in multiple myeloma (MM) because several susceptibility loci that have been identified to date are pleiotropic. Therefore, the aim of this study was to identify novel pleiotropic variants involved in MM risk using 28 684 independent single nucleotide polymorphisms (SNPs) from GWAS Catalog that reached a significant association (P < 5 × 10-8 ) with their respective trait. The selected SNPs were analyzed in 2434 MM cases and 3446 controls from the International Lymphoma Epidemiology Consortium (InterLymph). The 10 SNPs showing the strongest associations with MM risk in InterLymph were selected for replication in an independent set of 1955 MM cases and 1549 controls from the International Multiple Myeloma rESEarch (IMMEnSE) consortium and 418 MM cases and 147 282 controls from the FinnGen project. The combined analysis of the three studies identified an association between DNAJB4-rs34517439-A and an increased risk of developing MM (OR = 1.22, 95%CI 1.13-1.32, P = 4.81 × 10-7 ). rs34517439-A is associated with a modified expression of the FUBP1 gene, which encodes a multifunctional DNA and RNA-binding protein that it was observed to influence the regulation of various genes involved in cell cycle regulation, among which various oncogenes and oncosuppressors. In conclusion, with a pleiotropic scan approach we identified DNAJB4-rs34517439 as a potentially novel MM risk locus.

A polygenic risk score for multiple myeloma risk prediction.

There is overwhelming epidemiologic evidence that the risk of multiple myeloma (MM) has a solid genetic background. Genome-wide association studies (GWAS) have identified 23 risk loci that contribute to the genetic susceptibility of MM, but have low individual penetrance. Combining the SNPs in a polygenic risk score (PRS) is a possible approach to improve their usefulness. Using 2361 MM cases and 1415 controls from the International Multiple Myeloma rESEarch (IMMEnSE) consortium, we computed a weighted and an unweighted PRS. We observed associations with MM risk with OR = 3.44, 95% CI 2.53-4.69, p = 3.55 × 10-15 for the highest vs. lowest quintile of the weighted score, and OR = 3.18, 95% CI 2.1 = 34-4.33, p = 1.62 × 10-13 for the highest vs. lowest quintile of the unweighted score. We found a convincing association of a PRS generated with 23 SNPs and risk of MM. Our work provides additional validation of previously discovered MM risk variants and of their combination into a PRS, which is a first step towards the use of genetics for risk stratification in the general population.

First international workshop of the ATM and cancer risk group (4-5 December 2019).

The first International Workshop of the ATM and Cancer Risk group focusing on the role of Ataxia-Telangiectasia Mutated (ATM) gene in cancer was held on December 4 and 5, 2019 at Institut Curie in Paris, France. It was motivated by the fact that germline ATM pathogenic variants have been found to be associated with different cancer types. However, due to the lack of precise age-, sex-, and site-specific risk estimates, no consensus on management guidelines for variant carriers exists, and the clinical utility of ATM variant testing is uncertain. The meeting brought together epidemiologists, geneticists, biologists and clinicians to review current knowledge and on-going challenges related to ATM and cancer risk. This report summarizes the meeting sessions content that covered the latest results in family-based and population-based studies, the importance of accurate variant classification, the effect of radiation exposures for ATM variant carriers, and the characteristics of ATM-deficient tumors. The report concludes that ATM variant carriers outside of the context of Ataxia-Telangiectasia may benefit from effective cancer risk management and therapeutic strategies and that efforts to set up large-scale studies in the international framework to achieve this goal are necessary.

TUMOSPEC: A Nation-Wide Study of Hereditary Breast and Ovarian Cancer Families with a Predicted Pathogenic Variant Identified through Multigene Panel Testing.

Assessment of age-dependent cancer risk for carriers of a predicted pathogenic variant (PPV) is often hampered by biases in data collection, with a frequent under-representation of cancer-free PPV carriers. TUMOSPEC was designed to estimate the cumulative risk of cancer for carriers of a PPV in a gene that is usually tested in a hereditary breast and ovarian cancer context. Index cases are enrolled consecutively among patients who undergo genetic testing as part of their care plan in France. First- and second-degree relatives and cousins of PPV carriers are invited to participate whether they are affected by cancer or not, and genotyped for the familial PPV. Clinical, family and epidemiological data are collected, and all data including sequencing data are centralized at the coordinating centre. The three-year feasibility study included 4431 prospective index cases, with 19.1% of them carrying a PPV. When invited by the coordinating centre, 65.3% of the relatives of index cases (5.7 relatives per family, on average) accepted the invitation to participate. The study logistics were well adapted to clinical and laboratory constraints, and collaboration between partners (clinicians, biologists, coordinating centre and participants) was smooth. Hence, TUMOSPEC is being pursued, with the aim of optimizing clinical management guidelines specific to each gene.

Diagnostic chest X-rays and breast cancer risk among women with a hereditary predisposition to breast cancer unexplained by a BRCA1 or BRCA2 mutation.

BACKGROUND: Diagnostic ionizing radiation is a risk factor for breast cancer (BC). BC risk increases with increased dose to the chest and decreases with increased age at exposure, with possible effect modification related to familial or genetic predisposition. While chest X-rays increase the BC risk of BRCA1/2 mutation carriers compared to non-carriers, little is known for women with a hereditary predisposition to BC but who tested negative for a BRCA1 or BRCA2 (BRCA1/2) mutation. METHODS: We evaluated the effect of chest X-rays from diagnostic medical procedures in a dataset composed of 1552 BC cases identified through French family cancer clinics and 1363 unrelated controls. Participants reported their history of X-ray exposures in a detailed questionnaire and were tested for 113 DNA repair genes. Logistic regression and multinomial logistic regression models were used to assess the association with BC. RESULTS: Chest X-ray exposure doubled BC risk. A 3% increased BC risk per additional exposure was observed. Being 20 years old or younger at first exposure or being exposed before first full-term pregnancy did not seem to modify this risk. Birth after 1960 or carrying a rare likely deleterious coding variant in a DNA repair gene other than BRCA1/2 modified the effect of chest X-ray exposure. CONCLUSION: Ever/never chest X-ray exposure increases BC risk 2-fold regardless of age at first exposure and, by up to 5-fold when carrying 3 or more rare variants in a DNA repair gene. Further studies are needed to evaluate other DNA repair genes or variants to identify those which could modify radiation sensitivity. Identification of subpopulations that are more or less susceptible to ionizing radiation is important and potentially clinically relevant.

Smoking cessation using preference-based tools: a mixed method pilot study of a novel intervention among smokers with low socioeconomic position.

BACKGROUND: Compared to smokers with favorable socio-economic position (SEP), those with low SEP are less likely to have a successful smoking cessation attempt. Tailored approaches are therefore needed, and general practitioners could help reaching and assisting usually hard-to-reach population. METHOD: STOP (Sevrage Tabagique à l'aide d'Outils dédiés selon la Préférence) is a pilot study, examining the feasibility, acceptability and potentiality of a smoking cessation intervention centered on smoker's preference. Smokers with low SEP, wishing to quit, were recruited in six healthcare centers in the Greater Paris area. They were asked to choose between different types of nicotine replacement therapy (NRT) products and/or e-cigarette with liquids delivered free of charge to aid their smoking cessation attempt. We describe the characteristics of recruited participants, their perception of smoking cessation aids, and the evolution of their smoking status 4 to 6 weeks after recruitment. RESULTS: We recruited 49 participants, of which 29% chose an e-cigarette, 29% chose NRT and 42% chose both an e-cigarette and NRT. The intervention was shown to be acceptable by participants and health professionals. Among the 24 participants followed for at least one month, 14 (28% of all participants) stopped smoking, and 9 (18%) considerably reduced their consumption. CONCLUSION: The STOP intervention is feasible and acceptable, even if more efforts should be made to limit lost-to-follow-up. This preference-based intervention also shows interesting prospect in helping smokers with low SEP quit smoking. We will test the efficacy of this preference-based intervention in a randomized controlled trial.

A new hybrid record linkage process to make epidemiological databases interoperable: application to the GEMO and GENEPSO studies involving BRCA1 and BRCA2 mutation carriers.

BACKGROUND: Linking independent sources of data describing the same individuals enable innovative epidemiological and health studies but require a robust record linkage approach. We describe a hybrid record linkage process to link databases from two independent ongoing French national studies, GEMO (Genetic Modifiers of BRCA1 and BRCA2), which focuses on the identification of genetic factors modifying cancer risk of BRCA1 and BRCA2 mutation carriers, and GENEPSO (prospective cohort of BRCAx mutation carriers), which focuses on environmental and lifestyle risk factors. METHODS: To identify as many as possible of the individuals participating in the two studies but not registered by a shared identifier, we combined probabilistic record linkage (PRL) and supervised machine learning (ML). This approach (named "PRL + ML") combined together the candidate matches identified by both approaches. We built the ML model using the gold standard on a first version of the two databases as a training dataset. This gold standard was obtained from PRL-derived matches verified by an exhaustive manual review. Results The Random Forest (RF) algorithm showed a highest recall (0.985) among six widely used ML algorithms: RF, Bagged trees, AdaBoost, Support Vector Machine, Neural Network. Therefore, RF was selected to build the ML model since our goal was to identify the maximum number of true matches. Our combined linkage PRL + ML showed a higher recall (range 0.988-0.992) than either PRL (range 0.916-0.991) or ML (0.981) alone. It identified 1995 individuals participating in both GEMO (6375 participants) and GENEPSO (4925 participants). CONCLUSIONS: Our hybrid linkage process represents an efficient tool for linking GEMO and GENEPSO. It may be generalizable to other epidemiological studies involving other databases and registries.

Evaluation of the association of heterozygous germline variants in NTHL1 with breast cancer predisposition: an international multi-center study of 47,180 subjects.

Bi-allelic loss-of-function (LoF) variants in the base excision repair (BER) gene NTHL1 cause a high-risk hereditary multi-tumor syndrome that includes breast cancer, but the contribution of heterozygous variants to hereditary breast cancer is unknown. An analysis of 4985 women with breast cancer, enriched for familial features, and 4786 cancer-free women revealed significant enrichment for NTHL1 LoF variants. Immunohistochemistry confirmed reduced NTHL1 expression in tumors from heterozygous carriers but the NTHL1 bi-allelic loss characteristic mutational signature (SBS 30) was not present. The analysis was extended to 27,421 breast cancer cases and 19,759 controls from 10 international studies revealing 138 cases and 93 controls with a heterozygous LoF variant (OR 1.06, 95% CI: 0.82-1.39) and 316 cases and 179 controls with a missense variant (OR 1.31, 95% CI: 1.09-1.57). Missense variants selected for deleterious features by a number of in silico bioinformatic prediction tools or located within the endonuclease III functional domain showed a stronger association with breast cancer. Somatic sequencing of breast cancers from carriers indicated that the risk associated with NTHL1 appears to operate through haploinsufficiency, consistent with other described low-penetrance breast cancer genes. Data from this very large international multicenter study suggests that heterozygous pathogenic germline coding variants in NTHL1 may be associated with low- to moderate- increased risk of breast cancer.

Genetically determined telomere length and multiple myeloma risk and outcome.

Telomeres are involved in processes like cellular growth, chromosomal stability, and proper segregation to daughter cells. Telomere length measured in leukocytes (LTL) has been investigated in different cancer types, including multiple myeloma (MM). However, LTL measurement is prone to heterogeneity due to sample handling and study design (retrospective vs. prospective). LTL is genetically determined; genome-wide association studies identified 11 SNPs that, combined in a score, can be used as a genetic instrument to measure LTL and evaluate its association with MM risk. This approach has been already successfully attempted in various cancer types but never in MM. We tested the "teloscore" in 2407 MM patients and 1741 controls from the International Multiple Myeloma rESEarch (IMMeNSE) consortium. We observed an increased risk for longer genetically determined telomere length (gdTL) (OR = 1.69; 95% CI 1.36-2.11; P = 2.97 × 10-6 for highest vs. lowest quintile of the score). Furthermore, in a subset of 1376 MM patients we tested the relationship between the teloscore and MM patients survival, observing a better prognosis for longer gdTL compared with shorter gdTL (HR = 0.93; 95% CI 0.86-0.99; P = 0.049). In conclusion, we report convincing evidence that longer gdTL is a risk marker for MM risk, and that it is potentially involved in increasing MM survival.

Role of DNA Repair Variants and Diagnostic Radiology Exams in Differentiated Thyroid Cancer Risk: A Pooled Analysis of Two Case-Control Studies.

BACKGROUND: Given the increased use and diversity of diagnostic procedures, it is important to understand genetic susceptibility to radiation-induced thyroid cancer. METHODS: On the basis of self-declared diagnostic radiology examination records in addition to existing literature, we estimated the radiation dose delivered to the thyroid gland from diagnostic procedures during childhood and adulthood in two case-control studies conducted in France. A total of 1,071 differentiated thyroid cancer (DTC) cases and 1,188 controls from the combined studies were genotyped using a custom-made Illumina OncoArray DNA chip. We focused our analysis on variants in genes involved in DNA damage response and repair pathways, representing a total of 5,817 SNPs in 571 genes. We estimated the OR per milli-Gray (OR/mGy) of the radiation dose delivered to the thyroid gland using conditional logistic regression. We then used an unconditional logistic regression model to assess the association between DNA repair gene variants and DTC risk. We performed a meta-analysis of the two studies. RESULTS: The OR/mGy was 1.02 (95% confidence interval, 1.00-1.03). We found significant associations between DTC and rs7164173 in CHD2 (P = 5.79 × 10-5), rs6067822 in NFATc2 (P = 9.26 × 10-5), rs1059394 and rs699517 both in ENOSF1/THYS, rs12702628 in RPA3, and an interaction between rs7068306 in MGMT and thyroid radiation doses (P = 3.40 × 10-4). CONCLUSIONS: Our results suggest a role for variants in CDH2, NFATc2, ENOSF1/THYS, RPA3, and MGMT in DTC risk. IMPACT: CDH2, NFATc2, ENOSF1/THYS, and RPA3 have not previously been shown to be associated with DTC risk.

Gene network and biological pathways associated with susceptibility to differentiated thyroid carcinoma.

Variants identified in earlier genome-wide association studies (GWAS) on differentiated thyroid carcinoma (DTC) explain about 10% of the overall estimated genetic contribution and could not provide complete insights into biological mechanisms involved in DTC susceptibility. Integrating systems biology information from model organisms, genome-wide expression data from tumor and matched normal tissue and GWAS data could help identifying DTC-associated genes, and pathways or functional networks in which they are involved. We performed data mining of GWAS data of the EPITHYR consortium (1551 cases and 1957 controls) using various pathways and protein-protein interaction (PPI) annotation databases and gene expression data from The Cancer Genome Atlas. We identified eight DTC-associated genes at known loci 2q35 (DIRC3), 8p12 (NRG1), 9q22 (FOXE1, TRMO, HEMGN, ANP32B, NANS) and 14q13 (MBIP). Using the EW_dmGWAS approach we found that gene networks related to glycogenolysis, glycogen metabolism, insulin metabolism and signal transduction pathways associated with muscle contraction were overrepresented with association signals (false discovery rate adjusted p-value < 0.05). Additionally, suggestive association of 21 KEGG and 75 REACTOME pathways with DTC indicate a link between DTC susceptibility and functions related to metabolism of cholesterol, amino sugar and nucleotide sugar metabolism, steroid biosynthesis, and downregulation of ERBB2 signaling pathways. Together, our results provide novel insights into biological mechanisms contributing to DTC risk.

Fine-mapping of two differentiated thyroid carcinoma susceptibility loci at 2q35 and 8p12 in Europeans, Melanesians and Polynesians.

Differentiated thyroid carcinoma (DTC) incidence is characterized by wide ethnic and geographic variations, with high incidence rates observed in Oceanian populations. Genome-wide association studies (GWAS) identified mainly four DTC susceptibility loci at 9q22.33, 14q13.3, 2q35 and 8p12. Here we performed fine-mapping of the 2q35 and 8p12 loci in the population of the EPITHYR consortium that includes Europeans, Melanesians and Polynesians to identify likely causal variants for DTC risk. We conducted a colocalization analysis using eQTLs data to determine the SNPs with the highest probability of causality. At 2q35, we highlighted rs16857609 located in DIRC3. This SNP has a high probability of causality in the three populations, and a significant association in Europeans (OR = 1.4, p = 1.9 x 10-10). It is also associated with expression of DIRC3 and of the nearby gene IGFBP5 in thyroid tumour cells. At 8p12, we identified rs7844425 which was significantly associated with DTC in Europeans (OR = 1.32, p = 7.6 x 10-8) and rs2439304, which was highlighted by the colocalization analysis but only moderately associated with DTC in our dataset (OR = 1.2, p = 0.001). These SNPs are linked to the expression of NRG1 in thyroid tissue. Hence, our study identified novel variants at 2q35 and 8p12 to be prioritized for further functional studies.

Boosting GWAS using biological networks: A study on susceptibility to familial breast cancer.

Genome-wide association studies (GWAS) explore the genetic causes of complex diseases. However, classical approaches ignore the biological context of the genetic variants and genes under study. To address this shortcoming, one can use biological networks, which model functional relationships, to search for functionally related susceptibility loci. Many such network methods exist, each arising from different mathematical frameworks, pre-processing steps, and assumptions about the network properties of the susceptibility mechanism. Unsurprisingly, this results in disparate solutions. To explore how to exploit these heterogeneous approaches, we selected six network methods and applied them to GENESIS, a nationwide French study on familial breast cancer. First, we verified that network methods recovered more interpretable results than a standard GWAS. We addressed the heterogeneity of their solutions by studying their overlap, computing what we called the consensus. The key gene in this consensus solution was COPS5, a gene related to multiple cancer hallmarks. Another issue we observed was that network methods were unstable, selecting very different genes on different subsamples of GENESIS. Therefore, we proposed a stable consensus solution formed by the 68 genes most consistently selected across multiple subsamples. This solution was also enriched in genes known to be associated with breast cancer susceptibility (BLM, CASP8, CASP10, DNAJC1, FGFR2, MRPS30, and SLC4A7, P-value = 3 × 10-4). The most connected gene was CUL3, a regulator of several genes linked to cancer progression. Lastly, we evaluated the biases of each method and the impact of their parameters on the outcome. In general, network methods preferred highly connected genes, even after random rewirings that stripped the connections of any biological meaning. In conclusion, we present the advantages of network-guided GWAS, characterize their shortcomings, and provide strategies to address them. To compute the consensus networks, implementations of all six methods are available at https://github.com/hclimente/gwas-tools.