Horizontal gene transfer among genomes: The complexity hypothesis

Increasingly, studies of genes and genomes are indicating that considerable horizontal transfer has occurred between prokaryotes. Extensive horizontal transfer has occurred for operational genes (those involved in housekeeping), whereas informational genes (those involved in transcription, translation, and related processes) are seldomly horizontally transferred. Through phylogenetic analysis of six complete prokaryotic genomes and the identification of 312 sets of orthologous genes present in all six genomes, we tested two theories describing the temporal flow of horizontal transfer. We show that operational genes have been horizontally transferred continuously since the divergence of the prokaryotes, rather than having been exchanged in one, or a few, massive events that occurred early in the evolution of prokaryotes. In agreement with earlier studies, we found that differences in rates of evolution between operational and informational genes are minimal, suggesting that factors other than rate of evolution are responsible for the observed differences in horizontal transfer. We propose that a major factor in the more frequent horizontal transfer of operational genes is that informational genes are typically members of large, complex systems, whereas operational genes are not, thereby making horizontal transfer of informational gene products less probable (the complexity hypothesis).

Ectopic expression of the minimal whiE polyketide synthase generates a library of aromatic polyketides of diverse sizes and shapes

The single recombinant expressing the Streptomyces coelicolor minimal whiE (spore pigment) polyketide synthase (PKS) is uniquely capable of generating a large array of well more than 30 polyketides, many of which, so far, are novel to this recombinant. The characterized polyketides represent a diverse set of molecules that differ in size (chain length) and shape (cyclization pattern). This combinatorial biosynthetic library is, by far, the largest and most complex of its kind described to date and indicates that the minimal whiE PKS does not independently control polyketide chain length nor dictate the first cyclization event. Rather, the minimal PKS enzyme complex must rely on the stabilizing effects of additional subunits (i.e., the cyclase whiE-ORFVI) to ensure that the chain reaches the full 24 carbons and cyclizes correctly. This dramatic loss of control implies that the growing polyketide chain does not remain enzyme bound, resulting in the spontaneous cyclization of the methyl terminus. Among the six characterized dodecaketides, four different first-ring cyclization regiochemistries are represented, including C7/C12, C8/C13, C10/C15, and C13/C15. The dodecaketide TW93h possesses a unique 2,4-dioxaadamantane ring system and represents a new structural class of polyketides with no related structures isolated from natural or engineered organisms, thus supporting the claim that engineered biosynthesis is capable of producing novel chemotypes.

Multiple-complete-digest restriction fragment mapping: Generating sequence-ready maps for large-scale DNA sequencing

Multiple-complete-digest mapping is a DNA mapping technique based on complete-restriction-digest fingerprints of a set of clones that provides highly redundant coverage of the mapping target. The maps assembled from these fingerprints order both the clones and the restriction fragments. Maps are coordinated across three enzymes in the examples presented. Starting with yeast artificial chromosome contigs from the 7q31.3 and 7p14 regions of the human genome, we have produced cosmid-based maps spanning more than one million base pairs. Each yeast artificial chromosome is first subcloned into cosmids at a redundancy of ×15–30. Complete-digest fragments are electrophoresed on agarose gels, poststained, and imaged on a fluorescent scanner. Aberrant clones that are not representative of the underlying genome are rejected in the map construction process. Almost every restriction fragment is ordered, allowing selection of minimal tiling paths with clone-to-clone overlaps of only a few thousand base pairs. These maps demonstrate the practicality of applying the experimental and software-based steps in multiple-complete-digest mapping to a target of significant size and complexity. We present evidence that the maps are sufficiently accurate to validate both the clones selected for sequencing and the sequence assemblies obtained once these clones have been sequenced by a “shotgun” method.

Toward an understanding of the biochemical and pharmacological complexity of the saliva of a hematophagous sand fly Lutzomyia longipalpis

