Rat Mammary Epithelial Cell Transplantation into the Interscapular White Fat Pad.

As early as the 1970s, researchers have successfully transplanted mammary epithelial cells into the interscapular white fat pad of rats. Grafting mammary epithelium using transplantation techniques takes advantage of the hormonal environment provided by the adolescent rodent host. These studies are ideally suited to explore the impact of various biological manipulations on mammary gland development and dissect many aspects of mammary gland biology. A common, but limiting, feature is that transplanted epithelial cells are strongly influenced by the surrounding stroma and outcompeted by endogenous epithelium; to utilize native mammary tissue, the abdominal-inguinal white fat pad must be cleared to remove host mammary epithelium prior to the transplantation. A major obstacle when using the rat model organism is that clearing the developing mammary tree in post-weaned rats is not efficient. When transplanted into gland-free fat pads, donor epithelial cells can repopulate the cleared host fat pad and form a functional mammary gland. The interscapular fat pad is an alternative location for these grafts. A major advantage is that it lacks ductal structures yet provides the normal stroma that is necessary to promote epithelial outgrowth and is easily accessible in the rat. Another major advantage of this technique is that it is minimally invasive, because it eliminates the need to cauterize and remove the growing endogenous mammary tree. Additionally, the interscapular fat pad contains a medial blood vessel that can be used to separate sites for grafting. Because the endogenous glands remain intact, this technique can also be used for studies comparing the endogenous mammary gland to the transplanted gland. This paper describes the method of mammary epithelial cell transplantation into the interscapular white fat pad of rats.

Deletion of Cdkn1b in ACI rats leads to increased proliferation and pregnancy-associated changes in the mammary gland due to perturbed systemic endocrine environment.

Mammary epithelial progenitors are the normal cell-of-origin of breast cancer. We previously defined a population of p27+ quiescent hormone-responsive progenitor cells in the normal human breast whose frequency associates with breast cancer risk. Here, we describe that deletion of the Cdkn1b gene encoding the p27 cyclin-dependent kinase inhibitor in the estrogen-induced mammary tumor-susceptible ACI rat strain leads to a decrease in the relative frequencies of Cd49b+ mammary luminal epithelial progenitors and pregnancy-related differentiation. We show by comprehensive gene expression profiling of purified progenitor and differentiated mammary epithelial cell populations that p27 deletion has the most pronounced effects on luminal progenitors. Cdkn1b-/- females have decreased fertility, but rats that are able to get pregnant had normal litter size and were able to nurse their pups implying that loss of p27 in ACI rats does not completely abrogate ovarian function and lactation. Reciprocal mammary gland transplantation experiments indicate that the p27-loss-induced changes in mammary epithelial cells are not only caused by alterations in their intrinsic properties, but are likely due to altered hormonal signaling triggered by the perturbed systemic endocrine environment observed in Cdkn1b-/- females. We also observed a decrease in the frequency of mammary epithelial cells positive for progesterone receptor (Pr) and FoxA1, known direct transcriptional targets of the estrogen receptor (Erα), and an increase in phospho-Stat5 positive cells commonly induced by prolactin (Prl). Characterization of genome-wide Pr chromatin binding revealed distinct binding patterns in mammary epithelial cells of Cdkn1b+/+ and Cdkn1b-/- females and enrichment in genes with known roles in Notch, ErbB, leptin, and Erα signaling and regulation of G1-S transition. Our data support a role for p27 in regulating the pool size of hormone-responsive luminal progenitors that could impact breast cancer risk.

Deletion of the murine ortholog of the 8q24 gene desert has anti-cancer effects in transgenic mammary cancer models.

