Mitochondrial phosphoenolpyruvate carboxykinase (PEPCK-M) and serine biosynthetic pathway genes are co-ordinately increased during anabolic agent-induced skeletal muscle growth.

We aimed to identify novel molecular mechanisms for muscle growth during administration of anabolic agents. Growing pigs (Duroc/(Landrace/Large-White)) were administered Ractopamine (a beta-adrenergic agonist; BA; 20 ppm in feed) or Reporcin (recombinant growth hormone; GH; 10 mg/48 hours injected) and compared to a control cohort (feed only; no injections) over a 27-day time course (1, 3, 7, 13 or 27-days). Longissimus Dorsi muscle gene expression was analyzed using Agilent porcine transcriptome microarrays and clusters of genes displaying similar expression profiles were identified using a modified maSigPro clustering algorithm. Anabolic agents increased carcass (p = 0.002) and muscle weights (Vastus Lateralis: p < 0.001; Semitendinosus: p = 0.075). Skeletal muscle mRNA expression of serine/one-carbon/glycine biosynthesis pathway genes (Phgdh, Psat1 and Psph) and the gluconeogenic enzyme, phosphoenolpyruvate carboxykinase-M (Pck2/PEPCK-M), increased during treatment with BA, and to a lesser extent GH (p < 0.001, treatment x time interaction). Treatment with BA, but not GH, caused a 2-fold increase in phosphoglycerate dehydrogenase (PHGDH) protein expression at days 3 (p < 0.05) and 7 (p < 0.01), and a 2-fold increase in PEPCK-M protein expression at day 7 (p < 0.01). BA treated pigs exhibit a profound increase in expression of PHGDH and PEPCK-M in skeletal muscle, implicating a role for biosynthetic metabolic pathways in muscle growth.

The successful application of systems approaches in plant biology.

Plant biology has had longstanding successes from applying mathematical models to plant systems. Of the >160 models that have been developed to date, a closer study is made of crop models and more recent plant models. The latter focus on hormone response networks, metabolism, circadian clock, biomechanics of growth and new organ development. The multicellular and multiscale models have provided important and novel insights into the regulation of hormone distribution, tissue development and environmental sensing. Finally, the challenges faced when establishing multidisciplinary teams are introduced along with successful management strategies and techniques to overcome them.

Application of experimentally verified transcription factor binding sites models for computational analysis of ChIP-Seq data.

BACKGROUND: ChIP-Seq is widely used to detect genomic segments bound by transcription factors (TF), either directly at DNA binding sites (BSs) or indirectly via other proteins. Currently, there are many software tools implementing different approaches to identify TFBSs within ChIP-Seq peaks. However, their use for the interpretation of ChIP-Seq data is usually complicated by the absence of direct experimental verification, making it difficult both to set a threshold to avoid recognition of too many false-positive BSs, and to compare the actual performance of different models. RESULTS: Using ChIP-Seq data for FoxA2 binding loci in mouse adult liver and human HepG2 cells we compared FoxA binding-site predictions for four computational models of two fundamental classes: pattern matching based on existing training set of experimentally confirmed TFBSs (oPWM and SiteGA) and de novo motif discovery (ChIPMunk and diChIPMunk). To properly select prediction thresholds for the models, we experimentally evaluated affinity of 64 predicted FoxA BSs using EMSA that allows safely distinguishing sequences able to bind TF. As a result we identified thousands of reliable FoxA BSs within ChIP-Seq loci from mouse liver and human HepG2 cells. It was found that the performance of conventional position weight matrix (PWM) models was inferior with the highest false positive rate. On the contrary, the best recognition efficiency was achieved by the combination of SiteGA & diChIPMunk/ChIPMunk models, properly identifying FoxA BSs in up to 90% of loci for both mouse and human ChIP-Seq datasets. CONCLUSIONS: The experimental study of TF binding to oligonucleotides corresponding to predicted sites increases the reliability of computational methods for TFBS-recognition in ChIP-Seq data analysis. Regarding ChIP-Seq data interpretation, basic PWMs have inferior TFBS recognition quality compared to the more sophisticated SiteGA and de novo motif discovery methods. A combination of models from different principles allowed identification of proper TFBSs.

