Search

Tumor necrosis factor alpha (TNFα), a pleiotropic cytokine, helps maintain a balance between proliferation and apoptosis in normal cells. This balance is often sacrificed in a diseased cell, such as that of a cancer, by preferring survival phenotype over apoptosis. Restoring this balance requires a detailed understanding of the causal intracellular mechanisms that govern TNFαstimulated apoptotic response. In this study, we use a systems biology approach to unravel the interplay between the intracellular signaling markers that orchestrate apoptosis levels. Our approach deciphered the synergism between the early intracellular markers phosphorylated JNK (pJNK) and phosphorylated AKT (pAKT) that modulate the activation of Caspase3, an important apoptotic regulator. We demonstrate that this synergism depends critically on the survival pathway signaling mediated by NFκB which plays a dominant role in controlling the extent of the overall apoptotic response. By systematic inhibition of the signaling markers, we establish that the dynamic cross-talk between the pJNK and pAKT transients directs the apoptosis phenotype via accumulated Caspase3 response. Interestingly, superposition of the semi-quantitative correlation between apoptosis and Caspase3 transient levels on the proposed TNFαnetwork model permits quantification of the dynamic apoptotic response under different stimulation conditions. Thus, the predictive model can be leveraged towards arriving at useful insights that can identify potential targeted therapeutic strategies for altering apoptotic response.

Detection of levels of intracellular phospho-proteins is key to analyzing the dynamics of signal transduction in cellular systems. Cell-to-cell variability in the form of differences in protein level in each cell affects signaling and is implicated in prognosis of many diseases. Quantitative analysis of such variability necessitate measuring the protein levels at single-cell resolution. Single-cell intracellular protein abundance detection in statistically significant number of adherent cells for short time sampling points post stimulation using classical flow cytometry (FCM) technique has thus far been a challenge due to the detrimental effects of cell detachment methods on the cellular machinery. We systematically show that cell suspension obtained by noninvasive temperature-sensitive detachment of adherent cells is amenable to high-throughput phospho-ERK1/2 protein detection at single-cell level using FCM in these short time sampling points. We demonstrate this on three adherent cell lines, viz., HeLa, A549, and MCF7, from distinct lineages having characteristically different elasticity at 37 °C. In particular, we use a right combination of multiplexing via fluorescent cell barcoding (FCB) and intracellular antibody staining for simultaneous detection of phospho-ERK1/2 (pERK) stimulated by epidermal growth factor (EGF) in multiple samples. Based on systematic characterization using Alexa 350 dye, we arrive at two conditions that must be satisfied for correct implementation of FCB. Our study reveals that the temperature-sensitive detachment of HeLa cells correctly captures the expected pronounced bimodal pERK distribution as an early response to EGF, which the enzymatic treatment methods fail to detect. © 2018 International Society for Advancement of Cytometry.

Reference genes are indispensable for normalizing mRNA levels across samples in real-time quantitative PCR. Their expression levels vary under different experimental conditions and because of several inherent characteristics. Appropriate reference gene selection is thus critical for gene-expression studies. This study aimed at selecting optimal reference genes for gene-expression analysis of reticulocytes and at validating them in hereditary spherocytosis (HS) and β-thalassemia intermedia (βTI) patients. Seven reference genes (PGK1, MPP1, HPRT1, ACTB, GAPDH, RN18S1, and SDHA) were selected because of published reports. Real-time quantitative PCR was performed on reticulocytes in 20 healthy volunteers, 15 HS patients, and 10 βTI patients. Threshold cycle values were compared with fold-change method and RefFinder software. The stable reference genes recommended by RefFinder were validated with SLC4A1 and flow cytometric eosin-5′-maleimide binding assay values in HS patients and HBG2 and high performance liquid chromatography–derived percentage of hemoglobin F in βTI. Comprehensive ranking predicted MPP1 and GAPDH as optimal reference genes for reticulocytes that were not affected in HS and βTI. This was further confirmed on validation with eosin-5′-maleimide results and percentage of hemoglobin F in HS and βTI patients, respectively. Hence, MPP1 and GAPDH are good reference genes for reticulocyte expression studies compared with ACTB and RN18S1, the two most commonly used reference genes.

