Colorectal Cancer Stem Cell Subtypes Orchestrate Distinct Tumor Microenvironments.

Several classification systems have been developed to define tumor subtypes in colorectal cancer (CRC). One system proposes that tumor heterogeneity derives in part from distinct cancer stem cell populations that co-exist as admixtures of varying proportions. However, the lack of single cell resolution has prohibited a definitive identification of these types of stem cells and therefore any understanding of how each influence tumor phenotypes. Here were report the isolation and characterization of two cancer stem cell subtypes from the SW480 CRC cell line. We find these cancer stem cells are oncogenic versions of the normal Crypt Base Columnar (CBC) and Regenerative Stem Cell (RSC) populations from intestinal crypts and that their gene signatures are consistent with the "Admixture" and other CRC classification systems. Using publicly available single cell RNA sequencing (scRNAseq) data from CRC patients, we determine that RSC and CBC cancer stem cells are commonly co-present in human CRC. To characterize influences on the tumor microenvironment, we develop subtype-specific xenograft models and we define their tumor microenvironments at high resolution via scRNAseq. RSCs create differentiated, inflammatory, slow growing tumors. CBCs create proliferative, undifferentiated, invasive tumors. With this enhanced resolution, we unify current CRC patient classification schema with TME phenotypes and organization.

Circadian control of tumor immunosuppression affects efficacy of immune checkpoint blockade.

The circadian clock is a critical regulator of immunity, and this circadian control of immune modulation has an essential function in host defense and tumor immunosurveillance. Here we use a single-cell RNA sequencing approach and a genetic model of colorectal cancer to identify clock-dependent changes to the immune landscape that control the abundance of immunosuppressive cells and consequent suppression of cytotoxic CD8+ T cells. Of these immunosuppressive cell types, PD-L1-expressing myeloid-derived suppressor cells (MDSCs) peak in abundance in a rhythmic manner. Disruption of the epithelial cell clock regulates the secretion of cytokines that promote heightened inflammation, recruitment of neutrophils and the subsequent development of MDSCs. We also show that time-of-day anti-PD-L1 delivery is most effective when synchronized with the abundance of immunosuppressive MDSCs. Collectively, these data indicate that circadian gating of tumor immunosuppression informs the timing and efficacy of immune checkpoint inhibitors.

Leveraging inter-individual transcriptional correlation structure to infer discrete signaling mechanisms across metabolic tissues.

Inter-organ communication is a vital process to maintain physiologic homeostasis, and its dysregulation contributes to many human diseases. Given that circulating bioactive factors are stable in serum, occur naturally, and are easily assayed from blood, they present obvious focal molecules for therapeutic intervention and biomarker development. Recently, studies have shown that secreted proteins mediating inter-tissue signaling could be identified by 'brute force' surveys of all genes within RNA-sequencing measures across tissues within a population. Expanding on this intuition, we reasoned that parallel strategies could be used to understand how individual genes mediate signaling across metabolic tissues through correlative analyses of gene variation between individuals. Thus, comparison of quantitative levels of gene expression relationships between organs in a population could aid in understanding cross-organ signaling. Here, we surveyed gene-gene correlation structure across 18 metabolic tissues in 310 human individuals and 7 tissues in 103 diverse strains of mice fed a normal chow or high-fat/high-sucrose (HFHS) diet. Variation of genes such as FGF21, ADIPOQ, GCG, and IL6 showed enrichments which recapitulate experimental observations. Further, similar analyses were applied to explore both within-tissue signaling mechanisms (liver PCSK9) and genes encoding enzymes producing metabolites (adipose PNPLA2), where inter-individual correlation structure aligned with known roles for these critical metabolic pathways. Examination of sex hormone receptor correlations in mice highlighted the difference of tissue-specific variation in relationships with metabolic traits. We refer to this resource as gene-derived correlations across tissues (GD-CAT) where all tools and data are built into a web portal enabling users to perform these analyses without a single line of code (gdcat.org). This resource enables querying of any gene in any tissue to find correlated patterns of genes, cell types, pathways, and network architectures across metabolic organs.

