Enhancer release and retargeting activates disease-susceptibility genes.

The functional engagement between an enhancer and its target promoter ensures precise gene transcription1. Understanding the basis of promoter choice by enhancers has important implications for health and disease. Here we report that functional loss of a preferred promoter can release its partner enhancer to loop to and activate an alternative promoter (or alternative promoters) in the neighbourhood. We refer to this target-switching process as 'enhancer release and retargeting'. Genetic deletion, motif perturbation or mutation, and dCas9-mediated CTCF tethering reveal that promoter choice by an enhancer can be determined by the binding of CTCF at promoters, in a cohesin-dependent manner-consistent with a model of 'enhancer scanning' inside the contact domain. Promoter-associated CTCF shows a lower affinity than that at chromatin domain boundaries and often lacks a preferred motif orientation or a partnering CTCF at the cognate enhancer, suggesting properties distinct from boundary CTCF. Analyses of cancer mutations, data from the GTEx project and risk loci from genome-wide association studies, together with a focused CRISPR interference screen, reveal that enhancer release and retargeting represents an overlooked mechanism that underlies the activation of disease-susceptibility genes, as exemplified by a risk locus for Parkinson's disease (NUCKS1-RAB7L1) and three loci associated with cancer (CLPTM1L-TERT, ZCCHC7-PAX5 and PVT1-MYC).

TREX1 is required for microglial cholesterol homeostasis and oligodendrocyte terminal differentiation in human neural assembloids.

Three Prime Repair Exonuclease 1 (TREX1) gene mutations have been associated with Aicardi-Goutières Syndrome (AGS) - a rare, severe pediatric autoimmune disorder that primarily affects the brain and has a poorly understood etiology. Microglia are brain-resident macrophages indispensable for brain development and implicated in multiple neuroinflammatory diseases. However, the role of TREX1 - a DNase that cleaves cytosolic nucleic acids, preventing viral- and autoimmune-related inflammatory responses - in microglia biology remains to be elucidated. Here, we leverage a model of human embryonic stem cell (hESC)-derived engineered microglia-like cells, bulk, and single-cell transcriptomics, optical and transmission electron microscopy, and three-month-old assembloids composed of microglia and oligodendrocyte-containing organoids to interrogate TREX1 functions in human microglia. Our analyses suggest that TREX1 influences cholesterol metabolism, leading to an active microglial morphology with increased phagocytosis in the absence of TREX1. Notably, regulating cholesterol metabolism with an HMG-CoA reductase inhibitor, FDA-approved atorvastatin, rescues these microglial phenotypes. Functionally, TREX1 in microglia is necessary for the transition from gliogenic intermediate progenitors known as pre-oligodendrocyte precursor cells (pre-OPCs) to precursors of the oligodendrocyte lineage known as OPCs, impairing oligodendrogenesis in favor of astrogliogenesis in human assembloids. Together, these results suggest routes for therapeutic intervention in pathologies such as AGS based on microglia-specific molecular and cellular mechanisms.

Nuclear receptor-enhanced transcription requires motor- and LSD1-dependent gene networking in interchromatin granules.

While the transcriptional machinery has been extensively dissected at the molecular level, little is known about regulation of chromosomal organization in the three-dimensional space of the nucleus to achieve integrated transcriptional responses to diverse signaling events. Here, we report that ligand induces rapid interchromosomal interactions among subsets of estrogen receptor alpha-bound transcription units, with a dramatic reorganization of nuclear territories requiring nuclear actin/myosin-I transport machinery, dynein light chain 1 (DLC1), and a specific subset of transcriptional coactivators and chromatin remodeling complexes. We establish a requirement for the histone lysine demethylase, LSD1, in directing specific interchromosomal interaction loci to distinct interchromatin granules, long thought to be "storage" sites for splicing machinery, and demonstrate that these three-dimensional motor-dependent interactions are required to achieve enhanced transcription of specific estrogen-receptor target genes. These findings reveal roles for the modulation of nuclear architecture in orchestrating regulated gene-expression programs in the mammalian nucleus.

Transcriptional regulation of the sodium channel gene (SCN5A) by GATA4 in human heart.