Sequence variants in the bovine gonadotrophin releasing hormone receptor gene and their associations with fertility.

Seven sequence variants (SVs) have been identified in exon 1 and in the promoter region upstream of the bovine gonadotrophin releasing hormone (GnRH) receptor gene, at nucleotides g.-331A>G, g.-108T>C, g.+206G>A, g.+260C>T, g.+341C>T, g.+383C>T and g.+410C>T relative to the translation start site. The SVs at nucleotides g.-108, g.260, g.341 and g.410 and those at g.206 and g.383 formed two groups with complete linkage disequilibrium within groups, but incomplete linkage disequilibrium between groups, and none of the SVs altered receptor amino acid sequence. The g.-108T>C allelic variants were associated with an approximately 0.4 day reduction in predicted transmitting ability for days to first service. None of the allelic variants affected the pattern of circulating LH following administration of GnRH. The g.260C>T alteration introduced a new transcription factor binding site in a region of DNA with relatively low nucleosome formation potential. The data suggest that selection for animals carrying the g.-108T>C group of alterations will improve fertility in the dairy cow.

Effective transcription factor binding site prediction using a combination of optimization, a genetic algorithm and discriminant analysis to capture distant interactions.

BACKGROUND: Reliable transcription factor binding site (TFBS) prediction methods are essential for computer annotation of large amount of genome sequence data. However, current methods to predict TFBSs are hampered by the high false-positive rates that occur when only sequence conservation at the core binding-sites is considered. RESULTS: To improve this situation, we have quantified the performance of several Position Weight Matrix (PWM) algorithms, using exhaustive approaches to find their optimal length and position. We applied these approaches to bio-medically important TFBSs involved in the regulation of cell growth and proliferation as well as in inflammatory, immune, and antiviral responses (NF-kappaB, ISGF3, IRF1, STAT1), obesity and lipid metabolism (PPAR, SREBP, HNF4), regulation of the steroidogenic (SF-1) and cell cycle (E2F) genes expression. We have also gained extra specificity using a method, entitled SiteGA, which takes into account structural interactions within TFBS core and flanking regions, using a genetic algorithm (GA) with a discriminant function of locally positioned dinucleotide (LPD) frequencies. To ensure a higher confidence in our approach, we applied resampling-jackknife and bootstrap tests for the comparison, it appears that, optimized PWM and SiteGA have shown similar recognition performances. Then we applied SiteGA and optimized PWMs (both separately and together) to sequences in the Eukaryotic Promoter Database (EPD). The resulting SiteGA recognition models can now be used to search sequences for BSs using the web tool, SiteGA. Analysis of dependencies between close and distant LPDs revealed by SiteGA models has shown that the most significant correlations are between close LPDs, and are generally located in the core (footprint) region. A greater number of less significant correlations are mainly between distant LPDs, which spanned both core and flanking regions. When SiteGA and optimized PWM models were applied together, this substantially reduced false positives at least at higher stringencies. CONCLUSION: Based on this analysis, SiteGA adds substantial specificity even to optimized PWMs and may be considered for large-scale genome analysis. It adds to the range of techniques available for TFBS prediction, and EPD analysis has led to a list of genes which appear to be regulated by the above TFs.

[Method SiteGA for the recognition of transcription factor binding sites].

A new approach, SiteGA, for the prediction of functional transcription factor binding sites has been developed. The approach is based on the detection of locally positioned dinucleotides by the genetic algorithm and discriminant analysis. The approach has been applied to recognize transcription factor binding sites involved in the regulation of immune responses and cell growth (AP-1, IRF1, ISGF3, NFkappaB, STAT1), obesity and lipid metabolism (HNF4, PPAR, SREBP), and the expression of steroidogenesis genes (SF-1). SiteGA is far superior in accuracy to the traditionally used method of position weight matrices. The approach was implemented in the web tool, SiteGA http://wwwmgs2. bionet.nsc.ru/mgs/programs/sitega.

Constructing an enzyme-centric view of metabolism.

