An epistatic interaction controls the latency of a transgene-induced mammary tumor

Roles of mitochondrial genetics in cancer metastasis

The contributions of mitochondria to cancer have been recognized for decades. However, the focus on the metabolic role of mitochondria and the diminutive size of the mitochondrial genome compared to the nuclear genome have hindered discovery of the roles of mitochondrial genetics in cancer. This review summarizes recent data demonstrating the contributions of mitochondrial DNA (mtDNA) copy-number variants (CNVs), somatic mutations, and germline polymorphisms to cancer initiation, progression, and metastasis. The goal is to summarize accumulating data to establish a framework for exploring the contributions of mtDNA to neoplasia and metastasis.

Transcriptome analysis of heterogeneity in mouse model of metastatic breast cancer

Background

Cancer metastasis is a complex process involving the spread of malignant cells from a primary tumor to distal organs. Understanding this cascade at a mechanistic level could provide critical new insights into the disease and potentially reveal new avenues for treatment. Transcriptome profiling of spontaneous cancer models is an attractive method to examine the dynamic changes accompanying tumor cell spread. However, such studies are complicated by the underlying heterogeneity of the cell types involved. The purpose of this study was to examine the transcriptomes of metastatic breast cancer cells using the well-established MMTV-PyMT mouse model.

Methods

Organ-derived metastatic cell lines were harvested from 10 female MMTV-PyMT mice. Cancer cells were isolated and sorted based on the expression of CD44^low/EpCAM^high or CD44^high/EpCAM^high surface markers. RNA from each cell line was extracted and sequenced using the NextSeq 500 Illumina platform. Tissue-specific genes were compared across the different metastatic and primary tumor samples. Reads were mapped to the mouse genome using STAR, and gene expression was quantified using RSEM. Single-cell RNA-seq (scRNA-seq) was performed on select samples using the ddSeq platform by BioRad and analyzed using Seurat v3.2.3. Monocle2 was used to infer pseudo-time progression.

Results

Comparison of RNA sequencing data across all cell populations produced distinct gene clusters. Differential gene expression patterns related to CD44 expression, organ tropism, and immunomodulatory signatures were observed. scRNA-seq identified expression profiles based on tissue-dependent niches and clonal heterogeneity. These cohorts of data were narrowed down to identify subsets of genes with high expression and known metastatic propensity. Dot plot analyses further revealed clusters expressing cancer stem cell and cancer dormancy markers. Changes in relevant genes were investigated across pseudo-time and tissue origin using Monocle2. These data revealed transcriptomes that may contribute to sub-clonal evolution and treatment evasion during cancer progression.

Conclusions

We performed a comprehensive transcriptome analysis of tumor heterogeneity and organ tropism during breast cancer metastasis. These data add to our understanding of metastatic progression and highlight targets for breast cancer treatment. These markers could also be used to image the impact of tumor heterogeneity on metastases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13058-021-01468-x.

Molecular evolutionary insights from PRLR in mammals

Prolactin (PRL) is a pleiotropic neurohormone secreted by the mammalian pituitary gland into the blood, thus reaching many tissues and organs beyond the brain. PRL binds to its receptor, PRLR, eliciting a molecular signaling cascade. This system modulates essential mammalian behaviors and promotes notable modifications in the reproductive female tissues and organs. Here, we explore how the intracellular domain of PRLR (PRLR-ICD) modulates the expression of the PRLR gene. Despite differences in the reproductive strategies between eutherian and metatherian mammals, there is no clear distinction between PRLR-ICD functional motifs. However, we found selection signatures that showed differences between groups, with many conserved functional elements strongly maintained through purifying selection across the class Mammalia. We observed a few residues under relaxed selection, the levels of which were more pronounced in Eutheria and particularly striking in primates (Simiiformes), which could represent a pre-adaptive genetic element protected from purifying selection. Alternative, new motifs, such as YLDP (318-321) and others with residues Y283 and Y290, may already be functional. These motifs would have been co-opted in primates as part of a complex genetic repertoire related to some derived adaptive phenotypes, but these changes would have no impact on the primordial functions that characterize the mammals as a whole and that are related to the PRL-PRLR system.

Roles of the mitochondrial genetics in cancer metastasis: not to be ignored any longer

Mitochondrial DNA (mtDNA) encodes for only a fraction of the proteins that are encoded within the nucleus, and therefore has typically been regarded as a lesser player in cancer biology and metastasis. Accumulating evidence, however, supports an increased role for mtDNA impacting tumor progression and metastatic susceptibility. Unfortunately, due to this delay, there is a dearth of data defining the relative contributions of specific mtDNA polymorphisms (SNP), which leads to an inability to effectively use these polymorphisms to guide and enhance therapeutic strategies and diagnosis. In addition, evidence also suggests that differences in mtDNA impact not only the cancer cells but also the cells within the surrounding tumor microenvironment, suggesting a broad encompassing role for mtDNA polymorphisms in regulating the disease progression. mtDNA may have profound implications in the regulation of cancer biology and metastasis. However, there are still great lengths to go to understand fully its contributions. Thus, herein, we discuss the recent advances in our understanding of mtDNA in cancer and metastasis, providing a framework for future functional validation and discovery.

Genetic variation in sensitivity to estrogens and breast cancer risk

