Experimental evolution links post-transcriptional regulation to Leishmania fitness gain.

The protozoan parasite Leishmania donovani causes fatal human visceral leishmaniasis in absence of treatment. Genome instability has been recognized as a driver in Leishmania fitness gain in response to environmental change or chemotherapy. How genome instability generates beneficial phenotypes despite potential deleterious gene dosage effects is unknown. Here we address this important open question applying experimental evolution and integrative systems approaches on parasites adapting to in vitro culture. Phenotypic analyses of parasites from early and late stages of culture adaptation revealed an important fitness tradeoff, with selection for accelerated growth in promastigote culture (fitness gain) impairing infectivity (fitness costs). Comparative genomics, transcriptomics and proteomics analyses revealed a complex regulatory network associated with parasite fitness gain, with genome instability causing highly reproducible, gene dosage-independent and -dependent changes. Reduction of flagellar transcripts and increase in coding and non-coding RNAs implicated in ribosomal biogenesis and protein translation were not correlated to dosage changes of the corresponding genes, revealing a gene dosage-independent, post-transcriptional mechanism of regulation. In contrast, abundance of gene products implicated in post-transcriptional regulation itself correlated to corresponding gene dosage changes. Thus, RNA abundance during parasite adaptation is controled by direct and indirect gene dosage changes. We correlated differential expression of small nucleolar RNAs (snoRNAs) with changes in rRNA modification, providing first evidence that Leishmania fitness gain in culture may be controlled by post-transcriptional and epitranscriptomic regulation. Our findings propose a novel model for Leishmania fitness gain in culture, where differential regulation of mRNA stability and the generation of modified ribosomes may potentially filter deleterious from beneficial gene dosage effects and provide proteomic robustness to genetically heterogenous, adapting parasite populations. This model challenges the current, genome-centric approach to Leishmania epidemiology and identifies the Leishmania transcriptome and non-coding small RNome as potential novel sources for the discovery of biomarkers that may be associated with parasite phenotypic adaptation in clinical settings.

Germline CYBB mutations that selectively affect macrophages in kindreds with X-linked predisposition to tuberculous mycobacterial disease.

Germline mutations in CYBB, the human gene encoding the gp91(phox) subunit of the phagocyte NADPH oxidase, impair the respiratory burst of all types of phagocytes and result in X-linked chronic granulomatous disease (CGD). We report here two kindreds in which otherwise healthy male adults developed X-linked recessive Mendelian susceptibility to mycobacterial disease (MSMD) syndromes. These patients had previously unknown mutations in CYBB that resulted in an impaired respiratory burst in monocyte-derived macrophages but not in monocytes or granulocytes. The macrophage-specific functional consequences of the germline mutation resulted from cell-specific impairment in the assembly of the NADPH oxidase. This 'experiment of nature' indicates that CYBB is associated with MSMD and demonstrates that the respiratory burst in human macrophages is a crucial mechanism for protective immunity to tuberculous mycobacteria.

Papuan mitochondrial genomes and the settlement of Sahul.

New Guineans represent one of the oldest locally continuous populations outside Africa, harboring among the greatest linguistic and genetic diversity on the planet. Archeological and genetic evidence suggest that their ancestors reached Sahul (present day New Guinea and Australia) by at least 55,000 years ago (kya). However, little is known about this early settlement phase or subsequent dispersal and population structuring over the subsequent period of time. Here we report 379 complete Papuan mitochondrial genomes from across Papua New Guinea, which allow us to reconstruct the phylogenetic and phylogeographic history of northern Sahul. Our results support the arrival of two groups of settlers in Sahul within the same broad time window (50-65 kya), each carrying a different set of maternal lineages and settling Northern and Southern Sahul separately. Strong geographic structure in northern Sahul remains visible today, indicating limited dispersal over time despite major climatic, cultural, and historical changes. However, following a period of isolation lasting nearly 20 ky after initial settlement, environmental changes postdating the Last Glacial Maximum stimulated diversification of mtDNA lineages and greater interactions within and beyond Northern Sahul, to Southern Sahul, Wallacea and beyond. Later, in the Holocene, populations from New Guinea, in contrast to those of Australia, participated in early interactions with incoming Asian populations from Island Southeast Asia and continuing into Oceania.

A novel rare c.-39C>T mutation in the PROS1 5'UTR causing PS deficiency by creating a new upstream translation initiation codon.

