Gene expression in hypothalamus and brown adipose tissue of mice divergently selected for heat loss

Comprehensive Analysis of NAFLD and the Therapeutic Target Identified

Objective: Non-alcoholic fatty liver disease (NAFLD) is a serious health threat worldwide. The aim of this study was to comprehensively describe the metabolic and immunologic characteristics of NAFLD, and to explore potential therapeutic drug targets for NAFLD. Methods: Six NAFLD datasets were downloaded from the Gene Expression Omnibus (GEO) database, including GSE48452, GSE63067, GSE66676, GSE89632, GSE24807, and GSE37031. The datasets we then used to identify and analyze genes that were differentially expressed in samples from patients with NAFLD and normal subjects, followed by analysis of the metabolic and immunologic characteristics of patients with NAFLD. We also identified potential therapeutic drugs for NAFLD using the Connectivity Map (CMAP) database. Moreover, we constructed a prediction model using minimum depth random forest analysis and screened for potential therapeutic targets. Finally, therapeutic targets were verified in a fatty liver model stimulated by palmitic acid (PA). Results: A total of 1,358 differentially expressed genes (DEGs) were obtained, which were mainly enriched in carbohydrate metabolism, lipid metabolism, and other metabolic pathways. Immune infiltration analysis showed that memory B cells, regulatory T cells and M1 macrophage were significantly up-regulated, while T cells follicular helper were down regulated in NAFLD. These may provide a reference for the immune-metabolism interaction in the pathogenesis of NAFLD. Digoxin and helveticoside were identified as potential therapeutic drugs for NAFLD via the CMAP database. In addition, a five-gene prediction model based on minimum depth random forest analysis was constructed, and the receiver operating characteristic (ROC) curves of both training and validation set reached 1. The five candidate therapeutic targets were ENO3, CXCL10, INHBE, LRRC31, and OPTN. Moreover, the efficiency of hepatocyte adipogenesis decreased after OPTN knockout, confirming the potential use of OPTN as a new therapeutic target for NAFLD. Conclusion: This study provides a deeper insight into the molecular pathogenesis of NAFLD. We used five key genes to construct a diagnostic model with a strong predictive effect. Therefore, these five key genes may play an important role in the diagnosis and treatment of NAFLD, particularly those with increased OPTN expression.

Diet-induced obesity causes ghrelin resistance in arcuate NPY/AgRP neurons

Circulating ghrelin is decreased in obesity, and peripheral ghrelin does not induce food intake in obese mice. We investigated whether ghrelin resistance was a centrally mediated phenomenon involving dysregulated neuropeptide Y (NPY) and agouti-related peptide (AgRP) circuits. We show that diet-induced obesity (DIO) (12 wk) suppresses the neuroendocrine ghrelin system by decreasing acylated and total plasma ghrelin, decreasing ghrelin and Goat mRNA in the stomach, and decreasing expression of hypothalamic GHSR. Peripheral (ip) or central (intracerebroventricular) ghrelin injection was able to induce food intake and arcuate nucleus Fos immunoreactivity in chow-fed but not high-fat diet-fed mice. DIO decreased expression of Npy and Agrp mRNA, and central ghrelin was unable to promote expression of these genes. Ghrelin did not induce AgRP or NPY secretion in hypothalamic explants from DIO mice. Injection of NPY intracerebroventricularly increased food intake in both chow-fed and high-fat diet-fed mice, indicating that downstream NPY/AgRP neural targets are intact and that defective NPY/AgRP function is a primary cause of ghrelin resistance. Ghrelin resistance in DIO is not confined to the NPY/AgRP neurons, because ghrelin did not stimulate growth hormone secretion in DIO mice. Collectively, our data suggests that DIO causes ghrelin resistance by reducing NPY/AgRP responsiveness to plasma ghrelin and suppressing the neuroendocrine ghrelin axis to limit further food intake. Ghrelin has a number of functions in the brain aside from appetite control, including cognitive function, mood regulation, and protecting against neurodegenerative diseases. Thus, central ghrelin resistance may potentiate obesity-related cognitive decline, and restoring ghrelin sensitivity may provide therapeutic outcomes for maintaining healthy aging.

Exploration of genes showing intramuscular fat deposition-associated expression changes in musculus longissimus muscle

