Exogenous mitochondrial transfer increases energy expenditure and attenuates adiposity gains in mice with diet-induced obesity

Obesity is associated with chronic multi-system bioenergetic stress that may be improved by increasing the number of healthy mitochondria available across organ systems. However, treatments capable of increasing mitochondrial content are generally limited to endurance exercise training paradigms, which are not always sustainable long-term, let alone feasible for many patients with obesity. Recent studies have shown that local transfer of exogenous mitochondria from healthy donor tissues can improve bioenergetic outcomes and alleviate the effects of tissue injury in recipients with organ specific disease. Thus, the aim of this project was to determine the feasibility of systemic mitochondrial transfer for improving energy balance regulation in the setting of diet-induced obesity. We found that transplantation of mitochondria from lean mice into mice with diet-induced obesity attenuated adiposity gains by increasing energy expenditure and promoting the mobilization and oxidation of lipids. Additionally, mice that received exogenous mitochondria demonstrated improved glucose uptake, greater insulin responsiveness, and complete reversal of hepatic steatosis. These changes were, in part, driven by adaptations occurring in white adipose tissue. Together, these findings are proof-of-principle that mitochondrial transplantation is an effective therapeutic strategy for limiting the deleterious metabolic effects of diet-induced obesity in mice.

Pleiotropic influence of DNA methylation QTLs on physiological and ageing traits

DNA methylation is influenced by genetic and non-genetic factors. Here, we chart quantitative trait loci (QTLs) that modulate levels of methylation at highly conserved CpGs using liver methylome data from mouse strains belonging to the BXD family. A regulatory hotspot on chromosome 5 had the highest density of trans-acting methylation QTLs (trans-meQTLs) associated with multiple distant CpGs. We refer to this locus as meQTL.5a. Trans-modulated CpGs showed age-dependent changes and were enriched in developmental genes, including several members of the MODY pathway (maturity onset diabetes of the young). The joint modulation by genotype and ageing resulted in a more 'aged methylome' for BXD strains that inherited the DBA/2J parental allele at meQTL.5a. Further, several gene expression traits, body weight, and lipid levels mapped to meQTL.5a, and there was a modest linkage with lifespan. DNA binding motif and protein-protein interaction enrichment analyses identified the hepatic nuclear factor, Hnf1a (MODY3 gene in humans), as a strong candidate. The pleiotropic effects of meQTL.5a could contribute to variations in body size and metabolic traits, and influence CpG methylation and epigenetic ageing that could have an impact on lifespan.

Association of DNA methylation signatures with premature ageing and cardiovascular death in patients with end-stage kidney disease: a pilot epigenome-wide association study

Patients with end-stage kidney disease (ESKD) display features of premature aging. There is strong evidence that changes in DNA methylation (DNAm) contribute to age-related pathologies; however, little is known about their association with premature aging and cardiovascular mortality in patients with ESKD. We assayed genome-wide DNAm in a pilot case-control study of 60 hemodialysis patients with (n=30, cases) and without (n=30, controls) a fatal cardiovascular event. DNAm was profiled on the Illumina EPIC BeadChip. Four established DNAm clocks (i.e., Horvath-, Hannum-, Pheno-, and GrimAge) were used to estimate epigenetic age (DNAmAge). Epigenetic age acceleration (EAA) was derived as the residuals of regressing DNAmAge on chronological age (chroAge), and its association with cardiovascular death was examined using multivariable conditional logistic regression. An epigenome-wide association study (EWAS) was performed to identify differentially methylated CpGs associated with cardiovascular death. All clocks performed well at predicting chroAge (correlation between DNAmAges and chroAge of r=0.76-0.89), with GrimAge showing the largest deviation from chroAge (a mean of +21.3 years). There was no significant association of EAAs with cardiovascular death. In the EWAS, a CpG (cg22305782) in the FBXL19 gene had the strongest association with cardiovascular death with significantly lower DNAm in cases vs. controls (P_FDR=2.0x10^-6). FBXL19 is involved in cell apoptosis, inflammation, and adipogenesis. Overall, we observed more accelerated aging in patients with ESKD, although there was no significant association of EAAs with cardiovascular death. EWAS suggests a potential novel DNAm biomarker for premature cardiovascular mortality in ESKD.

Polyglucosan body density in the aged mouse hippocampus is controlled by a novel modifier locus on chromosome 1

Aging can be associated with the accumulation of hypobranched glycogen molecules (polyglucosan bodies, PGBs), particularly in astrocytes of the hippocampus. While PGBs have a detrimental effect on cognition in diseases such as adult polyglucosan body disease and Lafora disease, the underlying mechanism and clinical relevance of age-related PGB accumulation remains unknown. Here, we have investigated the genetic basis and functional impact of age-related PGB accumulation in 32 fully sequenced BXD-type strains of mice which exhibit a 400-fold variation in PGB burden in 16-18 month old females. We mapped a major locus controlling PGB density in the hippocampus to chromosome 1 at 72-75 Mb (linkage of 4.9 -logP), which we defined as the Pgb1 locus. To identify potentially causal gene variants within Pgb1, we generated extensive hippocampal transcriptome datasets and identified two strong candidate genes for which mRNA correlates with PGB density-Smarcal1 and Usp37. In addition, both Smarcal1 and Usp37 contain non-synonymous allele variations likely to impact protein function. A phenome-wide association analysis highlighted a trans-regulatory effect of the Pgb1 locus on expression of Hp1bp3, a gene known to play a role in age-related changes in learning and memory. To investigate the potential impact of PGBs on cognition, we performed conditioned fear memory testing on strains displaying varying degrees of PGB burden, and a phenome-wide association scan of ~12,000 traits. Importantly, we did not find any evidence suggesting a negative impact of PGB burden on cognitive capacity. Taken together, we have identified a major modifier locus controlling PGB burden in the hippocampus and shed light on the genetic architecture and clinical relevance of this strikingly heterogeneous hippocampal phenotype.

Universal DNA methylation age across mammalian tissues

Aging, often considered a result of random cellular damage, can be accurately estimated using DNA methylation profiles, the foundation of pan-tissue epigenetic clocks. Here, we demonstrate the development of universal pan-mammalian clocks, using 11,754 methylation arrays from our Mammalian Methylation Consortium, which encompass 59 tissue types across 185 mammalian species. These predictive models estimate mammalian tissue age with high accuracy (r > 0.96). Age deviations correlate with human mortality risk, mouse somatotropic axis mutations and caloric restriction. We identified specific cytosines with methylation levels that change with age across numerous species. These sites, highly enriched in polycomb repressive complex 2-binding locations, are near genes implicated in mammalian development, cancer, obesity and longevity. Our findings offer new evidence suggesting that aging is evolutionarily conserved and intertwined with developmental processes across all mammals.

