To identify genetic and environmental factors contributing to the pathogenesis of non-alcoholic fatty liver disease, we examined liver steatosis and related clinical and molecular traits in more than 100 unique inbred mouse strains, which were fed a diet rich in fat and carbohydrates. A >30-fold variation in hepatic TG accumulation was observed among the strains. Genome-wide association studies revealed three loci associated with hepatic TG accumulation. Utilizing transcriptomic data from the liver and adipose tissue, we identified several high-confidence candidate genes for hepatic steatosis, including Gde1, a glycerophosphodiester phosphodiesterase not previously implicated in triglyceride metabolism. We confirmed the role of Gde1 by in vivo hepatic over-expression and shRNA knockdown studies. We hypothesize that Gde1 expression increases TG production by contributing to the production of glycerol-3-phosphate. Our multi-level data, including transcript levels, metabolite levels, and gut microbiota composition, provide a framework for understanding genetic and environmental interactions underlying hepatic steatosis.

Evolutionary studies are often limited by missing data that are critical to understanding the history of selection. Selection experiments, which reproduce rapid evolution under controlled conditions, are excellent tools to study how genomes evolve under selection. Here we present a genomic dissection of the Longshanks selection experiment, in which mice were selectively bred over 20 generations for longer tibiae relative to body mass, resulting in 13% longer tibiae in two replicates. We synthesized evolutionary theory, genome sequences and molecular genetics to understand the selection response and found that it involved both polygenic adaptation and discrete loci of major effect, with the strongest loci tending to be selected in parallel between replicates. We show that selection may favor de-repression of bone growth through inactivating two limb enhancers of an inhibitor, Nkx3-2. Our integrative genomic analyses thus show that it is possible to connect individual base-pair changes to the overall selection response.

Hepatic metastases are a poor prognostic factor of colorectal carcinoma (CRC) and new strategies to reduce the risk of liver CRC colonization are highly needed. Herein, we used mouse models of hepatic metastatization to demonstrate that the continuous infusion of therapeutic doses of interferon-alpha (IFNα) controls CRC invasion by acting on hepatic endothelial cells (HECs). Mechanistically, IFNα promoted the development of a vascular antimetastatic niche characterized by liver sinusoidal endothelial cells (LSECs) defenestration extracellular matrix and glycocalyx deposition, thus strengthening the liver vascular barrier impairing CRC trans-sinusoidal migration, without requiring a direct action on tumor cells, hepatic stellate cells, hepatocytes, or liver dendritic cells (DCs), Kupffer cells (KCs) and liver capsular macrophages (LCMs). Moreover, IFNα endowed LSECs with efficient cross-priming potential that, along with the early intravascular tumor burden reduction, supported the generation of antitumor CD8+ T cells and ultimately led to the establishment of a protective long-term memory T cell response. These findings provide a rationale for the use of continuous IFNα therapy in perioperative settings to reduce CRC metastatic spreading to the liver.

Genetic variation conferring resistance and susceptibility to carcinogen-induced tumorigenesis is frequently studied in mice. We have now turned this idea to melanoma using the collaborative cross (CC), a resource of mouse strains designed to discover genes for complex diseases. We studied melanoma-prone transgenic progeny across seventy CC genetic backgrounds. We mapped a strong quantitative trait locus for rapid onset spontaneous melanoma onset to Prkdc, a gene involved in detection and repair of DNA damage. In contrast, rapid onset UVR-induced melanoma was linked to the ribosomal subunit gene Rrp15. Ribosome biogenesis was upregulated in skin shortly after UVR exposure. Mechanistically, variation in the 'usual suspects' by which UVR may exacerbate melanoma, defective DNA repair, melanocyte proliferation, or inflammatory cell infiltration, did not explain melanoma susceptibility or resistance across the CC. Instead, events occurring soon after exposure, such as dysregulation of ribosome function, which alters many aspects of cellular metabolism, may be important.

Sequence variation in enhancers that control cell-type-specific gene transcription contributes significantly to phenotypic variation within human populations. However, it remains difficult to predict precisely the effect of any given sequence variant on enhancer function due to the complexity of DNA sequence motifs that determine transcription factor (TF) binding to enhancers in their native genomic context. Using F1-hybrid cells derived from crosses between distantly related inbred strains of mice, we identified thousands of enhancers with allele-specific TF binding and/or activity. We find that genetic variants located within the central region of enhancers are most likely to alter TF binding and enhancer activity. We observe that the AP-1 family of TFs (Fos/Jun) are frequently required for binding of TEAD TFs and for enhancer function. However, many sequence variants outside of core motifs for AP-1 and TEAD also impact enhancer function, including sequences flanking core TF motifs and AP-1 half sites. Taken together, these data represent one of the most comprehensive assessments of allele-specific TF binding and enhancer function to date and reveal how sequence changes at enhancers alter their function across evolutionary timescales.