Genome scale reconstruction of a Salmonella metabolic model comparison of similarity and differences with a commensal Escherichia coli strain

Salmonella are closely related to commensal Escherichia coli but have gained virulence factors enabling them to behave as enteric pathogens. Less well studied are the similarities and differences that exist between the metabolic properties of these organisms that may contribute toward niche adaptation of Salmonella pathogens. To address this, we have constructed a genome scale Salmonella metabolic model (iMA945). The model comprises 945 open reading frames or genes, 1964 reactions, and 1036 metabolites. There was significant overlap with genes present in E. coli MG1655 model iAF1260. In silico growth predictions were simulated using the model on different carbon, nitrogen, phosphorous, and sulfur sources. These were compared with substrate utilization data gathered from high throughput phenotyping microarrays revealing good agreement. Of the compounds tested, the majority were utilizable by both Salmonella and E. coli. Nevertheless a number of differences were identified both between Salmonella and E. coli and also within the Salmonella strains included. These differences provide valuable insight into differences between a commensal and a closely related pathogen and within different pathogenic strains opening new avenues for future explorations.

Proteomic analyses of a Listeria monocytogenes mutant lacking σB identify new components of the σB regulon and highlight a role for σB in the utilization of glycerol

In Listeria monocytogenes the alternative sigma factor σB plays important roles in both virulence and stress tolerance. In this study a proteomic approach was used to define components of the σB regulon in L. monocytogenes 10403S (serotype 1/2a). Using two-dimensional gel electrophoresis and the recently developed isobaric tags for relative and absolute quantitation technique, the protein expression profiles of the wild type and an isogenic ΔsigB deletion strain were compared. Overall, this study identified 38 proteins whose expression was σB dependent; 17 of these proteins were found to require the presence of σB for full expression, while 21 were expressed at a higher level in the ΔsigB mutant background. The data obtained with the two proteomic approaches showed limited overlap (four proteins were identified by both methods), a finding that highlights the complementarity of the two technologies. Overall, the proteomic data reaffirmed a role for σB in the general stress response and highlighted a probable role for σB in metabolism, especially in the utilization of alternative carbon sources. Proteomic and physiological data revealed the involvement of σB in glycerol metabolism. Five newly identified members of the σB regulon were shown to be under direct regulation of σB using reverse transcription-PCR (RT-PCR), while random amplification of cDNA ends-PCR was used to map four σB-dependent promoters upstream from lmo0796, lmo1830, lmo2391, and lmo2695. Using RT-PCR analysis of known and newly identified σB-dependent genes, as well as proteomic analyses, σB was shown to play a major role in the stationary phase of growth in complex media.

Genistein protects prostate cells against hydrogen peroxide-induced DNA damage and induces expression of genes involved in the defence against oxidative stress

Prostate cancer is one of the most frequent cancer types in Western societies and predominately occurs in the elderly male. The strong age-related increase of prostate cancer is associated with a progressive accumulation of oxidative DNA damage which is presumably supported by a decline of the cellular antioxidative defence during ageing. Risk of developing prostate cancer is much lower in many Asian countries where soy food is an integral part of diet. Therefore, isoflavones from soy were suggested to have chemopreventive activities in prostate cells. Here, we have investigated the hypothesis that the soy-isoflavone genistein could protect DNA of LAPC-4 prostate cells from oxidative stress-related damage by enhancing the expression of antioxidative genes and proteins. A 24 h preincubation with genistein (1-30 microM) protected cells from hydrogen peroxide-induced DNA damage, as determined by the comet assay. Analysis of two cDNA macroarrays, each containing 96 genes of biotransformation and stress response, revealed a modulated expression of 3 genes at 1 microM and of 19 genes at 10 microM genistein. Real-time PCR confirmed the induction of three genes encoding products with antioxidant activities, namely glutathione reductase (2.7-fold), microsomal glutathione S-transferase 1 (1.9-fold) and metallothionein 1X (6.3-fold), at 1-30 microM genistein. 17Beta-estradiol, in contrast, decreased the expression of metallothionein 1X at 0.3 microM (2.0-fold), possibly pointing to an estrogen receptor-mediated regulation of this gene. Immunocytochemical staining revealed an induction of metallothionein proteins at 30 microM genistein, while their intracellular localization was unaffected. Metallothioneins were previously found to protect cells from hydrogen peroxide-induced DNA damage. Hence, our findings indicate that genistein protects prostate cells from oxidative stress-related DNA damage presumably by inducing the expression of antioxidative products, such as metallothioneins. Genistein, therefore, might counteract the age-related decline of important antioxidative defence systems which in turn maintain DNA integrity.

