Time-calibrated phylogenomics of the porcine epidemic diarrhea virus: genome-wide insights into the spatio-temporal dynamics

Porcine epidemic diarrhea virus (PEDV), the causative agent of porcine epidemic diarrhea (PED), has led to tremendous economic losses in the global swine industry. Although the phylogeny of PEDV has been investigated extensively at the molecular level, there was no time-calibrated phylogenomic study on the virus. To improve insight into this topic, we analyzed 138 published genome sequences using the Bayesian coalescent analyses as well as Bayesian inferences and maximum likelihood methods. All of the global PEDV isolates were divided into six groups, except for one unclassified isolate. Of the six groups, Groups 1–5 comprised pandemic viruses while the remaining Group 6 contained classical isolates. Interestingly, the two clades, both pandemic and classical, consisted of clade-specific amino acid sequences in five genes: ORF1a, ORF1b, S, ORF3, and N. Within the pandemic clade, Group 1 and Group 2 originated from North America, whereas Group 3–Group 5 were derived from Asia. In Group 2, there was a common origin of S INDEL isolates. Within each group, there was no apparent association between temporal or geographic origin and heterogeneity of PEDVs. Our findings also showed that the PEDV virus evolved at a rate of 3.38?×?10^?4 substitutions/site/year, and the most recent common ancestor of the virus emerged 75.9 years ago. Our Bayesian skyline plot analysis indicated that the PEDV had maintained constant effective population size excluding only a short period, around 2012, when a valley shaped decline in the effective number of infections occurred.     

Time-Calibrated Phylogenomics of the Classical Swine Fever Viruses: Genome-Wide Bayesian Coalescent Approach

The phylogeny of classical swine fever virus (CSFV), the causative agent of classical swine fever (CSF), has been investigated extensively. However, no evolutionary research has been performed using the whole CSFV genome. In this study, we used 37 published genome sequences to investigate the time-calibrated phylogenomics of CSFV. In phylogenomic trees based on Bayesian inference (BI) and Maximum likelihood (ML), the 37 isolates were categorized into five genetic types (1.1, 1.2, 2.1, 2.3, and 3.4). Subgenotype 1.1 is divided into 3 groups and 1 unclassified isolate, 2.1 into 4 groups, 2.3 into 2 groups and 1 unclassified isolate, and subgenotype 1.2 and 3.4 consisted of one isolate each. We did not observe an apparent temporal or geographical relationship between isolates. Of the 14 genomic regions, NS4B showed the most powerful phylogenetic signal. Results of this evolutionary study using Bayesian coalescent approach indicate that CSFV has evolved at a rate of 13×.010^-4 substitutions per site per year. The most recent common ancestor of CSFV appeared 2770.2 years ago, which was about 8000 years after pig domestication. The effective population size of CSFV underwent a slow increase until the 1950s, after which it has remained constant.     

Molecular evolution of poxviruses

Previous restriction fragment length polymorphism analysis divided variola virus (VARV) strains into two subtypes, one of which included West African and South American isolates. This allowed a dating to be introduced for the first time in estimation of the VARV evolution rate. The results were used to analyze the molecular evolution of the total family Poxviridae. Comparisons of the known nucleotide sequences were performed for the extended conserved central genome region in 42 orthopoxvirus strains and for the eight genes of multisubunit RNA polymerase in 65 viruses belonging to various genera of the family Poxviridae. Using the Bayesian dating method, the mutation accumulation rate of poxviruses was estimated at (1.7–8.8) × 10^?6 nucleotide substitutions per site per year. Computations showed that the modern poxvirus genera started diverging from an ancestral virus more than 200 thousand years ago and that an ancestor of the genus Orthopoxvirus emerged 131 ± 45 thousand years ago. The other genera of mammalian poxviruses with a low GC content diverged approximately 110–90 thousand years ago. The independent evolution of VARV started 3.4 ± 0.8 thousand years ago. It was shown with the example of VARV and the monkeypox virus (MPXV) that divergent evolution of these orthopoxviruses started and the West African subtypes of VARV and MPXV were formed as geographical conditions changed to allow isolation of West African animals from other African regions.     

