A comprehensive BAC resource

The Human Genome Project has generated extensive map and sequence data for a large number of Bacterial Artificial Chromosome (BAC) clones. In order to maximize the efficient use of the data and to minimize the redundant work for the research community, The Institute for Genomic Research (TIGR) comprehensive BAC resource (cBACr) (http://www.tigr.org/tdb/BacResource/BAC_resource_intro.html) was built as an expansion of the TIGR human BAC ends database. This resource collects, integrates and reports the information on library, maps, sequence, annotation and functions for each human and mouse BAC. The current database contains 635 016 human BACs and 265 617 mouse BACs that were characterized by various approaches, among which 22 705 human clones and 1000 mouse clones have sequence and annotation data.

ICB database: the gyrB database for identification and classification of bacteria

The Identification and Classification of Bacteria (ICB) database (http:/www.mbio.co.jp/icb) contains currently available information about the DNA gyrase subunit B (gyrB) gene in bacteria. The database is designed to provide the scientific community with a reference point for using gyrB as an evolutionary and taxonomic marker. Nucleic and amino acid sequence data are currently available for over 850 strains, along with alignments at several different taxonomic levels and an exhaustive review of primer selection and background information.

REBASE—restriction enzymes and methylases

REBASE contains comprehensive information about restriction enzymes, DNA methylases and related proteins such as nicking enzymes, specificity subunits and control proteins. It contains published and unpublished references, recognition and cleavage sites, isoschizomers, commercial availability, methy­lation sensitivity, crystal data and sequence data. Homing endonucleases are also included. Most recently, extensive information about the methy­lation sensitivity of restriction enzymes has been added and a new feature contains complete analyses of the putative restriction systems in the sequenced bacterial and archaeal genomes. The data is distributed via email, ftp (ftp.neb.com) and the Web (http://rebase.neb.com).

Toward a plant genomics initiative: Thoughts on the value of cross-species and cross-genera comparisons in the grasses

Comparative genomics offers unparalleled opportunities to integrate historically distinct disciplines, to link disparate biological kingdoms, and to bridge basic and applied science. Cross-species, cross-genera, and cross-kingdom comparisons are proving key to understanding how genes are structured, how gene structure relates to gene function, and how changes in DNA have given rise to the biological diversity on the planet. The application of genomics to the study of crop species offers special opportunities for innovative approaches for combining sequence information with the vast reservoirs of historical information associated with crops and their evolution. The grasses provide a particularly well developed system for the development of tools to facilitate comparative genetic interpretation among members of a diverse and evolutionarily successful family. Rice provides advantages for genomic sequencing because of its small genome and its diploid nature, whereas each of the other grasses provides complementary genetic information that will help extract meaning from the sequence data. Because of the importance of the cereals to the human food chain, developments in this area can lead directly to opportunities for improving the health and productivity of our food systems and for promoting the sustainable use of natural resources.

Maize as a model for the evolution of plant nuclear genomes

The maize genome is replete with chromosomal duplications and repetitive DNA. The duplications resulted from an ancient polyploid event that occurred over 11 million years ago. Based on DNA sequence data, the polyploid event occurred after the divergence between sorghum and maize, and hence the polyploid event explains some of the difference in DNA content between these two species. Genomic rearrangement and diploidization followed the polyploid event. Most of the repetitive DNA in the maize genome is retrotransposable elements, and they comprise 50% of the genome. Retrotransposon multiplication has been relatively recent—within the last 5–6 million years—suggesting that the proliferation of retrotransposons has also contributed to differences in DNA content between sorghum and maize. There are still unanswered questions about repetitive DNA, including the distribution of repetitive DNA throughout the genome, the relative impacts of retrotransposons and chromosomal duplication in plant genome evolution, and the hypothesized correlation of duplication events with transposition. Population genetic processes also affect the evolution of genomes. We discuss how centromeric genes should, in theory, contain less genetic diversity than noncentromeric genes. In addition, studies of diversity in the wild relatives of maize indicate that different genes have different histories and also show that domestication and intensive breeding have had heterogeneous effects on genetic diversity across genes.

A novel application of gene arrays: Escherichia coli array provides insight into the biology of the obligate endosymbiont of tsetse flies

Symbiotic associations with microorganisms are pivotal in many insects. Yet, the functional roles of obligate symbionts have been difficult to study because it has not been possible to cultivate these organisms in vitro. The medically important tsetse fly (Diptera: Glossinidae) relies on its obligate endosymbiont, Wigglesworthia glossinidia, a member of the Enterobacteriaceae, closely related to Escherichia coli, for fertility and possibly nutrition. We show here that the intracellular Wigglesworthia has a reduced genome size smaller than 770 kb. In an attempt to understand the composition of its genome, we used the gene arrays developed for E. coli. We were able to identify 650 orthologous genes in Wigglesworthia corresponding to ≈85% of its genome. The arrays were also applied for expression analysis using Wigglesworthia cDNA and 61 gene products were detected, presumably coding for some of its most abundant products. Overall, genes involved in cell processes, DNA replication, transcription, and translation were found largely retained in the small genome of Wigglesworthia. In addition, genes coding for transport proteins, chaperones, biosynthesis of cofactors, and some amino acids were found to comprise a significant portion, suggesting an important role for these proteins in its symbiotic life. Based on its expression profile, we predict that Wigglesworthia may be a facultative anaerobic organism that utilizes ammonia as its major source of nitrogen. We present an application of E. coli gene arrays to obtain broad genome information for a closely related organism in the absence of complete genome sequence data.

Two highly homologous ribonuclease genes expressed in mouse eosinophils identify a larger subgroup of the mammalian ribonuclease superfamily.

Two putative ribonucleases have been isolated from the secondary granules of mouse eosinophils. Degenerate oligonucleotide primers inferred from peptide sequence data were used in reverse transcriptase-PCR reactions of bone marrow-derived cDNA. The resulting PCR product was used to screen a C57BL/6J bone marrow cDNA library, and comparisons of representative clones showed that these genes and encoded proteins are highly homologous (96% identity at the nucleotide level; 92/94% identical/similar at the amino acid level). The mouse proteins are only weakly homologous (approximately 50% amino acid identity) with the human eosinophil-associated ribonucleases (i.e., eosinophil-derived neurotoxin and eosinophil cationic protein) and show no sequence bias toward either human protein. Phylogenetic analyses established that the human and mouse loci shared an ancestral gene, but that independent duplication events have occurred since the divergence of primates and rodents. The duplication event generating the mouse genes was estimated to have occurred < 5 x 10(6) years ago (versus 30 to 40 x 10(6) years ago in primates). The identification of independent duplication events in two extant mammalian orders suggests a selective advantage to having multiple eosinophil granule ribonucleases. Southern blot analyses in the mouse demonstrated the existence of three additional highly homologous genes (i.e., five genes total) as well as several more divergent family members. The potential significance of this observation is the implication of a larger gene subfamily in primates (i.e., humans).

Comparison of circumsporozoite proteins from avian and mammalian malarias: biological and phylogenetic implications.

The circumsporozoite (CS) protein of malaria parasites (Plasmodium) covers the surface of sporozoites that invade hepatocytes in mammalian hosts and macrophages in avian hosts. CS genes have been characterized from many Plasmodium that infect mammals; two domains of the corresponding proteins, identified initially by their conservation (region I and region II), have been implicated in binding to hepatocytes. The CS gene from the avian parasite Plasmodium gallinaceum was characterized to compare these functional domains to those of mammalian Plasmodium and for the study of Plasmodium evolution. The P. gallinaceum protein has the characteristics of CS proteins, including a secretory signal sequence, central repeat region, regions of charged amino acids, and an anchor sequence. Comparison with CS signal sequences reveals four distinct groupings, with P. gallinaceum most closely related to the human malaria Plasmodium falciparum. The 5-amino acid sequence designated region I, which is identical in all mammalian CS and implicated in hepatocyte invasion, is different in the avian protein. The P. gallinaceum repeat region consists of 9-amino acid repeats with the consensus sequence QP(A/V)GGNGG(A/V). The conserved motif designated region II-plus, which is associated with targeting the invasion of liver cells, is also conserved in the avian protein. Phylogenetic analysis of the aligned Plasmodium CS sequences yields a tree with a topology similar to the one obtained using sequence data from the small subunit rRNA gene. The phylogeny using the CS gene supports the proposal that the human malaria P. falciparum is significantly more related to avian parasites than to other parasites infecting mammals, although the biology of sporozoite invasion is different between the avian and mammalian species.

The miniaturized nuclear genome of eukaryotic endosymbiont contains genes that overlap, genes that are cotranscribed, and the smallest known spliceosomal introns.

Chlorarachniophyte algae contain a complex, multi-membraned chloroplast derived from the endosymbiosis of a eukaryotic alga. The vestigial nucleus of the endosymbiont, called the nucleomorph, contains only three small linear chromosomes with a haploid genome size of 380 kb and is the smallest known eukaryotic genome. Nucleotide sequence data from a subtelomeric fragment of chromosome III were analyzed as a preliminary investigation of the coding capacity of this vestigial genome. Several housekeeping genes including U6 small nuclear RNA (snRNA), ribosomal proteins S4 and S13, a core protein of the spliceosome [small nuclear ribonucleoprotein (snRNP) E], and a cip-like protease (clpP) were identified. Expression of these genes was confirmed by combinations of Northern blot analysis, in situ hybridization, immunocytochemistry, and cDNA analysis. The protein-encoding genes are typically eukaryotic in overall structure and their messenger RNAs are polyadenylylated. A novel feature is the abundance of 18-, 19-, or 20-nucleotide introns; the smallest spliceosomal introns known. Two of the genes, U6 and S13, overlap while another two genes, snRNP E and clpP, are cotranscribed in a single mRNA. The overall gene organization is extraordinarily compact, making the nucleomorph a unique model for eukaryotic genomics.

Bootstrap confidence levels for phylogenetic trees.

Evolutionary trees are often estimated from DNA or RNA sequence data. How much confidence should we have in the estimated trees? In 1985, Felsenstein [Felsenstein, J. (1985) Evolution 39, 783-791] suggested the use of the bootstrap to answer this question. Felsenstein's method, which in concept is a straightforward application of the bootstrap, is widely used, but has been criticized as biased in the genetics literature. This paper concerns the use of the bootstrap in the tree problem. We show that Felsenstein's method is not biased, but that it can be corrected to better agree with standard ideas of confidence levels and hypothesis testing. These corrections can be made by using the more elaborate bootstrap method presented here, at the expense of considerably more computation.

Cloning and expression of human deoxyguanosine kinase cDNA.

A human cDNA sequence homologous to human deoxycytidine kinase (dCK; EC 2.7.1.74) was identified in the GenBank sequence data base. The longest open reading frame encoded a protein that was 48% identical to dCK at the amino acid level. The cDNA was expressed in Escherichia coli and shown to encode a protein with the same substrate specificity as described for the mitochondrial deoxyguanosine kinase (dGK; EC 2.7.1.113). The N terminus of the deduced amino acid sequence had properties characteristic for a mitochondrial translocation signal, and cleavage at a putative mitochondrial peptidase cleavage site would give a mature protein size of 28 kDa. Northern blot analysis determined the length of dGK mRNA to 1.3 kbp with no cross-hybridization to the 2.8-kbp dCK mRNA. dGK mRNA was detected in all tissues investigated with the highest expression levels in muscle, brain, liver, and lymphoid tissues. Alignment of the dGK and herpes simplex virus type 1 thymidine kinase amino acid sequences showed that five regions, including the substrate-binding pocket and the ATP-binding glycine loop, were also conserved in dGK. To our knowledge, this is the first report of a cloned mitochondrial nucleoside kinase and the first demonstration of a general sequence homology between two mammalian deoxyribonucleoside kinases. Our findings suggest that dCK and dGK are evolutionarily related, as well as related to the family of herpes virus thymidine kinases.

Discontinuous epitopes of hepatitis B surface antigen derived from a filamentous phage peptide library.

The structure of the small hepatitis B virus surface antigen (HBsAg) was investigated by epitope mapping of four anti-HBsAg monoclonal antibodies (mAbs). Amino acid sequences of epitopes were derived from affinity-enrichment experiments (biopanning) using a filamentous phage peptide library. The library consists of 10(9) different clones bearing a 30-residue peptide fused to gene III. Sequence homologies between peptides obtained from panning the library against the antibodies and the native HBsAg sequence allowed for precise description of the binding regions. Three of four mAbs were found to bind to distinct discontinuous epitopes between amino acid residues 101 and 207 of HBsAg. The fourth mAb was demonstrated to bind to residues 121-124. The sequence data are supported by ELISA assays demonstrating the binding of the HBsAg-specific peptides on filamentous phage to mAbs. The sequence data were used to map the surface of HBsAg and to derive a topological model for the alpha-carbon trace of the 101-207 region of HBsAg. The approach should be useful for other proteins for which the crystal structure is not available but a representative set of mAbs can be obtained.

Trinucleotide repeats and long homopeptides in genes and proteins associated with nervous system disease and development.

Several human neurological disorders are associated with proteins containing abnormally long runs of glutamine residues. Strikingly, most of these proteins contain two or more additional long runs of amino acids other than glutamine. We screened the current human, mouse, Drosophila, yeast, and Escherichia coli protein sequence data bases and identified all proteins containing multiple long homopeptides. This search found multiple long homopeptides in about 12% of Drosophila proteins but in only about 1.7% of human, mouse, and yeast proteins and none among E. coli proteins. Most of these sequences show other unusual sequence features, including multiple charge clusters and excessive counts of homopeptides of length > or = two amino acid residues. Intriguingly, a large majority of the identified Drosophila proteins are essential developmental proteins and, in particular, most play a role in central nervous system development. Almost half of the human and mouse proteins identified are homeotic homologs. The role of long homopeptides in fine-tuning protein conformation for multiple functional activities is discussed. The relative contributions of strand slippage and of dynamic mutation are also addressed. Several new experiments are proposed.

The 145-kDa protein induced to associate with Shc by multiple cytokines is an inositol tetraphosphate and phosphatidylinositol 3,4,5-triphosphate 5-phosphatase.

A 145-kDa tyrosine-phosphorylated protein that becomes associated with Shc in response to multiple cytokines has been purified from the murine hemopoietic cell line B6SUtA1. Amino acid sequence data were used to clone the cDNA encoding this protein from a B6SUtA1 library. The predicted amino acid sequence encodes a unique protein containing an N-terminal src homology 2 domain, two consensus sequences that are targets for phosphotyrosine binding domains, a proline-rich region, and two motifs highly conserved among inositol polyphosphate 5-phosphatases. Cell lysates immunoprecipitated with antiserum to this protein exhibited both phosphatidylinositol 3,4,5-trisphosphate and inositol 1,3,4,5-tetrakisphosphate polyphosphate 5-phosphatase activity. This novel signal transduction intermediate may serve to modulate both Ras and inositol signaling pathways. Based on its properties, we suggest the 145-kDa protein be called SHIP for SH2-containing inositol phosphatase.

A Fas-associated protein factor, FAF1, potentiates Fas-mediated apoptosis.

Fas, a member of the tumor necrosis factor receptor family, can induce apoptosis when activated by Fas ligand binding or anti-Fas antibody crosslinking. Genetic studies have shown that a defect in Fas-mediated apoptosis resulted in abnormal development and function of the immune system in mice. A point mutation in the cytoplasmic domain of Fas (a single base change from T to A at base 786), replacing isoleucine with asparagine, abolishes the signal transducing property of Fas. Mice homozygous for this mutant allele (lprcg/lprcg mice) develop lymphadenopathy and a lupus-like autoimmune disease. Little is known about the mechanism of signal transduction in Fas-mediated apoptosis. In this study, we used the two-hybrid screen in yeast to isolate a Fas-associated protein factor, FAF1, which specifically interacts with the cytoplasmic domain of wild-type Fas but not the lprcg-mutated Fas protein. This interaction occurs not only in yeast but also in mammalian cells. When transiently expressed in L cells, FAF1 potentiated Fas-induced apoptosis. A search of available DNA and protein sequence data banks did not reveal significant homology between FAF1 and known proteins. Therefore, FAF1 is an unusual protein that binds to the wild type but not the inactive point mutant of Fas. FAF1 potentiates Fas-induced cell killing and is a candidate signal transducing molecule in the regulation of apoptosis.