Assessment Method for a Power Analysis to Identify Differentially Expressed Pathways

Gene expression data can provide a very rich source of information for elucidating the biological function on the pathway level if the experimental design considers the needs of the statistical analysis methods. The purpose of this paper is to provide a comparative analysis of statistical methods for detecting the differentially expression of pathways (DEP). In contrast to many other studies conducted so far, we use three novel simulation types, producing a more realistic correlation structure than previous simulation methods. This includes also the generation of surrogate data from two large-scale microarray experiments from prostate cancer and ALL. As a result from our comprehensive analysis of 41,004 parameter configurations, we find that each method should only be applied if certain conditions of the data from a pathway are met. Further, we provide method-specific estimates for the optimal sample size for microarray experiments aiming to identify DEP in order to avoid an underpowered design. Our study highlights the sensitivity of the studied methods on the parameters of the system. © 2012 Tripahti and Emmert-Streib.

Sistemas de informação para DNA microarrays

O projecto de sequenciação do genoma humano veio abrir caminho para o surgimento de novas áreas transdisciplinares de investigação, como a biologia computacional, a bioinformática e a bioestatística. Um dos resultados emergentes desde advento foi a tecnologia de DNA microarrays, que permite o estudo do perfil da expressão de milhares de genes, quando sujeitos a perturbações externas. Apesar de ser uma tecnologia relativamente consolidada, continua a apresentar um conjunto vasto de desafios, nomeadamente do ponto de vista computacional e dos sistemas de informação. São exemplos a optimização dos procedimentos de tratamento de dados bem como o desenvolvimento de metodologias de interpretação semi-automática dos resultados. O principal objectivo deste trabalho consistiu em explorar novas soluções técnicas para agilizar os procedimentos de armazenamento, partilha e análise de dados de experiências de microarrays. Com esta finalidade, realizou-se uma análise de requisitos associados às principais etapas da execução de uma experiência, tendo sido identificados os principais défices, propostas estratégias de melhoramento e apresentadas novas soluções. Ao nível da gestão de dados laboratoriais, é proposto um LIMS (Laboratory Information Management System) que possibilita a gestão de todos os dados gerados e dos procedimentos realizados. Este sistema integra ainda uma solução que permite a partilha de experiências, de forma a promover a participação colaborativa de vários investigadores num mesmo projecto, mesmo usando LIMS distintos. No contexto da análise de dados, é apresentado um modelo que facilita a integração de algoritmos de processamento e de análise de experiências no sistema desenvolvido. Por fim, é proposta uma solução para facilitar a interpretação biológica de um conjunto de genes diferencialmente expressos, através de ferramentas que integram informação existente em diversas bases de dados biomédicas.

Modulation des réactions alloimmunitaires par les cytokines maîtresses IFN-γ et TGF-β

L’injection de cellules immunologiquement compétentes à un hôte histo-incompatible amène une réaction qui peut se traduire par la maladie du greffon-contre-l’hôte (GVHD). La GVHD demeure une barrière importante à une utilisation plus répandue de la greffe allogénique de cellules hématopoïétiques (AHCT), pourtant un traitement efficace pour traiter de nombreuses maladies. Une meilleure compréhension des mécanismes qui sous-tendent cette pathologie pourrait en faciliter le traitement et la prévention. L’Interféron-gamma (IFN-γ) et le Transforming Growth Factor-béta (TGF-β) sont deux cytokines maîtresses de l’immunité impliquées dans la fonction et l’homéostasie des cellules greffées. Nous démontrons chez la souris que l’IFN-γ limite la reconstitution lympho-hématopoïétique de façon dose-dépendante en mobilisant des mécanismes d’apoptose et en inhibant la prolifération cellulaire. Le TGF-β est quant à lui généralement connu comme un immunosuppresseur qui contrôle l’immunité en utilisant plusieurs voies de signalisation. Le rôle relatif de ces voies en AHCT est inconnu. Nous avons étudié une de ces voies en greffant des cellules provenant de donneurs déficients pour le gène SMAD3 (SMAD3-KO), un médiateur central de la voie canonique du TGF-β, à des souris histo-incompatibles. Bien que l’absence de SMAD3 ne cause aucune maladie chez nos souris donneuses, l’injection de cellules SMAD3-KO amène une GVHD du colon sévère chez le receveur. Cette atteinte est caractérisée par une différenciation Th1 et une infiltration massive de granulocytes témoignant d’un rôle central de SMAD3 dans la physiologie des lymphocytes T CD4 et des cellules myéloïdes. Nous avons focalisé ensuite nos efforts sur le rôle de SMAD3 chez les lymphocytes T CD4 en sachant que SMAD3 était actif chez les lymphocytes T CD4 tolérants. Nous avons découvert que SMAD3 était rapidement inactivé après une activation des cellules T, suggérant que l’inactivation de SMAD3 était fonctionnellement importante pour briser l’état de tolérance. Des études de micro-puces d’ADNc nous ont montré que SMAD3 contrôlait en effet l’expression de nombreux transcrits de gènes connus comme étant reliés à la tolérance et/ou à des processus biologiques dont les rôles dans le maintien de la tolérance sont plausibles.