DNA methylation networks underlying mammalian traits

Using DNA methylation profiles (n = 15,456) from 348 mammalian species, we constructed phyloepigenetic trees that bear marked similarities to traditional phylogenetic ones. Using unsupervised clustering across all samples, we identified 55 distinct cytosine modules, of which 30 are related to traits such as maximum life span, adult weight, age, sex, and human mortality risk. Maximum life span is associated with methylation levels in HOXL subclass homeobox genes and developmental processes and is potentially regulated by pluripotency transcription factors. The methylation state of some modules responds to perturbations such as caloric restriction, ablation of growth hormone receptors, consumption of high-fat diets, and expression of Yamanaka factors. This study reveals an intertwined evolution of the genome and epigenome that mediates the biological characteristics and traits of different mammalian species.

Epigenetic analysis in a murine genetic model of Gulf War illness

Of the nearly 1 million military personnel who participated in the 1990-1991 Gulf War, between 25% and 35% became ill with what now is referred to as Gulf War Illness (GWI) by the Department of Defense. Symptoms varied from gastrointestinal distress to lethargy, memory loss, inability to concentrate, depression, respiratory, and reproductive problems. The symptoms have persisted for 30 years in those afflicted but the basis of the illness remains largely unknown. Nerve agents and other chemical exposures in the war zone have been implicated but the long-term effects of these acute exposures have left few if any identifiable signatures. The major aim of this study is to elucidate the possible genomic basis for the persistence of symptoms, especially of the neurological and behavioral effects. To address this, we performed a whole genome epigenetic analysis of the proposed cause of GWI, viz., exposure to organophosphate neurotoxicants combined with high circulating glucocorticoids in two inbred mouse strains, C57BL/6J and DBA/2J. The animals received corticosterone in their drinking water for 7 days followed by injection of diisopropylfluorophosphate, a nerve agent surrogate. Six weeks after DFP injection, the animals were euthanized and medial prefrontal cortex harvested for genome-wide DNA methylation analysis using high-throughput sequencing. We observed 67 differentially methylated genes, notably among them, Ttll7, Akr1c14, Slc44a4, and Rusc2, all related to different symptoms of GWI. Our results support proof of principle of genetic differences in the chronic effects of GWI-related exposures and may reveal why the disease has persisted in many of the now aging Gulf War veterans.

Characterization of xenobiotic and steroid disposition potential of human placental tissue and cell lines (BeWo, JEG-3, JAR, and HTR-8/SVneo) by quantitative proteomics

The placenta is a fetal organ that performs critical functions to maintain pregnancy and support fetal development, including metabolism and transport of xenobiotics and steroids between the maternal-fetal unit. In vitro placenta models are used to study xenobiotic and steroid disposition, but how well these models recapitulate the human placenta is not well understood. We first characterized the abundance of proteins involved in xenobiotic and steroid disposition in human placental tissue. In pooled human placenta, the following xenobiotic and steroid disposition proteins were detected (highest to lowest), i) enzymes: glutathione S-transferase P, carbonyl reductase 1, aldo-keto reductase 1B1, hydroxysteroid dehydrogenases (HSD3B1 and HSD11B1), aromatase, epoxide hydrolase 1 (EPHX1) and steryl-sulfatase, and ii) transporters: monocarboxylate transporters (MCT1 and 4), organic anion transporting polypeptide 2B1, organic anion transporter 4, and breast cancer resistance protein (BCRP). Then, the tissue proteomics data was compared with four placental cell lines (BeWo, JEG-3, JAR, and HTR-8/SVneo). The differential global proteomics analysis revealed that the tissue and cell lines shared 1420 cytosolic and 1186 membrane proteins. Although extravillous trophoblast and cytotrophoblast marker proteins were detected in all cell lines, only BeWo and JEG-3 cells expressed the syncytiotrophoblast marker, chorionic somatomammotropin hormone 1. BeWo and JEG-3 cells expressed most target proteins including aromatase, HSDs, EPHX1, MCT1, and BCRP. JEG-3 cells treated with phthalates commonly detected in human biofluids, showed dysregulation of steroid pathways. The data presented here showed that BeWo and JEG-3 cells are closer to the placental tissue for studying xenobiotic and steroid disposition. Significance Statement This is the first study to compare proteomics data of human placental tissue and cell lines (BeWo, JAR, JEG-3, and HTR-8/SVneo). The placental cell line and tissue proteomes are vastly different, but BeWo and JEG-3 cells showed greater resemblance to the tissue in the expression of xenobiotic and steroid disposition proteins. These data will assist researchers to select an optimum cell model for mechanistic investigations on xenobiotic and steroid disposition in the placenta.

Changes in Dietary Intake of Methionine, Folate/Folic Acid and Vitamin B12 and Survival in Postmenopausal Women with Breast Cancer: A Prospective Cohort Study

BACKGROUND

Previous experimental studies showed that limiting methionine in the diet of animals or in cell culture media suppresses mammary cancer cell proliferation or metastasis. However, no previous study has investigated the associations of changes in methionine intake with survival among breast cancer survivors. We aimed to examine the association between changes in dietary intake of methionine, folate/folic acid, and vitamin B12 from before to after diagnosis of breast cancer, and mortality among breast cancer survivors.

METHODS

We included 1553 postmenopausal women from the Women's Health Initiative who were diagnosed with invasive breast cancer and completed a food frequency questionnaire both before and after breast cancer diagnosis. Multivariable Cox proportional hazards regression models were used to estimate adjusted hazard ratios (HRs) and 95% confidence (CIs) of all-cause and breast cancer mortality associated with changes in methionine intake and changes in folate/folic acid and vitamin B12 intake.

RESULTS

Relative to pre-diagnosis, 28% of women decreased methionine intake by ≥20%, 30% of women increased methionine intake by ≥20%, and 42% of women had a relatively stable methionine intake (±19.9%) following breast cancer diagnosis. During a mean 16.1 years of follow up, there were 772 deaths in total, including 195 deaths from breast cancer. Compared to women with relatively stable methionine intake, women with decreased methionine intake had lower risks of all-cause (HR 0.78, 95% CI 0.62-0.97) and breast cancer mortality (HR 0.58, 95% CI 0.37-0.91) in fully adjusted models. In contrast, increased methionine intake or changes in folate/folic acid or vitamin B12 intake were not associated with all-cause or breast cancer mortality.