Cloning, expression, and characterization of human cytosolic aminopeptidase P: a single manganese(II)-dependent enzyme

The mammalian bradykinin-degrading enzyme aminopeptidase P (AP-P; E. C. 3.4.11.9) is a metal-dependent enzyme and is a member of the peptidase clan MG. AP-P exists as membrane-bound and cytosolic forms, which represent distinct gene products. A partially truncated clone encoding the cytosolic form was obtained from a human pancreatic cDNA library and the 5' region containing the initiating Met was obtained by 5' rapid accumulation of cDNA ends (RACE). The open reading frame encodes a protein of 623 amino acids with a calculated molecular mass of 69,886 Da. The full-length cDNA with a C-terminal hexahistidine tag was expressed in Escherichia coli and COS-1 cells and migrated on SDS-PAGE with a molecular mass of 71 kDa. The expressed cytosolic AP-P hydrolyzed the X-Pro bond of bradykinin and substance P but did not hydrolyze Gly-Pro-hydroxyPro. Hydrolysis of bradykinin was inhibited by 1,10-phenanthroline and by the specific inhibitor of the membrane-bound form of mammalian AP-P, apstatin. Inductively coupled plasma atomic emission spectroscopy of AP-P expressed in E. coli revealed the presence of 1 mol of manganese/mol of protein and insignificant amounts of cobalt, iron, and zinc. The enzymatic activity of AP-P was promoted in the presence of Mn(II), and this activation was increased further by the addition of glutathione. The only other metal ion to cause slight activation of the enzyme was Co(II), with Ca(II), Cu(II), Mg(II), Ni(II), and Zn(II) all being inhibitory. Removal of the metal ion from the protein was achieved by treatment with 1,10-phenanthroline. The metal-free enzyme was reactivated by the addition of Mn(II) and, partially, by Fe(II). Neither Co(II) nor Zn(II) reactivated the metal-free enzyme. On the basis of these data we propose that human cytosolic AP-P is a single metal ion-dependent enzyme and that manganese is most likely the metal ion used in vivo.

The early transcriptomic response to interleukin 1β and interleukin 33 in rat neonatal cardiomyocytes

In the heart, inflammatory cytokines including interleukin (IL) 1β are implicated in regulating adaptive and maladaptive changes, whereas IL33 negatively regulates cardiomyocyte hypertrophy and promotes cardioprotection. These agonists signal through a common co-receptor but, in cardiomyocytes, IL1β more potently activates mitogen-activated protein kinases and NFκB, pathways that regulate gene expression. We compared the effects of external application of IL1β and IL33 on the cardiomyocyte transcriptome. Neonatal rat cardiomyocytes were exposed to IL1β or IL33 (0.5, 1 or 2h). Transcriptomic profiles were determined using Affymetrix rat genome 230 2.0 microarrays and data were validated by quantitative PCR. IL1β induced significant changes in more RNAs than IL33 and, generally, to a greater degree. It also had a significantly greater effect in downregulating mRNAs and in regulating mRNAs associated with selected pathways. IL33 had a greater effect on a small, select group of specific transcripts. Thus, differences in intensity of intracellular signals can deliver qualitatively different responses. Quantitatively different responses in production of receptor agonists and transcription factors may contribute to qualitative differences at later times resulting in different phenotypic cellular responses.

Isolation and characterization of an AINTEGUMENTA-like gene in different coconut (Cocos nucifera L.) varieties from Sri Lanka

