A small nucleolar RNA requirement for site-specific ribose methylation of rRNA in Xenopus

Vertebrate cells contain a large number of small nucleolar RNA (snoRNA) species, the vast majority of which bind fibrillarin. Most of the fibrillarin-associated snoRNAs can form 10- to 21-nt duplexes with rRNA and are thought to guide 2′-O-methylation of selected nucleotides in rRNA. These include mammalian UHG (U22 host gene)-encoded U25–U31 snoRNAs. We have characterized two novel human snoRNA species, U62 and U63, which similarly exhibit 15- (with one interruption) and 12-nt complementarities and are therefore predicted to direct 2′-O-methylation of A590 in 18S and A4531 in 28S rRNA, respectively. To establish the function of antisense snoRNAs in vertebrates, we exploited the Xenopus oocyte system. Cloning of the Xenopus U25–U31 snoRNA genes indicated that they are encoded within multiple homologs of mammalian UHG. Depletion of U25 from the Xenopus oocyte abolished 2′-O-methylation of G1448 in 18S rRNA; methylation could be restored by injecting either the Xenopus or human U25 transcript into U25-depleted oocytes. Comparison of Xenopus and human U25 sequences revealed that only boxes C, D, and D′, as well as the 18S rRNA complement, were invariant, suggesting that they may be the only elements required for U25 snoRNA stability and function.

Nucleotide sequence of the Kaposi sarcoma-associated herpesvirus (HHV8)

The genome of the Kaposi sarcoma-associated herpesvirus (KSHV or HHV8) was mapped with cosmid and phage genomic libraries from the BC-1 cell line. Its nucleotide sequence was determined except for a 3-kb region at the right end of the genome that was refractory to cloning. The BC-1 KSHV genome consists of a 140.5-kb-long unique coding region flanked by multiple G+C-rich 801-bp terminal repeat sequences. A genomic duplication that apparently arose in the parental tumor is present in this cell culture-derived strain. At least 81 ORFs, including 66 with homology to herpesvirus saimiri ORFs, and 5 internal repeat regions are present in the long unique region. The virus encodes homologs to complement-binding proteins, three cytokines (two macrophage inflammatory proteins and interleukin 6), dihydrofolate reductase, bcl-2, interferon regulatory factors, interleukin 8 receptor, neural cell adhesion molecule-like adhesin, and a D-type cyclin, as well as viral structural and metabolic proteins. Terminal repeat analysis of virus DNA from a KS lesion suggests a monoclonal expansion of KSHV in the KS tumor.

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.