CONCLUSIONS

Among breast cancer survivors, decreased methionine intake after breast cancer diagnosis was associated with lower risk of all-cause and breast cancer mortality.

Genetic loci and metabolic states associated with murine epigenetic aging

Changes in DNA methylation (DNAm) are linked to aging. Here, we profile highly conserved CpGs in 339 predominantly female mice belonging to the BXD family for which we have deep longevity and genomic data. We use a 'pan-mammalian' microarray that provides a common platform for assaying the methylome across mammalian clades. We computed epigenetic clocks and tested associations with DNAm entropy, diet, weight, metabolic traits, and genetic variation. We describe the multifactorial variance of methylation at these CpGs and show that high-fat diet augments the age-related changes. Entropy increases with age. The progression to disorder, particularly at CpGs that gain methylation over time, was predictive of genotype-dependent life expectancy. The longer-lived BXD strains had comparatively lower entropy at a given age. We identified two genetic loci that modulate epigenetic age acceleration (EAA): one on chromosome (Chr) 11 that encompasses the Erbb2/Her2 oncogenic region, and the other on Chr19 that contains a cytochrome P450 cluster. Both loci harbor genes associated with EAA in humans, including STXBP4, NKX2-3, and CUTC. Transcriptome and proteome analyses revealed correlations with oxidation-reduction, metabolic, and immune response pathways. Our results highlight concordant loci for EAA in humans and mice, and demonstrate a tight coupling between the metabolic state and epigenetic aging.

Gene-by-environment modulation of lifespan and weight gain in the murine BXD family

How lifespan and body weight vary as a function of diet and genetic differences is not well understood. Here we quantify the impact of differences in diet on lifespan in a genetically diverse family of female mice, split into matched isogenic cohorts fed a low-fat chow diet (CD, n = 663) or a high-fat diet (HFD, n = 685). We further generate key metabolic data in a parallel cohort euthanized at four time points. HFD feeding shortens lifespan by 12%: equivalent to a decade in humans. Initial body weight and early weight gains account for longevity differences of roughly 4-6 days per gram. At 500 days, animals on a HFD typically gain four times as much weight as control, but variation in weight gain does not correlate with lifespan. Classic serum metabolites, often regarded as health biomarkers, are not necessarily strong predictors of longevity. Our data indicate that responses to a HFD are substantially modulated by gene-by-environment interactions, highlighting the importance of genetic variation in making accurate individualized dietary recommendations.

Body weight and high-fat diet are associated with epigenetic aging in female members of the BXD murine family

DNA methylation (DNAm) is shaped by genetic and environmental factors and modulated by aging. Here, we examine interrelations between epigenetic aging, body weight (BW), and life span in 12 isogenic strains from the BXD family of mice that exhibit over twofold variation in longevity. Genome-wide DNAm was assayed in 70 liver specimens from predominantly female cases, 6-25 months old, that were maintained on normal chow or high-fat diet (HFD). We defined subsets of CpG regions associated with age, BW at young adulthood, and strain-by-diet-dependent life span. These age-associated differentially methylated CpG regions (age-DMRs) featured distinct genomic characteristics, with DNAm gains over time occurring in sites such as promoters and exons that have high CpG density and low average methylation. CpG regions associated with BW were enriched in introns, tended to have lower methylation in mice with higher BW, and were inversely correlated with gene expression (i.e., higher mRNA levels in mice with higher BW). CpG regions associated with life span were linked to genes involved in life span modulation, including the telomerase reverse transcriptase gene, Tert, which had both lower methylation and higher expression in long-lived strains. An epigenetic clock defined from age-DMRs revealed accelerated aging in mice belonging to strains with shorter life spans. Both higher BW and the HFD were associated with accelerated epigenetic aging. Our results highlight the age-accelerating effect of heavier BW. Furthermore, we demonstrate that the measure of epigenetic aging derived from age-DMRs can predict genotype and diet-induced differences in life span among female BXD members.

Effect of short-term prescription opioids on DNA methylation of the OPRM1 promoter

Background

A long-term opioid use has been associated with hypermethylation of the opioid receptor mu 1 (OPRM1) promoter. Very little is currently known about the early epigenetic response to therapeutic opioids. Here, we examine whether we can detect DNA methylation changes associated with a few days’ use of prescribed opioids. Genome-wide DNA methylation was assayed in a cohort of 33 opioid-naïve participants who underwent standard dental surgery followed by opioid self-administration. Saliva samples were collected before surgery (visit 1), and at two postsurgery visits at 2.7 ± 1.5 days (visit 2), and 39 ± 10 days (visit 3) after the discontinuation of opioid analgesics.

Results

The perioperative methylome underwent significant changes over the three visits that were primarily due to postoperative inflammatory response and cell heterogeneity. To specifically examine the effect of opioids, we started with a candidate gene approach and evaluated 10 CpGs located in the OPRM1 promoter. There was a significant cross-sectional variability in opioid use, and for participants who self-administered the prescribed drugs, the total dosage ranged from 5–210 morphine milligram equivalent (MME). Participants were categorized by cumulative dosage into three groups: < 25 MME, 25–90 MME, and ≥ 90 MME. Using mixed-effects modeling, 4 CpGs had significant positive associations with opioid dose at two-tailed p value < 0.05, and overall, 9 of the 10 OPRM1 promoter CpGs showed the predicted higher methylation in the higher dose groups relative to the lowest dose group. After adjustment for age, cellular heterogeneity, and past tobacco use, the promoter mean methylation also had positive associations with cumulative MME (regression coefficient = 0.0002, one-tailed p value = 0.02) and duration of opioid use (regression coefficient = 0.003, one-tailed p value = 0.001), but this effect was significant only for visit 3. A preliminary epigenome-wide association study identified a significant CpG in the promoter of the RAS-related signaling gene, RASL10A, that may be predictive of opioid dosage.

Conclusion

The present study provides evidence that the hypermethylation of the OPRM1 promoter is in response to opioid use and that epigenetic differences in OPRM1 and other sites are associated with a short-term use of therapeutic opioids.

An Application of the Bayesian Periodicity Test to Identify Diurnal Rhythm Genes in the Brain

