Genome-wide SNP detection, validation, and development of an 8K SNP array for apple

As high-throughput genetic marker screening systems are essential for a range of genetics studies and plant breeding applications, the International RosBREED SNP Consortium (IRSC) has utilized the Illumina Infinium® II system to develop a medium- to high-throughput SNP screening tool for genome-wide evaluation of allelic variation in apple (Malus×domestica) breeding germplasm. For genome-wide SNP discovery, 27 apple cultivars were chosen to represent worldwide breeding germplasm and re-sequenced at low coverage with the Illumina Genome Analyzer II. Following alignment of these sequences to the whole genome sequence of 'Golden Delicious', SNPs were identified using SoapSNP. A total of 2,113,120 SNPs were detected, corresponding to one SNP to every 288 bp of the genome. The Illumina GoldenGate® assay was then used to validate a subset of 144 SNPs with a range of characteristics, using a set of 160 apple accessions. This validation assay enabled fine-tuning of the final subset of SNPs for the Illumina Infinium® II system. The set of stringent filtering criteria developed allowed choice of a set of SNPs that not only exhibited an even distribution across the apple genome and a range of minor allele frequencies to ensure utility across germplasm, but also were located in putative exonic regions to maximize genotyping success rate. A total of 7867 apple SNPs was established for the IRSC apple 8K SNP array v1, of which 5554 were polymorphic after evaluation in segregating families and a germplasm collection. This publicly available genomics resource will provide an unprecedented resolution of SNP haplotypes, which will enable marker-locus-trait association discovery, description of the genetic architecture of quantitative traits, investigation of genetic variation (neutral and functional), and genomic selection in apple.

Dopamine transporter genotype predicts behavioural and neural measures of response inhibition

The ability to inhibit unwanted actions is a heritable executive function that may confer risk to disorders such as attention deficit hyperactivity disorder (ADHD). Converging evidence from pharmacology and cognitive neuroscience suggests that response inhibition is instantiated within frontostriatal circuits of the brain with patterns of activity that are modulated by the catecholamines dopamine and noradrenaline. A total of 405 healthy adult participants performed the stop-signal task, a paradigmatic measure of response inhibition that yields an index of the latency of inhibition, termed the stop-signal reaction time (SSRT). Using this phenotype, we tested for genetic association, performing high-density single-nucleotide polymorphism mapping across the full range of autosomal catecholamine genes. Fifty participants also underwent functional magnetic resonance imaging to establish the impact of associated alleles on brain and behaviour. Allelic variation in polymorphisms of the dopamine transporter gene (SLC6A3: rs37020; rs460000) predicted individual differences in SSRT, after corrections for multiple comparisons. Furthermore, activity in frontal regions (anterior frontal, superior frontal and superior medial gyri) and caudate varied additively with the T-allele of rs37020. The influence of genetic variation in SLC6A3 on the development of frontostriatal inhibition networks may represent a key risk mechanism for disorders of behavioural inhibition.

Single nucleotide polymorphisms and mRNA expression of VEGF-A in papillary thyroid carcinoma: Potential markers for aggressive phenotypes

BACKGROUND AND OBJECTIVES Polymorphisms of the VEGF gene are known to affect the biological behaviour of cancers but have seldom been studied in thyroid cancer. The aim of the current study is to evaluate the prevalence and relevance of VEGF-A polymorphisms and mRNA expression in papillary thyroid carcinoma (PTC). MATERIALS AND METHODS Genomic DNA and total RNA were isolated from paraffin-embedded tissue from 91 PTC (51 conventional PTC and 40 follicular variant) and 78 control thyroid tissues. Three DNA polymorphisms (+936C > T, +405C > G and -141A > C) in the 3' and 5' untranslated region (3'-UTR, 5'-UTR) of VEGF-A were studied using PCR and RFLP. Also, the mRNA expression of VEGF-A in these tissues was studied by real-time PCR. RESULTS Distribution of polymorphisms in the 5'-UTR (VEGF-VEGF -141A > C and +405C > G) and 3'-UTR (VEGF +936C > T) were all significantly different in PTC and benign thyroid tissue (p = 0.0001, 0.001 and 0.028 respectively). The VEGF -141 C allele was more common in PTC with lymph node metastases (p = 0.026). VEGF + 405 Galleles andVEGF +936 CC genotype were more common in PTC of advanced pathological staging (p = 0.018 and 0.017 respectively). Also, increased expression of VEGF-A mRNA was noted in PTC compared to control (p = 0.009). Within the group of patients with conventional PTC, those with lymph nodal metastases had a higher level of VEGF-A mRNA expression than other patients (p = 0.0003). CONCLUSION These findings suggest that VEGF polymorphisms and mRNA expression may predict the aggressiveness behaviour of thyroid cancer.