Loss of CTRP10 results in female obesity with preserved metabolic health.

Obesity is a major risk factor for type 2 diabetes, dyslipidemia, cardiovascular disease, and hypertension. Intriguingly, there is a subset of metabolically healthy obese (MHO) individuals who are seemingly able to maintain a healthy metabolic profile free of metabolic syndrome. The molecular underpinnings of MHO, however, are not well understood. Here, we report that CTRP10/C1QL2-deficient mice represent a unique female model of MHO. CTRP10 modulates weight gain in a striking and sexually dimorphic manner. Female, but not male, mice lacking CTRP10 develop obesity with age on a low-fat diet while maintaining an otherwise healthy metabolic profile. When fed an obesogenic diet, female Ctrp10 knockout (KO) mice show rapid weight gain. Despite pronounced obesity, Ctrp10 KO female mice do not develop steatosis, dyslipidemia, glucose intolerance, insulin resistance, oxidative stress, or low-grade inflammation. Obesity is largely uncoupled from metabolic dysregulation in female KO mice. Multi-tissue transcriptomic analyses highlighted gene expression changes and pathways associated with insulin-sensitive obesity. Transcriptional correlation of the differentially expressed gene (DEG) orthologous in humans also show sex differences in gene connectivity within and across metabolic tissues, underscoring the conserved sex-dependent function of CTRP10. Collectively, our findings suggest that CTRP10 negatively regulates body weight in females, and that loss of CTRP10 results in benign obesity with largely preserved insulin sensitivity and metabolic health. This female MHO mouse model is valuable for understanding sex-biased mechanisms that uncouple obesity from metabolic dysfunction.

Systems genetics approaches for understanding complex traits with relevance for human disease.

Quantitative traits are often complex because of the contribution of many loci, with further complexity added by environmental factors. In medical research, systems genetics is a powerful approach for the study of complex traits, as it integrates intermediate phenotypes, such as RNA, protein, and metabolite levels, to understand molecular and physiological phenotypes linking discrete DNA sequence variation to complex clinical and physiological traits. The primary purpose of this review is to describe some of the resources and tools of systems genetics in humans and rodent models, so that researchers in many areas of biology and medicine can make use of the data.

Feeding neurons integrate metabolic and reproductive states in mice.

Balance between metabolic and reproductive processes is important for survival, particularly in mammals that gestate their young. How the nervous system coordinates this balance is an active area of study. Herein, we demonstrate that somatostatin (SST) neurons of the tuberal hypothalamus alter feeding in a manner sensitive to metabolic and reproductive states in mice. Whereas chemogenetic activation of SST neurons increased food intake across sexes, ablation decreased food intake only in female mice during proestrus. This ablation effect was only apparent in animals with low body mass. Fat transplantation and bioinformatics analysis of SST neuronal transcriptomes revealed white adipose as a key modulator of these effects. These studies indicate that SST hypothalamic neurons integrate metabolic and reproductive cues by responding to varying levels of circulating estrogens to modulate feeding differentially based on energy stores. Thus, gonadal steroid modulation of neuronal circuits can be context dependent and gated by metabolic status.

CTRP13 ablation improves systemic glucose and lipid metabolism.