Breast cancer risk is intimately intertwined with exposure to estrogens. While more than 160 breast cancer risk loci have been identified in humans, genetic interactions with estrogen exposure remain to be established. Strains of rodents exhibit striking differences in their responses to endogenous ovarian estrogens (primarily 17β-estradiol). Similar genetic variation has been observed for synthetic estrogen agonists (ethinyl estradiol) and environmental chemicals that mimic the actions of estrogens (xenoestrogens). This review of literature highlights the extent of variation in responses to estrogens among strains of rodents and compiles the genetic loci underlying pathogenic effects of excessive estrogen signaling. Genetic linkage studies have identified a total of the 35 quantitative trait loci (QTL) affecting responses to 17β-estradiol or diethylstilbestrol in five different tissues. However, the QTL appear to act in a tissue-specific manner with 9 QTL affecting the incidence or latency of mammary tumors induced by 17β-estradiol or diethylstilbestrol. Mammary gland development during puberty is also exquisitely sensitive to the actions of endogenous estrogens. Analysis of mammary ductal growth and branching in 43 strains of inbred mice identified 20 QTL. Regions in the human genome orthologous to the mammary development QTL harbor loci associated with breast cancer risk or mammographic density. The data demonstrate extensive genetic variation in regulation of estrogen signaling in rodent mammary tissues that alters susceptibility to tumors. Genetic variants in these pathways may identify a subset of women who are especially sensitive to either endogenous estrogens or environmental xenoestrogens and render them at increased risk of breast cancer.

Angiogenesis QTL on Mouse Chromosome 8 Colocalizes with Differential β-Defensin Expression

Identification of genetic factors that modify complex traits is often complicated by gene-environment interactions that contribute to the observed phenotype. In model systems, the phenotypic outcomes quantified are typically traits that maximize observed variance, which in turn, should maximize the detection of quantitative trait loci (QTL) in subsequent mapping studies. However, when the observed trait is dependent on multiple interacting factors, it can complicate genetic analysis, reducing the likelihood that the modifying mutation will ultimately be found. Alternatively, by focusing on intermediate phenotypes of a larger condition, we can reduce a model's complexity, which will, in turn, limit the number of QTL that contribute to variance. We used a novel method to follow angiogenesis in mice that reduces environmental variance by measuring endothelial cell growth from culture of isolated skin biopsies that varies depending on the genetic source of the tissue. This method, in combination with a backcross breeding strategy, is intended to reduce genetic complexity and limit the phenotypic effects to fewer modifier loci. We determined that our approach was an efficient means to generate recombinant progeny and used this cohort to map a novel s.c. angiogenesis QTL to proximal mouse chromosome (Chr.) 8 with suggestive QTL on Chr. 2 and 7. Global mRNA expression analysis of samples from parental reference strains revealed β-defensins as potential candidate genes for future study.

Genetic background may contribute to PAM50 gene expression breast cancer subtype assignments

Recent advances in genome wide transcriptional analysis have provided greater insights into the etiology and heterogeneity of breast cancer. Molecular signatures have been developed that stratify the conventional estrogen receptor positive or negative categories into subtypes that are associated with differing clinical outcomes. It is thought that the expression patterns of the molecular subtypes primarily reflect cell-of-origin or tumor driver mutations. In this study however, using a genetically engineered mouse mammary tumor model we demonstrate that the PAM50 subtype signature of tumors driven by a common oncogenic event can be significantly influenced by the genetic background on which the tumor arises. These results have important implications for interpretation of "snapshot" expression profiles, as well as suggesting that incorporation of genetic background effects may allow investigation into phenotypes not initially anticipated in individual mouse models of cancer.

Epistatic Control of Mammary Cancer Susceptibility in Mice may Depend on the Dietary Environment

Recent studies have linked a high fat diet to the development of breast cancer, but any genetic basis for this association is poorly understood. We investigated this association with an epistatic analysis of seven cancer traits in a segregating population of mice with metastatic mammary cancer that were fed either a control or a high-fat diet. We used an interval mapping approach with single nucleotide polymorphisms to scan all 19 autosomes, and discovered a number of diet-independent epistatic interactions of quantitative trait loci (QTLs) affecting these traits. More importantly, we also discovered significant epistatic by diet interactions affecting some of the traits that suggested these epistatic effects varied depending on the dietary environment. An analysis of these interactions showed some were due to epistasis that occurred in mice fed only the control diet or only the high-fat diet whereas other interactions were generated by differential effects of epistasis in the two dietary environments. Some of the epistatic QTLs appeared to colocalize with cancer QTLs mapped in other mouse populations and with candidate genes identified from eQTLs previously mapped in this population, but others represented novel modifying loci affecting these cancer traits. It was concluded that these diet-dependent epistatic QTLs contribute to a genetic susceptibility of dietary effects on breast cancer, and their identification may eventually lead to a better understanding that will be needed for the design of more effective treatments for this disease.

Bromodomain-Containing Protein 4: A Dynamic Regulator of Breast Cancer Metastasis through Modulation of the Extracellular Matrix

Metastasis is an extremely complex process that accounts for most cancer-related deaths. Malignant primary tumors can be removed surgically, but the cells that migrate, invade, and proliferate at distant organs are often the cells that prove most difficult to target therapeutically. There is growing evidence that host factors outside of the primary tumors are of major importance in the development of metastasis. Recently, we have shown that the bromodomain-containing protein 4 or bromodomain 4 (Brd4) functions as an inherited susceptibility gene for breast cancer progression and metastasis. In this paper, we will discuss that host genetic background on which a tumor arises can significantly alter the biology of the subsequent metastatic disease, and we will focus on the role of Brd4 in regulating metastasis susceptibility.

Genetic heterogeneity of skin microvasculature

