Stochastic variation of transcript abundance in C57BL/6J mice

Interindividual differences of dietary fat-inducible Mest in white adipose tissue of C57BL/6J mice are not heritable

OBJECTIVE

Differences in white adipose tissue (WAT) expression of mesoderm-specific transcript (Mest) in C57BL6/J mice fed a high-fat diet (HFD) are concomitant with and predictive for the development of obesity. However, the basis for differences in WAT Mest among mice is unknown. This study investigated whether HFD-inducible WAT Mest, as well as susceptibility to obesity, is transmissible from parents to offspring.

METHODS

WAT biopsies of mice fed an HFD for 2 weeks identified parents with low and high WAT Mest for breeding. Obesity phenotypes, WAT Mest, hepatic gene expression, and serum metabolites were determined in offspring fed an HFD for 2 weeks.

RESULTS

Offspring showed no heritability of obesity or WAT Mest phenotypes from parents but did show hepatic and serum metabolite changes consistent with their WAT Mest. Importantly, retired male breeders showed WAT Mest expression congruent with initial WAT biopsies even though HFD exposure occurred early in life.

CONCLUSIONS

Disparity of HFD-induced Mest in mice is not heritable but, rather, is reestablished during each generation and remains fixed from an early age to adulthood. Short-term HFD feeding reveals variation of WAT Mest expression within isogenic mice that is positively associated with the development of obesity.

Comparative urine metabolomics of mice treated with non-toxic and toxic oral doses of (-)-epigallocatechin-3-gallate

The green tea polyphenol, (-)-epigallocatechin-3-gallate (EGCG), has been studied for its potential positive health effects, but human and animal model studies have reported potential toxicity at high oral bolus doses. This study used liquid chromatography-mass spectrometry-based metabolomics to compare the urinary EGCG metabolite profile after administration of a single non-toxic (100 mg kg-1) or toxic (750 mg kg-1) oral bolus dose to male C57BL6/J mice to better understand how EGCG metabolism varies with dose. EGCG metabolites, including methyl, glucuronide, sulfate, and glucoside conjugates, were tentatively identified based on their mass to charge (m/z) ratio and fragment ion patterns. Partial least squares discriminant analysis (PLS-DA) results showed clear separation of the urine metabolite profiles between treatment groups. The most differentiating metabolites in the negative and positive ion modes were provisionally identified as di-glucuronidated EGCG quinone and di-glucuronidated EGCG, respectively. The presence of EGCG oxidation products at toxic dose is consistent with studies showing that EGCG toxicity is associated with oxidative stress. Relative amounts of methylated metabolites increased with dose to a lesser extent than glucuronide and sulfate metabolites, indicating that methylation is more prominent at low doses, whereas glucuronidation and sulfation may be more important at higher doses. One limitation of the current work is that the lack of commercially-available EGCG metabolite standards prevented absolute metabolite quantification and identification. Despite this limitation, these findings provide a basis for better understanding the dose-dependent changes in EGCG metabolism and advance studies on how these differences may contribute to the toxicity of high doses of EGCG.

Circadian clocks of the kidney: function, mechanism, and regulation

An intrinsic cellular circadian clock is located in nearly every cell of the body. The peripheral circadian clocks within the cells of the kidney contribute to the regulation of a variety of renal processes. In this review, we summarize what is currently known regarding the function, mechanism, and regulation of kidney clocks. Additionally, the effect of extrarenal physiological processes, such as endocrine and neuronal signals, on kidney function is also reviewed. Circadian rhythms in renal function are an integral part of kidney physiology, underscoring the importance of considering time of day as a key biological variable. The field of circadian renal physiology is of tremendous relevance, but with limited physiological and mechanistic information on the kidney clocks this is an area in need of extensive investigation.

MicroRNAs in Learning and Memory and Their Impact on Alzheimer’s Disease