OBJECTIVE: Tissue crosstalk mediated by secreted hormones underlies the integrative control of metabolism. We previously showed that CTRP13/C1QL3, a secreted protein of the C1q family, can improve glucose metabolism and insulin action in vitro and reduce food intake and body weight in mice when centrally delivered. A role for CTRP13 in regulating insulin secretion in isolated islets has also been demonstrated. It remains unclear, however, whether the effects of CTRP13 on cultured cells and in mice reflect the physiological function of the protein. Here, we use a loss-of-function mouse model to address whether CTRP13 is required for metabolic homeostasis. METHODS: WT and Ctrp13 knockout (KO) mice fed a standard chow or a high-fat diet were subjected to comprehensive metabolic phenotyping. Transcriptomic analyses were carried out on visceral and subcutaneous fat, liver, and skeletal muscle to identify pathways altered by CTRP13 deficiency. RNA-seq data was further integrated with the Metabolic Syndrome in Man (METSIM) cohort data. Adjusted regression analysis was used to demonstrate that genetic variation of CTRP13 expression accounts for a significant proportion of variance between differentially expressed genes (DEGs) in adipose tissue and metabolic traits in humans. RESULTS: Contrary to expectation, chow-fed Ctrp13-KO male mice had elevated physical activity, lower body weight, and improved lipid handling. On a high-fat diet (HFD), Ctrp13-KO mice of either sex were consistently more active and leaner. Loss of CTRP13 reduced hepatic glucose output and improved glucose tolerance, insulin sensitivity, and triglyceride clearance, though with notable sex differences. Consistent with the lean phenotype, transcriptomic analyses revealed a lower inflammatory profile in visceral fat and liver. Reduced hepatic steatosis was correlated with the suppression of lipid synthesis and enhanced lipid catabolism gene expression. Visceral fat had the largest number of DEGs and mediation analyses on the human orthologs of the DEGs suggested the potential causal contribution of CTRP13 to human metabolic syndrome. CONCLUSIONS: Our results suggest that CTRP13 is a negative metabolic regulator, and its deficiency improves systemic metabolic profiles. Our data also suggest the reduction in circulating human CTRP13 levels seen in obesity and diabetes may reflect a compensatory physiologic response to counteract insulin resistance.

Genetic architecture of heart mitochondrial proteome influencing cardiac hypertrophy.

Mitochondria play an important role in both normal heart function and disease etiology. We report analysis of common genetic variations contributing to mitochondrial and heart functions using an integrative proteomics approach in a panel of inbred mouse strains called the Hybrid Mouse Diversity Panel (HMDP). We performed a whole heart proteome study in the HMDP (72 strains, n=2-3 mice) and retrieved 848 mitochondrial proteins (quantified in ≥50 strains). High-resolution association mapping on their relative abundance levels revealed three trans-acting genetic loci on chromosomes (chr) 7, 13 and 17 that regulate distinct classes of mitochondrial proteins as well as cardiac hypertrophy. DAVID enrichment analyses of genes regulated by each of the loci revealed that the chr13 locus was highly enriched for complex-I proteins (24 proteins, P=2.2E-61), the chr17 locus for mitochondrial ribonucleoprotein complex (17 proteins, P=3.1E-25) and the chr7 locus for ubiquinone biosynthesis (3 proteins, P=6.9E-05). Follow-up high resolution regional mapping identified NDUFS4, LRPPRC and COQ7 as the candidate genes for chr13, chr17 and chr7 loci, respectively, and both experimental and statistical analyses supported their causal roles. Furthermore, a large cohort of Diversity Outbred mice was used to corroborate Lrpprc gene as a driver of mitochondrial DNA (mtDNA)-encoded gene regulation, and to show that the chr17 locus is specific to heart. Variations in all three loci were associated with heart mass in at least one of two independent heart stress models, namely, isoproterenol-induced heart failure and diet-induced obesity. These findings suggest that common variations in certain mitochondrial proteins can act in trans to influence tissue-specific mitochondrial functions and contribute to heart hypertrophy, elucidating mechanisms that may underlie genetic susceptibility to heart failure in human populations.

Conserved multi-tissue transcriptomic adaptations to exercise training in humans and mice.

Physical activity is associated with beneficial adaptations in human and rodent metabolism. We studied over 50 complex traits before and after exercise intervention in middle-aged men and a panel of 100 diverse strains of female mice. Candidate gene analyses in three brain regions, muscle, liver, heart, and adipose tissue of mice indicate genetic drivers of clinically relevant traits, including volitional exercise volume, muscle metabolism, adiposity, and hepatic lipids. Although ∼33% of genes differentially expressed in skeletal muscle following the exercise intervention are similar in mice and humans independent of BMI, responsiveness of adipose tissue to exercise-stimulated weight loss appears controlled by species and underlying genotype. We leveraged genetic diversity to generate prediction models of metabolic trait responsiveness to volitional activity offering a framework for advancing personalized exercise prescription. The human and mouse data are publicly available via a user-friendly Web-based application to enhance data mining and hypothesis development.