Genome scale metabolic model provides an overview of an organism's metabolic capability. These genome-specific metabolic reconstructions are based on identification of gene to protein to reaction (GPR) associations and, in turn, on homology with annotated genes from other organisms. Cyanobacteria are photosynthetic prokaryotes which have diverged appreciably from their nonphotosynthetic counterparts. They also show significant evolutionary divergence from plants, which are well studied for their photosynthetic apparatus. We argue that context-specific sequence and domain similarity can add to the repertoire of the GPR associations and significantly expand our view of the metabolic capability of cyanobacteria. We took an approach that combines the results of context-specific sequence-to-sequence similarity search with those of sequence-to-profile searches. We employ PSI-BLAST for the former, and CDD, Pfam, and COG for the latter. An optimization algorithm was devised to arrive at a weighting scheme to combine the different evidences with KEGG-annotated GPRs as training data. We present the algorithm in the form of software "Systematic, Homology-based Automated Re-annotation for Prokaryotes (SHARP)." We predicted 3,781 new GPR associations for the 10 prokaryotes considered of which eight are cyanobacteria species. These new GPR associations fall in several metabolic pathways and were used to annotate 7,718 gaps in the metabolic network. These new annotations led to discovery of several pathways that may be active and thereby providing new directions for metabolic engineering of these species for production of useful products. Metabolic model developed on such a reconstructed network is likely to give better phenotypic predictions.

Unicellular diazotrophic cyanobacteria such as Cyanothece sp. ATCC 51142 (henceforth Cyanothece), temporally separate the oxygen sensitive nitrogen fixation from oxygen evolving photosynthesis not only under diurnal cycles (LD) but also in continuous light (LL). However, recent reports demonstrate that the oscillations in LL occur with a shorter cycle time of similar to 11 h. We find that indeed, majority of the genes oscillate in LL with this cycle time. Genes that are upregulated at a particular time of day under diurnal cycle also get upregulated at an equivalent metabolic phase under LL suggesting tight coupling of various cellular events with each other and with the cell's metabolic status. A number of metabolic processes get upregulated in a coordinated fashion during the respiratory phase under LL including glycogen degradation, glycolysis, oxidative pentose phosphate pathway, and tricarboxylic acid cycle. These precede nitrogen fixation apparently to ensure sufficient energy and anoxic environment needed for the nitrogenase enzyme. Photosynthetic phase sees upregulation of photosystem II, carbonate transport, carbon concentrating mechanism, RuBisCO, glycogen synthesis and light harvesting antenna pigment biosynthesis. In Synechococcus elongates PCC 7942, a non-nitrogen fixing cyanobacteria, expression of a relatively smaller fraction of genes oscillates under LL condition with the major periodicity being 24 h. In contrast, the entire cellular machinery of Cyanothece orchestrates coordinated oscillation in anticipation of the ensuing metabolic phase in both LD and LL. These results may have important implications in understanding the timing of various cellular events and in engineering cyanobacteria for biofuel production.

This study investigates the influence of mixotrophy on physiology and metabolism by analysis of global gene expression in unicellular diazotrophic cyanobacterium Cyanothece sp. ATCC 51142 ( henceforth Cyanothece 51142). It was found that Cyanothece 51142 continues to oscillate between photosynthesis and respiration in continuous light under mixotrophy with cycle time of similar to 13 h. Mixotrophy is marked by an extended respiratory phase compared with photoautotrophy. It can be argued that glycerol provides supplementary energy for nitrogen fixation, which is derived primarily from the glycogen reserves during photoautotrophy. The genes of NDH complex, cytochrome c oxidase and ATP synthase are significantly overexpressed in mixotrophy during the day compared to autotrophy with synchronous expression of the bidirectional hydrogenase genes possibly to maintain redox balance. However, nitrogenase complex remains exclusive to nighttime metabolism concomitantly with uptake hydrogenase. This study throws light on interrelations between metabolic pathways with implications in design of hydrogen producer strains. (C) 2015 Elsevier Ltd. All rights reserved.

Cyanobacteria, a group of photosynthetic prokaryotes, oscillate between day and night time metabolisms with concomitant oscillations in gene expression in response to light/dark cycles (LD). The oscillations in gene expression have been shown to sustain in constant light (LL) with a free running period of 24 h in a model cyanobacterium Synechococcus elongatus PCC 7942. However, equivalent oscillations in metabolism are not reported under LL in this non-nitrogen fixing cyanobacterium. Here we focus on Cyanothece sp. ATCC 51142, a unicellular, nitrogen-fixing cyanobacterium known to temporally separate the processes of oxygenic photosynthesis and oxygen-sensitive nitrogen fixation. In a recent report, metabolism of Cyanothece 51142 has been shown to oscillate between photosynthetic and respiratory phases under LL with free running periods that are temperature dependent but significantly shorter than the circadian period. Further, the oscillations shift to circadian pattern at moderate cell densities that are concomitant with slower growth rates. Here we take this understanding forward and demonstrate that the utradian rhythm under LL sustains at much higher cell densities when grown under turbulent regimes that simulate flashing light effect. Our results suggest that the ultradian rhythm in metabolism may be needed to support higher carbon and nitrogen requirements of rapidly growing cells under LL. With a comprehensive Real time PCR based gene expression analysis we account for key regulatory interactions and demonstrate the interplay between clock genes and the genes of key metabolic pathways. Further, we observe that several genes that peak at dusk in Synechococcus peak at dawn in Cyanothece and vice versa. The circadian rhythm of this organism appears to be more robust with peaking of genes in anticipation of the ensuing photosynthetic and respiratory metabolic phases.

