Variation in histone configurations correlates with gene expression across nine inbred strains of mice.

The Diversity Outbred (DO) mice and their inbred founders are widely used models of human disease. However, although the genetic diversity of these mice has been well documented, their epigenetic diversity has not. Epigenetic modifications, such as histone modifications and DNA methylation, are important regulators of gene expression and, as such, are a critical mechanistic link between genotype and phenotype. Therefore, creating a map of epigenetic modifications in the DO mice and their founders is an important step toward understanding mechanisms of gene regulation and the link to disease in this widely used resource. To this end, we performed a strain survey of epigenetic modifications in hepatocytes of the DO founders. We surveyed four histone modifications (H3K4me1, H3K4me3, H3K27me3, and H3K27ac), as well as DNA methylation. We used ChromHMM to identify 14 chromatin states, each of which represents a distinct combination of the four histone modifications. We found that the epigenetic landscape is highly variable across the DO founders and is associated with variation in gene expression across strains. We found that epigenetic state imputed into a population of DO mice recapitulated the association with gene expression seen in the founders, suggesting that both histone modifications and DNA methylation are highly heritable mechanisms of gene expression regulation. We illustrate how DO gene expression can be aligned with inbred epigenetic states to identify putative cis-regulatory regions. Finally, we provide a data resource that documents strain-specific variation in the chromatin state and DNA methylation in hepatocytes across nine widely used strains of laboratory mice.

Patch grafting, strategies for transplantation of organoids into solid organs such as liver.

Epithelial cell therapies have been at an impasse because of inefficient methods of transplantation to solid organs. Patch grafting strategies were established enabling transplantation of ≥107th organoids/patch of porcine GFP+ biliary tree stem/progenitors into livers of wild type hosts. Grafts consisted of organoids embedded in soft (~100 Pa) hyaluronan hydrogels, both prepared in serum-free Kubota's Medium; placed against target sites; covered with a silk backing impregnated with more rigid hyaluronan hydrogels (~700 Pa); and use of the backing to tether grafts with sutures or glue to target sites. Hyaluronan coatings (~200-300 Pa) onto the serosal surface of the graft served to minimize adhesions with neighboring organs. The organ's clearance of hyaluronans enabled restoration of tissue-specific paracrine and systemic signaling, resulting in return of normal hepatic histology, with donor parenchymal cells uniformly integrated amidst host cells and that had differentiated to mature hepatocytes and cholangiocytes. Grafts containing donor mature hepatocytes, partnered with endothelia, and in the same graft biomaterials as for stem/progenitor organoids, did not engraft. Engraftment occurred if porcine liver-derived mesenchymal stem cells (MSCs) were co-transplanted with donor mature cells. RNA-seq analyses revealed that engraftment correlated with expression of matrix-metalloproteinases (MMPs), especially secreted isoforms that were found expressed strongly by organoids, less so by MSCs, and minimally, if at all, by adult cells. Engraftment with patch grafting strategies occurred without evidence of emboli or ectopic cell distribution. It was successful with stem/progenitor organoids or with cells with a source(s) of secreted MMP isoforms and offers significant potential for enabling cell therapies for solid organs.

A Thalamic Orphan Receptor Drives Variability in Short-Term Memory.

Working memory is a form of short-term memory that involves maintaining and updating task-relevant information toward goal-directed pursuits. Classical models posit persistent activity in prefrontal cortex (PFC) as a primary neural correlate, but emerging views suggest additional mechanisms may exist. We screened ∼200 genetically diverse mice on a working memory task and identified a genetic locus on chromosome 5 that contributes to a substantial proportion (17%) of the phenotypic variance. Within the locus, we identified a gene encoding an orphan G-protein-coupled receptor, Gpr12, which is sufficient to drive substantial and bidirectional changes in working memory. Molecular, cellular, and imaging studies revealed that Gpr12 enables high thalamus-PFC synchrony to support memory maintenance and choice accuracy. These findings identify an orphan receptor as a potent modifier of short-term memory and supplement classical PFC-based models with an emerging thalamus-centric framework for the mechanistic understanding of working memory.

