David J. Samuelson

Rat Mcs5a is a compound quantitative trait locus with orthologous human loci that associate with breast cancer risk

Breast cancer risk is a polygenic trait. To identify breast cancer modifier alleles that have a high population frequency and low penetrance we used a comparative genomics approach. Quantitative trait loci (QTL) were initially identified by linkage analysis in a rat mammary carcinogenesis model followed by verification in congenic rats carrying the specific QTL allele under study. The Mcs5a locus was identified by fine-mapping Mcs5 in a congenic model. Here we characterize the Mcs5a locus, which when homozygous for the Wky allele, reduces mammary cancer risk by 50%. The Mcs5a locus is a compound QTL with at least two noncoding interacting elements: Mcs5a1 and Mcs5a2. The resistance phenotype is only observed in rats carrying at least one copy of the Wky allele of each element on the same chromosome. Mcs5a1 is located within the ubiquitin ligase Fbxo10, whereas Mcs5a2 includes the 5′ portion of Frmpd1. Resistant congenic rats show a down-regulation of Fbxo10 in the thymus and an up-regulation of Frmpd1 in the spleen. The association of the Mcs5a1 and Mcs5a2 human orthologs with breast cancer was tested in two population-based breast cancer case-control studies (≈12,000 women). The minor alleles of rs6476643 (MCS5A1) and rs2182317 (MCS5A2) were independently associated with breast cancer risk. The minor allele of rs6476643 increases risk, whereas the rs2182317 minor allele decreases risk. Both alleles have a high population frequency and a low penetrance toward breast cancer risk.

Fine Mapping Reveals Multiple Loci and a Possible Epistatic Interaction within the Mammary Carcinoma Susceptibility Quantitative Trait Locus, Mcs5

To identify high-frequency, low-penetrance breast cancer modifier genes, we have developed a rat genetic model that uses the Wistar-Kyoto (WKy) inbred strain, resistant to developing 7,12-dimethylbenz[a]anthracene-induced mammary carcinogenesis, as a congenic donor and the susceptible Wistar-Furth (WF) strain as the recipient. Here, data from congenic rat lines containing smaller WKy genomic intervals of the Mcs5 quantitative trait locus region are presented to fine map three independently acting Mcs5 subloci. WKy-homozygous females from congenic lines defining Mcs5a, Mcs5b, and Mcs5c averaged, respectively, 4.0 +/- 0.4, 11.6 +/- 0.6, and 3.5 +/- 0.4 mammary carcinomas per rat. These phenotypic values are statistically different from the WF-homozygous phenotype value of 8.0 +/- 0.4, which is the baseline phenotype used for these experiments. We identified a likely Mcs5a x Mcs5b epistatic interaction that results in masking the increased susceptibility effect of the Mcs5b WKy allele by the Mcs5a WKy allele. We also provide evidence for a Mcs5a x Mcs5c interaction that is synergistic to decrease mammary carcinoma susceptibility below the additive effects of WKy alleles at each locus independently. The Mcs5 subloci are currently localized to 1.0, 7.5, and 4.5 Mb of rat chromosome 5, and the orthologous regions are on human chromosome 9 and mouse chromosome 4. These loci will provide unbiased candidate gene loci for evaluation in human case-control association studies.

The non-protein coding breast cancer susceptibility locus Mcs5a acts in a non-mammary cell-autonomous fashion through the immune system and modulates T-cell homeostasis and functions

