Genome Sequencing and Phylogenetic Analysis of 39 Human Parainfluenza Virus Type 1 Strains Isolated from 1997–2010

Thirty-nine human parainfluenza type 1 (HPIV-1) genomes were sequenced from samples collected in Milwaukee, Wisconsin from 1997–2010. Following sequencing, phylogenetic analyses of these sequences plus any publicly available HPIV-1 sequences (from GenBank) were performed. Phylogenetic analysis of the whole genomes, as well as individual genes, revealed that the current HPIV-1 viruses group into three different clades. Previous evolutionary studies of HPIV-1 in Milwaukee revealed that there were two genotypes of HPIV-1 co-circulating in 1991 (previously described as HPIV-1 genotypes C and D). The current study reveals that there are still two different HPIV-1 viruses co-circulating in Milwaukee; however, both groups of HPIV-1 viruses are derived from genotype C indicating that genotype D may no longer be in circulation in Milwaukee. Analyses of genetic diversity indicate that while most of the genome is under purifying selection some regions of the genome are more tolerant of mutation. In the 40 HPIV-1 genomes sequenced in this study, the nucleotide sequence of the L gene is the most conserved while the sequence of the P gene is the most variable. Over the entire protein coding region of the genome, 81 variable amino acid residues were observed and as with nucleotide diversity, the P protein seemed to be the most tolerant of mutation (and contains the greatest proportion of non-synonymous to synonymous substitutions) while the M protein appears to be the least tolerant of amino acid substitution.     

Genome-Wide Comparisons of Phylogenetic Similarities between Partial Genomic Regions and the Full-Length Genome in Hepatitis E Virus Genotyping

Besides the complete genome, different partial genomic sequences of Hepatitis E virus (HEV) have been used in genotyping studies, making it difficult to compare the results based on them. No commonly agreed partial region for HEV genotyping has been determined. In this study, we used a statistical method to evaluate the phylogenetic performance of each partial genomic sequence from a genome wide, by comparisons of evolutionary distances between genomic regions and the full-length genomes of 101 HEV isolates to identify short genomic regions that can reproduce HEV genotype assignments based on full-length genomes. Several genomic regions, especially one genomic region at the 3′-terminal of the papain-like cysteine protease domain, were detected to have relatively high phylogenetic correlations with the full-length genome. Phylogenetic analyses confirmed the identical performances between these regions and the full-length genome in genotyping, in which the HEV isolates involved could be divided into reasonable genotypes. This analysis may be of value in developing a partial sequence-based consensus classification of HEV species.     

Discovery and characterization of mammalian endogenous parvoviruses

Public databases of nucleotide sequences contain exponentially increasing amounts of sequence data from mammalian genomes. Through the use of large-scale bioinformatic screening for sequences homologous to exogenous mammalian viruses, we found several sequences related to human and animal parvoviruses (PVs) in the Parvovirus and Dependovirus genera within genomes of several mammals, including rats, wallabies, opossums, guinea pigs, hedgehogs, African elephants, and European rabbits. However, phylogenetic analysis of these endogenous parvovirus (EnPV) sequences demonstrated substantial genetic divergence from exogenous mammalian PVs characterized to date. Entire nonstructural and capsid gene sequences of a novel EnPV were amplified and genetically characterized from rat (Rattus norvegicus) genomic DNA. Rat EnPV sequences were most closely related to members of the genus Parvovirus, with >70% and 65% amino acid identities to nonstructural and capsid proteins of canine parvovirus, respectively. Integration of EnPV into chromosome 5 of rats was confirmed by PCR cloning and sequence analysis of the viral and chromosomal junctions. Using inverse PCR, we determined that the rat genome contains a single copy of rat EnPV. Considering mammalian phylogeny, we estimate that EnPV integrated into the rat genome less than 30 million years ago. Comparative phylogenetic analysis done using all known representative exogenous parvovirus (ExPV) and EnPV sequences showed two major genetic groups of EnPVs, one genetically more similar to genus Parvovirus and the other genetically more similar to the genus Dependovirus. The full extent of the genetic diversity of parvoviruses that have undergone endogenization during evolution of mammals and other vertebrates will be recognized only once complete genomic sequences from a wider range of classes, orders, and species of animals become available.     

