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.

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.

Identification of a cryptic gene associated with an insertion sequence not previously identified in Bacillus thuringiensis.

After screening several Bt strains with a cryII toxin probe, clones from two strains were found to contain a cryptic cryIIB gene associated with an insertion sequence element belonging to the IS2/IS3 family. The lack of expression of this gene appears to result from mutation of the upstream orf2 gene which has been shown to be necessary for cryII expression.

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.

Models for the structure and function of the Bacillus thuringiensis delta-endotoxins determined by compilational analysis.

The protein delta-endotoxins of Bacillus thuringiensis are a commercially and environmentally important class of highly specific insecticides. From an alignment of their sequences, certain structural and functional domains can be inferred which may shed light on the mode of action of these toxins.

Secondary structure analysis identifies a putative mouse protein demonstrating similarity to the repeat units found in CDC4, the G protein beta subunits and related proteins.

The predicted protein product of an anonymous clone isolated from a cDNA library prepared from 12 day post coitum (p.c) embryonic mouse heart tissue demonstrated the same segmental repeats previously identified in the cell division control protein, CDC4 and the G protein beta 1 subunit. A search of the protein database subsequently identified three other classes of protein containing the repeat. Secondary structure analyses performed on the repeat sequences revealed a high degree of conservation suggesting that the repeat motif performs a specific function in a diverse range of proteins.

The elucidation of protein function from its amino acid sequence.

This review gives an outline of how computers may be used to determine the function of a protein, when only its primary structure is known. The current programming methods are outlined in general terms before their detailed application is discussed, and the common ways of predicting protein structure are also introduced. Identification is usually by database searching and sequence alignment, though a collection of motifs relating sequence to function are also described.

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.

The elucidation of protein function by sequence motif analysis.

Protein sequence motifs are acquiring increasing prominence in the area of sequence analysis. This review describes the current methods of their construction and their use in the determination of protein function, and offers guidelines on interpreting data obtained. An appendix is attached which refers to 200 motifs of various kinds.

The herpes simplex virus type 2 equivalent of the herpes simplex virus type 1 US7 gene and its flanking sequences.

Nucleotide sequencing studies (D. J. McGeoch, A. Dolan, S. Donald, and F. Rixon, 1985, J. Mol. Biol. 181, 1-14) have indicated that herpes simplex virus type 1 (HSV-1) has a coding sequence, referred to as US7, between the genes for the glycoproteins D and E (gD and gE). Northern blot analysis and nucleotide sequencing have been carried out to show that the type 2 virus (HSV-2) has an equivalent to the US7 gene. A comparison with the HSV-1 sequence has revealed some surprising similarities and differences. At the nucleotide level, HSV-2 has inserted a large sequence into the gE promoter, retained a large palindrome present in the coding sequence but not some tandem repeats, and deleted a region beside those repeats. At the amino acid level, the putative transmembrane sequence has been remarkably well conserved, and hydrophobic moment analysis indicates that it could be interacting with polar species within the plane of the membrane. Immediately after the deletion in the HSV-2 sequence, there is an N-glycosylation signal, and HSV-2 has one more such signal than HSV-1. The longest conserved sequence at the nucleotide level codes for a region of polypeptide that is strongly predicted to fold into alpha-helix. Implications of these analyses to the structure and possible function of these molecules are discussed.

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.

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.

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

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.

The construction of Bacillus thuringiensis strains expressing novel entomocidal delta-endotoxin combinations.

Using our recently reported method of electroporation to transform Bacillus thuringiensis [Bone & Ellar (1989) FEMS Microbiol. Lett. 58, 171-178], cloned B. thuringiensis entomocidal delta-endotoxin genes have been introduced into several native B. thuringiensis strains. In many cases the resulting transformants expressed both their native toxins and the cloned toxin, producing strains with broader toxicity spectra. The introduction of the var. tenebrionis toxin gene into B. thuringiensis var. israelensis resulted in a strain with activity against Pieris brassicae (cabbage white butterfly), an activity which neither parent strain possesses. We discuss further the possibility of synergism and also the problems associated with introducing cloned DNA by this method.