Background Carotid artery intima/media thickness (IMT) is a hallmark trait associated with future cardiovascular events. The goal of this study was to map new genes that regulate carotid IMT by genome-wide association. Methods and Results We induced IMT by ligation procedure of the left carotid artery in 30 inbred mouse strains. Histologic reconstruction revealed significant variation in left carotid artery intima, media, adventitia, external elastic lamina volumes, intima-to-media ratio, and (intima+media)/external elastic lamina percent ratio in inbred mice. The carotid remodeling trait was regulated by distinct genomic signatures with a dozen common single-nucleotide polymorphisms associated with left carotid artery intima volume, intima-to-media ratio, and (intima+media)/external elastic lamina percent ratio. Among genetic loci on mouse chromosomes 1, 4, and 12, there was natriuretic peptide receptor 2 (Npr2), a strong candidate gene. We observed that only male, not female, mice heterozygous for a targeted Npr2 deletion (Npr2+/-) exhibited defective carotid artery remodeling compared with Npr2 wild-type (Npr2+/+) littermates. Fibrosis in carotid IMT was significantly increased in Npr2+/- males compared with Npr2+/- females or Npr2+/+ mice. We also detected decreased Npr2 expression in human atherosclerotic plaques, similar to that seen in studies in Npr2+/- mice. Conclusions We found that components of carotid IMT were regulated by distinct genetic factors. We also showed a critical role for Npr2 in genetic regulation of vascular fibrosis associated with defective carotid remodeling.

A traditional approach to investigate the genetic basis of complex diseases is to identify genes with a global change in expression between diseased and healthy individuals. However, population heterogeneity may undermine the effort to uncover genes with significant but individual contribution to the spectrum of disease phenotypes within a population. Here we investigate individual changes of gene expression when inducing hypertrophy and heart failure in 100 + strains of genetically distinct mice from the Hybrid Mouse Diversity Panel (HMDP). We find that genes whose expression fold-change correlates in a statistically significant way with the severity of the disease are either up or down-regulated across strains, and therefore missed by a traditional population-wide analysis of differential gene expression. Furthermore, those "fold-change" genes are enriched in human cardiac disease genes and form a dense co-regulated module strongly interacting with the cardiac hypertrophic signaling network in the human interactome. We validate our approach by showing that the knockdown of Hes1, predicted as a strong candidate, induces a dramatic reduction of hypertrophy by 80-90% in neonatal rat ventricular myocytes. Our results demonstrate that individualized approaches are crucial to identify genes underlying complex diseases as well as to develop personalized therapies.

Liver fibrosis is a multifactorial trait that develops in response to chronic liver injury. Our aim was to characterize the genetic architecture of carbon tetrachloride (CCl4)-induced liver fibrosis using the Hybrid Mouse Diversity Panel, a panel of more than 100 genetically distinct mouse strains optimized for genome-wide association studies and systems genetics.

Elevated plasma cholesterol and type 2 diabetes (T2D) are associated with coronary artery disease (CAD). Individuals treated with cholesterol-lowering statins have increased T2D risk, while individuals with hypercholesterolemia have reduced T2D risk. We explore the relationship between lipid and glucose control by constructing network models from the STARNET study with sequencing data from seven cardiometabolic tissues obtained from CAD patients during coronary artery by-pass grafting surgery. By integrating gene expression, genotype, metabolomic, and clinical data, we identify a glucose and lipid determining (GLD) regulatory network showing inverse relationships with lipid and glucose traits. Master regulators of the GLD network also impact lipid and glucose levels in inverse directions. Experimental inhibition of one of the GLD network master regulators, lanosterol synthase (LSS), in mice confirms the inverse relationships to glucose and lipid levels as predicted by our model and provides mechanistic insights.

Lipocalin-2 (LCN2) is a secreted protein involved in innate immunity and has also been associated with several cardiometabolic traits in both mouse and human studies. However, the causal relationship of LCN2 to these traits is unclear, and most studies have examined only males.

Cardiovascular diseases (CVD) and type 2 diabetes (T2D) are closely interrelated complex diseases likely sharing overlapping pathogenesis driven by aberrant activities in gene networks. However, the molecular circuitries underlying the pathogenic commonalities remain poorly understood. We sought to identify the shared gene networks and their key intervening drivers for both CVD and T2D by conducting a comprehensive integrative analysis driven by five multi-ethnic genome-wide association studies (GWAS) for CVD and T2D, expression quantitative trait loci (eQTLs), ENCODE, and tissue-specific gene network models (both co-expression and graphical models) from CVD and T2D relevant tissues. We identified pathways regulating the metabolism of lipids, glucose, and branched-chain amino acids, along with those governing oxidation, extracellular matrix, immune response, and neuronal system as shared pathogenic processes for both diseases. Further, we uncovered 15 key drivers including HMGCR, CAV1, IGF1 and PCOLCE, whose network neighbors collectively account for approximately 35% of known GWAS hits for CVD and 22% for T2D. Finally, we cross-validated the regulatory role of the top key drivers using in vitro siRNA knockdown, in vivo gene knockout, and two Hybrid Mouse Diversity Panels each comprised of >100 strains. Findings from this in-depth assessment of genetic and functional data from multiple human cohorts provide strong support that common sets of tissue-specific molecular networks drive the pathogenesis of both CVD and T2D across ethnicities and help prioritize new therapeutic avenues for both CVD and T2D.