BACKGROUND: The gene desert on human chromosomal band 8q24 harbors multiple genetic variants associated with common cancers, including breast cancer. The locus, including the gene desert and its flanking genes, MYC, PVT1 and FAM84B, is also frequently amplified in human breast cancer. We generated a megadeletion (MD) mouse model lacking 430-Kb of sequence orthologous to the breast cancer-associated region in the gene desert. The goals were to examine the effect of the deletion on mammary cancer development and on transcript level regulation of the candidate genes within the locus. METHODS: The MD allele was engineered using the MICER system in embryonic stem cells and bred onto 3 well-characterized transgenic models for breast cancer, namely MMTV-PyVT, MMTV-neu and C3(1)-TAg. Mammary tumor growth, latency, multiplicity and metastasis were compared between homozygous MD and wild type mice carrying the transgenes. A reciprocal mammary gland transplantation assay was conducted to distinguish mammary cell-autonomous from non-mammary cell-autonomous anti-cancer effects. Gene expression analysis was done using quantitative real-time PCR. Chromatin interactions were evaluated by 3C. Gene-specific patient outcome data were analysed using the METABRIC and TCGA data sets through the cBioPortal website. RESULTS: Mice homozygous for the MD allele are viable, fertile, lactate sufficiently to nourish their pups, but maintain a 10% lower body weight mainly due to decreased adiposity. The deletion interferes with mammary tumorigenesis in mouse models for luminal and basal breast cancer. In the MMTV-PyVT model the mammary cancer-reducing effects of the allele are mammary cell-autonomous. We found organ-specific effects on transcript level regulation, with Myc and Fam84b being downregulated in mammary gland, prostate and mammary tumor samples. Through analysis using the METABRIC and TCGA datasets, we provide evidence that MYC and FAM84B are frequently co-amplified in breast cancer, but in contrast with MYC, FAM84B is frequently overexpressed in the luminal subtype, whereas MYC activity affect basal breast cancer outcomes. CONCLUSION: Deletion of a breast cancer-associated non-protein coding region affects mammary cancer development in 3 transgenic mouse models. We propose Myc as a candidate susceptibility gene, regulated by the gene desert locus, and a potential role for Fam84b in modifying breast cancer development.

The Non-coding Mammary Carcinoma Susceptibility Locus, Mcs5c, Regulates Pappa Expression via Age-Specific Chromatin Folding and Allele-Dependent DNA Methylation.

In understanding the etiology of breast cancer, the contributions of both genetic and environmental risk factors are further complicated by the impact of breast developmental stage. Specifically, the time period ranging from childhood to young adulthood represents a critical developmental window in a woman's life when she is more susceptible to environmental hazards that may affect future breast cancer risk. Although the effects of environmental exposures during particular developmental Windows of Susceptibility (WOS) are well documented, the genetic mechanisms governing these interactions are largely unknown. Functional characterization of the Mammary Carcinoma Susceptibility 5c, Mcs5c, congenic rat model of breast cancer at various stages of mammary gland development was conducted to gain insight into the interplay between genetic risk factors and WOS. Using quantitative real-time PCR, chromosome conformation capture, and bisulfite pyrosequencing we have found that Mcs5c acts within the mammary gland to regulate expression of the neighboring gene Pappa during a critical mammary developmental time period in the rat, corresponding to the human young adult WOS. Pappa has been shown to positively regulate the IGF signaling pathway, which is required for proper mammary gland/breast development and is of increasing interest in breast cancer pathogenesis. Mcs5c-mediated regulation of Pappa appears to occur through age-dependent and mammary gland-specific chromatin looping, as well as genotype-dependent CpG island shore methylation. This represents, to our knowledge, the first insight into cellular mechanisms underlying the WOS phenomenon and demonstrates the influence developmental stage can have on risk locus functionality. Additionally, this work represents a novel model for further investigation into how environmental factors, together with genetic factors, modulate breast cancer risk in the context of breast developmental stage.

The gene desert mammary carcinoma susceptibility locus Mcs1a regulates Nr2f1 modifying mammary epithelial cell differentiation and proliferation.

