PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of gdataGuide for AuthorsAboutExplore this JournalGenomics Data
 
Genom Data. 2017 June; 12: 24–27.
Published online 2017 February 24. doi:  10.1016/j.gdata.2017.02.009
PMCID: PMC5331152

Hypothalamic transcriptome analysis of congenic-derived F2 mice (chromosome 17:3–45 Mb) exhibiting preferential carbohydrate (versus fat) intake

Table thumbnail

2. Experimental design, materials and methods

2.1. Experimental design

Along with the general drive to eat, macronutrient-specific appetites are likely encoded by unique molecular changes in the hypothalamus, which occur largely as the result of unknown interoceptive signaling [1]. We conducted a transcriptome analysis of total RNA extracted from the whole hypothalamus of B6·CAST-17.1 congenic F2 mice that were either homozygous cast/cast or b6/b6 in the chromosome 17 congenic interval (Chr 17:3.19–45.73 Mb), on an otherwise B6 genome (see strain description below). A complete description of samples used and files that are available on GEO is found in Table 1.

Table 1
Sample information.

2.2. Mice

Our development of the B6·CAST-17.1 congenic strain has been described previously [2], [3], [4]. Here, the 42.5 Mb cast/cast congenic donor segment was delimited by SNP markers rs49640908 (proximal) and rs48762654 (distal). To fine-map the Mnic1 QTL using genetic recombination, a large congenic-by-recipient F2 population was generated [5]. For transcriptome analysis, only non-recombinant congenic F2 littermates were selected, i.e., those possessing a homozygous genotype of either cast/cast or b6/b6 throughout the 42.5 Mb chromosome 17 congenic region, and a b6/b6 genotype across the rest of the genome. All animal procedures were approved by the Pennington Biomedical Research Center Institutional Animal Care and Use Committee.

2.3. Animal genotyping

Genomic DNA was isolated from spleen and a custom SNP panel (Illumina)was designed to genotype the animals; for details, see [5]. Genotypes were coded for the chromosome 17 congenic segment based on 301 SNP markers.

2.4. Phenotyping

Food intake was determined in singly-housed adult male mice using a two-choice macronutrient diet protocol. For 10 days mice were given a choice between two diets: a fat + protein versus a carbohydrate + protein mixture, each containing protein (casein) (22% of energy) with the balance of calories contributed by either fat (vegetable shortening) or carbohydrate (corn starch-sucrose) (78% of energy). Both diets contained minerals, vitamins and cellulose. Diet intake including all spillage were measured daily to 0.1 g. The experimental diet composition and complete details of the phenotyping procedures have been previously described [5], [6].

2.5. Sample selection for gene expression analyses

This experiment was carried out using B6·CAST-17.1 congenic-derived F2 mice that possessed a non-recombined, chromosome 17 donor interval (3.19–45.73 Mb). Based on high-density SNP marker analysis, we determined the absence of recombination in this congenic interval, resulting in genotypes of either cast/cast (n = 44) or b6/b6 (n = 47) across the 42.5 Mb region, and that of b6/b6 throughout the remaining genome. Animals were first subjected to a 10 d macronutrient selection test (see Section 2.4), followed by an extended wash-out period on rodent chow. The macronutrient-rich diets were then reinitiated for 48 h prior to euthanasia and tissue harvest. The 48 h time point for tissue collection was chosen on the basis of temporal variation in genetic linkage, i.e., the absence of genetic linkage for the carbohydrate-rich diet on day 1, and the presence of linkage on days 2 and following [4]. Metabolic signals arising from food ingestion may act to influence the animals' choice beginning on day 2 of diet exposure. Remarkably, the proportion of carbohydrate + protein versus fat + protein diet selected (kcal%) between the 10 d and 2 d macronutrient diet selection tests was highly correlated (r = 0.86, P < 0.01), emphasizing the stability of the phenotype (preferential carbohydrate intake). The GEO repository contains the data obtained from individuals that displayed the most extreme values for self-selected intake of the carbohydrate + protein diet mixture (versus fat + protein diet) in the 10 d test (10 d sum of kcal): cast/cast congenic F2 mice in the upper quartile (n = 12) and b6/b6 congenic F2 s in the lowest quartile (n = 12) (Fig. 1). Due to a smaller RNA integrity number (RIN), two samples from the 25th quartile for each genotype were replaced with samples associated with the next highest (cast/cast) or lowest (b6/b6) phenotypic value. Total RNA was isolated from whole hypothalamus, collected from individual animals. RNA sequencing analysis of gene expression was performed on these 24 samples and genes that were differentially expressed between the two genotypes were identified.

