The incompatible interaction between Phytophthora parasitica var. nicotianae race 0 and tobacco is suppressed in transgenic plants expressing antisense lipoxygenase sequences

Nicotiana tabacum 46-8 cultivar displays an incompatible interaction with race 0 of Phytophthora parasitica var. nicotianae (Ppn), a fungal pathogen of most tobacco cultivars. At the plant level, incompatibility is characterized by the induction of lipoxygenase (LOX, EC = 1.13.11.12) activity and localized hypersensitive cell death before defense gene activation. To evaluate the involvement of LOX in the onset of plant defense, tobacco 46-8 plants were genetically engineered using full-length or partial-length antisense (AS) tobacco LOX cDNA constructs. AS expression strongly reduced elicitor- and pathogen-induced LOX activity. Eight independent AS-LOX lines were selected and assayed for their response to Ppn. After root or stem inoculation with race 0, all AS-LOX lines but one displayed a compatible phenotype whereas control transformed plants, not containing the AS-LOX cassette, showed the typical incompatible reaction. The presence of the fungus in transgenic lines was demonstrated by PCR amplification of a Ppn-specific genomic sequence. A linear relationship was found between the extent of LOX suppression and the size of the lesion caused by the fungus. The AS-LOX plants also showed enhanced susceptibility toward the compatible fungus Rhizoctonia solani. The results demonstrate the strong involvement of LOX in the establishment of incompatibility in plant–microorganism interactions, consistent with its role in the defense of host plants.

Mode matches and their locations in the hydrophobic free energy sequences of peptide ligands and their receptor eigenfunctions

Patterns in sequences of amino acid hydrophobic free energies predict secondary structures in proteins. In protein folding, matches in hydrophobic free energy statistical wavelengths appear to contribute to selective aggregation of secondary structures in “hydrophobic zippers.” In a similar setting, the use of Fourier analysis to characterize the dominant statistical wavelengths of peptide ligands’ and receptor proteins’ hydrophobic modes to predict such matches has been limited by the aliasing and end effects of short peptide lengths, as well as the broad-band, mode multiplicity of many of their frequency (power) spectra. In addition, the sequence locations of the matching modes are lost in this transformation. We make new use of three techniques to address these difficulties: (i) eigenfunction construction from the linear decomposition of the lagged covariance matrices of the ligands and receptors as hydrophobic free energy sequences; (ii) maximum entropy, complex poles power spectra, which select the dominant modes of the hydrophobic free energy sequences or their eigenfunctions; and (iii) discrete, best bases, trigonometric wavelet transformations, which confirm the dominant spectral frequencies of the eigenfunctions and locate them as (absolute valued) moduli in the peptide or receptor sequence. The leading eigenfunction of the covariance matrix of a transmembrane receptor sequence locates the same transmembrane segments seen in n-block-averaged hydropathy plots while leaving the remaining hydrophobic modes unsmoothed and available for further analyses as secondary eigenfunctions. In these receptor eigenfunctions, we find a set of statistical wavelength matches between peptide ligands and their G-protein and tyrosine kinase coupled receptors, ranging across examples from 13.10 amino acids in acid fibroblast growth factor to 2.18 residues in corticotropin releasing factor. We find that the wavelet-located receptor modes in the extracellular loops are compatible with studies of receptor chimeric exchanges and point mutations. A nonbinding corticotropin-releasing factor receptor mutant is shown to have lost the signatory mode common to the normal receptor and its ligand. Hydrophobic free energy eigenfunctions and their transformations offer new quantitative physical homologies in database searches for peptide-receptor matches.

Backtracking leukemia to birth: Identification of clonotypic gene fusion sequences in neonatal blood spots

Epidemiological evidence has suggested that some pediatric leukemias may be initiated in utero and, for some pairs of identical twins with concordant leukemia, this possibility has been strongly endorsed by molecular studies of clonality. Direct evidence for a prenatal origin can only be derived by prospective or retrospective detection of leukemia-specific molecular abnormalities in fetal or newborn samples. We report a PCR-based method that has been developed to scrutinize neonatal blood spots (Guthrie cards) for the presence of numerically infrequent leukemic cells at birth in individuals who subsequently developed leukemia. We demonstrate that unique or clonotypic MLL-AF4 genomic fusion sequences are present and detectable in neonatal blood spots from individuals who were diagnosed with acute lymphoblastic leukemia at ages 5 months to 2 years and, therefore, have arisen during fetal hematopoiesis in utero. This result provides unequivocal evidence for a prenatal initiation of acute leukemia in young patients. The method should be applicable to other fusion genes in children with common subtypes of leukemia and will be of value in attempts to unravel the natural history and etiology of this major subtype of pediatric cancer.