Genome-wide association studies have revealed that many low-penetrance breast cancer susceptibility loci are located in non-protein coding genomic regions; however, few have been characterized. In a comparative genetics approach to model such loci in a rat breast cancer model, we previously identified the mammary carcinoma susceptibility locus Mcs1a. We now localize Mcs1a to a critical interval (277 Kb) within a gene desert. Mcs1a reduces mammary carcinoma multiplicity by 50% and acts in a mammary cell-autonomous manner. We developed a megadeletion mouse model, which lacks 535 Kb of sequence containing the Mcs1a ortholog. Global gene expression analysis by RNA-seq revealed that in the mouse mammary gland, the orphan nuclear receptor gene Nr2f1/Coup-tf1 is regulated by Mcs1a. In resistant Mcs1a congenic rats, as compared with susceptible congenic control rats, we found Nr2f1 transcript levels to be elevated in mammary gland, epithelial cells, and carcinoma samples. Chromatin looping over ∼820 Kb of sequence from the Nr2f1 promoter to a strongly conserved element within the Mcs1a critical interval was identified. This element contains a 14 bp indel polymorphism that affects a human-rat-mouse conserved COUP-TF binding motif and is a functional Mcs1a candidate. In both the rat and mouse models, higher Nr2f1 transcript levels are associated with higher abundance of luminal mammary epithelial cells. In both the mouse mammary gland and a human breast cancer global gene expression data set, we found Nr2f1 transcript levels to be strongly anti-correlated to a gene cluster enriched in cell cycle-related genes. We queried 12 large publicly available human breast cancer gene expression studies and found that the median NR2F1 transcript level is consistently lower in 'triple-negative' (ER-PR-HER2-) breast cancers as compared with 'receptor-positive' breast cancers. Our data suggest that the non-protein coding locus Mcs1a regulates Nr2f1, which is a candidate modifier of differentiation, proliferation, and mammary cancer risk.

EBSeq: an empirical Bayes hierarchical model for inference in RNA-seq experiments.

MOTIVATION: Messenger RNA expression is important in normal development and differentiation, as well as in manifestation of disease. RNA-seq experiments allow for the identification of differentially expressed (DE) genes and their corresponding isoforms on a genome-wide scale. However, statistical methods are required to ensure that accurate identifications are made. A number of methods exist for identifying DE genes, but far fewer are available for identifying DE isoforms. When isoform DE is of interest, investigators often apply gene-level (count-based) methods directly to estimates of isoform counts. Doing so is not recommended. In short, estimating isoform expression is relatively straightforward for some groups of isoforms, but more challenging for others. This results in estimation uncertainty that varies across isoform groups. Count-based methods were not designed to accommodate this varying uncertainty, and consequently, application of them for isoform inference results in reduced power for some classes of isoforms and increased false discoveries for others. RESULTS: Taking advantage of the merits of empirical Bayesian methods, we have developed EBSeq for identifying DE isoforms in an RNA-seq experiment comparing two or more biological conditions. Results demonstrate substantially improved power and performance of EBSeq for identifying DE isoforms. EBSeq also proves to be a robust approach for identifying DE genes. AVAILABILITY AND IMPLEMENTATION: An R package containing examples and sample datasets is available at http://www.biostat.wisc.edu/kendzior/EBSEQ/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Effective flow cytometric phenotyping of cells using minimal amounts of antibody.

Here we introduce a modified antibody staining method that uses up to 80% less antibody for flow cytometry. We demonstrate this method for the detection of antigens expressed at high, moderate, or low levels in mouse and rat lymphocytes as well as mouse mammary epithelial cells. We obtained reproducibly accurate results for the detection of up to seven parameters for activation induced-proliferation, cell cycle analysis, and phenotyping of cell-surface and intracellular antigens.

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.

Quantification of epithelial cell differentiation in mammary glands and carcinomas from DMBA- and MNU-exposed rats.