Autosomal dominant inherited Protein S deficiency (PSD) (MIM 612336) is a rare disorder caused by rare mutations, mainly located in the coding sequence of the structural PROS1 gene, and associated with an increased risk of venous thromboembolism. To identify the molecular defect underlying PSD observed in an extended French pedigree with seven PSD affected members in whom no candidate deleterious PROS1 mutation was detected by Sanger sequencing of PROS1 exons and their flanking intronic regions or via an multiplex ligation-dependent probe amplification (MLPA) approach, a whole genome sequencing strategy was adopted. This led to the identification of a never reported C to T substitution at c.-39 from the natural ATG codon of the PROS1 gene that completely segregates with PSD in the whole family. This substitution ACG→ATG creates a new start codon upstream of the main ATG. We experimentally demonstrated in HeLa cells that the variant generates a novel overlapping upstream open reading frame (uORF) and inhibits the translation of the wild-type PS. This work describes the first example of 5'UTR PROS1 mutation causing PSD through the creation of an uORF, a mutation that is not predicted to be deleterious by standard annotation softwares, and emphasizes the need for better exploration of such type of non-coding variations in clinical genomics.

Analysis with the exome array identifies multiple new independent variants in lipid loci.

It has been hypothesized that low frequency (1-5% minor allele frequency (MAF)) and rare (<1% MAF) variants with large effect sizes may contribute to the missing heritability in complex traits. Here, we report an association analysis of lipid traits (total cholesterol, LDL-cholesterol, HDL-cholesterol triglycerides) in up to 27 312 individuals with a comprehensive set of low frequency coding variants (ExomeChip), combined with conditional analysis in the known lipid loci. No new locus reached genome-wide significance. However, we found a new lead variant in 26 known lipid association regions of which 16 were >1000-fold more significant than the previous sentinel variant and not in close LD (six had MAF <5%). Furthermore, conditional analysis revealed multiple independent signals (ranging from 1 to 5) in a third of the 98 lipid loci tested, including rare variants. Addition of our novel associations resulted in between 1.5- and 2.5-fold increase in the proportion of heritability explained for the different lipid traits. Our findings suggest that rare coding variants contribute to the genetic architecture of lipid traits.

Assessment of the prognostic role of a 94-single nucleotide polymorphisms risk score in early breast cancer in the SIGNAL/PHARE prospective cohort: no correlation with clinico-pathological characteristics and outcomes.

BACKGROUND: Genome-wide association studies (GWAS) have to date identified 94 genetic variants (single nucleotide polymorphisms (SNPs)) associated with risk of developing breast cancer. A score based on the combined effect of the 94 risk alleles can be calculated to measure the global risk of breast cancer. We aimed to test the hypothesis that the 94-SNP-based risk score is associated with clinico-pathological characteristics, breast cancer subtypes and outcomes in early breast cancer. METHODS: A 94-SNP risk score was calculated in 8703 patients in the PHARE and SIGNAL prospective case cohorts. This score is the total number of inherited risk alleles based on 94 selected SNPs. Clinical data and outcomes were prospectively registered. Genotyping was obtained from a GWAS. RESULTS: The median 94-SNP risk score in 8703 patients with early breast cancer was 77.5 (range: 58.1-97.6). The risk score was not associated with usual prognostic and predictive factors (age; tumor, node, metastasis (TNM) status; Scarff-Bloom-Richardson grade; inflammatory features; estrogen receptor status; progesterone receptor status; human epidermal growth factor receptor 2 (HER2) status) and did not correlate with breast cancer subtypes. The 94-SNP risk score did not predict outcomes represented by overall survival or disease-free survival. CONCLUSIONS: In a prospective case cohort of 8703 patients, a risk score based on 94 SNPs was not associated with breast cancer characteristics, cancer subtypes, or patients' outcomes. If we hypothesize that prognosis and subtypes of breast cancer are determined by constitutional genetic factors, our results suggest that a score based on breast cancer risk-associated SNPs is not associated with prognosis. TRIAL REGISTRATION: PHARE cohort: NCT00381901 , Sept. 26, 2006 - SIGNAL cohort: INCa RECF1098, Jan. 28, 2009.

Major Loci on Chromosomes 8q and 3q Control Interferon γ Production Triggered by Bacillus Calmette-Guerin and 6-kDa Early Secretory Antigen Target, Respectively, in Various Populations.

