Microsatellite instability at U2AF-binding polypyrimidic tract sites perturbs alternative splicing during colorectal cancer initiation.

BACKGROUND: Microsatellite instability (MSI) due to mismatch repair deficiency (dMMR) is common in colorectal cancer (CRC). These cancers are associated with somatic coding events, but the noncoding pathophysiological impact of this genomic instability is yet poorly understood. Here, we perform an analysis of coding and noncoding MSI events at the different steps of colorectal tumorigenesis using whole exome sequencing and search for associated splicing events via RNA sequencing at the bulk-tumor and single-cell levels. RESULTS: Our results demonstrate that MSI leads to hundreds of noncoding DNA mutations, notably at polypyrimidine U2AF RNA-binding sites which are endowed with cis-activity in splicing, while higher frequency of exon skipping events are observed in the mRNAs of MSI compared to non-MSI CRC. At the DNA level, these noncoding MSI mutations occur very early prior to cell transformation in the dMMR colonic crypt, accounting for only a fraction of the exon skipping in MSI CRC. At the RNA level, the aberrant exon skipping signature is likely to impair colonic cell differentiation in MSI CRC affecting the expression of alternative exons encoding protein isoforms governing cell fate, while also targeting constitutive exons, making dMMR cells immunogenic in early stage before the onset of coding mutations. This signature is characterized by its similarity to the oncogenic U2AF1-S34F splicing mutation observed in several other non-MSI cancer. CONCLUSIONS: Overall, these findings provide evidence that a very early RNA splicing signature partly driven by MSI impairs cell differentiation and promotes MSI CRC initiation, far before coding mutations which accumulate later during MSI tumorigenesis.

Human genetic structure in Northwest France provides new insights into West European historical demography.

The demographical history of France remains largely understudied despite its central role toward understanding modern population structure across Western Europe. Here, by exploring publicly available Europe-wide genotype datasets together with the genomes of 3234 present-day and six newly sequenced medieval individuals from Northern France, we found extensive fine-scale population structure across Brittany and the downstream Loire basin and increased population differentiation between the northern and southern sides of the river Loire, associated with higher proportions of steppe vs. Neolithic-related ancestry. We also found increased allele sharing between individuals from Western Brittany and those associated with the Bell Beaker complex. Our results emphasise the need for investigating local populations to better understand the distribution of rare (putatively deleterious) variants across space and the importance of common genetic legacy in understanding the sharing of disease-related alleles between Brittany and people from western Britain and Ireland.

Exome sequencing in undiagnosed congenital myopathy reveals new genes and refines genes-phenotypes correlations.

BACKGROUND: Congenital myopathies are severe genetic diseases with a strong impact on patient autonomy and often on survival. A large number of patients do not have a genetic diagnosis, precluding genetic counseling and appropriate clinical management. Our objective was to find novel pathogenic variants and genes associated with congenital myopathies and to decrease diagnostic odysseys and dead-end. METHODS: To identify pathogenic variants and genes implicated in congenital myopathies, we established and conducted the MYOCAPTURE project from 2009 to 2018 to perform exome sequencing in a large cohort of 310 families partially excluded for the main known genes. RESULTS: Pathogenic variants were identified in 156 families (50%), among which 123 families (40%) had a conclusive diagnosis. Only 44 (36%) of the resolved cases were linked to a known myopathy gene with the corresponding phenotype, while 55 (44%) were linked to pathogenic variants in a known myopathy gene with atypical signs, highlighting that most genetic diagnosis could not be anticipated based on clinical-histological assessments in this cohort. An important phenotypic and genetic heterogeneity was observed for the different genes and for the different congenital myopathy subtypes, respectively. In addition, we identified 14 new myopathy genes not previously associated with muscle diseases (20% of all diagnosed cases) that we previously reported in the literature, revealing novel pathomechanisms and potential therapeutic targets. CONCLUSIONS: Overall, this approach illustrates the importance of massive parallel gene sequencing as a comprehensive tool for establishing a molecular diagnosis for families with congenital myopathies. It also emphasizes the contribution of clinical data, histological findings on muscle biopsies, and the availability of DNA samples from additional family members to the diagnostic success rate. This study facilitated and accelerated the genetic diagnosis of congenital myopathies, improved health care for several patients, and opened novel perspectives for either repurposing of existing molecules or the development of novel treatments.