The potential of integrating human and mouse discovery platforms to advance our understanding of cardiometabolic diseases.

Cardiometabolic diseases encompass a range of interrelated conditions that arise from underlying metabolic perturbations precipitated by genetic, environmental, and lifestyle factors. While obesity, dyslipidaemia, smoking, and insulin resistance are major risk factors for cardiometabolic diseases, individuals still present in the absence of such traditional risk factors, making it difficult to determine those at greatest risk of disease. Thus, it is crucial to elucidate the genetic, environmental, and molecular underpinnings to better understand, diagnose, and treat cardiometabolic diseases. Much of this information can be garnered using systems genetics, which takes population-based approaches to investigate how genetic variance contributes to complex traits. Despite the important advances made by human genome-wide association studies (GWAS) in this space, corroboration of these findings has been hampered by limitations including the inability to control environmental influence, limited access to pertinent metabolic tissues, and often, poor classification of diseases or phenotypes. A complementary approach to human GWAS is the utilisation of model systems such as genetically diverse mouse panels to study natural genetic and phenotypic variation in a controlled environment. Here, we review mouse genetic reference panels and the opportunities they provide for the study of cardiometabolic diseases and related traits. We discuss how the post-GWAS era has prompted a shift in focus from discovery of novel genetic variants to understanding gene function. Finally, we highlight key advantages and challenges of integrating complementary genetic and multi-omics data from human and mouse populations to advance biological discovery.

Dysregulated systemic metabolism in a Down syndrome mouse model.

OBJECTIVE: Trisomy 21 is one of the most complex genetic perturbations compatible with postnatal survival. Dosage imbalance arising from the triplication of genes on human chromosome 21 (Hsa21) affects multiple organ systems. Much of Down syndrome (DS) research, however, has focused on addressing how aneuploidy dysregulates CNS function leading to cognitive deficit. Although obesity, diabetes, and associated sequelae such as fatty liver and dyslipidemia are well documented in the DS population, only limited studies have been conducted to determine how gene dosage imbalance affects whole-body metabolism. Here, we conduct a comprehensive and systematic analysis of key metabolic parameters across different physiological states in the Ts65Dn trisomic mouse model of DS. METHODS: Ts65Dn mice and euploid littermates were subjected to comprehensive metabolic phenotyping under basal (chow-fed) state and the pathophysiological state of obesity induced by a high-fat diet (HFD). RNA sequencing of liver, skeletal muscle, and two major fat depots were conducted to determine the impact of aneuploidy on tissue transcriptome. Pathway enrichments, gene-centrality, and key driver estimates were performed to provide insights into tissue autonomous and non-autonomous mechanisms contributing to the dysregulation of systemic metabolism. RESULTS: Under the basal state, chow-fed Ts65Dn mice of both sexes had elevated locomotor activity and energy expenditure, reduced fasting serum cholesterol levels, and mild glucose intolerance. Sexually dimorphic deterioration in metabolic homeostasis became apparent when mice were challenged with a high-fat diet. While obese Ts65Dn mice of both sexes exhibited dyslipidemia, male mice also showed impaired systemic insulin sensitivity, reduced mitochondrial activity, and elevated fibrotic and inflammatory gene signatures in the liver and adipose tissue. Systems-level analysis highlighted conserved pathways and potential endocrine drivers of adipose-liver crosstalk that contribute to dysregulated glucose and lipid metabolism. CONCLUSIONS: A combined alteration in the expression of trisomic and disomic genes in peripheral tissues contribute to metabolic dysregulations in Ts65Dn mice. These data lay the groundwork for understanding the impact of aneuploidy on in vivo metabolism.

Lipopolysaccharide-induced chronic inflammation increases female serum gonadotropins and shifts the pituitary transcriptomic landscape.

