Sample matching by inferred agonal stress in gene expression analyses of the brain

BackgroundGene expression patterns in the brain are strongly influenced by the severity and duration of physiological stress at the time of death. This agonal effect, if not well controlled, can lead to spurious findings and diminished statistical power in case-control comparisons. While some recent studies match samples by tissue pH and clinically recorded agonal conditions, we found that these indicators were sometimes at odds with observed stress-related gene expression patterns, and that matching by these criteria still sometimes results in identifying case-control differences that are primarily driven by residual agonal effects. This problem is analogous to the one encountered in genetic association studies, where self-reported race and ethnicity are often imprecise proxies for an individual's actual genetic ancestry.ResultsWe developed an Agonal Stress Rating (ASR) system that evaluates each sample's degree of stress based on gene expression data, and used ASRs in post hoc sample matching or covariate analysis. While gene expression patterns are generally correlated across different brain regions, we found strong region-region differences in empirical ASRs in many subjects that likely reflect inter-individual variabilities in local structure or function, resulting in region-specific vulnerability to agonal stress.ConclusionVariation of agonal stress across different brain regions differs between individuals, revealing a new level of complexity for gene expression studies of brain tissues. The Agonal Stress Ratings quantitatively assess each sample's extent of regulatory response to agonal stress, and allow a strong control of this important confounder.

[1]  Carl D Langefeld,et al.  Ordered subset analysis in genetic linkage mapping of complex traits , 2004, Genetic epidemiology.

[2]  Ilya Shmulevich,et al.  In silico microdissection of microarray data from heterogeneous cell populations , 2005, BMC Bioinformatics.

[3]  R. Myers,et al.  Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data , 2005, Nucleic acids research.

[4]  J. Wang-Rodriguez,et al.  In silico dissection of cell-type-associated patterns of gene expression in prostate cancer. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[5]  Hugues Bersini,et al.  Separation of samples into their constituents using gene expression data , 2001, ISMB.

[6]  N. Cairns,et al.  Quantifying mRNA in postmortem human brain: influence of gender, age at death, postmortem interval, brain pH, agonal state and inter-lobe mRNA variance. , 2003, Brain research. Molecular brain research.

[7]  P. Nelson,et al.  Project normal: Defining normal variance in mouse gene expression , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[8]  K. Roeder,et al.  The power of genomic control. , 2000, American journal of human genetics.

[9]  David J Balding,et al.  Logistic regression protects against population structure in genetic association studies. , 2005, Genome research.

[10]  Pui-Yan Kwok,et al.  Population structure, admixture, and aging-related phenotypes in African American adults: the Cardiovascular Health Study. , 2005, American journal of human genetics.

[11]  P. Donnelly,et al.  Inference of population structure using multilocus genotype data. , 2000, Genetics.

[12]  Mark D Shriver,et al.  Control of confounding of genetic associations in stratified populations. , 2003, American journal of human genetics.

[13]  E G Jones,et al.  Dysregulation of the fibroblast growth factor system in major depression. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[14]  R. Tibshirani,et al.  Diagnosis of multiple cancer types by shrunken centroids of gene expression , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[15]  R. Williams,et al.  Gm3;5,13,14 and type 2 diabetes mellitus: an association in American Indians with genetic admixture. , 1988, American journal of human genetics.

[16]  Huda Akil,et al.  Effect of agonal and postmortem factors on gene expression profile: quality control in microarray analyses of postmortem human brain , 2004, Biological Psychiatry.

[17]  D. Reich,et al.  Detecting association in a case‐control study while correcting for population stratification , 2001, Genetic epidemiology.

[18]  Douglas A. Hosack,et al.  Identifying biological themes within lists of genes with EASE , 2003, Genome Biology.

[19]  R. Ross,et al.  Brain pH has a significant impact on human postmortem hippocampal gene expression profiles , 2006, Brain Research.

[20]  K. Roeder,et al.  Genomic Control for Association Studies , 1999, Biometrics.

[21]  Debashis Ghosh,et al.  Mixture models for assessing differential expression in complex tissues using microarray data , 2004, Bioinform..