Signaling pathways that promote adipose tissue thermogenesis are well characterized, but the limiters of energy expenditure are largely unknown. Here, we show that ablation of the anti-inflammatory cytokine IL-10 improves insulin sensitivity, protects against diet-induced obesity, and elicits the browning of white adipose tissue. Mechanistic studies define bone marrow cells as the source of the IL-10 signal and adipocytes as the target cell type mediating these effects. IL-10 receptor alpha is highly enriched in mature adipocytes and is induced in response to differentiation, obesity, and aging. Assay for transposase-accessible chromatin sequencing (ATAC-seq), ChIP-seq, and RNA-seq reveal that IL-10 represses the transcription of thermogenic genes in adipocytes by altering chromatin accessibility and inhibiting ATF and C/EBPβ recruitment to key enhancer regions. These findings expand our understanding of the relationship between inflammatory signaling pathways and adipose tissue function and provide insight into the physiological control of thermogenesis that could inform future therapy.

The upregulated expression of heparin binding EGF-like growth factor (HB-EGF) in the vessel and circulation is associated with risk of cardiovascular disease. In this study, we tested the effects of HB-EGF targeting using HB-EGF-specific antisense oligonucleotide (ASO) on the development of aortic aneurysm in a mouse aneurysm model.

Mouse models have been used extensively to study human coronary artery disease (CAD) or atherosclerosis and to test therapeutic targets. However, whether mouse and human share similar genetic factors and pathogenic mechanisms of atherosclerosis has not been thoroughly investigated in a data-driven manner. We conducted a cross-species comparison study to better understand atherosclerosis pathogenesis between species by leveraging multiomics data. Specifically, we compared genetically driven and thus CAD-causal gene networks and pathways, by using human GWAS of CAD from the CARDIoGRAMplusC4D consortium and mouse GWAS of atherosclerosis from the Hybrid Mouse Diversity Panel (HMDP) followed by integration with functional multiomics human (STARNET and GTEx) and mouse (HMDP) databases. We found that mouse and human shared >75% of CAD causal pathways. Based on network topology, we then predicted key regulatory genes for both the shared pathways and species-specific pathways, which were further validated through the use of single cell data and the latest CAD GWAS. In sum, our results should serve as a much-needed guidance for which human CAD-causal pathways can or cannot be further evaluated for novel CAD therapies using mouse models.

To elucidate the contributions of specific lipid species to metabolic traits, we integrated global hepatic lipid data with other omics measures and genetic data from a cohort of about 100 diverse inbred strains of mice fed a high-fat/high-sucrose diet for 8 weeks. Association mapping, correlation, structure analyses, and network modeling revealed pathways and genes underlying these interactions. In particular, our studies lead to the identification of Ifi203 and Map2k6 as regulators of hepatic phosphatidylcholine homeostasis and triacylglycerol accumulation, respectively. Our analyses highlight mechanisms for how genetic variation in hepatic lipidome can be linked to physiological and molecular phenotypes, such as microbiota composition.

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.

Coronary artery disease (CAD) is a complex, multifactorial disease caused, in particular, by inflammation and cholesterol metabolism. At the molecular level, the role of tissue-specific signaling pathways leading to CAD is still largely unexplored. This study relied on two main resources: (1) genes with impact on atherosclerosis/CAD, and (2) liver-specific transcriptome analyses from human and mouse studies. The transcription factor activating transcription factor 3 (ATF3) was identified as a key regulator of a liver network relevant to atherosclerosis and linked to inflammation and cholesterol metabolism. ATF3 was predicted to be a direct and indirect (via MAF BZIP Transcription Factor F (MAFF)) regulator of low-density lipoprotein receptor (LDLR). Chromatin immunoprecipitation DNA sequencing (ChIP-seq) data from human liver cells revealed an ATF3 binding motif in the promoter regions of MAFF and LDLR. siRNA knockdown of ATF3 in human Hep3B liver cells significantly upregulated LDLR expression (p < 0.01). Inflammation induced by lipopolysaccharide (LPS) stimulation resulted in significant upregulation of ATF3 (p < 0.01) and subsequent downregulation of LDLR (p < 0.001). Liver-specific expression data from human CAD patients undergoing coronary artery bypass grafting (CABG) surgery (STARNET) and mouse models (HMDP) confirmed the regulatory role of ATF3 in the homeostasis of cholesterol metabolism. This study suggests that ATF3 might be a promising treatment candidate for lowering LDL cholesterol and reducing cardiovascular risk.