Marbling, as defined by the amount of intramuscular fat, is an economically important trait in beef cattle. Intramuscular fat deposition is postulated to arise mainly from a series of adipogenic events in intramuscular adipocyte-lineage cells and in the physiological or anatomical milieux surrounding them. This study was designed to investigate gene-expression patterns associated with fat deposition in musculus longissimus muscle, including adipocyte-lineage cells and part of the milieux. Differential-display PCR (ddPCR) was used to examine expression differences between low-marbled and high-marbled steer groups at 8, 10, 12 and 14 months of age, encompassing the time that marbling starts to appear. Seventy-four of 2114 total bands on ddPCR gel-bands were significant (P < 0.05) for the group effect, the interaction effect between group and age, or both the group and the interaction effects. Sequence analysis of 72 of these bands revealed 77 genes, including 35 annotated genes and 42 novel sequences. Among the 35 annotated genes, 6 (BTG2, PDHB, SORBS1, TRDN, TTN and MGP) have been related to changes in intramuscular fat deposition, possibly by exerting effects on adipocyte-lineage cells or on the milieux surrounding them.

Characterization of QTL with major effects on fatness and growth on mouse chromosome 2

OBJECTIVE

To isolate and characterize a region on mouse chromosome 2 harboring quantitative trait loci with large influences on growth and fatness.

RESEARCH METHODS AND PROCEDURES

A congenic line [M16i.B6-(D2Mit306-D2Mit52); MB2] was created using the polygenic obese M16i selection line as the recipient for an approximately 38-centimorgan region from C57BL/6J. Males and females from M16i and MB2 were compared for body weight, body composition, feed consumption, and additional traits at 6, 15, and 24 weeks. Interactions of genotype and environment (low and high dietary fat) were investigated. Males (8 weeks) were evaluated for fatty acid profiles in liver and for transcriptional profiles in liver and adipose.

RESULTS

Consequences of replacing M16i alleles with C57BL/6J alleles in MB2 were maximized at 15 weeks. MB2 mice were up to 15% lighter than M16i at this age, with no differences in feed consumption. As a percentage of body weight, MB2 had dramatically less epididymal (males) or perimetrial (females) fat (1.17% vs. 2.79% pooled across sex) and lower total lipids (16.1% vs. 23.3%) than M16i. Decreased adiposity in MB2 was not dependent on gender or diet. MB2 mice also had significant decreases in levels of leptin, insulin, and glucose, decreased de novo synthesis of hepatic fatty acid, and transcriptional changes for many genes both within, and external to, the congenic region.

DISCUSSION

Results confirm the presence and large effects of mouse chromosome 2 quantitative trait loci and further define their phenotypic consequences related to energy balance. The MB2 congenic line is a powerful resource for eventual identification of pathways and mutations within genes regulating predisposition to growth and obesity.

Evaluation of gene expression in pigs selected for enhanced reproduction using differential display PCR: II. Anterior pituitary

The objective of this study was to identify differentially expressed genes in the anterior pituitary (AP) of sows selected for enhanced reproductive phenotypes. Selection in the Index (I) line was based on an index of ovulation rate and embryo survival, whereas random selection was used in the Control (C) line. Average numbers of fully formed piglets at birth were 12.5 +/- 1.5 and 9.9 +/- 2.0 for Line I and C sows used in this study, respectively. In order to induce luteolysis and synchronize follicle development, sows were injected (i.m.) with 2 mL of prostaglandin F2alpha analog between d 12 and 14 of the estrous cycle. Tissue was harvested 2 d (d2) or 4 d (d4) after injection, resulting in four experimental groups: Cd2 (n = 6), Cd4 (n = 4), Id2 (n = 6), and Id4 (n = 7). Differential display PCR (ddPCR) was used to search for transcriptional changes between selection lines in the AP, using samples within line but pooled across days. Northern hybridization was used to confirm ddPCR results. For ddPCR, two pools were used from each line (C and I). Three genes were confirmed to be differentially expressed between Lines I and C: G-beta like protein, ferritin heavy-chain, and follicle stimulating hormone beta subunit, whereas many other expressed sequence tags were observed to be differentially expressed but still require confirmation. Our findings indicate that long-term selection to increase ovulation rate and decrease embryo mortality has altered transcriptional patterns in the anterior pituitary, most likely as correlated responses.

Physical activity and resting metabolic rate

