Evolution and diversity of the Microviridae viral family through a collection of 81 new complete genomes assembled from virome reads

Recent studies suggest that members of the Microviridae (a family of ssDNA bacteriophages) might play an important role in a broad spectrum of environments, as they were found in great number among the viral fraction from seawater and human gut samples. 24 completely sequenced Microviridae have been described so far, divided into three distinct groups named Microvirus, Gokushovirinae and Alpavirinae, this last group being only composed of prophages. In this study, we present the analysis of 81 new complete Microviridae genomes, assembled from viral metagenomes originating from various ecosystems. The phylogenetic analysis of the core genes highlights the existence of four groups, confirming the three sub-families described so far and exhibiting a new group, named Pichovirinae. The genomic organizations of these viruses are strikingly coherent with their phylogeny, the Pichovirinae being the only group of this family with a different organization of the three core genes. Analysis of the structure of the major capsid protein reveals the presence of mushroom-like insertions conserved within all the groups except for the microviruses. In addition, a peptidase gene was found in 10 Microviridae and its analysis indicates a horizontal gene transfer that occurred several times between these viruses and their bacterial hosts. This is the first report of such gene transfer in Microviridae. Finally, searches against viral metagenomes revealed the presence of highly similar sequences in a variety of biomes indicating that Microviridae probably have both an important role in these ecosystems and an ancient origin.     

Experimental evolution of viruses: Microviridae as a model system

φX174 was developed as a model system for experimental studies of evolution because of its small genome size and ease of cultivation. It has been used extensively to address statistical questions about the dynamics of adaptive evolution. Molecular changes seen during experimental evolution of φX174 under a variety of conditions were compiled from 10 experiments comprising 58 lineages, where whole genomes were sequenced. A total of 667 substitutions was seen. Parallel evolution was rampant, with over 50 per cent of substitutions occurring at sites with three or more events. Comparisons of experimentally evolved sites to variation seen among wild phage suggest that at least some of the adaptive mechanisms seen in the laboratory are relevant to adaptation in nature. Elucidation of these mechanisms is aided by the availability of capsid and pro-capsid structures for φX174 and builds on years of genetic studies of the phage life history.     

Possible finger structure in gene A protein of Microviridae

Microvirid phages alpha 3 and phi K encode for A protein which functions in initiation of the viral DNA synthesis. By nucleotide sequencing analysis, we have found that each gene A protein has 'finger motif structure' which conserves two cysteine and histidine residues similar to the consensus sequence deduced from more than thirty finger motifs reported from many eukaryotic regulatory proteins. In closely related phages, phi X174 and G4, we have detected the same structures in their gene A proteins also.     

Shiga toxins even when different are encoded at identical positions in the genomes of related temperate bacteriophages

The nucleotide sequence of an 11,142-bp region including the stx ₂ operon in the genome of the temperate bacteriophage 933W in the EDL933 strain of Escherichia coli O157 was determined and compared to the respective regions derived from other lambdoid bacteriophages. In phage 933W, a region of ORFs interlinked by overlapping start-stop codons (ATGA) was detected preceding the toxin gene. These ORFs show a high degree of sequence identity to genes of the nin region of phage λ. Immediately downstream of these nin genes we identified an ORF that may code for an antiterminator similar to the λ Q protein. It is concluded that toxin expression is directly associated with the initiation of cell lysis. Downstream of the stx ₂ operon we identified an ORF that is homologous to the holin gene S of bacteriophage PA-2. PCR primers were designed, which, based on a comparison of the phage sequences, appeared to be common to both stx ₁ - and stx ₂ -harbouring phages. However, only seven of the 22 STEC strains investigated from serogroups O157, O26, O103 and O111 yielded the expected PCR amplification product. The data reported here may be useful in developing new strategies for inhibiting the expression of Stx and for developing universal diagnostic primers for use in tracking the origin and evolution of Shiga toxins and the phages that carry them. Received: 24 February 1999 / Accepted: 7 September 1999     

High-throughput sequencing of six bamboo chloroplast genomes: phylogenetic implications for temperate woody bamboos (Poaceae: Bambusoideae)

Background

Bambusoideae is the only subfamily that contains woody members in the grass family, Poaceae. In phylogenetic analyses, Bambusoideae, Pooideae and Ehrhartoideae formed the BEP clade, yet the internal relationships of this clade are controversial. The distinctive life history (infrequent flowering and predominance of asexual reproduction) of woody bamboos makes them an interesting but taxonomically difficult group. Phylogenetic analyses based on large DNA fragments could only provide a moderate resolution of woody bamboo relationships, although a robust phylogenetic tree is needed to elucidate their evolutionary history. Phylogenomics is an alternative choice for resolving difficult phylogenies.

Methodology/Principal Findings

