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.     

Identification of high and low responders to lipopolysaccharide in normal subjects: an unbiased approach to identify modulators of innate immunity

LPS stimulates a vigorous inflammatory response from circulating leukocytes that varies greatly from individual to individual. The goal of this study was to use an unbiased approach to identify differences in gene expression that may account for the high degree of interindividual variability in inflammatory responses to LPS in the normal human population. We measured LPS-induced cytokine production ex vivo in whole blood from 102 healthy human subjects and identified individuals who consistently showed either very high or very low responses to LPS (denoted lps(high) and lps(low), respectively). Comparison of gene expression profiles between the lps(high) and lps(low) individuals revealed 80 genes that were differentially expressed in the presence of LPS and 21 genes that were differentially expressed in the absence of LPS (p < 0.005, ANOVA). Expression of a subset of these genes was confirmed using real-time RT-PCR. Functional relevance for one gene confirmed to be expressed at a higher level in lps(high), adipophilin, was inferred when reduction in adipophilin mRNA by small interfering RNA in the human monocyte-like cell line THP-1 resulted in a modest but significant reduction in LPS-induced MCP-1 mRNA expression. These data illustrate a novel approach to the identification of factors that determine interindividual variability in innate immune inflammatory responses and identify adipophilin as a novel potential regulator of LPS-induced MCP-1 production in human monocytes.     

Specific site methylation in the 5'-flanking region of CYP1A2 interindividual differences in human livers

Human cytochrome P4501A2 (CYP1A2) is involved in the metabolism of a large number of common drugs and is responsible for the metabolic activation of numerous promutagens and procarcinogens. Large interindividual differences exist in the expression of this enzyme. This variability has important implications for drug efficacy and cancer susceptibility. In this sudy, the methylation status of the CCGG site (bp -2759) located adjacent to an AP-1 site in the 5'-flanking region of the CYP1A2 gene was assessed in liver samples from a pool of 55 human donors. DNA methylation is an important epigenetic mechanism controlling gene expression and may be one of the molecular mechanisms underlying the interindividual variation. Analysis was conducted using Hpa II digestion and PCR. Results showed that individual samples varied in the methylation status at this site. The site was found to be hypermethylated in approximately 60% of the samples. To compare methylation status with level of CYP1A2 expression, results were analyzed by median test. In 44% of the samples with expression levels above the median the CCGG site was hypermethylated. However, for those samples with levels below the median hypermethylation of the site was found in 73% of the samples. The difference was statistically significant (p<0.05). These findings indicate that CpG methylation may be involved in controlling the expression of CYP1A2, with hypermethylation reducing expression. Moreover, this approach can be useful in assessing the role of site-specific DNA methylation in interindividual variation of CYP1A2. Analysis of other CpG sites in potentially important regulatory elements of the CYP1A2 gene will continue.     

Hepatocytes--the choice to investigate drug metabolism and toxicity in man: in vitro variability as a reflection of in vivo

The pharmaceutical industry is committed to marketing safer drugs with fewer side effects, predictable pharmacokinetic properties and quantifiable drug-drug interactions. Drug metabolism is a major determinant of drug clearance and interindividual pharmacokinetic differences, and an indirect determinant of the clinical efficacy and toxicity of drugs. Progressive advances in the knowledge of metabolic routes and enzymes responsible for drug biotransformation have contributed to understanding the great metabolic variations existing in human beings. Phenotypic as well genotypic differences in the expression of the enzymes involved in drug metabolism are the main causes of this variability. However, only a minor part of phenotypic variability in man is attributable to gene polymorphisms, thus making the definition of a normal liver complex. At present, the use of human in vitro hepatic models at early preclinical stages means that the process of selecting drug candidates is becoming much more rational. Cultured human hepatocytes are considered to be the closest model to human liver. However, the fact that hepatocytes are located in a microenvironment that differs from that of the cell in the liver raises the question: to what extent does drug metabolism variability observed in vitro actually reflect that of the liver in vivo? By comparing the metabolism of a model compound both in vitro and in vivo in the same individual, a good correlation between the in vitro and in vivo relative abundance of oxidized metabolites and the hydrolysis of the compound was observed. Thus, it is reasonable to consider that the variability observed in human hepatocytes reflects the existing phenotypic heterogeneity of the P450 expression in human liver.     

Significance analysis and statistical dissection of variably methylated regions

It has recently been proposed that variation in DNA methylation at specific genomic locations may play an important role in the development of complex diseases such as cancer. Here, we develop 1- and 2-group multiple testing procedures for identifying and quantifying regions of DNA methylation variability. Our method is the first genome-wide statistical significance calculation for increased or differential variability, as opposed to the traditional approach of testing for mean changes. We apply these procedures to genome-wide methylation data obtained from biological and technical replicates and provide the first statistical proof that variably methylated regions exist and are due to interindividual variation. We also show that differentially variable regions in colon tumor and normal tissue show enrichment of genes regulating gene expression, cell morphogenesis, and development, supporting a biological role for DNA methylation variability in cancer.     