The saliva of blood-sucking arthropods contains powerful pharmacologically active substances and may be a vaccine target against some vector-borne diseases. Subtractive cloning combined with biochemical approaches was used to discover activities in the salivary glands of the hematophagous fly Lutzomyia longipalpis. Sequences of nine full-length cDNA clones were obtained, five of which are possibly associated with blood-meal acquisition, each having cDNA similarity to: (i) the bed bug Cimex lectularius apyrase, (ii) a 5′-nucleotidase/phosphodiesterase, (iii) a hyaluronidase, (iv) a protein containing a carbohydrate-recognition domain (CRD), and (v) a RGD-containing peptide with no significant matches to known proteins in the blast databases. Following these findings, we observed that the salivary apyrase activity of L. longipalpis is indeed similar to that of Cimex apyrase in its metal requirements. The predicted isoelectric point of the putative apyrase matches the value found for Lutzomyia salivary apyrase. A 5′-nucleotidase, as well as hyaluronidase activity, was found in the salivary glands, and the CRD-containing cDNA matches the N-terminal sequence of the HPLC-purified salivary anticlotting protein. A cDNA similar to α-amylase was discovered and salivary enzymatic activity demonstrated for the first time in a blood-sucking arthropod. Full-length clones were also found coding for three proteins of unknown function matching, respectively, the N-terminal sequence of an abundant salivary protein, having similarity to the CAP superfamily of proteins and the Drosophila yellow protein. Finally, two partial sequences are reported that match possible housekeeping genes. Subtractive cloning will considerably enhance efforts to unravel the salivary pharmacopeia of blood-sucking arthropods.

Transfected Drosophila cells as a probe for defining the minimal requirements for stimulating unprimed CD8+ T cells

Stimulation of naive T cells by antigen-presenting cells (APC) is thought to involve two qualitatively different signals: signal one results from T-cell receptor (TCR) recognition of antigenic peptides bound to major histocompatibility complex (MHC) molecules, whereas signal two reflects contact with one or more costimulatory molecules. The requirements for stimulating naive T cells were studied with MHC class I-restricted CD8+ T cells from a T-cell receptor transgenic line, with defined peptides as antigen and transfected Drosophila cells as APC. Three main findings are reported. First, stimulation of naive T cells via signal one alone (MHC plus peptide) was essentially nonimmunogenic; thus T cells cultured with peptides presented by MHC class I-transfected Drosophila APC lacking costimulatory molecules showed little or no change in their surface phenotype. Second, cotransfection of two costimulatory molecules, B7-1 and intercellular adhesion molecule 1 (ICAM-1), converted class I+ Drosophila cells to potent APC capable of inducing strong T-proliferative responses and cytokine (interleukin 2) production. Third, B7-1 and ICAM-1 acted synergistically, indicating that signal two is complex; synergy between B7-1 and ICAM-1 varied from moderate to extreme and was influenced by both the dose and affinity of the peptide used and the parameter of T-cell activation studied. Transfected Drosophila cells are thus a useful tool for examining the minimal APC requirements for naive T cells.

Fundamental patterns underlying gene expression profiles: Simplicity from complexity

Analysis of previously published sets of DNA microarray gene expression data by singular value decomposition has uncovered underlying patterns or “characteristic modes” in their temporal profiles. These patterns contribute unequally to the structure of the expression profiles. Moreover, the essential features of a given set of expression profiles are captured using just a small number of characteristic modes. This leads to the striking conclusion that the transcriptional response of a genome is orchestrated in a few fundamental patterns of gene expression change. These patterns are both simple and robust, dominating the alterations in expression of genes throughout the genome. Moreover, the characteristic modes of gene expression change in response to environmental perturbations are similar in such distant organisms as yeast and human cells. This analysis reveals simple regularities in the seemingly complex transcriptional transitions of diverse cells to new states, and these provide insights into the operation of the underlying genetic networks.

Minimal and optimal mechanisms for GroE-mediated protein folding