Angiogenesis, the formation of new blood vessels from existing vasculature, is a complex process that is essential for normal embryonic development. Current models for experimental evaluation of angiogenesis often use tissue from large vessels like the aorta and umbilical vein, which are phenotypically distinct from microvasculature. We demonstrate that the utilization of skin to measure microvascular angiogenesis in embryonic and adult tissues is an efficient way to quantify microvasculature angiogenesis. We validate this approach and demonstrate its added value by showing significant differences in angiogenesis in monogenic and polygenic mouse models. We discovered that the pattern of angiogenic response among inbred mouse strains in this ex vivo assay differs from the strain distributions of previous in vivo angiogenesis assays. The difference between the ex vivo and in vivo assays may be related to systemic factors present in whole animals. Expression analysis of cultured skin biopsies from strains of mice with opposing angiogenic response was performed to identify pathways that contribute to differential angiogenic response. Increased expression of negative regulators of angiogenesis in C57Bl/6J mice was associated with lower growth rates.

Lessons in signaling and tumorigenesis from polyomavirus middle T antigen

The small DNA tumor viruses have provided a very long-lived source of insights into many aspects of the life cycle of eukaryotic cells. In recent years, the emphasis has been on cancer-related signaling. Here we review murine polyomavirus middle T antigen, its mechanisms, and its downstream pathways of transformation. We concentrate on the MMTV-PyMT transgenic mouse, one of the most studied models of breast cancer, which permits the examination of in situ tumor progression from hyperplasia to metastasis.

Gene expression profiles and breast cancer metastasis: a genetic perspective

The majority of cancer mortality is attributed to metastasis, which is the spread of tumor cells to a secondary site. Several studies have demonstrated that the genetic background on which a tumor arises has a major effect on both metastatic efficiency and on predictive gene expression profiles. These observations suggest that there is variability in metastasis frequency between individuals and that some individuals could be more prone to secondary tumor formation and development than others. Thus, genetic background might have important clinical implications in metastasis detection, management and prevention.

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(alpha)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.

Mouse modifier genes in mammary tumorigenesis and metastasis

Tumorigenesis and metastasis are complex multistep processes. In addition to the numerous somatic mutations that facilitate cancer progression, there is abundant evidence that an individual's genetic background not only contributes to overall cancer risk, but also specifically influences metastatic potential. The handful of human susceptibility genes that have been identified thus far do not fully account for hereditary cancer risk, and the discovery of additional susceptibility loci using population based studies is complex, time-consuming and expensive. Therefore, we and others have used a variety of mouse models to identify novel candidate susceptibility genes. Here we review how these mouse models have contributed to our understanding of the role of genetic background in modifying tumorigenesis and metastasis susceptibility.

Genotype X diet interactions in mice predisposed to mammary cancer: II. Tumors and metastasis

High dietary fat intake and obesity may increase the risk of susceptibility to certain forms of cancer. To study the interactions of dietary fat, obesity, and metastatic mammary cancer, we created a population of F(2) mice cosegregating obesity QTL and the MMTV-PyMT transgene. We fed the F(2) mice either a very high-fat or a matched-control-fat diet, and we measured growth, body composition, age at mammary tumor onset, tumor number and severity, and formation of pulmonary metastases. SNP genotyping across the genome facilitated analyses of QTL and QTL x diet interaction effects. Here we describe effects of diet on mammary tumor and metastases phenotypes, mapping of tumor/metastasis modifier genes, and the interaction between dietary fat levels and effects of cancer modifiers. Results demonstrate that animals fed a high-fat diet are not only more likely to experience decreased mammary cancer latency but increased tumor growth and pulmonary metastases occurrence over an equivalent time. We identified 25 modifier loci for mammary cancer and pulmonary metastasis, likely representing 13 unique loci after accounting for pleiotropy, and novel QTL x diet interactions at a majority of these loci. These findings highlight the importance of accurately modeling not only the human cancer characteristics in mice but also the environmental exposures of human populations.

Genotype X diet interactions in mice predisposed to mammary cancer. I. Body weight and fat

High dietary fat intake and obesity may increase susceptibility to certain forms of cancer. To study the interactions of dietary fat, obesity, and metastatic mammary cancer, we created a population of F(2) mice cosegregating obesity QTL and the MMTV-PyMT transgene. We fed the F(2) mice either a very-high-fat or a matched-control-fat diet and measured growth, body composition, age at mammary tumor onset, tumor number and severity, and formation of pulmonary metastases. SNP genotyping across the genome facilitated analyses of QTL and QTL x diet interaction effects. Here we describe development of the F(2) population (n = 615) which resulted from a cross between the polygenic obesity model M16i and FVB/NJ-TgN (MMTV-PyMT)(634Mul), effects of diet on growth and body composition, and QTL and QTL x diet and/or gender interaction effects for growth and obesity-related phenotypes. We identified 38 QTL for body composition traits that were significant at the genome-wide 0.05 level, likely representing nine distinct loci after accounting for pleiotropic effects. QTL x diet and/or gender interactions were present at 15 of these QTL, indicating that such interactions play a significant role in defining the genetic architecture of complex traits such as body weight and obesity.

Identification of novel modifier loci of Apc Min affecting mammary tumor development

Genetic background affects the susceptibility to mammary tumor development in Apc(Min/+) mice. Here we report the identification of four novel modifier loci that influence different aspects of mammary tumor development in Apc(Min/+) mice. Analysis of tumor development in a backcross of (FVBB6 Apc(Min/+)) x B6 Apc(Min/+) mice has identified a modifier on chromosome 9 that significantly affects tumor multiplicity, and a modifier on chromosome 4 that significantly affects tumor latency and affects tumor number with suggestive significance. This modifier was also identified in a backcross involving 129X1/SvJ and B6 Apc(Min/+) mice. A modifier on chromosome 18 specifically affects tumor latency but not tumor number. Kaplan-Meier analysis suggests there is at least an additive interaction affecting tumor latency between the loci on chromosomes 4 and 18. We also identified a modifier locus on chromosome 6 that interacts with the loci on chromosome 4 and chromosome 9 to affect tumor number. These results suggest that multiple genetic loci control different aspects of mammary tumor development. None of these modifiers is associated with intestinal tumor susceptibility, which indicates that these modifiers act on tumor development in a tissue-specific manner.