Aberrant expression of the sodium channel gene (SCN5A) has been proposed to disrupt cardiac action potential and cause human cardiac arrhythmias, but the mechanisms of SCN5A gene regulation and dysregulation still remain largely unexplored. To gain insight into the transcriptional regulatory networks of SCN5A, we surveyed the promoter and first intronic regions of the SCN5A gene, predicting the presence of several binding sites for GATA transcription factors (TFs). Consistent with this prediction, chromatin immunoprecipitation (ChIP) and sequential ChIP (Re-ChIP) assays show co-occupancy of cardiac GATA TFs GATA4 and GATA5 on promoter and intron 1 SCN5A regions in fresh-frozen human left ventricle samples. Gene reporter experiments show GATA4 and GATA5 synergism in the activation of the SCN5A promoter, and its dependence on predicted GATA binding sites. GATA4 and GATA6 mRNAs are robustly expressed in fresh-frozen human left ventricle samples as measured by highly sensitive droplet digital PCR (ddPCR). GATA5 mRNA is marginally but still clearly detected in the same samples. Importantly, GATA4 mRNA levels are strongly and positively correlated with SCN5A transcript levels in the human heart. Together, our findings uncover a novel mechanism of GATA TFs in the regulation of the SCN5A gene in human heart tissue. Our studies suggest that GATA5 but especially GATA4 are main contributors to SCN5A gene expression, thus providing a new paradigm of SCN5A expression regulation that may shed new light into the understanding of cardiac disease.

Reprogramming transcription by distinct classes of enhancers functionally defined by eRNA.

Mammalian genomes are populated with thousands of transcriptional enhancers that orchestrate cell-type-specific gene expression programs, but how those enhancers are exploited to institute alternative, signal-dependent transcriptional responses remains poorly understood. Here we present evidence that cell-lineage-specific factors, such as FoxA1, can simultaneously facilitate and restrict key regulated transcription factors, exemplified by the androgen receptor (AR), to act on structurally and functionally distinct classes of enhancer. Consequently, FoxA1 downregulation, an unfavourable prognostic sign in certain advanced prostate tumours, triggers dramatic reprogramming of the hormonal response by causing a massive switch in AR binding to a distinct cohort of pre-established enhancers. These enhancers are functional, as evidenced by the production of enhancer-templated non-coding RNA (eRNA) based on global nuclear run-on sequencing (GRO-seq) analysis, with a unique class apparently requiring no nucleosome remodelling to induce specific enhancer-promoter looping and gene activation. GRO-seq data also suggest that liganded AR induces both transcription initiation and elongation. Together, these findings reveal a large repository of active enhancers that can be dynamically tuned to elicit alternative gene expression programs, which may underlie many sequential gene expression events in development, cell differentiation and disease progression.

PHF8 mediates histone H4 lysine 20 demethylation events involved in cell cycle progression.

While reversible histone modifications are linked to an ever-expanding range of biological functions, the demethylases for histone H4 lysine 20 and their potential regulatory roles remain unknown. Here we report that the PHD and Jumonji C (JmjC) domain-containing protein, PHF8, while using multiple substrates, including H3K9me1/2 and H3K27me2, also functions as an H4K20me1 demethylase. PHF8 is recruited to promoters by its PHD domain based on interaction with H3K4me2/3 and controls G1-S transition in conjunction with E2F1, HCF-1 (also known as HCFC1) and SET1A (also known as SETD1A), at least in part, by removing the repressive H4K20me1 mark from a subset of E2F1-regulated gene promoters. Phosphorylation-dependent PHF8 dismissal from chromatin in prophase is apparently required for the accumulation of H4K20me1 during early mitosis, which might represent a component of the condensin II loading process. Accordingly, the HEAT repeat clusters in two non-structural maintenance of chromosomes (SMC) condensin II subunits, N-CAPD3 and N-CAPG2 (also known as NCAPD3 and NCAPG2, respectively), are capable of recognizing H4K20me1, and ChIP-Seq analysis demonstrates a significant overlap of condensin II and H4K20me1 sites in mitotic HeLa cells. Thus, the identification and characterization of an H4K20me1 demethylase, PHF8, has revealed an intimate link between this enzyme and two distinct events in cell cycle progression.

Decoding a signature-based model of transcription cofactor recruitment dictated by cardinal cis-regulatory elements in proximal promoter regions.

