Comments on duck circovirus (DuCV) genotype definition

A standardised methodology has been used to define genotypes based on pairwise sequence comparisons (PASC). PASC is a widely accepted method in virus taxonomy, which is based on the histogram of pairwised differences among sequences. Recently, Zhang et al. (2013) concluded that the average p-distance of duck circovirus (DuCV) between genotypes 1 and 2 was 0.170, and subtype distance thresholds were 0.032 in DuCV-1 and 0.018 in DuCV-2, respectively. However, there might be some concerns on the methodology application to define the genotype of DuCV. Taking into account the concerns mentioned above, our authors conducted the PASC analyses of 54 capsid gene (ORF2) and genomic sequences including all the sequences from Zhang et al. (2013). Our results confirmed the existence of two DuCV genotypes (1 and 2) and, we suggest that DuCV ORF2 and genomic distance genotype thresholds were 0.061 and 0.038, respectively.     

The marsupial MHC: The Tammar Wallaby,Macropus eugenii,contains an expressed DNA-like gene on chromosome 1

Isolates of cauliflower mosaic virus (CaMV) differ in host range and symptomatology. Knowledge of their sequence relationships should assist in identifying nucleotide sequences responsible for isolate-specific characters. Complete nucleotide sequences of the DNAs of eight isolates of CaMV were aligned and the aligned sequences were used to analyze phylogenetic relationships by maximum likelihood, bootstrapped parsimony, and distance methods. Isolates found in North America clustered separately from those isolated from other parts of the world. Additional isolates, for which partial sequences were available, were incorporated into phylogenetic analysis of the sequences of genome segments corresponding to individual protein coding regions or the large intergenic region of CaMV DNA. The analysis revealed several instances where the position of an isolate on a tree for one coding region did not agree with the position of the isolate on the tree for the complete genome or with its position on trees for other coding regions. Examination of the distribution of shared residue types of phylogenetically informative positions in anomalous regions suggested that most of the anomalies were due to recombination events during the evolution of the isolates. Application of an algorithm that searches for segments of significant length that are identical between pairs of isolates or contain a significantly high concentration of polymorphisms suggested two additional recombination events between progenitors of the isolates studied and an event between the XinJing isolate and a CaMV not represented in the data set. An earlier phylogenetic origin for CaMV than for carnation etched ring virus, the caulimovirus used as outgroup in these analyses, was deduced from the position of the outgroup with North American isolates in some trees, but with non-North American isolates in other trees. Correspondence to: U. Melcher     

全基因组序列的系统发育与重组分析鉴定我国马铃薯Y病毒株系

为明确我国马铃薯Y病毒(Potato virus Y,PVY)株系组成及其在不同寄主上的分布特征,文章对31条PVY中国分离物的基因组编码区进行了系统发育及重组分析。系统发育分析结果显示,最大似然法(maximum likelihood,ML)和邻接法(neighbor-joining,NJ)重建的系统发育拓扑结构基本一致,GF_YL20、FZ10、DF、SD1、1116和9703-3等6个分离物均未与已知株系相聚,其他株系来源相同的分离物以较高的自展值相聚成簇,显示出较强的株系特异性。重组分析显示,31个PVY分离物中,共有28个PVY分离物均检测到具有显著的重组信号,这些分离物主要包括PVY~(N-Wi)、PVY~(NTN)和PVY~(NTN-NW)三种株系。进一步分析显示,分离物1104、DF、SD1、1116和9703-3为未报道过的PVY新重组分离物。株系统计结果显示PVY重组株系已上升为优势株系,生产上需要密切注意这些株系,尤其新重组分离物的发展动态。以上结果表明,综合应用PVY全基因组序列的系统发育和重组分析,可以实现PVY株系的准确鉴定。     

PCR detection and sequence analysis of duck circovirus in sick Muscovy ducks