We have analyzed the effects of different components of the GroE chaperonin system on protein folding by using a nonpermissive substrate (i.e., one that has very low spontaneous refolding yield) for which rate data can be acquired. In the absence of GroES and nucleotides, the rate of GroEL-mediated refolding of heat- and DTT-denatured mitochondrial malate dehydrogenase was extremely low, but some three times higher than the spontaneous rate. This GroEL-mediated rate was increased 17-fold by saturating concentrations of ATP, 11-fold by ADP and GroES, and 465-fold by ATP and GroES. Optimal refolding activity was observed when the dissociation of GroES from the chaperonin complex was dramatically reduced. Although GroEL minichaperones were able to bind denatured mitochondrial malate dehydrogenase, they were ineffective in enhancing the refolding rate. The spectrum of mechanisms for GroE-mediated protein folding depends on the nature of the substrate. The minimal mechanism for permissive substrates (i.e., having significant yields of spontaneous refolding), requires only binding to the apical domain of GroEL. Slow folding rates of nonpermissive substrates are limited by the transitions between high- and low-affinity states of GroEL alone. The optimal mechanism, which requires holoGroEL, physiological amounts of GroES, and ATP hydrolysis, is necessary for the chaperonin-mediated folding of nonpermissive substrates at physiologically relevant rates under conditions in which retention of bound GroES prevents the premature release of aggregation-prone folding intermediates from the chaperonin complex. The different mechanisms are described in terms of the structural features of mini- and holo-chaperones.

Methylation of the minimal promoter of an embryonic globin gene silences transcription in primary erythroid cells

Methylation of cytosines in the dinucleotide CpG has been shown to suppress transcription of a number of tissue-specific genes, yet the precise mechanism is not fully understood. The vertebrate globin genes were among the first examples in which an inverse correlation was shown between CpG methylation and transcription. We studied the methylation pattern of the 235-bp ρ-globin gene promoter in genomic DNA from primary chicken erythroid cells using the sodium bisulfite conversion technique and found all CpGs in the promoter to be methylated in erythroid cells from adult chickens in which the ρ-globin gene is silent but unmethylated in 5-day (primitive) embryonic red cells in which the gene is transcribed. To elucidate further the mechanism of methylation-induced silencing, an expression construct consisting of 235 bp of 5′ promoter sequence of the ρ-globin gene along with a strong 5′ erythroid enhancer driving a chloramphenicol acetyltransferase reporter gene, ρ-CAT, was transfected into primary avian erythroid cells derived from 5-day embryos. Methylation of just the 235-bp ρ-globin gene promoter fragment at every CpG resulted in a 20- to 30-fold inhibition of transcription, and this effect was not overridden by the presence of potent erythroid-specific enhancers. The ability of the 235-bp ρ-globin gene promoter to bind to a DNA Methyl Cytosine binding Protein Complex (MeCPC) was tested in electrophoretic mobility shift assays utilizing primary avian erythroid cell nuclear extract. The results were that fully methylated but not unmethylated 235-bp ρ-globin gene promoter fragment could compete efficiently for MeCPC binding. These results are a direct demonstration that site-specific methylation of a globin gene promoter at the exact CpGs that are methylated in vivo can silence transcription in homologous primary erythroid cells. Further, these data implicate binding of MeCPC to the promoter in the mechanism of silencing.

Defining the minimal length of sequence homology required for selective gene isolation by TAR cloning

The transformation-associated recombination (TAR) cloning technique allows selective and accurate isolation of chromosomal regions and genes from complex genomes. The technique is based on in vivo recombination between genomic DNA and a linearized vector containing homologous sequences, or hooks, to the gene of interest. The recombination occurs during transformation of yeast spheroplasts that results in the generation of a yeast artificial chromosome (YAC) containing the gene of interest. To further enhance and refine the TAR cloning technology, we determined the minimal size of a specific hook required for gene isolation utilizing the Tg.AC mouse transgene as a targeted region. For this purpose a set of vectors containing a B1 repeat hook and a Tg.AC-specific hook of variable sizes (from 20 to 800 bp) was constructed and checked for efficiency of transgene isolation by a radial TAR cloning. When vectors with a specific hook that was ≥60 bp were utilized, ∼2% of transformants contained circular YACs with the Tg.AC transgene sequences. Efficiency of cloning dramatically decreased when the TAR vector contained a hook of 40 bp or less. Thus, the minimal length of a unique sequence required for gene isolation by TAR is ∼60 bp. No transgene-positive YAC clones were detected when an ARS element was incorporated into a vector, demonstrating that the absence of a yeast origin of replication in a vector is a prerequisite for efficient gene isolation by TAR cloning.

