In Vivo Evidence for Serine Biosynthesis-Defined Sensitivity of Lung Metastasis, but Not of Primary Breast Tumors, to mTORC1 Inhibition.

In tumors, nutrient availability and metabolism are known to be important modulators of growth signaling. However, it remains elusive whether cancer cells that are growing out in the metastatic niche rely on the same nutrients and metabolic pathways to activate growth signaling as cancer cells within the primary tumor. We discovered that breast-cancer-derived lung metastases, but not the corresponding primary breast tumors, use the serine biosynthesis pathway to support mTORC1 growth signaling. Mechanistically, pyruvate uptake through Mct2 supported mTORC1 signaling by fueling serine biosynthesis-derived α-ketoglutarate production in breast-cancer-derived lung metastases. Consequently, expression of the serine biosynthesis enzyme PHGDH was required for sensitivity to the mTORC1 inhibitor rapamycin in breast-cancer-derived lung tumors, but not in primary breast tumors. In summary, we provide in vivo evidence that the metabolic and nutrient requirements to activate growth signaling differ between the lung metastatic niche and the primary breast cancer site.

Resf1 is a compound G4 quadruplex-associated tumor suppressor for triple negative breast cancer.

Patients with ER-negative breast cancer have the worst prognosis of all breast cancer subtypes, often experiencing rapid recurrence or progression to metastatic disease shortly after diagnosis. Given that metastasis is the primary cause of mortality in most solid tumors, understanding metastatic biology is crucial for effective intervention. Using a mouse systems genetics approach, we previously identified 12 genes associated with metastatic susceptibility. Here, we extend those studies to identify Resf1, a poorly characterized gene, as a novel metastasis susceptibility gene in ER- breast cancer. Resf1 is a large, unstructured protein with an evolutionarily conserved intron-exon structure, but with poor amino acid conservation. CRISPR or gene trap mouse models crossed to the Polyoma Middle-T antigen genetically engineered mouse model (MMTV-PyMT) demonstrated that reduction of Resf1 resulted in a significant increase in tumor growth, a shortened overall survival time, and increased incidence and number of lung metastases, consistent with patient data. Furthermore, an analysis of matched tail and primary tissues revealed loss of the wildtype copy in tumor tissue, consistent with Resf1 being a tumor suppressor. Mechanistic analysis revealed a potential role of Resf1 in transcriptional control through association with compound G4 quadruplexes in expressed sequences, particularly those associated with ribosomal biogenesis. These results suggest that loss of Resf1 enhances tumor progression in ER- breast cancer through multiple alterations in both transcriptional and translational control.

Host CLIC4 expression in the tumor microenvironment is essential for breast cancer metastatic competence.

The TGF-ß-regulated Chloride Intracellular Channel 4 (CLIC4) is an essential participant in the formation of breast cancer stroma. Here, we used data available from the TCGA and METABRIC datasets to show that CLIC4 expression was higher in breast cancers from younger women and those with early-stage metastatic disease. Elevated CLIC4 predicted poor outcome in breast cancer patients and was linked to the TGF-ß pathway. However, these associations did not reveal the underlying biological contribution of CLIC4 to breast cancer progression. Constitutive ablation of host Clic4 in two murine metastatic breast cancer models nearly eliminated lung metastases without reducing primary tumor weight, while tumor cells ablated of Clic4 retained metastatic capability in wildtype hosts. Thus, CLIC4 was required for host metastatic competence. Pre- and post-metastatic proteomic analysis identified circulating pro-metastatic soluble factors that differed in tumor-bearing CLIC4-deficient and wildtype hosts. Vascular abnormalities and necrosis increased in primary tumors from CLIC4-deficient hosts. Transcriptional profiles of both primary tumors and pre-metastatic lungs of tumor-bearing CLIC4-deficient hosts were consistent with a microenvironment where inflammatory pathways were elevated. Altogether, CLIC4 expression in human breast cancers may serve as a prognostic biomarker; therapeutic targeting of CLIC4 could reduce primary tumor viability and host metastatic competence.

The genomic landscape of metastasis in treatment-naïve breast cancer models.