Comparaison des méthodes d'analyse de l'expression différentielle basée sur la dépendance des niveaux d'expression

La technologie des microarrays demeure à ce jour un outil important pour la mesure de l'expression génique. Au-delà de la technologie elle-même, l'analyse des données provenant des microarrays constitue un problème statistique complexe, ce qui explique la myriade de méthodes proposées pour le pré-traitement et en particulier, l'analyse de l'expression différentielle. Toutefois, l'absence de données de calibration ou de méthodologie de comparaison appropriée a empêché l'émergence d'un consensus quant aux méthodes d'analyse optimales. En conséquence, la décision de l'analyste de choisir telle méthode plutôt qu'une autre se fera la plupart du temps de façon subjective, en se basant par exemple sur la facilité d'utilisation, l'accès au logiciel ou la popularité. Ce mémoire présente une approche nouvelle au problème de la comparaison des méthodes d'analyse de l'expression différentielle. Plus de 800 pipelines d'analyse sont appliqués à plus d'une centaine d'expériences sur deux plateformes Affymetrix différentes. La performance de chacun des pipelines est évaluée en calculant le niveau moyen de co-régulation par l'entremise de scores d'enrichissements pour différentes collections de signatures moléculaires. L'approche comparative proposée repose donc sur un ensemble varié de données biologiques pertinentes, ne confond pas la reproductibilité avec l'exactitude et peut facilement être appliquée à de nouvelles méthodes. Parmi les méthodes testées, la supériorité de la sommarisation FARMS et de la statistique de l'expression différentielle TREAT est sans équivoque. De plus, les résultats obtenus quant à la statistique d'expression différentielle corroborent les conclusions d'autres études récentes à propos de l'importance de prendre en compte la grandeur du changement en plus de sa significativité statistique.

Histone H2B-R95A mutant identifies the pheromone pathway that signals cell cycle arrest during rapamycin response

La rapamycine est un immunosuppresseur utilisé pour traiter plusieurs types de maladies dont le cancer du rein. Son fonctionnement par l’inhibition de la voie de Tor mène à des changements dans des processus physiologiques, incluant le cycle cellulaire. Chez Saccharomyces cerevisiae, la rapamycine conduit à une altération rapide et globale de l’expression génique, déclenchant un remodelage de la chromatine. Nous proposons que les modifications des histones peuvent jouer un rôle crucial dans le remodelage de la chromatine en réponse à la rapamycine. Notre objectif principal est d’identifier d’une banque de mutants d’histone les variantes qui vont échouer à répondre à la rapamycine dans une tentative de réaliser une caractérisation des modifications d’histone critiques pour la réponse à cette drogue. Ainsi, nous avons réalisé un criblage d’une banque de mutants d’histone et identifié plusieurs mutants d‘histone dont la résistance à la rapamycine a été altérée. Nous avons caractérisé une de ces variantes d’histone, à savoir H2B, qui porte une substitution de l’alanine en arginine en position 95 (H2B-R95A) et démontré que ce mutant est extrêmement résistant à la rapamycine, et non à d’autres drogues. Des immunoprécipitations ont démontré que H2B-R95A est défectueux pour former un complexe avec Spt16, un facteur essentiel pour la dissociation de H2A et H2B de la chromatine, permetant la réplication et la transcription par les ADN et ARN polymérases, respectivement. Des expériences de ChIP-Chip et de micropuce ont démontré que l’arginine 95 de H2B est requise pour recruter Spt16 afin de permettre l’expression d’une multitude de gènes, dont certains font partie de la voie des phéromones. Des évidences seront présentées pour la première fois démontrant que la rapamycine peut activer la voie des phéromones et qu’une défectuosité dans cette voie cause la résistante à cette drogue.