Genome-wide maps of DNase I hypersensitive sites (DHSs) reveal that most human promoters contain perpetually active cis-regulatory elements between -150 bp and +50 bp (-150/+50 bp) relative to the transcription start site (TSS). Transcription factors (TFs) recruit cofactors (chromatin remodelers, histone/protein-modifying enzymes, and scaffold proteins) to these elements in order to organize the local chromatin structure and coordinate the balance of post-translational modifications nearby, contributing to the overall regulation of transcription. However, the rules of TF-mediated cofactor recruitment to the -150/+50 bp promoter regions remain poorly understood. Here, we provide evidence for a general model in which a series of cis-regulatory elements (here termed 'cardinal' motifs) prefer acting individually, rather than in fixed combinations, within the -150/+50 bp regions to recruit TFs that dictate cofactor signatures distinctive of specific promoter subsets. Subsequently, human promoters can be subclassified based on the presence of cardinal elements and their associated cofactor signatures. In this study, furthermore, we have focused on promoters containing the nuclear respiratory factor 1 (NRF1) motif as the cardinal cis-regulatory element and have identified the pervasive association of NRF1 with the cofactor lysine-specific demethylase 1 (LSD1/KDM1A). This signature might be distinctive of promoters regulating nuclear-encoded mitochondrial and other particular genes in at least some cells. Together, we propose that decoding a signature-based, expanded model of control at proximal promoter regions should lead to a better understanding of coordinated regulation of gene transcription.

Presenilins regulate neurotrypsin gene expression and neurotrypsin-dependent agrin cleavage via cyclic AMP response element-binding protein (CREB) modulation.

Presenilins, the catalytic components of the γ-secretase complex, are upstream regulators of multiple cellular pathways via regulation of gene transcription. However, the underlying mechanisms and the genes regulated by these pathways are poorly characterized. In this study, we identify Tequila and its mammalian ortholog Prss12 as genes negatively regulated by presenilins in Drosophila larval brains and mouse embryonic fibroblasts, respectively. Prss12 encodes the serine protease neurotrypsin, which cleaves the heparan sulfate proteoglycan agrin. Altered neurotrypsin activity causes serious synaptic and cognitive defects; despite this, the molecular processes regulating neurotrypsin expression and activity are poorly understood. Using γ-secretase drug inhibitors and presenilin mutants in mouse embryonic fibroblasts, we found that a mature γ-secretase complex was required to repress neurotrypsin expression and agrin cleavage. We also determined that PSEN1 endoproteolysis or processing of well known γ-secretase substrates was not essential for this process. At the transcriptional level, PSEN1/2 removal induced cyclic AMP response element-binding protein (CREB)/CREB-binding protein binding, accumulation of activating histone marks at the neurotrypsin promoter, and neurotrypsin transcriptional and functional up-regulation that was dependent on GSK3 activity. Upon PSEN1/2 reintroduction, this active epigenetic state was replaced by a methyl CpG-binding protein 2 (MeCP2)-containing repressive state and reduced neurotrypsin expression. Genome-wide analysis revealed hundreds of other mouse promoters in which CREB binding is similarly modulated by the presence/absence of presenilins. Our study thus identifies Tequila and neurotrypsin as new genes repressed by presenilins and reveals a novel mechanism used by presenilins to modulate CREB signaling based on controlling CREB recruitment.

Mechanisms establishing TLR4-responsive activation states of inflammatory response genes.

Precise control of the innate immune response is required for resistance to microbial infections and maintenance of normal tissue homeostasis. Because this response involves coordinate regulation of hundreds of genes, it provides a powerful biological system to elucidate the molecular strategies that underlie signal- and time-dependent transitions of gene expression. Comprehensive genome-wide analysis of the epigenetic and transcription status of the TLR4-induced transcriptional program in macrophages suggests that Toll-like receptor 4 (TLR4)-dependent activation of nearly all immediate/early- (I/E) and late-response genes results from a sequential process in which signal-independent factors initially establish basal levels of gene expression that are then amplified by signal-dependent transcription factors. Promoters of I/E genes are distinguished from those of late genes by encoding a distinct set of signal-dependent transcription factor elements, including TATA boxes, which lead to preferential binding of TBP and basal enrichment for RNA polymerase II immediately downstream of transcriptional start sites. Global nuclear run-on (GRO) sequencing and total RNA sequencing further indicates that TLR4 signaling markedly increases the overall rates of both transcriptional initiation and the efficiency of transcriptional elongation of nearly all I/E genes, while RNA splicing is largely unaffected. Collectively, these findings reveal broadly utilized mechanisms underlying temporally distinct patterns of TLR4-dependent gene activation required for homeostasis and effective immune responses.