A quantitative genotype algorithm reflecting H5N1 Avian influenza niches

MOTIVATION: Computational genotyping analyses are critical for characterizing molecular evolutionary footprints, thus providing important information for designing the strategies of influenza prevention and control. Most of the current methods that are available are based on multiple sequence alignment and phylogenetic tree construction, which are time consuming and limited by the number of taxa. Arbitrarily defining genotypes further complicates the interpretation of genotyping results. METHODS: In this study, we describe a quantitative influenza genotyping algorithm based on the theory of quasispecies. First, the complete composition vector (CCV) was utilized to calculate the pairwise evolutionary distance between genotypes. Next, Hierarchical Bayesian Modeling using the Gibbs Sampling algorithm was applied to identify the segment genotype threshold, which is used to identify influenza segment genotype through a modularity calculation. The viral genotype was defined by combining eight segment genotypes based on the genetic reassortment feature of influenza A viruses. RESULTS: We applied this method for H5N1 avian influenza viruses and identified 107 niches among 283 viruses with a complete genome set. The diversity of viral genotypes, and their correlation with geographic locations suggests that these viruses form local niches after being introduced to a new ecological environment through poultry trade or bird migration. This novel method allows us to define genotypes in a robust, quantitative as well as hierarchical manner. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Comparative genomics of foot-and-mouth disease virus

Here we present complete genome sequences, including a comparative analysis, of 103 isolates of foot-and-mouth disease virus (FMDV) representing all seven serotypes and including the first complete sequences of the SAT1 and SAT3 genomes. The data reveal novel highly conserved genomic regions, indicating functional constraints for variability as well as novel viral genomic motifs with likely biological relevance. Previously undescribed invariant motifs were identified in the 5' and 3' untranslated regions (UTR), as was tolerance for insertions/deletions in the 5' UTR. Fifty-eight percent of the amino acids encoded by FMDV isolates are invariant, suggesting that these residues are critical for virus biology. Novel, conserved sequence motifs with likely functional significance were identified within proteins L(pro), 1B, 1D, and 3C. An analysis of the complete FMDV genomes indicated phylogenetic incongruities between different genomic regions which were suggestive of interserotypic recombination. Additionally, a novel SAT virus lineage containing nonstructural protein-encoding regions distinct from other SAT and Euroasiatic lineages was identified. Insights into viral RNA sequence conservation and variability and genetic diversity in nature will likely impact our understanding of FMDV infections, host range, and transmission.     

The Complete Genome Sequence of Natrinema sp. J7-2, a Haloarchaeon Capable of Growth on Synthetic Media without Amino Acid Supplements

Natrinema sp. J7-2 is an extreme haloarchaeon capable of growing on synthetic media without amino acid supplements. Here we report the complete genome sequence of Natrinema sp. J7-2 which is composed of a 3,697,626-bp chromosome and a 95,989-bp plasmid pJ7-I. This is the first complete genome sequence of a member of the genus Natrinema. We demonstrate that Natrinema sp. J7-2 can use gluconate, glycerol, or acetate as the sole carbon source and that its genome encodes complete metabolic pathways for assimilating these substrates. The biosynthetic pathways for all 20 amino acids have been reconstructed, and we discuss a possible evolutionary relationship between the haloarchaeal arginine synthetic pathway and the bacterial lysine synthetic pathway. The genome harbors the genes for assimilation of ammonium and nitrite, but not nitrate, and has a denitrification pathway to reduce nitrite to N(2)O. Comparative genomic analysis suggests that most sequenced haloarchaea employ the TrkAH system, rather than the Kdp system, to actively uptake potassium. The genomic analysis also reveals that one of the three CRISPR loci in the Natrinema sp. J7-2 chromosome is located in an integrative genetic element and is probably propagated via horizontal gene transfer (HGT). Finally, our phylogenetic analysis of haloarchaeal genomes provides clues about evolutionary relationships of haloarchaea.     

