A multi-level multi-scale approach to study essential genes in Mycobacterium tuberculosis

Background: The set of indispensable genes that are required by an organism to grow and sustain life are termed as essential genes. There is a strong interest in identification of the set of essential genes, particularly in pathogens, not only for a better understanding of the pathogen biology, but also for identifying drug targets and the minimal gene set for the organism. Essentiality is inherently a systems property and requires consideration of the system as a whole for their identification. The available experimental approaches capture some aspects but each method comes with its own limitations. Moreover, they do not explain the basis for essentiality in most cases. A powerful prediction method to recognize this gene pool including rationalization of the known essential genes in a given organism would be very useful. Here we describe a multi-level multi-scale approach to identify the essential gene pool in a deadly pathogen, Mycobacterium tuberculosis. Results: The multi-level workflow analyses the bacterial cell by studying (a) genome-wide gene expression profiles to identify the set of genes which show consistent and significant levels of expression in multiple samples of the same condition, (b) indispensability for growth by using gene expression integrated flux balance analysis of a genome-scale metabolic model, (c) importance for maintaining the integrity and flow in a protein-protein interaction network and (d) evolutionary conservation in a set of genomes of the same ecological niche. In the gene pool identified, the functional basis for essentiality has been addressed by studying residue level conservation and the sub-structure at the ligand binding pockets, from which essential amino acid residues in that pocket have also been identified. 283 genes were identified as essential genes with high-confidence. An agreement of about 73.5% is observed with that obtained from the experimental transposon mutagenesis technique. A large proportion of the identified genes belong to the class of intermediary metabolism and respiration. Conclusions: The multi-scale, multi-level approach described can be generally applied to other pathogens as well. The essential gene pool identified form a basis for designing experiments to probe their finer functional roles and also serve as a ready shortlist for identifying drug targets.

The N-terminal region containing the zinc finger domain of tobacco streak virus coat protein is essential for the formation of virus-like particles

Tobacco streak virus (TSV), a member of the genus Ilarvirus (family Bromoviridae), has a tripartite genome and forms quasi-isometric virions. All three viral capsids, encapsidating RNA 1, RNA 2 or RNA 3 and subgenomic RNA 4, are constituted of a single species of coat protein (CP). Formation of virus-like particles (VLPs) could be observed when the TSV CP gene was cloned and the recombinant CP (rCP) was expressed in E. coli. TSV VLPs were found to be stabilized by Zn2+ ions and could be disassembled in the presence of 500 mM CaCl2. Mutational analysis corroborated previous studies that showed that an N-terminal arginine-rich motif was crucial for RNA binding; however, the results presented here demonstrate that the presence of RNA is not a prerequisite for assembly of TSV VLPs. Instead, the N-terminal region containing the zinc finger domain preceding the arginine-rich motif is essential for assembly of these VLPs.

A MatLAB Based Virtual Phantom for 2D Electrical Impedance Tomography (MatVP2DEIT): Studying the Medical Electrical Impedance Tomography Reconstruction in Computer