The duck circovirus (DuCV) infection in sick ducks from Fujian Province was investigated. The liver samples of 43 sick Muscovy ducks with infectious serositis were collected from 12 duck farms in Fujian Province.Based on the published sequences of DuCV, two primers were designed for the detection of DuCV and four pairs of primers were designed to amplify four overlapping fragments that cover the complete genome of DuCV. The specific PCR products were amplified from positive samples. The fragments were then cloned into pMD18-T vector and sequenced, and the full length genomic sequence of the FJ0601 isolate of DuCV was obtained. PCR analysis showed that the proportion of ducks which were positive for circovirus was 79% and 10 out of the 12 farms were positive. Sequence analysis showed that the complete genome of DuCV-FJ0601 was 1988 bp and possessed features common to the family Circoviridae which included a stem-loop structure and the Rep protein motifs. Homology analysis showed that FJ0601 isolate of DuCV had 97.3%～97.5% nucleotide sequence identity to all the four Taiwan isolates (TC1/2002, TC2/2002, TC3/2002, TC4/2002), 82.9% identity to the America (33753-52) isolate and 82.3% identity to the Germany isolate. Phylogenetic analysis with Clustal W, however,showed that FJ0601 isolate of DuCV was on a common branch with Taiwan isolates, and Germany and America isolates belonged to the other branch.     

Complete-genome phylogenetic approach to varicella-zoster virus evolution: genetic divergence and evidence for recombination

Recent studies of varicella-zoster virus (VZV) DNA sequence variation, involving large numbers of globally distributed clinical isolates, suggest that this virus has diverged into at least three distinct genotypes designated European (E), Japanese (J), and mosaic (M). In the present study, we determined and analyzed the complete genomic sequences of two M VZV strains and compared them to the sequences of three E strains and two J strains retrieved from GenBank (including the Oka vaccine preparation, V-Oka). Except for a few polymorphic tandem repeat regions, the whole genome, representing approximately 125,000 nucleotides, is highly conserved, presenting a genetic similarity between the E and J genotypes of approximately 99.85%. These analyses revealed that VZV strains distinctly segregate into at least four genotypes (E, J, M1, and M2) in phylogenetic trees supported by high bootstrap values. Separate analyses of informative sites revealed that the tree topology was dependent on the region of the VZV genome used to determine the phylogeny; collectively, these results indicate the observed strain variation is likely to have resulted, at least in part, from interstrain recombination. Recombination analyses suggest that strains belonging to the M1 and M2 genotypes are mosaic recombinant strains that originated from ancestral isolates belonging to the E and J genotypes through recombination on multiple occasions. Furthermore, evidence of more recent recombination events between M1 and M2 strains is present in six segments of the VZV genome. As such, interstrain recombination in dually infected cells seems to figure prominently in the evolutionary history of VZV, a feature it has in common with other herpesviruses. In addition, we report here six novel genomic targets located in open reading frames 51 to 58 suitable for genotyping of clinical VZV isolates.     

Recombination between sequences of hepatitis B virus from different genotypes

A comparison of 25 hepatitis B virus (HBV) isolates for which complete genome sequences are available revealed two that occupied different positions in phylogenetic trees reconstructed from different open reading frames. Further analysis indicated that this incongruence was the result of recombination between viruses of different genomic and antigenic types. Both putative recombinants originated from geographic regions where multiple genotypes are known to cocirculate. A search of the sequence databases showed evidence of similar intergenotypic recombinants. These observations indicate that recombination between divergent strains may represent an important source of genetic variation in HBV.     

Construction of a virtual Mycobacterium tuberculosis consensus genome and its application to data from a next generation sequencer

Background