Rat mammary carcinogenesis models have been used extensively to study breast cancer initiation, progression, prevention, and intervention. Nevertheless, quantitative molecular data on epithelial cell differentiation in mammary glands of untreated and carcinogen-exposed rats is limited. Here, we describe the characterization of rat mammary epithelial cells (RMECs) by multicolor flow cytometry using antibodies against cell surface proteins CD24, CD29, CD31, CD45, CD49f, CD61, Peanut Lectin, and Thy-1, intracellular proteins CK14, CK19, and FAK, along with phalloidin and Hoechst staining. We identified the luminal and basal/myoepithelial populations and actively dividing RMECs. In inbred rats susceptible to mammary carcinoma development, we quantified the changes in differentiation of the RMEC populations at 1, 2, and 4 weeks after exposure to mammary carcinogens DMBA and MNU. DMBA exposure did not alter the percentage of basal or luminal cells, but upregulated CD49f (Integrin α6) expression and increased cell cycle activity. MNU exposure resulted in a temporary disruption of the luminal/basal ratio and no CD49f upregulation. When comparing DMBA- or MNU-induced mammary carcinomas, the RMEC differentiation profiles are indistinguishable. The carcinomas compared with mammary glands from untreated rats, showed upregulation of CD29 (Integrin ß1) and CD49f expression, increased FAK (focal adhesion kinase) activation especially in the CD29hi population, and decreased CD61 (Integrin ß3) expression. This study provides quantitative insight into the protein expression phenotypes underlying RMEC differentiation. The results highlight distinct RMEC differentiation etiologies of DMBA and MNU exposure, while the resulting carcinomas have similar RMEC differentiation profiles. The methodology and data will enhance rat mammary carcinogenesis models in the study of the role of epithelial cell differentiation in breast cancer.

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.

INTRODUCTION: Mechanisms 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. METHODS: We 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. RESULTS: We 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. CONCLUSIONS: These 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.

A mutation in Myo15 leads to Usher-like symptoms in LEW/Ztm-ci2 rats.

The LEW/Ztm-ci2 rat is an animal model for syndromal deafness that arose from a spontaneous mutation. Homozygous animals show locomotor abnormalities like lateralized circling behavior. Additionally, an impaired vision can be observed in some animals through behavioral studies. Syndromal deafness as well as retinal degeneration are features of the Usher syndrome in humans. In the present study, the mutation was identified as a base substitution (T->C) in exon 56 of Myo15, leading to an amino acid exchange from leucine (Leu) to proline (Pro) within the carboxy-terminal MyTH4 domain in the proteins' tail region. Myo15 mRNA was expressed in the retina as demonstrated for the first time with the help of in-situ hybridization and PCR. To characterize the visual phenotype, rats were examined by scotopic and photopic electroretinography and, additionally, histological analyses of the retinas were conducted. The complete loss of sight was detected along with a severe degeneration of photoreceptor cells. Interestingly, the manifestation of the disease does not solely depend on the mutation, but also on environmental factors. Since the LEW/Ztm-ci2 rat features the entire range of symptoms of the human Usher syndrome we think that this strain is an appropriate model for this disease. Our findings display that mutations in binding domains of myosin XV do not only cause non-syndromic hearing loss but can also lead to syndromic disorders including retinal dysfunction.

ENU mutagenesis to generate genetically modified rat models.

The rat is one of the most preferred model organisms in biomedical research and has been extremely useful for linking physiology and pathology to the genome. However, approaches to genetically modify specific genes in the rat germ line remain relatively scarce. To date, the most efficient approach for generating genetically modified rats has been the target-selected N-ethyl-N-nitrosourea (ENU) mutagenesis-based technology. Here, we describe the detailed protocols for ENU mutagenesis and mutant retrieval in the rat model organism.

Lack of DNA mismatch repair protein MSH6 in the rat results in hereditary non-polyposis colorectal cancer-like tumorigenesis.

To understand genetic instability in relation to tumorigenesis, experimental animal models have proven very useful. The DNA mismatch repair (MMR) machinery safeguards genomic integrity by repairing mismatches, insertion or deletion loops and responding to genotoxic agents. Here, we describe the functional characterization of a novel rat mutant model in which the MMR gene Msh6 has been genetically inactivated by N-ethyl-N-nitrosourea-driven target-selected mutagenesis. This model shows a robust mutator phenotype that is reflected by microsatellite instability and an increased germ line point mutation frequency. Consequently, these rats develop a spectrum of tumors with a high similarity to atypical hereditary non-polyposis colorectal cancer in humans. The MSH6 knockout rat complements existing models for studying genetic instable tumorigenesis as it provides experimental opportunities that are not available or suboptimal in current models.

Genetically engineered rat models for breast cancer.

Rat mammary cancers generally resemble many features of human breast cancer. With the recent developments in rat genetic engineering, the rat has become an excellent model system to study aspects of the molecular etiology of breast cancer. In this review, we describe the efforts to generate genetically engineered rat models for breast cancer.