BACKGROUND: Interferon γ (IFN-γ) release assays (IGRAs) provide an in vitro measurement of antimycobacterial immunity that is widely used as a test for Mycobacterium tuberculosis infection. IGRA outcomes are highly heritable in various populations, but the nature of the involved genetic factors remains unknown. METHODS: We conducted a genome-wide linkage analysis of IGRA phenotypes in families from a tuberculosis household contact study in France and a replication study in families from South Africa to confirm the loci identified. RESULTS: We identified a major locus on chromosome 8q controlling IFN-γ production in response to stimulation with live bacillus Calmette-Guerin (BCG; LOD score, 3.81; P = 1.40 × 10(-5)). We also detected a second locus, on chromosome 3q, that controlled IFN-γ levels in response to stimulation with 6-kDa early secretory antigen target, when accounting for the IFN-γ production shared with that induced by BCG (LOD score, 3.72; P = 1.8 × 10(-5)). Both loci were replicated in South African families, where tuberculosis is hyperendemic. These loci differ from those previously identified as controlling the response to the tuberculin skin test (TST1 and TST2) and the production of TNF-α (TNF1). CONCLUSIONS: The identification of 2 new linkage signals in populations of various ethnic origins living in different M. tuberculosis exposure settings provides new clues about the genetic control of human antimycobacterial immunity.

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.

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.

Accelerated genome sequencing with controlled costs for infants in intensive care units: a feasibility study in a French hospital network.

Obtaining a rapid etiological diagnosis for infants with early-onset rare diseases remains a major challenge. These diseases often have a severe presentation and unknown prognosis, and the genetic causes are very heterogeneous. In a French hospital network, we assessed the feasibility of performing accelerated trio-genome sequencing (GS) with limited additional costs by integrating urgent requests into the routine workflow. In addition to evaluating our capacity for such an approach, this prospective multicentre pilot study was designed to identify pitfalls encountered during its implementation. Over 14 months, we included newborns and infants hospitalized in neonatal or paediatric intensive care units with probable genetic disease and in urgent need for etiological diagnosis to guide medical care. The duration of each step and the pitfalls were recorded. We analysed any deviation from the planned schedule and identified obstacles. Trio-GS was performed for 37 individuals, leading to a molecular diagnosis in 18/37 (49%), and 21/37 (57%) after reanalysis. Corrective measures and protocol adaptations resulted in a median duration of 42 days from blood sampling to report. Accelerated trio-GS is undeniably valuable for individuals in an urgent care context. Such a circuit should coexist with a rapid or ultra-rapid circuit, which, although more expensive, can be used in particularly urgent cases. The drop in GS costs should result in its generalized use for diagnostic purposes and lead to a reduction of the costs of rapid GS.

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.

Mutations and variants of ONECUT1 in diabetes.

Genes involved in distinct diabetes types suggest shared disease mechanisms. Here we show that One Cut Homeobox 1 (ONECUT1) mutations cause monogenic recessive syndromic diabetes in two unrelated patients, characterized by intrauterine growth retardation, pancreas hypoplasia and gallbladder agenesis/hypoplasia, and early-onset diabetes in heterozygous relatives. Heterozygous carriers of rare coding variants of ONECUT1 define a distinctive subgroup of diabetic patients with early-onset, nonautoimmune diabetes, who respond well to diabetes treatment. In addition, common regulatory ONECUT1 variants are associated with multifactorial type 2 diabetes. Directed differentiation of human pluripotent stem cells revealed that loss of ONECUT1 impairs pancreatic progenitor formation and a subsequent endocrine program. Loss of ONECUT1 altered transcription factor binding and enhancer activity and NKX2.2/NKX6.1 expression in pancreatic progenitor cells. Collectively, we demonstrate that ONECUT1 controls a transcriptional and epigenetic machinery regulating endocrine development, involved in a spectrum of diabetes, encompassing monogenic (recessive and dominant) as well as multifactorial inheritance. Our findings highlight the broad contribution of ONECUT1 in diabetes pathogenesis, marking an important step toward precision diabetes medicine.

Prediction of Breast Cancer Treatment-Induced Fatigue by Machine Learning Using Genome-Wide Association Data.

BACKGROUND: We aimed at predicting fatigue after breast cancer treatment using machine learning on clinical covariates and germline genome-wide data. METHODS: We accessed germline genome-wide data of 2799 early-stage breast cancer patients from the Cancer Toxicity study (NCT01993498). The primary endpoint was defined as scoring zero at diagnosis and higher than quartile 3 at 1 year after primary treatment completion on European Organization for Research and Treatment of Cancer quality-of-life questionnaires for Overall Fatigue and on the multidimensional questionnaire for Physical, Emotional, and Cognitive fatigue. First, we tested univariate associations of each endpoint with clinical variables and genome-wide variants. Then, using preselected clinical (false discovery rate < 0.05) and genomic (P < .001) variables, a multivariable preconditioned random-forest regression model was built and validated on a hold-out subset to predict fatigue. Gene set enrichment analysis identified key biological correlates (MetaCore). All statistical tests were 2-sided. RESULTS: Statistically significant clinical associations were found only with Emotional and Cognitive Fatigue, including receipt of chemotherapy, anxiety, and pain. Some single nucleotide polymorphisms had some degree of association (P < .001) with the different fatigue endpoints, although there were no genome-wide statistically significant (P < 5.00 × 10-8) associations. Only for Cognitive Fatigue, the predictive ability of the genomic multivariable model was statistically significantly better than random (area under the curve = 0.59, P = .01) and marginally improved with clinical variables (area under the curve = 0.60, P = .005). Single nucleotide polymorphisms found to be associated (P < .001) with Cognitive Fatigue belonged to genes linked to inflammation (false discovery rate adjusted P = .03), cognitive disorders (P = 1.51 × 10-12), and synaptic transmission (P = 6.28 × 10-8). CONCLUSIONS: Genomic analyses in this large cohort of breast cancer survivors suggest a possible genetic role for severe Cognitive Fatigue that warrants further exploration.