Although Mycobacterium tuberculosis isolates are consisted of several different lineages and the epidemiology analyses are usually assessed relative to a particular reference genome, M. tuberculosis H37Rv, which might introduce some biased results. Those analyses are essentially based genome sequence information of M. tuberculosis and could be performed in sillico in theory, with whole genome sequence (WGS) data available in the databases and obtained by next generation sequencers (NGSs). As an approach to establish higher resolution methods for such analyses, whole genome sequences of the M. tuberculosis complexes (MTBCs) strains available on databases were aligned to construct virtual reference genome sequences called the consensus sequence (CS), and evaluated its feasibility in in sillico epidemiological analyses.

Results

The consensus sequence (CS) was successfully constructed and utilized to perform phylogenetic analysis, evaluation of read mapping efficacy, which is crucial for detecting single nucleotide polymorphisms (SNPs), and various MTBC typing methods virtually including spoligotyping, VNTR, Long sequence polymorphism and Beijing typing. SNPs detected based on CS, in comparison with H37Rv, were utilized in concatemer-based phylogenetic analysis to determine their reliability relative to a phylogenetic tree based on whole genome alignment as the gold standard. Statistical comparison of phylogenic trees based on CS with that of H37Rv indicated the former showed always better results that that of later. SNP detection and concatenation with CS was advantageous because the frequency of crucial SNPs distinguishing among strain lineages was higher than those of H37Rv. The number of SNPs detected was lower with the consensus than with the H37Rv sequence, resulting in a significant reduction in computational time. Performance of each virtual typing was satisfactory and accorded with those published when those are available.

Conclusions

These results indicated that virtual CS constructed from genome sequence data is an ideal approach as a reference for MTBC studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1368-9) contains supplementary material, which is available to authorized users.     

Isolation and identification of swine influenza recombinant A/Swine/Shandong/1/2003(H9N2) virus

Ten influenza virus isolates were obtained from infected pigs from different places in Shandong province showing clinical symptoms from October 2002 to January 2003. All 10 isolates were identified in China's National Influenza Research Center as influenza A virus of H9N2 subtype. The complete genome of one isolate, designated A/Swine/Shandong/1/2003(H9N2), was sequenced and compared with sequences available in GenBank. The results of analyses indicated that the sequence of A/Swine/Shandong/1/2003(H9N2) was similar to those of several chicken influenza viruses and duck influenza viruses recently prevalent in South China. According to phylogenetic analysis of the complete gene sequences, A/Swine/Shandong/1/2003(H9N2) possibly originated from the reassortment of chicken influenza viruses and duck influenza viruses. It was found that the amino acid sequence at the HA cleavage site in Sw/SD/1/2003 is R-S-L-R-G, differing clearly from that of other H9N2 subtype isolates of swine influenza and avian influenza, which is R-S-S-R-G.     

Complete nucleotide sequence analysis of <Emphasis Type="Italic">Cymbidium mosaic virus</Emphasis> Indian isolate: further evidence for natural recombination among potexviruses

The complete nucleotide sequence of an Indian strain of Cymbidium mosaic virus (CymMV) was determined and compared with other potexviruses. Phylogenetic analyses on the basis of RNA-dependent RNA polymerase (RdRp), triple gene block protein and coat protein (CP) amino acid sequences revealed that CymMV is closely related to the Narcissus mosaic virus (NMV), Scallion virus X (SVX), Pepino mosaic virus (PepMV) and Potato aucuba mosaic virus (PAMV). Different sets of primers were used for the amplification of different regions of the genome through RT-PCR and the amplified genes were cloned in a suitable vector. The full genome of the Indian isolate of CymMV from Phaius tankervilliae shares 96–97% similarity with isolates reported from other countries. It was found that the CP gene of CymMV shares a high similarity with each other and other potexviruses. One of the Indian isolates seems to be a recombinant formed by the intermolecular recombination of two other CymMV isolates. The phylogenetic analyses, Recombination Detection Program (RDP2) analyses and sequence alignment survey provided evidence for the occurrence of a recombination between an Indian isolate (AM055720) as the major parent, and a Korean type-2 isolate (AF016914) as the minor parent. Recombination was also observed between a Singapore isolate (U62963) as the major parent, and a Taiwan CymMV (AY571289) as the minor parent.     