Time scale of Poxvirus evolution

Unlike in vertebrates and RNA viruses, the molecular clock has not been estimated so far for DNA viruses. The extended conserved central region (102 kb) of the orthopoxvirus genome and the DNA polymerase gene (3 kb) were analyzed in viruses representing several genera of the family Poxviridae. Analysis was based on the known dating of the variola virus (VARV) transfer from Western Africa to South America and previous data on the phylogenetic relatedness of modern West African and South American isolates of VARV. The mutation accumulation rate was for the first time estimated for these DNA viruses at (0.9–1.2) × 10^?6 substitutions per site per year. It was assumed that poxviruses diverged from an ancestor approximately 500,000 years ago to form the recent species and that the ancestor of the genus Orthopoxvirus emerged approximately 300,000 years ago and gave origin to the modern species approximately 14,000 years ago.     

Molecular characterization of a highly pathogenetic porcine reproductive and respiratory syndrome virus variant in Hubei, China

Porcine reproductive and respiratory syndrome virus (PRRSV) has been recognized as one of the most important pathogens of pigs throughout the world. In 2006, more than 10 provinces of China have experienced an epizootic outbreak of pig diseases characterized by high fever, reddened skin and high morbidity and mortality. From June 2006 to April 2007, we have investigated some clinical samples in Hubei province by RT-PCR and cloned several major genes, N, GP5 and NSP2 gene, shown in this study. Phylogenetic analysis of these genes revealed that the highly pathogenic PRRSV variant, ZB, was responsible for 2006 emergent outbreak of pig disease in Hubei province similar with those variants isolated from other provinces in China in 2006, and belongs to the NA-type PRRSV. In the PRRSV variants, the N and GP5 shear about 90% identity with prototypic ATCC VR-2332 and some typical NA-type Chinese isolates, except the 2850bp NSP2 gene (only shares 65% identity with ATCC VR-2332). But they all shear more than and 97% identity with other highly pathogenetic Chinese PRRSV strains. Additionally, there are extensive amino acid (aa) mutations in the GP5 protein and 2 deletions in the Nsp2 protein when compared with the previous isolates. Most of the variants found in 2006 epizootic outbreak of pig diseases in China were the farthest variants from the typical NA-type PRRSV in phylogenetic distance, and these diversities may be responsible for the differences in the pathogenicity observed between these variants and original Chinese PRRSV strains.     

Complete genome sequence of porcine reproductive and respiratory syndrome virus isolated from piglet stool samples

WUH4 is a highly pathogenic North American porcine reproductive and respiratory syndrome virus (PRRSV). Unlike previous PRRSV isolates, which were mainly recovered from sera or tissues, WUH4 was isolated from a piglet stool sample. Here we announce its complete genome sequence.     

Tracing the spatio-temporal dynamics of endangered fin whales (<Emphasis Type="Italic">Balaenoptera physalus</Emphasis>) within baleen whale (Mysticeti) lineages: a mitogenomic perspective