Learning and memory formation rely on the precise spatiotemporal regulation of gene expression, such as microRNA (miRNA)-associated silencing, to fine-tune gene expression for the induction and maintenance of synaptic plasticity. Much progress has been made in presenting direct evidence of miRNA regulation in learning and memory. Here, we summarize studies that have manipulated miRNA expression using various approaches in rodents, with changes in cognitive performance. Some of these are involved in well-known mechanisms, such as the CREB-dependent signaling pathway, and some of their roles are in fear- and stress-related disorders, particularly cognitive impairment. We also summarize extensive studies on miRNAs correlated with pathogenic tau and amyloid-β that drive the processes of Alzheimer’s disease (AD). Although altered miRNA profiles in human patients with AD and in mouse models have been well studied, little is known about their clinical applications and therapeutics. Studies on miRNAs as biomarkers still show inconsistencies, and more challenges need to be confronted in standardizing blood-based biomarkers for use in AD.

Activation of a hippocampal CREB-pCREB-miRNA-MEF2 axis modulates individual variation of spatial learning and memory capability

Phenotypic variation is a fundamental prerequisite for cell and organism evolution by natural selection. Whereas the role of stochastic gene expression in phenotypic diversity of genetically identical cells is well studied, not much is known regarding the relationship between stochastic gene expression and individual behavioral variation in animals. We demonstrate that a specific miRNA (miR-466f-3p) is upregulated in the hippocampus of a portion of individual inbred mice upon a Morris water maze task. Significantly, miR-466f-3p positively regulates the neuron morphology, function and spatial learning, and memory capability of mice. Mechanistically, miR-466f-3p represses translation of MEF2A, a negative regulator of learning/memory. Finally, we show that varied upregulation of hippocampal miR-466f-3p results from randomized phosphorylation of hippocampal cyclic AMP (cAMP)-response element binding (CREB) in individuals. This finding of modulation of spatial learning and memory via a randomized hippocampal signaling axis upon neuronal stimulation represents a demonstration of how variation in tissue gene expression lead to varied animal behavior.

Human gene expression variability and its dependence on methylation and aging

BACKGROUND

Phenotypic variability of human populations is partly the result of gene polymorphism and differential gene expression. As such, understanding the molecular basis for diversity requires identifying genes with both high and low population expression variance and identifying the mechanisms underlying their expression control. Key issues remain unanswered with respect to expression variability in human populations. The role of gene methylation as well as the contribution that age, sex and tissue-specific factors have on expression variability are not well understood.

RESULTS

Here we used a novel method that accounts for sampling error to classify human genes based on their expression variability in normal human breast and brain tissues. We find that high expression variability is almost exclusively unimodal, indicating that variance is not the result of segregation into distinct expression states. Genes with high expression variability differ markedly between tissues and we find that genes with high population expression variability are likely to have age-, but not sex-dependent expression. Lastly, we find that methylation likely has a key role in controlling expression variability insofar as genes with low expression variability are likely to be non-methylated.

CONCLUSIONS

We conclude that gene expression variability in the human population is likely to be important in tissue development and identity, methylation, and in natural biological aging. The expression variability of a gene is an important functional characteristic of the gene itself and the classification of a gene as one with Hyper-Variability or Hypo-Variability in a human population or in a specific tissue should be useful in the identification of important genes that functionally regulate development or disease.

Predictability of human differential gene expression

The identification of genes that are differentially expressed provides a molecular foothold onto biological questions of interest. Whether some genes are more likely to be differentially expressed than others, and to what degree, has never been assessed on a global scale. Here, we reanalyze more than 600 studies and find that knowledge of a gene’s prior probability of differential expression (DE) allows for accurate prediction of DE hit lists, regardless of the biological question. This result suggests redundancy in transcriptomics experiments that both informs gene set interpretation and highlights room for growth within the field.