Complete Genome Sequencing of a Velogenic Viscerotropic Avian Paramyxovirus 1 Isolated from Pheasants (Pucrasia macrolopha) in Lahore,Pakistan

We report the complete genome sequence of avian paramyxovirus 1 (APMV-1) isolated from an acute and highly contagious outbreak in pheasants (Pucrasia macrolopha) in Lahore, Pakistan. Biological and serological characterization showed a velogenic strain of APMV-1, which was further confirmed by the sequence analysis of the cleavage site in the fusion protein. Complete genome sequencing and phylogenetic analysis indicated that this isolate belonged to genotype VII, specifically to subgenotype VIIa, and clustered closely with isolates characterized from Indonesia. Notably, the isolate showed significant differences from previously characterized APMV-1 from Pakistani commercial and rural chicken.     

Mosaic genome structure of simian immunodeficiency virus from west African green monkeys.

Elucidation of the phylogenetic origins of simian and human immunodeficiency viruses (SIV and HIV) is fundamental to the understanding of HIV pathogenesis and the spread of AIDS worldwide. In this study, we molecularly characterized multiple SIVAGM isolates from four different African green monkey species (vervet, grivet, sabaeus and tantalus monkeys). Phylogenetic analysis of partial (1 kb) env sequences indicated that all SIVAGM strains cluster together, and that they fall into four distinct sequence sub-groups according to their species of origin. However, alignment of long terminal repeat sequences revealed that SIVs from West African sabaeus monkeys contain a structural feature (a duplication of the transactivation response element) thus far only found in otherwise highly divergent lentiviruses infecting sooty mangabeys (SIVSM) and humans (HIV-2). To determine whether there were additional similarities with the SIVSM/HIV-2 group, a full-length replication competent sabaeus provirus was cloned and sequenced. In phylogenetic trees derived from the central and 3' coding regions, the sabaeus virus clustered with SIVAGM isolates from other African green monkey species. However, in trees derived from the 3' half of gag and the adjacent 5' region of pol, the sabaeus virus grouped with the SIVSM/HIV-2 lineage. These results indicated that the sabaeus virus comprised a mosaic genome which must have resulted from recombination of divergent lentiviruses in the distant past. A second, independent sabaeus isolate exhibited similar phylogenetic relationships, suggesting that all West African green monkey viruses share this complex evolutionary history. Taken together, these results indicate that African green monkeys have been infected with SIVAGM for very long periods of time, and that recombination and cross-species transmission in the wild have contributed to the genetic complexity of primate lentiviruses.     

Molecular characterization of Indian Sugarcane streak mosaic virus isolates reveals recombination and negative selection in the P1 gene

Sugarcane streak mosaic virus (SCSMV), a member of the genus Poacevirus is an important viral pathogen affecting sugarcane production in India. The P1 gene of ten Indian isolates was sequenced and compared with previously reported SCSMV isolates. Comparative sequence analysis revealed a high level of diversity in the P1 gene (83–98% nucleotide sequence identity; 87–100% amino acid sequence identity), and the Indian SCSMV isolates were found to be the most variable (up to 9% diversity at the amino acid level). Phylogenetic tree analysis showed clustering of 17 SCSMV isolates into two groups: group I included isolates from India (except SCSMV-TPT) and Pakistan, and group II consisted of isolates from Japan, Indonesia, Thailand and SCSMV-TPT. The results obtained from phylogenetic study were further supported by the different in silico analysis viz. SNPs (single nucleotide polymorphism), INDELs (insertion and deletion) and evolutionary distance analysis. A significant proportion of recombination sites were observed at the N terminal region of P1 gene. Analysis of selection pressure indicated that the P1 gene of the Indian SCSMV isolates is under strong negative or purifying selection. It is likely that recombination identified in Indian SCSMV isolates, along with strong purifying selection, enhances the speed of elimination of deleterious mutations in the P1 gene. The evolutionary processes (recombination and selection pressure) together contributed to the observed genetic diversity and population structure of Indian SCSMV isolates.     

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.     