MOTIVATION: The current paradigm for viewing metabolism, such as the Boehringer Chart or KEGG, takes a metabolite-centric view that is not ideal for genomics analysis because the same enzyme can appear in multiple places. Therefore an enzyme-centric view is also required. RESULTS: We have eliminated synonymous compound names taken from the ENZYME database ensuring that it is computationally parseable at all levels. Based on these results, we have written a software to create enzyme-centric graphs from reaction data, and we have created a second dataset with hub molecules removed, allowing a greater depth of information to be extracted from these graphs. We also present a detailed analysis of the various stages of the reconditioning process and the characteristics of the subgraphs resulting from the application of our software to the revised datasets. AVAILABILITY: Complete datasets and supplementary material may be downloaded from http://helix.ex.ac.uk/metabolism. The software for the creation of enzyme-centric graphs from reaction data is available on request from the authors.

The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models.

MOTIVATION: Molecular biotechnology now makes it possible to build elaborate systems models, but the systems biology community needs information standards if models are to be shared, evaluated and developed cooperatively. RESULTS: We summarize the Systems Biology Markup Language (SBML) Level 1, a free, open, XML-based format for representing biochemical reaction networks. SBML is a software-independent language for describing models common to research in many areas of computational biology, including cell signaling pathways, metabolic pathways, gene regulation, and others. AVAILABILITY: The specification of SBML Level 1 is freely available from http://www.sbml.org/

Simultaneous modelling of metabolic, genetic and product-interaction networks.

The creation of cell models from annotated genome information, as well as additional data from other databases, requires both a format and medium for its distribution. Standards are described for the representation of the data in the form of Document Type Definitions (DTDs) for XML files. Separate DTDs are detailed for genetic, metabolic and gene product-interaction networks, which can be used to hold information on individual subsystems, or which may be combined to create a whole cell DTD. In the execution of this work, a fifth DTD was also created for a metabolite thesaurus, which allows incorporation of metabolite synonyms and generic nomenclature data into the models. A gene-regulation classification scheme was also created, to facilitate incorporation of gene regulatory information in an efficient manner. The work is described with particular reference to the metabolic network of Escherichia coli, which contains 808 individual enzymes. The assignment of confidence levels to these data, through the use of Gene Ontology evidence codes, is highlighted. In silico investigations may now be performed using the mathematical simulation workbench, DBsolve, which incorporates the facility to introduce data directly from XML.

A historical perspective on gene/protein functional assignment.

Sequence determination and analysis began on proteins in the 1950s, with RNA starting about a decade later and DNA a similar period later still. Hence many of the concepts for function prediction were first developed by looking at amino acid sequences. Over time these methods have become much more sophisticated, allowing better discrimination of only weak similarities. The most recent developments concern an examination of contextual information, such as operon structure, metabolic reconstruction or co-expression profiles.

Comparison of functional annotation schemes for genomes.

In this paper we survey a number of functional classification schemes applicable to genomes. We present the concepts of depth, breadth and resolution as descriptors of the schemes' scope and architecture and compare selected classifications according to these criteria. We also generate a 'Combined Scheme' against which we map six classifications which we believe are representative of the range currently available. The mapping allows the generation of 'FuncWheels', which are graphical representations of hierarchical classification schemes. They are used to illustrate similarities and differences in functional space coverage. This survey highlights many issues related to the design and implementation of gene product functional classifications, which are discussed in the light of emerging 'second-generation' schemes.

Mathematical simulation and analysis of cellular metabolism and regulation.

MOTIVATION: A better understanding of the biological phenomena observed in cells requires the creation and analysis of mathematical models of cellular metabolism and physiology. The formulation and study of such models must also be simplified as far as possible to cope with the increasing complexity demanded and exponential accumulation of the metabolic reconstructions computed from sequenced genomes. RESULTS: A mathematical simulation workbench, DBsolve, has been developed to simplify the derivation and analysis of mathematical models. It combines: (i) derivation of large-scale mathematical models from metabolic reconstructions and other data sources; (ii) solving and parameter continuation of non-linear algebraic equations (NAEs), including metabolic control analysis; (iii) solving the non-linear stiff systems of ordinary differential equations (ODEs); (iv) bifurcation analysis of ODEs; (v) parameter fitting to experimental data or functional criteria based on constrained optimization. The workbench has been successfully used for dynamic metabolic modeling of some typical biochemical networks (Dolgacheva et al., Biochemistry (Moscow), 6, 1063-1068, 1996; Goldstein and Goryanin, Mol. Biol. (Moscow), 30, 976-983, 1996), including microbial glycolytic pathways, signal transduction pathways and receptor-ligand interactions. AVAILABILITY: DBsolve 5. 00 is freely available from http://websites.ntl.com/ approximately igor.goryanin. CONTACT: gzz78923@ggr.co.uk