To explore the spatio-temporal dynamics of endangered fin whales (Balaenoptera physalus) within the baleen whale (Mysticeti) lineages, we analyzed 148 published mitochondrial genome sequences of baleen whales. We used a Bayesian coalescent approach as well as Bayesian inferences and maximum likelihood methods. The results showed that the fin whales had a single maternal origin, and that there is a significant correlation between geographic location and evolution of global fin whales. The most recent common female ancestor of this species lived approximately 9.88 million years ago (Mya). Here, North Pacific fin whales first appeared about 7.48 Mya, followed by a subsequent divergence in Southern Hemisphere approximately 6.63 Mya and North Atlantic about 4.42 Mya. Relatively recently, approximately 1.76 and 1.42 Mya, there were two additional occurrences of North Pacific populations; one originated from the Southern Hemisphere and the other from an uncertain location. The evolutionary rate of this species was 1.002?×?10^?3 substitutions/site/My. Our Bayesian skyline plot illustrates that the fin whale population has the rapid expansion event since ~?2.5 Mya, during the Quaternary glaciation stage. Additionally, this study indicates that the fin whale has a sister group relationship with humpback whale (Meganoptera novaeangliae) within the baleen whale lineages. Of the 16 genomic regions, NADH5 showed the most powerful signal for baleen whale phylogenetics. Interestingly, fin whales have 16 species-specific amino acid residues in eight mitochondrial genes: NADH2, COX2, COX3, ATPase6, ATPase8, NADH4, NADH5, and Cytb.     

Bayesian Phylogeography of Crimean-Congo Hemorrhagic Fever Virus in Europe

Crimean-Congo hemorrhagic fever (CCHF) is a zoonosis mainly transmitted by ticks that causes severe hemorrhagic fever and has a mortality rate of 5-60%. The first outbreak of CCHF occurred in the Crimean peninsula in 1944-45 and it has recently emerged in the Balkans and eastern Mediterranean. In order to reconstruct the origin and pathway of the worldwide dispersion of the virus at global and regional (eastern European) level, we investigated the phylogeography of the infection by analysing 121 publicly available CCHFV S gene sequences including two recently characterised Albanian isolates. The spatial and temporal phylogeny was reconstructed using a Bayesian Markov chain Monte Carlo approach, which estimated a mean evolutionary rate of 2.96 x 10^-4 (95%HPD=1.6 and 4.7 x 10^-4) substitutions/site/year for the analysed fragment. All of the isolates segregated into seven highly significant clades that correspond to the known geographical clades: in particular the two new isolates from northern Albania clustered significantly within the Europe 1 clade. Our phylogeographical reconstruction suggests that the global CCHFV clades originated about one thousand years ago from a common ancestor probably located in Africa. The virus then spread to Asia in the XV century and entered Europe on at least two occasions: the first in the early 1800s, when a still circulating but less or non-pathogenic virus emerged in Greece and Turkey, and the second in the early 1900s, when a pathogenic CCHFV strain began to spread in eastern Europe. The most probable location for the origin of this European clade 1 was Russia, but Turkey played a central role in spreading the virus throughout Europe. Given the close proximity of the infected areas, our data suggest that the movement of wild and domestic ungulates from endemic areas was probably the main cause of the dissemination of the virus in eastern Europe.     

Variations in Abundance,Genome Size,Morphology, and Functional Role of the Virioplankton in Lakes Annecy and Bourget over a 1-Year Period