Characterization of breast cancer by array comparative genomic hybridization

Cancer progression is due to the accumulation of recurrent genomic alterations that induce growth advantage and clonal expansion. Most of these genomic changes can be detected using the array comparative genomic hybridization (CGH) technique. The accurate classification of these genomic alterations is expected to have an important impact on translational and basic research. Here we review recent advances in CGH technology used in the characterization of different features of breast cancer. First, we present bioinformatics methods that have been developed for the analysis of CGH arrays; next, we discuss the use of array CGH technology to classify tumor stages and to identify and stratify subgroups of patients with different prognoses and clinical behaviors. We finish our review with a discussion of how CGH arrays are being used to identify oncogenes, tumor suppressor genes, and breast cancer susceptibility genes.

Genetic mapping in mice identifies DMBT1 as a candidate modifier of mammary tumors and breast cancer risk

Low-penetrance breast cancer susceptibility alleles seem to play a significant role in breast cancer risk but are difficult to identify in human cohorts. A genetic screen of 176 N2 backcross progeny of two Trp53(+/-) strains, BALB/c and C57BL/6, which differ in their susceptibility to mammary tumors, identified a modifier of mammary tumor susceptibility in an approximately 25-Mb interval on mouse chromosome 7 (designated SuprMam1). Relative to heterozygotes, homozygosity for BALB/c alleles of SuprMam1 significantly decreased mammary tumor latency from 70.7 to 61.1 weeks and increased risk twofold (P = 0.002). Dmbt1 (deleted in malignant brain tumors 1) was identified as a candidate modifier gene within the SuprMam1 interval because it was differentially expressed in mammary tissues from BALB/c-Trp53(+/-) and C57BL/6-Trp53(+/-) mice. Dmbt1 mRNA and protein was reduced in mammary glands of the susceptible BALB/c mice. Immunohistochemical staining demonstrated that DMBT1 protein expression was also significantly reduced in normal breast tissue from women with breast cancer (staining score, 1.8; n = 46) compared with cancer-free controls (staining score, 3.9; n = 53; P < 0.0001). These experiments demonstrate the use of Trp53(+/-) mice as a sensitized background to screen for low-penetrance modifiers of cancer. The results identify a novel mammary tumor susceptibility locus in mice and support a role for DMBT1 in suppression of mammary tumors in both mice and women.

Effects of FVB/NJ and C57Bl/6J strain backgrounds on mammary tumor phenotype in inducible nitric oxide synthase deficient mice

The ability to genetically manipulate mice has led to rapid progress in our understanding of the roles of different gene products in human disease. Transgenic mice have often been created in the FVB/NJ (FVB) strain due to its high fecundity, while gene-targeted mice have been developed in the 129/SvJ-C57Bl/6J strains due to the capacity of 129/SvJ embryonic stem cells to facilitate germline transmission. Gene-targeted mice are commonly backcrossed into the C57Bl/6J (B6) background for comparison with existing data. Genetic modifiers have been shown to modulate mammary tumor latency in mouse models of breast cancer and it is commonly known that the FVB strain is susceptible to mammary tumors while the B6 strain is more resistant. Since gene-targeted mice in the B6 background are frequently bred into the polyomavirus middle T (PyMT) mouse model of breast cancer in the FVB strain, we have sought to understand the impact of the different genetic backgrounds on the resulting phenotype. We bred mice deficient in the inducible nitric oxide synthase (iNOS) until they were congenic in the PyMT model in the FVB and B6 strains. Our results reveal that the large difference in mean tumor latencies in the two backgrounds of 53 and 92 days respectively affect the ability to discern smaller differences in latency due to the Nos2 genetic mutation. Furthermore, the longer latency in the B6 strain enables a more detailed analysis of tumor formation indicating that individual tumor development is not stoichastic, but is initiated in the #1 glands and proceeds in early and late phases. NO production affects tumors that develop early suggesting an association of iNOS-induced NO with a more aggressive tumor phenotype, consistent with human clinical data positively correlating iNOS expression with breast cancer progression. An examination of lung metastases, which are significantly reduced in PyMT/iNOS^−/− mice compared with PyMT/iNOS^+/+ mice only in the B6 background, is concordant with these findings. Our data suggest that PyMT in the B6 background provides a useful model for the study of inflammation-induced breast cancer.

Parallel analysis of transcript and translation profiles: identification of metastasis-related signal pathways differentially regulated by drug and genetic modifications

Tumor metastasis is a complex multistep process normally involving dysregulation of multiple signal transduction pathways. In this study, we developed a novel approach to efficiently define dysreguated pathways associated with metastasis by comparing global gene and protein expressions of two distinct metastasis-suppressed models. Consequently, we identified common features shared by the two models which are potentially associated with metastasis. The efficiency of metastasis from the highly aggressive polyoma middle T-induced mouse mammary tumors was suppressed by either prolonged caffeine exposure or by breeding the animal to a low metastatic mouse strain. Molecular profiles of the primary tumors from both metastasis-suppressed classes were then derived to identify molecules and pathways that might underlie a common mechanism of metastasis. A number of differentially regulated genes and proteins were identified, including genes encoding basement membrane components, which were inversely related to metastatic efficiency. In addition, the analysis revealed that the Stat signal transduction pathways were potentially associated with metastasis inhibition, as demonstrated by enhanced Stat1 activation, and decreased Stat5 phosphorylation in both genetic and pharmacological modification models. Tumor cells of low-metastatic genotypes also demonstrated anti-apoptotic properties. The common changes of these pathways in all of the metastasis-suppressed systems suggest that they may be critical components in the metastatic cascade, at least in this model system. Our data demonstrate that analysis of common changes in genes and proteins in a metastatic-related context greatly decrease the complexity of data analysis, and may serve as a screening tool to identify biological important factors from large scale data.

