The virion glycoproteins of Ebola viruses are encoded in two reading frames and are expressed through transcriptional editing.

In late 1994 and early 1995, Ebola (EBO) virus dramatically reemerged in Africa, causing human disease in the Ivory Coast and Zaire. Analysis of the entire glycoprotein genes of these viruses and those of other EBO virus subtypes has shown that the virion glycoprotein (130 kDa) is encoded in two reading frames, which are linked by transcriptional editing. This editing results in the addition of an extra nontemplated adenosine within a run of seven adenosines near the middle of the coding region. The primary gene product is a smaller (50-70 kDa), nonstructural, secreted glycoprotein, which is produced in large amounts and has an unknown function. Phylogenetic analysis indicates that EBO virus subtypes are genetically diverse and that the recent Ivory Coast isolate represents a new (fourth) subtype of EBO virus. In contrast, the EBO virus isolate from the 1995 outbreak in Kikwit, Zaire, is virtually identical to the virus that caused a similar epidemic in Yambuku, Zaire, almost 20 years earlier. This genetic stability may indicate that EBO viruses have coevolved with their natural reservoirs and do not change appreciably in the wild.

The non-coding RNAs as riboregulators

The non-coding RNAs database (http://biobases.ibch.poznan.pl/ncRNA/) contains currently available data on RNAs, which do not have long open reading frames and act as riboregulators. Non-coding RNAs are involved in the specific recognition of cellular nucleic acid targets through complementary base pairing to control cell growth and differentiation. Some of them are connected with several well known developmental and neuro­behavioral disorders. We have divided them into four groups. This paper is a short introduction to the database and presents its latest, updated edition.

DNA repeats identify novel virulence genes in Haemophilus influenzae.

The whole genome sequence (1.83 Mbp) of Haemophilus influenzae strain Rd was searched to identify tandem oligonucleotide repeat sequences. Loss or gain of one or more nucleotide repeats through a recombination-independent slippage mechanism is known to mediate phase variation of surface molecules of pathogenic bacteria, including H. influenzae. This facilitates evasion of host defenses and adaptation to the varying microenvironments of the host. We reasoned that iterative nucleotides could identify novel genes relevant to microbe-host interactions. Our search of the Rd genome sequence identified 9 novel loci with multiple (range 6-36, mean 22) tandem tetranucleotide repeats. All were found to be located within putative open reading frames and included homologues of hemoglobin-binding proteins of Neisseria, a glycosyltransferase (IgtC gene product) of Neisseria, and an adhesin of Yersinia. These tetranucleotide repeat sequences were also shown to be present in two other epidemiologically different H. influenzae type b strains, although the number and distribution of repeats was different. Further characterization of the IgtC gene showed that it was involved in phenotypic switching of a lipopolysaccharide epitope and that this variable expression was associated with changes in the number of tetranucleotide repeats. Mutation of IgtC resulted in attenuated virulence of H. influenzae in an infant rat model of invasive infection. These data indicate the rapidity, economy, and completeness with which whole genome sequences can be used to investigate the biology of pathogenic bacteria.

Gene disruptions using P transposable elements: an integral component of the Drosophila genome project.

Biologists require genetic as well as molecular tools to decipher genomic information and ultimately to understand gene function. The Berkeley Drosophila Genome Project is addressing these needs with a massive gene disruption project that uses individual, genetically engineered P transposable elements to target open reading frames throughout the Drosophila genome. DNA flanking the insertions is sequenced, thereby placing an extensive series of genetic markers on the physical genomic map and associating insertions with specific open reading frames and genes. Insertions from the collection now lie within or near most Drosophila genes, greatly reducing the time required to identify new mutations and analyze gene functions. Information revealed from these studies about P element site specificity is being used to target the remaining open reading frames.

The nucleotide sequence of chromosome I from Saccharomyces cerevisiae.

Chromosome I from the yeast Saccharomyces cerevisiae contains a DNA molecule of approximately 231 kbp and is the smallest naturally occurring functional eukaryotic nuclear chromosome so far characterized. The nucleotide sequence of this chromosome has been determined as part of an international collaboration to sequence the entire yeast genome. The chromosome contains 89 open reading frames and 4 tRNA genes. The central 165 kbp of the chromosome resembles other large sequenced regions of the yeast genome in both its high density and distribution of genes. In contrast, the remaining sequences flanking this DNA that comprise the two ends of the chromosome and make up more than 25% of the DNA molecule have a much lower gene density, are largely not transcribed, contain no genes essential for vegetative growth, and contain several apparent pseudogenes and a 15-kbp redundant sequence. These terminally repetitive regions consist of a telomeric repeat called W', flanked by DNA closely related to the yeast FLO1 gene. The low gene density, presence of pseudogenes, and lack of expression are consistent with the idea that these terminal regions represent the yeast equivalent of heterochromatin. The occurrence of such a high proportion of DNA with so little information suggests that its presence gives this chromosome the critical length required for proper function.

Tight frames of k-plane ridgelets and the problem of representing objects that are smooth away from d-dimensional singularities in Rn

For each pair (n, k) with 1 ≤ k < n, we construct a tight frame (ρλ : λ ∈ Λ) for L2 (Rn), which we call a frame of k-plane ridgelets. The intent is to efficiently represent functions that are smooth away from singularities along k-planes in Rn. We also develop tools to help decide whether k-plane ridgelets provide the desired efficient representation. We first construct a wavelet-like tight frame on the X-ray bundle χn,k—the fiber bundle having the Grassman manifold Gn,k of k-planes in Rn for base space, and for fibers the orthocomplements of those planes. This wavelet-like tight frame is the pushout to χn,k, via the smooth local coordinates of Gn,k, of an orthonormal basis of tensor Meyer wavelets on Euclidean space Rk(n−k) × Rn−k. We then use the X-ray isometry [Solmon, D. C. (1976) J. Math. Anal. Appl. 56, 61–83] to map this tight frame isometrically to a tight frame for L2(Rn)—the k-plane ridgelets. This construction makes analysis of a function f ∈ L2(Rn) by k-plane ridgelets identical to the analysis of the k-plane X-ray transform of f by an appropriate wavelet-like system for χn,k. As wavelets are typically effective at representing point singularities, it may be expected that these new systems will be effective at representing objects whose k-plane X-ray transform has a point singularity. Objects with discontinuities across hyperplanes are of this form, for k = n − 1.

Three members of a novel small gene-family from Arabidopsis thaliana able to complement functionally an Escherichia coli mutant defective in PAPS reductase activity encode proteins with a thioredoxin-like domain and “APS reductase” activity

Three different cDNAs, Prh-19, Prh-26, and Prh-43 [3′-phosphoadenosine-5′-phosphosulfate (PAPS) reductase homolog], have been isolated by complementation of an Escherichia coli cysH mutant, defective in PAPS reductase activity, to prototrophy with an Arabidopsis thaliana cDNA library in the expression vector λYES. Sequence analysis of the cDNAs revealed continuous open reading frames encoding polypeptides of 465, 458, and 453 amino acids, with calculated molecular masses of 51.3, 50.5, and 50.4 kDa, respectively, that have strong homology with fungal, yeast, and bacterial PAPS reductases. However, unlike microbial PAPS reductases, each PRH protein has an N-terminal extension, characteristic of a plastid transit peptide, and a C-terminal extension that has amino acid and deduced three-dimensional homology to thioredoxin proteins. Adenosine 5′-phosphosulfate (APS) was shown to be a much more efficient substrate than PAPS when the activity of the PRH proteins was tested by their ability to convert 35S-labeled substrate to acid-volatile 35S-sulfite. We speculate that the thioredoxin-like domain is involved in catalytic function, and that the PRH proteins may function as novel “APS reductase” enzymes. Southern hybridization analysis showed the presence of a small multigene family in the Arabidopsis genome. RNA blot hybridization with gene-specific probes revealed for each gene the presence of a transcript of ≈1.85 kb in leaves, stems, and roots that increased on sulfate starvation. To our knowledge, this is the first report of the cloning and characterization of plant genes that encode proteins with APS reductase activity and supports the suggestion that APS can be utilized directly, without activation to PAPS, as an intermediary substrate in reductive sulfate assimilation.

Linking genome and proteome by mass spectrometry: Large-scale identification of yeast proteins from two dimensional gels

The function of many of the uncharacterized open reading frames discovered by genomic sequencing can be determined at the level of expressed gene products, the proteome. However, identifying the cognate gene from minute amounts of protein has been one of the major problems in molecular biology. Using yeast as an example, we demonstrate here that mass spectrometric protein identification is a general solution to this problem given a completely sequenced genome. As a first screen, our strategy uses automated laser desorption ionization mass spectrometry of the peptide mixtures produced by in-gel tryptic digestion of a protein. Up to 90% of proteins are identified by searching sequence data bases by lists of peptide masses obtained with high accuracy. The remaining proteins are identified by partially sequencing several peptides of the unseparated mixture by nanoelectrospray tandem mass spectrometry followed by data base searching with multiple peptide sequence tags. In blind trials, the method led to unambiguous identification in all cases. In the largest individual protein identification project to date, a total of 150 gel spots—many of them at subpicomole amounts—were successfully analyzed, greatly enlarging a yeast two-dimensional gel data base. More than 32 proteins were novel and matched to previously uncharacterized open reading frames in the yeast genome. This study establishes that mass spectrometry provides the required throughput, the certainty of identification, and the general applicability to serve as the method of choice to connect genome and proteome.

HDA1 and RPD3 are members of distinct yeast histone deacetylase complexes that regulate silencing and transcription

Increased histone acetylation has been correlated with increased transcription, and regions of heterochromatin are generally hypoacetylated. In investigating the cause-and-effect relationship between histone acetylation and gene activity, we have characterized two yeast histone deacetylase complexes. Histone deacetylase-A (HDA) is an ≈350-kDa complex that is highly sensitive to the deacetylase inhibitor trichostatin A. Histone deacetylase-B (HDB) is an ≈600-kDa complex that is much less sensitive to trichostatin A. The HDA1 protein (a subunit of the HDA activity) shares sequence similarity to RPD3, a factor required for optimal transcription of certain yeast genes. RPD3 is associated with the HDB activity. HDA1 also shares similarity to three new open reading frames in yeast, designated HOS1, HOS2, and HOS3. We find that both hda1 and rpd3 deletions increase acetylation levels in vivo at all sites examined in both core histones H3 and H4, with rpd3 deletions having a greater impact on histone H4 lysine positions 5 and 12. Surprisingly, both hda1 and rpd3 deletions increase repression at telomeric loci, which resemble heterochromatin with rpd3 having a greater effect. In addition, rpd3 deletions retard full induction of the PHO5 promoter fused to the reporter lacZ. These data demonstrate that histone acetylation state has a role in regulating both heterochromatic silencing and regulated gene expression.

The interphotoreceptor retinoid binding protein gene in therian mammals: Implications for higher level relationships and evidence for loss of function in the marsupial mole

The subclass Theria of Mammalia includes marsupials (infraclass Metatheria) and placentals (infraclass Eutheria). Within each group, interordinal relationships remain unclear. One limitation of many studies is incomplete ordinal representation. Here, we analyze DNA sequences for part of exon 1 of the interphotoreceptor retinoid binding protein gene, including 10 that are newly reported, for representatives of all therian orders. Among placentals, the most robust clades are Cetartiodactyla, Paenungulata, and an expanded African clade that includes paenungulates, tubulidentates, and macroscelideans. Anagalida, Archonta, Altungulata, Hyracoidea + Perissodactyla, Ungulata, and the “flying primate” hypothesis are rejected by statistical tests. Among marsupials, the most robust clade includes all orders except Didelphimorphia. The phylogenetic placement of the monito del monte and the marsupial mole remains unclear. However, the marsupial mole sequence contains three frameshift indels and numerous stop codons in all three reading frames. Given that the interphotoreceptor retinoid binding protein gene is a single-copy gene that functions in the visual cycle and that the marsupial mole is blind with degenerate eyes, this finding suggests that phenotypic degeneration of the eyes is accompanied by parallel changes at the molecular level as a result of relaxed selective constraints.

VIDA: a virus database system for the organization of animal virus genome open reading frames

VIDA is a new virus database that organizes open reading frames (ORFs) from partial and complete genomic sequences from animal viruses. Currently VIDA includes all sequences from GenBank for Herpesviridae, Coronaviridae and Arteriviridae. The ORFs are organized into homologous protein families, which are identified on the basis of sequence similarity relationships. Conserved sequence regions of potential functional importance are identified and can be retrieved as sequence alignments. We use a controlled taxonomical and functional classification for all the proteins and protein families in the database. When available, protein structures that are related to the families have also been included. The database is available for online search and sequence information retrieval at http://www.biochem.ucl.ac.uk/bsm/virus_database/VIDA.html.

Multiple DNA and protein sequence alignment based on segment-to-segment comparison.

In this paper, a new way to think about, and to construct, pairwise as well as multiple alignments of DNA and protein sequences is proposed. Rather than forcing alignments to either align single residues or to introduce gaps by defining an alignment as a path running right from the source up to the sink in the associated dot-matrix diagram, we propose to consider alignments as consistent equivalence relations defined on the set of all positions occurring in all sequences under consideration. We also propose constructing alignments from whole segments exhibiting highly significant overall similarity rather than by aligning individual residues. Consequently, we present an alignment algorithm that (i) is based on segment-to-segment comparison instead of the commonly used residue-to-residue comparison and which (ii) avoids the well-known difficulties concerning the choice of appropriate gap penalties: gaps are not treated explicity, but remain as those parts of the sequences that do not belong to any of the aligned segments. Finally, we discuss the application of our algorithm to two test examples and compare it with commonly used alignment methods. As a first example, we aligned a set of 11 DNA sequences coding for functional helix-loop-helix proteins. Though the sequences show only low overall similarity, our program correctly aligned all of the 11 functional sites, which was a unique result among the methods tested. As a by-product, the reading frames of the sequences were identified. Next, we aligned a set of ribonuclease H proteins and compared our results with alignments produced by other programs as reported by McClure et al. [McClure, M. A., Vasi, T. K. & Fitch, W. M. (1994) Mol. Biol. Evol. 11, 571-592]. Our program was one of the best scoring programs. However, in contrast to other methods, our protein alignments are independent of user-defined parameters.

Human cytomegalovirus US3 impairs transport and maturation of major histocompatibility complex class I heavy chains.

The human cytomegalovirus (HCMV) early glycoprotein products of the US11 and US2 open reading frames cause increased turnover of major histocompatibility complex (MHC) class I heavy chains. Since US2 is homologous to another HCMV gene (US3), we hypothesized that the US3 gene product also may affect MHC class I expression. In cells constitutively expressing the HCMV US3 gene, MHC class I heavy chains formed a stable complex with beta 2-microglobulin. However, maturation of the N-linked glycan of MHC class I heavy chains was impaired in US3+ cells. The glycoprotein product of US3 (gpUS3) occurs mostly in a high-mannose form and coimmunoprecipitates with beta 2-microglobulin associated class I heavy chains. Mature class I molecules were detected at steady state on the surface of US3+ cells, as in control cells. Substantial perinuclear accumulation of heavy chains was observed in US3+ cells. The data suggest that gpUS3 impairs egress of MHC class I heavy chains from the endoplasmic reticulum.

A chromosomal locus required for copper resistance, competitive fitness, and cytochrome c biogenesis in Pseudomonas fluorescens.

A chromosomal locus required for copper resistance and competitive fitness was cloned from a strain of Pseudomonas fluorescens isolated from copper-contaminated agricultural soil. Sequence analysis of this locus revealed six open reading frames with homology to genes involved in cytochrome c biogenesis in other bacteria, helC, cycJ, cycK, tipB, cycL, and cycH, with the closest similarity being to the aeg-46.5(yej) region of the Escherichia coli chromosome. The proposed functions of these genes in other bacteria include the binding, transport, and coupling of heme to apocytochrome c in the periplasm of these Gram-negative bacteria. Putative heme-binding motifs were present in the predicted products of cycK and cycL, and TipB contained a putative disulfide oxidoreductase active site proposed to maintain the heme-binding site of the apocytochrome in a reduced state for ligation of heme. Tn3-gus mutagenesis showed that expression of the genes was constitutive but enhanced by copper, and confirmed that the genes function both in copper resistance and production of active cytochrome c. However, two mutants in cycH were copper-sensitive and oxidase-positive, suggesting that the functions of these genes, rather than cytochrome c oxidase itself, were required for resistance to copper.

Genetic rearrangements in the rfb regions of Vibrio cholerae O1 and O139.

The recent emergence of a pathogenic new non-O1 serotype (O139) of Vibrio cholerae has led to numerous studies in an attempt to identify the origins of this new strain. Our studies indicate that O139 strains have clear differences in the surface polysaccharides when compared with O1 strains: the lipopolysaccharide can be described as semi-rough. Southern hybridization with the O1 rfb region demonstrates that O139 strains no longer contain any of the rfb genes required for the synthesis of the O1 O-antigen or its modification and also lack at least 6 kb of additional contiguous DNA. However, O139 strains have retained rfaD and have a single open reading frame closely related to three small open reading frames of the O1 rfb region. This region is closely related to the H-repeat of Escherichia coli and to the transposases of a number of insertion sequence elements and has all the features of an insertion sequence element that has been designated VcIS1. Transposon insertion mutants defective in O139 O-antigen (and capsule) biosynthesis map to the same fragment as VcIS1. Preliminary sequence data of complementing clones indicate that this DNA encodes a galactosyl-transferase and other enzymes for the utilization of galactose in polysaccharide biosynthesis. We propose a mechanism by which both the Ogawa serotype of O1 strains and the O139 serotype strains may have evolved.