The Library

Statistical modelling of transcript profiles of differentially regulated genes

Tools

Eastwood, Daniel C., Mead, A. (Andrew), Sergeant, Martin J. and Burton, Kerry S.. (2008) Statistical modelling of transcript profiles of differentially regulated genes. BMC Molecular Biology, Vol.9 (No.66). ISSN 1471-2199

Preview

PDF
WRAP_Mead_1471-2199-9-66.pdf - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader
Download (750Kb)

Official URL: http://dx.doi.org/10.1186/1471-2199-9-66

Abstract

Background: The vast quantities of gene expression profiling data produced in microarray studies, and the more precise quantitative PCR, are often not statistically analysed to their full potential. Previous studies have summarised gene expression profiles using simple descriptive statistics, basic analysis of variance (ANOVA) and the clustering of genes based on simple models fitted to their expression profiles over time. We report the novel application of statistical non-linear regression modelling techniques to describe the shapes of expression profiles for the fungus Agaricus bisporus, quantified by PCR, and for E. coli and Rattus norvegicus, using microarray technology. The use of parametric non-linear regression models provides a more precise description of expression profiles, reducing the "noise" of the raw data to produce a clear "signal" given by the fitted curve, and describing each profile with a small number of biologically interpretable parameters. This approach then allows the direct comparison and clustering of the shapes of response patterns between genes and potentially enables a greater exploration and interpretation of the biological processes driving gene expression. Results: Quantitative reverse transcriptase PCR-derived time-course data of genes were modelled. "Splitline" or "broken-stick" regression identified the initial time of gene up-regulation, enabling the classification of genes into those with primary and secondary responses. Five-day profiles were modelled using the biologically-oriented, critical exponential curve, y(t) = A + (B + Ct)Rt + ε. This non-linear regression approach allowed the expression patterns for different genes to be compared in terms of curve shape, time of maximal transcript level and the decline and asymptotic response levels. Three distinct regulatory patterns were identified for the five genes studied. Applying the regression modelling approach to microarray-derived time course data allowed 11% of the Escherichia coli features to be fitted by an exponential function, and 25% of the Rattus norvegicus features could be described by the critical exponential model, all with statistical significance of p < 0.05. Conclusion: The statistical non-linear regression approaches presented in this study provide detailed biologically oriented descriptions of individual gene expression profiles, using biologically variable data to generate a set of defining parameters. These approaches have application to the modelling and greater interpretation of profiles obtained across a wide range of platforms, such as microarrays. Through careful choice of appropriate model forms, such statistical regression approaches allow an improved comparison of gene expression profiles, and may provide an approach for the greater understanding of common regulatory mechanisms between genes.