Measuring the Efficiency of Purging by non-random Mating in Human Populations.

Human populations harbor a high concentration of deleterious genetic variants. Here, we tested the hypothesis that non-random mating practices affect the distribution of these variants, through exposure in the homozygous state, leading to their purging from the population gene pool. To do so, we produced whole-genome sequencing data for two pairs of Asian populations exhibiting different alliance rules and rates of inbreeding, but with similar effective population sizes. The results show that populations with higher rates of inbred matings do not purge deleterious variants more efficiently. Purging therefore has a low efficiency in human populations, and different mating practices lead to a similar mutational load.

Multi-ancestry polygenic risk scores for venous thromboembolism.

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality, with large disparities in incidence rates between Black and White Americans. Polygenic risk scores (PRSs) limited to variants discovered in genome-wide association studies in European-ancestry samples can identify European-ancestry individuals at high risk of VTE. However, there is limited evidence on whether high-dimensional PRS constructed using more sophisticated methods and more diverse training data can enhance the predictive ability and their utility across diverse populations. We developed PRSs for VTE using summary statistics from the International Network against Venous Thrombosis (INVENT) consortium genome-wide association studies meta-analyses of European- (71 771 cases and 1 059 740 controls) and African-ancestry samples (7482 cases and 129 975 controls). We used LDpred2 and PRS-CSx to construct ancestry-specific and multi-ancestry PRSs and evaluated their performance in an independent European- (6781 cases and 103 016 controls) and African-ancestry sample (1385 cases and 12 569 controls). Multi-ancestry PRSs with weights tuned in European-ancestry samples slightly outperformed ancestry-specific PRSs in European-ancestry test samples (e.g. the area under the receiver operating curve [AUC] was 0.609 for PRS-CSx_combinedEUR and 0.608 for PRS-CSxEUR [P = 0.00029]). Multi-ancestry PRSs with weights tuned in African-ancestry samples also outperformed ancestry-specific PRSs in African-ancestry test samples (PRS-CSxAFR: AUC = 0.58, PRS-CSx_combined AFR: AUC = 0.59), although this difference was not statistically significant (P = 0.34). The highest fifth percentile of the best-performing PRS was associated with 1.9-fold and 1.68-fold increased risk for VTE among European- and African-ancestry subjects, respectively, relative to those in the middle stratum. These findings suggest that the multi-ancestry PRS might be used to improve performance across diverse populations to identify individuals at highest risk for VTE.

Plasma levels of complement components C5 and C9 are associated with thrombin generation.

BACKGROUND: The thrombin generation assay (TGA) evaluates the potential of plasma to generate thrombin over time, providing a global picture of an individual's hemostatic balance. OBJECTIVES: This study aimed to identify novel biological determinants of thrombin generation using a multiomics approach. METHODS: Associations between TGA parameters and plasma levels of 377 antibodies targeting 236 candidate proteins for cardiovascular risk were tested using multiple linear regression analysis in 770 individuals with venous thrombosis from the Marseille Thrombosis Association (MARTHA) study. Proteins associated with at least 3 TGA parameters were selected for validation in an independent population of 536 healthy individuals (Etablissement Français du Sang Alpes-Méditerranée [EFS-AM]). Proteins with strongest associations in both groups underwent additional genetic analyses and in vitro experiments. RESULTS: Eighteen proteins were associated (P < 1.33 × 10⁻4) with at least 3 TGA parameters in MARTHA, among which 13 demonstrated a similar pattern of associations in EFS-AM. Complement proteins C5 and C9 had the strongest associations in both groups. Ex vivo supplementation of platelet-poor plasma with purified C9 protein had a significant dose-dependent effect on TGA parameters. No effect was observed with purified C5. Several single nucleotide polymorphisms associated with C5 and C9 plasma levels were identified, with the strongest association for the C5 missense variant rs17611, which was associated with a decrease in C5 levels, endogenous thrombin potential, and peak in MARTHA. No association of this variant with TGA parameters was observed in EFS-AM. CONCLUSION: This study identified complement proteins C5 and C9 as potential determinants of thrombin generation. Further studies are warranted to establish causality and elucidate the underlying mechanisms.

MYLK*FLNB and DOCK1*LAMA2 gene-gene interactions associated with rheumatoid arthritis in the focal adhesion pathway.