Development of an efficient tissue culture protocol in coconut is hampered by numerous technical constraints. Thus a greater understanding of the fundamental aspects of embryogenesis is essential. The role of AINTEGUMENTA-like genes in embryogenesis has been elucidated not only in model plants but also in economically important crops. A coconut gene, CnANT, that encodes two APETALA2 (AP2) domains and a conserved linker region similar to those of the BABY BOOM transcription factor was cloned, characterized, and its tissue specific expression was examined. The full-length cDNA of 1,780 bp contains a 1,425-bp open reading frame that encodes a putative peptide of 474 amino acids. The genomic DNA sequence includes 2,317 bp and consists of nine exons interrupted by eight introns. The exon/intron organization of CnANT is similar to that of homologous genes in other plant species. Analysis of differential tissue expression by real-time polymerase chain reaction indicated that CnANT is expressed more highly in in vitro grown tissues than in other vegetative tissues. Sequence comparison of the genomic sequence of CnANT in different coconut varieties revealed one single nucleotide polymorphism and one indel in the first exon and first intron, respectively, which differentiate the Tall group of trees from Dwarfs. The indel sequence, which can be considered a simple sequence repeats marker, was successfully used to distinguish the Tall and Dwarf groups as well as to develop a marker system, which may be of value in the identification of parental varieties that are used in coconut breeding programs in Sri Lanka.

Construction of a Festuca pratensis BAC library for map-based cloning in Festulolium substitution lines

Introgression in Festulolium is a potentially powerful tool to isolate genes for a large number of traits which differ between Festuca pratensis Huds. and Lolium perenne L. Not only are hybrids between the two species fertile, but the two genomes can be distinguished by genomic in situ hybridisation and a high frequency of recombination occurs between homoeologous chromosomes and chromosome segments. By a programme of introgression and a series of backcrosses, L. perenne lines have been produced which contain small F. pratensis substitutions. This material is a rich source of polymorphic markers targeted towards any trait carried on the F. pratensis substitution not observed in the L. perenne background. We describe here the construction of an F. pratensis BAC library, which establishes the basis of a map-based cloning strategy in L. perenne. The library contains 49,152 clones, with an average insert size of 112 kbp, providing coverage of 2.5 haploid genome equivalents. We have screened the library for eight amplified fragment length polymorphism (AFLP) derived markers known to be linked to an F. pratensis gene introgressed into L. perenne and conferring a staygreen phenotype as a consequence of a mutation in primary chlorophyll catabolism. While for four of the markers it was possible to identify bacterial artificial chromosome (BAC) clones, the other four AFLPs were too repetitive to enable reliable identification of locus-specific BACs. Moreover, when the four BACs were partially sequenced, no obvious coding regions could be identified. This contrasted to BACs identified using cDNA sequences, when multiple genes were identified on the same BAC.

Error, reproducibility and sensitivity: a pipeline for data processing of Agilent oligonucleotide expression arrays

Background: Expression microarrays are increasingly used to obtain large scale transcriptomic information on a wide range of biological samples. Nevertheless, there is still much debate on the best ways to process data, to design experiments and analyse the output. Furthermore, many of the more sophisticated mathematical approaches to data analysis in the literature remain inaccessible to much of the biological research community. In this study we examine ways of extracting and analysing a large data set obtained using the Agilent long oligonucleotide transcriptomics platform, applied to a set of human macrophage and dendritic cell samples. Results: We describe and validate a series of data extraction, transformation and normalisation steps which are implemented via a new R function. Analysis of replicate normalised reference data demonstrate that intrarray variability is small (only around 2 of the mean log signal), while interarray variability from replicate array measurements has a standard deviation (SD) of around 0.5 log(2) units (6 of mean). The common practise of working with ratios of Cy5/Cy3 signal offers little further improvement in terms of reducing error. Comparison to expression data obtained using Arabidopsis samples demonstrates that the large number of genes in each sample showing a low level of transcription reflect the real complexity of the cellular transcriptome. Multidimensional scaling is used to show that the processed data identifies an underlying structure which reflect some of the key biological variables which define the data set. This structure is robust, allowing reliable comparison of samples collected over a number of years and collected by a variety of operators. Conclusions: This study outlines a robust and easily implemented pipeline for extracting, transforming normalising and visualising transcriptomic array data from Agilent expression platform. The analysis is used to obtain quantitative estimates of the SD arising from experimental (non biological) intra- and interarray variability, and for a lower threshold for determining whether an individual gene is expressed. The study provides a reliable basis for further more extensive studies of the systems biology of eukaryotic cells.

