Identification of a novel selD homolog from Eukaryotes, Bacteria, and Archaea: Is there an autoregulatory mechanism in selenocysteine metabolism?

Escherichia coli selenophosphate synthetase (SPS, the selD gene product) catalyzes the production of monoselenophosphate, the selenium donor compound required for synthesis of selenocysteine (Sec) and seleno-tRNAs. We report the molecular cloning of human and mouse homologs of the selD gene, designated Sps2, which contains an in-frame TGA codon at a site corresponding to the enzyme’s putative active site. These sequences allow the identification of selD gene homologs in the genomes of the bacterium Haemophilus influenzae and the archaeon Methanococcus jannaschii, which had been previously misinterpreted due to their in-frame TGA codon. Sps2 mRNA levels are elevated in organs previously implicated in the synthesis of selenoproteins and in active sites of blood cell development. In addition, we show that Sps2 mRNA is up-regulated upon activation of T lymphocytes and have mapped the Sps2 gene to mouse chromosome 7. Using the mouse gene isolated from the hematopoietic cell line FDCPmixA4, we devised a construct for protein expression that results in the insertion of a FLAG tag sequence at the N terminus of the SPS2 protein. This strategy allowed us to document the readthrough of the in-frame TGA codon and the incorporation of 75Se into SPS2. These results suggest the existence of an autoregulatory mechanism involving the incorporation of Sec into SPS2 that might be relevant to blood cell biology. This mechanism is likely to have been present in ancient life forms and conserved in a variety of living organisms from all domains of life.

DAF-16 recruits the CREB-binding protein coactivator complex to the insulin-like growth factor binding protein 1 promoter in HepG2 cells

Insulin negatively regulates expression of the insulin-like growth factor binding protein 1 (IGFBP-1) gene by means of an insulin-responsive element (IRE) that also contributes to glucocorticoid stimulation of this gene. We find that the Caenorhabditis elegans protein DAF-16 binds the IGFBP-1⋅IRE with specificity similar to that of the forkhead (FKH) factor(s) that act both to enhance glucocorticoid responsiveness and to mediate the negative effect of insulin at this site. In HepG2 cells, DAF-16 and its mammalian homologs, FKHR, FKHRL1, and AFX, activate transcription through the IGFBP-1⋅IRE; this effect is inhibited by the viral oncoprotein E1A, but not by mutants of E1A that fail to interact with the coactivator p300/CREB-binding protein (CBP). We show that DAF-16 and FKHR can interact with both the KIX and E1A/SRC interaction domains of p300/CBP, as well as the steroid receptor coactivator (SRC). A C-terminal deletion mutant of DAF-16 that is nonfunctional in C. elegans fails to bind the KIX domain of CBP, fails to activate transcription through the IGFBP-1⋅IRE, and inhibits activation of the IGFBP-1 promoter by glucocorticoids. Thus, the interaction of DAF-16 homologs with the KIX domain of CBP is essential to basal and glucocorticoid-stimulated transactivation. Although AFX interacts with the KIX domain of CBP, it does not interact with SRC and does not respond to glucocorticoids or insulin. Thus, we conclude that DAF-16 and FKHR act as accessory factors to the glucocorticoid response, by recruiting the p300/CBP/SRC coactivator complex to an FKH factor site in the IGFBP-1 promoter, which allows the cell to integrate the effects of glucocorticoids and insulin on genes that carry this site.

Toll-related receptors and the control of antimicrobial peptide expression in Drosophila

Insects defend themselves against infectious microorganisms by synthesizing potent antimicrobial peptides. Drosophila has appeared in recent years as a favorable model to study this innate host defense. A genetic analysis of the regulation of the antifungal peptide drosomycin has demonstrated a key role for the transmembrane receptor Toll, which prompted the search for mammalian homologs. Two of these, Toll-like receptor (TLR)2 and TLR4, recently were shown to play a critical role in innate immunity against bacteria. Here we describe six additional Toll-related genes (Toll-3 to Toll-8) in Drosophila in addition to 18-wheeler. Two of these genes, Toll-3 and Toll-4, are expressed at a low level. Toll-6, -7, and -8, on the other hand, are expressed at high levels during embryogenesis and molting, suggesting that, like Toll and 18w, they perform developmental functions. Finally, Toll-5 is expressed only in larvae and adults. By using chimeric constructs, we have tested the capacity of the signaling Toll/IL-1R homology domains of these receptors to activate antimicrobial peptide promoters and found that only Toll and Toll-5 can activate the drosomycin promoter in transfected cells, thus demonstrating specificity at the level of the Toll/IL-1R homology domain. In contrast, none of these constructs activated antibacterial peptide promoters, suggesting that Toll-related receptors are not involved in the regulation of antibacterial peptide expression. This result was independently confirmed by the demonstration that a dominant-negative version of the kinase Pelle can block induction of drosomycin by the cytokine Spaetzle, but does not affect induction of the antibacterial peptide attacin by lipopolysaccharide.