Metastasis remains the principle cause of mortality for breast cancer and presents a critical challenge because secondary lesions are often refractory to conventional treatments. While specific genetic alterations are tightly linked to primary tumor development and progression, the role of genetic alteration in the metastatic process is not well-understood. The theory of tumor evolution postulated by Peter Nowell in 1976 has yet to be proven in the context of metastasis. Therefore, in order to investigate how somatic evolution contributes to breast cancer metastasis, we performed exome, whole genome, and RNA sequencing of matched metastatic and primary tumors from pre-clinical mouse models of breast cancer. Here we show that in a treatment-naïve setting, recurrent single nucleotide variants and copy number variation, but not gene fusion events, play key metastasis-driving roles in breast cancer. For instance, we identified recurrent mutations in Kras, a known driver of colorectal and lung tumorigenesis that has not been previously implicated in breast cancer metastasis. However, in a set of in vivo proof-of-concept experiments we show that the Kras G12D mutation is sufficient to significantly promote metastasis using three syngeneic allograft models. The work herein confirms the existence of metastasis-driving mutations and presents a novel framework to identify actionable metastasis-targeted therapies.

Aicardi-Goutières syndrome gene Rnaseh2c is a metastasis susceptibility gene in breast cancer.

Breast cancer is the second leading cause of cancer-related deaths in the United States, with the majority of these deaths due to metastatic lesions rather than the primary tumor. Thus, a better understanding of the etiology of metastatic disease is crucial for improving survival. Using a haplotype mapping strategy in mouse and shRNA-mediated gene knockdown, we identified Rnaseh2c, a scaffolding protein of the heterotrimeric RNase H2 endoribonuclease complex, as a novel metastasis susceptibility factor. We found that the role of Rnaseh2c in metastatic disease is independent of RNase H2 enzymatic activity, and immunophenotyping and RNA-sequencing analysis revealed engagement of the T cell-mediated adaptive immune response. Furthermore, the cGAS-Sting pathway was not activated in the metastatic cancer cells used in this study, suggesting that the mechanism of immune response in breast cancer is different from the mechanism proposed for Aicardi-Goutières Syndrome, a rare interferonopathy caused by RNase H2 mutation. These results suggest an important novel, non-enzymatic role for RNASEH2C during breast cancer progression and add Rnaseh2c to a panel of genes we have identified that together could determine patients with high risk for metastasis. These results also highlight a potential new target for combination with immunotherapies and may contribute to a better understanding of the etiology of Aicardi-Goutières Syndrome autoimmunity.

The Circadian Rhythm Gene Arntl2 Is a Metastasis Susceptibility Gene for Estrogen Receptor-Negative Breast Cancer.

Breast cancer mortality is primarily due to metastasis rather than primary tumors, yet relatively little is understood regarding the etiology of metastatic breast cancer. Previously, using a mouse genetics approach, we demonstrated that inherited germline polymorphisms contribute to metastatic disease, and that these single nucleotide polymorphisms (SNPs) could be used to predict outcome in breast cancer patients. In this study, a backcross between a highly metastatic (FVB/NJ) and low metastatic (MOLF/EiJ) mouse strain identified Arntl2, a gene encoding a circadian rhythm transcription factor, as a metastasis susceptibility gene associated with progression, specifically in estrogen receptor-negative breast cancer patients. Integrated whole genome sequence analysis with DNase hypersensitivity sites reveals SNPs in the predicted promoter of Arntl2. Using CRISPR/Cas9-mediated substitution of the MOLF promoter, we demonstrate that the SNPs regulate Arntl2 transcription and affect metastatic burden. Finally, analysis of SNPs associated with ARNTL2 expression in human breast cancer patients revealed reproducible associations of ARNTL2 expression quantitative trait loci (eQTL) SNPs with disease-free survival, consistent with the mouse studies.

Post-transcriptional Control of Tumor Cell Autonomous Metastatic Potential by CCR4-NOT Deadenylase CNOT7.