Haplotype block structure is conserved across mammals.

Genetic variation in genomes is organized in haplotype blocks, and species-specific block structure is defined by differential contribution of population history effects in combination with mutation and recombination events. Haplotype maps characterize the common patterns of linkage disequilibrium in populations and have important applications in the design and interpretation of genetic experiments. Although evolutionary processes are known to drive the selection of individual polymorphisms, their effect on haplotype block structure dynamics has not been shown. Here, we present a high-resolution haplotype map for a 5-megabase genomic region in the rat and compare it with the orthologous human and mouse segments. Although the size and fine structure of haplotype blocks are species dependent, there is a significant interspecies overlap in structure and a tendency for blocks to encompass complete genes. Extending these findings to the complete human genome using haplotype map phase I data reveals that linkage disequilibrium values are significantly higher for equally spaced positions in genic regions, including promoters, as compared to intergenic regions, indicating that a selective mechanism exists to maintain combinations of alleles within potentially interacting coding and regulatory regions. Although this characteristic may complicate the identification of causal polymorphisms underlying phenotypic traits, conservation of haplotype structure may be employed for the identification and characterization of functionally important genomic regions.

Rat genetics: the next episode.

More than a thousand quantitative trait loci (QTLs) relevant to many aspects of complex human disease have been identified in a wide range of rat inbred lines over the past few decades. With the complete rat genome available, it is now time for the next stage in rat genetic research: the identification and causal confirmation of underlying polymorphisms and genes. Recent developments in single nucleotide polymorphism-based genetic markers and technologies to manipulate the rat genome will undoubtedly be important tools in this next episode.

Rats go genomic.

A report on the meeting 'Rat Genomics and Models', Cold Spring Harbor, USA, 8-11 December 2005.

Generation of gene knockouts and mutant models in the laboratory rat by ENU-driven target-selected mutagenesis.

OBJECTIVE: The rat is one of the most important model organisms for biomedical and pharmacological research. However, the generation of novel models for studying specific aspects of human diseases largely depends on selection for specific traits using existing rat strains, thereby solely depending on naturally occurring variation. This study aims to provide the tools to manipulate the rat genome in a more directed way. METHODS: We developed robust, automated, and scaleable reverse genetic methodology based on ENU (N-ethyl-N-nitrosourea)-driven target-selected mutagenesis. Optimal mutagenesis conditions have been determined in three different rat strains and a universal, rapid, and cost-effective dideoxy resequencing-based screening setup was established for mutation discovery. The effectiveness of the approach is illustrated by the identification of 120 induced mutations in a set of genes of interest, including six that result in unique rat knockout models due to the introduction of premature stop codons. In addition, 56 mutations were found that change amino acids, including critical residues in transmembrane domains of receptors and channels. CONCLUSIONS: The approach described here allows for the systematic generation of knockout and protein function altering alleles in the rat. The resulting induced rat models will be powerful tools for studying many aspects of a wide variety of human diseases.

Efficient single nucleotide polymorphism discovery in laboratory rat strains using wild rat-derived SNP candidates.

BACKGROUND: The laboratory rat (Rattus norvegicus) is an important model for studying many aspects of human health and disease. Detailed knowledge on genetic variation between strains is important from a biomedical, particularly pharmacogenetic point of view and useful for marker selection for genetic cloning and association studies. RESULTS: We show that Single Nucleotide Polymorphisms (SNPs) in commonly used rat strains are surprisingly well represented in wild rat isolates. Shotgun sequencing of 814 Kbp in one wild rat resulted in the identification of 485 SNPs as compared with the Brown Norway genome sequence. Genotyping 36 commonly used inbred rat strains showed that 84% of these alleles are also polymorphic in a representative set of laboratory rat strains. CONCLUSION: We postulate that shotgun sequencing in a wild rat sample and subsequent genotyping in multiple laboratory or domesticated strains rather than direct shotgun sequencing of multiple strains, could be the most efficient SNP discovery approach. For the rat, laboratory strains still harbor a large portion of the haplotypes present in wild isolates, suggesting a relatively recent common origin and supporting the idea that rat inbred strains, in contrast to mouse inbred strains, originate from a single species, R. norvegicus.