Tid1/Rdh54 promotes colocalization of Rad51 and Dmc1 during meiotic recombination

Two RecA homologs, Rad51 and Dmc1, assemble as cytologically visible complexes (foci) at the same sites on meiotic chromosomes. Time course analysis confirms that co-foci appear and disappear as the single predominant form. A large fraction of co-foci are eliminated in a red1 mutant, which is expected as a characteristic of the interhomolog-specific recombination pathway. Previous studies suggested that normal Dmc1 loading depends on Rad51. We show here that a mutation in TID1/RDH54, encoding a RAD54 homolog, reduces Rad51-Dmc1 colocalization relative to WT. A rad54 mutation, in contrast, has relatively little effect on RecA homolog foci except when strains also contain a tid1/rdh54 mutation. The role of Tid1/Rdh54 in coordinating RecA homolog assembly may be very direct, because Tid1/Rdh54 is known to physically bind both Dmc1 and Rad51. Also, Dmc1 foci appear early in a tid1/rdh54 mutant. Thus, Tid1 may normally act with Rad51 to promote ordered RecA homolog assembly by blocking Dmc1 until Rad51 is present. Finally, whereas double-staining foci predominate in WT nuclei, a subset of nuclei with expanded chromatin exhibit individual Rad51 and Dmc1 foci side-by-side, suggesting that a Rad51 homo-oligomer and a Dmc1 homo-oligomer assemble next to one another at the site of a single double-strand break (DSB) recombination intermediate.

Assessing sequence comparison methods with reliable structurally identified distant evolutionary relationships

Pairwise sequence comparison methods have been assessed using proteins whose relationships are known reliably from their structures and functions, as described in the scop database [Murzin, A. G., Brenner, S. E., Hubbard, T. & Chothia C. (1995) J. Mol. Biol. 247, 536–540]. The evaluation tested the programs blast [Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. (1990). J. Mol. Biol. 215, 403–410], wu-blast2 [Altschul, S. F. & Gish, W. (1996) Methods Enzymol. 266, 460–480], fasta [Pearson, W. R. & Lipman, D. J. (1988) Proc. Natl. Acad. Sci. USA 85, 2444–2448], and ssearch [Smith, T. F. & Waterman, M. S. (1981) J. Mol. Biol. 147, 195–197] and their scoring schemes. The error rate of all algorithms is greatly reduced by using statistical scores to evaluate matches rather than percentage identity or raw scores. The E-value statistical scores of ssearch and fasta are reliable: the number of false positives found in our tests agrees well with the scores reported. However, the P-values reported by blast and wu-blast2 exaggerate significance by orders of magnitude. ssearch, fasta ktup = 1, and wu-blast2 perform best, and they are capable of detecting almost all relationships between proteins whose sequence identities are >30%. For more distantly related proteins, they do much less well; only one-half of the relationships between proteins with 20–30% identity are found. Because many homologs have low sequence similarity, most distant relationships cannot be detected by any pairwise comparison method; however, those which are identified may be used with confidence.

Genomic evidence for two functionally distinct gene classes