Accumulating evidence supports the role of an aberrant transcriptome as a driver of metastatic potential. Deadenylation is a general regulatory node for post-transcriptional control by microRNAs and other determinants of RNA stability. Previously, we demonstrated that the CCR4-NOT scaffold component Cnot2 is an inherited metastasis susceptibility gene. In this study, using orthotopic metastasis assays and genetically engineered mouse models, we show that one of the enzymatic subunits of the CCR4-NOT complex, Cnot7, is also a metastasis modifying gene. We demonstrate that higher expression of Cnot7 drives tumor cell autonomous metastatic potential, which requires its deadenylase activity. Furthermore, metastasis promotion by CNOT7 is dependent on interaction with CNOT1 and TOB1. CNOT7 ribonucleoprotein-immunoprecipitation (RIP) and integrated transcriptome wide analyses reveal that CNOT7-regulated transcripts are enriched for a tripartite 3'UTR motif bound by RNA-binding proteins known to complex with CNOT7, TOB1, and CNOT1. Collectively, our data support a model of CNOT7, TOB1, CNOT1, and RNA-binding proteins collectively exerting post-transcriptional control on a metastasis suppressive transcriptional program to drive tumor cell metastasis.

An Integrated Genome-Wide Systems Genetics Screen for Breast Cancer Metastasis Susceptibility Genes.

Metastasis remains the primary cause of patient morbidity and mortality in solid tumors and is due to the action of a large number of tumor-autonomous and non-autonomous factors. Here we report the results of a genome-wide integrated strategy to identify novel metastasis susceptibility candidate genes and molecular pathways in breast cancer metastasis. This analysis implicates a number of transcriptional regulators and suggests cell-mediated immunity is an important determinant. Moreover, the analysis identified novel or FDA-approved drugs as potentially useful for anti-metastatic therapy. Further explorations implementing this strategy may therefore provide a variety of information for clinical applications in the control and treatment of advanced neoplastic disease.

A multi-megabase copy number gain causes maternal transmission ratio distortion on mouse chromosome 2.

Significant departures from expected Mendelian inheritance ratios (transmission ratio distortion, TRD) are frequently observed in both experimental crosses and natural populations. TRD on mouse Chromosome (Chr) 2 has been reported in multiple experimental crosses, including the Collaborative Cross (CC). Among the eight CC founder inbred strains, we found that Chr 2 TRD was exclusive to females that were heterozygous for the WSB/EiJ allele within a 9.3 Mb region (Chr 2 76.9 - 86.2 Mb). A copy number gain of a 127 kb-long DNA segment (designated as responder to drive, R2d) emerged as the strongest candidate for the causative allele. We mapped R2d sequences to two loci within the candidate interval. R2d1 is located near the proximal boundary, and contains a single copy of R2d in all strains tested. R2d2 maps to a 900 kb interval, and the number of R2d copies varies from zero in classical strains (including the mouse reference genome) to more than 30 in wild-derived strains. Using real-time PCR assays for the copy number, we identified a mutation (R2d2WSBdel1) that eliminates the majority of the R2d2WSB copies without apparent alterations of the surrounding WSB/EiJ haplotype. In a three-generation pedigree segregating for R2d2WSBdel1, the mutation is transmitted to the progeny and Mendelian segregation is restored in females heterozygous for R2d2WSBdel1, thus providing direct evidence that the copy number gain is causal for maternal TRD. We found that transmission ratios in R2d2WSB heterozygous females vary between Mendelian segregation and complete distortion depending on the genetic background, and that TRD is under genetic control of unlinked distorter loci. Although the R2d2WSB transmission ratio was inversely correlated with average litter size, several independent lines of evidence support the contention that female meiotic drive is the cause of the distortion. We discuss the implications and potential applications of this novel meiotic drive system.

An integrated systems genetics screen reveals the transcriptional structure of inherited predisposition to metastatic disease.

