Gut microbiome plays a fundamental role in several aspects of host health and diseases. There has been an exponential surge in the use of animal models that can mimic different phenotypes of the human intestinal ecosystem. However, data on host species-specific signatures of gut microbiome and its metabolites like short-chain fatty acids (SCFAs; i.e., acetate, propionate, and butyrate) and lactate in these models and their similarities/differences from humans remain limited, due to high variability in protocols and analyses. Here, we analyze the fecal microbiota composition and the fecal levels of SCFAs and lactate in three of the most-widely used animal models, i.e., mice, rats, and non-human primates (NHPs) and compare them with human subjects, using data generated on a single platform with same protocols. The data show several species-specific similarities and differences in the gut microbiota and fecal organic acids between these species groups. Based on β-diversity, the gut microbiota in humans seems to be closer to NHPs than to mice and rats; however, among rodents, mice microbiota appears to be closer to humans than rats. The phylum-level analyses demonstrate higher Firmicutes–Bacteroidetes ratio in humans and NHPs vs. mice and rats. Human microbiota is dominated by Bacteroides followed by Ruminococcaceae and Clostridiales. Mouse gut is predominated by members of the family S24-7 followed by those from the order Clostridiales, whereas rats and NHPs have higher abundance of Prevotella compared with mice and humans. Also, fecal levels of lactate are higher in mice and rats vs. NHPs and humans, while acetate is highest in human feces. These data of host species-specific gut microbiota signatures in some of the most widely used animal models in context to the human microbiota might reflect disparities in host factors, e.g., diets, genetic origin, gender and age, and hence call for prospective studies investigating the features of gut microbiome in such animal models by controlling for these host elements.

Despite the obesity crisis in the United States, the underlying genetics are poorly understood. Our lab previously identified Keratinocyte-associated protein 3, Krtcap3, as a candidate gene for adiposity through a genome-wide association study in outbred rats, where increased liver expression of Krtcap3 correlated with decreased fat mass. Here we seek to confirm that Krtcap3 expression affects adiposity traits. To do so, we developed an in vivo whole-body Krtcap3 knock-out (KO) rat model. Wild-type (WT) and KO rats were placed onto a high-fat (HFD) or low-fat diet (LFD) at 6 weeks of age and were maintained on diet for 13 weeks, followed by assessments of metabolic health. We hypothesized that Krtcap3-KO rats will have increased adiposity and a worsened metabolic phenotype relative to WT. We found that KO male and female rats have significantly increased body weight versus WT, with the largest effect in females on a HFD. KO females also ate more and had greater adiposity, but were more insulin sensitive than WT regardless of diet condition. Although KO males weighed more than WT under both diet conditions, there were no differences in eating behavior or fat mass. Interestingly, KO males on a HFD were more insulin resistant than WT. This study confirms that Krtcap3 plays a role in body weight regulation and demonstrates genotype- and sex-specific effects on food intake, adiposity, and insulin sensitivity. Future studies will seek to better understand these sex differences, the role of diet, and establish a mechanism for Krtcap3 in obesity.

We previously identified Keratinocyte-associated protein 3, Krtcap3, as an obesity-related gene in female rats where a whole-body Krtcap3 knock-out (KO) led to increased adiposity compared to wild-type (WT) controls when fed a high-fat diet (HFD). We sought to replicate this work to better understand the function of Krtcap3 but were unable to reproduce the adiposity phenotype. In the current work, WT female rats ate more compared to WT in the prior study, with corresponding increases in body weight and fat mass, while there were no changes in these measures in KO females between the studies. The prior study was conducted before the COVID-19 pandemic, while the current study started after initial lock-down orders and was completed during the pandemic with a generally less stressful environment. We hypothesize that the environmental changes impacted stress levels and may explain the failure to replicate our results. Analysis of corticosterone (CORT) at euthanasia showed a significant study by genotype interaction where WT had significantly higher CORT relative to KO in Study 1, with no differences in Study 2. These data suggest that decreasing Krtcap3 expression may alter the environmental stress response to influence adiposity. We also found that KO rats in both studies, but not WT, experienced a dramatic increase in CORT after their cage mate was removed, suggesting a separate connection to social behavioral stress. Future work is necessary to confirm and elucidate the finer mechanisms of these relationships, but these data indicate the possibility of Krtcap3 as a novel stress gene.

Obesity is influenced by both genetics and diet and has wide ranging comorbidities, including anxiety and depressive disorders. In this study, outbred Heterogeneous Stock (HS) rats were fed high fat diet (HFD) or low fat diet (LFD) and tested for metabolic and behavioral outcomes. Male HS rats were fed 60% HFD or 10% LFD ad lib for 12-21 weeks. Rats were weighed weekly. We assessed multiple metabolic (glucose tolerance, fasting insulin, fat pad weights) and behavioral (elevated plus maze - EPM, open field test - OFT, splash test - SpT, and forced swim test - FST) measures. We found negative effects of HFD on metabolic outcomes, including increased body weight and fat pad weights, decreased glucose tolerance, and increased fasting insulin. We also found negative effects of HFD on coping response to stress and/or despair- and anxiety-like behaviors, including increased closed arm time and decreased open arm time in the EPM, increased movement and rest episodes and decreased rearing episodes in the OFT, increased SpT grooming time, and increased FST immobility. Positive correlations were found between FST immobility and visceral fat pad weights. We also found diet-dependent correlations between EPM, OFT and SpT measures. HFD negatively affected metabolic and behavioral health in HS rats. Changes in multiple coping/despair-and anxiety-like behaviors, independent of locomotion, indicate a role of HFD on behavioral health beyond just fat mass. This work establishes the HS rat as a model to study gene by diet interactions affecting obesity and behavioral health.