microRNA-146a-5p association with the cardiometabolic disease risk factor TMAO.

Trimethylamine-N-oxide (TMAO), a microbial choline metabolism byproduct that is processed in the liver and excreted into circulation, is associated with increased atherosclerotic lesion formation and cardiovascular disease risk. Genetic regulators of TMAO levels are largely unknown. In the present study, we used 288 mice from a genetically heterogeneous mouse population [Diversity Outbred (DO)] to determine hepatic microRNA associations with TMAO in the context of an atherogenic diet. We also validated findings in two additional animal models of atherosclerosis: liver-specific insulin receptor knockout mice fed a chow diet (LIRKO) and African green monkeys fed high-fat/high-cholesterol diet. Small RNA-sequencing analysis in DO mice, LIRKO mice, and African green monkeys identified only one hepatic microRNA (miR-146a-5p) that is aberrantly expressed across all three models. Moreover, miR-146a-5p levels are associated with circulating TMAO after atherogenic diet in each of these models. We also performed high-resolution genetic mapping and identified a novel quantitative trait locus on Chromosome 12 for TMAO levels. This interval includes two genes, Numb and Dlst, which are inversely correlated with both miR-146a and TMAO and are predicted targets of miR-146a. Both of these genes have been validated as direct targets of miR-146a, though in other cellular contexts. This is the first report to our knowledge of a link between miR-146 and TMAO. Our findings suggest that miR-146-5p, as well as one or more genes at the Chromosome 12 QTL (possibly Numb or Dlst), is strongly linked to TMAO levels and likely involved in the control of atherosclerosis.

Colonic epithelial miR-31 associates with the development of Crohn's phenotypes.

BACKGROUND: Crohn's disease (CD) is highly heterogeneous, due in large part to variability in cellular processes that underlie the natural history of CD, thereby confounding effective therapy. There is a critical need to advance understanding of the cellular mechanisms that drive CD heterogeneity. METHODS: We performed small RNA sequencing of adult colon tissue from CD and NIBD controls. Colonic epithelial cells and immune cells were isolated from colonic tissues, and microRNA-31 (miR-31) expression was measured. miR-31 expression was measured in colonoid cultures generated from controls and patients with CD. We performed small RNA-sequencing of formalin-fixed paraffin-embedded colon and ileum biopsies from treatment-naive pediatric patients with CD and controls and collected data on disease features and outcomes. RESULTS: Small RNA-sequencing and microRNA profiling in the colon revealed 2 distinct molecular subtypes, each with different clinical associations. Notably, we found that miR-31 expression was a driver of these 2 subtypes and, further, that miR-31 expression was particularly pronounced in epithelial cells. Colonoids revealed that miR-31 expression differences are preserved in this ex vivo system. In adult patients, low colonic miR-31 expression levels at the time of surgery were associated with worse disease outcome as measured by need for an end ileostomy and recurrence of disease in the neoterminal ileum. In pediatric patients, lower miR-31 expression at the time of diagnosis was associated with future development of fibrostenotic ileal CD requiring surgeryCONCLUSIONS. These findings represent an important step forward in designing more effective clinical trials and developing personalized CD therapies. FUNDING: This work was supported by CCF Career Development Award (SZS), R01-ES024983 from NIEHS (SZS and TSF), 1R01DK104828-01A1 from NIDDK (SZS and TSF), P01-DK094779-01A1 from NIDDK (SZS), P30-DK034987 from NIDDK (SZS), 1-16-ACE-47 ADA Pathway Award (PS), UNC Nutrition Obesity Research Center Pilot & Feasibility Grant P30DK056350 (PS), CCF PRO-KIIDS NETWORK (SZS and PS), UNC CGIBD T32 Training Grant from NIDDK (JBB), T32 Training Grant (5T32GM007092-42) from NIGMS (MH), and SHARE from the Helmsley Trust (SZS). The UNC Translational Pathology Laboratory is supported, in part, by grants from the National Cancer Institute (3P30CA016086) and the UNC University Cancer Research Fund (UCRF) (PS).

Comprehensive analysis of The Cancer Genome Atlas reveals a unique gene and non-coding RNA signature of fibrolamellar carcinoma.