The Effect of Vaccination on the Evolution and Population Dynamics of Avian Paramyxovirus-1

Newcastle Disease Virus (NDV) is a pathogenic strain of avian paramyxovirus (aPMV-1) that is among the most serious of disease threats to the poultry industry worldwide. Viral diversity is high in aPMV-1; eight genotypes are recognized based on phylogenetic reconstruction of gene sequences. Modified live vaccines have been developed to decrease the economic losses caused by this virus. Vaccines derived from avirulent genotype II strains were developed in the 1950s and are in use globally, whereas Australian strains belonging to genotype I were developed as vaccines in the 1970s and are used mainly in Asia. In this study, we evaluated the consequences of attenuated live virus vaccination on the evolution of aPMV-1 genotypes. There was phylogenetic incongruence among trees based on individual genes and complete coding region of 54 full length aPMV-1 genomes, suggesting that recombinant sequences were present in the data set. Subsequently, five recombinant genomes were identified, four of which contained sequences from either genotype I or II. The population history of vaccine-related genotype II strains was distinct from other aPMV-1 genotypes; genotype II emerged in the late 19^th century and is evolving more slowly than other genotypes, which emerged in the 1960s. Despite vaccination efforts, genotype II viruses have experienced constant population growth to the present. In contrast, other contemporary genotypes showed population declines in the late 1990s. Additionally, genotype I and II viruses, which are circulating in the presence of homotypic vaccine pressure, have unique selection profiles compared to nonvaccine-related strains. Collectively, these data show that vaccination with live attenuated viruses has changed the evolution of aPMV-1 by maintaining a large effective population size of a vaccine-related genotype, allowing for coinfection and recombination of vaccine and wild type strains, and by applying unique selective pressures on viral glycoproteins.     

Genotypic and pathotypic characterization of Newcastle disease viruses from India

Newcastle disease virus (NDV) is an avian paramyxovirus that causes significant economic losses to the poultry industry in most parts of the world. The susceptibility of a wide variety of avian species coupled with synanthropic bird reservoirs has contributed to the vast genomic diversity of this virus as well as diagnostic failures. Since the first panzootic in 1926, Newcastle disease (ND) became enzootic in India with recurrent outbreaks in multiple avian species. The genetic characteristics of circulating strains in India, however, are largely unknown. To understand the nature of NDV genotypes in India, we characterized two representative strains isolated 13 years apart from a chicken and a pigeon by complete genome sequence analysis and pathotyping. The viruses were characterized as velogenic by pathogenicity indices devised to distinguish these strains. The genome length was 15,186 nucleotides (nt) and consisted of six non-overlapping genes, with conserved and complementary 3' leader and 5' trailer regions, conserved gene starts, gene stops, and intergenic sequences similar to those in avian paramyxovirus 1 (APMV-1) strains. Matrix gene sequence analysis grouped the pigeon isolate with APMV-1 strains. Phylogeny based on the fusion (F), and hemagglutinin (HN) genes and complete genome sequence grouped these viruses into genotype IV. Genotype IV strains are considered to have "died out" after the first panzootic (1926-1960) of ND. But, our results suggest that there is persistence of genotype IV strains in India.     

Genome analysis with the conditional multinomial distribution profile

The focus of the research is on the analysis of genome sequences. Based on the inter-nucleotide distance sequence, we propose the conditional multinomial distribution profile for the complete genomic sequence. These profiles can be used to define a very simple, computationally efficient, alignment-free, distance measure that reflects the evolutionary relationships between genomic sequences. We use this distance measure to classify chromosomes according to species of origin, to build the phylogenetic tree of 24 complete genome sequences of coronaviruses. Our results demonstrate the new method is powerful and efficient.     

新分离的副粘病毒Tianjin株的全基因组序列分析