Differential expression (DE) is commonly used to explore molecular mechanisms of biological conditions. While many studies report significant results between their groups of interest, the degree to which results are specific to the question at hand is not generally assessed, potentially leading to inaccurate interpretation. This could be particularly problematic for metaanalysis where replicability across datasets is taken as strong evidence for the existence of a specific, biologically relevant signal, but which instead may arise from recurrence of generic processes. To address this, we developed an approach to predict DE based on an analysis of over 600 studies. A predictor based on empirical prior probability of DE performs very well at this task (mean area under the receiver operating characteristic curve, ∼0.8), indicating that a large fraction of DE hit lists are nonspecific. In contrast, predictors based on attributes such as gene function, mutation rates, or network features perform poorly. Genes associated with sex, the extracellular matrix, the immune system, and stress responses are prominent within the “DE prior.” In a series of control studies, we show that these patterns reflect shared biology rather than technical artifacts or ascertainment biases. Finally, we demonstrate the application of the DE prior to data interpretation in three use cases: (i) breast cancer subtyping, (ii) single-cell genomics of pancreatic islet cells, and (iii) metaanalysis of lung adenocarcinoma and renal transplant rejection transcriptomics. In all cases, we find hallmarks of generic DE, highlighting the need for nuanced interpretation of gene phenotypic associations.

The hepatic innate immune response is lobe-specific in a murine model endotoxemia

The liver plays a central role in the innate immune response to endotoxemia. While previous studies have demonstrated lobe-specific transcriptional responses to various insults, whether this is true in response to endotoxemia is unknown. We sought to assess whether there were significant intra- and inter-lobe differences in the murine hepatic innate immune transcriptional response to endotoxemia. Adult male ICR mice were exposed to i.p. LPS (5 mg/kg, 30 min, 60 min, 5 h) and primary (Tnf, Cxcl1, Nfkbia, Tnfiap3) and secondary (Il6, Nos2) innate immune response gene expression was assessed in the left medial, right medial, left lateral, and right lateral lobes, and the papillary and caudate processes. The expression of all innate immune response genes increased following i.p. LPS challenge. When tested at the early time points (30 and 60 min), the left medial lobe and caudate process consistently demonstrated the highest induction of gene expression. Most inter-lobe differences were attenuated at later time points (5 h). To improve reproducibility of the study of endotoxemia induced by i.p. LPS challenge, inclusion of appropriate methodological details regarding collection of hepatic tissue should be included when reporting scientific results in published manuscripts.

Single subject transcriptome analysis to identify functionally signed gene set or pathway activity

Analysis of single-subject transcriptome response data is an unmet need of precision medicine, made challenging by the high dimension, dynamic nature and difficulty in extracting meaningful signals from biological or stochastic noise. We have proposed a method for single subject analysis that uses a mixture model for transcript fold-change clustering from isogenically paired samples, followed by integration of these distributions with Gene Ontology Biological Processes (GO-BP) to reduce dimension and identify functional attributes. We then extended these methods to develop functional signing metrics for gene set process regulation by incorporating biological repressor relationships encoded in GO-BP as negatively_regulates edges. Results revealed reproducible and biologically meaningful signals from analysis of a single subject's response, opening the door to future transcriptomic studies where subject and resource availability are currently limiting. We used inbred mouse strains fed different diets to provide isogenic biological replicates, permitting rigorous validation of our method. We compared significant genotype-specific GO-BP term results for overlap and rank order across three replicate pairs per genotype, and cross-methods to reference standards (limma+FET, SAM+FET, and GSEA). All single-subject analytics findings were robust and highly reproducible (median area under the ROC curve=0.96, n=24 genotypes × 3 replicates), providing confidence and validation of this approach for analyses in single subjects. R code is available online at http://www.lussiergroup.org/publications/PathwayActivity.

Design and Sampling Plan Optimization for RT-qPCR Experiments in Plants: A Case Study in Blueberry