Biological systems are extremely dynamic and many aspects of cellular processes show rhythmic circadian patterns. Extracting such information from large expression data is challenging. In this work, we present a modified application of the Empirical Bayes periodicity test to identify genes with diurnal rhythmic behavior in two brain regions. The hypothalamus and amygdala gene expression data were generated from 100 BXD recombinant inbred mice during the day hours. Brain samples were collected over the course of two days. We first filtered the transcripts based on rank correlation at matched time points between day-1 and day-2. We then applied the proposed test of periodicity to identify diurnal rhythm genes in the full cohort and gender-specific sub-cohorts. In hypothalamus, at a Benjamini-Hochberg false discovery rate (BH-FDR) of 0.01, we identified 15 transcripts with cyclic behavior in the full cohort, none, and 53 transcripts in the female and male cohort, respectively. Similarly, in amygdala, we identified 58 diurnal rhythm genes in the full cohort, and 1 and 28 in the female and male cohorts, respectively. In conclusion, we present a modified version of the empirical Bayes periodicity test to detect periodic expression patterns. Our results demonstrate that this approach can capture cyclic patterns from relatively noisy expression data sets.

Cognitive resilience among APOE ε4 carriers in the oldest old

OBJECTIVES

Relatively few APOE ε4+ carriers survive to old age (age 80+) without cognitive impairment (CI); thus, little is known about distinguishing characteristics of resilient APOE ε4+ carriers. Herein, we describe the sociodemographic characteristics of a large sample of resilient APOE ε4+ women from the Women's Health Initiative Memory Study (WHIMS) and compare them to noncarriers and APOE ε4+ women who developed CI before age 80.

METHODS

Women were recruited for clinical trials evaluating postmenopausal hormone therapy and incidence of dementia. During posttrial follow-up, cognitive status was adjudicated annually. Among 5716 women, we compared groups by APOE ε4 status using logistic regression, covarying for treatment, demographics, lifestyle, cardiovascular and physical function, well-being, and self-rated general health.

RESULTS

Among 557 APOE ε4+ women, those who survived to age 80+ without CI had higher baseline self-rated general health (odds ratio [OR]: 1.02; 95% confidence interval [CI], 1.01-1.04) and cognitive scores (OR: 1.18; 95% CI, 1.12-1.25) than those who did not reach age 80 without CI. Baseline high total cholesterol and low-density lipoprotein (LDL) levels were similar across APOE ε4+ groups but were higher compared with APOE ε4- women. Among women who survived to 80+ without CI, more APOE ε4+ women had a history of high total cholesterol (P = .003) and LDL cholesterol (OR: 1.01; 95% CI, 1.00-1.01). There were no differences in hypertension, diabetes, or other vascular risk factors in APOE ε4+ women compared with noncarriers.

CONCLUSIONS

Results highlight the importance of baseline cognitive function and general health for late-life cognition among ε4+ women.

Neutrophils-related host factors associated with severe disease and fatality in patients with influenza infection

Severe influenza infection has no effective treatment available. One of the key barriers to developing host-directed therapy is a lack of reliable prognostic factors needed to guide such therapy. Here, we use a network analysis approach to identify host factors associated with severe influenza and fatal outcome. In influenza patients with moderate-to-severe diseases, we uncover a complex landscape of immunological pathways, with the main changes occurring in pathways related to circulating neutrophils. Patients with severe disease display excessive neutrophil extracellular traps formation, neutrophil-inflammation and delayed apoptosis, all of which have been associated with fatal outcome in animal models. Excessive neutrophil activation correlates with worsening oxygenation impairment and predicted fatal outcome (AUROC 0.817-0.898). These findings provide new evidence that neutrophil-dominated host response is associated with poor outcomes. Measuring neutrophil-related changes may improve risk stratification and patient selection, a critical first step in developing host-directed immune therapy.

Longitudinal study of leukocyte DNA methylation and biomarkers for cancer risk in older adults

Background

Changes in DNA methylation over the course of life may provide an indicator of risk for cancer. We explored longitudinal changes in CpG methylation from blood leukocytes, and likelihood of future cancer diagnosis.

Methods

Peripheral blood samples were obtained at baseline and at follow-up visit from 20 participants in the Health, Aging and Body Composition prospective cohort study. Genome-wide CpG methylation was assayed using the Illumina Infinium Human MethylationEPIC (HM850K) microarray.

Results

Global patterns in DNA methylation from CpG-based analyses showed extensive changes in cell composition over time in participants who developed cancer. By visit year 6, the proportion of CD8+ T-cells decreased (p-value = 0.02), while granulocytes cell levels increased (p-value = 0.04) among participants diagnosed with cancer compared to those who remained cancer-free (cancer-free vs. cancer-present: 0.03 ± 0.02 vs. 0.003 ± 0.005 for CD8+ T-cells; 0.52 ± 0.14 vs. 0.66 ± 0.09 for granulocytes). Epigenome-wide analysis identified three CpGs with suggestive p-values ≤10^− 5 for differential methylation between cancer-free and cancer-present groups, including a CpG located in MTA3, a gene linked with metastasis. At a lenient statistical threshold (p-value ≤3 × 10^− 5), the top 10 cancer-associated CpGs included a site near RPTOR that is involved in the mTOR pathway, and the candidate tumor suppressor genes REC8, KCNQ1, and ZSWIM5. However, only the CpG in RPTOR (cg08129331) was replicated in an independent data set. Analysis of within-individual change from baseline to Year 6 found significant correlations between the rates of change in methylation in RPTOR, REC8 and ZSWIM5, and time to cancer diagnosis.

Conclusion

The results show that changes in cellular composition explains much of the cross-sectional and longitudinal variation in CpG methylation. Additionally, differential methylation and longitudinal dynamics at specific CpGs could provide powerful indicators of cancer development and/or progression. In particular, we highlight CpG methylation in the RPTOR gene as a potential biomarker of cancer that awaits further validation.

Electronic supplementary material

The online version of this article (10.1186/s40364-019-0161-3) contains supplementary material, which is available to authorized users.

Fetal-Not Maternal-APOL1 Genotype Associated with Risk for Preeclampsia in Those with African Ancestry