Evidence of Austronesian Genetic Lineages in East Africa and South Arabia: Complex Dispersal from Madagascar and Southeast Asia.

The Austronesian dispersal across the Indonesian Ocean to Madagascar and the Comoros has been well documented, but in an unexplained anomaly, few to no traces have been found of the Austronesian expansion in East Africa or the Arabian Peninsula. To revisit this peculiarity, we surveyed the Western Indian Ocean rim populations to identify potential Austronesian genetic ancestry. We generated full mitochondrial DNA genomes and genome-wide genotyping data for these individuals and compared them with the Banjar, the Indonesian source population of the westward Austronesian dispersal. We find strong support for Asian genetic contributions to maternal lineages and autosomal variation in modern day Somalia and Yemen. Surprisingly, this input reveals two apparently different geographic origins and timings of admixture for the Austronesian contact; one at a very early phase (likely associated with the early Austronesian dispersals), and a later movement dating to the end of nineteenth century. These Austronesian gene flows come, respectively, from Madagascar and directly from an unidentified location in Island Southeast Asia. This result reveals a far more complex dynamic of Austronesian dispersals through the Western Indian Ocean than has previously been understood and suggests that Austronesian movements within the Indian Ocean may have been part of a lengthy process, probably continuing well into the modern era.

Genome-Wide Methylation Analysis Identifies Specific Epigenetic Marks In Severely Obese Children.

Obesity is a heterogeneous disease with many different subtypes. Epigenetics could contribute to these differences. The aim of this study was to investigate genome-wide DNA methylation searching for methylation marks associated with obesity in children and adolescents. We studied DNA methylation profiles in whole blood cells from 40 obese children and controls using Illumina Infinium HumanMethylation450 BeadChips. After correction for cell heterogeneity and multiple tests, we found that compared to lean controls, 31 CpGs are differentially methylated in obese patients. A greatest proportion of these CpGs is hypermethylated in obesity and located in CpG shores regions. We next focused on severely obese children and identified 151 differentially methylated CpGs among which 10 with a difference in methylation greater than 10%. The top pathways enriched among the identified CpGs included the "IRS1 target genes" and several pathways in cancer diseases. This study represents the first effort to search for differences in methylation in obesity and severe obesity, which may help understanding these different forms of obesity and their complications.

A GWAS in uveal melanoma identifies risk polymorphisms in the CLPTM1L locus.

Uveal melanoma, a rare malignant tumor of the eye, is predominantly observed in populations of European ancestry. A genome-wide association study of 259 uveal melanoma patients compared to 401 controls all of European ancestry revealed a candidate locus at chromosome 5p15.33 (region rs421284: OR = 1.7, CI 1.43-2.05). This locus was replicated in an independent set of 276 cases and 184 controls. In addition, risk variants from this region were positively associated with higher expression of CLPTM1L. In conclusion, the CLPTM1L region contains risk alleles for uveal melanoma susceptibility, suggesting that CLPTM1L could play a role in uveal melanoma oncogenesis.

dUTPase (DUT) Is Mutated in a Novel Monogenic Syndrome With Diabetes and Bone Marrow Failure.