IntroductionMechanisms underlying low-penetrance, common, non-protein coding variants in breast cancer risk loci are largely undefined. We showed previously that the non-protein coding mammary carcinoma susceptibility locus Mcs5a/MCS5A modulates breast cancer risk in rats and women. The Mcs5a allele from the Wistar-Kyoto (WKy) rat strain consists of two genetically interacting elements that have to be present on the same chromosome to confer mammary carcinoma resistance. We also found that the two interacting elements of the resistant allele are required for the downregulation of transcript levels of the Fbxo10 gene specifically in T-cells. Here we describe mechanisms through which Mcs5a may reduce mammary carcinoma susceptibility.MethodsWe performed mammary carcinoma multiplicity studies with three mammary carcinoma-inducing treatments, namely 7,12-dimethylbenz(a)anthracene (DMBA) and N-nitroso-N-methylurea (NMU) carcinogenesis, and mammary ductal infusion of retrovirus expressing the activated HER2/neu oncogene. We used mammary gland and bone marrow transplantation assays to assess the target tissue of Mcs5a activity. We used immunophenotyping assays on well-defined congenic rat lines carrying susceptible and resistant Mcs5a alleles to identify changes in T-cell homeostasis and function associated with resistance.ResultsWe show that Mcs5a acts beyond the initial step of mammary epithelial cell transformation, during early cancer progression. We show that Mcs5a controls susceptibility in a non-mammary cell-autonomous manner through the immune system. The resistant Mcs5a allele was found to be associated with an overabundance of gd T-cell receptor (TCR)+ T-cells as well as a CD62L (L-selectin)-high population of all T-cell classes. In contrast to in mammary carcinoma, gdTCR+ T-cells are the predominant T-cell type in the mammary gland and were found to be overabundant in the mammary epithelium of Mcs5a resistant congenic rats. Most of them simultaneously expressed the CD4, CD8, and CD161α markers. In cultured T-cells of Mcs5a resistant congenic rats we found increased mitogen-induced proliferation and production of Th1 cytokines IFNg, IL-2, and Tumor Necrosis Factor (TNF), but not Th2 cytokines IL-4 and IL-6, or Th17 cytokine IL-17 when compared with susceptible control rats.ConclusionsThese data support a hypothesis that Mcs5a displays a non-mammary cell-autonomous mechanism of action to modulate breast cancer risk through the immune system. The resistant Mcs5a allele is associated with alterations in T-cell homeostasis and functions, and overabundance of γδTCR+ T-cells in carcinogen-exposed mammary epithelium.

An insulator loop resides between the synthetically interacting elements of the human/rat conserved breast cancer susceptibility locus MCS5A/Mcs5a

Many low-penetrance breast cancer susceptibility loci are found to be located in non-protein-coding regions, suggesting their involvement in gene expression regulation. We identified the human/rat-conserved breast cancer susceptibility locus MCS5A/Mcs5a. This locus has been shown to act in a non-mammary cell-autonomous fashion through the immune system. The resistant Mcs5a allele from the Wistar–Kyoto (WKy) rat strain consists of two non-protein-coding genetic elements that must be located on the same chromosome to elicit the phenotype. In this study, we show the presence of a conserved higher order chromatin structure in MCS5A/Mcs5a located in between the synthetically interacting genetic elements. The looped elements are shown to be bound by CTCF and cohesin. We identify the downregulation of Fbxo10 expression in T cells as a strong candidate mechanism through which the interacting genetic elements of the resistant Mcs5a allele modulate mammary carcinoma susceptibility. Finally, we show that the human MCS5A polymorphisms associated with breast cancer risk are located at both sides of the looped structure and functionally interact to downregulate transcriptional activity, similar to rat Mcs5a. We propose a mechanistic model for MCS5a/Mcs5a in which a CTCF-mediated insulator loop encompassing the TOMM5/Tomm5 gene, resides in between and brings into closer physical proximity the synthetically and functionally interacting resistant genetic variants.

Rat Mcs1b Is Concordant to the Genome-Wide Association-Identified Breast Cancer Risk Locus at Human 5q11.2 and MIER3 is a Candidate Cancer Susceptibility Gene

Low-penetrance alleles associated with breast cancer risk have been identified in population-based studies. Most risk loci contain either no or multiple potential candidate genes. Rat mammary carcinoma susceptibility 1b (Mcs1b) is a quantitative trait locus on RN02 that confers decreased susceptibility when Copenhagen (COP)-resistant alleles are introgressed into a Wistar Furth (WF)-susceptible genome. Five WF.COP congenic lines containing COP RN02 segments were compared. One line developed an average of 3.4 ± 2.0 and 5.5 ± 3.6 mammary carcinomas per rat ± SD when females were Mcs1b-resistant homozygous and Mcs1b heterozygous, respectively. These phenotypes were significantly different from susceptible genotype littermates (7.8 ± 3.1 mean mammary carcinomas per rat ± SD, P = 0.0001 and P = 0.0413, respectively). All other congenic lines tested were susceptible. Thus, Mcs1b was narrowed to 1.8 Mb of RN02 between genetic markers ENSRNOSNP2740854 and g2UL2-27. Mammary gland-graft carcinoma susceptibility assays were used to determine that donor (P = 0.0019), but not recipient Mcs1b genotype (P = 0.9381), was associated with ectopic mammary carcinoma outcome. Rat Mcs1b contains sequence orthologous to human 5q11.2, a breast cancer susceptibility locus identified in multiple genome-wide association studies. Human/rat MAP3K1/Map3k1 and mesoderm induction early response (MIER; MIER3)/MIER3 are within these orthologous segments. We identified MIER3 as a candidate Mcs1b gene based on 4.5-fold higher mammary gland levels of MIER3 transcripts in susceptible compared with Mcs1b-resistant females. These data suggest that the human 5q11.2 breast cancer risk allele marked by rs889312 is mammary gland autonomous, and MIER3 is a candidate breast cancer susceptibility gene.