Black Americans are at increased risk for preeclampsia. Genetic variants in apolipoprotein L1 (APOL1) account for much of the increased risk for kidney disease in blacks. APOL1 is expressed in human placenta and transgenic mice expressing APOL1 develop preeclampsia. We evaluated the role of APOL1 variants in human preeclampsia. We determined maternal and fetal APOL1 genotypes in black women with preeclampsia in two populations. At Einstein Montefiore Center (EMC) Affiliated Hospitals, we studied 121 pregnancies in black women with preeclampsia. At University of Tennessee Health Science Center (UTHSC), we studied 93 pregnancies in black women with preeclampsia and 793 pregnancies without preeclampsia. We measured serum markers of preeclampsia soluble fms-like tyrosine kinase 1 (sFlt-1), placental growth factor (PlGF), and soluble endoglin (sEng). Fetal APOL1 high-risk (HR) genotype was associated with preeclampsia, with odds ratios at EMC and UTHSC of 1.84 (95% CI 1.11, 2.93) and 1.92 (95% CI 1.05, 3.49), respectively. Maternal APOL1 HR genotype was not associated with preeclampsia. Mothers with the fetal APOL1 HR genotype had more cerebral or visual disturbances (63% versus 37%, p = 0.04). In addition, fetal APOL1 HR genotype was associated with a higher sFLT-1/PlGF ratio at birth (p = 0.04). Fetal APOL1 high-risk genotype increases the risk for preeclampsia, likely by adversely affecting placental function. Further research is needed to assess whether APOL1 genetic testing can predict preeclampsia and improve pregnancy outcomes.

Genetic cartography of longevity in humans and mice: Current landscape and horizons

Aging is a complex and highly variable process. Heritability of longevity among humans and other species is low, and this finding has given rise to the idea that it may be futile to search for DNA variants that modulate aging. We argue that the problem in mapping longevity genes is mainly one of low power and the genetic and environmental complexity of aging. In this review we highlight progress made in mapping genes and molecular networks associated with longevity, paying special attention to work in mice and humans. We summarize 40 years of linkage studies using murine cohorts and 15 years of studies in human populations that have exploited candidate gene and genome-wide association methods. A small but growing number of gene variants contribute to known longevity mechanisms, but a much larger set have unknown functions. We outline these and other challenges and suggest some possible solutions, including more intense collaboration between research communities that use model organisms and human cohorts. Once hundreds of gene variants have been linked to differences in longevity in mammals, it will become feasible to systematically explore gene-by-environmental interactions, dissect mechanisms with more assurance, and evaluate the roles of epistasis and epigenetics in aging. A deeper understanding of complex networks-genetic, cellular, physiological, and social-should position us well to improve healthspan.

Profiling DNA methylation differences between inbred mouse strains on the Illumina Human Infinium MethylationEPIC microarray

The Illumina Infinium MethylationEPIC provides an efficient platform for profiling DNA methylation in humans at over 850,000 CpGs. Model organisms such as mice do not currently benefit from an equivalent array. Here we used this array to measure DNA methylation in mice. We defined probes targeting conserved regions and performed differential methylation analysis and compared between the array-based assay and affinity-based DNA sequencing of methyl-CpGs (MBD-seq) and reduced representation bisulfite sequencing. Mouse samples consisted of 11 liver DNA from two strains, C57BL/6J (B6) and DBA/2J (D2), that varied widely in age. Linear regression was applied to detect differential methylation. In total, 13,665 probes (1.6% of total probes) aligned to conserved CpGs. Beta-values (β-value) for these probes showed a distribution similar to that in humans. Overall, there was high concordance in methylation signal between the EPIC array and MBD-seq (Pearson correlation r = 0.70, p-value < 0.0001). However, the EPIC probes had higher quantitative sensitivity at CpGs that are hypo- (β-value < 0.3) or hypermethylated (β-value > 0.7). In terms of differential methylation, no EPIC probe detected a significant difference between age groups at a Benjamini-Hochberg threshold of 10%, and the MBD-seq performed better at detecting age-dependent change in methylation. However, the top most significant probe for age (cg13269407; uncorrected p-value = 1.8 x 10-5) is part of the clock CpGs used to estimate the human epigenetic age. For strain, 219 EPIC probes detected significant differential methylation (FDR cutoff 10%) with ~80% CpGs associated with higher methylation in D2. This higher methylation profile in D2 compared to B6 was also replicated by the MBD-seq data. To summarize, we found only a small subset of EPIC probes that target conserved sites. However, for this small subset the array provides a reliable assay of DNA methylation and can be effectively used to measure differential methylation in mice.

Genetic Analysis of Mitochondrial Ribosomal Proteins and Cognitive Aging in Postmenopausal Women

Genes encoding mitochondrial ribosomal proteins (MRPs) have been linked to aging and longevity in model organisms (i.e., mice, Caenorhabditis elegans). Here we evaluated if the MRPs have conserved effects on aging traits in humans. We utilized data from 4,504 participants of the Women's Health Initiative Memory Study (WHIMS) who had both longitudinal cognitive data and genetic data. Two aging phenotypes were considered: (1) gross lifespan (time to all-cause mortality), and (2) cognitive aging (longitudinal rate of change in modified mini-mental state scores). We tested genetic association with variants in 78 members of the MRP gene family. Genetic association tests were done at the single nucleotide polymorphism (SNP) level, and at gene-set level using two distinct procedures (GATES and MAGMA). We included SNPs in APOE and adjusted the tests for the APOE-ε4 allele, a known risk factor for dementia. The strongest association signal is for the known cognitive aging SNP, rs429358, in APOE (p-value = 5 × 10⁻²⁸ for cognitive aging; p-value = 0.03 for survival). We found no significant association between the MRPs and survival time. For cognitive aging, we detected SNP level association for rs189661478 in MRPL23 (p-value < 9 × 10⁻⁶). Furthermore, the gene-set analysis showed modest but significant association between the MRP family and cognitive aging. In conclusion, our results indicate a potential pathway-level association between the MRPs and cognitive aging that is independent of the APOE locus. We however did not detect association between the MRPs and lifespan.

Patients experiencing statin-induced myalgia exhibit a unique program of skeletal muscle gene expression following statin re-challenge