Transcriptomic analysis in metabolically active tissues allows a systems genetics approach to identify causal genes and networks involved in metabolic disease. Outbred heterogeneous stock (HS) rats are used for genetic mapping of complex traits, but to-date, a systems genetics analysis of metabolic tissues has not been done. We investigated whether adiposity-associated genes and gene co-expression networks in outbred heterogeneous stock (HS) rats overlap those found in humans. We analyzed RNAseq data from adipose tissue of 415 male HS rats, correlated these transcripts with body weight (BW) and compared transcriptome signatures to two human cohorts: the "African American Genetics of Metabolism and Expression" and "Metabolic Syndrome in Men". We used weighted gene co-expression network analysis to identify adiposity-associated gene networks and mediation analysis to identify genes under genetic control whose expression drives adiposity. We identified 554 orthologous "consensus genes" whose expression correlates with BW in the rat and with body mass index (BMI) in both human cohorts. Consensus genes fell within eight co-expressed networks and were enriched for genes involved in immune system function, cell growth, extracellular matrix organization and lipid metabolic processes. We identified 19 consensus genes for which genetic variation may influence BW via their expression, including those involved in lipolysis (e.g., Hcar1), inflammation (e.g., Rgs1), adipogenesis (e.g., Tmem120b) or no previously known role in obesity (e.g., St14, Msa4a6). Strong concordance between HS rat and human BW/BMI associated transcripts demonstrates translational utility of the rat model, while identification of novel genes expands our knowledge of the genetics underlying obesity.

Despite the successes of human genome-wide association studies, the causal genes underlying most metabolic traits remain unclear. We used outbred heterogeneous stock (HS) rats, coupled with expression data and mediation analysis, to identify quantitative trait loci (QTLs) and candidate gene mediators for adiposity, glucose tolerance, serum lipids, and other metabolic traits. Physiological traits were measured in 1519 male HS rats, with liver and adipose transcriptomes measured in over 410 rats. Genotypes were imputed from low coverage whole genome sequence. Linear mixed models were used to detect physiological and expression QTLs (pQTLs and eQTLs, respectively), employing both SNP- and haplotype-based models for pQTL mapping. Genes with cis-eQTLs that overlapped pQTLs were assessed as causal candidates through mediation analysis. We identified 14 SNP-based pQTLs and 19 haplotype-based pQTLs, of which 10 were in common. Using mediation, we identified the following genes as candidate mediators of pQTLs: Grk5 for a fat pad weight pQTL on Chr1, Krtcap3 for fat pad weight and serum lipids pQTLs on Chr6, Ilrun for a fat pad weight pQTL on Chr20 and Rfx6 for a whole pancreatic insulin content pQTL on Chr20. Furthermore, we verified Grk5 and Ktrcap3 using gene knock-down/out models, thereby shedding light on novel regulators of obesity.

Objective: Despite the obesity crisis in the United States, the underlying genetics are poorly understood. Our lab previously identified Keratinocyte-associated protein 3, Krtcap3, as a candidate gene for adiposity where increased expression of Krtcap3 correlated with decreased fat mass. Here we seek to confirm that Krtcap3 expression affects adiposity traits. Methods: We developed an in vivo whole-body Krtcap3 knock-out (KO) rat model. Wild-type (WT) and KO rats were placed onto a high-fat or low-fat diet at six weeks of age and were maintained on diet for 13 weeks, followed by assessments of metabolic health. We hypothesized that Krtcap3-KO rats will have increased adiposity and a worsened metabolic phenotype relative to WT. Results: We found that KO male and female rats have significantly increased body weight versus WT. KO females ate more, had more fat mass, but were also more insulin sensitive than WT. Alternatively, KO males weighed more and were more insulin resistant than WT, with no differences in eating or fat mass. Conclusions: This study validates Krtcap3 in body weight regulation and demonstrates sex-specific effects on food intake, adiposity, and insulin sensitivity. Future studies will investigate how Krtcap3 is acting and seek to better understand these sex differences.

Integration of retroviral elements into the host genome is a phenomena observed among many classes of retroviruses. Much information concerning integration of retroviral elements has been documented based on in vitro analysis or expression of selectable markers. To identify possible Tf1 integration events within silent regions of the S. pombe genome, we focused on performing an in vivo genome-wide analysis of Tf1 integration events from the nonselective phase of the retrotransposition assay. We analyzed 1000 individual colonies streaked from four independent Tf1 transposed patches under nonselection conditions. Our analysis detected a population of G418S/neo+ Tf1 integration events that would have been overlooked during the selective phase of the assay. Further RNA analysis from the G418S/neo+ clones revealed 50% of clones expressing the neo selectable marker. Our data reveals Tf1’s ability to insert within silent regions of S. pombe’s genome.