Analyses of complete genomes indicate that a massive prokaryotic gene transfer (or transfers) preceded the formation of the eukaryotic cell. In comparisons of the entire set of Methanococcus jannaschii genes with their orthologs from Escherichia coli, Synechocystis 6803, and the yeast Saccharomyces cerevisiae, it is shown that prokaryotic genomes consist of two different groups of genes. The deeper, diverging informational lineage codes for genes which function in translation, transcription, and replication, and also includes GTPases, vacuolar ATPase homologs, and most tRNA synthetases. The more recently diverging operational lineage codes for amino acid synthesis, the biosynthesis of cofactors, the cell envelope, energy metabolism, intermediary metabolism, fatty acid and phospholipid biosynthesis, nucleotide biosynthesis, and regulatory functions. In eukaryotes, the informational genes are most closely related to those of Methanococcus, whereas the majority of operational genes are most closely related to those of Escherichia, but some are closest to Methanococcus or to Synechocystis.

Genetic footprinting with mariner-based transposition in Pseudomonas aeruginosa

The complete DNA sequence of Pseudomonas aeruginosa provides an opportunity to apply functional genomics to a major human pathogen. A comparative genomics approach combined with genetic footprinting was used as a strategy to identify genes required for viability in P. aeruginosa. Use of a highly efficient in vivo mariner transposition system in P. aeruginosa facilitated the analysis of candidate genes of this class. We have developed a rapid and efficient allelic exchange system by using the I-SceI homing endonuclease in conjunction with in vitro mariner mutagenesis to generate mutants within targeted regions of the P. aeruginosa chromosome for genetic footprinting analyses. This technique for generating transposon insertion mutants should be widely applicable to other organisms that are not naturally transformable or may lack well developed in vivo transposition systems. We tested this system with three genes in P. aeruginosa that have putative essential homologs in Haemophilus influenzae. We show that one of three H. influenzae essential gene homologs is needed for growth in P. aeruginosa, validating the practicality of this comparative genomics strategy to identify essential genes in P. aeruginosa.

Molecular characterization and assembly of the needle complex of the Salmonella typhimurium type III protein secretion system

Many bacterial pathogens of plants and animals have evolved a specialized protein-secretion system termed type III to deliver bacterial proteins into host cells. These proteins stimulate or interfere with host cellular functions for the pathogen's benefit. The Salmonella typhimurium pathogenicity island 1 encodes one of these systems that mediates this bacterium's ability to enter nonphagocytic cells. Several components of this type III secretion system are organized in a supramolecular structure termed the needle complex. This structure is made of discrete substructures including a base that spans both membranes and a needle-like projection that extends outward from the bacterial surface. We demonstrate here that the type III secretion export apparatus is required for the assembly of the needle substructure but is dispensable for the assembly of the base. We show that the length of the needle segment is determined by the type III secretion associated protein InvJ. We report that InvG, PrgH, and PrgK constitute the base and that PrgI is the main component of the needle of the type III secretion complex. PrgI homologs are present in type III secretion systems from bacteria pathogenic for animals but are absent from bacteria pathogenic for plants. We hypothesize that the needle component may establish the specificity of type III secretion systems in delivering proteins into either plant or animal cells.

Expression of homeobox genes shows chelicerate arthropods retain their deutocerebral segment

Expression patterns of six homeobox containing genes in a model chelicerate, the oribatid mite Archegozetes longisetosus, were examined to establish homology of chelicerate and insect head segments and to investigate claims that the chelicerate deutocerebral segment has been reduced or lost. engrailed (en) expression, which has been used to demonstrate the presence of segments in insects, fails to demonstrate a reduced deutocerebral segment. Expression patterns of the chelicerate homologs of the Drosophila genes Antennapedia (Antp), Sex combs reduced (Scr), Deformed (Dfd), proboscipedia (pb), and orthodenticle (otd) confirm direct correspondence of head segments. The chelicerate deutocerebral segment has not been reduced or lost. We make further inferences concerning the evolution of heads and Hox genes in arthropods.

Independent deletions of a pathogen-resistance gene in Brassica and Arabidopsis