Human MCS5A1 candidate breast cancer susceptibility gene FBXO10 is induced by cellular stress and correlated with lens epithelium-derived growth factor (LEDGF).

Genetic variation and candidate genes associated with breast cancer susceptibility have been identified. Identifying molecular interactions between associated genetic variation and cellular proteins may help to better understand environmental risk. Human MCS5A1 breast cancer susceptibility associated SNP rs7042509 is located in F‐box protein 10 (FBXO10). An orthologous Rattus norvegicus DNA‐sequence that contains SNV ss262858675 is located in rat Mcs5a1, which is part of a mammary carcinoma susceptibility locus controlling tumor development in a non‐mammary cell‐autonomous manner via an immune cell‐mediated mechanism. Higher Fbxo10 expression in T cells is associated with Mcs5a increased susceptibility alleles. A common DNA–protein complex bound human and rat sequences containing MCS5A1/Mcs5a1 rs7042509/ss262858675 in electrophoretic mobility shift assays (EMSAs). Lens epithelium‐derived growth factor (LEDGF), a stress‐response protein, was identified as a candidate to bind both human and rat sequences using DNA‐pulldown and mass spectrometry. LEDGF binding was confirmed by LEDGF‐antibody EMSA and chromatin immunoprecipitation (ChIP). Ectopic expression of LEDGF/p75 increased luciferase activities of co‐transfected reporters containing both human and rat orthologs. Over‐expressed LEDGF/p75 increased endogenous FBXO10 mRNA levels in Jurkat cells, a human T‐cell line, implying LEDGF may be involved in increasing FBXO10 transcript levels. Oxidative and thermal stress of Jurkat cells increased FBXO10 and LEDGF expression, further supporting a hypothesis that LEDGF binds to a regulatory region of FBXO10 and increases expression during conditions favoring carcinogenesis. We conclude that FBXO10, a candidate breast cancer susceptibility associated gene, is induced by cellular stress and LEDGF may play a role in expression of this gene.

Physical confirmation and mapping of overlapping rat mammary carcinoma susceptibility QTLs, Mcs2 and Mcs6.

Only a portion of the estimated heritability of breast cancer susceptibility has been explained by individual loci. Comparative genetic approaches that first use an experimental organism to map susceptibility QTLs are unbiased methods to identify human orthologs to target in human population-based genetic association studies. Here, overlapping rat mammary carcinoma susceptibility (Mcs) predicted QTLs, Mcs6 and Mcs2, were physically confirmed and mapped to identify the human orthologous region. To physically confirm Mcs6 and Mcs2, congenic lines were established using the Wistar-Furth (WF) rat strain, which is susceptible to developing mammary carcinomas, as the recipient (genetic background) and either Wistar-Kyoto (WKy, Mcs6) or Copenhagen (COP, Mcs2), which are resistant, as donor strains. By comparing Mcs phenotypes of WF.WKy congenic lines with distinct segments of WKy chromosome 7 we physically confirmed and mapped Mcs6 to ∼33 Mb between markers D7Rat171 and gUwm64-3. The predicted Mcs2 QTL was also physically confirmed using segments of COP chromosome 7 introgressed into a susceptible WF background. The Mcs6 and Mcs2 overlapping genomic regions contain multiple annotated genes, but none have a clear or well established link to breast cancer susceptibility. Igf1 and Socs2 are two of multiple potential candidate genes in Mcs6. The human genomic region orthologous to rat Mcs6 is on chromosome 12 from base positions 71,270,266 to 105,502,699. This region has not shown a genome-wide significant association to breast cancer risk in pun studies of breast cancer susceptibility.

Mcs5c: a Mammary Carcinoma Susceptibility Locus Located in a Gene Desert that Associates with Tenascin C Expression