We sampled the surface waters (2–50 m) of two deep peri-alpine lakes over a 1-year period in order to examine (1) the abundance, vertical distribution, genome size, and morphology structures of the virioplankton; (2) the virus-mediated bacterial mortality; and (3) the specific genome size range of double-stranded DNA (dsDNA) phytoplankton viruses. Virus-like particle (VLP) concentrations varied between 4.16?×?10⁷ (January) and 2.08?×?10⁸ part mL^?1 (May) in Lake Bourget and between 2.7?×?10⁷ (June) and 8.39?×?10⁷ part mL^?1 (November) in Lake Annecy. Our flow cytometry analysis revealed at least three viral groups (referred to as virus-like particles 1, 2, and 3) that exhibited distinctive dynamics suggestive of different host types. Phage-induced bacterial mortality varied between 6.1 % (June) and 33.2 % (October) in Lake Bourget and between 7.4 % (June) and 52.6 % (November) in Lake Annecy, suggesting that viral lysis may be a key cause of mortality of the bacterioplankton. Virioplankton genome size ranged from 27 to 486 kb in Lake Bourget, while it reached 620 kb in Lake Annecy for which larger genome sizes were recorded. Our analysis of pulsed field gel electrophoresis bands using different PCR primers targeting both cyanophages and algal viruses showed that (1) dsDNA viruses infecting phytoplankton may range from 65 to 486 kb, and (2) both cyanophage and algal “diversity” were higher in Lake Annecy. Lakes Annecy and Bourget also differed regarding the proportions of both viral families (with the dominance of myoviruses vs. podoviruses) and infected bacterial morphotypes (short rods vs. elongated rods), in each of these lakes, respectively. Overall, our results reveal that (1) viruses displayed distinct temporal and vertical distribution, dynamics, community structure in terms of genome size and morphology, and viral activity in the two lakes; (2) the Myoviridae seemed to be the main cause of bacterial mortality in both lakes and this group seemed to be related to VLP2; and (3) phytoplankton viruses may have a broader range of genome size than previously thought. This study adds to growing evidence that viruses are diverse and play a significant role in freshwater microbial dynamics and more globally lake functioning. It highlights the importance of further considering this biological compartment for a better understanding of plankton ecology in peri-alpine lakes.     

Molecular evolution of the tick-borne encephalitis and Powassan viruses

Issues associated with newly emerging viruses, their genetic diversity, and viral evolution in modern environments are currently attracting growing attention. In this study, a phylogenetic analysis was performed and the evolution rate was evaluated for such pathogenic flaviviruses endemic to Russia as tick-borne encephalitis virus (TBEV) and Powassan virus (PV). The analysis involved 47 nucleotide sequences of the TBEV genome region encoding protein E and 17 sequences of the PV NS5-encoding region. The nucleotide substitution rate was estimated as 1.4 × 10⁻⁴ and 5.4 × 10⁻⁵ substitutions per site per year for the E protein-encoding region of the TBEV genome and for the NS5 genome region of PV, respectively. The ratio of non-synonymous to synonymous nucleotide substitutions (dN/dS) in viral sequences was calculated as 0.049 for TBEV and 0.098 for PV. The highest dN/dS values of 0.201–0.220 were found in the subcluster of Russian and Canadian PV strains, and the lowest value of 0.024 was observed in the cluster of Russian and Chinese strains of the Far Eastern TBEV genotype. Evaluation of time intervals between the events of viral evolution showed that the European subtype of TBEV diverged from the common TBEV ancestor approximately 2750 years ago, while the Siberian and Far Eastern subtypes emerged approximately 2250 years ago. The PV was introduced into its natural foci of the Russian Primorskii krai only approximately 70 years ago; these strains were very close to Canadian PV strains. The pattern of PV evolution in North America was similar to the evolution of the Siberian and Far Eastern TBEV subtypes in Asia. The moments of divergence between major genetic groups of TBEV and PV coincide with historical periods of climate warming and cooling, suggesting that climate change was a key factor in the evolution of flaviviruses in past millennia.     

Molecular evolution of poxviruses

Previous restriction fragment length polymorphism analysis divided variola virus (VARV) strains into two subtypes, one of which included West African and South American isolates. This allowed a dating to be introduced for the first time in estimation of the VARV evolution rate. The results were used to analyze the molecular evolution of the total family Poxviridae. Comparisons of the known nucleotide sequences were performed for the extended conserved central genome region in 42 orthopoxvirus strains and for the eight genes of multisubunit RNA polymerase in 65 viruses belonging to various genera of the family Poxviridae. Using the Bayesian dating method, the mutation accumulation rate of poxviruses was estimated at (1.7-8.8) x 10(-6) nucleotide substitutions per site per year. Computations showed that the modem poxvirus genera started diverging from an ancestral virus more than 200 thousand years ago and that an ancestor of the genus Orthopoxvirus emerged 131 +/- 45 thousand years ago. The other genera of mammalian poxviruses with a low GC content diverged approximately 110-90 thousand years ago. The independent evolution of VARV started 3.4 +/- 0.8 thousand years ago. It was shown with the example of VARV and the monkeypox virus (MPXV) that divergent evolution of these orthopoxviruses started and the West African subtypes of VARV and MPXV were formed as geographical conditions changed to allow isolation of West African animals from other African regions.     