Changes in gene expression in Arabidopsis shoots during phosphate starvation and the potential for developing smart plants

Our aim was to generate and prove the concept of "smart" plants to monitor plant phosphorus (P) status in Arabidopsis. Smart plants can be genetically engineered by transformation with a construct containing the promoter of a gene up-regulated specifically by P starvation in an accessible tissue upstream of a marker gene such as beta-glucuronidase (GUS). First, using microarrays, we identified genes whose expression changed more than 2.5-fold in shoots of plants growing hydroponically when P, but not N or K, was withheld from the nutrient solution. The transient changes in gene expression occurring immediately (4 h) after P withdrawal were highly variable, and many nonspecific, shock-induced genes were up-regulated during this period. However, two common putative cis-regulatory elements (a PHO-like element and a TATA box-like element) were present significantly more often in the promoters of genes whose expression increased 4 h after the withdrawal of P compared with their general occurrence in the promoters of all genes represented on the microarray. Surprisingly, the expression of only four genes differed between shoots of P-starved and -replete plants 28 h after P was withdrawn. This lull in differential gene expression preceded the differential expression of a new group of 61 genes 100 h after withdrawing P. A literature survey indicated that the expression of many of these "late" genes responded specifically to P starvation. Shoots had reduced P after 100 h, but growth was unaffected. The expression of SQD1, a gene involved in the synthesis of sulfolipids, responded specifically to P starvation and was increased 100 h after withdrawing P. Leaves of Arabidopsis bearing a SQD1::GUS construct showed increased GUS activity after P withdrawal, which was detectable before P starvation limited growth. Hence, smart plants can monitor plant P status. Transferring this technology to crops would allow precision management of P fertilization, thereby maintaining yields while reducing costs, conserving natural resources, and preventing pollution.

Gene expression changes in phosphorus deficient potato (Solanum tuberosum L.) leaves and the potential for diagnostic gene expression markers

Background: There are compelling economic and environmental reasons to reduce our reliance on inorganic phosphate (Pi) fertilisers. Better management of Pi fertiliser applications is one option to improve the efficiency of Pi fertiliser use, whilst maintaining crop yields. Application rates of Pi fertilisers are traditionally determined from analyses of soil or plant tissues. Alternatively, diagnostic genes with altered expression under Pi limiting conditions that suggest a physiological requirement for Pi fertilisation, could be used to manage Pifertiliser applications, and might be more precise than indirect measurements of soil or tissue samples. Results: We grew potato (Solanum tuberosum L.) plants hydroponically, under glasshouse conditions, to control their nutrient status accurately. Samples of total leaf RNA taken periodically after Pi was removed from the nutrient solution were labelled and hybridised to potato oligonucleotide arrays. A total of 1,659 genes were significantly differentially expressed following Pi withdrawal. These included genes that encode proteins involved in lipid, protein, and carbohydrate metabolism, characteristic of Pi deficient leaves and included potential novel roles for genes encoding patatin like proteins in potatoes. The array data were analysed using a support vector machine algorithm to identify groups of genes that could predict the Pi status of the crop. These groups of diagnostic genes were tested using field grown potatoes that had either been fertilised or unfertilised. A group of 200 genes could correctly predict the Pi status of field grown potatoes. Conclusions: This paper provides a proof-of-concept demonstration for using microarrays and class prediction tools to predict the Pi status of a field grown potato crop. There is potential to develop this technology for other biotic and abiotic stresses in field grown crops. Ultimately, a better understanding of crop stresses may improve our management of the crop, improving the sustainability of agriculture.

Cardiac protein kinases: the cardiomyocyte kinome and differential kinase expression in human failing hearts