Rheumatoid arthritis (RA) is a chronic, systemic autoimmune disease caused by a combination of genetic and environmental factors. Rare variants with low predicted effects in genes participating in the same biological function might be involved in developing complex diseases such as RA. From whole-exome sequencing (WES) data, we identified genes containing rare non-neutral variants with complete penetrance and no phenocopy in at least one of nine French multiplex families. Further enrichment analysis highlighted focal adhesion as the most significant pathway. We then tested if interactions between the genes participating in this function would increase or decrease the risk of developing RA disease. The model-based multifactor dimensionality reduction (MB-MDR) approach was used to detect epistasis in a discovery sample (19 RA cases and 11 healthy individuals from 9 families and 98 unrelated CEU controls from the International Genome Sample Resource). We identified 9 significant interactions involving 11 genes (MYLK, FLNB, DOCK1, LAMA2, RELN, PIP5K1C, TNC, PRKCA, VEGFB, ITGB5, and FLT1). One interaction (MYLK*FLNB) increasing RA risk and one interaction decreasing RA risk (DOCK1*LAMA2) were confirmed in a replication sample (200 unrelated RA cases and 91 GBR unrelated controls). Functional and genomic data in RA samples or relevant cell types argue the key role of these genes in RA.

Genome-wide investigation of exogenous female hormones, genetic variation, and venous thromboembolism risk.

BACKGROUND: Increased risk of venous thromboembolism (VTE) is a life-threatening side effect for users of oral contraceptives (OCs) or hormone therapy (HT). OBJECTIVES: To investigate the potential for genetic predisposition to VTE in OC or HT users, we conducted a gene-by-environment case-only meta-analysis of genome-wide association studies (GWAS). METHODS: Use or nonuse of OCs (7 studies) or HT (8 studies) at the time of the VTE event was determined by pharmacy records or self-report. A synergy index (SI) was modeled for each variant in each study and submultiplicative/supramultiplicative gene-by-environment interactions were estimated. The SI parameters were first meta-analyzed across OC and HT studies and subsequently meta-analyzed to obtain an overall estimate. The primary analysis was agnostic GWAS and interrogated all imputed genotypes using a P value threshold of <5.0 × 10-8; secondary analyses were candidate-based. RESULTS: The VTE case-only OC meta-analysis included 2895 OC users and 6607 nonusers; the case-only HT meta-analysis included 2434 HT users and 12 793 nonusers. In primary GWAS meta-analyses, no variant reached genome-wide significance, but the smallest P value approached statistical significance: rs9386463 (P = 5.03 × 10-8). We tested associations for 138 candidate variants and identified 2 that exceeded statistical significance (0.05/138 = 3.62 × 10-4): F5 rs6025 (P = 1.87 × 10-5; SI, 1.29; previously observed) and F11 rs2036914 (P = 2.0 × 10-4; SI, 0.91; new observation). CONCLUSION: The candidate variant approach to identify submultiplictive/supramultiplicative associations between genetic variation and OC and HT use identified a new association with common genetic variation in F11, while the agnostic interrogations did not yield new discoveries.

Positive selection in the genomes of two Papua New Guinean populations at distinct altitude levels.

Highlanders and lowlanders of Papua New Guinea have faced distinct environmental stress, such as hypoxia and environment-specific pathogen exposure, respectively. In this study, we explored the top genomics regions and the candidate driver SNPs for selection in these two populations using newly sequenced whole-genomes of 54 highlanders and 74 lowlanders. We identified two candidate SNPs under selection - one in highlanders, associated with red blood cell traits and another in lowlanders, which is associated with white blood cell count - both potentially influencing the heart rate of Papua New Guineans in opposite directions. We also observed four candidate driver SNPs that exhibit linkage disequilibrium with an introgressed haplotype, highlighting the need to explore the possibility of adaptive introgression within these populations. This study reveals that the signatures of positive selection in highlanders and lowlanders of Papua New Guinea align closely with the challenges they face, which are specific to their environments.

TAD boundary deletion causes PITX2-related cardiac electrical and structural defects.