Targeted phasing of 2-200 kilobase DNA fragments with a short-read sequencer and a single-tube linked-read library method.

In the human genome, heterozygous sites refer to genomic positions with a different allele or nucleotide variant on the maternal and paternal chromosomes. Resolving these allelic differences by chromosomal copy, also known as phasing, is achievable on a short-read sequencer when using a library preparation method that captures long-range genomic information. TELL-Seq is a library preparation that captures long-range genomic information with the aid of molecular identifiers (barcodes). The same barcode is used to tag the reads derived from the same long DNA fragment within a range of up to 200 kilobases (kb), generating linked-reads. This strategy can be used to phase an entire genome. Here, we introduce a TELL-Seq protocol developed for targeted applications, enabling the phasing of enriched loci of varying sizes, purity levels, and heterozygosity. To validate this protocol, we phased 2-200 kb loci enriched with different methods: CRISPR/Cas9-mediated excision coupled with pulse-field electrophoresis for the longest fragments, CRISPR/Cas9-mediated protection from exonuclease digestion for mid-size fragments, and long PCR for the shortest fragments. All selected loci have known clinical relevance: BRCA1, BRCA2, MLH1, MSH2, MSH6, APC, PMS2, SCN5A-SCN10A, and PKI3CA. Collectively, the analyses show that TELL-Seq can accurately phase 2-200 kb targets using a short-read sequencer.

Opposing LSD1 complexes function in developmental gene activation and repression programmes.

Precise control of transcriptional programmes underlying metazoan development is modulated by enzymatically active co-regulatory complexes, coupled with epigenetic strategies. One thing that remains unclear is how specific members of histone modification enzyme families, such as histone methyltransferases and demethylases, are used in vivo to simultaneously orchestrate distinct developmental gene activation and repression programmes. Here, we report that the histone lysine demethylase, LSD1--a component of the CoREST-CtBP co-repressor complex--is required for late cell-lineage determination and differentiation during pituitary organogenesis. LSD1 seems to act primarily on target gene activation programmes, as well as in gene repression programmes, on the basis of recruitment of distinct LSD1-containing co-activator or co-repressor complexes. LSD1-dependent gene repression programmes can be extended late in development with the induced expression of ZEB1, a Krüppel-like repressor that can act as a molecular beacon for recruitment of the LSD1-containing CoREST-CtBP co-repressor complex, causing repression of an additional cohort of genes, such as Gh, which previously required LSD1 for activation. These findings suggest that temporal patterns of expression of specific components of LSD1 complexes modulate gene regulatory programmes in many mammalian organs.

Targeted Phasing of 2-200 Kilobase DNA Fragments with a Short-Read Sequencer and a Single-Tube Linked-Read Library Method.

In the human genome, heterozygous sites are genomic positions with different alleles inherited from each parent. On average, there is a heterozygous site every 1-2 kilobases (kb). Resolving whether two alleles in neighboring heterozygous positions are physically linked-that is, phased-is possible with a short-read sequencer if the sequencing library captures long-range information. TELL-Seq is a library preparation method based on millions of barcoded micro-sized beads that enables instrument-free phasing of a whole human genome in a single PCR tube. TELL-Seq incorporates a unique molecular identifier (barcode) to the short reads generated from the same high-molecular-weight (HMW) DNA fragment (known as 'linked-reads'). However, genome-scale TELL-Seq is not cost-effective for applications focusing on a single locus or a few loci. Here, we present an optimized TELL-Seq protocol that enables the cost-effective phasing of enriched loci (targets) of varying sizes, purity levels, and heterozygosity. Targeted TELL-Seq maximizes linked-read efficiency and library yield while minimizing input requirements, fragment collisions on microbeads, and sequencing burden. To validate the targeted protocol, we phased seven 180-200 kb loci enriched by CRISPR/Cas9-mediated excision coupled with pulse-field electrophoresis, four 20 kb loci enriched by CRISPR/Cas9-mediated protection from exonuclease digestion, and six 2-13 kb loci amplified by PCR. The selected targets have clinical and research relevance (BRCA1, BRCA2, MLH1, MSH2, MSH6, APC, PMS2, SCN5A-SCN10A, and PKI3CA). These analyses reveal that targeted TELL-Seq provides a reliable way of phasing allelic variants within targets (2-200 kb in length) with the low cost and high accuracy of short-read sequencing.