Multilocus inheritance determines predisposition to α-radiation induced bone tumourigenesis in mice

In a recent study, we presented evidence for genetic predisposition governing radiation osteosarcomagenesis in mice. Following the incorporation of the bone-seeking alpha emitter 227Th, approximately 25% of the variance in osteosarcoma incidence was determined by inherited genetic factors. We have now mapped 5 susceptibility loci in crosses between the more susceptible BALB/c and the more resistant CBA/Ca strains. The major QTL on chromosome 14 overlaps with a locus that was already found in our previous study, using different strains of mice. Here, we investigate the effect by which the major susceptibility locus and 4 minor modifier loci interact to influence osteosarcoma predisposition. Following incorporation of the bone-seeking isotope, 100% of mice that harbour high-risk genotypes at all 5 susceptibility loci develop osteosarcoma with an average of 472 days latency times. In 10 mice inheriting exclusively low-risk genotypes only 1 osteosarcoma was found, arising after 733 days latency time. Inheritance of distinct combinations of BALB/c and CBA/Ca alleles at the susceptibility loci confer more extreme phenotypes in terms of susceptibility or resistance than observed in either of the two parental inbred strains. From the present study, we demonstrate that additive effects of multiple alleles, each making only a minor phenotypic contribution, can combine and significantly alter tumour risk. This mechanism can be of particular importance in genetically heterogeneous populations such as man.

Genetic identification of distinct loci controlling mammary tumor multiplicity, latency, and aggressiveness in the rat

The rat is considered an excellent model for studying human breast cancer. Therefore, understanding the genetic basis of susceptibility to mammary cancer in this species is of great interest. Previous studies based on crosses involving the susceptible strain WF (crossed with the resistant strains COP or WKY) and focusing on tumor multiplicity as the susceptibility phenotype led to the identification of several loci that control chemically induced mammary cancer. The present study was aimed to determine whether other loci can be identified by analyzing crosses derived from another susceptible strain on the one hand, and by including phenotypes other than tumor multiplicity on the other hand. A backcross was generated between the susceptible SPRD-Cu3 strain and the resistant WKY strain. Female progeny were genotyped with microsatellite markers covering all rat autosomes, treated with a single dose of DMBA, and phenotyped with respect to tumor latency, tumor multiplicity, and tumor aggressiveness. Seven loci controlling mammary tumor development were detected. Different loci control tumor multiplicity, latency, and aggressiveness. While some of these loci colocalize with loci identified in crosses involving the susceptible strain WF, new loci have been uncovered, indicating that the use of distinct susceptible and resistant strain pairs will help in establishing a comprehensive inventory of mammary cancer susceptibility loci.

Syngeneic mouse mammary carcinoma cell lines: two closely related cell lines with divergent metastatic behavior

Two cell lines, Met-1(fvb2) and DB-7(fvb2), with different metastatic potential, were derived from mammary carcinomas in FVB/N-Tg(MMTV-PyVmT) and FVB/N-Tg(MMTV-PyVmT ( Y315F/Y322F )) mice, transplanted into syngeneic FVB/N hosts and characterized. The lines maintain a stable morphological and biological phenotype after multiple rounds of in vitro culture and in vivo transplantation. The Met-1(fvb2) line derived from a FVB/N-Tg(MMTV-PyVmT) tumor exhibits invasive growth and 100% metastases when transplanted into the females FVB/N mammary fat pad. The DB-7(fvb2) line derived from the FVB/N-Tg(MMTV-PyVmT ( Y315F/Y322F )) with a "double base" modification at Y315F/Y322F exhibits more rapid growth when transplanted into the mammary fat pad, but a lower rate of metastasis (17%). The Met1(fvb2) cells show high activation of AKT, while DB-7(fvb2) cells show very low levels of AKT activation. The DNA content and gene expression levels of both cell lines are stable over multiple generations. Therefore, these two cell lines provide a stable, reproducible resource for the study of metastasis modulators, AKT molecular pathway interactions, and gene target and marker discovery.

Sipa1 is a candidate for underlying the metastasis efficiency modifier locus Mtes1

We previously identified loci in the mouse genome that substantially influence the metastatic efficiency of mammary tumors. Here, we present data supporting the idea that the signal transduction molecule, Sipa1, is a candidate for underlying the metastasis efficiency modifier locus Mtes1. Analysis of candidate genes identified a nonsynonymous amino acid polymorphism in Sipa1 that affects the Sipa1 Rap-GAP function. Spontaneous metastasis assays using cells ectopically expressing Sipa1 or cells with knocked-down Sipa1 expression showed that metastatic capacity was correlated with cellular Sipa1 levels. We examined human expression data and found that they were consistent with the idea that Sipa1 concentration has a role in metastasis. Taken together, these data suggest that the Sipa1 polymorphism is one of the genetic polymorphisms underlying the Mtes1 locus. This report is also the first demonstration, to our knowledge, of a constitutional genetic polymorphism affecting tumor metastasis.

Modifiers of mammary tumor progression and metastasis on mouse chromosomes 7, 9, and 17