A population-specific HTR2B stop codon predisposes to severe impulsivity

Impulsivity, describing action without foresight, is an important feature of several psychiatric diseases, suicidality and violent behaviour. The complex origins of impulsivity hinder identification of the genes influencing it and the diseases with which it is associated. Here we perform exon-focused sequencing of impulsive individuals in a founder population, targeting fourteen genes belonging to the serotonin and dopamine domain. A stop codon in HTR2B was identified that is common (minor allele frequency > 1%) but exclusive to Finnish people. Expression of the gene in the human brain was assessed, as well as the molecular functionality of the stop codon, which was associated with psychiatric diseases marked by impulsivity in both population and family-based analyses. Knockout of Htr2b increased impulsive behaviours in mice, indicative of predictive validity. Our study shows the potential for identifying and tracing effects of rare alleles in complex behavioural phenotypes using founder populations, and indicates a role for HTR2B in impulsivity.

G9a histone methyltransferase contributes to imprinting in the mouse placenta

Whereas DNA methylation is essential for genomic imprinting, the importance of histone methylation in the allelic expression of imprinted genes is unclear. Imprinting control regions (ICRs), however, are marked by histone H3-K9 methylation on their DNA-methylated allele. In the placenta, the paternal silencing along the Kcnq1 domain on distal chromosome 7 also correlates with the presence of H3-K9 methylation, but imprinted repression at these genes is maintained independently of DNA methylation. To explore which histone methyltransferase (HMT) could mediate the allelic H3-K9 methylation on distal chromosome 7, and at ICRs, we generated mouse conceptuses deficient for the SET domain protein G9a. We found that in the embryo and placenta, the differential DNA methylation at ICRs and imprinted genes is maintained in the absence of G9a. Accordingly, in embryos, imprinted gene expression was unchanged at the domains analyzed, in spite of a global loss of H3-K9 dimethylation (H3K9me2). In contrast, the placenta-specific imprinting of genes on distal chromosome 7 is impaired in the absence of G9a, and this correlates with reduced levels of H3K9me2 and H3K9me3. These findings provide the first evidence for the involvement of an HMT and suggest that histone methylation contributes to imprinted gene repression in the trophoblast.

Global dissemination of a multidrug resistant Escherichia coli clone

Escherichia coli sequence type 131 (ST131) is a globally disseminated, multidrug resistant (MDR) clone responsible for a high proportion of urinary tract and bloodstream infections. The rapid emergence and successful spread of E. coli ST131 is strongly associated with several factors, including resistance to fluoroquinolones, high virulence gene content, the possession of the type 1 fimbriae FimH30 allele, and the production of the CTX-M-15 extended spectrum β-lactamase (ESBL). Here, we used genome sequencing to examine the molecular epidemiology of a collection of E. coli ST131 strains isolated from six distinct geographical locations across the world spanning 2000–2011. The global phylogeny of E. coli ST131, determined from whole-genome sequence data, revealed a single lineage of E. coli ST131 distinct from other extraintestinal E. coli strains within the B2 phylogroup. Three closely related E. coli ST131 sublineages were identified, with little association to geographic origin. The majority of single-nucleotide variants associated with each of the sublineages were due to recombination in regions adjacent to mobile genetic elements (MGEs). The most prevalent sublineage of ST131 strains was characterized by fluoroquinolone resistance, and a distinct virulence factor and MGE profile. Four different variants of the CTX-M ESBL–resistance gene were identified in our ST131 strains, with acquisition of CTX-M-15 representing a defining feature of a discrete but geographically dispersed ST131 sublineage. This study confirms the global dispersal of a single E. coli ST131 clone and demonstrates the role of MGEs and recombination in the evolution of this important MDR pathogen.

The cytochrome P-450 isoenzyme CYP2E1 in the biological processing of industrial chemicals : consequences for occupational and environmental medicine