Analysis of Clonal Composition in Human iPSC and ESC and Derived 2D and 3D Differentiated Cultures.

Human induced pluripotent and embryonic stem cell cultures (hiPSC/hESC) are phenotypically heterogeneous and prone to clonal deviations during subculturing and differentiation. Clonal deviations often emerge unnoticed, but they can change the biology of the cell culture with a negative impact on experimental reproducibility. Here, we describe a computational workflow to profile the bulk clonal composition in a hiPSC/hESC culture that can also be used to infer clonal deviations. This workflow processes data obtained with two versions of the same method. The two versions-epigenetic and transcriptomic-rely on a mechanism of stochastic H3K4me3 deposition during hiPSC/hESC derivation. This mechanism generates a signature of ten or more H3K4me3-enriched clustered protocadherin (PCDH) promoters distinct in every single cell. The aggregate of single-cell signatures provides an identificatory feature in every hiPSC/hESC line. This feature is stably transmitted to the cell progeny of the culture even after differentiation unless there is a clonal deviation event that changes the internal balance of single-cell signatures. H3K4me3 signatures can be profiled by chromatin immunoprecipitation and next-generation sequencing (ChIP-seq). Alternatively, an equivalent PCDH-expression version can be profiled by RNA-seq in PCDH-expressing hiPSC/hESC-derived cells (such as neurons, astrocytes, and cardiomyocytes; and, in long-term cultures, such as cerebral organoids). Notably, our workflow can also distinguish genetically identical hiPSC/hESC lines derived from the same patient or generated in the same editing process. Together, we propose a method to improve data sharing and reproducibility in the hiPSC and hESC fields.

Enhancing nuclear receptor-induced transcription requires nuclear motor and LSD1-dependent gene networking in interchromatin granules.

Although the role of liganded nuclear receptors in mediating coactivator/corepressor exchange is well-established, little is known about the potential regulation of chromosomal organization in the 3-dimensional space of the nucleus in achieving integrated transcriptional responses to diverse signaling events. Here, we report that ligand induces rapid interchromosomal interactions among specific subsets of estrogen receptor alpha-bound transcription units, with a dramatic reorganization of nuclear territories, which depends on the actions of nuclear actin/myosin-I machinery and dynein light chain 1. The histone lysine demethylase, LSD1, is required for these ligand-induced interactive loci to associate with distinct interchromatin granules, long thought to serve as "storage" sites for the splicing machinery, some critical transcription elongation factors, and various chromatin remodeling complexes. We demonstrate that this 2-step nuclear rearrangement is essential for achieving enhanced, coordinated transcription of nuclear receptor target genes.

Analysis of Brugada syndrome loci reveals that fine-mapping clustered GWAS hits enhances the annotation of disease-relevant variants.

Genome-wide association studies (GWASs) are instrumental in identifying loci harboring common single-nucleotide variants (SNVs) that affect human traits and diseases. GWAS hits emerge in clusters, but the focus is often on the most significant hit in each trait- or disease-associated locus. The remaining hits represent SNVs in linkage disequilibrium (LD) and are considered redundant and thus frequently marginally reported or exploited. Here, we interrogate the value of integrating the full set of GWAS hits in a locus repeatedly associated with cardiac conduction traits and arrhythmia, SCN5A-SCN10A. Our analysis reveals 5 common 7-SNV haplotypes (Hap1-5) with 2 combinations associated with life-threatening arrhythmia-Brugada syndrome (the risk Hap1/1 and protective Hap2/3 genotypes). Hap1 and Hap2 share 3 SNVs; thus, this analysis suggests that assuming redundancy among clustered GWAS hits can lead to confounding disease-risk associations and supports the need to deconstruct GWAS data in the context of haplotype composition.

Super-Enhancer Redistribution as a Mechanism of Broad Gene Dysregulation in Repeatedly Drug-Treated Cancer Cells.