Site-Specific Integration in Mammalian Cells Mediated by a New Hybrid Baculovirus–Adeno-Associated Virus Vector

Baculovirus can transiently transduce primary human and rat hepatocytes, as well as a subset of stable cell lines. To prolong transgene expression, we have developed new hybrid vectors which associate key elements from adeno-associated virus (AAV) with the elevated transducing capacity of baculovirus. The hybrid vectors contain a transgene cassette composed of the β-galactosidase (β-Gal) reporter gene and the hygromycin resistance (Hyg^r) gene flanked by the AAV inverted terminal repeats (ITRs), which are necessary for AAV replication and integration in the host genome. Constructs were derived both with and without the AAV rep gene under the p5 and p19 promoters cloned in different positions with respect to the baculovirus polyheidrin promoter. A high-titer preparation of baculovirus-AAV (Bac-AAV) chimeric virus containing the ITR–Hyg^r–β-Gal sequence was obtained with insect cells only when the rep gene was placed in an antisense orientation to the polyheidrin promoter. Infection of 293 cells with Bac-AAV virus expressing the rep gene results in a 10- to 50-fold increase in the number of Hyg^r stable cell clones. Additionally, rep expression determined the localization of the transgene cassette in the aavs1 site in approximately 41% of cases as detected by both Southern blotting and fluorescent in situ hybridization analysis. Moreover, site-specific integration of the ITR-flanked DNA was also detected by PCR amplification of the ITR-aavs1 junction in transduced human fibroblasts. These data indicate that Bac-AAV hybrid vectors can allow permanent, nontoxic gene delivery of DNA constructs for ex vivo treatment of primary human cells.     

Recombination in the Evolution of Enterovirus C Species Sub-Group that Contains Types CVA-21, CVA-24, EV-C95, EV-C96 and EV-C99

Genetic recombination is considered to be a very frequent phenomenon among enteroviruses (Family Picornaviridae, Genus Enterovirus). However, the recombination patterns may differ between enterovirus species and between types within species. Enterovirus C (EV-C) species contains 21 types. In the capsid coding P1 region, the types of EV-C species cluster further into three sub-groups (designated here as A–C). In this study, the recombination pattern of EV-C species sub-group B that contains types CVA-21, CVA-24, EV-C95, EV-C96 and EV-C99 was determined using partial 5′UTR and VP1 sequences of enterovirus strains isolated during poliovirus surveillance and previously published complete genome sequences. Several inter-typic recombination events were detected. Furthermore, the analyses suggested that inter-typic recombination events have occurred mainly within the distinct sub-groups of EV-C species. Only sporadic recombination events between EV-C species sub-group B and other EV-C sub-groups were detected. In addition, strict recombination barriers were inferred for CVA-21 genotype C and CVA-24 variant strains. These results suggest that the frequency of inter-typic recombinations, even within species, may depend on the phylogenetic position of the given viruses.     

Complete Genomic Sequence of a New Muscovy Duck-Origin Reovirus from China

The complete genomic sequence of a new Muscovy duck-origin reovirus (N-MDRV), strain J18 from China, was determined. The virus has a tricistronic S1 genome segment that is distinct from the originally described MDRV, which possesses a bicistronic S4 genome segment. Pairwise comparisons and phylogenetic analyses suggest that N-MDRV J18 is a new isolate within the species Avian orthoreovirus.     

Frequency and dynamics of recombination within different species of human enteroviruses

