@article {949, title = {Transcriptome profiling of the intoxication response of Tenebrio molitor larvae to Bacillus thuringiensis Cry3Aa protoxin.}, journal = {PloS one}, volume = {7}, year = {2012}, month = {2012}, pages = {e34624}, abstract = {Bacillus thuringiensis (Bt) crystal (Cry) proteins are effective against a select number of insect pests, but improvements are needed to increase efficacy and decrease time to mortality for coleopteran pests. To gain insight into the Bt intoxication process in Coleoptera, we performed RNA-Seq on cDNA generated from the guts of Tenebrio molitor larvae that consumed either a control diet or a diet containing Cry3Aa protoxin. Approximately 134,090 and 124,287 sequence reads from the control and Cry3Aa-treated groups were assembled into 1,318 and 1,140 contigs, respectively. Enrichment analyses indicated that functions associated with mitochondrial respiration, signalling, maintenance of cell structure, membrane integrity, protein recycling/synthesis, and glycosyl hydrolases were significantly increased in Cry3Aa-treated larvae, whereas functions associated with many metabolic processes were reduced, especially glycolysis, tricarboxylic acid cycle, and fatty acid synthesis. Microarray analysis was used to evaluate temporal changes in gene expression after 6, 12 or 24 h of Cry3Aa exposure. Overall, microarray analysis indicated that transcripts related to allergens, chitin-binding proteins, glycosyl hydrolases, and tubulins were induced, and those related to immunity and metabolism were repressed in Cry3Aa-intoxicated larvae. The 24 h microarray data validated most of the RNA-Seq data. Of the three intoxication intervals, larvae demonstrated more differential expression of transcripts after 12 h exposure to Cry3Aa. Gene expression examined by three different methods in control vs. Cry3Aa-treated larvae at the 24 h time point indicated that transcripts encoding proteins with chitin-binding domain 3 were the most differentially expressed in Cry3Aa-intoxicated larvae. Overall, the data suggest that T. molitor larvae mount a complex response to Cry3Aa during the initial 24 h of intoxication. Data from this study represent the largest genetic sequence dataset for T. molitor to date. Furthermore, the methods in this study are useful for comparative analyses in organisms lacking a sequenced genome.}, keywords = {Administration, Animals, Bacterial Proteins, Base Sequence, Biosynthetic Pathways, Complementary, DNA, Endotoxins, Energy Metabolism, Gene Expression Profiling, Hemolysin Proteins, Larva, Microarray Analysis, Molecular Sequence Data, Oral, Sequence Analysis, Tenebrio, Time Factors, Transcriptome}, issn = {1932-6203}, doi = {10.1371/journal.pone.0034624}, author = {Oppert, Brenda and Dowd, Scot E and Bouffard, Pascal and Li, Lewyn and Ana Conesa and Lorenzen, Marc{\'e} D and Toutges, Michelle and Marshall, Jeremy and Huestis, Diana L and Fabrick, Jeff and Oppert, Cris and Jurat-Fuentes, Juan Luis} } @article {580, title = {Functional assessment of time course microarray data.}, journal = {BMC Bioinformatics}, volume = {10 Suppl 6}, year = {2009}, month = {2009 Jun 16}, pages = {S9}, abstract = {

MOTIVATION: Time-course microarray experiments study the progress of gene expression along time across one or several experimental conditions. Most developed analysis methods focus on the clustering or the differential expression analysis of genes and do not integrate functional information. The assessment of the functional aspects of time-course transcriptomics data requires the use of approaches that exploit the activation dynamics of the functional categories to where genes are annotated.

METHODS: We present three novel methodologies for the functional assessment of time-course microarray data. i) maSigFun derives from the maSigPro method, a regression-based strategy to model time-dependent expression patterns and identify genes with differences across series. maSigFun fits a regression model for groups of genes labeled by a functional class and selects those categories which have a significant model. ii) PCA-maSigFun fits a PCA model of each functional class-defined expression matrix to extract orthogonal patterns of expression change, which are then assessed for their fit to a time-dependent regression model. iii) ASCA-functional uses the ASCA model to rank genes according to their correlation to principal time expression patterns and assess functional enrichment on a GSA fashion. We used simulated and experimental datasets to study these novel approaches. Results were compared to alternative methodologies.

RESULTS: Synthetic and experimental data showed that the different methods are able to capture different aspects of the relationship between genes, functions and co-expression that are biologically meaningful. The methods should not be considered as competitive but they provide different insights into the molecular and functional dynamic events taking place within the biological system under study.

}, keywords = {Computer Simulation, Gene Expression Profiling, Oligonucleotide Array Sequence Analysis, Time Factors}, issn = {1471-2105}, doi = {10.1186/1471-2105-10-S6-S9}, author = {Nueda, Maria Jos{\'e} and Sebasti{\'a}n, Patricia and Tarazona, Sonia and Garcia-Garcia, Francisco and Dopazo, Joaquin and Ferrer, Alberto and Conesa, Ana} }