Statins, the 3-hydroxy-3-methyl-glutaryl (HMG)-CoA reductase inhibitors, are widely prescribed for treatment of hypercholesterolemia. Although statins are generally well tolerated, up to ten percent of statin-treated patients experience myalgia symptoms, defined as muscle pain without elevated creatinine phosphokinase (CPK) levels. Myalgia is the most frequent reason for discontinuation of statin therapy. The mechanisms underlying statin myalgia are not clearly understood. To elucidate changes in gene expression associated with statin myalgia, we compared profiles of gene expression in skeletal muscle biopsies from patients with statin myalgia who were undergoing statin re-challenge (cases) versus those of statin-tolerant controls. A robust separation of case and control cohorts was revealed by Principal Component Analysis of differentially expressed genes (DEGs). To identify putative gene expression and metabolic pathways that may be perturbed in skeletal muscles of patients with statin myalgia, we subjected DEGs to Ingenuity Pathways (IPA) and DAVID (Database for Annotation, Visualization and Integrated Discovery) analyses. The most prominent pathways altered by statins included cellular stress, apoptosis, cell senescence and DNA repair (TP53, BARD1, Mre11 and RAD51); activation of pro-inflammatory immune response (CXCL12, CST5, POU2F1); protein catabolism, cholesterol biosynthesis, protein prenylation and RAS-GTPase activation (FDFT1, LSS, TP53, UBD, ATF2, H-ras). Based on these data we tentatively conclude that persistent myalgia in response to statins may emanate from cellular stress underpinned by mechanisms of post-inflammatory repair and regeneration. We also posit that this subset of individuals is genetically predisposed to eliciting altered statin metabolism and/or increased end-organ susceptibility that lead to a range of statin-induced myopathies. This mechanistic scenario is further bolstered by the discovery that a number of single nucleotide polymorphisms (e.g., SLCO1B1, SLCO2B1 and RYR2) associated with statin myalgia and myositis were observed with increased frequency among patients with statin myalgia.

Genetic divergence in the transcriptional engram of chronic alcohol abuse: A laser-capture RNA-seq study of the mouse mesocorticolimbic system

Genetic factors that influence the transition from initial drinking to dependence remain enigmatic. Recent studies have leveraged chronic intermittent ethanol (CIE) paradigms to measure changes in brain gene expression in a single strain at 0, 8, 72 h, and even 7 days following CIE. We extend these findings using LCM RNA-seq to profile expression in 11 brain regions in two inbred strains - C57BL/6J (B6) and DBA/2J (D2) - 72 h following multiple cycles of ethanol self-administration and CIE. Linear models identified differential expression based on treatment, region, strain, or interactions with treatment. Nearly 40% of genes showed a robust effect (FDR < 0.01) of region, and hippocampus CA1, cortex, bed nucleus stria terminalis, and nucleus accumbens core had the highest number of differentially expressed genes after treatment. Another 8% of differentially expressed genes demonstrated a robust effect of strain. As expected, based on similar studies in B6, treatment had a much smaller impact on expression; only 72 genes (p < 0.01) are modulated by treatment (independent of region or strain). Strikingly, many more genes (415) show a strain-specific and largely opposite response to treatment and are enriched in processes related to RNA metabolism, transcription factor activity, and mitochondrial function. Over 3 times as many changes in gene expression were detected in D2 compared to B6, and weighted gene co-expression network analysis (WGCNA) module comparison identified more modules enriched for treatment effects in D2. Substantial strain differences exist in the temporal pattern of transcriptional neuroadaptation to CIE, and these may drive individual differences in risk of addiction following excessive alcohol consumption.

Hippocampal proteomics defines pathways associated with memory decline and resilience in normal aging and Alzheimer's disease mouse models

Alzheimer's disease (AD), the most common form of dementia in the elderly, has no cure. Thus, the identification of key molecular mediators of cognitive decline in AD remains a top priority. As aging is the most significant risk factor for AD, the goal of this study was to identify altered proteins and pathways associated with the development of normal aging and AD memory deficits, and identify unique proteins and pathways that may contribute to AD-specific symptoms. We used contextual fear conditioning to diagnose 8-month-old 5XFAD and non-transgenic (Ntg) mice as having either intact or impaired memory, followed by liquid chromatography-tandem mass spectrometry (LC-MS/MS) to quantify hippocampal membrane proteins across groups. Subsequent analysis detected 113 proteins differentially expressed relative to memory status (intact vs impaired) in Ntg mice and 103 proteins in 5XFAD mice. Thirty-six proteins, including several involved in neuronal excitability and synaptic plasticity (e.g., GRIA1, GRM3, and SYN1), were altered in both normal aging and AD. Pathway analysis highlighted HDAC4 as a regulator of observed protein changes in both genotypes and identified the REST epigenetic regulatory pathway and G_i intracellular signaling as AD-specific pathways involved in regulating the onset of memory deficits. Comparing the hippocampal membrane proteome of Ntg versus AD, regardless of cognitive status, identified 138 differentially expressed proteins, including confirmatory proteins APOE and CLU. Overall, we provide a novel list of putative targets and pathways with therapeutic potential, including a set of proteins associated with cognitive status in normal aging mice or gene mutations that cause AD.

Joint mouse-human phenome-wide association to test gene function and disease risk

Phenome-wide association is a novel reverse genetic strategy to analyze genome-to-phenome relations in human clinical cohorts. Here we test this approach using a large murine population segregating for ∼5 million sequence variants, and we compare our results to those extracted from a matched analysis of gene variants in a large human cohort. For the mouse cohort, we amassed a deep and broad open-access phenome consisting of ∼4,500 metabolic, physiological, pharmacological and behavioural traits, and more than 90 independent expression quantitative trait locus (QTL), transcriptome, proteome, metagenome and metabolome data sets—by far the largest coherent phenome for any experimental cohort (www.genenetwork.org). We tested downstream effects of subsets of variants and discovered several novel associations, including a missense mutation in fumarate hydratase that controls variation in the mitochondrial unfolded protein response in both mouse and Caenorhabditis elegans, and missense mutations in Col6a5 that underlies variation in bone mineral density in both mouse and human.

Phenome-wide association is a novel method that links sequence variants to a spectrum of phenotypes and diseases. Here the authors generate detailed mouse genetic and phenome data which links their phenome-wide association study (PheWAS) of mouse to corresponding PheWAS in human.

Ancestry dependent DNA methylation and influence of maternal nutrition

There is extensive variation in DNA methylation between individuals and ethnic groups. These differences arise from a combination of genetic and non-genetic influences and potential modifiers include nutritional cues, early life experience, and social and physical environments. Here we compare genome-wide DNA methylation in neonatal cord blood from African American (AA; N = 112) and European American (EA; N = 91) participants of the CANDLE Study (Conditions Affecting Neurocognitive Development and Learning in Early Childhood). Our goal is to determine if there are replicable ancestry-specific methylation patterns that may implicate risk factors for diseases that have differential prevalence between populations. To identify the most robust ancestry-specific CpG sites, we replicate our results in lymphoblastoid cell lines from Yoruba African and CEPH European panels of HapMap. We also evaluate the influence of maternal nutrition—specifically, plasma levels of vitamin D and folate during pregnancy—on methylation in newborns. We define stable ancestry-dependent methylation of genes that include tumor suppressors and cell cycle regulators (e.g., APC, BRCA1, MCC). Overall, there is lower global methylation in African ancestral groups. Plasma levels of 25-hydroxy vitamin D are also considerably lower among AA mothers and about 60% of AA and 40% of EA mothers have concentrations below 20 ng/ml. Using a weighted correlation analysis, we define a network of CpG sites that is jointly modulated by ancestry and maternal vitamin D. Our results show that differences in DNA methylation patterns are remarkably stable and maternal micronutrients can exert an influence on the child epigenome.

