Cross-talk between one-carbon metabolism and xenobiotic metabolism: implications on oxidative DNA damage and susceptibility to breast cancer

The aim of this case–control study is to explore the role of aberrations in xenobiotic metabolism in inducing oxidative DNA damage and altering the susceptibility to breast cancer. Cytochrome P4501A1 (CYP1A1) m1 (OR: 1.41, 95% CI 1.08–1.84), CYP1A1 m4 (OR: 5.13, 95% CI 2.68–9.81), Catecholamine-O-methyl transferase (COMT) H108L (OR: 1.49, 95% CI 1.16–1.92), and glutathione S-transferase (GST) T1 null (OR: 1.68, 95% CI 1.09–2.59) variants showed association with breast cancer risk. Reduced folate carrier 1 (RFC1) 80A/CYP1A1 m1/CYP1A1 m4 and RFC1 80A/thymidylate synthase (TYMS) 5′-UTR 2R/methionine synthase (MTR) 2756G/COMT 108L genetic combinations were found to inflate breast cancer risk under the conditions of low dietary folate (345 ± 110 vs. 379 ± 139 μg/day) and low plasma folate (6.81 ± 1.25 vs. 7.09 ± 1.26 ng/ml) by increasing plasma 8-oxo-2′-deoxyguanosine (8-oxodG). This increase in 8-oxodG is attributed to low methionine (49.38 ± 23.74 vs. 53.90 ± 23.85 μmol/l); low glutathione (378 ± 242 vs. 501 ± 126 μmol/l) and GSTT1 null variant; and hypermethylation of CpG island of extracellular-superoxide dismutase (EC-SOD) (92.78 ± 11.49 vs. 80.45 ± 9.86%), which impair O-methylation of catechol estrogens to methoxy estrogens, conjugation of glutathione to semiquinones/quinones and free radical scavenging respectively. Our results suggest cross-talk between one-carbon metabolism and xenobiotic metabolism influencing oxidative DNA damage and susceptibility to breast cancer.     

Epistatic interactions between four rad loci in yeast

Haploid yeast strains carrying mutations in two or more of four ad genes were contrusted by tetrad dissection, and the UV survival of these strains was measured. It was found that (with one exception) double mutant strains were not significantly more sensitive than the most sensitive single mutants, for strains involving mutant loci rad 1, rad 3 and rad 4. The exception was the double mutant rad 1–5 rad 4-4, but another double mutant involving different alleles of the the same loci did not show an enhanced UV sensitivity. Triple and quadruple mutants also failed to show a significantly increased UV sensitivity with respect to the single mutants. The results indicate that all these four mutant loci confer UV sensitivity by the same mechanism, and it is suggested that the wild-type alleles mediate excision-repair of UV-induced DNA lesions. Enhanced sensitivity of the genotype rad 1–5 rad 4-4 is attributed to leakiness of these alleles.     

Genetic loci controlling breast cancer susceptibility in the Wistar-Kyoto rat

In this study, the Wistar-Kyoto (WKy) rat was genetically characterized for loci that modify susceptibility to mammary carcinogenesis. We used a genetic backcross between resistant WKy and susceptible Wistar-Furth (WF) rats as a panel for linkage mapping to genetically identify mammary carcinoma susceptibility (Mcs) loci underlying the resistance of the WKy rat. Rats were phenotyped for DMBA-induced mammary carcinomas and genotyped using microsatellite markers. To detect quantitative trait loci (QTL), we analyzed the genome scan data under both parametric and nonparametric distributional assumptions and used permutation tests to calculate significance thresholds. A generalized linear model analysis was also performed to test for interactions between significant QTL. This methodology was extended to identify interactions between the significant QTL and other genome locations. Chromosomes 5, 7, 10, and 14 were found to contain significant QTL, termed Mcs5, Mcs6, Mcs7, and Mcs8, respectively. The WKy alleles of Mcs5, -6, and -8 are associated with mammary carcinoma resistance; the WKy allele of Mcs7 is associated with an increased incidence of mammary cancer. In addition, we identified an interaction between Mcs8 and a region on chromosome 6 termed Mcsm1 (modifier of Mcs), which had no significant main effect on mammary cancer susceptibility in this genetic analysis.     