While 3D chromatin organization in topologically associating domains (TADs) and loops mediating regulatory element-promoter interactions is crucial for tissue-specific gene regulation, the extent of their involvement in human Mendelian disease is largely unknown. Here, we identify 7 families presenting a new cardiac entity associated with a heterozygous deletion of 2 CTCF binding sites on 4q25, inducing TAD fusion and chromatin conformation remodeling. The CTCF binding sites are located in a gene desert at 1 Mb from the Paired-like homeodomain transcription factor 2 gene (PITX2). By introducing the ortholog of the human deletion in the mouse genome, we recapitulate the patient phenotype and characterize an opposite dysregulation of PITX2 expression in the sinoatrial node (ectopic activation) and ventricle (reduction), respectively. Chromatin conformation assay performed in human induced pluripotent stem cell-derived cardiomyocytes harboring the minimal deletion identified in family#1 reveals a conformation remodeling and fusion of TADs. We conclude that TAD remodeling mediated by deletion of CTCF binding sites causes a new autosomal dominant Mendelian cardiac disorder.

A comparative genomics approach reveals a local genetic signature of Leishmania tropica in Morocco.

In Morocco, cutaneous leishmaniasis (CL) caused by Leishmania (L.) tropica is an important health problem. Despite the high incidence of CL in the country, the genomic heterogeneity of these parasites is still incompletely understood. In this study, we sequenced the genomes of 14 Moroccan isolates of L. tropica collected from confirmed cases of CL to investigate their genomic heterogeneity. Comparative genomics analyses were conducted by applying the recently established Genome Instability Pipeline (GIP), which allowed us to conduct phylogenomic and principal components analyses (PCA), and to assess genomic variations at the levels of the karyotype, gene copy number, single nucleotide polymorphisms (SNPs) and small insertions/deletions (INDELs) variants. Read-depth analyses revealed a mostly disomic karyotype, with the exception of the stable tetrasomy of chromosome 31. In contrast, we identified important gene copy number variations across all isolates, which affect known virulence genes and thus were probably selected in the field. SNP-based cluster analysis of the 14 isolates revealed a core group of 12 strains that formed a tight cluster and shared 45.1 % (87 751) of SNPs, as well as two strains (M3015, Ltr_16) that clustered separately from each other and the core group, suggesting the circulation of genetically highly diverse strains in Morocco. Phylogenetic analysis, which compared our 14 L. tropica isolates against 40 published genomes of L. tropica from a diverse array of locations, confirmed the genetic difference of our Moroccan isolates from all other isolates examined. In conclusion, our results indicate potential regional variations in SNP profiles that may differentiate Moroccan L. tropica from other L. tropica strains circulating in endemic countries in the Middle East. Our report paves the way for future research with a larger number of strains that will allow correlation of diverse phenotypes (resistance to treatments, virulence) and origins (geography, host species, year of isolation) to defined genomic signals such as gene copy number variations or SNP profiles that may represent interesting biomarker candidates.

Upstream open reading frame-introducing variants in patients with primary familial brain calcification.

More than 50% of patients with primary familial brain calcification (PFBC), a rare neurological disorder, remain genetically unexplained. While some causative genes are yet to be identified, variants in non-coding regions of known genes may represent a source of missed diagnoses. We hypothesized that 5'-Untranslated Region (UTR) variants introducing an AUG codon may initiate mRNA translation and result in a loss of function in some of the PFBC genes. After reannotation of exome sequencing data of 113 unrelated PFBC probands, we identified two upstream AUG-introducing variants in the 5'UTR of PDGFB. One, NM_002608.4:c.-373C>G, segregated with PFBC in the family. It was predicted to create an upstream open reading frame (ORF). The other one, NM_002608.4:c.-318C>T, was found in a simplex case. It was predicted to result in an ORF overlapping the natural ORF with a frameshift. In a GFP reporter assay, both variants were associated with a dramatic decrease in GFP levels, and, after restoring the reading frame with the GFP sequence, the c.-318C>T variant was associated with a strong initiation of translation as measured by western blotting. Overall, we found upstream AUG-introducing variants in the 5'UTR of PDGFB in 2/113 (1.7%) undiagnosed PFBC cases. Such variants thus represent a source of putative pathogenic variants.

De Novo Variants Found in Three Distinct Schizophrenia Populations Hit a Common Core Gene Network Related to Microtubule and Actin Cytoskeleton Gene Ontology Classes.