Metastasis is the result of stochastic genomic and epigenetic events leading to gene expression profiles that drive tumor dissemination. Here we exploit the principle that metastatic propensity is modified by the genetic background to generate prognostic gene expression signatures that illuminate regulators of metastasis. We also identify multiple microRNAs whose germline variation is causally linked to tumor progression and metastasis. We employ network analysis of global gene expression profiles in tumors derived from a panel of recombinant inbred mice to identify a network of co-expressed genes centered on Cnot2 that predicts metastasis-free survival. Modulating Cnot2 expression changes tumor cell metastatic potential in vivo, supporting a functional role for Cnot2 in metastasis. Small RNA sequencing of the same tumor set revealed a negative correlation between expression of the Mir216/217 cluster and tumor progression. Expression quantitative trait locus analysis (eQTL) identified cis-eQTLs at the Mir216/217 locus, indicating that differences in expression may be inherited. Ectopic expression of Mir216/217 in tumor cells suppressed metastasis in vivo. Finally, small RNA sequencing and mRNA expression profiling data were integrated to reveal that miR-3470a/b target a high proportion of network transcripts. In vivo analysis of Mir3470a/b demonstrated that both promote metastasis. Moreover, Mir3470b is a likely regulator of the Cnot2 network as its overexpression down-regulated expression of network hub genes and enhanced metastasis in vivo, phenocopying Cnot2 knockdown. The resulting data from this strategy identify Cnot2 as a novel regulator of metastasis and demonstrate the power of our systems-level approach in identifying modifiers of metastasis.

Host genetics influence tumour metastasis.

The complexity of the metastatic process has made it difficult to gain a full understanding of the origins of this most lethal aspect of cancer. Many factors probably have an important role, including somatic mutation, epigenetic modulations, interactions with normal stroma, and environmental stimuli. Additionally, recent evidence implies a significant role for germline polymorphisms in cancer progression. The existence of inherited metastasis risk factors (or prospective metastatic biomarkers) has potentially significant implications for our models of metastasis, clinical prognosis and the development of tailored treatment. Further investigations into the inherited components of metastasis might help resolve many of the questions that remain about tumour progression.

Cadm1 is a metastasis susceptibility gene that suppresses metastasis by modifying tumor interaction with the cell-mediated immunity.

Metastasis is a complex process utilizing both tumor-cell-autonomous properties and host-derived factors, including cellular immunity. We have previously shown that germline polymorphisms can modify tumor cell metastatic capabilities through cell-autonomous mechanisms. However, how metastasis susceptibility genes interact with the tumor stroma is incompletely understood. Here, we employ a complex genetic screen to identify Cadm1 as a novel modifier of metastasis. We demonstrate that Cadm1 can specifically suppress metastasis without affecting primary tumor growth. Unexpectedly, Cadm1 did not alter tumor-cell-autonomous properties such as proliferation or invasion, but required the host's adaptive immune system to affect metastasis. The metastasis-suppressing effect of Cadm1 was lost in mice lacking T cell-mediated immunity, which was partially phenocopied by depleting CD8(+) T cells in immune-competent mice. Our data show a novel function for Cadm1 in suppressing metastasis by sensitizing tumor cells to immune surveillance mechanisms, and this is the first report of a heritable metastasis susceptibility gene engaging tumor non-autonomous factors.

Allelic variation and differential expression of the mSIN3A histone deacetylase complex gene Arid4b promote mammary tumor growth and metastasis.

Accumulating evidence suggests that breast cancer metastatic progression is modified by germline polymorphism, although specific modifier genes have remained largely undefined. In the current study, we employ the MMTV-PyMT transgenic mouse model and the AKXD panel of recombinant inbred mice to identify AT-rich interactive domain 4B (Arid4b; NM_194262) as a breast cancer progression modifier gene. Ectopic expression of Arid4b promoted primary tumor growth in vivo as well as increased migration and invasion in vitro, and the phenotype was associated with polymorphisms identified between the AKR/J and DBA/2J alleles as predicted by our genetic analyses. Stable shRNA-mediated knockdown of Arid4b caused a significant reduction in pulmonary metastases, validating a role for Arid4b as a metastasis modifier gene. ARID4B physically interacts with the breast cancer metastasis suppressor BRMS1, and we detected differential binding of the Arid4b alleles to histone deacetylase complex members mSIN3A and mSDS3, suggesting that the mechanism of Arid4b action likely involves interactions with chromatin modifying complexes. Downregulation of the conserved Tpx2 gene network, which is comprised of many factors regulating cell cycle and mitotic spindle biology, was observed concomitant with loss of metastatic efficiency in Arid4b knockdown cells. Consistent with our genetic analysis and in vivo experiments in our mouse model system, ARID4B expression was also an independent predictor of distant metastasis-free survival in breast cancer patients with ER+ tumors. These studies support a causative role of ARID4B in metastatic progression of breast cancer.

