Are There 1031 Virus Particles on Earth, or More, or Fewer?

The number of virus particles on Earth is frequently reported in the scientific literature and in general-interest publications as being on the order of 1031, with some confusion about whether this is a high or low estimate. This number is often given without a source, although it should be attributed to a paper by Hendrix et al. published in 1999 (R. W. Hendrix, M. C. Smith, R. N. Burns, M. E. Ford, and G. F. Hatfull, Proc Natl Acad Sci U S A 96:2192-2197, 1999, https://doi.org/10.1073/pnas.96.5.2192). As with any oft-repeated statistic, it is informative to know how it has been derived and whether it should be revised in the light of new evidence. I review the history of the 1031 estimate and use more recent assessments of the number of bacterial and viral particles in various habitats to conclude that the best estimate of the number of virus particles on Earth ("the Hendrix product") remains close to 1031 and is unlikely to be either much less or much more than that.

Functional specialization and evolution of leader proteinases in the family Closteroviridae.

Members of the Closteroviridae and Potyviridae families of the plant positive-strand RNA viruses encode one or two papain-like leader proteinases. In addition to a C-terminal proteolytic domain, each of these proteinases possesses a nonproteolytic N-terminal domain. We compared functions of the several leader proteinases using a gene swapping approach. The leader proteinase (L-Pro) of Beet yellows virus (BYV; a closterovirus) was replaced with L1 or L2 proteinases of Citrus tristeza virus (CTV; another closterovirus), P-Pro proteinase of Lettuce infectious yellows virus (LIYV; a crinivirus), and HC-Pro proteinase of Tobacco etch virus (a potyvirus). Each foreign proteinase efficiently processed the chimeric BYV polyprotein in vitro. However, only L1 and P-Pro, not L2 and HC-Pro, were able to rescue the amplification of the chimeric BYV variants. The combined expression of L1 and L2 resulted in an increased RNA accumulation compared to that of the parental BYV. Remarkably, this L1-L2 chimera exhibited reduced invasiveness and inability to move from cell to cell. Similar analyses of the BYV hybrids, in which only the papain-like domain of L-Pro was replaced with those derived from L1, L2, P-Pro, and HC-Pro, also revealed functional specialization of these domains. In subcellular-localization experiments, distinct patterns were observed for the leader proteinases of BYV, CTV, and LIYV. Taken together, these results demonstrated that, in addition to a common proteolytic activity, the leader proteinases of closteroviruses possess specialized functions in virus RNA amplification, virus invasion, and cell-to-cell movement. The phylogenetic analysis suggested that functionally distinct L1 and L2 of CTV originated by a gene duplication event.

Quod erat demonstrandum? The mystery of experimental validation of apparently erroneous computational analyses of protein sequences.

BACKGROUND: Computational predictions are critical for directing the experimental study of protein functions. Therefore it is paradoxical when an apparently erroneous computational prediction seems to be supported by experiment. RESULTS: We analyzed six cases where application of novel or conventional computational methods for protein sequence and structure analysis led to non-trivial predictions that were subsequently supported by direct experiments. We show that, on all six occasions, the original prediction was unjustified, and in at least three cases, an alternative, well-supported computational prediction, incompatible with the original one, could be derived. The most unusual cases involved the identification of an archaeal cysteinyl-tRNA synthetase, a dihydropteroate synthase and a thymidylate synthase, for which experimental verifications of apparently erroneous computational predictions were reported. Using sequence-profile analysis, multiple alignment and secondary-structure prediction, we have identified the unique archaeal 'cysteinyl-tRNA synthetase' as a homolog of extracellular polygalactosaminidases, and the 'dihydropteroate synthase' as a member of the beta-lactamase-like superfamily of metal-dependent hydrolases. CONCLUSIONS: In each of the analyzed cases, the original computational predictions could be refuted and, in some instances, alternative strongly supported predictions were obtained. The nature of the experimental evidence that appears to support these predictions remains an open question. Some of these experiments might signify discovery of extremely unusual forms of the respective enzymes, whereas the results of others could be due to artifacts.

Conserved phosphoprotein interaction motif is functionally interchangeable between ataxin-7 and arrestins.

Olivopontocerebellar atrophy with retinal degeneration is a hereditary neurodegenerative disorder that belongs to the subtype II of the autosomal dominant cerebellar ataxias and is characterized by early-onset cerebellar and macular degeneration preceded by diagnostically useful tritan colorblindness. The gene mutated in the disease (SCA7) has been mapped to chromosome 3p12-13.5, and positional cloning identified the cause of the disease as CAG repeat expansion in this gene. The SCA7 gene product, ataxin-7, is an 897 amino acid protein with an expandable polyglutamine tract close to its N-terminus. No clues to ataxin-7 function have been obtained from sequence database searches. Here we report that ataxin-7 has a motif of ca. 50 amino acids, related to the phosphate-binding site of arrestins. To test the relevance of this sequence similarity, we introduced the putative ataxin-7 phosphate-binding site into visual arrestin and beta-arrestin. Both chimeric arrestins retain receptor-binding affinity and show characteristic high selectivity for phosphorylated activated forms of rhodopsin and beta-adrenergic receptor, respectively. Although the insertion of a Gly residue (absent in arrestins but present in the putative phosphate-binding site of ataxin-7) disrupts the function of visual arrestin-ataxin-7 chimera, it enhances the function of beta-arrestin-ataxin-7 chimera. Taken together, our data suggest that the arrestin-like site in the ataxin-7 sequence is a functional phosphate-binding site. The presence of the phosphate-binding site in ataxin-7 suggests that this protein may be involved in phosphorylation-dependent binding to its protein partner(s) in the cell.

Large-scale taxonomic profiling of eukaryotic model organisms: a comparison of orthologous proteins encoded by the human, fly, nematode, and yeast genomes.

Comparisons of DNA and protein sequences between humans and model organisms, including the yeast Saccharomyces cerevisiae, the nematode Caenorhabditis elegans, and the fruit fly Drosophila melanogaster, are a significant source of information about the function of human genes and proteins in both normal and disease states. Important questions regarding cross-species sequence comparison remain unanswered, including (1) the fraction of the metabolic, signaling, and regulatory pathways that is shared by humans and the various model organisms; and (2) the validity of functional inferences based on sequence homology. We addressed these questions by analyzing the available fractions of human, fly, nematode, and yeast genomes for orthologous protein-coding genes, applying strict criteria to distinguish between candidate orthologous and paralogous proteins. Forty-two quartets of proteins could be identified as candidate orthologs. Twenty-four Drosophila protein sequences were more similar to their human orthologs than the corresponding nematode proteins. Analysis of sequence substitutions and evolutionary distances in this data set revealed that most C. elegans genes are evolving more rapidly than Drosophila genes, suggesting that unequal evolutionary rates may contribute to the differences in similarity to human protein sequences. The available fraction of Drosophila proteins appears to lack representatives of many protein families and domains, reflecting the relative paucity of genomic data from this species.

Comparison of archaeal and bacterial genomes: computer analysis of protein sequences predicts novel functions and suggests a chimeric origin for the archaea.

Protein sequences encoded in three complete bacterial genomes, those of Haemophilus influenzae, Mycoplasma genitalium and Synechocystis sp., and the first available archaeal genome sequence, that of Methanococcus jannaschii, were analysed using the BLAST2 algorithm and methods for amino acid motif detection. Between 75% and 90% of the predicted proteins encoded in each of the bacterial genomes and 73% of the M. jannaschii proteins showed significant sequence similarity to proteins from other species. The fraction of bacterial and archaeal proteins containing regions conserved over long phylogenetic distances is nearly the same and close to 70%. Functions of 70-85% of the bacterial proteins and about 70% of the archaeal proteins were predicted with varying precision. This contrasts with the previous report that more than half of the archaeal proteins have no homologues and shows that, with more sensitive methods and detailed analysis of conserved motifs, archaeal genomes become as amenable to meaningful interpretation by computer as bacterial genomes. The analysis of conserved motifs resulted in the prediction of a number of previously undetected functions of bacterial and archaeal proteins and in the identification of novel protein families. In spite of the generally high conservation of protein sequences, orthologues of 25% or less of the M. jannaschii genes were detected in each individual completely sequenced genome, supporting the uniqueness of archaea as a distinct domain of life. About 53% of the M. jannaschii proteins belong to families of paralogues, a fraction similar to that in bacteria with larger genomes, such as Synechocystis sp. and Escherichia coli, but higher than that in H. influenzae, which has approximately the same number of genes as M. jannaschii. Certain groups of proteins, e.g. molecular chaperones and DNA repair enzymes, thought to be ubiquitous and represented in the minimal gene set derived by bacterial genome comparison, are missing in M. jannaschii, indicating massive non-orthologous displacement of genes responsible for essential functions. An unexpectedly large fraction of the M. jannaschii gene products, 44%, shows significantly higher similarity to bacterial than to eukaryotic proteins, compared with 13% that have eukaryotic proteins as their closest homologues (the rest of the proteins show approximately the same level of similarity to bacterial and eukaryotic homologues or have no homologues). Proteins involved in translation, transcription, replication and protein secretion are most closely related to eukaryotic proteins, whereas metabolic enzymes, metabolite uptake systems, enzymes for cell wall biosynthesis and many uncharacterized proteins appear to be 'bacterial'. A similar prevalence of proteins of apparent bacterial origin was observed among the currently available sequences from the distantly related archaeal genus, Sulfolobus. It is likely that the evolution of archaea included at least one major merger between ancestral cells from the bacterial lineage and the lineage leading to the eukaryotic nucleocytoplasm.

Positionally cloned human disease genes: patterns of evolutionary conservation and functional motifs.

Positional cloning has already produced the sequences of more than 70 human genes associated with specific diseases. In addition to their medical importance, these genes are of interest as a set of human genes isolated solely on the basis of the phenotypic effect of the respective mutations. We analyzed the protein sequences encoded by the positionally cloned disease genes using an iterative strategy combining several sensitive computer methods. Comparisons to complete sequence databases and to separate databases of nematode, yeast, and bacterial proteins showed that for most of the disease gene products, statistically significant sequence similarities are detectable in each of the model organisms. Only the nematode genome encodes apparent orthologs with conserved domain architecture for the majority of the disease genes. In yeast and bacterial homologs, domain organization is typically not conserved, and sequence similarity is limited to individual domains. Generally, human genes complement mutations only in orthologous yeast genes. Most of the positionally cloned genes encode large proteins with several globular and nonglobular domains, the functions of some or all of which are not known. We detected conserved domains and motifs not described previously in a number of proteins encoded by disease genes and predicted functions for some of them. These predictions include an ATP-binding domain in the product of hereditary nonpolyposis colon cancer gene (a MutL homolog), which is conserved in the HS90 family of chaperone proteins, type II DNA topoisomerases, and histidine kinases, and a nuclease domain homologous to bacterial RNase D and the 3'-5' exonuclease domain of DNA polymerase I in the Werner syndrome gene product.

Discrete regions of the sensor protein virA determine the strain-specific ability of Agrobacterium to agroinfect maize.

The ability of Agrobacterium strains to infect transformation-recalcitrant maize plants has been shown to be determined mainly by the virA locus, implicating vir gene induction as the major factor influencing maize infection. In this report, we further explore the roles of vir induction-associated bacterial factors in maize infection using the technique of agroinfection. The Ti plasmid and virA source are shown to be important in determining the ability of a strain to infect maize, and the monosaccharide binding protein ChvE is absolutely required for maize agroinfection. The linker domain of VirAC58 from an agroinfection-competent strain, C58, is sufficient to convert VirAA6 of a nonagroinfecting strain, A348,to agroinfection competence. The periplasmic domain of VirAC58 is also able to confer a moderate level of agroinfection competence to VirAA6. In addition, the VirAA6 protein from A348 is agroinfection competent when removed from its cognate Ti plasmid background and placed in a pTiC58 background. The presence of a pTiA6-encoded, VirAA6-specific inhibitor is hypothesized and examined.

Complete genome sequences of cellular life forms: glimpses of theoretical evolutionary genomics.

The availability of complete genome sequences of cellular life forms creates the opportunity to explore the functional content of the genomes and evolutionary relationships between them at a new qualitative level. With the advent of these sequences, the construction of a minimal gene set sufficient for sustaining cellular life and reconstruction of the genome of the last common ancestor of bacteria, eukaryotes, and archaea become realistic, albeit challenging, research projects. A version of the minimal gene set for modern-type cellular life derived by comparative analysis of two bacterial genomes, those of Haemophilus influenzae and Mycoplasma genitalium, consists of approximately 250 genes. A comparison of the protein sequences encoded in these genes with those of the proteins encoded in the complete yeast genome suggests that the last common ancestor of all extant life might have had an RNA genome.

Sequence analysis of eukaryotic developmental proteins: ancient and novel domains.

Most of the genes involved in the development of multicellular eukaryotes encode large, multidomain proteins. To decipher the major trends in the evolution of these proteins and make functional predictions for uncharacterized domains, we applied a strategy of sequence database search that includes construction of specialized data sets and iterative subsequence masking. This computational approach allowed us to detect previously unnoticed but potentially important sequence similarities. Developmental gene products are enriched in predicted nonglobular regions as compared to unbiased sets of eukaryotic and bacterial proteins. Developmental genes that act intracellularly, primarily at the level of transcription regulation, typically code for proteins containing highly conserved DNA-binding domains, most of which appear to have evolved before the radiation of bacteria and eukaryotes. We identified bacterial homologues, namely a protein family that includes the Escherichia coli universal stress protein UspA, for the MADS-box transcription regulators previously described only in eukaryotes. We also show that the FUS6 family of eukaryotic proteins contains a putative DNA-binding domain related to bacterial helix-turn-helix transcription regulators. Developmental proteins that act extracellularly are less conserved and often do not have bacterial homologues. Nevertheless, several provocative similarities between different groups of such proteins were detected.

The SPO1 gene product required for meiosis in yeast has a high similarity to phospholipase B enzymes.

The SPO1 gene of Saccharomyces cerevisiae has been cloned and sequenced. The Spo1 protein reveals significant similarity with fungal phospholipase B (PLB) enzymes. Features of the SPO1 gene sequence are presented.

A minimal gene set for cellular life derived by comparison of complete bacterial genomes.

The recently sequenced genome of the parasitic bacterium Mycoplasma genitalium contains only 468 identified protein-coding genes that have been dubbed a minimal gene complement [Fraser, C.M., Gocayne, J.D., White, O., Adams, M.D., Clayton, R.A., et al. (1995) Science 270, 397-403]. Although the M. genitalium gene complement is indeed the smallest among known cellular life forms, there is no evidence that it is the minimal self-sufficient gene set. To derive such a set, we compared the 468 predicted M. genitalium protein sequences with the 1703 protein sequences encoded by the other completely sequenced small bacterial genome, that of Haemophilus influenzae. M. genitalium and H. influenzae belong to two ancient bacterial lineages, i.e., Gram-positive and Gram-negative bacteria, respectively. Therefore, the genes that are conserved in these two bacteria are almost certainly essential for cellular function. It is this category of genes that is most likely to approximate the minimal gene set. We found that 240 M. genitalium genes have orthologs among the genes of H. influenzae. This collection of genes falls short of comprising the minimal set as some enzymes responsible for intermediate steps in essential pathways are missing. The apparent reason for this is the phenomenon that we call nonorthologous gene displacement when the same function is fulfilled by nonorthologous proteins in two organisms. We identified 22 nonorthologous displacements and supplemented the set of orthologs with the respective M. genitalium genes. After examining the resulting list of 262 genes for possible functional redundancy and for the presence of apparently parasite-specific genes, 6 genes were removed. We suggest that the remaining 256 genes are close to the minimal gene set that is necessary and sufficient to sustain the existence of a modern-type cell. Most of the proteins encoded by the genes from the minimal set have eukaryotic or archaeal homologs but seven key proteins of DNA replication do not. We speculate that the last common ancestor of the three primary kingdoms had an RNA genome. Possibilities are explored to further reduce the minimal set to model a primitive cell that might have existed at a very early stage of life evolution.

Evidence for selection in evolution of alpha satellite DNA: the central role of CENP-B/pJ alpha binding region.

Conservation of DNA segments performing sequence-related functions is a landmark of selection and functional significance. Phylogenetic variability of alpha satellite and apparent absence of conserved regions calls its functional significance into question, even though sequence-specific alpha satellite-binding proteins pJ alpha and CENP-B have been discovered. Moreover, the function of pJ alpha is obscure and CENP-B binding satellite DNA, which is thought to participate in centromere formation, is found only in few species and not necessarily in all chromosomes. Analysis of alpha satellite evolution allows us to recognize the order in this variability. Here we report a new alpha satellite suprachromosomal family, which together with the four defined earlier, covers all known alpha satellite sequences. Although each family has its characteristic types of monomers, they all descend from two prototypes, A and B. We show that most differences between prototypes are concentrated in a short region (positions 35 to 51), which exists in two alternative states: it matches a binding site for pJ alpha in type A and the one for CENP-B in type B. Lower primates have only type A monomers whereas great apes have both A and B. The new family is formed by monomeric types almost identical to A and B prototypes, thus representing a living relic of alpha satellite. Analysis of these data shows that selection-driven evolution, rather than random fixation of mutations, formed the distinction between A and B types. To our knowledge, this is the first evidence for selection in any of the known satellite DNAs.

Organization of the 3'-terminal half of beet yellow stunt virus genome and implications for the evolution of closteroviruses.

The 3'-terminal half of the beet yellow stunt virus (BYSV) genome 10,545 nt, has been cloned and sequenced. The sequenced portion of the BYSV genome encompasses 10 open reading frames (ORFs) and 241 nt of the 3' untranslated region. The sequence spans, in the 5' to 3' direction, the C-terminal region of the replication-associated polyprotein gene (ORF 1a) which includes the set of motifs typical of helicases (HEL), the entire 53-kDa polymerase (RdRp) gene (ORF 1b), and genes encoding 30-kDa (ORF 2), 6-kDa (ORF 3), 66-kDa (ORF 4), 61-kDa (ORF 5), 25-kDa (ORF 6), 23.7-kDa (coat protein, CP) (ORF 7), 18-kDa (ORF 8), and 22-kDa (ORF 9) proteins. The double-stranded RNA "replicative form" of the BYSV was demonstrated to have a nontemplate G residue at the 3' terminus of the (+) strand. The RdRp of BYSV is presumably expressed via a +1 ribosomal frameshift. The five-gene module conserved among closteroviruses was identified in BYSV; it includes a gene array coding for a 6-kDa small hydrophobic protein, a 66-kDa homolog of the cellular HSP70 heat shock proteins, a 61-kDa protein, and a 25-kDa diverged copy of the CP followed by the CP gene itself. Phylogenetic analysis of the replication-associated HEL and RdRp domains as well as proteins from the five-gene module demonstrated the closest relationship between BYSV and two other closteroviruses, beet yellows (BYV) and citrus tristeza (CTV) viruses. Like CTV, the BYSV genome contains a 30-kDa protein gene between the RdRp and the 6-kDa protein genes, and like BYV it has only two genes downstream of the CP gene. The organization of the BYSV genome appears to be intermediate between BYV and CTV, which suggests that these three viruses might represent three distinct but probably close stages in the closterovirus evolution.