Electrical Impedance Tomography (EIT) is a computerized medical imaging technique which reconstructs the electrical impedance images of a domain under test from the boundary voltage-current data measured by an EIT electronic instrumentation using an image reconstruction algorithm. Being a computed tomography technique, EIT injects a constant current to the patient's body through the surface electrodes surrounding the domain to be imaged (Omega) and tries to calculate the spatial distribution of electrical conductivity or resistivity of the closed conducting domain using the potentials developed at the domain boundary (partial derivative Omega). Practical phantoms are essentially required to study, test and calibrate a medical EIT system for certifying the system before applying it on patients for diagnostic imaging. Therefore, the EIT phantoms are essentially required to generate boundary data for studying and assessing the instrumentation and inverse solvers a in EIT. For proper assessment of an inverse solver of a 2D EIT system, a perfect 2D practical phantom is required. As the practical phantoms are the assemblies of the objects with 3D geometries, the developing of a practical 2D-phantom is a great challenge and therefore, the boundary data generated from the practical phantoms with 3D geometry are found inappropriate for assessing a 2D inverse solver. Furthermore, the boundary data errors contributed by the instrumentation are also difficult to separate from the errors developed by the 3D phantoms. Hence, the errorless boundary data are found essential to assess the inverse solver in 2D EIT. In this direction, a MatLAB-based Virtual Phantom for 2D EIT (MatVP2DEIT) is developed to generate accurate boundary data for assessing the 2D-EIT inverse solvers and the image reconstruction accuracy. MatVP2DEIT is a MatLAB-based computer program which simulates a phantom in computer and generates the boundary potential data as the outputs by using the combinations of different phantom parameters as the inputs to the program. Phantom diameter, inhomogeneity geometry (shape, size and position), number of inhomogeneities, applied current magnitude, background resistivity, inhomogeneity resistivity all are set as the phantom variables which are provided as the input parameters to the MatVP2DEIT for simulating different phantom configurations. A constant current injection is simulated at the phantom boundary with different current injection protocols and boundary potential data are calculated. Boundary data sets are generated with different phantom configurations obtained with the different combinations of the phantom variables and the resistivity images are reconstructed using EIDORS. Boundary data of the virtual phantoms, containing inhomogeneities with complex geometries, are also generated for different current injection patterns using MatVP2DEIT and the resistivity imaging is studied. The effect of regularization method on the image reconstruction is also studied with the data generated by MatVP2DEIT. Resistivity images are evaluated by studying the resistivity parameters and contrast parameters estimated from the elemental resistivity profiles of the reconstructed phantom domain. Results show that the MatVP2DEIT generates accurate boundary data for different types of single or multiple objects which are efficient and accurate enough to reconstruct the resistivity images in EIDORS. The spatial resolution studies show that, the resistivity imaging conducted with the boundary data generated by MatVP2DEIT with 2048 elements, can reconstruct two circular inhomogeneities placed with a minimum distance (boundary to boundary) of 2 mm. It is also observed that, in MatVP2DEIT with 2048 elements, the boundary data generated for a phantom with a circular inhomogeneity of a diameter less than 7% of that of the phantom domain can produce resistivity images in EIDORS with a 1968 element mesh. Results also show that the MatVP2DEIT accurately generates the boundary data for neighbouring, opposite reference and trigonometric current patterns which are very suitable for resistivity reconstruction studies. MatVP2DEIT generated data are also found suitable for studying the effect of the different regularization methods on reconstruction process. Comparing the reconstructed image with an original geometry made in MatVP2DEIT, it would be easier to study the resistivity imaging procedures as well as the inverse solver performance. Using the proposed MatVP2DEIT software with modified domains, the cross sectional anatomy of a number of body parts can be simulated in PC and the impedance image reconstruction of human anatomy can be studied.

IS A DEEP ONE-CELL MERIDIONAL CIRCULATION ESSENTIAL FOR THE FLUX TRANSPORT SOLAR DYNAMO?

The solar activity cycle is successfully modeled by the flux transport dynamo, in which the meridional circulation of the Sun plays an important role. Most of the kinematic dynamo simulations assume a one-cell structure of the meridional circulation within the convection zone, with the equatorward return flow at its bottom. In view of the recent claims that the return flow occurs at a much shallower depth, we explore whether a meridional circulation with such a shallow return flow can still retain the attractive features of the flux transport dynamo (such as a proper butterfly diagram, the proper phase relation between the toroidal and poloidal fields). We consider additional cells of the meridional circulation below the shallow return flow-both the case of multiple cells radially stacked above one another and the case of more complicated cell patterns. As long as there is an equatorward flow in low latitudes at the bottom of the convection zone, we find that the solar behavior is approximately reproduced. However, if there is either no flow or a poleward flow at the bottom of the convection zone, then we cannot reproduce solar behavior. On making the turbulent diffusivity low, we still find periodic behavior, although the period of the cycle becomes unrealistically large. In addition, with a low diffusivity, we do not get the observed correlation between the polar field at the sunspot minimum and the strength of the next cycle, which is reproduced when diffusivity is high. On introducing radially downward pumping, we get a more reasonable period and more solar-like behavior even with low diffusivity.

THE KNICKKOPF DOMON DOMAIN IS ESSENTIAL FOR CUTICLE DIFFERENTIATION IN Drosophila melanogaster