The qPCR assay has become a routine technology in plant biotechnology and agricultural research. It is unlikely to be technically improved, but there are still challenges which center around minimizing the variability in results and transparency when reporting technical data in support of the conclusions of a study. There are a number of aspects of the pre- and post-assay workflow that contribute to variability of results. Here, through the study of the introduction of error in qPCR measurements at different stages of the workflow, we describe the most important causes of technical variability in a case study using blueberry. In this study, we found that the stage for which increasing the number of replicates would be the most beneficial depends on the tissue used. For example, we would recommend the use of more RT replicates when working with leaf tissue, while the use of more sampling (RNA extraction) replicates would be recommended when working with stems or fruits to obtain the most optimal results. The use of more qPCR replicates provides the least benefit as it is the most reproducible step. By knowing the distribution of error over an entire experiment and the costs at each step, we have developed a script to identify the optimal sampling plan within the limits of a given budget. These findings should help plant scientists improve the design of qPCR experiments and refine their laboratory practices in order to conduct qPCR assays in a more reliable-manner to produce more consistent and reproducible data.

Diurnal Variation Has Effect on Differential Gene Expression Analysis in the Hippocampus of the Pilocarpine-Induced Model of Mesial Temporal Lobe Epilepsy

The molecular mechanisms underlying epileptogenesis have been widely investigated by differential gene expression approach, especially RT-qPCR methodology. However, controversial findings highlight the occurrence of unpredictable sources of variance in the experimental designs. Here, we investigated if diurnal rhythms of transcript's levels may impact on differential gene expression analysis in hippocampus of rats with experimental epilepsy. For this, we have selected six core clock genes (Per1, Per3, Bmal1, Clock, Cry1 and Cry2), whose rhythmic expression pattern in hippocampus had been previously reported. Initially, we identified Tubb2a/Rplp1 and Tubb2a/Ppia as suitable normalizers for circadian studies in hippocampus of rats maintained to 12:12 hour light:dark (LD) cycle. Next, we confirmed the temporal profiling of Per1, Per3, Bmal1, Cry1 and Cry2 mRNA levels in the hippocampus of naive rats by both Acrophase and CircWave statistical tests for circadian analysis. Finally, we showed that temporal differences of sampling can change experimental results for Per1, Per3, Bmal1, Cry1 and Cry2, but not for Clock, which was consistently decreased in rats with epilepsy in all comparison to the naive group. In conclusion, our study demonstrates it is mandatory to consider diurnal oscillations, in order to avoid erroneous conclusions in gene expression analysis in hippocampus of rats with epilepsy. Investigators, therefore, should be aware that genes with circadian expression could be out of phase in different animals of experimental and control groups. Moreover, our results indicate that a sub-expression of Clock may be involved in epileptogenicity, although the functional significance of this remains to be investigated.

Gene expression in teratogenic exposures: a new approach to understanding individual risk

The phenomenon of partial or incomplete penetrance is common to many paradigms of exposure to teratogens, where only some of the exposed individuals exhibit developmental defects. We here argue that the most widely used experimental approaches in reproductive toxicology do not take partial penetrance into account, and are thus likely to miss differences between affected and unaffected individuals that contribute to susceptibility for teratogenesis. We propose that focus on the variation between exposed individuals could help to discover factors that may play a causative role for abnormal developmental processes that occur with incomplete penetrance.

Dietary oleanolic acid mediates circadian clock gene expression in liver independently of diet and animal model but requires apolipoprotein A1