副粘病毒Tianjin株是一株对普通棉耳狨猴具有高致病性,并可能与人类下呼吸道感染密切相关的毒株.为了明确其基因结构、变异特点及种系进化地位,采用RT PCR、测序和拼接,获得了副粘病毒Tianjin株全基因组序列,与GenBank登录的副粘病毒科7个属和尚未分类的28株病毒及7株仙台病毒代表株,进行同源性比较及系统进化分析.结果表明,副粘病Tianjin株属于副粘病毒科、副粘病毒亚科、呼吸道病毒属,与仙台病毒关系最近.基因组全长及组成规律与仙台病毒相似,只是L基因末尾A15240C变异而使L蛋白增加了一个谷氨酸残基.副粘病毒Tianjin株存在440个独特的核苷酸变异位点,导致110个氨基酸残基的改变,系统进化上构成独立的分支.副粘病毒Tianjin株在基因组序列、宿主亲嗜性和致病性等方面与已知仙台病毒存在较大的差异,可能代表仙台病毒的一个新基因型.     

Evidence of two Lyssavirus phylogroups with distinct pathogenicity and immunogenicity

The genetic diversity of representative members of the Lyssavirus genus (rabies and rabies-related viruses) was evaluated using the gene encoding the transmembrane glycoprotein involved in the virus-host interaction, immunogenicity, and pathogenicity. Phylogenetic analysis distinguished seven genotypes, which could be divided into two major phylogroups having the highest bootstrap values. Phylogroup I comprises the worldwide genotype 1 (classic Rabies virus), the European bat lyssavirus (EBL) genotypes 5 (EBL1) and 6 (EBL2), the African genotype 4 (Duvenhage virus), and the Australian bat lyssavirus genotype 7. Phylogroup II comprises the divergent African genotypes 2 (Lagos bat virus) and 3 (Mokola virus). We studied immunogenic and pathogenic properties to investigate the biological significance of this phylogenetic grouping. Viruses from phylogroup I (Rabies virus and EBL1) were found to be pathogenic for mice when injected by the intracerebral or the intramuscular route, whereas viruses from phylogroup II (Mokola and Lagos bat viruses) were only pathogenic by the intracerebral route. We showed that the glycoprotein R333 residue essential for virulence was naturally replaced by a D333 in the phylogroup II viruses, likely resulting in their attenuated pathogenicity. Moreover, cross-neutralization distinguished the same phylogroups. Within each phylogroup, the amino acid sequence of the glycoprotein ectodomain was at least 74% identical, and antiglycoprotein virus-neutralizing antibodies displayed cross-neutralization. Between phylogroups, the identity was less than 64.5% and the cross-neutralization was absent, explaining why the classical rabies vaccines (phylogroup I) cannot protect against lyssaviruses from phylogroup II. Our tree-axial analysis divided lyssaviruses into two phylogroups that more closely reflect their biological characteristics than previous serotypes and genotypes.     

Previously unknown and highly divergent ssDNA viruses populate the oceans

Single-stranded DNA (ssDNA) viruses are economically important pathogens of plants and animals, and are widespread in oceans; yet, the diversity and evolutionary relationships among marine ssDNA viruses remain largely unknown. Here we present the results from a metagenomic study of composite samples from temperate (Saanich Inlet, 11 samples; Strait of Georgia, 85 samples) and subtropical (46 samples, Gulf of Mexico) seawater. Most sequences (84%) had no evident similarity to sequenced viruses. In total, 608 putative complete genomes of ssDNA viruses were assembled, almost doubling the number of ssDNA viral genomes in databases. These comprised 129 genetically distinct groups, each represented by at least one complete genome that had no recognizable similarity to each other or to other virus sequences. Given that the seven recognized families of ssDNA viruses have considerable sequence homology within them, this suggests that many of these genetic groups may represent new viral families. Moreover, nearly 70% of the sequences were similar to one of these genomes, indicating that most of the sequences could be assigned to a genetically distinct group. Most sequences fell within 11 well-defined gene groups, each sharing a common gene. Some of these encoded putative replication and coat proteins that had similarity to sequences from viruses infecting eukaryotes, suggesting that these were likely from viruses infecting eukaryotic phytoplankton and zooplankton.     

Complete sequences of two highly divergent european isolates of TT virus