Tumor progression, the growth and dissemination of primary tumor to secondary sites, is of critical clinical importance since the vast majority of patients succumb to metastatic disease rather than to the primary tumor. Many factors are likely to influence this process, including the primary oncogenic events, environmental exposures and stress and progressive stochastic mutations. Previously, our laboratory demonstrated that an additional factor, the genetic background on which tumors arose, had a significant effect on metastatic efficiency. Using a highly metastatic transgene-induced mammary tumor model, a locus modulating metastatic efficiency, Mtes1, was localized on proximal mouse Chromosome 19. In addition, a number of additional suggestive loci were observed on several other chromosomes. To confirm the presence of these additional loci before initiating cloning strategies, chromosomal substitution strains have been constructed and assayed for modification of the cancer phenotypes. Using the chromosomal substitution strains, an additional modifier modulating tumor latency was confirmed, as well as three new modifier genes that alter the kinetics of tumor progression. Identification and analysis of these loci will likely present interesting and novel information about cancer heterogeneity in the human population.

Caffeine suppresses metastasis in a transgenic mouse model: a prototype molecule for prophylaxis of metastasis

A significant fraction of cancer patients have occult disseminated tumors at the time of primary diagnosis, which usually progress to become clinically relevant lesions. Since the majority of cancer mortality is associated with metastatic disease, the ability to inhibit the growth of the secondary tumors would significantly reduce cancer-related morbidity and mortality. We have investigated whether caffeine, which has been shown to suppress tumor cell invasiveness and experimental metastasis, can suppress metastasis in a spontaneous transgene-induced mammary tumor model. Chronic exposure to caffeine prior to the appearance of palpable mammary tumors significantly reduced both tumor burden and metastatic colonization. However, when caffeine exposure began after the appearance of frank tumors, caffeine suppressed metastasis without changing primary tumor burden. The means by which caffeine suppressed metastatic activity may be associated with inhibition of malignant transformation of mammary epithelial cells, inhibition of conversion of dormant tumor cells to micrometastases, micrometastases to macrometastases, or inhibition of tumor cell adhesion and motility. Gene and protein expression patterns resulting from caffeine treatment showed that metastasis suppression may be associated with up-regulation the mRNA expression of multiple extracellular matrix genes, including Fbln1, Bgn, Sparc, Fbn1, Loxl1, Colla1, Col3a1, Col5a1, ColS5a2, ColSa3, Col6a1, Col6a2, and Col6a3. These data suggested that caffeine or other methyl xanthine derivatives may improve the clinical outcome in patients prior to and following the diagnosis of metastatic disease, and could potentially reduce the morbidity and mortality associated with disseminated tumors.

Ron receptor signaling augments mammary tumor formation and metastasis in a murine model of breast cancer

The tyrosine kinase receptor Ron has been implicated in several types of cancer, including overexpression in human breast cancer. This is the first report describing the effect of Ron signaling on tumorigenesis and metastasis in a mouse model of breast cancer. Mice with a targeted deletion of the Ron tyrosine kinase signaling domain (TK-/-) were crossed to mice expressing the polyoma virus middle T antigen (pMT) under the control of the mouse mammary tumor virus promoter. Both pMT-expressing wild-type control (pMT+/- TK+/+) and pMT+/- TK-/- mice developed mammary tumors and lung metastases. However, a significant decrease in mammary tumor initiation and growth was found in the pMT+/- TK-/- mice compared with controls. An examination of mammary tumors showed that there was a significant decrease in microvessel density, significantly decreased cellular proliferation, and a significant increase in terminal deoxynucleotidyl transferase-mediated nick end labeling-positive staining in mammary tumor cells from the pMT+/- TK-/- mice compared with the pMT+/- TK+/+ mice. Biochemical analyses on mammary tumor lysates showed that whereas both the pMT-expressing TK+/+ and TK-/- tumors have increased Ron expression compared with normal mammary glands, the pMT-expressing TK-/- tumors have deficits in mitogen-activated protein kinase and AKT activation. These results indicate that Ron signaling synergizes with pMT signaling to induce mammary tumor formation, growth, and metastasis. This effect may be mediated in part through the regulation of angiogenesis and through proliferative and cell survival pathways regulated by mitogen-activated protein kinase and AKT.

Genetic loci that control the angiogenic response to basic fibroblast growth factor

Angiogenesis is controlled by a balance between stimulatory growth factors and endogenous inhibitors. We propose that the balance of stimulators and inhibitors, as well as the general sensitivity of the endothelium to these factors, varies from individual to individual. Indeed, we have found that individual mouse strains have dramatically different responses to growth factor-induced neovascularization. Quantitative trait loci (QTLs), which influence the extent of angiogenesis induced by vascular endothelial growth factor (VEGF), were previously identified by our laboratory. Since genetic susceptibility may vary according to the angiogenic stimulator, we have undertaken a similar mapping approach to identify QTLs that influence basic fibroblast growth factor (FGF2) induced neovascularization in the BXD series of recombinant inbred mouse strains. Composite and multiple interval mapping identified areas of chromosomes 4, 13, 15, and 18. These new angiogenesis QTLs, named AngFq1 through AngFq4 (for angiogenesis due to FGF2), are different from previously identified VEGF QTLs. The mapped regions contain several genes involved in the angiogenic process including matrix metalloproteinase 16, eph receptor A7, angiopoetin 1, endothelial lipase, and autotaxin. Differences in these regions may influence individual susceptibility to angiogenesis related diseases such as cancer, macular degeneration, atherosclerosis, and arthritis.

Bax regulates c-Myc-induced mammary tumour apoptosis but not proliferation in MMTV-c-myc transgenic mice