Pea-barley intercropping and short-term subsequent crop effects across European organic cropping conditions

Grain legumes are known to increase the soil mineral nitrogen (N) content, reduce the infection pressure of soil borne pathogens, and hence enhance subsequent cereals yields. Replicated field experiments were performed throughout W. Europe (Denmark, United Kingdom, France, Germany and Italy) to asses the effect of intercropping pea and barley on the N supply to subsequent wheat in organic cropping systems. Pea and barley were grown either as sole crops at the recommended plant density (P100 and B100, respectively) or in replacement (P50B50) or additive (P100B50) intercropping designs. In the replacement design the total relative plant density is kept constant, while the additive design uses the optimal sole crop density for pea supplementing with 'extra' barley plants. The pea and barley crops were followed by winter wheat with and without N application. Additional experiments in Denmark and the United Kingdom included subsequent spring wheat with grass-clover as catch crops. The experiment was repeated over the three cropping seasons of 2003, 2004 and 2005. Irrespective of sites and intercrop design pea-barley intercropping improved the plant resource utilization (water, light, nutrients) to grain N yield with 25-30% using the Land Equivalent ratio. In terms of absolute quantities, sole cropped pea accumulated more N in the grains as compared to the additive design followed by the replacement design and then sole cropped barley. The post harvest soil mineral N content was unaffected by the preceding crops. Under the following winter wheat, the lowest mineral N content was generally found in early spring. Variation in soil mineral N content under the winter wheat between sites and seasons indicated a greater influence of regional climatic conditions and long-term cropping history than annual preceding crop and residue quality. Just as with the soil mineral N, the subsequent crop response to preceding crop was negligible. Soil N balances showed general negative values in the 2-year period, indicating depletion of N independent of preceding crop and cropping strategy. It is recommended to develop more rotational approaches to determine subsequent crop effects in organic cropping systems, since preceding crop effects, especially when including legumes, can occur over several years of cropping.

Modelling simultaneously water content and dry matter dynamics of wheat grains

A model was devised to describe simultaneously the grain masses of water and dry matter against thermal time during grain filling and maturation of winter wheat. The model accounted for a linear increase in water mass of duration anthesis-m(1) (end of rapid water assimilation phase) and rate a, followed by a more stable water mass until in,, after which water mass declined rapidly at rate e. Grain dry matter was described as a linear increase of rate bgf until a maximum size (maxgf) was attained at m(2).The model was fitted to plot data from weekly samples of grains taken from replicated field experiments investigating effects of grain position (apical or medial), fungicide (five contrasting treatments), sowing date (early or late), cultivar (Malacca or Shamrock) and season (2001/2002 and 2002/2003) on grain filling. The model accounted for between 83 and 99% of the variation ( 2) when fitted to data from individual plots, and between 97 and 99% when fitted to treatment means. Endosperm cell number of grains from early-sown plots in the first season were also counted. Differences in maxgf between grain positions and also between cultivars were mostly the result of effects on bgf and were empirically associated with water mass at nil. Fungicide application controlled S. tritici and powdery mildew infection, delayed flag leaf senescence, increased water mass at m(1) (wm(1)), and also increased m(2), bgf and maxgf. Fungicide effects on water mass were detected before fungicide effects on dry matter, but comparison of the effects of individual fungicide treatments showed no evidence that effects on wm(1), nor on endosperm cell numbers at about m(1), were required for fungicide effects on maxgf, (c) 2005 Elsevier B.V. All rights reserved.