Here we present the complete nucleotide sequences of six woody bamboo chloroplast (cp) genomes using Illumina sequencing. These genomes are similar to those of other grasses and rather conservative in evolution. We constructed a phylogeny of Poaceae from 24 complete cp genomes including 21 grass species. Within the BEP clade, we found strong support for a sister relationship between Bambusoideae and Pooideae. In a substantial improvement over prior studies, all six nodes within Bambusoideae were supported with ≥0.95 posterior probability from Bayesian inference and 5/6 nodes resolved with 100% bootstrap support in maximum parsimony and maximum likelihood analyses. We found that repeats in the cp genome could provide phylogenetic information, while caution is needed when using indels in phylogenetic analyses based on few selected genes. We also identified relatively rapidly evolving cp genome regions that have the potential to be used for further phylogenetic study in Bambusoideae.

Conclusions/Significance

The cp genome of Bambusoideae evolved slowly, and phylogenomics based on whole cp genome could be used to resolve major relationships within the subfamily. The difficulty in resolving the diversification among three clades of temperate woody bamboos, even with complete cp genome sequences, suggests that these lineages may have diverged very rapidly.     

M-MuLV-induced leukemogenesis: Integration and structure of recombinant proviruses in tumors

M-MuLV-specific DNA probes were used to establish the state of integration and amplification of recombinant proviral sequences in Moloney virus-induced tumors of Balb/Mo, Balb/c and 129 mice. The somatically acquired viral sequences contain both authentic M-MuLV genomes and recombinants of M-MuLV with endogenous viral sequences. All reintegrated genomes carry long terminal repeat (LTR) sequences at both termini of their genome. In the preleukemic stage a large population of cells exhibiting a random distribution of reintegrated M-MuLV genomes are seen, but during outgrowth of the tumor, selection of cells occurs leaving one or a few clonal descendants in the outgrown tumor. In this latter stage recombinant genomes can be detected. Although these recombinants constitute a heterogeneous group of proviruses, characteristic molecular markers are conserved among many individual proviral recombinants, lending credence to the notion that a certain recombinant structure is a prerequisite for the onset of neoplasia. The structure of these recombinants shows close structural similarities to the previously described mink cell focus-inducing (MCF)-type viruses.     

Structural basis of sequestration of the anti-Shine-Dalgarno sequence in the Bacteroidetes ribosome

Genomic studies have indicated that certain bacterial lineages such as the Bacteroidetes lack Shine-Dalgarno (SD) sequences, and yet with few exceptions ribosomes of these organisms carry the canonical anti-SD (ASD) sequence. Here, we show that ribosomes purified from Flavobacterium johnsoniae, a representative of the Bacteroidetes, fail to recognize the SD sequence of mRNA in vitro. A cryo-electron microscopy structure of the complete 70S ribosome from F. johnsoniae at 2.8 Å resolution reveals that the ASD is sequestered by ribosomal proteins bS21, bS18 and bS6, explaining the basis of ASD inhibition. The structure also uncovers a novel ribosomal protein—bL38. Remarkably, in F. johnsoniae and many other Flavobacteriia, the gene encoding bS21 contains a strong SD, unlike virtually all other genes. A subset of Flavobacteriia have an alternative ASD, and in these organisms the fully complementary sequence lies upstream of the bS21 gene, indicative of natural covariation. In other Bacteroidetes classes, strong SDs are frequently found upstream of the genes for bS21 and/or bS18. We propose that these SDs are used as regulatory elements, enabling bS21 and bS18 to translationally control their own production.     

Molecular cloning of MuLV proviruses integrated into the genome of mouse erythroleukemia cells. I. Characteristics of endogenous proviruses

The library of genes was obtained from erythroleukemic AKR cells (C-1), that were maintained as suspension culture. Thirty four clones that had homology with 60-70S RNA of Rauscher Leukemia virus (RLV) were separated from this library. The restriction mapping was carried out with 14 clones, that contained most extensive proviral sequences. One clone (107) contains proviral sequences that are derived from one of the components of the RLV complex. The other 13 clones contain sequences of endogenous xenotropic viruses. The endogenous retroviral sequences obtained differ in restrictive maps from proviruses of ecotropic and xenotropic infectious endogenous MuLV and, apparently, might be attributed as non-inducible infectious xenotropic MuLV of class III. Some of the cloned retroviral sequences had symmetrical structure, that is typical for integrated proviruses, i. e. these sequences were separated from flanking cellular ones by long terminal repeats. All investigated retroviral sequences are deletion mutants of MuLV proviruses. It was shown that the inner regions of proviruses diverged more than the long terminal repeats. The expression of the main inner MuLV polypeptide (p30) was detected in NIH 3T3 cells, transfected with DNA of some clones.     

Human lymphocyte glucocorticoid receptors reside mainly in the cytoplasm

Karyoplasts prepared from human lymphocytes by treatment with cytochalasin B or lymphocytes permeabilized with digitonin contain only 25-40% of the total cellular glucocorticoid binding activity. The presented results thus support the original concept that unfilled steroid hormone receptors are mainly cytoplasmic.