Integrated cross-species transcriptional network analysis of metastatic susceptibility.

Metastatic disease is the proximal cause of mortality for most cancers and remains a significant problem for the clinical management of neoplastic disease. Recent advances in global transcriptional analysis have enabled better prediction of individuals likely to progress to metastatic disease. However, minimal overlap between predictive signatures has precluded easy identification of key biological processes contributing to the prometastatic transcriptional state. To overcome this limitation, we have applied network analysis to two independent human breast cancer datasets and three different mouse populations developed for quantitative analysis of metastasis. Analysis of these datasets revealed that the gene membership of the networks is highly conserved within and between species, and that these networks predicted distant metastasis free survival. Furthermore these results suggest that susceptibility to metastatic disease is cell-autonomous in estrogen receptor-positive tumors and associated with the mitotic spindle checkpoint. In contrast, nontumor genetics and pathway activities-associated stromal biology are significant modifiers of the rate of metastatic spread of estrogen receptor-negative tumors. These results suggest that the application of network analysis across species may provide a robust method to identify key biological programs associated with human cancer progression.

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.

SMARCD1 is a "Goldilocks" metastasis modifier.

Breast cancer is the most frequently diagnosed cancer worldwide, constituting around 15% of all diagnosed cancers in 2023. The predominant cause of breast cancer-related mortality is metastasis to distant essential organs, and a lack of metastasis-targeted therapies perpetuates dismal outcomes for late-stage patients. However, through our use of meiotic genetics to study inherited transcriptional network regulation, we have identified a new class of "Goldilocks" genes that are promising candidates for the development of metastasis-targeted therapeutics. Building upon previous work that implicated the CCR4-NOT RNA deadenylase complex in metastasis, we now demonstrate that the RNA-binding proteins (RNA-BPs) NANOS1, PUM2, and CPSF4 also regulate metastatic potential. Using cell lines, 3D culture, mouse models, and clinical data, we pinpoint Smarcd1 mRNA as a key target of all three RNA-BPs. Strikingly, both high and low expression of Smarcd1 is associated with positive clinical outcomes, while intermediate expression significantly reduces the probability of survival. Applying the theory of "essential genes" from evolution, we identify an additional 50 genes that span several cellular processes and must be maintained within a discrete window of expression for metastasis to occur. In the case of Smarcd1, small perturbations in its expression level significantly reduce metastasis in laboratory mouse models and alter splicing programs relevant to the ER+/HER2-enriched breast cancer subtype. The identification of subtype-specific "Goldilocks" metastasis modifier genes introduces a new class of genes and potential catalogue of novel targets that, when therapeutically "nudged" in either direction, may significantly improve late-stage patient outcomes.

Loss of NAT10 disrupts enhancer organization via p300 mislocalization and suppresses transcription of genes necessary for metastasis progression.

Acetylation of protein and RNA represent a critical event for development and cancer progression. NAT10 is the only known RNA acetylase that catalyzes the N4-actylcytidine (ac4C) modification of RNAs. Here, we show that the loss of NAT10 significantly decreases lung metastasis in allograft and genetically engineered mouse models of breast cancer. NAT10 interacts with a mechanosensitive, metastasis susceptibility protein complex at the nuclear pore. In addition to its canonical role in RNA acetylation, we find that NAT10 interacts with p300 at gene enhancers. NAT10 loss is associated with p300 mislocalization into heterochromatin regions. NAT10 depletion disrupts enhancer organization, leading to alteration of gene transcription necessary for metastatic progression, including reduced myeloid cell-recruiting chemokines that results in a less metastasis-prone tumor microenvironment. Our study uncovers a distinct role of NAT10 in enhancer organization of metastatic tumor cells and suggests its involvement in the tumor-immune crosstalk dictating metastatic outcomes.