Nucleoprotein architecture and ColE1 dimer resolution: a hypothesis.

Dimers of plasmid ColE1 are converted to monomers by site-specific recombination, a process that requires 240 bp of DNA (cer) and four host-encoded proteins (XerC, XerD, ArgR and PepA). Here, we propose structures for nucleoprotein complexes involved in cer-Xer recombination based upon existing knowledge of the structures of component proteins and computational analyses of protein structure and DNA curvature. We propose that, in the nucleoprotein complex at a single cer site, a PepA hexamer acts as an adaptor, connecting the heterodimeric recombinase (XerCD) to an ArgR hexamer. This provides a protein core around which the cer site wraps, its exact path being defined by strong sequence-specific interactions with ArgR and XerCD, weak interactions with PepA and sequence-dependent flexibility of cer. The initial association of single-site complexes (pairing) is proposed to occur via an ArgR-PepA interaction. Pairing between sites in a plasmid dimer is stabilized by DNA supercoiling and is followed by a structural isomerization to form a recombination-proficient synaptic complex. We propose that paired structures formed between sites in trans are too short-lived to permit synaptic complex formation. There is thus an energetic barrier to inappropriate recombination reactions. Our proposals are consistent with a wide range of experimental observations.

Automated management of gene discovery projects.

UNLABELLED: The paper reports on programs, organization and Hypertext Mark-up Language (HTML) presentation provided by the Hierarchical Automated Gene Identification System (HAGIS) designed for managing sequencing projects. AVAILABILITY: On request from the authors. CONTACT: st23646@ggr.co.uk

A new dynamic tool to perform assembly of expressed sequence tags (ESTs).

MOTIVATION: Expressed Sequence Tags (ESTs) are short single-pass DNA sequences obtained from either ends of cDNA clones. To exploit these sequences efficiently, a dynamic Web-tool has been developed which uses these data to perform fast virtual cloning of cDNAs. RESULTS: Starting with a query sequence, the user is able to identify related ESTs and extend the sequence of interest step by step, possibly to a full-length transcript. Graphical views of the clustering are used to monitor the progress of a particular 'cloning' project. Potential open reading frames are detected by positional base preference, and hyperlinks to other Worldwide Web sites allows the user to retrieve information relevant to each EST in a cluster (e.g. sequence traces, clone size, plate position). Apart from cDNA cloning, this tool also provides a mechanism for collating gene families and polymorphism sites.

Lewis X biosynthesis in Helicobacter pylori. Molecular cloning of an alpha(1,3)-fucosyltransferase gene.

The lipopolysaccharide of certain strains of Helicobacter pylori was recently shown to contain the Lewis X (Lex) trisaccharide (Galbeta-1, 4-(Fucalpha(1,3))-GlcNAc). Lex is an oncofetal antigen which appears on human gastric epithelium, and its mimicry by carbohydrate structures on the surface of H. pylori may play an important part in the interaction of this pathogen with its host. Potential roles for bacterial Lex in mucosal adhesion, immune evasion, and autoantibody induction have been proposed (Moran, A. P., Prendergast, M. M., and Appelmelk, B. J. (1996) FEMS Immunol. Med. Microbiol. 16, 105-115). In mammals, the final step of Lex biosynthesis is the alpha(1,3)-fucosylation of GlcNAc in a terminal Galbeta(1-->4)GlcNAc unit, and a corresponding GDP-fucose:N-acetylglucosaminyl alpha(1,3) fucosyltransferase (alpha(1,3)-Fuc-T) activity was recently discovered in H. pylori extracts. We used part of a human alpha(1, 3)-Fuc-T amino acid sequence to search an H. pylori genomic data base for related sequences. Using a probe based upon weakly matching data base sequences, we retrieved clones from a plasmid library of H. pylori DNA. DNA sequence analysis of the library clones revealed a gene which we have named fucT, encoding a protein with localized homology to the human alpha(1,3)-Fuc-Ts. We have demonstrated that fucT encodes an active Fuc-T enzyme by expressing the gene in Escherichia coli. The recombinant enzyme shows a strong preference for type 2 (e.g. LacNAc) over type 1 (e.g. lacto-N-biose) acceptors in vitro. Certain residues in a short segment of the H. pylori protein are completely conserved throughout the alpha(1,3)-Fuc-T family, defining an alpha(1,3)-Fuc-T motif which may be of use in identifying new fucosyltransferase genes.