We describe a new syndrome characterized by early-onset diabetes associated with bone marrow failure, affecting mostly the erythrocytic lineage. Using whole-exome sequencing in a remotely consanguineous patient from a family with two affected siblings, we identified a single homozygous missense mutation (chr15.hg19:g.48,626,619A>G) located in the dUTPase (DUT) gene (National Center for Biotechnology Information Gene ID 1854), affecting both the mitochondrial (DUT-M p.Y142C) and the nuclear (DUT-N p.Y54C) isoforms. We found the same homozygous mutation in an unrelated consanguineous patient with diabetes and bone marrow aplasia from a family with two affected siblings, whereas none of the >60,000 subjects from the Exome Aggregation Consortium (ExAC) was homozygous for this mutation. This replicated observation probability was highly significant, thus confirming the role of this DUT mutation in this syndrome. DUT is a key enzyme for maintaining DNA integrity by preventing misincorporation of uracil into DNA, which results in DNA toxicity and cell death. We showed that DUT silencing in human and rat pancreatic ß-cells results in apoptosis via the intrinsic cell death pathway. Our findings support the importance of tight control of DNA metabolism for ß-cell integrity and warrant close metabolic monitoring of patients treated by drugs affecting dUTP balance.

Constitutional variants are not associated with HER2-positive breast cancer: results from the SIGNAL/PHARE clinical cohort.

Human epidermal growth factor receptor 2-positive breast cancer is a subtype of interest regarding its outcome and the impressive impact of human epidermal growth factor receptor 2 targeted therapy. Constitutional variants may be involved in the aetiology of human epidermal growth factor receptor 2-positive breast cancer, and we propose a case-case study to test the hypothesis that single nucleotide polymorphisms may be associated with human epidermal growth factor receptor 2 status. A Genome-Wide Association Study was used in a cohort of 9836 patients from the SIGNAL/PHARE study (NCT00381901-RECF1098). The main goal was to identify variants specifically related to human epidermal growth factor receptor 2-positive breast cancer. A two-staged genotyping strategy was carried out to cover as large a proportion of the genome as possible. All subjects were genotyped using the Illumina HumanCore Exome chip set. Principal Components Analysis and k-means were then used to characterize the ancestry of the participants. A random sample of subjects from the main "European" cluster was genotyped with the Omni5 chip set. These data were then used to impute missing genotypes from the remaining subjects genotyped only using the HumanCore Exome array. From the 9836 patients, a total of 8703 cases including 3230 patients with human epidermal growth factor receptor 2-positive breast cancer were analyzed. Despite having 80% power to detect an odds ratio of 1.23 in this population, no variant achieved genome-wide significance for association with the occurrence of human epidermal growth factor receptor 2-positive breast cancer vs. any other subtype of breast tumour. Our study was unable to identify constitutional polymorphisms that are strongly associated with human epidermal growth factor receptor 2-positive status among breast cancer patients.

Genome-wide association study identifies multiple risk loci for renal cell carcinoma.

Previous genome-wide association studies (GWAS) have identified six risk loci for renal cell carcinoma (RCC). We conducted a meta-analysis of two new scans of 5,198 cases and 7,331 controls together with four existing scans, totalling 10,784 cases and 20,406 controls of European ancestry. Twenty-four loci were tested in an additional 3,182 cases and 6,301 controls. We confirm the six known RCC risk loci and identify seven new loci at 1p32.3 (rs4381241, P=3.1 × 10-10), 3p22.1 (rs67311347, P=2.5 × 10-8), 3q26.2 (rs10936602, P=8.8 × 10-9), 8p21.3 (rs2241261, P=5.8 × 10-9), 10q24.33-q25.1 (rs11813268, P=3.9 × 10-8), 11q22.3 (rs74911261, P=2.1 × 10-10) and 14q24.2 (rs4903064, P=2.2 × 10-24). Expression quantitative trait analyses suggest plausible candidate genes at these regions that may contribute to RCC susceptibility.

GWAS in the SIGNAL/PHARE clinical cohort restricts the association between the FGFR2 locus and estrogen receptor status to HER2-negative breast cancer patients.

Genetic polymorphisms are associated with breast cancer risk. Clinical and epidemiological observations suggest that clinical characteristics of breast cancer, such as estrogen receptor or HER2 status, are also influenced by hereditary factors. To identify genetic variants associated with pathological characteristics of breast cancer patients, a Genome Wide Association Study was performed in a cohort of 9365 women from the French nationwide SIGNAL/PHARE studies (NCT00381901/RECF1098). Strong association between the FGFR2 locus and ER status of breast cancer patients was observed (ER-positive n=6211, ER-negative n=2516; rs3135718 OR=1.34 p=5.46×10-12). This association was limited to patients with HER2-negative tumors (ER-positive n=4267, ER-negative n=1185; rs3135724 OR=1.85 p=1.16×10-11). The FGFR2 locus is known to be associated with breast cancer risk. This study provides sound evidence for an association between variants in the FGFR2 locus and ER status among breast cancer patients, particularly among patients with HER2-negative disease. This refinement of the association between FGFR2 variants and ER-status to HER2-negative disease provides novel insight to potential biological and clinical influence of genetic polymorphisms on breast tumors.