The direct effects of physical activity interventions on energy expenditure are relatively small when placed in the context of total daily energy demands. Hence, the suggestion has been made that exercise produces energetic benefits in other components of the daily energy budget, thus generating a net effect on energy balance much greater than the direct energy cost of the exercise alone. Resting metabolic rate (RMR) is the largest component of the daily energy budget in most human societies and, therefore, any increases in RMR in response to exercise interventions are potentially of great importance. Animal studies have generally shown that single exercise events and longer-term training produce increases in RMR. This effect is observed in longer-term interventions despite parallel decreases in body mass and fat mass. Flight is an exception, as both single flights and long-term flight training induce reductions in RMR. Studies in animals that measure the effect of voluntary exercise regimens on RMR are less commonly performed and do not show the same response as that to forced exercise. In particular, they indicate that exercise does not induce elevations in RMR. Many studies of human subjects indicate a short-term elevation in RMR in response to single exercise events (generally termed the excess post-exercise O2 consumption; EPOC). This EPOC appears to have two phases, one lasting < 2 h and a smaller much more prolonged effect lasting up to 48 h. Many studies have shown that long-term training increases RMR, but many other studies have failed to find such effects. Data concerning long-term effects of training are potentially confounded by some studies not leaving sufficient time after the last exercise bout for the termination of the long-term EPOC. Long-term effects of training include increases in RMR due to increases in lean muscle mass. Extreme interventions, however, may induce reductions in RMR, in spite of the increased lean tissue mass, similar to the changes observed in animals in response to flight.

Evaluation of gene expression in pigs selected for enhanced reproduction using differential display PCR and human microarrays: I. Ovarian follicles1,2

Differential display PCR (ddPCR) and complementary DNA microarray analyses were used to evaluate gene expression differences in porcine ovarian follicles between a line of pigs selected for an index of ovulation rate and embryo survival (Line I) and its randomly selected control line (Line C). Follicles (4.0 to 7.0 mm) were dissected from ovaries of multiparous sows (n = 27) at either 2 or 4 d following PGF2alpha analog injection on d 12 to 14 of the estrous cycle. Using ddPCR, differentially expressed bands (n = 282) were excised from gels and 107 were sequenced, yielding 84 unique porcine follicle expressed sequence tags. Northern hybridization confirmed differential expression (between lines, days, or follicle sizes) for messenger RNA representing the calpain I light subunit, cytochrome C oxidase subunit III, cytochrome P450 aromatase, and cytochrome P450 side chain cleavage genes. For microarray analysis, two mRNA pools representing follicles (d 2; 4.50 to 4.75 mm) from Line I and Line C sows were hybridized to the Incyte UniGEM V1.0 human chip (approximately 7,000 gene probes). A second analysis was performed using mRNA from follicles (d 2; 4.50 to 5.00 mm) hybridized to the Incyte UniGEM V2.0 human chip (approximately 9,100 gene probes). A total of 33 and 21 genes were identified with significant expression differences using UniGEM V1.0 and V2.0, respectively (twofold or greater relative expression following adjustment for expression of control probes). However, there was little overlap between results of the two hybridizations. Expression differences between lines for two genes, follistatin and nuclear receptor subfamily 4, group A, member 1, were confirmed using Northern hybridization. These results demonstrate changes in follicular gene expression as the result of long-term selection for enhanced reproduction. These correlated responses may directly represent allelic variation utilized by selection (e.g., quantitative trait loci), or more likely, transcriptional changes in other genes that interact with reproductive QTL. This work represents one of the first applications of gene expression analysis to evaluate long-term selection response in livestock populations.

Complementary DNA macroarray analyses of differential gene expression in porcine fetal and postnatal muscle

To study differential gene expression in porcine skeletal muscle, a porcine complementary DNA (cDNA) macroarray was produced that contained 327 expressed sequence tags (EST) derived from whole embryo and adult skeletal muscle, and differential display PCR products from fetal and postnatal muscle. Total RNA from four muscle samples, 75- and 105-d fetal hind limb muscles, and 1- and 7-wk postnatal semitendinosus muscle was used to make radiolabeled targets for duplicate hybridization to the macroarray membranes in an initial screen for expression. All EST that gave clear signals (n = 238) were then re-arrayed, and hybridization was conducted with additional biological replication of samples in the 75-d and 1-wk ages. Signal intensity for each gene was normalized to signal intensity measured at control spots on each membrane, which consisted of total cDNA from liver, lung, spleen, and skeletal muscle. Both normalized ratio levels and a mixed linear model analyses were used to identify genes differentially expressed among the muscle samples. Results showed 28 genes had differences in expression level greater than twofold between the 75-d fetal and 1-wk muscle RNA samples. All 28 genes were also identified as genes with significantly different (P < 0.01) expression using a mixed linear model analysis. Nineteen of these 28 genes had significant matches (basic local alignment search tool [BLAST] score > 100; P < 0.01) to known genes, two matched genes encoding human hypothetical proteins, and seven had no significant matches to Genbank nonredundant and dbEST (database of expressed sequence tags) entries. These results were confirmed for representative genes with RNA blot analysis of seven developmental time points, including RNA from the same muscle samples tested previously in the macroarray. The RNA blot results confirmed the macroarray results for all selected genes, demonstrating that the macroarray technique used in this study is accurate and reproducible. An unknown muscle clone (M218) with a slightly less than twofold increase in expression from the 75-d to the 1-wk age (1 wk/75 d = 1.94; P = 0.0114) was also shown to differ between these two ages using RNA blot analysis, demonstrating the methods used to identify differentially expressed genes may be conservative. The association between expression patterns of vimentin and desmin was also investigated. Results indicate the switch in intermediate filament protein from vimentin to desmin occurs primarily at the level of transcription and/or RNA processing.