Fig. 1
Scatter plot.

2.6. RNA extraction and quality control

Total RNA from whole hypothalamus was extracted and purified using the AllPrep RNA/Protein kit (Qiagen, Valencia, CA) and re-suspended in Nanopure water. RNA quality was assessed via nanodrop (Thermo Scientific, Waltham, MA) and Agilent 2100 Bioanalyzer (Santa Clara, CA); all samples showed an RIN of ≥ 7.0 with the exception of two cast/cast samples (6.6, 6.8) and one b6/b6 sample (6.6).

2.7. cDNA library preparation, SAGE-sequencing & transcriptome analysis

Transcriptome profiling was performed by 3′-expression tag sequencing (SAGE) on an Applied Biosystems (AB) SOLiD 5500XL next generation sequencer. Briefly, sequencing libraries containing 27-bp, 3′ tags for all transcripts within a sample were constructed from hypothalamus using the SOLiD SAGE kit from the manufacturer (Life Technologies), at the PBRC Genomics Core Facility. Each library was then labeled with a unique barcode sequence. Sequence mapping was performed using a modified version of the SOLiD SAGE Analysis Software v1.10 (Life Technologies) and defined analysis parameters. Sequence reads were aligned to mouse RefSeq transcripts (genome build GRCm38/mm10) as the reference. Tag hits, i.e., successfully aligned reads, were normalized or adjusted for coverage according to DESeq by estimating the size factor (median of the ratios of observed counts) for each sample library, and dividing the sample counts by the corresponding size factor [7]. The mapping statistics are summarized in Table 1; only uniquely mapped sequence reads were included in the expression count for each RefSeq gene.

2.8. Differential expression analysis

A principal components analysis (PCA) of gene expression data identified directions or principal components comprising the largest variation in the data [8]. As shown in the two-dimensional PCA biplot (Fig. 3), this analysis revealed two cluster-like patterns of overall gene expression levels characteristic of the congenic cast/cast and b6/b6 F2 samples. These patterns demonstrate a large influence by genotype on the gene expression profile of hypothalamic cells.

Fig. 3
Scatter plot of PCA for 24 transcript profiles.

Differential expression analysis of count data was carried out by means of DESeq, an R/Bioconductor package [7]. Genes having both ≥ 1.5 fold change in expression and P < 0.01, were considered to be significantly differentially expressed in the congenic cast/cast F2 compared to congenic b6/b6 F2 samples. Two-way hierarchical clustering, via Ward's minimum variance criterion method, was applied to normalized and standardized expression data using the tools in JMP Genomics, Version 10 (SAS Institute Inc., Cary, NC) [5]. The SAGE-seq data were transformed into a scatter plot of the log2 ratio (fold change) versus the mean (Fig. 2). A symmetrical distribution for gene expression was obtained, validating the data quality.

Fig. 2
MA plot.

Consistent with the genetic model, the 42.5 Mb Chr 17 cast/cast congenic donor interval was enriched for differentially expressed genes. Using filtering criteria of FC ≥ 1.5 and P < 0.01, we identified 86 differentially expressed genes within this interval (Chr 17:3.19–45.73 Mb), while none outside this interval met these criteria. Of the 86 differentially expressed genes, 21 exhibited increased expression and 65 were decreased in congenic cast/cast F2s compared to congenic b6/b6 F2. When we employed a less stringent FC criterion of ≥ 1.2, our analysis revealed 1–3 differentially expressed genes on each of eleven other autosomes, for a total of 18 DE genes outside the boundaries of the cast/cast congenic donor region, defined as Chr 17:3.19–45.73 Mb by SNP genotyping [5]. Differentially expressed genes located outside this 42.5 Mb congenic region are postulated to be under the control of trans-acting genetic factors.