Transcriptomic analysis of rhizobium leguminosarum biovar viciae in symbiosis with host plants pisum sativum and Vicia cracca

Rhizobium leguminosarum bv. viciae forms nitrogen-fixing nodules on several legumes, including pea (Pisum sativum) and vetch (Vicia cracca), and has been widely used as a model to study nodule biochemistry. To understand the complex biochemical and developmental changes undergone by R. leguminosarum bv. viciae during bacteroid development, microarray experiments were first performed with cultured bacteria grown on a variety of carbon substrates (glucose, pyruvate, succinate, inositol, acetate, and acetoacetate) and then compared to bacteroids. Bacteroid metabolism is essentially that of dicarboxylate-grown cells (i.e., induction of dicarboxylate transport, gluconeogenesis and alanine synthesis, and repression of sugar utilization). The decarboxylating arm of the tricarboxylic acid cycle is highly induced, as is gamma-aminobutyrate metabolism, particularly in bacteroids from early (7-day) nodules. To investigate bacteroid development, gene expression in bacteroids was analyzed at 7, 15, and 21 days postinoculation of peas. This revealed that bacterial rRNA isolated from pea, but not vetch, is extensively processed in mature bacteroids. In early development (7 days), there were large changes in the expression of regulators, exported and cell surface molecules, multidrug exporters, and heat and cold shock proteins. fix genes were induced early but continued to increase in mature bacteroids, while nif genes were induced strongly in older bacteroids. Mutation of 37 genes that were strongly upregulated in mature bacteroids revealed that none were essential for nitrogen fixation. However, screening of 3,072 mini-Tn5 mutants on peas revealed previously uncharacterized genes essential for nitrogen fixation. These encoded a potential magnesium transporter, an AAA domain protein, and proteins involved in cytochrome synthesis.

Population growth rate and carrying capacity for springtails Folsomia candida exposed to ivermectin

Forecasting the effects of stressors on the dynamics of natural populations requires assessment of the joint effects of a stressor and population density on the population response. The effects can be depicted as a contour map in which the population response, here assessed by Population growth rate, varies with stress and density in the same way that the height of land above sea level varies with latitude and longitude. We present the first complete map of this type using as our model Folsomia candida exposed to five different concentrations of the widespread anthelmintic veterinary medicine ivermectin in replicated microcosm experiments lasting 49 days. The concentrations of ivermectin in yeast were 0.0, 6.8 28.83 66.4 and 210.0 mg/L wet weight. Increasing density and chemical concentration both significantly reduced the population growth rate of Folsomia candida, in part through effects on food consumption and fecundity. The interaction between density and ivermectin concentration was "less-than-additive," implying that at high density populations were able to compensate for the effects of the chemical. This result demonstrates that regulatory protocols carried out at low density (as in most past experiments) may seriously overestimate effects in the field, where densities are locally high and populations are resource limited (e.g., in feces of livestock treated with ivermectin).

Exposure of Escherichia coli and Salmonella enterica serovar Typhimurium to triclosan induces a species-specific response, including drug detoxification