TT virus is a virus distantly related to the Circoviridae family. We report here the complete genome characterization of two European human isolates (T3PB and TUPB) using a new and simple protocol for sequencing GC-rich genomic regions. Sequence analysis confirmed the existence of two major ORFs, of a CAV-like VP2 motif in ORF2 and of potential stem-loop structures in non-coding regions. Phylogenetic analyses based on complete genomic sequences of human isolates suggested that three different lineages exist at least. The first lineage includes genotypes 1, 2, and 3, and two other lineages include viruses related to the Japanese SANBAN and to the North American TUS01 isolates respectively. Sequence comparison made it possible to assign strain T3PB to genotype 3, and strain TUPB to the TUS01 group. Consequently, this study reports the first full-length sequence of a genotype 3 isolate and demonstrates that viruses belonging to the TUS01 lineage are present in the Old Word.     

Comparative Analysis of the Peanut Witches'-Broom Phytoplasma Genome Reveals Horizontal Transfer of Potential Mobile Units and Effectors

Phytoplasmas are a group of bacteria that are associated with hundreds of plant diseases. Due to their economical importance and the difficulties involved in the experimental study of these obligate pathogens, genome sequencing and comparative analysis have been utilized as powerful tools to understand phytoplasma biology. To date four complete phytoplasma genome sequences have been published. However, these four strains represent limited phylogenetic diversity. In this study, we report the shotgun sequencing and evolutionary analysis of a peanut witches''-broom (PnWB) phytoplasma genome. The availability of this genome provides the first representative of the 16SrII group and substantially improves the taxon sampling to investigate genome evolution. The draft genome assembly contains 13 chromosomal contigs with a total size of 562,473 bp, covering ∼90% of the chromosome. Additionally, a complete plasmid sequence is included. Comparisons among the five available phytoplasma genomes reveal the differentiations in gene content and metabolic capacity. Notably, phylogenetic inferences of the potential mobile units (PMUs) in these genomes indicate that horizontal transfer may have occurred between divergent phytoplasma lineages. Because many effectors are associated with PMUs, the horizontal transfer of these transposon-like elements can contribute to the adaptation and diversification of these pathogens. In summary, the findings from this study highlight the importance of improving taxon sampling when investigating genome evolution. Moreover, the currently available sequences are inadequate to fully characterize the pan-genome of phytoplasmas. Future genome sequencing efforts to expand phylogenetic diversity are essential in improving our understanding of phytoplasma evolution.     

Emergence of genotype I of Japanese encephalitis virus as the dominant genotype in Asia

Japanese encephalitis virus (JEV), a mosquito-borne zoonotic pathogen, is one of the major causes of viral encephalitis worldwide. Previous phylogenetic studies based on the envelope protein indicated that there are four genotypes, and surveillance data suggest that genotype I is gradually replacing genotype III as the dominant strain. Here we report an evolutionary analysis based on 98 full-length genome sequences of JEV, including 67 new samples isolated from humans, pigs, mosquitoes, midges. and bats in affected areas. To investigate the relationships between the genotypes and the significance of genotype I in recent epidemics, we estimated evolutionary rates, ages of common ancestors, and population demographics. Our results indicate that the genotypes diverged in the order IV, III, II, and I and that the genetic diversity of genotype III has decreased rapidly while that of genotype I has increased gradually, consistent with its emergence as the dominant genotype.     

The ITS1-5.8S-ITS2 sequence region in the Musaceae: structure, diversity and use in molecular phylogeny