Introduction: Female reproductive function depends on a choreographed sequence of hormonal secretion and action, where specific stresses such as inflammation exert profound disruptions. Specifically, acute LPS-induced inflammation inhibits gonadotropin production and secretion from the pituitary, thereby impacting the downstream production of sex hormones. These outcomes have only been observed in acute inflammatory stress and little is known about the mechanisms by which chronic inflammation affects reproduction. In this study we seek to understand the chronic effects of LPS on pituitary function and consequent luteinizing and follicle stimulating hormone secretion. Methods: A chronic inflammatory state was induced in female mice by twice weekly injections with LPS over 6 weeks. Serum gonadotropins were measured and bulk RNAseq was performed on the pituitaries from these mice, along with basic measurements of reproductive biology. Results: Surprisingly, serum luteinizing and follicle stimulating hormone was not inhibited and instead we found it was increased with repeated LPS treatments. Discussion: Analysis of bulk RNA-sequencing of murine pituitary revealed paracrine activation of TGFß pathways as a potential mechanism regulating FSH secretion in response to chronic LPS. These results provide a framework with which to begin dissecting the impacts of chronic inflammation on reproductive physiology.

Proteome-wide systems genetics identifies UFMylation as a regulator of skeletal muscle function.

Improving muscle function has great potential to improve the quality of life. To identify novel regulators of skeletal muscle metabolism and function, we performed a proteomic analysis of gastrocnemius muscle from 73 genetically distinct inbred mouse strains, and integrated the data with previously acquired genomics and >300 molecular/phenotypic traits via quantitative trait loci mapping and correlation network analysis. These data identified thousands of associations between protein abundance and phenotypes and can be accessed online (https://muscle.coffeeprot.com/) to identify regulators of muscle function. We used this resource to prioritize targets for a functional genomic screen in human bioengineered skeletal muscle. This identified several negative regulators of muscle function including UFC1, an E2 ligase for protein UFMylation. We show UFMylation is up-regulated in a mouse model of amyotrophic lateral sclerosis, a disease that involves muscle atrophy. Furthermore, in vivo knockdown of UFMylation increased contraction force, implicating its role as a negative regulator of skeletal muscle function.

Liver-heart cross-talk mediated by coagulation factor XI protects against heart failure.

Tissue-tissue communication by endocrine factors is a vital mechanism for physiologic homeostasis. A systems genetics analysis of transcriptomic and functional data from a cohort of diverse, inbred strains of mice predicted that coagulation factor XI (FXI), a liver-derived protein, protects against diastolic dysfunction, a key trait of heart failure with preserved ejection fraction. This was confirmed using gain- and loss-of-function studies, and FXI was found to activate the bone morphogenetic protein (BMP)-SMAD1/5 pathway in the heart. The proteolytic activity of FXI is required for the cleavage and activation of extracellular matrix-associated BMP7 in the heart, thus inhibiting genes involved in inflammation and fibrosis. Our results reveal a protective role of FXI in heart injury that is distinct from its role in coagulation.

Tryptophan metabolism is a physiological integrator regulating circadian rhythms.

OBJECTIVE: The circadian clock aligns physiology with the 24-hour rotation of Earth. Light and food are the main environmental cues (zeitgebers) regulating circadian rhythms in mammals. Yet, little is known about the interaction between specific dietary components and light in coordinating circadian homeostasis. Herein, we focused on the role of essential amino acids. METHODS: Mice were fed diets depleted of specific essential amino acids and their behavioral rhythms were monitored and tryptophan was selected for downstream analyses. The role of tryptophan metabolism in modulating circadian homeostasis was studied using isotope tracing as well as transcriptomic- and metabolomic- analyses. RESULTS: Dietary tryptophan depletion alters behavioral rhythms in mice. Furthermore, tryptophan metabolism was shown to be regulated in a time- and light- dependent manner. A multi-omics approach and combinatory diet/light interventions demonstrated that tryptophan metabolism modulates temporal regulation of metabolism and transcription programs by buffering photic cues. Specifically, tryptophan metabolites regulate central circadian functions of the suprachiasmatic nucleus and the core clock machinery in the liver. CONCLUSIONS: Tryptophan metabolism is a modulator of circadian homeostasis by integrating environmental cues. Our findings propose tryptophan metabolism as a potential point for pharmacologic intervention to modulate phenotypes associated with disrupted circadian rhythms.

Disruption of the circadian clock drives Apc loss of heterozygosity to accelerate colorectal cancer.

An alarming rise in young onset colorectal cancer (CRC) has been reported; however, the underlying molecular mechanism remains undefined. Suspected risk factors of young onset CRC include environmental aspects, such as lifestyle and dietary factors, which are known to affect the circadian clock. We find that both genetic disruption and environmental disruption of the circadian clock accelerate Apc-driven CRC pathogenesis in vivo. Using an intestinal organoid model, we demonstrate that clock disruption promotes transformation by driving Apc loss of heterozygosity, which hyperactivates Wnt signaling. This up-regulates c-Myc, a known Wnt target, which drives heightened glycolytic metabolism. Using patient-derived organoids, we show that circadian rhythms are lost in human tumors. Last, we identify that variance between core clock and Wnt pathway genes significantly predicts the survival of patients with CRC. Overall, our findings demonstrate a previously unidentified mechanistic link between clock disruption and CRC, which has important implications for young onset cancer prevention.

Sex differences in heart mitochondria regulate diastolic dysfunction.

Heart failure with preserved ejection fraction (HFpEF) exhibits a sex bias, being more common in women than men, and we hypothesize that mitochondrial sex differences might underlie this bias. As part of genetic studies of heart failure in mice, we observe that heart mitochondrial DNA levels and function tend to be reduced in females as compared to males. We also observe that expression of genes encoding mitochondrial proteins are higher in males than females in human cohorts. We test our hypothesis in a panel of genetically diverse inbred strains of mice, termed the Hybrid Mouse Diversity Panel (HMDP). Indeed, we find that mitochondrial gene expression is highly correlated with diastolic function, a key trait in HFpEF. Consistent with this, studies of a "two-hit" mouse model of HFpEF confirm that mitochondrial function differs between sexes and is strongly associated with a number of HFpEF traits. By integrating data from human heart failure and the mouse HMDP cohort, we identify the mitochondrial gene Acsl6 as a genetic determinant of diastolic function. We validate its role in HFpEF using adenoviral over-expression in the heart. We conclude that sex differences in mitochondrial function underlie, in part, the sex bias in diastolic function.

Clonally expanded CD8 T cells characterize amyotrophic lateral sclerosis-4.

Amyotrophic lateral sclerosis (ALS) is a heterogenous neurodegenerative disorder that affects motor neurons and voluntary muscle control1. ALS heterogeneity includes the age of manifestation, the rate of progression and the anatomical sites of symptom onset. Disease-causing mutations in specific genes have been identified and define different subtypes of ALS1. Although several ALS-associated genes have been shown to affect immune functions2, whether specific immune features account for ALS heterogeneity is poorly understood. Amyotrophic lateral sclerosis-4 (ALS4) is characterized by juvenile onset and slow progression3. Patients with ALS4 show motor difficulties by the time that they are in their thirties, and most of them require devices to assist with walking by their fifties. ALS4 is caused by mutations in the senataxin gene (SETX). Here, using Setx knock-in mice that carry the ALS4-causative L389S mutation, we describe an immunological signature that consists of clonally expanded, terminally differentiated effector memory (TEMRA) CD8 T cells in the central nervous system and the blood of knock-in mice. Increased frequencies of antigen-specific CD8 T cells in knock-in mice mirror the progression of motor neuron disease and correlate with anti-glioma immunity. Furthermore, bone marrow transplantation experiments indicate that the immune system has a key role in ALS4 neurodegeneration. In patients with ALS4, clonally expanded TEMRA CD8 T cells circulate in the peripheral blood. Our results provide evidence of an antigen-specific CD8 T cell response in ALS4, which could be used to unravel disease mechanisms and as a potential biomarker of disease state.