The dopamine monoxygenase N-terminal (DOMON) domain is found in extracellular proteins across several eukaryotic and prokaryotic taxa. It has been proposed that this domain binds to heme or sugar moieties. Here, we have analyzed the role of four highly conserved amino acids in the DOMON domain of the Drosophila melanogaster Knickkopf protein that is inserted into the apical plasma membrane and assists extracellular chitin organization. In principal, we generated Knickkopf versions with exchanged residues tryptophan(299,) methionine(333), arginine(401), or histidine(437), and scored for the ability of the respective engineered protein to normalize the knickkopf mutant phenotype. Our results confirm the absolute necessity of tryptophan(299,) methionine(333), and histidine(437) for Knickkopf function and stability, the latter two being predicted to be critical for heme binding. In contrast, arginine(401) is required for full efficiency of Knickkopf activity. Taken together, our genetic data support the prediction of these residues to mediate the function of Knickkopf during cuticle differentiation in insects. Hence, the DOMON domain is apparently an essential factor contributing to the construction of polysaccharide-based extracellular matrices.

Trm1p, a Zn(II)(2)Cys(6)-type transcription factor, is essential for the transcriptional activation of genes of methanol utilization pathway, in Pichia pastoris

The zinc finger transcription factors Mxr1p and Rop are key regulators of methanol metabolism in the methylotrophic yeast, Pichia pastoris, while Trm1p and Trm2p regulate methanol metabolism in Candida boidinii. Here, we demonstrate that Trm1p is essential for the expression of genes of methanol utilization (mut) pathway in P. pastoris as well. Expression of AOXI and other genes of mut pathway is severely compromised in P. pastoris Delta Trm1 strain resulting in impaired growth on media containing methanol as the sole source of carbon. Trm1p localizes to the nucleus of cells cultured on glucose or methanol. The zinc finger domain of Mxr1p but not Trm1p binds to AOXI promoter sequences in vitro, indicating that these two positive regulators act by different mechanisms. We conclude that both Trm1p and Mxr1p are essential for the expression of genes of mut pathway in P. pastoris and the mechanism of transcriptional regulation of mut pathway may be similar in P. pastoris and C. boidinii. (C) 2014 Elsevier Inc. All rights reserved.

Methylation Silencing of ULK2, an Autophagy Gene, Is Essential for Astrocyte Transformation and Tumor Growth

Glioblastoma (GBM) is the most aggressive type of brain tumor and shows very poor prognosis. Here, using genome-wide methylation analysis, we show that G-CIMP+ and G-CIMP-subtypes enrich distinct classes of biological processes. One of the hypermethylated genes in GBM, ULK2, an upstream autophagy inducer, was found to be down-regulated in GBM. Promoter hypermethylation of ULK2 was confirmed by bisulfite sequencing. GBM and glioma cell lines had low levels of ULK2 transcripts, which could be reversed upon methylation inhibitor treatment. ULK2 promoter methylation and transcript levels showed significant negative correlation. Ectopic overexpression of ULK2-induced autophagy, which further enhanced upon nutrient starvation or temozolomide chemotherapy. ULK2 also inhibited the growth of glioma cells, which required autophagy induction as kinase mutant of ULK2 failed to induce autophagy and inhibit growth. Furthermore, ULK2 induced autophagy and inhibited growth in Ras-transformed immortalized Baby Mouse Kidney (iBMK) ATG5(+/+) but not in autophagy-deficient ATG5(-/-) cells. Growth inhibition due to ULK2 induced high levels of autophagy under starvation or chemotherapy utilized apoptotic cell death but not at low levels of autophagy. Growth inhibition by ULK2 also appears to involve catalase degradation and reactive oxygen species generation. ULK2 overexpression inhibited anchorage independent growth, inhibited astrocyte transformation in vitro and tumor growth in vivo. Of all autophagy genes, we found ULK2 and its homologue ULK1 were only down-regulated in all grades of glioma. Thus these results altogether suggest that inhibition of autophagy by ULK1/2 down-regulation is essential for glioma development.

The polycomb group gene EMF2B is essential for maintenance of floral meristem determinacy in rice