Non-alcoholic fatty liver disease (NAFLD) is a liver manifestation of metabolic syndrome, and is estimated to affect one billion individuals worldwide. An increased intake of a high-fat diet (HFD) and sugar-sweetened beverages are risk-factors for NAFLD development, but how their combined intake promotes progression to a more severe form of liver injury is unknown. Here we show that fructose metabolism via ketohexokinase (KHK) C isoform increases endoplasmic reticulum (ER) stress in a dose dependent fashion, so when fructose is coupled with a HFD intake it leads to unresolved ER stress. Conversely, a liver-specific knockdown of KHK in C57BL/6J male mice consuming fructose on a HFD is adequate to improve the NAFLD activity score and exert a profound effect on the hepatic transcriptome. Overexpression of KHK-C in cultured hepatocytes is sufficient to induce ER stress in fructose free media. Upregulation of KHK-C is also observed in genetically obesity ob/ob, db/db and lipodystrophic FIRKO male mice, whereas KHK knockdown in these mice improves metabolic function. Additionally, in over 100 inbred strains of male or female mice hepatic KHK expression correlates positively with adiposity, insulin resistance, and liver triglycerides. Similarly, in 241 human subjects and their controls, hepatic Khk expression is upregulated in early, but not late stages of NAFLD. In summary, we describe a novel role of KHK-C in triggering ER stress, which offers a mechanistic understanding of how the combined intake of fructose and a HFD propagates the development of metabolic complications.

The molecular bases of how host genetic variation impacts the gut microbiome remain largely unknown. Here we used a genetically diverse mouse population and applied systems genetics strategies to identify interactions between host and microbe phenotypes including microbial functions, using faecal metagenomics, small intestinal transcripts and caecal lipids that influence microbe-host dynamics. Quantitative trait locus (QTL) mapping identified murine genomic regions associated with variations in bacterial taxa; bacterial functions including motility, sporulation and lipopolysaccharide production and levels of bacterial- and host-derived lipids. We found overlapping QTL for the abundance of Akkermansia muciniphila and caecal levels of ornithine lipids. Follow-up in vitro and in vivo studies revealed that A. muciniphila is a major source of these lipids in the gut, provided evidence that ornithine lipids have immunomodulatory effects and identified intestinal transcripts co-regulated with these traits including Atf3, which encodes for a transcription factor that plays vital roles in modulating metabolism and immunity. Collectively, these results suggest that ornithine lipids are potentially important for A. muciniphila-host interactions and support the role of host genetics as a determinant of responses to gut microbes.

The gut microbiome influences nutrient processing as well as host physiology. Plasma lipid levels have been associated with the microbiome, although the underlying mechanisms are largely unknown, and the effects of dietary lipids on the gut microbiome in humans are not well-studied. We used a compilation of four studies utilizing non-human primates (Chlorocebus aethiops and Macaca fascicularis) with treatments that manipulated plasma lipid levels using dietary and pharmacological techniques, and characterized the microbiome using 16S rDNA. High-fat diets significantly reduced alpha diversity (Shannon) and the Firmicutes/Bacteroidetes ratio compared to chow diets, even when the diets had different compositions and were applied in different orders. When analyzed for differential abundance using DESeq2, Bulleidia, Clostridium, Ruminococcus, Eubacterium, Coprocacillus, Lachnospira, Blautia, Coprococcus, and Oscillospira were greater in both chow diets while Succinivibrio, Collinsella, Streptococcus, and Lactococcus were greater in both high-fat diets (oleic blend or lard fat source). Dietary cholesterol levels did not affect the microbiome and neither did alterations of plasma lipid levels through treatments of miR-33 antisense oligonucleotide (anti-miR-33), Niemann-Pick C1-Like 1 (NPC1L1) antisense oligonucleotide (ASO), and inducible degrader of LDLR (IDOL) ASO. However, a liver X receptor (LXR) agonist shifted the microbiome and decreased bile acid levels. Fifteen genera increased with the LXR agonist, while seven genera decreased. Pseudomonas increased on the LXR agonist and was negatively correlated to deoxycholic acid, cholic acid, and total bile acids while Ruminococcus was positively correlated with taurolithocholic acid and taurodeoxycholic acid. Seven of the nine bile acids identified in the feces significantly decreased due to the LXR agonist, and total bile acids (nmol/g) was reduced by 62%. These results indicate that plasma lipid levels have, at most, a modest effect on the microbiome, whereas bile acids, derived in part from plasma lipids, are likely responsible for the indirect relationship between lipid levels and the microbiome.