Cisplatin is an antineoplastic drug administered at suboptimal and intermittent doses to avoid life-threatening effects. Although this regimen shortly improves symptoms in the short term, it also leads to more malignant disease in the long term. We describe a multilayered analysis ranging from chromatin to translation-integrating chromatin immunoprecipitation sequencing (ChIP-seq), global run-on sequencing (GRO-seq), RNA sequencing (RNA-seq), and ribosome profiling-to understand how cisplatin confers (pre)malignant features by using a well-established ovarian cancer model of cisplatin exposure. This approach allows us to segregate the human transcriptome into gene modules representing distinct regulatory principles and to characterize that the most cisplatin-disrupted modules are associated with underlying events of super-enhancer plasticity. These events arise when cancer cells initiate without ultimately ending the program of drug-stimulated death. Using a PageRank-based algorithm, we predict super-enhancer regulator ISL1 as a driver of this plasticity and validate this prediction by using CRISPR/dCas9-KRAB inhibition (CRISPRi) and CRISPR/dCas9-VP64 activation (CRISPRa) tools. Together, we propose that cisplatin reprograms cancer cells when inducing them to undergo near-to-death experiences.

Tissue-specific transcription reprogramming promotes liver metastasis of colorectal cancer.

Metastasis, the development of secondary malignant growths at a distance from a primary tumor, is the cause of death for 90% of cancer patients, but little is known about how metastatic cancer cells adapt to and colonize new tissue environments. Here, using clinical samples, patient-derived xenograft (PDX) samples, PDX cells, and primary/metastatic cell lines, we discovered that liver metastatic colorectal cancer (CRC) cells lose their colon-specific gene transcription program yet gain a liver-specific gene transcription program. We showed that this transcription reprogramming is driven by a reshaped epigenetic landscape of both typical enhancers and super-enhancers. Further, we identified that the liver-specific transcription factors FOXA2 and HNF1A can bind to the gained enhancers and activate the liver-specific gene transcription, thereby driving CRC liver metastasis. Importantly, similar transcription reprogramming can be observed in multiple cancer types. Our data suggest that reprogrammed tissue-specific transcription promotes metastasis and should be targeted therapeutically.

Chromatin establishes an immature version of neuronal protocadherin selection during the naive-to-primed conversion of pluripotent stem cells.

In the mammalian genome, the clustered protocadherin (cPCDH) locus provides a paradigm for stochastic gene expression with the potential to generate a unique cPCDH combination in every neuron. Here we report a chromatin-based mechanism that emerges during the transition from the naive to the primed states of cell pluripotency and reduces, by orders of magnitude, the combinatorial potential in the human cPCDH locus. This mechanism selectively increases the frequency of stochastic selection of a small subset of cPCDH genes after neuronal differentiation in monolayers, 10-month-old cortical organoids and engrafted cells in the spinal cords of rats. Signs of these frequent selections can be observed in the brain throughout fetal development and disappear after birth, except in conditions of delayed maturation such as Down's syndrome. We therefore propose that a pattern of limited cPCDH-gene expression diversity is maintained while human neurons still retain fetal-like levels of maturation.

Mnt-deficient mammary glands exhibit impaired involution and tumors with characteristics of myc overexpression.

The proto-oncogene c-Myc plays a central role in cell growth and the development of human tumors. c-Myc interacts with Max and Myc-Max complexes bind to E-box and related sequences to activate transcription. Max also interacts with Mnt but Mnt-Max complexes repress transcription when bound to these sequences. MNT maps to human chromosome 17p13.3, a region frequently deleted in various human tumors, including mammary gland tumors. Consistent with the possibility that Mnt functions as a Myc antagonist, Mnt-deficient fibroblasts exhibit many of the hallmark characteristics of cells that overexpress Myc, and conditional (Cre/Lox) inactivation of Mnt in mammary gland epithelium leads to adenocarcinomas. Here, we further characterize mammary gland tissue following conditional deletion of Mnt in the mammary gland. We show that loss of Mnt severely disrupts mammary gland involution and leads to hyperplastic ducts associated with reduced numbers of apoptotic cells. These findings suggest that loss of Mnt in mammary tissue has similarities to Myc overexpression. We tested this directly by using promoter array analysis and mRNA expression analysis by oligonucleotide arrays. We found that Mnt and c-Myc bound to similar promoters in tumors from MMTV-c-Myc transgenic mice, and mRNA expression patterns were similar between mammary tumors from MMTV-Cre/Mnt(KO/CKO) and MMTV-c-Myc transgenic mice. These results reveal an important role for Mnt in pregnancy-associated mammary gland development and suggest that mammary gland tumorigenesis in the absence of Mnt is analogous to that caused by Myc deregulation.