The influence of one-carbon metabolism on gene promoter methylation in a population-based breast cancer study

Abnormal methylation in gene promoters is a hallmark of the cancer genome; however, factors that may influence promoter methylation have not been well elucidated. As the one-carbon metabolism pathway provides the universal methyl donor for methylation reactions, perturbation of this pathway might influence DNA methylation and, ultimately, affect gene functions. Utilizing approximately 800 breast cancer tumor tissues from a large population-based study, we investigated the relationships between dietary and genetic factors involved in the one-carbon metabolism pathway and promoter methylation of a panel of 13 breast cancer-related genes. We found that CCND2, HIN1 and CHD1 were the most “dietary sensitive” genes, as methylation of their promoters was associated with intakes of at least two out of the eight dietary methyl factors examined. On the other hand, some micronutrients (i.e., B₂ and B₆) were more “epigenetically active” as their intake levels correlated with promoter methylation status in 3 out of the 13 breast cancer genes evaluated. Both positive (hypermethylation) and inverse (hypomethylation) associations with high micronutrient intake were observed. Unlike what we saw for dietary factors, we did not observe any clear patterns between one-carbon genetic polymorphisms and the promoter methylation status of the genes examined. Our results provide preliminary evidence that one-carbon metabolism may have the capacity to influence the breast cancer epigenome. Given that epigenetic alterations are thought to occur early in cancer development and are potentially reversible, dietary modifications may offer promising venues for cancer intervention and prevention.     

The influence of one-carbon metabolism on gene promoter methylation in a population-based breast cancer study

Abnormal methylation in gene promoters is a hallmark of the cancer genome; however, factors that may influence promoter methylation have not been well elucidated. As the one-carbon metabolism pathway provides the universal methyl donor for methylation reactions, perturbation of this pathway might influence DNA methylation and, ultimately, affect gene functions. Utilizing approximately 800 breast cancer tumor tissues from a large population-based study, we investigated the relationships between dietary and genetic factors involved in the one-carbon metabolism pathway and promoter methylation of a panel of 13 breast cancer-related genes. We found that CCND2, HIN1 and CHD1 were the most "dietary sensitive" genes, as methylation of their promoters was associated with intakes of at least two out of the eight dietary methyl factors examined. On the other hand, some micronutrients (i.e., B 2 and B 6) were more "epigenetically active" as their intake levels correlated with promoter methylation status in 3 out of the 13 breast cancer genes evaluated. Both positive (hypermethylation) and inverse (hypomethylation) associations with high micronutrient intake were observed. Unlike what we saw for dietary factors, we did not observe any clear patterns between one-carbon genetic polymorphisms and the promoter methylation status of the genes examined. Our results provide preliminary evidence that one-carbon metabolism may have the capacity to influence the breast cancer epigenome. Given that epigenetic alterations are thought to occur early in cancer development and are potentially reversible, dietary modifications may offer promising venues for cancer intervention and prevention.     

Use of monozygotic twins in search for breast cancer susceptibility loci.

We have used Swedish monozygotic twins concordant for breast cancer to study genetic changes associated with the development of breast cancer. Because loss of heterozygosity (LOH) at a specific genomic region may reflect the presence of a tumour suppressor gene, loss of the same allele in both of the twins concordant for breast cancer may pinpoint a tumour suppressor gene that confers a strong predisposition to breast cancer. DNA samples extracted from the matched tumour and normal tissues of nine twin pairs were analysed for allelic imbalance using a set of microsatellite markers on chromosomes 1, 13, 16 and 17, containing loci with known tumour suppressor genes. The two main regions, where more twin pairs than expected had lost the same allele, were located at 16qtel', including markers D16S393, D16S305 and D16S413, and at 17p13, distal to the p53 locus. Our results show that the monozygotic twin model can be used to suggest candidate regions of potential tumour suppressor genes, even with a limited number of twin pairs.     

Reprogramming of serine,glycine and one-carbon metabolism in cancer