Oleanolic acid is a triterpene widely distributed throughout the plant kingdom and present in virgin olive oil at a concentration of 57 mg/kg. To test the hypotheses that its long-term administration could modify hepatic gene expression in several animal models and that this could be influenced by the presence of APOA1-containing high-density lipoproteins (HDLs), diets including 0.01% oleanolic acid were provided to Apoe- and Apoa1-deficient mice and F344 rats. Hepatic transcriptome was analyzed in Apoe-deficient mice fed long-term semipurified Western diets differing in the oleanolic acid content. Gene expression changes, confirmed by reverse transcriptase quantitative polymerase chain reaction, were sought for their implication in hepatic steatosis. To establish the effect of oleanolic acid independently of diet and animal model, male rats were fed chow diet with or without oleanolic acid, and to test the influence of HDL, Apoa1-deficient mice consuming the latter diet were used. In Apoe-deficient mice, oleanolic acid intake increased hepatic area occupied by lipid droplets with no change in oxidative stress. Bmal1 and the other core component of the circadian clock, Clock, together with Elovl3, Tubb2a and Cldn1 expressions, were significantly increased, while Amy2a5, Usp2, Per3 and Thrsp were significantly decreased in mice receiving the compound. Bmal1 and Cldn1 expressions were positively associated with lipid droplets. Increased Clock and Bmal1 expressions were also observed in rats, but not in Apoa1-deficient mice. The core liver clock components Clock-Bmal1 are a target of oleanolic acid in two animal models independently of the diets provided, and this compound requires APOA1-HDL for its hepatic action.

The use of transcriptomics to unveil the role of nutrients in Mammalian liver

Liver is the organ primarily responding to diet, and it is crucial in determining plasma carbohydrate, protein, and lipid levels. In addition, it is mainly responsible for transformation of xenobiotics. For these reasons, it has been a target of transcriptomic analyses. In this review, we have covered the works dealing with the response of mammalian liver to different nutritional stimuli such as fasting/feeding, caloric restriction, dietary carbohydrate, cholesterol, fat, protein, bile acid, salt, vitamin, and oligoelement contents. Quality of fats or proteins has been equally addressed, and has the influence of minor dietary components. Other compounds, not purely nutritional as those represented by alcohol and food additives, have been included due to their relevance in processed food. The influence has been studied not only on mRNA but also on miRNA. The wide scope of the technology clearly reflects that any simple intervention has profound changes in many metabolic parameters and that there is a synergy in response when more compounds are included in the intervention. Standardized arrays to systematically test the same genes in all studies and analyzing data to establish patterns of response are required, particularly for RNA sequencing. Moreover, RNA is a valuable, easy-screening ally but always requires further confirmation.

Predicting the outcome of infectious diseases: variability among inbred mice as a new and powerful tool for biomarker discovery

Individuals respond differently to infectious diseases. Even among inbred mice that are presumed to be genetically identical, the response to a microbial pathogen is variable, which is generally thought to reflect experimental inconsistencies, technical errors, and stochastic processes. Here we describe the remarkable observation that the variability of Helicobacter pylori colonization density in the stomachs of experimentally infected C57BL/6J mice is tightly correlated with weight loss and viral load after a challenge with influenza virus, though H. pylori infection per se does not affect influenza and vice versa. Since these two infectious agents are found in different tissue compartments and are detected using unrelated methods, the correlation in microbial burden must represent a biological measure of disease susceptibility among genetically nearly identical individuals and not technical or stochastic factors. We hypothesize that inbred mice represent a powerful new tool for the identification of biomarkers to predict the outcome of infectious diseases.

Strain-specific hyperkyphosis and megaesophagus in Add1 null mice

The three adducin proteins (α, β, and γ) share extensive sequence, structural, and functional homology. Heterodimers of α- and β-adducin are vital components of the red cell membrane skeleton, which is required to maintain red cell elasticity and structural integrity. In addition to anemia, targeted deletion of the α-adducin gene (Add1) reveals unexpected, strain-dependent non-erythroid phenotypes. On an inbred 129 genetic background, Add1 null mice show abnormal inward curvature of the cervicothoracic spine with complete penetrance. More surprisingly, a subset of 129-Add1 null mice develop severe megaesophagus, while examination of peripheral nerves reveals a reduced number of axons in 129-Add1 null mice at four months of age. These unforeseen phenotypes, described here, reveal new functions for adducin and provide new models of mammalian disease.