Objectives: The use of triclosan within various environments has been linked to the development of multiple drug resistance (MDR) through the increased expression of efflux pumps such as AcrAB-ToIC. In this work, we investigate the effect of triclosan exposure in order to ascertain the response of two species to the presence of this widely used biocide. Methods: The transcriptomes of Salmonella enterica serovar Typhimurium SL1344 and Escherichia coli K-12 MG1655 after exposure to the MIC of triclosan (0.12 mg/L) were determined in microarray experiments. Phenotypic validation of the transcriptomic data included RT-PCR, ability to form a biofilm and motility assays. Results: Despite important differences in the triclosan-dependent transcriptomes of the two species, increased expression of efflux pump component genes was seen in both. Increased expression of soxS was observed in Salmonella Typhimurium, however, within E. coli, decreased expression was seen. Expression of fabBAGI in Salmonella Typhimurium was decreased, whereas in E. coli expression of fabABFH was increased. Increased expression of ompR and genes within this regulon (e.g. ompC, csgD and ssrA) was seen in the transcriptome of Salmonella Typhimurium. An unexpected response of E. coli was the differential expression of genes within operons involved in iron homeostasis; these included fhu, fep and ent. Conclusions: These data indicate that whilst a core response to triclosan exposure exists, the differential transcriptome of each species was different. This suggests that E. coli K-12 should not be considered the paradigm for the Enterobacteriaceae when exploring the effects of antimicrobial agents.

Nutrient sensing and signalling in plants: Potassium and phosphorus

Potassium and phosphorus are important macronutrients for crops but are often deficient in the field. Very little is known about how plants sense fluctuations in K and P and how information about K and P availability is integrated at the whole plant level into physiological and metabolic adaptations. This chapter reviews recent advances in discovering molecular responses of plants to K and P deficiency by microarray experiments. These studies provide us not only with a comprehensive picture of adaptive mechanisms, but also with a large number of transcriptional markers that can be used to identify upstream components of K and P signalling pathways. On the basis of the available information we discuss putative receptors and signals involved in the sensing and integration of K and P status both at the cellular and at the whole plant level. These involve membrane potential, voltage-dependent ion channels, intracellular Ca and pH, and transcription factors, as well as hormones and metabolites for systemic signalling. Genetic screens of reporter lines for transcriptional markers and metabolome analysis of K- and P-deficient plants are likely to further advance our knowledge in this area in the near future.

Shoot yield drives phosphorus use efficiency in Brassica oleracea and correlates with root architecture traits

The environmental and financial costs of using inorganic phosphate fertilizers to maintain crop yield and quality are high. Breeding crops that acquire and use phosphorus (P) more efficiently could reduce these costs. The variation in shoot P concentration (shoot-P) and various measures of P use efficiency (PUE) were quantified among 355 Brassica oleracea L. accessions, 74 current commercial cultivars, and 90 doubled haploid (DH) mapping lines from a reference genetic mapping population. Accessions were grown at two or more external P concentrations in glasshouse experiments; commercial and DH accessions were also grown in replicated field experiments. Within the substantial species-wide diversity observed for shoot-P and various measures of PUE in B. oleracea, current commercial cultivars have greater PUE than would be expected by chance. This may be a consequence of breeding for increased yield, which is a significant component of most measures of PUE, or early establishment. Root development and architecture correlate with PUE; in particular, lateral root number, length, and growth rate. Significant quantitative trait loci associated with shoot-P and PUE occur on chromosomes C3 and C7. These data provide information to initiate breeding programmes to improve PUE in B. oleracea.

The global effect of follicle-stimulating hormone and tumour necrosis factor α on gene expression in cultured bovine ovarian granulosa cells