Genetic factors have been estimated to account for at least 30% of a womans risk to develop breast cancer. We have developed a rat model using Wistar Furth (WF) and Wistar Kyoto (WKy) strains to genetically identify mammary cancer susceptibility loci. The WKy allele of the mammary carcinogenesis susceptibility locus Mcs5c, was previously shown to reduce carcinoma multiplicity after 7,12-dimethylbenz-[a]anthracene (DMBA) exposure. In this study, Mcs5c was fine-mapped using WF.WKy congenic lines. Mcs5c was located to a region of approximately 176 kb on rat chromosome 5. One of the Mcs5c congenic lines containing a narrow Mcs5c WKy interval displayed a 40% decrease in average carcinoma number compared with WF-homozygous congenic controls after mammary carcinogenesis induction using two different models. As genetically mapped, the Mcs5c locus is located in a gene desert and thus is devoid of genes and annotated RNAs; thus, a genetic element in Mcs5c was hypothesized to regulate the expression of genes outside the locus. Tenascin c (Tnc) was identified as a candidate gene due to its reduced expression in thymus and ovarian tissues of Mcs5c WKy-homozygous congenic females compared with WF-homozygous congenic controls. This allele-specific differential expression is environmentally controlled. Cancer Prev Res; 4(1); 97–106. ©2011 AACR.

Differential 12‐O‐Tetradecanoylphorbol‐13‐acetate‐induced activation of rat mammary carcinoma susceptibility Fbxo10 variant promoters via a PKC‐AP1 pathway

Rat mammary carcinoma susceptibility 5a1 (Mcs5a1), which is concordant to human MCS5A1 breast cancer risk locus, mediates susceptibility by a non‐mammary cell‐autonomous mechanism associated with T cell differential expression of F‐box protein 10 (Fbxo10). Human FBXO10, an evolutionarily conserved ubiquitin ligase gene, was shown to have a potential role in regulating cell death by controlling the degradation of Bcl‐2, a key protein involved in apoptosis. Breast cancer susceptibility is controlled by interactions between environmental and genetic factors; therefore, we sought to determine if breast cancer risk‐associated environmental chemicals interact with Mcs5a1 variants using luciferase reporter constructs containing 4.2 kb Fbxo10 promoters based on alleles of mammary cancer susceptible Wistar Furth (WF) and resistant Wistar Kyoto (WKY) rat strains. 12‐O‐Tetradecanoylphorbol‐13‐acetate (TPA) induced activation of a 4.2 kb WF Fbxo10 promoter region, but lower levels of activation of the homologous WKY Fbxo10 promoter region. Using general and specific protein kinase inhibitors, we identified a protein kinase C (PKC) pathway that mediated TPA activation. We narrowed the possible PKCs to a member of the atypical PKC isoforms, namely PKCµ. We also determined that activator protein 1 (AP1) family member c‐Fos mediated TPA activation of the 4.2 kb WF Fbxo10 promoter. TPA was shown to induce endogenous FBXO10 mRNA and FBXO10 protein in Jurkat cells, a human T cell line, with a maximal level of expression from 1.5 to 2.5 h after exposure. These results indicate that FBXO10/Fbxo10 expression is regulated by a PKC‐dependent pathway acting through c‐Fos, which binds AP1‐specific DNA elements in Mcs5a1.

Significant overlap between human genome-wide association-study nominated breast cancer risk alleles and rat mammary cancer susceptibility loci

IntroductionHuman population-based genome-wide association (GWA) studies identify low penetrance breast cancer risk alleles; however, GWA studies alone do not definitively determine causative genes or mechanisms. Stringent genome- wide statistical significance level requirements, set to avoid false-positive associations, yield many false-negative associations. Laboratory rats (Rattus norvegicus) are useful to study many aspects of breast cancer, including genetic susceptibility. Several rat mammary cancer associated loci have been identified using genetic linkage and congenic strain based-approaches. Here, we sought to determine the amount of overlap between GWA study nominated human breast and rat mammary cancer susceptibility loci.MethodsWe queried published GWA studies to identify two groups of SNPs, one that reached genome-wide significance and one comprised of SNPs failing a validation step and not reaching genome- wide significance. Human genome locations of these SNPs were compared to known rat mammary carcinoma susceptibility loci to determine if risk alleles existed in both species. Rat genome regions not known to associate with mammary cancer risk were randomly selected as control regions.ResultsSignificantly more human breast cancer risk GWA study nominated SNPs mapped at orthologs of rat mammary cancer loci than to regions not known to contain rat mammary cancer loci. The rat genome was useful to predict associations that had met human genome-wide significance criteria and weaker associations that had not.ConclusionsIntegration of human and rat comparative genomics may be useful to parse out false-negative associations in GWA studies of breast cancer risk.