Aims. Protein kinases are potential therapeutic targets for heart failure, but most studies of cardiac protein kinases derive from other systems, an approach that fails to account for specific kinases expressed in the heart and the contractile cardiomyocytes. We aimed to define the cardiomyocyte kinome (i.e. the protein kinases expressed in cardiomyocytes) and identify kinases with altered expression in human failing hearts. Methods and Results. Expression profiling (Affymetrix microarrays) detected >400 protein kinase mRNAs in rat neonatal ventricular myocytes (NVMs) and/or adult ventricular myocytes (AVMs), 32 and 93 of which were significantly upregulated or downregulated (>2-fold), respectively, in AVMs. Data for AGC family members were validated by qPCR. Proteomics analysis identified >180 cardiomyocyte protein kinases, with high relative expression of mitogen-activated protein kinase cascades and other known cardiomyocyte kinases (e.g. CAMKs, cAMP-dependent protein kinase). Other kinases are poorly-investigated (e.g. Slk, Stk24, Oxsr1). Expression of Akt1/2/3, BRaf, ERK1/2, Map2k1, Map3k8, Map4k4, MST1/3, p38-MAPK, PKCδ, Pkn2, Ripk1/2, Tnni3k and Zak was confirmed by immunoblotting. Relative to total protein, Map3k8 and Tnni3k were upregulated in AVMs vs NVMs. Microarray data for human hearts demonstrated variation in kinome expression that may influence responses to kinase inhibitor therapies. Furthermore, some kinases were upregulated (e.g. NRK, JAK2, STK38L) or downregulated (e.g. MAP2K1, IRAK1, STK40) in human failing hearts. Conclusions. This characterization of the spectrum of kinases expressed in cardiomyocytes and the heart (cardiomyocyte and cardiac kinomes) identified novel kinases, some of which are differentially expressed in failing human hearts and could serve as potential therapeutic targets.

Arabinose and protocatechuate catabolism genes are important for growth of Rhizobium leguminosarum biovar viciae in the pea rhizosphere

Background and aims: To form nitrogen-fixing nodules on pea roots, Rhizobium leguminosarum biovar viciae must be competitive in the rhizosphere. Our aim was to identify genes important for rhizosphere fitness. Methods: Signature-tagged mutants were screened using microarrays to identify mutants reduced for growth in pea rhizospheres. Candidate mutants were assessed relative to controls for growth in minimal medium, growth in pea rhizospheres and for infection of peas in mixed inoculants. Mutated genes were identified by DNA sequencing and confirmed by transduction. Results: Of 5508 signature-tagged mutants, microarrays implicated 50 as having decreased rhizosphere fitness. Growth tests identified six mutants with rhizosphere-specific phenotypes. The mutation in one of the genes (araE) was in an arabinose catabolism operon and blocked growth on arabinose. The mutation in another gene (pcaM), encoding a predicted solute binding protein for protocatechuate and hydroxybenzoate uptake, decreased growth on protocatechuate. Both mutants were decreased for nodule infection competitiveness with mixed inoculants, but nodulated peas normally when inoculated alone. Other mutants with similar phenotypes had mutations predicted to affect secondary metabolism. Conclusions: Catabolism of arabinose and protocatechuate in the pea rhizosphere is important for competitiveness of R.l. viciae. Other genes predicted to be involved in secondary metabolism are also important.

Relationships Among Murine CD11chigh Dendritic Cell Subsets as Revealed by Baseline Gene Expression Patterns