Arrayed transposase-binding sequences on the ends of transposon Tn5090/Tn402

The transposon Tn5090/Tn402 encodes a 559 amino acid transposase, TniA, with a DDE motif. Gel mobility shifting and cleavage protection analysis with DNase I and hydroxyl radical probes revealed that TniA binds to multiple repeat sequences on either terminus of Tn5090/Tn402. Four of these TniA-binding 19mers occurred on the left-hand (t) end and two on the right-hand (i) end. Hydroxyl radical cleavage protection demonstrated the presence of 3–6 bp contact sequences on one face of the DNA helix. The binding pattern and organisation of repeats suggested parallels between Tn5090/Tn402 and Mu, which controls its transpositional activity in the assembly step of a higher order transpososome complex. The complex terminal structure and genes of transposase and nucleotide-binding proteins in tandem are hallmarks of the handful of Mu-like elements that are known to date.

Complementarity of the 16S rRNA penultimate stem with sequences downstream of the AUG destabilizes the plastid mRNAs

Escherichia coli mRNA translation is facilitated by sequences upstream and downstream of the initiation codon, called Shine–Dalgarno (SD) and downstream box (DB) sequences, respectively. In E.coli enhancing the complementarity between the DB sequences and the 16S rRNA penultimate stem resulted in increased protein accumulation without a significant affect on mRNA stability. The objective of this study was to test whether enhancing the complementarity of plastid mRNAs downstream of the AUG (downstream sequence or DS) with the 16S rRNA penultimate stem (anti-DS or ADS region) enhances protein accumulation. The test system was the tobacco plastid rRNA operon promoter fused with the E.coli phage T7 gene 10 (T7g10) 5′-untranslated region (5′-UTR) and DB region. Translation efficiency was tested by measuring neomycin phosphotransferase (NPTII) accumulation in tobacco chloroplasts. We report here that the phage T7g10 5′-UTR and DB region promotes accumulation of NPTII up to ∼16% of total soluble leaf protein (TSP). Enhanced mRNA stability and an improved NPTII yield (∼23% of TSP) was obtained from a construct in which the T7g10 5′-UTR was linked with the NPTII coding region via a NheI site. However, replacing the T7g10 DB region with the plastid DS sequence reduced NPTII and mRNA levels to 0.16 and 28%, respectively. Reduced NPTII accumulation is in part due to accelerated mRNA turnover.

Identification of sample-specific sequences in mammalian cDNA and genomic DNA by the novel ligation-mediated subtraction (Limes)

The representational difference analysis (RDA) and other subtraction techniques are used to enrich sample-specific sequences by elimination of ubiquitous sequences existing in both the sample of interest (tester) and the subtraction partner (driver). While applying the RDA to genomic DNA of cutaneous lymphoma cells in order to identify tumor relevant alterations, we predominantly isolated repetitive sequences and artificial repeat-mediated fusion products of otherwise independent PCR fragments (PCR hybrids). Since these products severely interfered with the isolation of tester-specific fragments, we developed a considerably more robust and efficient approach, termed ligation-mediated subtraction (Limes). In first applications of Limes, genomic sequences and/or transcripts of genes involved in the regulation of transcription, such as transforming growth factor β stimulated clone 22 related gene (TSC-22R), cell death and cytokine production (caspase-1) or antigen presentation (HLA class II sequences), were found to be completely absent in a cutaneous lymphoma line. On the assumption that mutations in tumor-relevant genes can affect their transcription pattern, a protocol was developed and successfully applied that allows the identification of such sequences. Due to these results, Limes may substitute/supplement other subtraction/comparison techniques such as RDA or DNA microarray techniques in a variety of different research fields.

Cis-acting requirements in flanking DNA for the programmed elimination of mse2.9: a common mechanism for deletion of internal eliminated sequences from the developing macronucleus of Tetrahymena thermophila