The importance of the isoform CYP2E1 of the human cytochrome P-450 superfamily of enzymes for occupational and environmental medicine is derived from its unique substrate spectrum that includes a number of highly important high-production chemicals, such as aliphatic and aromatic hydrocarbons, solvents and industrial monomers (i.a. alkanes, alkenes, aromatic and halogenated hydrocarbons). Many polymorphic genes, such as CYP2E1, show considerable differences in allelic distribution between different human populations. The polymorphic nature of the human CYP2E1 gene is significant for inter-individual differences in toxicity of its substrates. Since the substrate spectrum of CYP2E1 includes many compounds of basic relevance to industrial toxicology, a rationale for metabolic interactions of different CYP2E1 substrates is provided. In-depth research into the inter-individual phenotypic differences of human CYP2E1 enzyme activities was enabled by the recognition that the 6-hydroxylation of the drug chlorzoxazone is mediated by CYP2E1. Studies on CYP2E1 phenotyping have pointed to inter-individual variations in enzyme activities. There are consistent ethnic differences in CYP2E1 enzyme expression, mostly demonstrated between European and Japanese populations, which point to a major impact of genetic factors. The most frequently studied genetic polymorphisms are the restriction fragment length polymorphisms PstI/RsaI (mutant allele: CYP2E1*5B) located in the 5′-flanking region of the gene, as well as the DraI polymorphism (mutant allele: CYP2E1*6) located in intron 6. These polymorphisms are partly related, as they form the common allele designated CYP2E1*5A. Striking inter-ethnic differences between Europeans and Asians appear with respect to the frequencies of the CYP2E1*5A allele (only approximately 5% of Europeans are heterozygous, but 37% of Asians are, whilst 6% of Asians are homozygous). Available studies indicate a wide variation in human CYP2E1 expression, which are very likely based on complex gene-environment interactions. Major inter-ethnic differences are apparent on the genotyping and the phenotyping levels. Selected cases are presented where inter-ethnic variations of CYP2E1 may provide likely explanations for unexplained findings concerning industrial chemicals that are CYP2E1 substrates. Possible consequences of differential inter-individual and inter-ethnic susceptibilities are related to individual expressions of clinical symptoms of chemical toxicity, to results of biological monitoring of exposed workers, and to the interpretation of results of epidemiological or molecular-epidemiological studies.

Relevance of the deletion polymorphisms of the glutathione S-transferases GSTT1 and GSTM1 in pharmacology and toxicology

Although cytosolic glutathione S-transferase (GST) enzymes occupy a key position in biological detoxification processes, two of the most relevant human isoenzymes, GSTT1-1 and GSTM1-1, are genetically deleted (non-functional alleles GSTT1*0 and GSTM1*0) in a high percentage of the human population, with major ethnic differences. The structures of the GSTT and GSTM gene areas explain the underlying genetic processes. GSTT1-1 is highly conserved during evolution and plays a major role in phase-II biotransformation of a number of drugs and industrial chemicals, e.g. cytostatic drugs, hydrocarbons and halogenated hydrocarbons. GSTM1-1 is particularly relevant in the deactivation of carcinogenic intermediates of polycyclic aromatic hydrocarbons. Several lines of evidence suggest that hGSTT1-1 and/or hGSTM1-1 play a role in the deactivation of reactive oxygen species that are likely to be involved in cellular processes of inflammation, ageing and degenerative diseases. There is cumulating evidence that combinations of the GSTM1*0 state with other genetic traits affecting the metabolism of carcinogens (CYP1A1, GSTP1) may predispose the aero-digestive tract and lung, especially in smokers, to a higher risk of cancer. The GSTM1*0 status appears also associated with a modest increase in the risk of bladder cancer, consistent with a GSTM1 interaction with carcinogenic tobacco smoke constituents. Both human GST deletions, although largely counterbalanced by overlapping substrate affinities within the GST superfamily, have consequences when the organism comes into contact with distinct man-made chemicals. This appears relevant in industrial toxicology and in drug metabolism.

Glutathione transferase isozyme genotypes in patients with prostate and bladder carcinoma