Evaluation of hypothalamic gene expression in mice divergently selected for heat loss

Mouse lines divergently selected for heat loss were evaluated for correlated responses in the hypothalamic transcriptome. High (MH) heat loss mice have approximately 50% greater heat loss, approximately 35% less body fat, approximately 20% greater feed intake, approximately 100% greater locomotor activity levels, and higher core body temperature compared with low (ML) heat loss mice. We evaluated hypothalamic expression between inbred lines derived from MH and ML lines (IH and IL, respectively) using cDNA microarrays and selected genes previously isolated in a large differential-display PCR experiment. Northern analysis was used to confirm differences, revealing higher hypothalamic mRNA expression of oxytocin (Oxt) and tissue inhibitor of metalloproteinase 2 (Timp-2) in the IH line. Real-time PCR assays were developed for Oxt, Timp-2, and ribosomal protein L3 (Rpl3, previously found to be upregulated in IL) and confirmed differential expression of these genes with potential physiological relevance in energy balance. These results provide information on correlated responses in the transcriptome of mice selected for high and low energy expenditure and reveal new information regarding genetic regulation of energy balance.

Differential expression of NAT1 translational repressor during development of bovine intramuscular adipocytes

This study was undertaken to test for differential gene expression in intramuscular adipocytes during fat deposition of feedlot steers. Angus x Hereford steers (n = 50) were fed a high-energy concentrate ration ad libitum for 20 (n = 5), 86 (n = 15), 121 (n = 15), and 146 days (n = 15) to obtain various degrees of intramuscular adipocyte development. Carcass traits were significantly different (P < 0.05) between the groups. Intramuscular adipose tissue was excised from the longissimus dorsi and snap frozen in liquid nitrogen. Pooled samples of total RNA representing each group were analyzed by differential-display polymerase chain reaction using 200 primer combinations comprising 20 arbitrary (5') and 10 anchor (3') oligonucleotides. Bands (n = 70) representing putative differences among treatment groups were excised, sequenced, and subjected to BLAST homology search. From these, 40 contained significant homology to known genes. One was of particular interest, the translational repressor NAT1 (novel APOBEC-1 target-1). NAT1 mRNA was quantified in individual animals to confirm differential expression among treatment groups. Results indicate that NAT1 message is more abundant (P < 0.05) in intramuscular adipocytes of younger/leaner animals.

Influence of ethanol on thermoregulation: mapping quantitative trait loci

The genetic basis for the effects of ethanol on thermoregulation was investigated by utilizing recombinant inbred mouse strains from C57BL/6J and DBA/2J progenitor strains. Changes in core body temperature (T(c)) and the degree of fluctuation of T(c) were monitored in male mice following the administration of ethanol in an environment with cyclic changes in ambient temperature (T(a)). Changes in T(c) were utilized to assess ethanol-induced effects on regulated T(c), whereas fluctuations in T(c) were utilized to assess thermoregulatory disruption. Ethanol was administered intraperitoneally at 1.5, 2.5, and 3.5 g/kg for all strains. Change in T(c) and increase in tail temperature were also evaluated at 2.5 g/kg ethanol in a constant T(a) of 26 degrees C. Associations between the measured physiological responses and previously mapped genetic markers were used to identify quantitative trait loci (QTLs). This established probable chromosome locations for a number of genes for the responses. To our knowledge, this is the first report of QTLs that underlie changes in regulation as well as the disruption of a physiological regulatory system.

Comparative mapping of RPL3, a gene overexpressed in multiple obesity models

The ribosomal protein 3 gene is differentially expressed in hypothalamus and brown adipose tissue between mouse lines divergently selected for heat loss, and in skeletal muscle of the ob/ob mouse model. Unfortunately, multiple Rpl3-processed pseudogenes have hampered mapping of the functional gene copy in mammalian species. Using PCR amplification with intronic primer binding, we have mapped Rpl3 to MMU15, and have also localized RPL3 to BTA5 in cattle. Comparative mapping implicates a previously mapped copy of RPL3 on HSA22 as the functional copy of human RPL3, while predictive mapping places the porcine homologue on SSC5.