Survival study of metallothionein-1 transgenic mice and respective controls (C57BL/6J): influence of a zinc-enriched environment

The role of metallothioneins (MTs) in aging is not completely understood. Several studies have shown evidence that these proteins could represent a defense system against oxidative damage, but survival studies on mice overexpressing MTs are poor. Here we describe a survival study performed on old MT-1-overexpressing mice (MT-TG) and their respective controls (C57BL/6J) fed a standard or zinc (Zn)-supplemented diet. MT-TG mice had significantly increased survival compared with control. Zn supplementation affects the survival curves of MT-TG and C57BL/6J mice differently. This study poses the basis for intervention based on gene therapy with MTs to enhance the health span of laboratory mice.

Simultaneous measurement of insulin sensitivity, insulin secretion, and the disposition index in conscious unhandled mice

Of the parameters that determine glucose disposal and progression to diabetes in humans: first-phase insulin secretion, glucose effectiveness (Sg), insulin sensitivity (Si), and the disposition index (DI), only Si can be reliably measured in conscious mice. To determine the importance of the other parameters in murine glucose homeostasis in lean and obese states, we developed the frequently sampled intravenous glucose tolerance test (FSIVGTT) for use in unhandled mice. We validated the conscious FSIVGTT against the euglycemic clamp for measuring Si in lean and obese mice. Insulin-resistant mice had increased first-phase insulin secretion, decreased Sg, and a reduced DI, qualitatively similar to humans. Intriguingly, although insulin secretion explained most of the variation in glucose disposal in lean mice, Sg and the DI more strongly predicted glucose disposal in obese mice. DI curves identified individual diet-induced obese (DIO) mice as having compensated or decompensated insulin secretion. Conscious FSIVGTT opens the door to apply mouse genetics to the determinants of in vivo insulin secretion, Sg, and DI, and further validates the mouse as a model of metabolic disease.

Sources of variance in baseline gene expression in the rodent liver

The use of gene expression profiling in both clinical and laboratory settings would be enhanced by better characterization of variation due to individual, environmental, and technical factors. Analysis of microarray data from untreated or vehicle-treated animals within the control arm of toxicogenomics studies has yielded useful information on baseline fluctuations in liver gene expression in the rodent. Here, studies which highlight contributions of different factors to gene expression variability in the rodent liver are discussed including a large meta-analysis of rat liver, which identified genes that vary in control animals in the absence of chemical treatment. Genes and their pathways that are the most and least variable were identified in a number of these studies. Life stage, fasting, sex, diet, circadian rhythm and liver lobe source can profoundly influence gene expression in the liver. Recognition of biological and technical factors that contribute to variability of background gene expression can help the investigator in the design of an experiment that maximizes sensitivity and reduces the influence of confounders that may lead to misinterpretation of genomic changes. The factors that contribute to variability in liver gene expression in rodents are likely analogous to those contributing to human interindividual variability in drug response and chemical toxicity. Identification of batteries of genes that are altered in a variety of background conditions could be used to predict responses to drugs and chemicals in appropriate models of the human liver.

Comparison of hypothalamic mRNA levels in mice euthanized by CO₂ inhalation and focused-beam microwave irradiation

Focused-beam microwave irradiation (FBMI) is a relatively new method for euthanasia of small mammals and is available to most researchers. Compared with CO₂ inhalation, this method of euthanasia has the advantage of preserving fast-degrading metabolites. But differences in brain RNA quantity and quality, gene expression and histology in mice euthanized by CO₂ inhalation versus FBMI have not been investigated. Here the authors report that a smaller quantity of RNA was isolated from brains of mice euthanized by FBMI compared with those of mice euthanized by CO₂ inhalation. They also found relative differences in the levels of the expression of some genes. These studies suggest that either method can be used for histological analysis or RNA isolation, but the authors caution against combining the techniques within a single study on gene expression.