Cloning and analysis of the first cry gene from Bacillus popilliae.

An 80-kDa parasporal crystal protein was detected in protein extracts of sporangia of Bacillus popilliae isolated from a diseased larva of the common cockchafer (Melolontha melolontha L.). Amino acid analysis of tryptic peptides revealed significant homology to the Cry2Aa endotoxins of Bacillus thuringiensis. The gene cryBP1 (cry18Aa1), which codes for the parasporal crystal protein, was found in a putative cry operon on the bacterial chromosome, which contains at least one further (smaller) open reading frame, orf1. The 706-amino-acid-long CryBP1 (Cry18Aa1) protein has a predicted molecular mass of 79 kDa and shows about 40% sequence identity to the Cry2 polypeptides of B. thuringiensis. In the light of published observations which suggest that the parasporal crystal proteins of B. popilliae are slightly toxic to their grub hosts, we propose the following survival strategy of B. popilliae. As an obligate pathogen of grubs, B. popilliae germinates in the gut of a grub and the parasporal crystal proteins are released and activated. The activated protein does not cause colloid osmotic lysis but instead damages the gut wall somehow to allow the vegetative cells to enter the hemolymph more easily. By becoming a parasite, B. popilliae can continue to proliferate efficiently while the living grub provides a food supply. This process is in contrast to that of B. thuringiensis, which rapidly kills the insect and is then limited to growth on the larval carcass.

pBR322 and other pSC101 derivatives have a gene fragment related to the mini-F plasmid resolvase.

Sequence comparisons have shown that part of pBR322 and other pSC101 derivatives contain a remnant resolvase gene closely related to that of plasmid mini-F. This observation illuminates part of the history of these plasmids.

Structural and functional studies of a synthetic peptide mimicking a proposed membrane inserting region of a Bacillus thuringiensis delta-endotoxin.

In order to study the mechanism of action of Bacillus thuringiensis delta-endotoxins, a synthetic 31-mer peptide corresponding to the sequence of a putative pore-forming segment of the CrylA(c) toxin was characterized structurally and functionally. The peptide maps onto the central helix (alpha 5) of the six-helix bundle of domain I of the crystal structure of the CryIIIA toxin. CD and NMR spectroscopic studies indicated that the peptide exists as an alpha-helix in methanol and a random coil in water. The peptide associated with liposomes at pH 4.7 and formed discrete, characterizable channels in planar lipid bilayers at low pHs. These channels had a conductance value of 60 picosiemens (pS). It is possible that this helix is a component of the transmembrane pore formed by B. thuringiensis delta-endotoxins in vivo.

Characterization of a Bacillus thuringiensis strain which is toxic to the housefly Musca domestica.

A Bacillus thuringiensis isolate has been discovered which is toxic to the common housefly (Musca domestica) as well as other Diptera and Lepidoptera. Crystal delta-endotoxins purified from this isolate killed 50% of Musca larvae at a concentration of 10.2 micrograms/ml, and beta-exotoxin was not detected. Sodium dodecyl polyacrylamide gel electrophoresis of the purified crystals revealed three protein species which were related to CryIA(b), CryIB and CryIIA toxins on the basis of immunoreactivity and amino-terminal sequence determination. Southern blot and DNA restriction analyses suggested that the strain has sequences related to one cryIA(b), one cryIIA, and two cryIIB genes.