Metabolic pathways leading to the synthesis, uptake, and usage of the nonessential amino acid serine are frequently amplified in cancer. Serine encounters diverse fates in cancer cells, including being charged onto tRNAs for protein synthesis, providing head groups for sphingolipid and phospholipid synthesis, and serving as a precursor for cellular glycine and one-carbon units, which are necessary for nucleotide synthesis and methionine cycle reloading. This review will focus on the participation of serine and glycine in the mitochondrial one-carbon (SGOC) pathway during cancer progression, with an emphasis on the genetic and epigenetic determinants that drive SGOC gene expression. We will discuss recently elucidated roles for SGOC metabolism in nucleotide synthesis, redox balance, mitochondrial function, and epigenetic modifications. Finally, therapeutic considerations for targeting SGOC metabolism in the clinic will be discussed.     

Genome-wide association study in east Asians identifies novel susceptibility loci for breast cancer

Genetic factors play an important role in the etiology of both sporadic and familial breast cancer. We aimed to discover novel genetic susceptibility loci for breast cancer. We conducted a four-stage genome-wide association study (GWAS) in 19,091 cases and 20,606 controls of East-Asian descent including Chinese, Korean, and Japanese women. After analyzing 690,947 SNPs in 2,918 cases and 2,324 controls, we evaluated 5,365 SNPs for replication in 3,972 cases and 3,852 controls. Ninety-four SNPs were further evaluated in 5,203 cases and 5,138 controls, and finally the top 22 SNPs were investigated in up to 17,423 additional subjects (7,489 cases and 9,934 controls). SNP rs9485372, near the TGF-β activated kinase (TAB2) gene in chromosome 6q25.1, showed a consistent association with breast cancer risk across all four stages, with a P-value of 3.8×10⁻¹² in the combined analysis of all samples. Adjusted odds ratios (95% confidence intervals) were 0.89 (0.85–0.94) and 0.80 (0.75–0.86) for the A/G and A/A genotypes, respectively, compared with the genotype G/G. SNP rs9383951 (P = 1.9×10⁻⁶ from the combined analysis of all samples), located in intron 5 of the ESR1 gene, and SNP rs7107217 (P = 4.6×10⁻⁷), located at 11q24.3, also showed a consistent association in each of the four stages. This study provides strong evidence for a novel breast cancer susceptibility locus represented by rs9485372, near the TAB2 gene (6q25.1), and identifies two possible susceptibility loci located in the ESR1 gene and 11q24.3, respectively.     

Heterogeneity of breast cancer associations with five susceptibility loci by clinical and pathological characteristics