Genes coding for 45S ribosomal RNA are organized in tandem arrays of up to several thousand copies and contain 18S, 5.8S and 26S rRNA units separated by internal transcribed spacers ITS1 and ITS2. While the rRNA units are evolutionary conserved, ITS show high level of interspecific divergence and have been used frequently in genetic diversity and phylogenetic studies. In this work we report on the structure and diversity of the ITS region in 87 representatives of the family Musaceae. We provide the first detailed information on ITS sequence diversity in the genus Musa and describe the presence of more than one type of ITS sequence within individual species. Both Sanger sequencing of amplified ITS regions and whole genome 454 sequencing lead to similar phylogenetic inferences. We show that it is necessary to identify putative pseudogenic ITS sequences, which may have negative effect on phylogenetic reconstruction at lower taxonomic levels. Phylogenetic reconstruction based on ITS sequence showed that the genus Musa is divided into two distinct clades--Callimusa and Australimusa and Eumusa and Rhodochlamys. Most of the intraspecific banana hybrids analyzed contain conserved parental ITS sequences, indicating incomplete concerted evolution of rDNA loci. Independent evolution of parental rDNA in hybrids enables determination of genomic constitution of hybrids using ITS. The observation of only one type of ITS sequence in some of the presumed interspecific hybrid clones warrants further study to confirm their hybrid origin and to unravel processes leading to evolution of their genomes.     

From gene trees to organismal phylogeny in prokaryotes: the case of the gamma-Proteobacteria

The rapid increase in published genomic sequences for bacteria presents the first opportunity to reconstruct evolutionary events on the scale of entire genomes. However, extensive lateral gene transfer (LGT) may thwart this goal by preventing the establishment of organismal relationships based on individual gene phylogenies. The group for which cases of LGT are most frequently documented and for which the greatest density of complete genome sequences is available is the gamma-Proteobacteria, an ecologically diverse and ancient group including free-living species as well as pathogens and intracellular symbionts of plants and animals. We propose an approach to multigene phylogeny using complete genomes and apply it to the case of the gamma-Proteobacteria. We first applied stringent criteria to identify a set of likely gene orthologs and then tested the compatibilities of the resulting protein alignments with several phylogenetic hypotheses. Our results demonstrate phylogenetic concordance among virtually all (203 of 205) of the selected gene families, with each of the exceptions consistent with a single LGT event. The concatenated sequences of the concordant families yield a fully resolved phylogeny. This topology also received strong support in analyses aimed at excluding effects of heterogeneity in nucleotide base composition across lineages. Our analysis indicates that single-copy orthologous genes are resistant to horizontal transfer, even in ancient bacterial groups subject to high rates of LGT. This gene set can be identified and used to yield robust hypotheses for organismal phylogenies, thus establishing a foundation for reconstructing the evolutionary transitions, such as gene transfer, that underlie diversity in genome content and organization.     

Coenzymes, viruses and the RNA world

The results of a detailed bioinformatic search for ribonucleotidyl coenzyme biosynthetic sequences in DNA- and RNA viral genomes are presented. No RNA viral genome sequence available as of April 2011 appears to encode for sequences involved in coenzyme biosynthesis. In both single- and double-stranded DNA viruses a diverse array of coenzyme biosynthetic genes has been identified, but none of the viral genomes examined here encodes for a complete pathway. Although our conclusions may be constrained by the unexplored diversity of viral genomes and the biases in the construction of viral genome databases, our results do not support the possibility that RNA viruses are direct holdovers from an ancient RNA/protein world. Extrapolation of our results to evolutionary epochs prior to the emergence of DNA genomes suggest that during those early stages living entities may have depended on discontinuous genetic systems consisting of multiple small-size RNA sequences.     

Complete nucleotide sequence of the <Emphasis Type="Italic">Cryptomeria japonica</Emphasis>D. Don. chloroplast genome and comparative chloroplast genomics: diversified genomic structure of coniferous species

Background  

The recent determination of complete chloroplast (cp) genomic sequences of various plant species has enabled numerous comparative analyses as well as advances in plant and genome evolutionary studies. In angiosperms, the complete cp genome sequences of about 70 species have been determined, whereas those of only three gymnosperm species, Cycas taitungensis, Pinus thunbergii, and Pinus koraiensis have been established. The lack of information regarding the gene content and genomic structure of gymnosperm cp genomes may severely hamper further progress of plant and cp genome evolutionary studies. To address this need, we report here the complete nucleotide sequence of the cp genome of Cryptomeria japonica, the first in the Cupressaceae sensu lato of gymnosperms, and provide a comparative analysis of their gene content and genomic structure that illustrates the unique genomic features of gymnosperms.