Genetic complexity of pathogen perception by plants: The example of Rcr3, a tomato gene required specifically by Cf-2

Genetic analysis of plant–pathogen interactions has demonstrated that resistance to infection is often determined by the interaction of dominant plant resistance (R) genes and dominant pathogen-encoded avirulence (Avr) genes. It was postulated that R genes encode receptors for Avr determinants. A large number of R genes and their cognate Avr genes have now been analyzed at the molecular level. R gene loci are extremely polymorphic, particularly in sequences encoding amino acids of the leucine-rich repeat motif. A major challenge is to determine how Avr perception by R proteins triggers the plant defense response. Mutational analysis has identified several genes required for the function of specific R proteins. Here we report the identification of Rcr3, a tomato gene required specifically for Cf-2-mediated resistance. We propose that Avr products interact with host proteins to promote disease, and that R proteins “guard” these host components and initiate Avr-dependent plant defense responses.

Photosystem II single crystals studied by EPR spectroscopy at 94 GHz: The tyrosine radical Y\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{D}^{{\bullet}}}}\end{equation*}\end{document}

Electron paramagnetic resonance (EPR) spectroscopy at 94 GHz is used to study the dark-stable tyrosine radical Y\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{D}^{{\bullet}}}}\end{equation*}\end{document} in single crystals of photosystem II core complexes (cc) isolated from the thermophilic cyanobacterium Synechococcus elongatus. These complexes contain at least 17 subunits, including the water-oxidizing complex (WOC), and 32 chlorophyll a molecules/PS II; they are active in light-induced electron transfer and water oxidation. The crystals belong to the orthorhombic space group P212121, with four PS II dimers per unit cell. High-frequency EPR is used for enhancing the sensitivity of experiments performed on small single crystals as well as for increasing the spectral resolution of the g tensor components and of the different crystal sites. Magnitude and orientation of the g tensor of Y\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{D}^{{\bullet}}}}\end{equation*}\end{document} and related information on several proton hyperfine tensors are deduced from analysis of angular-dependent EPR spectra. The precise orientation of tyrosine Y\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{D}^{{\bullet}}}}\end{equation*}\end{document} in PS II is obtained as a first step in the EPR characterization of paramagnetic species in these single crystals.

BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document} acts as a phosphate analog in proteins phosphorylated on aspartate: Structure of a BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document} complex with phosphoserine phosphatase

Protein phosphoaspartate bonds play a variety of roles. In response regulator proteins of two-component signal transduction systems, phosphorylation of an aspartate residue is coupled to a change from an inactive to an active conformation. In phosphatases and mutases of the haloacid dehalogenase (HAD) superfamily, phosphoaspartate serves as an intermediate in phosphotransfer reactions, and in P-type ATPases, also members of the HAD family, it serves in the conversion of chemical energy to ion gradients. In each case, lability of the phosphoaspartate linkage has hampered a detailed study of the phosphorylated form. For response regulators, this difficulty was recently overcome with a phosphate analog, BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document}, which yields persistent complexes with the active site aspartate of their receiver domains. We now extend the application of this analog to a HAD superfamily member by solving at 1.5-Å resolution the x-ray crystal structure of the complex of BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document} with phosphoserine phosphatase (PSP) from Methanococcus jannaschii. The structure is comparable to that of a phosphoenzyme intermediate: BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document} is bound to Asp-11 with the tetrahedral geometry of a phosphoryl group, is coordinated to Mg2+, and is bound to residues surrounding the active site that are conserved in the HAD superfamily. Comparison of the active sites of BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document}⋅PSP and BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document}⋅CeY, a receiver domain/response regulator, reveals striking similarities that provide insights into the function not only of PSP but also of P-type ATPases. Our results indicate that use of BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document} for structural studies of proteins that form phosphoaspartate linkages will extend well beyond response regulators.