During macronuclear development in the ciliated protozoan Tetrahymena thermophila, extensive DNA deletions occur, eliminating thousands of internal eliminated sequences (IESs). Using an rDNA-based transformation assay we have analyzed the role during DNA deletion of DNA flanking mse2.9, an IES within the second intron of a gene encoding an as yet incompletely characterized protein. We establish that a cis-acting sequence for mse2.9 deletion acts at a distance to specify deletion boundaries. A complex sequence element necessary for efficient and accurate mse2.9 deletion is located in the region 47–81 bp from the right side of mse2.9. The ability of a variety of IES flanking sequences to rescue a processing deficient mse2.9 construct indicates that some cis-acting signal is shared among different IESs. In addition, the short intronic sequence that flanks mse2.9 is able to direct efficient and accurate processing. Despite no obvious sequence similarity between mse2.9 and other IESs, we suggest that a common mechanism is used to delete different families of IESs in Tetrahymena.

DNA-bound transcription factor complexes analysed by mass-spectrometry: binding of novel proteins to the human c-fos SRE and related sequences

Transcription factors control eukaryotic polymerase II function by influencing the recruitment of multiprotein complexes to promoters and their subsequent integrated function. The complexity of the functional ‘transcriptosome’ has necessitated biochemical fractionation and subsequent protein sequencing on a grand scale to identify individual components. As a consequence, much is now known of the basal transcription complex. In contrast, less is known about the complexes formed at distal promoter elements. The c-fos SRE, for example, is known to bind Serum Response Factor (SRF) and ternary complex factors such as Elk-1. Their interaction with other factors at the SRE is implied but, to date, none have been identified. Here we describe the use of mass-spectrometric sequencing to identify six proteins, SRF, Elk-1 and four novel proteins, captured on SRE duplexes linked to magnetic beads. This approach is generally applicable to the characterisation of nucleic acid-bound protein complexes and the post-translational modification of their components.

BAliBASE (Benchmark Alignment dataBASE): enhancements for repeats, transmembrane sequences and circular permutations

BAliBASE is specifically designed to serve as an evaluation resource to address all the problems encountered when aligning complete sequences. The database contains high quality, manually constructed multiple sequence alignments together with detailed annotations. The alignments are all based on three-dimensional structural superpositions, with the exception of the transmembrane sequences. The first release provided sets of reference alignments dealing with the problems of high variability, unequal repartition and large N/C-terminal extensions and internal insertions. Here we describe version 2.0 of the database, which incorporates three new reference sets of alignments containing structural repeats, trans­membrane sequences and circular permutations to evaluate the accuracy of detection/prediction and alignment of these complex sequences. BAliBASE can be viewed at the web site http://www-igbmc.u-strasbg.fr/BioInfo/BAliBASE2/index.html or can be downloaded from ftp://ftp-igbmc.u-strasbg.fr/pub/BAliBASE2/.

The TIGR Gene Indices: analysis of gene transcript sequences in highly sampled eukaryotic species