The expression of the proto-oncogene c-myc is frequently deregulated, via multiple mechanisms, in human breast cancers. Deregulated expression of c-myc contributes to mammary epithelial cell transformation and is causally involved in mammary tumorigenesis in MMTV-c-myc transgenic mice. c-Myc is known to promote cellular proliferation, apoptosis, genomic instability and tumorigenesis in several distinct tissues, both in vivo and in vitro. Expression of the proapoptotic regulatory gene bax is reduced or absent in human breast cancers, and c-Myc has been shown to regulate the expression of Bax, as well as cooperate with Bax in controlling apoptosis in a fibroblast model. Additionally, loss of bax reduces c-Myc-induced apoptosis in lymphoid cells and increases c-Myc-mediated lymphomagenesis in vivo. In order to assess whether loss of bax could influence c-Myc-induced apoptosis and tumorigenesis in the mammary gland in vivo, we generated MMTV-c-myc transgenic mice in which neither, one, or both wild-type alleles of bax were eliminated. Haploid loss of bax in MMTV-c-myc transgenic mice resulted in significantly reduced mammary tumour apoptosis. As anticipated for an apoptosis-regulatory gene, loss of the wild-type bax alleles did not significantly alter cellular proliferation in either mammary adenocarcinomas or dysplastic mammary tissues. However, in contrast to c-Myc-mediated lymphomagenesis, loss of one or both alleles of bax in MMTV-c-myc transgenic mice did not significantly enhance mammary tumorigenesis, despite evidence that haploid loss of bax might modestly increase mammary tumour multiplicity. Our results demonstrate that Bax contributes significantly to c-Myc-induced apoptosis in mammary tumours. In addition, they suggest that in contrast to c-Myc-induced lymphomagenesis, mammary tumorigenesis induced by deregulated c-myc expression requires some amount of Bax expression.

Global expression profiling identifies signatures of tumor virulence in MMTV-PyMT-transgenic mice: correlation to human disease

FVB/N-Tg (MMTV-PyMT)(634Mul)-transgenic mice develop multifocal mammary tumors with a high incidence of pulmonary metastasis. We have demonstrated previously that mammary tumors derived from transgene-positive F1 progeny in particular inbred strains display altered latency, tumor growth rates, and metastatic rates when compared with the FVB/NJ homozygous parent. To identify genes with expression that might be critical in modifying the biological behavior of MMTV-PyMT tumors, we performed a detailed comparative analysis of expression profiles from mammary tumors arising in the parental FVB/NJ background and F1 progeny from crosses with I/LnJ, LP/J, MOLF/Ei, and NZB/B1NJ mice. Compared with normal mammary glands, gene expression profiles of tumors from all five strains exhibited up-regulation of genes involved in cell growth (e.g., Cks1 and CDC25C) and down-regulation of cell adhesion molecules, with many genes associated previously with human breast cancer such as STAT2, CD24 antigen, gelsolin, and lipocalin2. To identify genes with significant variation in expression between the five different genotypes, significance analysis of microarrays (SAM) and one-way ANOVA were used. Three definable groupings of tumors were identified: (a) tumors derived in the LP/J F1 and MOLF/Ei F1 strains in which tumor growth and dissemination are suppressed and latency prolonged; (b) the most aggressive tumors from the FVB/NJ parental strain and I/LnJ F1 genomic backgrounds; and (c) an intermediate virulence phenotype with tumors from NZB/B1NJ-F1 crosses. These array based assessments correlated well with a composite phenotype ranking using a "virulence" index. The gene expression signature that is associated with a high metastatic rate in the mouse contains the same 17 genes described recently as the signature gene set predictive of metastasis in human tumors (1) with 16 of the 17 genes exhibiting the same directional change in expression associated with human metastases. These results demonstrate that the genetic analysis of mouse models of tumorigenesis may be highly relevant to human cancer and that the metastatic phenotype of a tumor may be affected by the germline genetic configuration of the host.

Cancer susceptibility in the mouse: genetics, biology and implications for human cancer

Growing evidence that a large proportion of apparently non-hereditary sporadic cancers occur in genetically predisposed individuals has emphasized the need to identify the underlying susceptibility genes. Increasingly, it seems that the best approach to define the numerous genes that have small but cumulative effects is to first identify and map them in mice, and subsequently to study the role of their homologues in humans. Development of new gene-mapping resources and strategies in mice has, for the first time, allowed some of these genes to be identified. In future, this unique approach is likely to provide important insights into the pathways of tumour development and might ultimately lead to more effective individually targeted cancer-prevention strategies.

Mammary tumor latency is increased in mice lacking the inducible nitric oxide synthase

Nitric oxide (NO) and its metabolites are implicated in carcinogenesis and metastasis. Both stimulatory and inhibitory effects of NO have been reported in relation to breast cancer and its role in the development of malignancies and metastasis remains uncertain. We have used the polyomavirus middle T antigen (PyV-mT) targeted to the mouse mammary gland and bred into an inducible NO synthase (iNOS)-deficient C57Bl/6 strain to examine a role for nitric oxide in modulating tumors that develop in the complex environment of the whole animal. The development of hyperplasias was delayed to the extent that the earliest palpable tumors arose 2-4 weeks later in PyV-mT/iNOS(-/-) mice compared with PyV-mT/iNOS(+/+) mice, identifying a role for iNOS in early events in mammary tumor formation. Tumors that did develop in PyV-mT/iNOS(-/-) mice were characteristically well differentiated and had a cribriform pattern. Other tumors were myoepithelial adenocarcinomas with uniform nuclear size. In contrast, mice capable of iNOS activity typically developed solid nodular adenocarcinomas with a high mitotic index and pleomorphic nuclei. No significant effect of iNOS deficiency was found on vascular density in hyperplasias or tumors by examining CD31-positive vessels. The infiltration of lesions by macrophages, cells capable of significant NO production, remained unchanged in PyV-mT/iNOS(-/-) mice. Metastatic potential was retained by PyV-mT-transformed epithelium in the absence of iNOS, indicating that NO production by iNOS is not essential for this process. These results indicate a role for iNOS in tumorigenesis, particularly in the regulation of early events.

Novel murine mammary epithelial cell lines that form osteolytic bone metastases: effect of strain background on tumor homing