Sex-specific modulation of gene expression networks in murine hypothalamus

The hypothalamus contains nuclei and cell populations that are critical in reproduction and that differ significantly between the sexes in structure and function. To examine the molecular and genetic basis for these differences, we quantified gene expression in the hypothalamus of 39 pairs of adult male and female mice belonging to the BXD strains. This experimental design enabled us to define hypothalamic gene coexpression networks and provided robust estimates of absolute expression differences. As expected, sex has the strongest effect on the expression of genes on the X and Y chromosomes (e.g., Uty, Xist, Kdm6a). Transcripts associated with the endocrine system and neuropeptide signaling also differ significantly. Sex-differentiated transcripts often have well delimited expression within specific hypothalamic nuclei that have roles in reproduction. For instance, the estrogen receptor (Esr1) and neurokinin B (Tac2) genes have intense expression in the medial preoptic and arcuate nuclei and comparatively high expression in females. Despite the strong effect of sex on single transcripts, the global pattern of covariance among transcripts is well preserved, and consequently, males and females have well matched coexpression modules. However, there are sex-specific hub genes in functionally equivalent modules. For example, only in males is the Y-linked gene, Uty, a highly connected transcript in a network that regulates chromatin modification and gene transcription. In females, the X chromosome paralog, Kdm6a, takes the place of Uty in the same network. We also find significant effect of sex on genetic regulation and the same network in males and females can be associated with markedly different regulatory loci. With the exception of a few sex-specific modules, our analysis reveals a system in which sets of functionally related transcripts are organized into stable sex-independent networks that are controlled at a higher level by sex-specific modulators.

Complex control of GABA(A) receptor subunit mRNA expression: variation, covariation, and genetic regulation

GABA type-A receptors are essential for fast inhibitory neurotransmission and are critical in brain function. Surprisingly, expression of receptor subunits is highly variable among individuals, but the cause and impact of this fluctuation remains unknown. We have studied sources of variation for all 19 receptor subunits using massive expression data sets collected across multiple brain regions and platforms in mice and humans. Expression of Gabra1, Gabra2, Gabrb2, Gabrb3, and Gabrg2 is highly variable and heritable among the large cohort of BXD strains derived from crosses of fully sequenced parents—C57BL/6J and DBA/2J. Genetic control of these subunits is complex and highly dependent on tissue and mRNA region. Remarkably, this high variation is generally not linked to phenotypic differences. The single exception is Gabrb3, a locus that is linked to anxiety. We identified upstream genetic loci that influence subunit expression, including three unlinked regions of chromosome 5 that modulate the expression of nine subunits in hippocampus, and that are also associated with multiple phenotypes. Candidate genes within these loci include, Naaa, Nos1, and Zkscan1. We confirmed a high level of coexpression for subunits comprising the major channel—Gabra1, Gabrb2, and Gabrg2—and identified conserved members of this expression network in mice and humans. Gucy1a3, Gucy1b3, and Lis1 are novel and conserved associates of multiple subunits that are involved in inhibitory signaling. Finally, proximal and distal regions of the 3′ UTRs of single subunits have remarkably independent expression patterns in both species. However, corresponding regions of different subunits often show congruent genetic control and coexpression (proximal-to-proximal or distal-to-distal), even in the absence of sequence homology. Our findings identify novel sources of variation that modulate subunit expression and highlight the extraordinary capacity of biological networks to buffer 4–100 fold differences in mRNA levels.

A promoter polymorphism in the Per3 gene is associated with alcohol and stress response

The period homolog genes Per1, Per2 and Per3 are important components of the circadian clock system. In addition to their role in maintaining circadian rhythm, these genes have been linked to mood disorders, stress response and vulnerability to addiction and alcoholism. In this study, we combined high-resolution sequence analysis and quantitative trait locus (QTL) mapping of gene expression and behavioral traits to identify Per3 as a compelling candidate for the interaction between circadian rhythm, alcohol and stress response. In the BXD family of mouse strains, sequence variants in Per3 have marked effects on steady-state mRNA and protein levels. As a result, the transcript maps as a cis-acting expression QTL (eQTL). We found that an insertion/deletion (indel) variant in a putative stress response element in the promoter region of Per3 causes local control of transcript abundance. This indel results in differences in protein binding affinities between the two alleles through the Nrf2 transcriptional activator. Variation in Per3 is also associated with downstream differences in the expression of genes involved in circadian rhythm, alcohol, stress response and schizophrenia. We found that the Per3 locus is linked to stress/anxiety traits, and that the basal expression of Per3 is also correlated with several anxiety and addiction-related phenotypes. Treatment with alcohol results in increased expression of Per3 in the hippocampus, and this effect interacts with acute restraint stress. Our data provide strong evidence that variation in the Per3 transcript is causally associated with and also responsive to stress and alcohol.

Genetic modulation of striatal volume by loci on Chrs 6 and 17 in BXD recombinant inbred mice

Natural variation in the absolute and relative size of different parts of the human brain is substantial, with a range that often exceeds a factor of 2. Much of this variation is generated by the cumulative effects of sets of unknown gene variants that modulate the proliferation, growth and death of neurons and glial cells. Discovering and testing the functions of these genes should contribute significantly to our understanding of differences in brain development, behavior and disease susceptibility. We have exploited a large population of genetically well-characterized strains of mice (BXD recombinant inbred strains) to map gene variants that influence the volume of the dorsal striatum (caudate-putamen without nucleus accumbens). We used unbiased methods to estimate volumes bilaterally in a sex-balanced sample taken from the Mouse Brain Library (www.mbl.org). We generated a matched microarray data set to efficiently evaluate candidate genes (www.genenetwork.org). As in humans, volume of the striatum is highly heritable, with greater than twofold differences among strains. We mapped a locus that modulates striatal volume on chromosome (Chr) 6 at 88 +/- 5 Mb. We also uncovered an epistatic interaction between loci on Chr 6 and Chr 17 that modulates striatal volume. Using bioinformatic tools and the corresponding expression database, we have identified positional candidates in these quantitative trait locus intervals.