Schizophrenia (SZ) is a heterogeneous and debilitating psychiatric disorder with a strong genetic component. To elucidate functional networks perturbed in schizophrenia, we analysed a large dataset of whole-genome studies that identified SNVs, CNVs, and a multi-stage schizophrenia genome-wide association study. Our analysis identified three subclusters that are interrelated and with small overlaps: GO:0007017~Microtubule-Based Process, GO:00015629~Actin Cytoskeleton, and GO:0007268~SynapticTransmission. We next analysed three distinct trio cohorts of 75 SZ Algerian, 45 SZ French, and 61 SZ Japanese patients. We performed Illumina HiSeq whole-exome sequencing and identified de novo mutations using a Bayesian approach. We validated 88 de novo mutations by Sanger sequencing: 35 in French, 21 in Algerian, and 32 in Japanese SZ patients. These 88 de novo mutations exhibited an enrichment in genes encoding proteins related to GO:0051015~actin filament binding (p = 0.0011) using David, and enrichments in GO: 0003774~transport (p = 0.019) and GO:0003729~mRNA binding (p = 0.010) using Amigo. One of these de novo variant was found in CORO1C coding sequence. We studied Coro1c haploinsufficiency in a Coro1c+/- mouse and found defects in the corpus callosum. These results could motivate future studies of the mechanisms surrounding genes encoding proteins involved in transport and the cytoskeleton, with the goal of developing therapeutic intervention strategies for a subset of SZ cases.

Multi-ancestry polygenic risk scores for venous thromboembolism.

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality, with large disparities in incidence rates between Black and White Americans. Polygenic risk scores (PRSs) limited to variants discovered in genome-wide association studies in European-ancestry samples can identify European-ancestry individuals at high risk of VTE. However, there is limited evidence on whether high-dimensional PRS constructed using more sophisticated methods and more diverse training data can enhance the predictive ability and their utility across diverse populations. We developed PRSs for VTE using summary statistics from the International Network against Venous Thrombosis (INVENT) consortium GWAS meta-analyses of European- (71,771 cases and 1,059,740 controls) and African-ancestry samples (7,482 cases and 129,975 controls). We used LDpred2 and PRSCSx to construct ancestry-specific and multi-ancestry PRSs and evaluated their performance in an independent European- (6,261 cases and 88,238 controls) and African-ancestry sample (1,385 cases and 12,569 controls). Multi-ancestry PRSs with weights tuned in European- and African-ancestry samples, respectively, outperformed ancestry-specific PRSs in European- (PRSCSXEUR: AUC=0.61 (0.60, 0.61), PRSCSX_combinedEUR: AUC=0.61 (0.60, 0.62)) and African-ancestry test samples (PRSCSXAFR: AUC=0.58 (0.57, 0.6), PRSCSX_combined AFR: AUC=0.59 (0.57, 0.60)). The highest fifth percentile of the best-performing PRS was associated with 1.9-fold and 1.68-fold increased risk for VTE among European- and African-ancestry subjects, respectively, relative to those in the middle stratum. These findings suggest that the multi-ancestry PRS may be used to identify individuals at highest risk for VTE and provide guidance for the most effective treatment strategy across diverse populations.

Assessment of Mendelian and risk-factor genes in Alzheimer disease: A prospective nationwide clinical utility study and recommendations for genetic screening.

PURPOSE: To assess the likely pathogenic/pathogenic (LP/P) variants rates in Mendelian dementia genes and the moderate-to-strong risk factors rates in patients with Alzheimer disease (AD). METHODS: We included 700 patients in a prospective study and performed exome sequencing. A panel of 28 Mendelian and 6 risk-factor genes was interpreted and returned to patients. We built a framework for risk variant interpretation and risk gradation and assessed the detection rates among early-onset AD (EOAD, age of onset (AOO) ≤65 years, n = 608) depending on AOO and pedigree structure and late-onset AD (66 < AOO < 75, n = 92). RESULTS: Twenty-one patients carried a LP/P variant in a Mendelian gene (all with EOAD, 3.4%), 20 of 21 affected APP, PSEN1, or PSEN2. LP/P variant detection rates in EOAD ranged from 1.7% to 11.6% based on AOO and pedigree structure. Risk factors were found in 69.5% of the remaining 679 patients, including 83 (12.2%) being heterozygotes for rare risk variants, in decreasing order of frequency, in TREM2, ABCA7, ATP8B4, SORL1, and ABCA1, including 5 heterozygotes for multiple rare risk variants, suggesting non-monogenic inheritance, even in some autosomal-dominant-like pedigrees. CONCLUSION: We suggest that genetic screening should be proposed to all EOAD patients and should no longer be prioritized based on pedigree structure.