[22]  N. Risch,et al.  Estimation of individual admixture: Analytical and study design considerations , 2005, Genetic epidemiology.

[23]  Huda Akil,et al.  Systematic changes in gene expression in postmortem human brains associated with tissue pH and terminal medical conditions. , 2004, Human molecular genetics.

[24]  P. Dodd,et al.  Biochemical and molecular studies using human autopsy brain tissue , 2003, Journal of neurochemistry.

[25]  T. Speed,et al.  Summaries of Affymetrix GeneChip probe level data. , 2003, Nucleic acids research.

[26]  Elizabeth L. Ogburn,et al.  Demonstrating stratification in a European American population , 2005, Nature Genetics.

[27]  R. Myers,et al.  Mitochondrial related gene expression changes are sensitive to agonal-pH state , 2006, Molecular Psychiatry.

[28]  P. Donnelly,et al.  The effects of human population structure on large genetic association studies , 2004, Nature Genetics.

[29]  Pak Sham,et al.  Properties of Structured Association Approaches to Detecting Population Stratification , 2005, Human Heredity.

[30]  Michael J Bamshad,et al.  Human population genetic structure and inference of group membership. , 2003, American journal of human genetics.

[31]  Jun-Feng Wang,et al.  Downregulation in components of the mitochondrial electron transport chain in the postmortem frontal cortex of subjects with bipolar disorder. , 2006, Journal of psychiatry & neuroscience : JPN.

[32]  Eric R. Ziegel,et al.  Biostatistical Genetics and Genetic Epidemiology , 2002, Technometrics.

[33]  Jacques Corbeil,et al.  Comparative gene expression analysis of blood and brain provides concurrent validation of SELENBP1 up-regulation in schizophrenia , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[34]  Stanley J. Watson,et al.  Methodological considerations for gene expression profiling of human brain , 2007, Journal of Neuroscience Methods.

[35]  R. Ward,et al.  Informativeness of genetic markers for inference of ancestry. , 2003, American journal of human genetics.

[36]  J. Mesirov,et al.  Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. , 1999, Science.

[37]  L. Criswell,et al.  Markers informative for ancestry demonstrate consistent megabase-length linkage disequilibrium in the African American population , 2003, Human Genetics.

[38]  K. Konvička,et al.  Matching strategies for genetic association studies in structured populations. , 2004, American journal of human genetics.

[39]  R. Yolken,et al.  Mitochondrial dysfunction in schizophrenia: evidence for compromised brain metabolism and oxidative stress , 2004, Molecular Psychiatry.

[40]  D. Reich,et al.  Principal components analysis corrects for stratification in genome-wide association studies , 2006, Nature Genetics.

[41]  P. Donnelly,et al.  Association mapping in structured populations. , 2000, American journal of human genetics.

[42]  R. Yolken,et al.  The Stanley Foundation brain collection and Neuropathology Consortium , 2000, Schizophrenia Research.

[43]  Elena A. Erosheva,et al.  Grade of membership and latent structure models with application to disability survey data , 2002 .

[44]  Kazuya Iwamoto,et al.  Altered expression of mitochondria-related genes in postmortem brains of patients with bipolar disorder or schizophrenia, as revealed by large-scale DNA microarray analysis. , 2005, Human molecular genetics.

[45]  G A Satten,et al.  Accounting for unmeasured population substructure in case-control studies of genetic association using a novel latent-class model. , 2001, American journal of human genetics.

[46]  S. Gabriel,et al.  Assessing the impact of population stratification on genetic association studies , 2004, Nature Genetics.

[47]  J. Miller,et al.  Effect of pre- and postmortem variables on specific mRNA levels in human brain. , 1991, Brain research. Molecular brain research.

[48]  B. Winblad,et al.  The patients dying after long terminal phase have acidotic brains; implications for biochemical measurements on autopsy tissue , 2005, Journal of Neural Transmission.

[49]  Aleksey A. Nakorchevskiy,et al.  Expression deconvolution: A reinterpretation of DNA microarray data reveals dynamic changes in cell populations , 2003, Proceedings of the National Academy of Sciences of the United States of America.