A three-stage genome-wide association study recently identified single nucleotide polymorphisms (SNPs) in five loci (fibroblast growth receptor 2 (FGFR2), trinucleotide repeat containing 9 (TNRC9), mitogen-activated protein kinase 3 K1 (MAP3K1), 8q24, and lymphocyte-specific protein 1 (LSP1)) associated with breast cancer risk. We investigated whether the associations between these SNPs and breast cancer risk varied by clinically important tumor characteristics in up to 23,039 invasive breast cancer cases and 26,273 controls from 20 studies. We also evaluated their influence on overall survival in 13,527 cases from 13 studies. All participants were of European or Asian origin. rs2981582 in FGFR2 was more strongly related to ER-positive (per-allele OR (95%CI)=1.31 (1.27–1.36)) than ER-negative (1.08 (1.03–1.14)) disease (P for heterogeneity=10⁻¹³). This SNP was also more strongly related to PR-positive, low grade and node positive tumors (P=10⁻⁵, 10⁻⁸, 0.013, respectively). The association for rs13281615 in 8q24 was stronger for ER-positive, PR-positive, and low grade tumors (P=0.001, 0.011 and 10⁻⁴, respectively). The differences in the associations between SNPs in FGFR2 and 8q24 and risk by ER and grade remained significant after permutation adjustment for multiple comparisons and after adjustment for other tumor characteristics. Three SNPs (rs2981582, rs3803662, and rs889312) showed weak but significant associations with ER-negative disease, the strongest association being for rs3803662 in TNRC9 (1.14 (1.09–1.21)). rs13281615 in 8q24 was associated with an improvement in survival after diagnosis (per-allele HR=0.90 (0.83–0.97). The association was attenuated and non-significant after adjusting for known prognostic factors. Our findings show that common genetic variants influence the pathological subtype of breast cancer and provide further support for the hypothesis that ER-positive and ER-negative disease are biologically distinct. Understanding the etiologic heterogeneity of breast cancer may ultimately result in improvements in prevention, early detection, and treatment.     

Genetic identification of multiple loci that control breast cancer susceptibility in the rat.

We have used a rat model of induced mammary carcinomas in an effort to identify breast cancer susceptibility genes. Using genetic crosses between the carcinoma-resistant Copenhagen (COP) and carcinoma-sensitive Wistar-Furth rats, we have confirmed the identification of the Mcs1 locus that modulates tumor number. We have now also identified two additional loci, Mcs2 and Mcs3. These three loci map to chromosomes 2, 7, and 1, respectively, and interact additively to suppress mammary carcinoma development in the COP strain. They are responsible for a major portion of the tumor-resistant phenotype of the COP rat. No loss of heterozygosity was observed surrounding the three loci. A fourth COP locus, Mcs4, has also been identified on chromosome 8 and acts in contrast to increase the number of carcinomas. These results show that mammary carcinoma susceptibility in the COP rat is a polygenic trait. Interestingly, a polymorphism in the human genomic region homologous to the rat Mcs4 region is associated with an increased breast cancer risk in African-American women. The isolation of the Mcs genes may help elucidate novel mechanisms of carcinogenesis, provide information important for human breast cancer risk estimation, and also provide unique drug discovery targets for breast cancer prevention.     

Associations between genetic variation in one-carbon metabolism and LINE-1 DNA methylation in histologically normal breast tissues

Genome-wide DNA hypomethylation is an early event in the carcinogenic process. Percent methylation of long interspersed nucleotide element-1 (LINE-1) is a biomarker of genome-wide methylation and is a potential biomarker for breast cancer. Understanding factors associated with percent LINE-1 DNA methylation in histologically normal tissues could provide insight into early stages of carcinogenesis. In a cross-sectional study of 121 healthy women with no prior history of cancer who underwent reduction mammoplasty, we examined associations between plasma and breast folate, genetic variation in one-carbon metabolism, and percent LINE-1 methylation using multivariable regression models (adjusting for race, oral contraceptive use, and alcohol use). Results are expressed as the ratio of LINE-1 methylation relative to that of the referent group, with the corresponding 95% confidence intervals (CI). We found no significant associations between plasma or breast folate and percent LINE-1 methylation. Variation in MTHFR, MTR, and MTRR were significantly associated with percent LINE-1 methylation. Variant allele carriers of MTHFR A1289C had 4% lower LINE-1 methylation (Ratio 0.96, 95% CI 0.93–0.98), while variant allele carriers of MTR A2756G (Ratio 1.03, 95% CI 1.01–1.06) and MTRR A66G (Ratio 1.03, 95% CI 1.01–1.06) had 3% higher LINE-1 methylation, compared to those carrying the more common genotypes of these SNPs. DNA methylation of LINE-1 elements in histologically normal breast tissues is influenced by polymorphisms in genes in the one-carbon metabolism pathway. Future studies are needed to investigate the sociodemographic, environmental and additional genetic determinants of DNA methylation in breast tissues and the impact on breast cancer susceptibility.     

Epistatic and synergistic interactions between circadian clock mutations in Neurospora crassa

We identified a series of epistatic and synergistic interactions among the circadian clock mutations of Neurospora crassa that indicate possible physical interactions among the various clock components encoded by these genes. The period-6 (prd-6) mutation, a short-period temperature-sensitive clock mutation, is epistatic to both the prd-2 and prd-3 mutations. The prd-2 and prd-3 long-period mutations show a synergistic interaction in that the period length of the double mutant strain is considerably longer than predicted. In addition, the prd-2 prd-3 double mutant strain also exhibits overcompensation to changes in ambient temperature, suggesting a role in the temperature compensation machinery of the clock. The prd-2, prd-3, and prd-6 mutations also show significant interactions with the frq(7) long-period mutation. These results suggest that the gene products of prd-2, prd-3, and prd-6 play an important role in both the timing and temperature compensation mechanisms of the circadian clock and may interact with the FRQ protein.     

Epistatic partners of neurogenic genes modulate Drosophila olfactory behavior

The extent to which epistasis affects the genetic architecture of complex traits is difficult to quantify, and identifying variants in natural populations with epistatic interactions is challenging. Previous studies in Drosophila implicated extensive epistasis between variants in genes that affect neural connectivity and contribute to natural variation in olfactory response to benzaldehyde. In this study, we implemented a powerful screen to quantify the extent of epistasis as well as identify candidate interacting variants using 203 inbred wild‐derived lines with sequenced genomes of the Drosophila melanogaster Genetic Reference Panel (DGRP). We crossed the DGRP lines to P[GT1]‐element insertion mutants in Sema‐5c and neuralized (neur), two neurodevelopmental loci which affect olfactory behavior, and to their coisogenic wild‐type control. We observed significant variation in olfactory responses to benzaldehyde among F₁ genotypes and for the DGRP line by mutant genotype interactions for both loci, showing extensive nonadditive genetic variation. We performed genome‐wide association analyses to identify the candidate modifier loci. None of these polymorphisms were in or near the focal genes; therefore, epistasis is the cause of the nonadditive genetic variance. Candidate genes could be placed in interaction networks. Several candidate modifiers are associated with neural development. Analyses of mutants of candidate epistatic partners with neur (merry‐go‐round (mgr), prospero (pros), CG10098, Alhambra (Alh) and CG12535) and Sema‐5c (CG42540 and bruchpilot (brp)) showed aberrant olfactory responses compared with coisogenic controls. Thus, integrating genome‐wide analyses of natural variants with mutations at defined genomic locations in a common coisogenic background can unmask specific epistatic modifiers of behavioral phenotypes.     

Contrasting epistatic interactions between rat quantitative trait loci controlling mammary cancer development

We previously defined quantitative trait loci (QTLs) that control susceptibility to 7,12-dimethylbenz(α)anthracene-induced mammary carcinoma in SPRD-Cu3 (susceptible) and WKY (resistant) rats. Two of these QTLs, assigned to chromosomes (Chr) 10 and 18, control tumor growth rate and invasiveness. In this study we characterized a congenic strain in which a large segment of WKY Chr 10 was introduced in the SPRD-Cu3 genetic background and demonstrated that this chromosome segment controls this tumor trait. The WKY allele at this QTL (Mcsta1) reduces the growth rate of the fastest growing tumors by 26%. We also previously showed that two SPRD-Cu3-WKY congenic strains containing a WKY chromosome segment derived either from Chr 5 or from Chr 18 exhibit a reduction in tumor multiplicity (QTLs Msctm1 and Mcstm2, respectively) (with no reduction in tumor growth rate in the Chr 18 congenic). In this study we generated a double congenic strain, which contains the two WKY differential segments from Chr 5 and 18, to determine how these two segments interact with one another. Interestingly, two types of epistatic interactions were found: no additive effect was seen with respect to tumor multiplicity, while a reduction in tumor growth rate was observed. It thus appears that WKY alleles located on Chr 5 and Chr 8 interact epistatically in a contrasting manner to modulate tumor multiplicity (in a nonadditive manner) and growth rate (in a synergic manner). Tumor growth rate is thus influenced by two QTLs, on Chr 10 (Mcsta1) and on Chr 18 (Mcsta2), the action of the latter being dependent on the presence of the Chr5 QTL (Mcstm1). The expression level of positional and functional candidate genes was also analyzed. On Chr 5, Pla2g2a is subject to a syntenic control while expression of the Tp53 (Chr 10) and Pmai1/Noxa (Chr 18) genes appears to be controlled by several mammary cancer resistance QTLs.     

Relationships among biomarkers of one-carbon metabolism

One-carbon metabolism is a network of metabolic pathways, disruption of which has been associated with cancer and other pathological conditions. Biomarkers of these pathways include homocysteine (HCY), S-adenosylmethionine (SAM), and S-adenosylhomocysteine (SAH). A better understanding of the relationships between these biomarkers is needed for their utilization in research. This study investigated the relationships between fasting concentrations of plasma HCY, SAM, SAH and the ratio of SAM:SAH, and serum folate, vitamin B(12) and creatinine in a healthy adult population. A cross-sectional study recruited 678 volunteers; only subjects with complete data (n = 581) were included in this analysis. Correlations were used to examine bivariate relationships among the biomarkers and multivariate linear regression determined independent relationships with HCY, SAM and SAH treated as dependent variables in separate models. Multivariate logistic regression examined determinants of a low SAM:SAH ratio (defined as having a SAM:SAH ratio in the bottom quartile and SAH value in the top quartile). HCY correlated inversely with folate and vitamin B(12) and weakly correlated with SAH and creatinine. Both SAM and SAH correlated with creatinine but were independent of serum folate and vitamin B(12). In multivariate analyses, folate, vitamin B(12), creatinine, sex and age were associated with HCY; age and creatinine were determinants of SAM, and sex and creatinine determinants of SAH. Finally, male sex and increasing creatinine levels were associated with having a low SAM:SAH ratio. Findings suggest that HCY, SAM and SAH are relatively independent parameters and reflect distinct aspects of one-carbon metabolism.     

Comprehensive evaluation of one-carbon metabolism pathway gene variants and renal cell cancer risk

Introduction

Folate and one-carbon metabolism are linked to cancer risk through their integral role in DNA synthesis and methylation. Variation in one-carbon metabolism genes, particularly MTHFR, has been associated with risk of a number of cancers in epidemiologic studies, but little is known regarding renal cancer.

Methods

Tag single nucleotide polymorphisms (SNPs) selected to produce high genomic coverage of 13 gene regions of one-carbon metabolism (ALDH1L1, BHMT, CBS, FOLR1, MTHFR, MTR, MTRR, SHMT1, SLC19A1, TYMS) and the closely associated glutathione synthesis pathway (CTH, GGH, GSS) were genotyped for 777 renal cell carcinoma (RCC) cases and 1,035 controls in the Central and Eastern European Renal Cancer case-control study. Associations of individual SNPs (n = 163) with RCC risk were calculated using unconditional logistic regression adjusted for age, sex and study center. Minimum p-value permutation (Min-P) tests were used to identify gene regions associated with risk, and haplotypes were evaluated within these genes.

Results

The strongest associations with RCC risk were observed for SLC19A1 (P_min-P = 0.03) and MTHFR (P_min-P = 0.13). A haplotype consisting of four SNPs in SLC19A1 (rs12483553, rs2838950, rs2838951, and rs17004785) was associated with a 37% increased risk (p = 0.02), and exploratory stratified analysis suggested the association was only significant among those in the lowest tertile of vegetable intake.

Conclusions

To our knowledge, this is the first study to comprehensively examine variation in one-carbon metabolism genes in relation to RCC risk. We identified a novel association with SLC19A1, which is important for transport of folate into cells. Replication in other populations is required to confirm these findings.     

Epistatic interactions between ancestral genotype and beneficial mutations shape evolvability in Pseudomonas aeruginosa

The idea that interactions between mutations influence adaptation by driving populations to low and high fitness peaks on adaptive landscapes is deeply ingrained in evolutionary theory. Here, we investigate the impact of epistasis on evolvability by challenging populations of two Pseudomonas aeruginosa clones bearing different initial mutations (in rpoB conferring rifampicin resistance, and the type IV pili gene network) to adaptation to a medium containing l ‐serine as the sole carbon source. Despite being initially indistinguishable in fitness, populations founded by the two ancestral genotypes reached different fitness following 300 generations of evolution. Genome sequencing revealed that the difference could not be explained by acquiring mutations in different targets of selection; the majority of clones from both ancestors converged on one of the following two strategies: (1) acquiring mutations in either PA2449 (gcsR, an l ‐serine‐metabolism RpoN enhancer binding protein) or (2) protease genes. Additionally, populations from both ancestors converged on loss‐of‐function mutations in the type IV pili gene network, either due to ancestral or acquired mutations. No compensatory or reversion mutations were observed in RNA polymerase (RNAP) genes, in spite of the large fitness costs typically associated with mutations in rpoB. Although current theory points to sign epistasis as the dominant constraint on evolvability, these results suggest that the role of magnitude epistasis in constraining evolvability may be underappreciated. The contribution of magnitude epistasis is likely to be greatest under the biologically relevant mutation supply rates that make back mutations probabilistically unlikely.