Dissection of a QTL hotspot on mouse distal chromosome 1 that modulates neurobehavioral phenotypes and gene expression

A remarkably diverse set of traits maps to a region on mouse distal chromosome 1 (Chr 1) that corresponds to human Chr 1q21-q23. This region is highly enriched in quantitative trait loci (QTLs) that control neural and behavioral phenotypes, including motor behavior, escape latency, emotionality, seizure susceptibility (Szs1), and responses to ethanol, caffeine, pentobarbital, and haloperidol. This region also controls the expression of a remarkably large number of genes, including genes that are associated with some of the classical traits that map to distal Chr 1 (e.g., seizure susceptibility). Here, we ask whether this QTL-rich region on Chr 1 (Qrr1) consists of a single master locus or a mixture of linked, but functionally unrelated, QTLs. To answer this question and to evaluate candidate genes, we generated and analyzed several gene expression, haplotype, and sequence datasets. We exploited six complementary mouse crosses, and combed through 18 expression datasets to determine class membership of genes modulated by Qrr1. Qrr1 can be broadly divided into a proximal part (Qrr1p) and a distal part (Qrr1d), each associated with the expression of distinct subsets of genes. Qrr1d controls RNA metabolism and protein synthesis, including the expression of approximately 20 aminoacyl-tRNA synthetases. Qrr1d contains a tRNA cluster, and this is a functionally pertinent candidate for the tRNA synthetases. Rgs7 and Fmn2 are other strong candidates in Qrr1d. FMN2 protein has pronounced expression in neurons, including in the dendrites, and deletion of Fmn2 had a strong effect on the expression of few genes modulated by Qrr1d. Our analysis revealed a highly complex gene expression regulatory interval in Qrr1, composed of multiple loci modulating the expression of functionally cognate sets of genes.

Variation in mouse basolateral amygdala volume is associated with differences in stress reactivity and fear learning

A wealth of research identifies the amygdala as a key brain region mediating negative affect, and implicates amygdala dysfunction in the pathophysiology of anxiety disorders. Although there is a strong genetic component to anxiety disorders such as posttraumatic stress disorder (PTSD) there remains debate about whether abnormalities in amygdala function predispose to these disorders. In the present study, groups of C57BL/6 x DBA/2 (B x D) recombinant inbred strains of mice were selected for differences in volume of the basolateral amygdala complex (BLA). Strains with relatively small, medium, or large BLA volumes were compared for Pavlovian fear learning and memory, anxiety-related behaviors, depression-related behavior, and glucocorticoid responses to stress. Strains with relatively small BLA exhibited stronger conditioned fear responses to both auditory tone and contextual stimuli, as compared to groups with larger BLA. The small BLA group also showed significantly greater corticosterone responses to stress than the larger BLA groups. BLA volume did not predict clear differences in measures of anxiety-like behavior or depression-related behavior, other than greater locomotor inhibition to novelty in strains with smaller BLA. Neither striatal, hippocampal nor cerebellar volumes correlated significantly with any behavioral measure. The present data demonstrate a phenotype of enhanced fear conditioning and exaggerated glucocorticoid responses to stress associated with small BLA volume. This profile is reminiscent of the increased fear processing and stress reactivity that is associated with amygdala excitability and reduced amygdala volume in humans carrying loss of function polymorphisms in the serotonin transporter and monoamine oxidase A genes. Our study provides a unique example of how natural variation in amygdala volume associates with specific fear- and stress-related phenotypes in rodents, and further supports the role of amygdala dysfunction in anxiety disorders such as PTSD.

Genetic and Structural Analysis of the Basolateral Amygdala Complex in BXD Recombinant Inbred Mice

The amygdala integrates and coordinates emotional and autonomic responses. The genetics that underlie variation in amygdala structure may be coupled to variation in levels of aggression, fear, anxiety, and affiliated behaviors. We systematically quantified the volume and cell populations of the basolateral amygdala complex (BLAc) across 35 BXD recombinant inbred (RI) lines, the parental strains--C57BL/6J (B6) and DBA/2J (D2)--and F1 hybrids (n cases=199, bilateral analysis). Neuron number and volume vary 1.7- to 2-fold among strains (e.g., neuron number ranged from 88,000 to 170,000). Glial and endothelial populations ranged more widely (5- to 8-fold), in part because of higher technical error. A quantitative trait locus (QTL) for the BLAc size is located on chromosome (Chr) 8 near the Large gene. This locus may also influence volume of other regions including hippocampus and cerebellum. Cell populations in the BLAc appear to be modulated more weakly by loci on Chrs 11 and 13. Candidate genes were selected on the basis of correlation with BLAc traits, chromosomal location, single nucleotide polymorphism (SNP) density, and expression patterns in the Allen Brain Atlas. Neurod2, a gene shown to be significant for the formation of the BLAc by knockout studies, is among the candidates genes. Other candidates include Large, and Thra. Responses to drugs of abuse and locomotor activity were the most notable behavioral correlates of the BLAc traits.

Differential covariation in taste responsiveness to bitter stimuli in rats

Variation exists in the sensitivity of individual rodents and humans to different bitter tastants. An absence of uniform correlation in responsiveness to different bitter substances across individuals within a species suggests heterogeneity in the mechanisms underlying stimulus processing within this taste modality. Here, we examined taste responsiveness of individual rats to three bitter compounds (quinine hydrochloride, denatonium benzoate, and cycloheximide) in short-term lick tests to determine the magnitude of covariation among responses to these stimuli and infer commonalities in their receptor and neural mechanisms. Rats were tested with a given pair of bitter stimuli during three sessions comprising randomized trial blocks of six concentrations of each stimulus + deionized water. Psychophysical functions were generated for individual rats for respective stimulus pairs, and concentrations of each stimulus that produced equivalent lick suppression relative to water were correlated across animals. Behavioral taste responsiveness to quinine hydrochloride strongly covaried with responsiveness to denatonium benzoate (r = +0.82). Lick responsiveness to quinine was less robustly correlated with that to cycloheximide (r = +0.44), and denatonium and cycloheximide responses failed to correlate. These results imply substantial overlap in the bitter taste coding mechanisms for quinine and denatonium but some degree of independence in the mechanisms responsible for gustatory processing of cycloheximide. More generally, these data reinforce the notion that bitter taste processing is not a homogeneous event.