Genotype distributions for GSTP1, GSTM1, and GSTT1 were determined in 91 patients with prostatic carcinoma and 135 patients with bladder carcinoma and compared with those in 127 abdominal surgery patients without malignancies. None of the genotypes differed significantly with respect to age or sex among controls or cancer patients. In the group of prostatic carcinoma patients, GSTT1 null allele homozygotes were more prevalent (25% in carcinoma patients vs 13% in controls, Fisher P=0.02, χ2 P = 0.02, OR = 2.31, CI = 1.17-4.59) and the combined M1-/T1-null genotype was also more frequent (9% vs 3%, χ2 P= 0.02, Fisher P = 0.03). Homozygosity for the GSTM1 null allele was more frequent among bladder carcinoma patients (59% in bladder carcinoma patients vs 45% in controls, Fisher P = 0.03, χ2 P = 0.02, OR = 1.76, CI = 1.08-2.88). In contrast to a previous report, no significant increase in the frequency of the GSTP1b allele was found in the tumor patients. Except for the combined GSTM1-/T1-null genotype in prostatic carcinoma, none of the combined genotypes showed a significant association with either of the cancers. These findings suggest that specific single polymorphic GST genes, that is GSTM1 in the case of bladder cancer and GSTT1 in the case of prostatic carcinoma, are most relevant for the development of these urological malignancies among the general population in Central Europe.

Influence of polymorphisms of the human glutathione transferases and cytochrome P450 2E1 enzyme on the metabolism and toxicity of ethylene oxide and acrylonitrile

A cohort of 59 persons with industrial handling of low levels of acrylonitrile is being studied as part of a medical surveillance programme. Previously, an extended haemoglobin adduct monitoring (N-(cyanoethyl)valine and N-(hydroxyethyl)-valine) was performed regarding the glutathione transferases hGSTM1 and hGSTT1 polymorphisms but no influence of hGSTM1 or hGSTT1 polymorphisms on specific adduct levels was found. A compilation of case reports of human accidental poisonings had pointed to significant individual differences in human acrylonitrile metabolism and toxicity. Therefore, a re-evaluation of the industrial cohort included known polymorphisms of the glutathione transferases hGSTM3 and hGSTP1 as well as of the cytochrome P450 CYP2E1. A detailed statistical analysis revealed that exposed carriers of the allelic variants of hGSTP1, hGSTP1*B/hGSTP1*C, characterized by a single nucleotide polymorphism at nucleotide 313 which results in a change from Ile to Val at codon 104, had higher levels of the acrylonitrile-specific haemoglobin adduct N-(cyanoethyl)valine compared to the carriers of the codon 113 alleles hGSTP1*A and hGSTP1*D. The single nucleotide polymorphism at codon 113 of hGSTP1 (hGSTP1*A/hGSTP1*B versus hGSTP1*C/hGSTP1*D) did not show an effect, and also no influence was seen on specific haemoglobin adduct levels of the polymorphisms of hGSTM3 or CYP2E1. The data, therefore, point to a possible influence of a human enzyme polymorphism of the GSTP1 gene at codon 104 on the detoxication of acrylonitrile which calls for experimental toxicological investigation. The study also confirmed the impact of GSTT1 polymorphism on background N-(hydroxyethyl)-valine adduct levels in haemoglobin which are caused by endogenous ethylene oxide.

Markers of genetic susceptibility in human environmental hygiene and toxicology : the role of selected CYP, NAT and GST genes