The first phlebo‐like virus infecting plants: a case study on the adaptation of negative‐stranded RNA viruses to new hosts

A novel negative‐stranded (ns) RNA virus associated with a severe citrus disease reported more than 80 years ago has been identified. Transmission electron microscopy showed that this novel virus, tentatively named citrus concave gum‐associated virus, is flexuous and non‐enveloped. Notwithstanding, its two genomic RNAs share structural features with members of the genus Phlebovirus, which are enveloped arthropod‐transmitted viruses infecting mammals, and with a group of still unclassified phlebo‐like viruses mainly infecting arthropods. CCGaV genomic RNAs code for an RNA‐dependent RNA polymerase, a nucleocapsid protein and a putative movement protein showing structural and phylogenetic relationships with phlebo‐like viruses, phleboviruses and the unrelated ophioviruses, respectively, thus providing intriguing evidence of a modular genome evolution. Phylogenetic reconstructions identified an invertebrate‐restricted virus as the most likely ancestor of this virus, revealing that its adaptation to plants was independent from and possibly predated that of the other nsRNA plant viruses. These data are consistent with an evolutionary scenario in which trans‐kingdom adaptation occurred several times during the history of nsRNA viruses and followed different evolutionary pathways, in which genomic RNA segments were gained or lost. The need to create a new genus for this bipartite nsRNA virus and the impact of the rapid and specific detection methods developed here on citrus sanitation and certification are also discussed.     

Evolution of alphaviruses in the eastern equine encephalomyelitis complex.

Evolution of viruses in the eastern equine encephalomyelitis (EEE) complex was studied by analyzing RNA sequences and oligonucleotide fingerprints from isolates representing the North and South American antigenic varieties. By using homologous sequences of Venezuelan equine encephalomyelitis virus as an outgroup, phylogenetic trees revealed three main EEE virus monophyletic groups. A North American variety group included all isolates from North America and the Caribbean. One South American variety group included isolates from the Amazon basin in Brazil and Peru, while the other included strains from Argentina, Guyana, Ecuador, Panama, Trinidad, and Venezuela. No evidence of heterologous recombination was obtained when three separate regions of the EEE virus genome were analyzed independently. Estimates of the overall rate of EEE virus evolution (nucleotide substitution) were 1.6 x 10(-4) substitution per nucleotide per year for the North American group and 4.3 x 10(-4) for the Argentina-Panama South American group. Evolutionary rate estimates for the North American group increased over 10-fold (from about 2 x 10(-5) to 4 x 10(-4)) concurrent with divergence of two monophyletic groups during the early 1970s. The North and South American antigenic varieties diverged roughly 1,000 years ago, while the two main South American groups diverged about 450 years ago. Analysis of multiple strains isolated from an upstate New York transmission focus during the same years suggested that, in certain locations, EEE virus may be relatively isolated for short time periods.     

Diversification of rice yellow mottle virus and related viruses spans the history of agriculture from the neolithic to the present

The mechanisms of evolution of plant viruses are being unraveled, yet the timescale of their evolution remains an enigma. To address this critical issue, the divergence time of plant viruses at the intra- and inter-specific levels was assessed. The time of the most recent common ancestor (TMRCA) of Rice yellow mottle virus (RYMV; genus Sobemovirus) was calculated by a Bayesian coalescent analysis of the coat protein sequences of 253 isolates collected between 1966 and 2006 from all over Africa. It is inferred that RYMV diversified approximately 200 years ago in Africa, i.e., centuries after rice was domesticated or introduced, and decades before epidemics were reported. The divergence time of sobemoviruses and viruses of related genera was subsequently assessed using the age of RYMV under a relaxed molecular clock for calibration. The divergence time between sobemoviruses and related viruses was estimated to be approximately 9,000 years, that between sobemoviruses and poleroviruses approximately 5,000 years, and that among sobemoviruses approximately 3,000 years. The TMRCA of closely related pairs of sobemoviruses, poleroviruses, and luteoviruses was approximately 500 years, which is a measure of the time associated with plant virus speciation. It is concluded that the diversification of RYMV and related viruses has spanned the history of agriculture, from the Neolithic age to the present.     