Background Oocytes mature in ovarian follicles surrounded by granulosa cells. During follicle growth, granulosa cells replicate and secrete hormones, particularly steroids close to ovulation. However, most follicles cease growing and undergo atresia or regression instead of ovulating. To investigate the effects of stimulatory (follicle-stimulating hormone; FSH) and inhibitory (tumour necrosis factor alpha; TNFα) factors on the granulosa cell transcriptome, bovine ovaries were obtained from a local abattoir and pools of granulosa cells were cultured in vitro for six days under defined serum-free conditions with treatments present on days 3–6. Initially dose–response experiments (n = 4) were performed to determine the optimal concentrations of FSH (0.33 ng/ml) and TNFα (10 ng/ml) to be used for the microarray experiments. For array experiments cells were cultured under control conditions, with FSH, with TNFα, or with FSH plus TNFα (n = 4 per group) and RNA was harvested for microarray analyses. Results Statistical analysis showed primary clustering of the arrays into two groups, control/FSH and TNFα/TNFα plus FSH. The effect of TNFα on gene expression dominated that of FSH, with substantially more genes differentially regulated, and the pathways and genes regulated by TNFα being similar to those of FSH plus TNFα treatment. TNFα treatment reduced the endocrine activity of granulosa cells with reductions in expression of FST, INHA, INBA and AMH. The top-ranked canonical pathways and GO biological terms for the TNFα treatments included antigen presentation, inflammatory response and other pathways indicative of innate immune function and fibrosis. The two most significant networks also reflect this, containing molecules which are present in the canonical pathways of hepatic fibrosis/hepatic stellate cell activation and transforming growth factor β signalling, and these were up regulated. Upstream regulator analyses also predicted TNF, interferons γ and β1 and interleukin 1β. Conclusions In vitro, the transcriptome of granulosa cells responded minimally to FSH compared with the response to TNFα. The response to TNFα indicated an active process akin to tissue remodelling as would occur upon atresia. Additionally there was reduction in endocrine function and induction of an inflammatory response to TNFα that displays features similar to immune cells.

Protein variation in blood-dwelling schistosome worms generated by differential splicing of micro-exon gene transcripts

Schistosoma mansoni is a well-adapted blood-dwelling parasitic helminth, persisting for decades in its human host despite being continually exposed to potential immune attack. Here, we describe in detail micro-exon genes (MEG) in S. mansoni, some present in multiple copies, which represent a novel molecular system for creating protein variation through the alternate splicing of short (<= 36 bp) symmetric exons organized in tandem. Analysis of three closely related copies of one MEG family allowed us to trace several evolutionary events and propose a mechanism for micro-exon generation and diversification. Microarray experiments show that the majority of MEGs are up-regulated in life cycle stages associated with establishment in the mammalian host after skin penetration. Sequencing of RT-PCR products allowed the description of several alternate splice forms of micro-exon genes, highlighting the potential use of these transcripts to generate a complex pool of protein variants. We obtained direct evidence for the existence of such pools by proteomic analysis of secretions from migrating schistosomula and mature eggs. Whole-mount in situ hybridization and immunolocalization showed that MEG transcripts and proteins were restricted to glands or epithelia exposed to the external environment. The ability of schistosomes to produce a complex pool of variant proteins aligns them with the other major groups of blood parasites, but using a completely different mechanism. We believe that our data open a new chapter in the study of immune evasion by schistosomes, and their ability to generate variant proteins could represent a significant obstacle to vaccine development.

Global Gene Expression Analysis during Sporulation of the Aquatic Fungus Blastocladiella emersonii

The Blastocladiella emersonii life cycle presents a number of drastic biochemical and morphological changes, mainly during two cell differentiation stages: germination and sporulation. To investigate the transcriptional changes taking place during the sporulation phase, which culminates with the production of the zoospores, motile cells responsible for the dispersal of the fungus, microarray experiments were performed. Among the 3,773 distinct genes investigated, a total of 1,207 were classified as differentially expressed, relative to time zero of sporulation, at at least one of the time points analyzed. These results indicate that accurate transcriptional control takes place during sporulation, as well as indicating the necessity for distinct molecular functions throughout this differentiation process. The main functional categories overrepresented among upregulated genes were those involving the microtubule, the cytoskeleton, signal transduction involving Ca(2+), and chromosome organization. On the other hand, protein biosynthesis, central carbon metabolism, and protein degradation were the most represented functional categories among downregulated genes. Gene expression changes were also analyzed in cells sporulating in the presence of subinhibitory concentrations of glucose or tryptophan. Data obtained revealed overexpression of microtubule and cytoskeleton transcripts in the presence of glucose, probably causing the shape and motility problems observed in the zoospores produced under this condition. In contrast, the presence of tryptophan during sporulation led to upregulation of genes involved in oxidative stress, proteolysis, and protein folding. These results indicate that distinct physiological pathways are involved in the inhibition of sporulation due to these two classes of nutrient sources.