Episignatures in practice: independent evaluation of published episignatures for the molecular diagnostics of ten neurodevelopmental disorders.

Variants of uncertain significance (VUS) are a significant issue for the molecular diagnosis of rare diseases. The publication of episignatures as effective biomarkers of certain Mendelian neurodevelopmental disorders has raised hopes to help classify VUS. However, prediction abilities of most published episignatures have not been independently investigated yet, which is a prerequisite for an informed and rigorous use in a diagnostic setting. We generated DNA methylation data from 101 carriers of (likely) pathogenic variants in ten different genes, 57 VUS carriers, and 25 healthy controls. Combining published episignature information and new validation data with a k-nearest-neighbour classifier within a leave-one-out scheme, we provide unbiased specificity and sensitivity estimates for each of the signatures. Our procedure reached 100% specificity, but the sensitivities unexpectedly spanned a very large spectrum. While ATRX, DNMT3A, KMT2D, and NSD1 signatures displayed a 100% sensitivity, CREBBP-RSTS and one of the CHD8 signatures reached <40% sensitivity on our dataset. Remaining Cornelia de Lange syndrome, KMT2A, KDM5C and CHD7 signatures reached 70-100% sensitivity at best with unstable performances, suffering from heterogeneous methylation profiles among cases and rare discordant samples. Our results call for cautiousness and demonstrate that episignatures do not perform equally well. Some signatures are ready for confident use in a diagnostic setting. Yet, it is imperative to characterise the actual validity perimeter and interpretation of each episignature with the help of larger validation sample sizes and in a broader set of episignatures.

Inverse Modulation of Aurora Kinase A and Topoisomerase IIα in Normal and Tumor Breast Cells upon Knockdown of Mitochondrial ASncmtRNA.

Breast cancer is currently the most diagnosed form of cancer and the leading cause of death by cancer among females worldwide. We described the family of long non-coding mitochondrial RNAs (ncmtRNAs), comprised of sense (SncmtRNA) and antisense (ASncmtRNA) members. Knockdown of ASncmtRNAs using antisense oligonucleotides (ASOs) induces proliferative arrest and apoptotic death of tumor cells, but not normal cells, from various tissue origins. In order to study the mechanisms underlying this selectivity, in this study we performed RNAseq in MDA-MB-231 breast cancer cells transfected with ASncmtRNA-specific ASO or control-ASO, or left untransfected. Bioinformatic analysis yielded several differentially expressed cell-cycle-related genes, from which we selected Aurora kinase A (AURKA) and topoisomerase IIα (TOP2A) for RT-qPCR and western blot validation in MDA-MB-231 and MCF7 breast cancer cells, as well as normal breast epithelial cells (HMEC). We observed no clear differences regarding mRNA levels but both proteins were downregulated in tumor cells and upregulated in normal cells. Since these proteins play a role in genomic integrity, this inverse effect of ASncmtRNA knockdown could account for tumor cell downfall whilst protecting normal cells, suggesting this approach could be used for genomic protection under cancer treatment regimens or other scenarios.

A Case-Only Genome-Wide Interaction Study of Smoking and Bladder Cancer Risk: Results from the COBLAnCE Cohort.

Bladder cancer (BC) is the 6th most common cancer worldwide, with tobacco smoking considered as its main risk factor. Accumulating evidence has found associations between genetic variants and the risk of BC. Candidate gene-environment interaction studies have suggested interactions between cigarette smoking and NAT2/GSTM1 gene variants. Our objective was to perform a genome-wide association case-only study using the French national prospective COBLAnCE cohort (COhort to study BLAdder CancEr), focusing on smoking behavior. The COBLAnCE cohort comprises 1800 BC patients enrolled between 2012 and 2018. Peripheral blood samples collected at enrolment were genotyped using the Illumina Global Screening Array with a Multi-Disease drop-in panel. Genotyping data (9,719,614 single nucleotide polymorphisms (SNP)) of 1674, 1283, and 1342 patients were analyzed for smoking status, average tobacco consumption, and age at smoking initiation, respectively. A genome-wide association study (GWAS) was conducted adjusting for gender, age, and genetic principal components. The results suggest new candidate loci (4q22.1, 12p13.1, 16p13.3) interacting with smoking behavior for the risk of BC. Our results need to be validated in other case-control or cohort studies.