Plant disease resistance (R) genes confer race-specific resistance to pathogens and are genetically defined on the basis of intra-specific functional polymorphism. Little is known about the evolutionary mechanisms that generate this polymorphism. Most R loci examined to date contain alternate alleles and/or linked homologs even in disease-susceptible plant genotypes. In contrast, the resistance to Pseudomonas syringae pathovar maculicola (RPM1) bacterial resistance gene is completely absent (rpm1-null) in 5/5 Arabidopsis thaliana accessions that lack RPM1 function. The rpm1-null locus contains a 98-bp segment of unknown origin in place of the RPM1 gene. We undertook comparative mapping of RPM1 and flanking genes in Brassica napus to determine the ancestral state of the RPM1 locus. We cloned two B. napus RPM1 homologs encoding hypothetical proteins with ≈81% amino acid identity to Arabidopsis RPM1. Collinearity of genes flanking RPM1 is conserved between B. napus and Arabidopsis. Surprisingly, we found four additional B. napus loci in which the flanking marker synteny is maintained but RPM1 is absent. These B. napus rpm1-null loci have no detectable nucleotide similarity to the Arabidopsis rpm1-null allele. We conclude that RPM1 evolved before the divergence of the Brassicaceae and has been deleted independently in the Brassica and Arabidopsis lineages. These results suggest that functional polymorphism at R gene loci can arise from gene deletions.

Molecular evolution of two vertebrate aryl hydrocarbon (dioxin) receptors (AHR1 and AHR2) and the PAS family

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor through which halogenated aromatic hydrocarbons such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) cause altered gene expression and toxicity. The AHR belongs to the basic helix–loop–helix/Per-ARNT-Sim (bHLH-PAS) family of transcriptional regulatory proteins, whose members play key roles in development, circadian rhythmicity, and environmental homeostasis; however, the normal cellular function of the AHR is not yet known. As part of a phylogenetic approach to understanding the function and evolutionary origin of the AHR, we sequenced the PAS homology domain of AHRs from several species of early vertebrates and performed phylogenetic analyses of these AHR amino acid sequences in relation to mammalian AHRs and 24 other members of the PAS family. AHR sequences were identified in a teleost (the killifish Fundulus heteroclitus), two elasmobranch species (the skate Raja erinacea and the dogfish Mustelus canis), and a jawless fish (the lamprey Petromyzon marinus). Two putative AHR genes, designated AHR1 and AHR2, were found both in Fundulus and Mustelus. Phylogenetic analyses indicate that the AHR2 genes in these two species are orthologous, suggesting that an AHR gene duplication occurred early in vertebrate evolution and that multiple AHR genes may be present in other vertebrates. Database searches and phylogenetic analyses identified four putative PAS proteins in the nematode Caenorhabditis elegans, including possible AHR and ARNT homologs. Phylogenetic analysis of the PAS gene family reveals distinct clades containing both invertebrate and vertebrate PAS family members; the latter include paralogous sequences that we propose have arisen by gene duplication early in vertebrate evolution. Overall, our analyses indicate that the AHR is a phylogenetically ancient protein present in all living vertebrate groups (with a possible invertebrate homolog), thus providing an evolutionary perspective to the study of dioxin toxicity and AHR function.

Rapid evolution of the DNA-binding site in LAGLIDADG homing endonucleases

Sequence analysis of chloroplast and mitochondrial large subunit rRNA genes from over 75 green algae disclosed 28 new group I intron-encoded proteins carrying a single LAGLIDADG motif. These putative homing endonucleases form four subfamilies of homologous enzymes, with the members of each subfamily being encoded by introns sharing the same insertion site. We showed that four divergent endonucleases from the I-CreI subfamily cleave the same DNA substrates. Mapping of the 66 amino acids that are conserved among the members of this subfamily on the 3-dimensional structure of I-CreI bound to its recognition sequence revealed that these residues participate in protein folding, homodimerization, DNA recognition and catalysis. Surprisingly, only seven of the 21 I-CreI amino acids interacting with DNA are conserved, suggesting that I-CreI and its homologs use different subsets of residues to recognize the same DNA sequence. Our sequence comparison of all 45 single-LAGLIDADG proteins identified so far suggests that these proteins share related structures and that there is a weak pressure in each subfamily to maintain identical protein–DNA contacts. The high sequence variability we observed in the DNA-binding site of homologous LAGLIDADG endonucleases provides insight into how these proteins evolve new DNA specificity.

Efficient recognition of substrates and substrate analogs by the adenine glycosylase MutY requires the C-terminal domain