Enteroviruses are members of the family Picornaviridae that cause widespread infections in human and other mammalian populations. Enteroviruses are genetically and antigenically highly variable, and recombination within and between serotypes contributes to their genetic diversity. To investigate the dynamics of the recombination process, sequence phylogenies between three regions of the genome (VP4, VP1, and 3Dpol) were compared among species A and B enterovirus variants detected in a human population-based survey in Scotland between 2000 and 2001, along with contemporary virus isolates collected in the same geographical region. This analysis used novel bioinformatic methods to quantify phylogenetic compatibility and correlations with serotype assignments of evolutionary trees constructed for different regions of the enterovirus genome. Species B enteroviruses showed much more frequent, time-correlated recombination events than those found for species A, despite the equivalence in population sampling, concordant with a linkage analysis of previously characterized enterovirus sequences obtained over longer collection periods. An analysis of recombination among complete genome sequences by computation of a phylogenetic compatibility matrix (PCM) demonstrated sharply defined boundaries between the VP2/VP3/VP1 block and sequences to either side in phylogenetic compatibility. The PCM also revealed equivalent or frequently greater degrees of incompatibility between different parts within the nonstructural region (2A-3D), indicating the occurrence of extensive recombination events in the past evolution of this part of the genome. Together, these findings provide new insights into the dynamics of species A and B enterovirus recombination and evolution and into the contribution of structured sampling to documenting reservoirs, emergence, and spread of novel recombinant forms in human populations.     

Genetic Structure and Molecular Variability of Cucumber mosaic virus Isolates in the United States

Cucumber mosaic virus (CMV) has a worldwide distribution and the widest host range of any known plant virus. From 2000 to 2012, epidemics of CMV severely affected the production of snap bean (Phaseulos vulgaris L.) in the Midwest and Northeastern United States. Virus diversity leading to emergence of new strains is often considered a significant factor in virus epidemics. In addition to epidemics, new disease phenotypes arising from genetic exchanges or mutation can compromise effectiveness of plant disease management strategies. Here, we captured a snapshot of genetic variation of 32 CMV isolates collected from different regions of the U.S including new field as well as historic isolates. Nucleotide diversity (π) was low for U.S. CMV isolates. Sequence and phylogenetic analyses revealed that CMV subgroup I is predominant in the US and further showed that the CMV population is a mixture of subgroups IA and IB. Furthermore, phylogenetic analysis suggests likely reassortment between subgroups IA and IB within five CMV isolates. Based on phylogenetic and computational analysis, recombination between subgroups I and II as well as IA and IB in RNA 3 was detected. This is the first report of recombination between CMV subgroups I and II. Neutrality tests illustrated that negative selection was the major force operating upon the CMV genome, although some positively selected sites were detected for all encoded proteins. Together, these data suggest that different regions of the CMV genome are under different evolutionary constraints. These results also delineate composition of the CMV population in the US, and further suggest that recombination and reassortment among strain subgroups does occur but at a low frequency, and point towards CMV genomic regions that differ in types of selection pressure.     

Characteristics of Nucleotide Substitution in the Hepatitis C Virus Genome: Constraints on Sequence Change in Coding Regions at Both Ends of the Genome

Comparison of complete genome sequences for different variants of hepatitis C virus (HCV) reveals several different constraints on sequence change. Synonymous changes are suppressed in coding regions at both 5′ and 3′ ends of the genome. No evidence was found for the existence of alternative reading frames or for a lower mutation frequency in these regions. Instead, suppression may be due to constraints imposed by RNA secondary structures identified within the core and NS5b genes. Nonsynonymous substitutions are less frequent than synonymous ones except in the hypervariable region of E2 and, to a lesser extent, in E1, NS2, and NS5b. Transitions are more frequent than transversions, particularly at the third position of codons where the bias is 16:1. In addition, nucleotide substitutions may not occur symmetrically since there is a bias toward G or C at the third position of codons, while T ↔ C transitions were twice as frequent as A ↔ G transitions. These different biases do not affect the phylogenetic analysis of HCV variants but need to be taken into account in interpreting sequence change in longitudinal studies. Received: 9 September 1996 / Accepted: 20 April 1997