Data deposition to NCBI Genomes: This Whole Genome Shotgun project has been deposited at DDBJ/EMBL/GenBank under the accession AMXX00000000 (SMACv1.0, unscaffolded genome assembly). The version described in this paper is the first version (AMXX01000000). The scaffolded assembly (SMACv1.1) has been deposited at DDBJ/EMBL/GenBank under the accession AOUJ00000000, and is also the first version (AOUJ01000000). Strong biological interest in traits such as the acquisition and utilization of speech, cognitive abilities, and longevity catalyzed the utilization of two next-generation sequencing platforms to provide the first-draft de novo genome assembly for the large, new world parrot Ara macao (Scarlet Macaw). Despite the challenges associated with genome assembly for an outbred avian species, including 951,507 high-quality putative single nucleotide polymorphisms, the final genome assembly (>1.035 Gb) includes more than 997 Mb of unambiguous sequence data (excluding N's). Cytogenetic analyses including ZooFISH revealed complex rearrangements associated with two scarlet macaw macrochromosomes (AMA6, AMA7), which supports the hypothesis that translocations, fusions, and intragenomic rearrangements are key factors associated with karyotype evolution among parrots. In silico annotation of the scarlet macaw genome provided robust evidence for 14,405 nuclear gene annotation models, their predicted transcripts and proteins, and a complete mitochondrial genome. Comparative analyses involving the scarlet macaw, chicken, and zebra finch genomes revealed high levels of nucleotide-based conservation as well as evidence for overall genome stability among the three highly divergent species. Application of a new whole-genome analysis of divergence involving all three species yielded prioritized candidate genes and noncoding regions for parrot traits of interest (i.e., speech, intelligence, longevity) which were independently supported by the results of previous human GWAS studies. We also observed evidence for genes and noncoding loci that displayed extreme conservation across the three avian lineages, thereby reflecting their likely biological and developmental importance among birds.

Nitrogen fixing cyanobacteria are being increasingly explored for nitrogenase-dependent hydrogen production. Commercial success however will depend on the ability to grow these cultures at high cell densities. Photo-limitation at high cell densities leads to hindered photoautotrophic growth while turbulent conditions, which simulate flashing light effect, can lead to oxygen toxicity to the nitrogenase enzyme. Cyanothece sp. strain ATCC 51142, a known hydrogen producer, is reported to grow and fix nitrogen under moderately oxic conditions in shake flasks. In this study, we explore the growth and nitrogen fixing potential of this organism under turbulent conditions with volumetric oxygen mass transfer coefficient (K(L)a) values that are up to 20-times greater than in shake flasks. In a stirred vessel, the organism grows well in turbulent regime possibly due to a simulated flashing light effect with optimal growth at Reynolds number of approximately 35,000. A respiratory burst lasting for about 4h creates anoxic conditions intracellularly with near saturating levels of dissolved oxygen in the extracellular medium. This is concomitant with complete exhaustion of intracellular glycogen storage and upregulation of nifH and nifX, the genes encoding proteins of the nitrogenase complex. Further, the rhythmic oscillations in exhaust gas CO2 and O-2 profiles synchronize faithfully with those in biochemical parameters and gene expression thereby serving as an effective online monitoring tool. These results will have important implications in potential commercial success of nitrogenase-dependent hydrogen production by cyanobacteria. Biotechnol. Bioeng. 2013; 110:2371-2379. (c) 2013 Wiley Periodicals, Inc.

The white-tailed deer (Odocoileus virginianus) represents one of the most successful and widely distributed large mammal species within North America, yet very little nucleotide sequence information is available. We utilized massively parallel pyrosequencing of a reduced representation library (RRL) and a random shotgun library (RSL) to generate a complete mitochondrial genome sequence and identify a large number of putative single nucleotide polymorphisms (SNPs) distributed throughout the white-tailed deer nuclear and mitochondrial genomes. A SNP validation study designed to test specific classes of putative SNPs provides evidence for as many as 10,476 genome-wide SNPs in the current dataset. Based on cytogenetic evidence for homology between cow (Bos taurus) and white-tailed deer chromosomes, we demonstrate that a divergent genome may be used for estimating the relative distribution and density of de novo sequence contigs as well as putative SNPs for species without draft genome assemblies. Our approach demonstrates that bioinformatic tools developed for model or agriculturally important species may be leveraged to support next-generation research programs for species of biological, ecological and evolutionary importance. We also provide a functional annotation analysis for the de novo sequence contigs assembled from white-tailed deer pyrosequencing reads, a mitochondrial phylogeny involving 13,722 nucleotide positions for 10 unique species of Cervidae, and a median joining haplotype network as a putative representation of mitochondrial evolution in O. virginianus. The results of this study are expected to provide a detailed template enabling genome-wide sequence-based studies of threatened, endangered or conservationally important non-model organisms.

Life depends on the capacity of individual cells to respond effectively to cues about their changing internal and external environments. Cellular decision making and responses are orchestrated by complex molecular networks consisting of entities such as proteins or RNAs connected by interactions such as activation or synthesis. Information contained in primary databases and in the experimental literature relevant to these networks is so extensive and rapidly growing that it is increasingly difficult to integrate. As an aid to theoretical and experimental research, it is convenient to distill the inferences contained in the experimental literature and databases into knowledgebases that consist of annotated representations of biological pathways. Pathway building has been performed by individual groups studying a network of interest (e.g., Kitano's group who assembled an immune signaling pathway [1]) as well as by large bioinformatics consortia (e.g., the Reactome Project [2]) and commercial entities (e.g., Ingenuity Systems). Pathway building is the process of identifying and integrating the entities, interactions, and associated annotations, and populating the knowledgebase. Pathway construction can have either a data-driven objective (DDO) or a knowledge-driven objective (KDO). Data-driven pathway construction is used to generate relationship information of genes or proteins identified in a specific experiment such as a microarray study. Knowledge-driven pathway construction entails development of a detailed pathway knowledgebase for particular domains of interest, such as a cell type, disease, or system. To help researchers get their bearings in this field, in the subsequent sections we provide a brief, practical orientation to existing knowledgebases and to the methods of pathway construction and analysis.

Series MAPK enzymatic cascades, ubiquitously found in signaling networks, act as signal amplifiers and play a key role in processing information during signal transduction in cells. In activated cascades, cell-to-cell variability or noise is bound to occur and thereby strongly affects the cellular response. Commonly used linearization method (LM) applied to Langevin type stochastic model of the MAPK cascade fails to accurately predict intrinsic noise propagation in the cascade. We prove this by using extensive stochastic simulations for various ranges of biochemical parameters. This failure is due to the fact that the LM ignores the nonlinear effects on the noise. However, LM provides a good estimate of the extrinsic noise propagation. We show that the correct estimate of intrinsic noise propagation in signaling networks that contain at least one enzymatic step can be obtained only through stochastic simulations. Noise propagation in the cascade depends on the underlying biochemical parameters which are often unavailable. Based on a combination of global sensitivity analysis (GSA) and stochastic simulations, we developed a systematic methodology to characterize noise propagation in the cascade. GSA predicts that noise propagation in MAPK cascade is sensitive to the total number of upstream enzyme molecules and the total number of molecules of the two substrates involved in the cascade. We argue that the general systematic approach proposed and demonstrated on MAPK cascade must accompany noise propagation studies in biological networks.

Background In the yeast Saccharomyces cerevisiae, interactions between galactose, Gal3p, Gal80p, and Gal4p determine the transcriptional status of the genes required for the galactose utilization. Increase in the cellular galactose concentration causes the galactose molecules to bind onto Gal3p which, via Gal80p, activates Gal4p, which induces the GAL3 and GAL80 gene transcription. Recently, a linear time-invariant multi-input multi-output (MIMO) model of this GAL regulatory network has been proposed; the inputs being galactose and Gal4p, and the outputs being the active Gal4p and galactose utilization. Unfortunately, this model assumes the cell culture to be homogeneous, although it is not so in practice. We overcome this drawback by including more biochemical reactions, and derive a quadratic ordinary differential equation (ODE) based model. Results We show that the model, referred to above, does not exhibit bistability. We establish sufficiency conditions for the domain of attraction of an equilibrium point of our ODE model for the special case of full-state feedback controller. We observe that the GAL regulatory system of Kluyveromyces lactis exhibits an aberration of monotone nonlinearity and apply the Rantzer multipliers to establish a class of stabilizing controllers for this system. Conclusion Feedback in a GAL regulatory system can be used to enhance the cellular memory. We show that the system can be modeled as a quadratic nonlinear system for which the effect of feedback on the domain of attraction of the equilibrium point can be characterized using linear matrix inequality (LMI) conditions that are easily implementable in software. The benefit of this result is that a mathematically sound approach to the synthesis of full-state and partial-state feedback controllers to regulate the cellular memory is now possible, irrespective of the number of state-variables or parameters of interest.