The Escherichia coli DNA repair enzyme MutY plays an important role in the prevention of DNA mutations by removing misincorporated adenine residues from 7,8-dihydro-8-oxo-2′-deoxyguanosine:2′-deoxyadenosine (OG:A) mispairs. The N-terminal domain of MutY (Stop 225, Met1–Lys225) has a sequence and structure that is characteristic of a superfamily of base excision repair glycosylases; however, MutY and its homologs contain a unique C-terminal domain. Previous studies have shown that the C-terminal domain confers specificity for OG:A substrates over G:A substrates and exhibits homology to the d(OG)TPase MutT, suggesting a role in OG recognition. In order to provide additional information on the importance of the C-terminal domain in damage recognition, we have investigated the kinetic properties of a form lacking this domain (Stop 225) under multiple- and single-turnover conditions. In addition, the interaction of Stop 225 with a series of non-cleavable substrate and product analogs was evaluated using gel retardation assays and footprinting experiments. Under multiple-turnover conditions Stop 225 exhibits biphasic kinetic behavior with both OG:A and G:A substrates, likely due to rate-limiting DNA product release. However, the rate of turnover of Stop 225 was increased 2-fold with OG:A substrates compared to the wild-type enzyme. In contrast, the intrinsic rate for adenine removal by Stop 225 from both G:A and OG:A substrates is significantly reduced (10- to 25-fold) compared to the wild-type. The affinity of Stop 225 for substrate analogs was dramatically reduced, as was the ability to discriminate between substrate analogs paired with OG over G. Interestingly, similar hydroxyl radical and DMS footprinting patterns are observed for Stop 225 and wild-type MutY bound to DNA duplexes containing OG opposite an abasic site mimic or a non-hydrogen bonding A analog, suggesting that similar regions of the DNA are contacted by both enzyme forms. Importantly, Stop 225 has a reduced ability to prevent DNA mutations in vivo. This implies that the reduced adenine glycosylase activity translates to a reduced capacity of Stop 225 to prevent DNA mutations in vivo.

Structure of human DNMT2, an enigmatic DNA methyltransferase homolog that displays denaturant-resistant binding to DNA

DNMT2 is a human protein that displays strong sequence similarities to DNA (cytosine-5)-methyltransferases (m5C MTases) of both prokaryotes and eukaryotes. DNMT2 contains all 10 sequence motifs that are conserved among m5C MTases, including the consensus S-adenosyl-l-methionine-binding motifs and the active site ProCys dipeptide. DNMT2 has close homologs in plants, insects and Schizosaccharomyces pombe, but no related sequence can be found in the genomes of Saccharomyces cerevisiae or Caenorhabditis elegans. The crystal structure of a deletion mutant of DNMT2 complexed with S-adenosyl-l-homocysteine (AdoHcy) has been determined at 1.8 Å resolution. The structure of the large domain that contains the sequence motifs involved in catalysis is remarkably similar to that of M.HhaI, a confirmed bacterial m5C MTase, and the smaller target recognition domains of DNMT2 and M.HhaI are also closely related in overall structure. The small domain of DNMT2 contains three short helices that are not present in M.HhaI. DNMT2 binds AdoHcy in the same conformation as confirmed m5C MTases and, while DNMT2 shares all sequence and structural features with m5C MTases, it has failed to demonstrate detectable transmethylase activity. We show here that homologs of DNMT2, which are present in some organisms that are not known to methylate their genomes, contain a specific target-recognizing sequence motif including an invariant CysPheThr tripeptide. DNMT2 binds DNA to form a denaturant-resistant complex in vitro. While the biological function of DNMT2 is not yet known, the strong binding to DNA suggests that DNMT2 may mark specific sequences in the genome by binding to DNA through the specific target-recognizing motif.

Isolation of segments of homologous genes with only one conserved amino acid region via PCR

We present a method which allows the isolation of fragments from genes coding for homologous proteins via PCR when only one block of conserved amino acids is available. Sets of degenerated primers are defined by reverse translation of the conserved amino acids such that each set contains not more than 128 different sequences. The second primer binding site is provided by a special cassette that is designed such that it does not allow binding of the second primer prior to being copied by DNA synthesis. The cassette is ligated to partially-digested chromosomal DNA. The second primer is biotinylated to allow elimination of PCR products carrying degenerated primers on both sides via streptavidin binding. Fragments obtained after amplification and enrichment are cloned and sequenced. The feasibility of this method was demonstrated in a model experiment, where degenerated primers were deduced from six conserved amino acids within the family of homologs to the Escherichia coli Vsr protein.