The functional relationships and properties of different subtypes of dendritic cells (DC) remain largely undefined. To better characterize these cells, we used global gene analysis to determine gene expression patterns among murine CD11c(high) DC subsets. CD4(+), CD8alpha(+), and CD8alpha(-) CD4(-) (double negative (DN)) DC were purified from spleens of normal C57/BL6 mice and analyzed using Affymetrix microarrays. The CD4(+) and CD8alpha(+) DC subsets showed distinct basal expression profiles differing by >200 individual genes. These included known DC subset markers as well as previously unrecognized, differentially expressed CD Ags such as CD1d, CD5, CD22, and CD72. Flow cytometric analysis confirmed differential expression in nine of nine cases, thereby validating the microarray analysis. Interestingly, the microarray expression profiles for DN cells strongly resembled those of CD4(+) DC, differing from them by <25 genes. This suggests that CD4(+) and DN DC are closely related phylogenetically, whereas CD8alpha(+) DC represent a more distant lineage, supporting the historical distinction between CD8alpha(+) and CD8alpha(-) DC. However, staining patterns revealed that in contrast to CD4(+) DC, the DN subset is heterogeneous and comprises at least two subpopulations. Gene Ontology and literature mining analyses of genes expressed differentially among DC subsets indicated strong associations with immune response parameters as well as cell differentiation and signaling. Such associations offer clues to possible unique functions of the CD11c(high) DC subsets that to date have been difficult to define as rigid distinctions.

Transcriptomic analysis of Enterohaemorrhagic Escherichia coli O157:H7 in response to plant extracts

Enterohaemorrhagic Escherichia coli (EHEC) are a group of food and contact-borne pathogens responsible for haemorrhagic colitis. The bacteria can be transmitted by contaminated meat, but importantly, also by plants. The bacteria can use plants as an alternative host, where they associate with both the leaves and the roots. Colonisation in the rhizosphere of plants is thought to be the main habitat for colonisation. Four different plant species, commonly associated with EHEC outbreaks, were infected with EHEC O157:H7 isolates Sakai and TUV 93-0 over ten days to assess the colonisation potential of the bacteria in both the phyllosphere and rhizosphere of plants. The rhizosphere was found to sustain a higher population level of bacteria over time in comparison to the phyllosphere, yet both strains were unable to utilize root exudates for growth. Global gene expression changes of EHEC O157:H7 strain Sakai were measured in response to plant extracts such as leaf lysates, root exudates and leaf cell wall polysaccharides from spinach cultivar Amazon and lettuce cultivar Salinas. Microarrays analysis showed a significant change in expression of 17 % of genes on exposure to leaf lysates of spinach. A more specific response was seen to spinach leaf cell wall polysaccharides with only a 1.5 % change. In contrast, when exposed to lettuce leaf cell wall polysaccharides a higher change of 4.8 % was seen. Genes that were differentially expressed belonged to multiple functional groups, including metabolism, indicating the utilization of plant-specific polysaccharides. Several areas of further investigation have been determined from this project, including the importance of culturing bacterial strains at a relevant temperature, the proposed lack of the type III secretion system in plant colonization by EHEC O157:H7 and the utilization of plant components for growth and persistence in the plant environment.

Transcriptome analysis for discovering candidate genes involve in embryogenesis in coconut (Cocos nucifera L.) through 454 pyrosequencing

Coconut, Cocos nucifera L. is a major plantation crop, which ensures income for millions of people in the tropical region. Detailed molecular studies on zygotic embryo development would provide valuable clues for the identification of molecular markers to improve somatic embryogenesis. Since there is no ongoing genome project for this species, coconut expressed sequence tags (EST) would be an interesting technique to identify important coconut embryo specific genes as well as other functional genes in different biochemical pathways. The goal of this study was to analyse the ESTs by examining the transcriptome data of the different embryo tissue types together with one somatic tissue. Here, four cDNA libraries from immature embryo, mature embryo, microspore derived embryo and mature leaves were constructed. cDNA was sequenced by the Roche-454 GS-FLX system and assembled into 32621 putative unigenes and 155017 singletons. Of these unigenes, 18651 had significant sequence similarities to non-redundant protein database, from which 16153 were assigned to one or more gene ontology categories. Homologue genes, which are responsible for embryo development such as chitinase, beta-1,3-glucanase, ATP synthase CF0 subunit, thaumatin-like protein and metallothionein-like protein were identified among the embryo EST collection. Of the unigenes, 6694 were mapped into 139 KEGG pathways including carbohydrate metabolism, energy metabolism, lipid metabolism, amino acid metabolism and nucleotide metabolism. This collection of 454-derived EST data generated from different tissue types provides a significant resource for genome wide studies and gene discovery of coconut, a non-model species.