Fibrolamellar carcinoma (FLC) is a unique liver cancer primarily affecting young adults and characterized by a fusion event between DNAJB1 and PRKACA. By analyzing RNA-sequencing data from The Cancer Genome Atlas (TCGA) for >9,100 tumors across ~30 cancer types, we show that the DNAJB1-PRKACA fusion is specific to FLCs. We demonstrate that FLC tumors (n = 6) exhibit distinct messenger RNA (mRNA) and long intergenic non-coding RNA (lincRNA) profiles compared to hepatocellular carcinoma (n = 263) and cholangiocarcinoma (n = 36), the two most common liver cancers. We also identify a set of mRNAs (n = 16) and lincRNAs (n = 4), including LINC00473, that distinguish FLC from ~25 other liver and non-liver cancer types. We confirm this unique FLC signature by analysis of two independent FLC cohorts (n = 20 and 34). Lastly, we validate the overexpression of one specific gene in the FLC signature, carbonic anhydrase XII (CA12), at the protein level by western blot and immunohistochemistry. Both the mRNA and lincRNA signatures support a major role for protein kinase A (PKA) signaling in shaping the FLC gene expression landscape, and present novel candidate FLC oncogenes that merit further investigation.

Functional Transcriptomics in Diverse Intestinal Epithelial Cell Types Reveals Robust MicroRNA Sensitivity in Intestinal Stem Cells to Microbial Status.

Gut microbiota play an important role in regulating the development of the host immune system, metabolic rate, and at times, disease pathogenesis. The factors and mechanisms that mediate interactions between microbiota and the intestinal epithelium are not fully understood. We provide novel evidence that microbiota may control intestinal epithelial stem cell (IESC) proliferation in part through microRNAs (miRNAs). We demonstrate that miRNA profiles differ dramatically across functionally distinct cell types of the mouse jejunal intestinal epithelium and that miRNAs respond to microbiota in a highly cell type-specific manner. Importantly, we also show that miRNAs in IESCs are more prominently regulated by microbiota compared with miRNAs in any other intestinal epithelial cell subtype. We identify miR-375 as one miRNA that is significantly suppressed by the presence of microbiota in IESCs. Using a novel method to knockdown gene and miRNA expression ex vivo enteroids, we demonstrate that we can knock down gene expression in Lgr5+ IESCs. Furthermore, when we knock down miR-375 in IESCs, we observe significantly increased proliferative capacity. Understanding the mechanisms by which microbiota regulate miRNA expression in IESCs and other intestinal epithelial cell subtypes will elucidate a critical molecular network that controls intestinal homeostasis and, given the heightened interest in miRNA-based therapies, may offer novel therapeutic strategies in the treatment of gastrointestinal diseases associated with altered IESC function.

Quantitative trait locus mapping of genes that regulate HDL cholesterol in SM/J and NZB/B1NJ inbred mice.

To investigate the quantitative trait loci (QTL) regulating plasma cholesterol, the female progeny of an (SMxNZB/ B1NJ)xNZB/B1NJ backcross were fed an atherogenic diet. After 18 wk, plasma total cholesterol and high-density lipoprotein cholesterol (HDL-C) was measured. HDL-C concentrations were greater in NZB than in SM mice. For standard chow-fed mice, QTL were found near D5Mit370 and D18Mit34. For mice fed an atherogenic diet, a QTL was found near D5Mit239. The QTL for chow-fed and atherogenic-fed mice on chromosome 5 seem to be two different loci. We used a multitrait analysis to rule out pleiotropy in favor of a two-QTL hypothesis. Furthermore, the HDL-C in these strains was induced by the high-fat diet. For inducible HDL-C, one significant locus was found near D15Mit39. The gene for an HDL receptor, Srb1, maps close to the HDL-C QTL at D5Mit370, but the concentrations of Srb1 mRNA and SR-B1 protein and the gene sequence of NZB/B1NJ and SM/J did not support Srb1 as a candidate gene. With these QTL, we have identified chromosomal regions that affect lipoprotein profiles in these strains.