Inherited genetic traits co-determine the susceptibility of an individual to a toxic chemical. Special emphasis has been put on individual responses to environmental and industrial carcinogens, but other chronic diseases are of increasing interest. Polymorphisms of relevant xenobiotic metabolising enzymes may be used as toxicological susceptibility markers. A growing number of genes encoding enzymes involved in biotransformation of toxicants and in cellular defence against toxicant-induced damage to the cells has been identified and cloned, leading to increased knowledge of allelic variants of genes and genetic defects that may result in a differential susceptibility toward environmental toxicants. "Low penetrating" polymorphisms in metabolism genes tend to be much more common in the population than allelic variants of "high penetrating" cancer genes, and are therefore of considerable importance from a public health point of view. Positive associations between cancer and CYP1A1 alleles, in particular the *2C I462V allele, were found for tissues following the aerodigestive tract. Again, in most cases, the effect of the variant CYP1A1 allele becomes apparent or clearer in connection with the GSTM1 null allele. The CYP1B1 codon 432 polymorphism (CYP1B1*3) has been identified as a susceptibility factor in smoking-related head-and-neck squameous cell cancer. The impact of this polymorphic variant of CYP1B1 on cancer risk was also reflected by an association with the frequency of somatic mutations of the p53 gene. Combined genotype analysis of CYP1B1 and the glutathione transferases GSTM1 or GSTT1 has also pointed to interactive effects. Of particular interest for the industrial and environmental field is the isozyme CYP2E1. Several genotypes of this isozyme have been characterised which seem to be associated with different levels of expression of enzyme activity. The acetylator status for NAT2 can be determined by genotyping or by phenotyping. In the pathogenesis of human bladder cancer due to occupational exposure to "classical" aromatic amines (benzidine, 4-aminodiphenyl, 1-naphthylamine) acetylation by NAT2 is regarded as a detoxication step. Interestingly, the underlying European findings of a higher susceptibility of slow acetylators towards aromatic amines are in contrast to findings in Chinese workers occupationally exposed to aromatic amines which points to different mechanisms of susceptibility between European and Chinese populations. Regarding human bladder cancer, the hypothesis has been put forward that genetic polymorphism of GSTM1 might be linked with the occurrence of this tumour type. This supports the hypothesis that exposure to PAH might causally be involved in urothelial cancers. The human polymorphic GST catalysing conjugation of halomethanes, dihalomethanes, ethylene oxide and a number of other industrial compounds could be characterised as a class theta enzyme (GSTT1) by means of molecular biology. "Conjugator" and "non-conjugator" phenotypes are coincident with the presence and absence of the GSTT1 gene. There are wide variations in the frequencies of GSTT1 deletion (GSTT1 *0/0) among different ethnicities. Human phenotyping is facilitated by the GST activity towards methyl bromide or ethylene oxide in erythrocytes which is representative of the metabolic GSTT1 competence of the entire organism. Inter-individual variations in xenobiotic metabolism capacities may be due to polymorphisms of the genes coding for the enzymes themselves or of the genes coding for the receptors or transcription factors which regulate the expression of the enzymes. Also, polymorphisms in several regions of genes may cause altered ligand affinity, transactivation activity or expression levels of the receptor subsequently influencing the expression of the downstream target genes. Studies of individual susceptibility to toxicants and gene-environment interaction are now emerging as an important component of molecular epidemiology.

The hereditary transmission of the glutathione transferase hGSTT1-1 conjugator phenotype in a large family

The polymorphism of human glutathione transferase hGSTT1-1 is expressed in three phenotypes. Experimentally, individuals can be classified as non-conjugators, low conjugators and 'high' conjugators depending on the enzyme activity in blood towards methylene chloride using a gas chromatographic assay. Non-conjugators do not have a functional hGSTT1 gene; however, little is known about the molecular basis of the three conjugator phenotypes. The higher hGSTT1-1 activity in high conjugators may be the result of enzyme induction or be genetically determined. Twenty-nine members of a large family, including three generations were phenotyped and genotyped with respect to hGSTT1-1. The hGSTT1-1 enzyme activity of high conjugators was twice as high as that of low conjugators. The distribution of hGSTT1-1 phenotypes strongly indicates a Mendelian intermediary inheritance, in which a gene-dosage effect results in a doubled enzyme expression in the presence of two functional alleles. The Mendelian intermediary inheritance is further supported by the findings of a semiquantitative polymerase chain reaction method designed to distinguish the three genotypes of hGSTT1 for rapid screening of large study groups.

Evidence of associations between cytokine gene polymorphisms and quality of life in patients with cancer and their family caregivers