Conductances of ion channels and transporters controlling cardiac excitation may vary in a population of subjects with different cardiac gene expression patterns. However, the amount of variability and its origin are not quantitatively known. We propose a new conceptual approach to predict this variability that consists of finding combinations of conductances generating a normal intracellular Ca2+ transient without any constraint on the action potential. Furthermore, we validate experimentally its predictions using the Hybrid Mouse Diversity Panel, a model system of genetically diverse mouse strains that allows us to quantify inter-subject versus intra-subject variability. The method predicts that conductances of inward Ca2+ and outward K+ currents compensate each other to generate a normal Ca2+ transient in good quantitative agreement with current measurements in ventricular myocytes from hearts of different isogenic strains. Our results suggest that a feedback mechanism sensing the aggregate Ca2+ transient of the heart suffices to regulate ionic conductances.

The recent advent of 3D photography has created the potential for comprehensive facial evaluation. However, lack of practical true 3D analysis of the information collected from 3D images has been the factor limiting widespread utilization in orthodontics. Current evaluation of 3D facial soft tissue images relies on subjective visual evaluation and 2D distances to assess facial disharmony. The objectives of this project strive to map the surface and define boundaries of 3D facial soft tissue, modify mathematical functions to average multiple 3D facial images, and mathematically average 3D facial images allowing generation of color-coded surface deviation relative to a true average.

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.

Diabetic nephropathy (DN) is a major complication of diabetes and the leading cause of end-stage renal disease. DN is characterized by changes in kidney structure and function but the underlying genetic and molecular factors are poorly understood. We used a mouse diversity panel to explore the genetic basis of DN traits in mice carrying the Ins2 Akita mutation. Twenty-eight Akita strains were generated by breeding this panel to DBA/2.Akita mice. Male F1 diabetic and nondiabetic littermates were evaluated for DN-related traits. Urine albumin-to-creatinine ratios (ACRs), volume and cystatin C as well as blood urea nitrogen and lipoprotein levels varied significantly among the diabetic strains. For most Akita strains, ACR values increased 2- to 6-fold over euglycemic control values. However, six strains exhibited changes in ACR exceeding 10-fold with two strains (NOD/ShiLt and CBA) showing 50- to 83- fold increases. These increases are larger than previously reported among available DN mouse models establishing these strains as useful for additional studies of renal function. ACRs correlated with cystatin C (P = 0.0286), a measure of hyperfiltration and an interstitial tubular marker associated with DN onset in humans suggesting that tubule damage as well as podocyte-stress contributed to reduced kidney function assessed by ACR. Although large changes were seen for ACRs, severe nephropathology was absent. However, glomerular hypertrophy and collagen IV content were found to vary significantly among strains suggesting a genetic basis for early onset features of DN. Our results define the range of DN phenotypes that occur among common inbred strains of mice.

Recent studies indicate both clinical and mechanistic links between atherosclerotic heart disease and intestinal microbial metabolism of certain dietary nutrients producing trimethylamine N-oxide (TMAO). Here we test the hypothesis that gut microbial transplantation can transmit choline diet-induced TMAO production and atherosclerosis susceptibility. First, a strong association was noted between atherosclerotic plaque and plasma TMAO levels in a mouse diversity panel (n = 22 strains, r = 0.38; p = 0.0001). An atherosclerosis-prone and high TMAO-producing strain, C57BL/6J, and an atherosclerosis-resistant and low TMAO-producing strain, NZW/LacJ, were selected as donors for cecal microbial transplantation into apolipoprotein e null mice in which resident intestinal microbes were first suppressed with antibiotics. Trimethylamine (TMA) and TMAO levels were initially higher in recipients on choline diet that received cecal microbes from C57BL/6J inbred mice; however, durability of choline diet-dependent differences in TMA/TMAO levels was not maintained to the end of the study. Mice receiving C57BL/6J cecal microbes demonstrated choline diet-dependent enhancement in atherosclerotic plaque burden as compared with recipients of NZW/LacJ microbes. Microbial DNA analyses in feces and cecum revealed transplantation of donor microbial community features into recipients with differences in taxa proportions between donor strains that were transmissible to recipients and that tended to show coincident proportions with TMAO levels. Proportions of specific taxa were also identified that correlated with plasma TMAO levels in donors and recipients and with atherosclerotic lesion area in recipients. Atherosclerosis susceptibility may be transmitted via transplantation of gut microbiota. Gut microbes may thus represent a novel therapeutic target for modulating atherosclerosis susceptibility.