Polycomb Repressive Complex 2 (PRC2) represses the transcriptional activity of target genes through trimethylation of lysine 27 of histone H3. The functions of plant PRC2 have been chiefly described in Arabidopsis, but specific functions in other plant species, especially cereals, are still largely unknown. Here we characterize mutants in the rice EMF2B gene, an ortholog of the Arabidopsis EMBRYONIC FLOWER2 (EMF2) gene. Loss of EMF2B in rice results in complete sterility, and mutant flowers have severe floral organ defects and indeterminacy that resemble loss-of-function mutants in E-function floral organ specification genes. Transcriptome analysis identified the E-function genes OsMADS1, OsMADS6 and OsMADS34 as differentially expressed in the emf2b mutant compared with wild type. OsMADS1 and OsMADS6, known to be required for meristem determinacy in rice, have reduced expression in the emf2b mutant, whereas OsMADS34 which interacts genetically with OsMADS1 was ectopically expressed. Chromatin immunoprecipitation for H3K27me3 followed by quantitative (q)RT-PCR showed that all three genes are presumptive targets of PRC2 in the meristem. Therefore, in rice, and possibly other cereals, PRC2 appears to play a major role in floral meristem determinacy through modulation of the expression of E-function genes.

Beadex Function in the Motor Neurons Is Essential for Female Reproduction in Drosophila melanogaster

Drosophila melanogaster has served as an excellent model system for understanding the neuronal circuits and molecular mechanisms regulating complex behaviors. The Drosophila female reproductive circuits, in particular, are well studied and can be used as a tool to understand the role of novel genes in neuronal function in general and female reproduction in particular. In the present study, the role of Beadex, a transcription co-activator, in Drosophila female reproduction was assessed by generation of mutant and knock down studies. Null allele of Beadex was generated by transposase induced excision of P-element present within an intron of Beadex gene. The mutant showed highly compromised reproductive abilities as evaluated by reduced fecundity and fertility, abnormal oviposition and more importantly, the failure of sperm release from storage organs. However, no defect was found in the overall ovariole development. Tissue specific, targeted knock down of Beadex indicated that its function in neurons is important for efficient female reproduction, since its neuronal knock down led to compromised female reproductive abilities, similar to Beadex null females. Further, different neuronal class specific knock down studies revealed that Beadex function is required in motor neurons for normal fecundity and fertility of females. Thus, the present study attributes a novel and essential role for Beadex in female reproduction through neurons.

Plasmodium berghei glycine cleavage system T-protein is non-essential for parasite survival in vertebrate and invertebrate hosts

T-protein, an aminomethyltransferase, represents one of the four components of glycine cleavage system (GCS) and catalyzes the transfer of methylene group from H-protein intermediate to tetrahydrofolate (THF) forming N-5, N-10-methylene THF (CH2-THF) with the release of ammonia. The malaria parasite genome encodes T-, H- and L-proteins, but not P-protein which is a glycine decarboxylase generating the aminomethylene group. A putative GCS has been considered to be functional in the parasite mitochondrion despite the absence of a detectable P-protein homologue. In the present study, the mitochondrial localization of T-protein in the malaria parasite was confirmed by immunofluorescence and its essentiality in the entire parasite life cycle was studied by targeting the T-protein locus in Plasmodium berghei (Pb). PbT knock out parasites did not show any growth defect in asexual, sexual and liver stages indicating that the T-protein is dispensable for parasite survival in vertebrate and invertebrate hosts. The absence of P-protein homologue and the non-essentiality of T protein suggest the possible redundancy of GCS activity in the malaria parasite. Nevertheless, the H- and L-proteins of GCS could be essential for malaria parasite because of their involvement in alpha-lcetoacid dehydrogenase reactions. (C) 2014 Elsevier B.V. All rights reserved.

Genome-wide methylation profiling identifies an essential role of reactive oxygen species in pediatric glioblastoma multiforme and validates a methylome specific for H3 histone family 3A with absence of G-CIMP/isocitrate dehydrogenase 1 mutation