PURPOSE/OBJECTIVES: To identify latent classes of individuals with distinct quality-of-life (QOL) trajectories, to evaluate for differences in demographic characteristics between the latent classes, and to evaluate for variations in pro- and anti-inflammatory cytokine genes between the latent classes. DESIGN: Descriptive, longitudinal study. SETTING: Two radiation therapy departments located in a comprehensive cancer center and a community-based oncology program in northern California. SAMPLE: 168 outpatients with prostate, breast, brain, or lung cancer and 85 of their family caregivers (FCs). METHODS: Growth mixture modeling (GMM) was employed to identify latent classes of individuals based on QOL scores measured prior to, during, and for four months following completion of radiation therapy. Single nucleotide polymorphisms (SNPs) and haplotypes in 16 candidate cytokine genes were tested between the latent classes. Logistic regression was used to evaluate the relationships among genotypic and phenotypic characteristics and QOL GMM group membership. MAIN RESEARCH VARIABLES: QOL latent class membership and variations in cytokine genes. FINDINGS: Two latent QOL classes were found: higher and lower. Patients and FCs who were younger, identified with an ethnic minority group, had poorer functional status, or had children living at home were more likely to belong to the lower QOL class. After controlling for significant covariates, between-group differences were found in SNPs in interleukin 1 receptor 2 (IL1R2) and nuclear factor kappa beta 2 (NFKB2). For IL1R2, carrying one or two doses of the rare C allele was associated with decreased odds of belonging to the lower QOL class. For NFKB2, carriers with two doses of the rare G allele were more likely to belong to the lower QOL class. CONCLUSIONS: Unique genetic markers in cytokine genes may partially explain interindividual variability in QOL. IMPLICATIONS FOR NURSING: Determination of high-risk characteristics and unique genetic markers would allow for earlier identification of patients with cancer and FCs at higher risk for poorer QOL. Knowledge of these risk factors could assist in the development of more targeted clinical or supportive care interventions for those identified.

Association study of MTHFD1 coding polymorphisms R134K and R653Q with migraine susceptibility

Objective There is evidence that folate metabolism has a role in migraine pathophysiology, particularly in the migraine with aura subtype. In this study we investigate whether two non-synonymous single nucleotide polymorphisms (SNPs), rs1950902 (C401T; R134K) and rs2236225 (G1958A; R653Q), in MTHDF1 are associated with migraine in an Australian case-control population. Background Increased plasma levels of homocysteine (HCy), one of the metabolites produced in the folate pathway, has been found to be a risk factor for migraine. There is also a genetic link, as a common polymorphism (C667T) that reduces the catalytic activity of MTHFR, the enzyme that catalyses the formation of HCy, is associated with an increase in risk of the migraine with aura (MA) subtype. MTHFD1 is a crucial multifunctional enzyme that catalyses three separate reactions of the folate pathway and therefore variants in MTHFD1 may also influence migraine susceptibility. Methods The R134K and R653Q variants in MTHFD1 were genotyped in an Australian cohort of 520 unrelated migraineurs (162 were diagnosed with migraine without aura [MO] and 358 with MA) and 520 matched controls. Data were analysed for association with migraine and for interaction with the MTHFR C667T polymorphism. Results We find no significant differences in genotype or allele frequencies for either SNP between migraineurs and controls, or when either MO or MA cases were compared to controls. In addition these MTHFD1 polymorphisms did not appear to influence the risk of MA conferred by the MTHFR 667T allele. Conclusions We find no evidence for association of the MTHFD1 R134K and R653Q polymorphisms with migraine in our Australian case-control population. However, as folate metabolism appears to be important in migraine, particularly with respect to the aura component, future studies using high throughput methods to expand the number of SNPs in folate-related genes genotyped and investigation of interactions between SNPs may be justified.

Association of the SNP rs2623047 in the HSPG modification enzyme SULF1 with an Australian Caucasian Breast Cancer Cohort

Breast cancer is the second most common cancer worldwide and the most common cancer reported in women. This malignant tumour is characterised by a number of specific features including uncontrolled cell proliferation. It ranks fifth in the world as a cause of cancer death in women. Early diagnosis increases 5 year survival rates up to 95%. Heparan sulfate proteoglycans (HSPGs) are complex proteins composed of a core protein to which a number of highly sulfated side chains are synthesised by a highly co-ordinated process resulting in distinct sulfation patterns, which determine specific interations with cell-signaling partners including growth factors, their receptors, ligands and morphogens. The enzymes responsible for chain initiation, elongation and sulfation are critical for creating HS chain variability conferring biological functionality. This study investigated single nucleotide polymorphism in SULF1, the enzyme responsible for the 6-0 desulfation of heparan sulfate side chains. We investigated this SNP in an Australian Caucasian case-control breast cancer population and found a significant association between SULF1 and breast cancer at both the allelic and genotypic level (allele, p=0.016; genotype, p=0.032). Our results suggest the res2623047 SNP in SULF1 may impact breast cancer susceptibility. Specifically, the T allele of rs2623047 in SULF1 is associated with a increased risk of developing breast cancer in our cohort. The identification of markers including SULF1 may improve detection of this disease at its earliest stages improving patient treatment and prognosis.