Comparative analysis and classification of highly divergent mouse rDNA units based on their intergenic spacer (IGS) variability.

Ribosomal DNA (rDNA) repeat units are organized into tandem clusters in eukaryotic cells. In mice, these clusters are located on at least eight chromosomes and show extensive variation in the number of repeats between mouse genomes. To analyze intra- and inter-genomic variation of mouse rDNA repeats, we selectively isolated 25 individual rDNA units using Transformation-Associated Recombination (TAR) cloning. Long-read sequencing and subsequent comparative sequence analysis revealed that each full-length unit comprises an intergenic spacer (IGS) and a ∼13.4 kb long transcribed region encoding the three rRNAs, but with substantial variability in rDNA unit size, ranging from ∼35 to ∼46 kb. Within the transcribed regions of rDNA units, we found 209 variants, 70 of which are in external transcribed spacers (ETSs); but the rDNA size differences are driven primarily by IGS size heterogeneity, due to indels containing repetitive elements and some functional signals such as enhancers. Further evolutionary analysis categorized rDNA units into distinct clusters with characteristic IGS lengths; numbers of enhancers; and presence/absence of two common SNPs in promoter regions, one of which is located within promoter (p)RNA and may influence pRNA folding stability. These characteristic features of IGSs also correlated significantly with 5'ETS variant patterns described previously and associated with differential expression of rDNA units. Our results suggest that variant rDNA units are differentially regulated and open a route to investigate the role of rDNA variation on nucleolar formation and possible associations with pathology.

Transfer RNA acetylation regulates in vivo mammalian stress signaling.

Transfer RNA (tRNA) modifications are crucial for protein synthesis, but their position-specific physiological roles remain poorly understood. Here we investigate the impact of N4-acetylcytidine (ac4C), a highly conserved tRNA modification, using a Thumpd1 knockout mouse model. We find that loss of Thumpd1-dependent tRNA acetylation leads to reduced levels of tRNALeu, increased ribosome stalling, and activation of eIF2α phosphorylation. Thumpd1 knockout mice exhibit growth defects and sterility. Remarkably, concurrent knockout of Thumpd1 and the stress-sensing kinase Gcn2 causes penetrant postnatal lethality, indicating a critical genetic interaction. Our findings demonstrate that a modification restricted to a single position within type II cytosolic tRNAs can regulate ribosome-mediated stress signaling in mammalian organisms, with implications for our understanding of translation control as well as therapeutic interventions.

Integrative genomic identification of genes on 8p associated with hepatocellular carcinoma progression and patient survival.

BACKGROUND & AIMS: Hepatocellular carcinoma (HCC) is an aggressive malignancy; its mechanisms of development and progression are poorly understood. We used an integrative approach to identify HCC driver genes, defined as genes whose copy numbers associate with gene expression and cancer progression. METHODS: We combined data from high-resolution, array-based comparative genomic hybridization and transcriptome analysis of HCC samples from 76 patients with hepatitis B virus infection with data on patient survival times. Candidate genes were functionally validated using in vitro and in vivo models. RESULTS: Unsupervised analyses of array comparative genomic hybridization data associated loss of chromosome 8p with poor outcome (reduced survival time); somatic copy number alterations correlated with expression of 27.3% of genes analyzed. We associated expression levels of 10 of these genes with patient survival times in 2 independent cohorts (comprising 319 cases of HCC with mixed etiology) and 3 breast cancer cohorts (637 cases). Among the 10-gene signature, a cluster of 6 genes on 8p, (DLC1, CCDC25, ELP3, PROSC, SH2D4A, and SORBS3) were deleted in HCCs from patients with poor outcomes. In vitro and in vivo analyses indicated that the products of PROSC, SH2D4A, and SORBS3 have tumor-suppressive activities, along with the known tumor suppressor gene DLC1. CONCLUSIONS: We used an unbiased approach to identify 10 genes associated with HCC progression. These might be used in assisting diagnosis and to stage tumors based on gene expression patterns.