We have developed a series of novel mammary epithelial cell lines from tumors arising in strain 129 mice, with the ultimate goal of evaluating the role of host factors in the development of bone metastases. Mammary tumors were induced in mice with subcutaneously implanted medroxyprogesterone acetate (MPA) pellets followed by administration of DMBA by oral gavage. Mammary tumor development was efficient in the 129 strain and was independent of osteopontin (OPN) expression. Epithelial cell lines were isolated from these tumors; surprisingly, these cells did not form tumors upon inoculation into the mammary fat pad of syngeneic mice, even when MPA was present. One OPN-deficient cell line was selected for further study; full transformation of these cells required expression of both polyoma middle T and activated ras. These doubly transfected cells, 1029 GP+Er3, grew in soft agar, and formed hormone-independent tumors efficiently in the mammary fat pad that spontaneously metastasized to several soft tissue sites but not to the bone. Derivatives of these cells were isolated from tumors arising in the fat pad and from a lung metastasis (r3T and r3L, respectively): these cells formed tumors more rapidly in the fat pad than the parental GP+Er3 cells. Upon left ventricle injection, the r3T and r3L cells formed osteolytic bone metastases in 129 mice, with few metastases seen in other organs. These tumors filled the marrow cavity, and caused extensive destruction of both cortical and trabecular bone. Intriguingly, in an alternative syngeneic host, (129xC57B1/6) F1, osteolytic bone metastases were not seen on x-ray; instead extensive liver metastasis was present in these mice, indicating that genetic factors in these two strains regulate tumor cell homing and distribution during metastasis. These cell lines provide an important new tool in the study of bone metastasis, particularly in elucidating the role of host factors in the development of these lesions, as the 129 mouse strain is frequently used for genetic manipulations in the mouse.

A bioinformatics-based strategy identifies c-Myc and Cdc25A as candidates for the Apmt mammary tumor latency modifiers

The epistatically interacting modifier loci (Apmt1 and Apmt2) accelerate the polyoma Middle-T (PyVT)-induced mammary tumor. To identify potential candidate genes loci, a combined bioinformatics and genomics strategy was used. On the basis of the assumption that the loci were functioning in the same or intersecting pathways, a search of the literature databases was performed to identify molecular pathways containing genes from both candidate intervals. Among the genes identified by this method were the cell cycle-associated genes Cdc25A and c-Myc, both of which have been implicated in breast cancer. Genomic sequencing revealed noncoding polymorphism in both genes, in the promoter region of Cdc25A, and in the 3' UTR of c-Myc. Molecular and in vitro analysis showed that the polymorphisms were functionally significant. In vivo analysis was performed by generating compound PyVT/Myc double-transgenic animals to mimic the hypothetical model, and was found to recapitulate the age-of-onset phenotype. These data suggest that c-Myc and Cdc25A are Apmt1 and Apmt2, and suggest that, at least in certain instances, bioinformatics can be utilized to bypass congenic construction and subsequent mapping in conventional QTL studies.

Complexities of cancer research: mouse genetic models

Cancer susceptibility is a complex interaction of an individual's genetic composition and environmental exposures. Huge strides have been made in understanding cancer over the past 100 yr, from recognition of cancer as a genetic disease, to identification of specific carcinogens, isolation of oncogenes, and recognition of tumor suppressors. A tremendous amount of knowledge has accumulated about the etiology of cancer. Cancer genetics has played a significant role in these discoveries. Analysis of high-risk familial cancers has led to the discovery of new tumor suppressor genes and important cancer pathways. These families, however, represent only a small fraction of cancer in the general population. Most cancer is instead probably the result of an intricate interaction of polymorphic susceptibility genes with the sea of environmental exposures that humans experience. Although the central cadre of cancer genes is known, little is understood about the peripheral genes that likely comprise the polymorphic susceptibility loci. The challenge for cancer genetics is therefore to move forward from the mendelian genetics of the rare familial cancer syndromes into the field of quantitative trait loci, susceptibility factors, and modifier genes. By identifying the genes that modulate an individual's susceptibility to cancer after an environmental exposure, researchers will be able to gain important insights into human biology, cancer prevention, and cancer treatment. This article summarizes the current state of quantitative trait genetic analysis and the tools, both proven and theoretical, that may be used to unravel one of the great challenges in cancer genetics.

Recent advances in inducible expression in transgenic mice

In order to accurately analyze gene function in transgenic mice, as well as to generate credible murine models of human diseases, the ability to regulate temporal- and spatial-specific expression of target genes is absolutely critical. Pioneering work in inducible transgenics, begun in the 1980s and continuing to the present, has led to the development of a variety of different inducible systems dedicated to this goal, the shared basis of which is the accurate conditional expression of a given transgene. Recent advances in inducible transgene expression in mice are discussed.

Three loci modify growth of a transgene-induced mammary tumor: suppression of proliferation associated with decreased microvessel density

In earlier studies it was observed that the genetic background significantly affected the phenotype of a transgene-induced mammary tumor. Tumors arising in an (I/LnJ x PyMT) F1 hybrid background appeared earlier than in the FVB/N-TgN(MMTV-PyVT)(634Mul) parent, but accumulated less tumor mass, indicating a net decrease in tumor growth. Quantitative genetic mapping in a backcross identified three loci that were associated with the decreased proliferative capacity of the I/LnJ F1 tumors. Molecular analysis of the tumors suggests that these loci may act by restricting the tumor's ability to recruit microvessels. The three loci, designated Mmtg1-3, are unlinked to the angiogenic genes Fgf2, Flt1, Flk4, Flk1, Vegf, and Vegfc, as well as the precursors of the endogenous antiangiogenic molecules angiostatin and endostatin. The Mmtg loci may therefore provide novel targets for antiangiogenic therapeutic strategies.