3. Discussion

Here we describe our methods of hypothalamic transcriptome analysis to generate a list of differentially expressed genes in carbohydrate-preferring B6·CAST-17.1 cast/cast congenic F2 mice compared to fat-preferring B6-CAST-17.1 b6/b6 congenic F2, in which both genotypes possess an otherwise b6/b6 genome. The Chr 17 CAST congenic segment encompasses the fine-mapped Mnic QTL (Chr 17:26.08–45.12 Mb, 95% confidence interval) that is associated with the preferential carbohydrate intake of cast/cast animals in a macronutrient diet choice paradigm [5]. The hypothalamic transcriptome represented by this fine-mapped 19.0 Mb region contains 55 differentially expressed genes, of which nearly 50% were classified as metabolism genes by gene ontology analysis [5]. Seven of these genes are associated with lipids and/or fatty acids and showed decreased expression in fat-preferring b6/b6 congenic F2 (Decr2, Ppard, Agpat1, Tnxb, Neu1, Pla2g7, Cyp39a1). By contrast, two genes with increased expression in carbohydrate-preferring cast/cast congenic F2 mice have relevance to carbohydrate metabolism (Glo1, Neu1). Notably, three of these same genes (Agpat1, Pla2g7, Glo1) were found to exhibit significantly differential expression, in the same direction and with similar magnitude, in the hypothalamus of other congenic lines containing a CAST Chr 17 segment on a B6 background [2], [4].

Funding

This research was supported by the National Institutes of Health (NIH) (DK-53113 to BKR). The work utilized the Genomics Core Facilities (GCF) which is partially supported by Center of Biomedical Research Excellence (COBRE) NIH 8 P20 GM103528 and Nutrition Obesity Research Center (NORC) NIH P30DK072476. We thank Dr. Michael Salbaum, Director of the GCF at Pennington Biomedical for helpful advice on the SAGE experiment.

References

1. Berthoud H.R., Munzberg H., Richards B.K., Morrison C.D. Neural and metabolic regulation of macronutrient intake and selection. Proc. Nutr. Soc. 2012;71:390–400. [PubMed]
2. Kumar K.G., DiCarlo L.M., Volaufova J., Zuberi A.R., Richards B.K. Increased physical activity cosegregates with higher intake of carbohydrate and total calories in a subcongenic mouse strain. Mamm. Genome. 2010;21:52–63. [PubMed]
3. Kumar K.G., Poole A.C., York B., Volaufova J., Zuberi A., Richards B.K. Quantitative trait loci for carbohydrate and total energy intake on mouse chromosome 17: congenic strain confirmation and candidate gene analyses (Glo1, Glp1r) Am. J. Phys. Regul. Integr. Comp. Phys. 2007;292:R207–R216. [PubMed]
4. Kumar K.G., Smith Richards B.K. Transcriptional profiling of chromosome 17 quantitative trait loci for carbohydrate and total calorie intake in a mouse congenic strain reveals candidate genes and pathways. J Nutrigenet Nutrigenomics. 2008;1:155–171. [PubMed]
5. Gularte-Merida R., DiCarlo L.M., Robertson G., Simon J., Johnson W.D., Kappen C., Medrano J.F., Richards B.K. High-resolution mapping of a genetic locus regulating preferential carbohydrate intake, total kilocalories, and food volume on mouse chromosome 17. PLoS One. 2014;9 [PMC free article] [PubMed]
6. Smith Richards B.K., Belton B.N., Poole A.C., Mancuso J.J., Churchill G.A., Li R., Volaufova J., Zuberi A., York B. QTL analysis of self-selected macronutrient diet intake: fat, carbohydrate, and total kilocalories. Physiol. Genomics. 2002;11:205–217. [PubMed]
7. Anders S., Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106. [PubMed]
8. Ringner M. What is principal component analysis? Nat. Biotechnol. 2008;26:303–304. [PubMed]

Articles from Genomics Data are provided here courtesy of Elsevier