Background. Pediatric glioblastoma multiforme (GBM) is rare, and there is a single study, a seminal discovery showing association of histone H3.3 and isocitrate dehydrogenase (IDH) 1 mutation with a DNA methylation signature. The present study aims to validate these findings in an independent cohort of pediatric GBM, compare it with adult GBM, and evaluate the involvement of important functionally altered pathways. Methods. Genome-wide methylation profiling of 21 pediatric GBM cases was done and compared with adult GBM data (GSE22867). We performed gene mutation analysis of IDH1 and H3 histone family 3A (H3F3A), status evaluation of glioma cytosine-phosphate-guanine island methylator phenotype (G-CIMP), and Gene Ontology analysis. Experimental evaluation of reactive oxygen species (ROS) association was also done. Results. Distinct differences were noted between methylomes of pediatric and adult GBM. Pediatric GBM was characterized by 94 hypermethylated and 1206 hypomethylated cytosine-phosphate-guanine (CpG) islands, with 3 distinct clusters, having a trend to prognostic correlation. Interestingly, none of the pediatric GBM cases showed G-CIMP/IDH1 mutation. Gene Ontology analysis identified ROS association in pediatric GBM, which was experimentally validated. H3F3A mutants (36.4%; all K27M) harbored distinct methylomes and showed enrichment of processes related to neuronal development, differentiation, and cell-fate commitment. Conclusions. Our study confirms that pediatric GBM has a distinct methylome compared with that of adults. Presence of distinct clusters and an H3F3A mutation-specific methylome indicate existence of epigenetic subgroups within pediatric GBM. Absence of IDH1/G-CIMP status further indicates that findings in adult GBM cannot be simply extrapolated to pediatric GBM and that there is a strong need for identification of separate prognostic markers. A possible role of ROS in pediatric GBM pathogenesis is demonstrated for the first time and needs further evaluation.

Reservoir Operation with Fuzzy State Variables for Irrigation of Multiple Crops

This paper presents the development and application of a stochastic dynamic programming model with fuzzy state variables for irrigation of multiple crops. A fuzzy stochastic dynamic programming (FSDP) model is developed in which the reservoir storage and soil moisture of the crops are considered as fuzzy numbers, and the reservoir inflow is considered as a stochastic variable. The model is formulated with an objective of minimizing crop yield deficits, resulting in optimal water allocations to the crops by maintaining storage continuity and soil moisture balance. The standard fuzzy arithmetic method is used to solve all arithmetic equations with fuzzy numbers, and the fuzzy ranking method is used to compare two or more fuzzy numbers. The reservoir operation model is integrated with a daily-based water allocation model, which results in daily temporal variations of allocated water, soil moisture, and crop deficits. A case study of an existing Bhadra reservoir in Karnataka, India, is chosen for the model application. The FSDP is a more realistic model because it considers the uncertainty in discretization of state variables. The results obtained using the FSDP model are found to be more acceptable for the case study than those of the classical stochastic dynamic model and the standard operating model, in terms of 10-day releases from the reservoir and evapotranspiration deficit. (C) 2015 American Society of Civil Engineers.

Cybernetic Modeling Revisited: A Method for Inferring the Cybernetic Variables u(i) from Experimental Data

The cybernetic modeling framework for the growth of microorganisms provides for an elegant methodology to account for the unknown regulatory phenomena through the use of cybernetic variables for enzyme induction and activity. In this paper, we revisit the assumption of limited resources for enzyme induction (Sigma u(i) = 1) used in the cybernetic modeling framework by presenting a methodology for inferring the individual cybernetic variables u(i) from experimental data. We use this methodology to infer u(i) during the simultaneous consumption of glycerol and lactose by Escherichia coli and then model the fitness trade-offs involved in the recently discovered predictive regulation strategy of microorganisms.

A Viscoelastic Constitutive Model with Nonlinear Evolutionary Internal Variables

Based on the internal variable theory, a viscoelastic constitutive model of a highly deformable continuous medium is proposed. A set of second rank tensorial internal state variables corresponding to Biot's strain is introduced, and a nonlinear evolution

Essential and Non-Essential Work of Fracture of PI/SiO2 Hybrid Thin Films

Essential work of fracture (EWF) analysis is used to study the effect of the silica doping level on fracture toughness of polyimide/silica (PI/SiO2) hybrid films. By using double-edge-notched-tension (DENT) specimens with different ligament lengths, it seems that the introduction of silica additive can improve the specific essential work of fracture (w (e) ) of PI thin films, but the specific non-essential work of fracture (beta w (p) ) will decease significantly as the silica doping level increasing from 1 to 5 wt.%, and even lower than that of neat PI. The failure process of the fracture is investigated with online scanning electron microscope (SEM) observation and the parameters of non-essential work of fracture, beta and w (p) , are calculated based on finite element (FE) method.