Gene Expression Variability within and between Human Populations and Implications toward Disease Susceptibility

Variations in gene expression level might lead to phenotypic diversity across individuals or populations. Although many human genes are found to have differential mRNA levels between populations, the extent of gene expression that could vary within and between populations largely remains elusive. To investigate the dynamic range of gene expression, we analyzed the expression variability of ∼18, 000 human genes across individuals within HapMap populations. Although ∼20% of human genes show differentiated mRNA levels between populations, our results show that expression variability of most human genes in one population is not significantly deviant from another population, except for a small fraction that do show substantially higher expression variability in a particular population. By associating expression variability with sequence polymorphism, intriguingly, we found SNPs in the untranslated regions (5′ and 3′UTRs) of these variable genes show consistently elevated population heterozygosity. We performed differential expression analysis on a genome-wide scale, and found substantially reduced expression variability for a large number of genes, prohibiting them from being differentially expressed between populations. Functional analysis revealed that genes with the greatest within-population expression variability are significantly enriched for chemokine signaling in HIV-1 infection, and for HIV-interacting proteins that control viral entry, replication, and propagation. This observation combined with the finding that known human HIV host factors show substantially elevated expression variability, collectively suggest that gene expression variability might explain differential HIV susceptibility across individuals.     

Differential metabolic and multi-tissue transcriptomic responses to fructose consumption among genetically diverse mice

Understanding how individuals react differently to the same treatment is a major concern in precision medicine. Metabolic challenges such as the one posed by high fructose intake are important determinants of disease mechanisms. We embarked on studies to determine how fructose affects differential metabolic dysfunctions across genetically dissimilar mice, namely, C57BL/6 J (B6), DBA/2 J (DBA) and FVB/NJ (FVB), by integrating physiological and gene regulatory mechanisms. We report that fructose has strain-specific effects, involving tissue-specific gene regulatory cascades in hypothalamus, liver, and white adipose tissues. DBA mice showed the largest numbers of genes associated with adiposity, congruent with their highest susceptibility to adiposity gain and glucose intolerance across the three tissues. In contrast, B6 and FVB mainly exhibited cholesterol phenotypes, accompanying the largest number of adipose genes correlating with total cholesterol in B6, and liver genes correlating with LDL in FVB mice. Tissue-specific network modeling predicted strain-and tissue-specific regulators such as Fgf21 (DBA) and Lss (B6), which were subsequently validated in primary hepatocytes. Strain-specific fructose-responsive genes revealed susceptibility for human diseases such that genes in liver and adipose tissue in DBA showed strong enrichment for human type 2 diabetes and obesity traits. Liver and adipose genes in FVB were mostly related to lipid traits, and liver and adipose genes in B6 showed relevance to most cardiometabolic traits tested. Our results show that fructose induces gene regulatory pathways that are tissue specific and dependent on the genetic make-up, which may underlie interindividual variability in cardiometabolic responses to high fructose consumption.     

DNA methylation-dependent regulation of TrkA, TrkB, and TrkC genes in human hepatocellular carcinoma

The tropomyosin-related kinase (Trk) family of neurotrophin receptors, TrkA, TrkB and TrkC, has been implicated in the growth and survival of human cancers. Here we report that Trks are frequently overexpressed in hepatocellular carcinoma (HCC) from patients and human liver cancer cell lines. To unravel the underlying molecular mechanism(s) for this phenomenon, DNA methylation patterns of CpG islands in TrkA, TrkB, and TrkC genes were examined in normal and cancer cell lines derived from liver. A good correlation was observed between promoter hypermethylation and lower expression of TrkA, TrkB, and TrkC genes, which was supported by the data that inhibiting DNA methylation with 5-azacytidine restored expression of those genes in normal liver cell lines. Furthermore, Trks promoted the proliferation of HepG2 and induced expression of the metastatic regulator, Twist. These results suggest that Trks may contribute to growth and metastasis of liver cancer.     

Human Placental-Specific Epipolymorphism and its Association with Adverse Pregnancy Outcomes