While genome sequencing projects are advancing rapidly, EST sequencing and analysis remains a primary research tool for the identification and categorization of gene sequences in a wide variety of species and an important resource for annotation of genomic sequence. The TIGR Gene Indices (http://www.tigr.org/tdb/tgi.shtml) are a collection of species-specific databases that use a highly refined protocol to analyze EST sequences in an attempt to identify the genes represented by that data and to provide additional information regarding those genes. Gene Indices are constructed by first clustering, then assembling EST and annotated gene sequences from GenBank for the targeted species. This process produces a set of unique, high-fidelity virtual transcripts, or Tentative Consensus (TC) sequences. The TC sequences can be used to provide putative genes with functional annotation, to link the transcripts to mapping and genomic sequence data, to provide links between orthologous and paralogous genes and as a resource for comparative sequence analysis.

rSNP_Guide, a database system for analysis of transcription factor binding to target sequences: application to SNPs and site-directed mutations

rSNP_Guide is a novel curated database system for analysis of transcription factor (TF) binding to target sequences in regulatory gene regions altered by mutations. It accumulates experimental data on naturally occurring site variants in regulatory gene regions and site-directed mutations. This database system also contains the web tools for SNP analysis, i.e., active applet applying weight matrices to predict the regulatory site candidates altered by a mutation. The current version of the rSNP_Guide is supplemented by six sub-databases: (i) rSNP_DB, on DNA–protein interaction caused by mutation; (ii) SYSTEM, on experimental systems; (iii) rSNP_BIB, on citations to original publications; (iv) SAMPLES, on experimentally identified sequences of known regulatory sites; (v) MATRIX, on weight matrices of known TF sites; (vi) rSNP_Report, on characteristic examples of successful rSNP_Tools implementation. These databases are useful for the analysis of natural SNPs and site-directed mutations. The databases are available through the Web, http://wwwmgs.bionet.nsc.ru/mgs/systems/rsnp/.

The Homeodomain Resource: sequences, structures, DNA binding sites and genomic information

The Homeodomain Resource is an annotated collection of non-redundant protein sequences, three-dimensional structures and genomic information for the homeodomain protein family. Release 3.0 contains 795 full-length homeodomain-containing sequences, 32 experimentally-derived structures and 143 homeo­box loci implicated in human genetic disorders. Entries are fully hyperlinked to facilitate easy retrieval of the original records from source databases. A simple search engine with a graphical user interface is provided to query the component databases and assemble customized data sets. A new feature for this release is the addition of DNA recognition sites for all human homeodomain proteins described in the literature. The Homeodomain Resource is freely available through the World Wide Web at http://genome.nhgri.nih.gov/homeodomain.

trEST, trGEN and Hits: access to databases of predicted protein sequences

High throughput genome (HTG) and expressed sequence tag (EST) sequences are currently the most abundant nucleotide sequence classes in the public database. The large volume, high degree of fragmentation and lack of gene structure annotations prevent efficient and effective searches of HTG and EST data for protein sequence homologies by standard search methods. Here, we briefly describe three newly developed resources that should make discovery of interesting genes in these sequence classes easier in the future, especially to biologists not having access to a powerful local bioinformatics environment. trEST and trGEN are regularly regenerated databases of hypothetical protein sequences predicted from EST and HTG sequences, respectively. Hits is a web-based data retrieval and analysis system providing access to precomputed matches between protein sequences (including sequences from trEST and trGEN) and patterns and profiles from Prosite and Pfam. The three resources can be accessed via the Hits home page (http://hits.isb-sib.ch).

Estimation of divergence times from multiprotein sequences for a few mammalian species and several distantly related organisms

When many protein sequences are available for estimating the time of divergence between two species, it is customary to estimate the time for each protein separately and then use the average for all proteins as the final estimate. However, it can be shown that this estimate generally has an upward bias, and that an unbiased estimate is obtained by using distances based on concatenated sequences. We have shown that two concatenation-based distances, i.e., average gamma distance weighted with sequence length (d2) and multiprotein gamma distance (d3), generally give more satisfactory results than other concatenation-based distances. Using these two distance measures for 104 protein sequences, we estimated the time of divergence between mice and rats to be approximately 33 million years ago. Similarly, the time of divergence between humans and rodents was estimated to be approximately 96 million years ago. We also investigated the dependency of time estimates on statistical methods and various assumptions made by using sequence data from eubacteria, protists, plants, fungi, and animals. Our best estimates of the times of divergence between eubacteria and eukaryotes, between protists and other eukaryotes, and between plants, fungi, and animals were 3, 1.7, and 1.3 billion years ago, respectively. However, estimates of ancient divergence times are subject to a substantial amount of error caused by uncertainty of the molecular clock, horizontal gene transfer, errors in sequence alignments, etc.

Herpes simplex virus DNA packaging sequences adopt novel structures that are specifically recognized by a component of the cleavage and packaging machinery

The product of the herpes simplex virus type 1 UL28 gene is essential for cleavage of concatemeric viral DNA into genome-length units and packaging of this DNA into viral procapsids. To address the role of UL28 in this process, purified UL28 protein was assayed for the ability to recognize conserved herpesvirus DNA packaging sequences. We report that DNA fragments containing the pac1 DNA packaging motif can be induced by heat treatment to adopt novel DNA conformations that migrate faster than the corresponding duplex in nondenaturing gels. Surprisingly, these novel DNA structures are high-affinity substrates for UL28 protein binding, whereas double-stranded DNA of identical sequence composition is not recognized by UL28 protein. We demonstrate that only one strand of the pac1 motif is responsible for the formation of novel DNA structures that are bound tightly and specifically by UL28 protein. To determine the relevance of the observed UL28 protein–pac1 interaction to the cleavage and packaging process, we have analyzed the binding affinity of UL28 protein for pac1 mutants previously shown to be deficient in cleavage and packaging in vivo. Each of the pac1 mutants exhibited a decrease in DNA binding by UL28 protein that correlated directly with the reported reduction in cleavage and packaging efficiency, thereby supporting a role for the UL28 protein–pac1 interaction in vivo. These data therefore suggest that the formation of novel DNA structures by the pac1 motif confers added specificity on recognition of DNA packaging sequences by the UL28-encoded component of the herpesvirus cleavage and packaging machinery.