Item Type:	Journal Article
Subjects:	Q Science > QH Natural history > QH426 Genetics
Divisions:	Faculty of Science > Life Sciences (2010- ) > Warwick HRI (2004-2010)
Library of Congress Subject Headings (LCSH):	Gene mapping, Regression analysis
Journal or Publication Title:	BMC Molecular Biology
Publisher:	Biomed central
ISSN:	1471-2199
Date:	23 July 2008
Volume:	Vol.9
Number:	No.66
Identification Number:	10.1186/1471-2199-9-66
Status:	Peer Reviewed
Access rights to Published version:	Open Access
Description:	Final version (published as open access).
References:	1. Dopazo J, Zanders E, Dragoni I, Amphlett G, Falciani F: Methods and approaches in the analysis of gene expression data. J Immunol Methods 2001, 250(1–2):93-112. 2. Tamames J, Clark D, Herrero J, Dopazo J, Blaschke C, Fernandez JM, Oliveros JC, Valencia A: Bioinformatics methods for the analysis of expression arrays: data clustering and information extraction. J Biotechnol 2002, 98:269-283. 3. Reimers M: Statistical analysis of microarray data. Addict Biol 2005, 10:23-35. 4. Tai YC, Speed TP: A multivariate empirical Bayes statistic for replicated microarray time course data. Ann Stat 2006, 34(6):2387-2412. 5. Manfield IW, Jen CH, Pinney JW, Michalopoulos I, Bradford JR, Gilmartin PM, Westhead DR: Arabidopsis Co-expression Tool (ACT): web server tools for microarray-based gene expression analysis. Nucleic Acids Res 2006, 34:504-509. 6. Persson S, Wei H, Milne J, Page GR, Somerville CR: Identification of genes required for cellulose synthesis by regression analysis of public microarray data sets. Proc Natl Acad Sci USA 2005, 102(24):8633-8638. 7. Conesa A, Neuda MJ, Ferrer A, Talön M: maSigPro: A method to identify significantly differential expression profiles in timecourse microarray experiments. Bioinformatics 2006, 22(9):1096-1102. 8. Heard NA, Holmes CC, Stephens DA: A quantitative study of gene regulation involved in the immune response of Anopheline mosquitoes: An application of Bayesian hierarchical clustering of curves. J Am Stat Assoc 2006, 101(473):18-29. 9. Kahnin R, Vinciotti V, Mersinias V, Smith CP, Wit P: Statistical reconstruction of transcription factor activity using Michaelis- Menten kinetics. Biometrics 2007, 63(3):816-823. 10. Nachman I, Regev A, Friedman N: Inferring quantitative models of regulatory networks from expression data. Bioinformatics 2004, 20(suppl 1):i248-i256. 11. Liss B: Improved quantitative real-time RT-PCR for expression profiling of individual cells. Nucleic Acids Res 2002, 30(17):89-98. 12. Eastwood DC, Kingsnorth CS, Jones HJ, Burton KS: Genes with increased transcript levels following harvest of the sporophore of Agaricus bisporus have multiple physiological roles. Mycol Res 2001, 105(10):1223-123. 13. Kingsnorth CS, Eastwood DC, Burton KS: Cloning and post-harvest expression of serine proteinase transcripts in the cultivated mushroom Agaricus bisporus. Fungal Gen Biol 2001, 32:135-144. 14. Wagemaker MJM, Eastwood DC, Wellboren W, Burton KS, Drift C Van Der, Jetten MSM, Van Griensven LJLD, Op Den Camp HJM: Argininosuccinate synthetase and argininosuccinate lyase: two ornithine cycle enzymes from Agaricus bisporus. Mycol Res 2007, 111:493-502. 15. Miyazaki Y, Nakamura M, Babasaki K: Molecular cloning of developmentally specific genes by representational difference analysis during fruiting body formation in the basidiomycete Lentinula edodes. Fungal Gen Biol 2005, 42:493-505. 16. Lee S-H, Kim B-G, Kim K-J, Lee J-S, Yun D-W, Hahn J-H, Kim G-H, Lee K-H, Suh D-S, Kwon S-T, Lee C-S, Yoo Y-B: Comparative analysis of sequences expressed during liquid-cultured mycelia and fruit body stages of Pleurotus ostreatus. Fungal Gen Biol 2002, 35:115-134. 17. Yamada M, Sakuraba S, Shibata K, Taguchi G, Inatomi S, Okazaki M, Shimosaka M: Isolation and analysis of genes specifically expressed during fruiting body development in the basidiomycete Flammulina velutipes by fluorescence differential display. FEMS Microbiol Lett 2006, 254:165-172. 18. Kues U: Life history and developmental processes in the basidiomycete Coprinus cinereus. Microbiol Mol Biol Rev 2000, 64(2):316-353. 19. Kamada T: Molecular genetics of sexual development in the mushroom Coprinus cinereus. BioEssays 2002, 24:449-459. 20. Sreenivasaprasad S, Eastwood DC, Browning N, Lewis SMJ, Burton KS: Differential expression of a putative riboflavin-aldehydeforming enzyme (raf) gene during development and postharvest storage and in different tissue of the sporophore in Agaricus bisporus. Appl Microbiol Biotechnol 2006, 70:470-476. 21. Blanchard JL, Wholey W-Y, Conlon EM, Pomposiello PJ: Rapid changes in gene expression dynamics in response to superoxide reveal SoxRS-dependent and independent transcriptional networks. PLoS ONE 2007, 2(11):e1186. 22. Jin JY, Almon RR, Dubois DC, Jusko J: Modelling of corticosteroid pharmacogenomics in rat liver using gene microarrays. J Pharmacol Exp Ther 2003, 307:93-109. 23. Bascampte J, Rodriguez MA: Habitat patchiness and plant species richness. Ecol Lett 2001, 4:417-420. 24. Carroll JE, Wilcox WF: Effects of humidity on the development of grapevine powdery mildew. Phytopathology 2003, 93:1137-1144. 25. Price RI, Walters MJ, Retallack RW, Henderson NK, Kerr D, Henzell S, Dhaliwal S, Prince RL: Impact of the analysis of a bone density reference range on determination of the T-score. J Clin Densitom 2003, 6(1):51-62. 26. Ikner A, Shiozaki K: Yeast signalling pathways in oxidative stress response. Mutat Res-Fund Mol Mech Mutagen 2005, 569(1– 2):13-27. 27. Umar MH, van Griensven LJLD: Morphological studies on the life span, development stages, senescence and death of fruiting bodies of Agaricus bisporus. Mycol Res 1997, 101:1409-1422. 28. Braaksma A, van Doorn AA, Kieft H, van Aelist AC: Morphometric analysis of ageing mushrooms (Agaricus bisporus) during post-harvest development. Postharvest Biol Technol 1998, 13:71-79. 29. Ge H, Roeder RG: The high-mobility group protein HMG1 can reversibly inhibit class-II gene-transcription by interaction with the TATA-binding protein. J Biol Chem 1994, 269(25):17136-17140. 30. Stros M, Ozaki T, Bacikova A, Kageyama H, Nakagawara A: HMGB1 and HMGB2 cell-specifically down-regulate the p53-and p- 73-dependant sequence-specific transactivation from the human Bax gene promoter. J Biol Chem 2002, 277(9):7157-7164. 31. Hammond JBW, Nichols R: Changes in respiration and carbohydrates during the post-harvest storage of mushrooms (Agaricus bisporus). J Sci Food Agric 1975, 26:835-842. 32. Stekel D: Microarray Bioinformatics Cambridge University Press, Cambridge, UK; 2003. 33. Rangel C, Angus J, Ghahramani Z, Lioumi M, Sotheran E, Gaiba A, Wild DL, Falciani F: Modelling T-cell activation using gene expression profiling and state-space models. Bioinformatics 2004, 20(9):1361-1372. 34. Beal MJ, Falciani F, Ghahramani Z, Rangel C, Wild DL: A Bayesian approach to reconstructing genetic regulatory networks with hiddenfactors. Bioinformatics 2005, 21(3):349-356. 35. Sambrook J, Russell DW: Molecular cloning: a laboratory manual 3rd edition. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY; 2001. 36. Rosen KM, Villa-Komaroff L: An alternative method for the visualization of RNA in formaldehyde agarose gels. Focus 1993, 12(2):23-24. 37. Goidin D, Mamessier A, Statquet M.-J, Schmitt D, Berthier-Vergnes O: Ribosomal 18S RNA prevails over Glyceraldehyde-3- phosphate dehydrogenase and β-actin genes as internal standards for quantitative comparison of mRNA levels in invasive and noninvasive human melanoma cell subpopulations. Anal Biochem 2000, 295(1):17-21. 38. Lekanne Deprez RH, Fijnvandraat AC, Ruijter JM, Moorman AFM: Sensitivity and accuracy of quantitative real-time polymerase chain reaction using SYBR green I depends on cDNA synthesis conditions. Anal Biochem 2002, 307(1):63-69. 39. de Groot PWJ, Schaap PJ, van Griensven LJLD, Visser J: Isolation of developmentally regulated genes from the edible mushroom Agaricus bisporus. Microbiology 1997, 143:1993-2001.
URI:	http://wrap.warwick.ac.uk/id/eprint/204