猪圆环病毒2型及猪繁殖与呼吸综合症病毒的快速检测

According to the published genome sequences of Porcine circovirus type 2(PCV2)and Porcine reproductive and respiratory syndrome virus(PRRSV),primers were designed and PCR,RT-PCR were set up for the detection of PCV2 and PRRSV,respectively,With the established methods,38 clinical samples from the respiratory disease pigs were detected for the presence of PCV2 and PRRSV. The results demonstrated that 22 samples were positive for PCV2,27 samples were positive for PRRSV and among the above positive samples,18 samples were positive for both viruses,The data obtained in the present study indicated that PCV2 and PRRSV maybe play an important role in the course of the development of respiratiory diseases.     

Microbiological Identification and Analysis of Swine Lungs Collected from Carcasses in Swine Farms, China

The primary objective of this 3 years study was to determine the prevalence of porcine pathogens of the lungs of swine in swine farms in southern China. A total of 5,420 samples were collected from 200 swine farms. The bacterium that was most commonly isolated was Streptococcus suis, with 10.24 % of the samples being positive, 114 lungs (2.1 %) were positive for pseudorabies virus and 263 (4.85 %) were positive for classical swine fever virus; much lower than positive for PRRSV (15.1 %, p = 0.023) and PCV2 (13.8 %, p = 0.038). lungs that were positive for PRRSV and/or PCV-2 have significantly increased odds of being positive for any of the S. suis (9.79 vs. 0.44 %, p = 0.003).     

3D genomics imposes evolution of the domain model of eukaryotic genome organization

The hypothesis that the genome is composed of a patchwork of structural and functional domains (units) that may be either active or repressed was proposed almost 30 years ago. Here, we examine the evolution of the domain model of eukaryotic genome organization in view of the expansion of genome-scale techniques in the twenty-first century that have provided us with a wealth of information on genome organization, folding, and functioning.     

A serological survey of selected pathogens in wild boar (<Emphasis Type="Italic">Sus scrofa</Emphasis>) in northern Turkey

During the hunting season in March 2012, a total of 93 blood samples were collected from wild boars (Sus scrofa) shot in the area of northern Turkey (Samsun and Gumushane provinces). These blood samples were examined by enzyme immunoassay (ELISA) for the presence of antibodies to classical swine fever virus (CSFV), Aujeszky’s disease virus (ADV), porcine reproductive and respiratory syndrome virus (PRRSV), porcine respiratory coronavirus (PRCV), swine influenza virus (SIV), porcine parvovirus (PPV), swine vesicular disease virus (SVDV), hepatitis E virus (HEV), African swine fever virus (ASFV), porcine rotavirus (PRV), transmissible gastroenteritis virus (TGEV) and bovine viral diarrhoea virus (BVDV). Out of 93 serum samples examined, 65 (69.9 %) were positive for PRV, 22 (23.7 %) were positive for ADV, 5 (5.4 %) were positive for BVDV, 4 (4.3 %) were positive for PPV and 2 (2.2 %) were positive for PRRSV. All sera were negative for ASFV, SVDV, HEV, SIV, PRCV, TGEV and CSFV. The results, recorded for the first time in Turkey, supported the hypothesis that wild boar act as a potential reservoir of selected viruses and thus have a role in the epidemiology of these diseases.