Interindividual variation in DNA-methylation level is widespread in the human genome, despite its critical role in regulating gene expression. The nature of this variation, including its tissue-specific nature, and the role it may play in human phenotypic variation and disease is still poorly characterized. The placenta plays a critical role in regulating fetal growth and development in ways that have lifelong effects on health. To identify genes with a high degree of interindividual DNA methylation variation in the human placenta, we surveyed the human genome using the Illumina GoldenGate Methylation Cancer panel targeting 1505 CpG sites of 807 genes. While many sites show a continuous pattern of methylation levels, WNT2, TUSC3 and EPHB4 were identified to have a polymorphic “on-or-off” pattern of DNA methylation variation at their promoter region which was confirmed by pyrosequencing. Methylation of these genes can be found in 7%–25% of over 100 placentas tested. The methylation state at the promoter of these genes is concordant with mRNA allelic expression. In three informative cases TUSC3 was observed to be methylated on the maternal allele, and it is thus possible this represents a polymorphically imprinted gene. Furthermore, TUSC3 promoter methylation showed evidence for association with preeclampsia. A biological significance of these methylation allelic polymorphisms (MAPs) to human placental diversity is further implied by their placental specificity and absence in mouse. An extended study of blood suggests that MAPs may also be found in other tissues, implicating their utility for tissue-specific association with complex disorders. The identification of such “epipolymorphism” in other tissues and their use in association studies, should improve our understanding of interindividual phenotypic variability and complex disease susceptibility.     

Functional Study of Carboxylesterase 1 Protein Isoforms

Carboxylesterase 1 (CES1) is a primary human hepatic hydrolase involved in hydrolytic biotransformation of numerous medications. Considerable interindividual variability in CES1 expression and activity has been consistently reported. Four isoforms of the CES1 protein are produced by alternative splicing (AS). In the current study, the activity and expression of each CES1 isoform are examined using transfected cell lines, and CES1 isoform composition and its impact on CES1 activity in human livers are determined. In transfected cells, isoforms 3 and 4 show mRNA and protein expressions comparable to isoforms 1 and 2, but have significantly impaired activity when hydrolyzing enalapril and clopidogrel. In individual human liver samples, isoforms 1 and 2 are the major forms, contributing 73–90% of the total CES1 protein expression. In addition, the protein expression ratios of isoforms 1 and 2 to isoforms 3 and 4 are positively associated with CES1 activity in the liver, suggesting that CES1 isoform composition is a factor contributing to the variability in hepatic CES1 function. Further investigations of the regulation of CES1 AS would improve the understanding of CES1 variability and help develop a strategy to optimize the pharmacotherapy of many CES1 substrate medications.     

Mixed-model reanalysis of primate data suggests tissue and species biases in oligonucleotide-based gene expression profiles

An emerging issue in evolutionary genetics is whether it is possible to use gene expression profiling to identify genes that are associated with morphological, physiological, or behavioral divergence between species and whether these genes have undergone positive selection. Some of these questions were addressed in a recent study (Enard et al. 2002) of the difference in gene expression among human, chimp, and orangutan, which suggested an accelerated rate of divergence in gene expression in the human brain relative to liver. Reanalysis of the Affymetrix data set using analysis of variance methods to quantify the contributions of individuals and species to variation in expression of 12,600 genes indicates that as much as one-quarter of the genome shows divergent expression between primate species at the 5% level. The magnitude of fold change ranges from 1.2-fold up to 8-fold. Similar conclusions apply to reanalysis of Enard et al. 2002 parallel murine data set. However, biases inherent to short oligonucleotide microarray technology may account for some of the tissue and species effects. At high significance levels, more differences were observed in the liver than in the brain in each of the pairwise species comparisons, so it is not clear that expression divergence is accelerated in the human brain. Further, there is an apparent bias toward upregulation of gene expression in the brain in both primates and mice, whereas genes are equally likely to be up- or downregulated in the liver when these species diverge. A small subset of genes that are candidates for adaptive divergence may be identified on the basis of a high ratio of interspecific to intraspecific divergence.     

A comparative analysis of the transcriptome and signal pathways in hepatic differentiation of human adipose mesenchymal stem cells

The specific features of the plasticity of adult stem cells are largely unknown. Recently, we demonstrated the hepatic differentiation of human adipose tissue-derived mesenchymal stem cells (AT-MSCs). To identify the genes responsible for hepatic differentiation, we examined the gene expression profiles of AT-MSC-derived hepatocytes (AT-MSC-Hepa) using several microarray methods. The resulting sets of differentially expressed genes (1639 clones) were comprehensively analyzed to identify the pathways expressed in AT-MSC-Hepa. Clustering analysis revealed a striking similarity of gene clusters between AT-MSC-Hepa and the whole liver, indicating that AT-MSC-Hepa were similar to liver with regard to gene expression. Further analysis showed that enriched categories of genes and signaling pathways such as complementary activation and the blood clotting cascade in the AT-MSC-Hepa were relevant to liver-specific functions. Notably, decreases in Twist and Snail expression indicated that mesenchymal-to-epithelial transition occurred in the differentiation of AT-MSCs into hepatocytes. Our data show a similarity between AT-MSC-Hepa and the liver, suggesting that AT-MSCs are modulated by their environmental conditions, and that AT-MSC-Hepa may be useful in basic studies of liver function as well as in the development of stem cell-based therapy.