Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • mSphere of Influence: Commentaries from Early Career Microbiologists
    • Archive
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mSphere
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSphere
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • mSphere of Influence: Commentaries from Early Career Microbiologists
    • Archive
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mSphere
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Clinical Science and Epidemiology

Phylogenetic Analyses Suggest that Factors Other Than the Capsid Protein Play a Role in the Epidemic Potential of GII.2 Norovirus

Kentaro Tohma, Cara J. Lepore, Lauren A. Ford-Siltz, Gabriel I. Parra
W. Paul Duprex, Editor
Kentaro Tohma
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Cara J. Lepore
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Lauren A. Ford-Siltz
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gabriel I. Parra
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Gabriel I. Parra
W. Paul Duprex
Boston University School of Medicine
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Martin Chan
Chinese University of Hong Kong
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rowena Bull
University of New South Wales
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jan Vinje
Centers for Disease Control and Prevention
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSphereDirect.00187-17
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

Article Figures & Data

Figures

  • Tables
  • Supplemental Material
  • FIG 1
    • Open in new tab
    • Download powerpoint
    FIG 1

    Evolutionary dynamics of the major capsid protein of GII.2 noroviruses. (a) Maximum likelihood phylogenetic tree of VP1-encoding (nt) sequences from GII.2 strains. The numbers on the ancestral nodes of the clusters represent the node support value calculated by the approximate likelihood-ratio test using PhyML. (b) Root-to-tip divergence plot of VP1 nt sequences of GII.2 strains. The x axis indicates isolation year, and the y axis shows the root-to-tip divergence on the maximum likelihood phylogenetic tree. The black line indicates a linear regression line of the root-to-tip divergence and isolation year. A total of 151 sequences were included for both analyses. Each strain is represented by a circle and colored according to its polymerase genotype. Strains whose polymerase genotype was not available are indicated in black. The genomic region used in the analyses spanned nt 5085 to 6419 relative to Snow Mountain virus (GenBank accession number AY134748). (c) The kernel density plot indicates the posterior estimates of substitution rate (substitutions/site/year) of VP1 nt sequences. To calculate the substitution rate of the VP1 associated with different polymerases, strains were clustered using the VP1 phylogenetic tree. Thus, the P2 cluster is shaded in red, the P16 cluster in green, and the P2-P16 mix in yellow. The overall rate is indicated by black shading.

  • FIG 2
    • Open in new tab
    • Download powerpoint
    FIG 2

    Amino acid diversification of the major capsid protein of GII.2 noroviruses. (a) Heat map indicating the accumulation of aa mutations in the partial VP1 sequences of GII.2 strains (445 amino acids) over time. The x axis indicates the pairwise aa mutations, and the y axis indicates isolation year. The number of pairwise comparisons is represented by the color gradient in the heat map. (b) The plot indicates the substitution rates (substitutions/site/year) of VP1 nucleotide sequences at each codon position. The mean and the 95% highest posterior density intervals are indicated. The substitution rates were calculated using the same sequence matrix and groups as those represented in Fig. 1.

  • FIG 3
    • Open in new tab
    • Download powerpoint
    FIG 3

    Positive evolution of the major capsid protein of GII.2 noroviruses. Maximum likelihood trees of VP1 nucleotide (nt) sequences of GII.2 indicate the branches under possible positive selection. The branches with evidence of positive selection (empirical Bayes factor of >10, P < 0.05; estimated using the mixed-effect model of evolution [MEME] method) at codon positions 344 (a), 345 (b), 384 (c), and 385 (d) are represented by red branch lines. The color shade indicates the phylogenetic clustering based on the polymerase genotype as indicated in Fig. 1. Reemerging GII.P16-GII.2 strains are indicated in blue.

  • FIG 4
    • Open in new tab
    • Download powerpoint
    FIG 4

    Evolutionary dynamics of the GII.2 polymerase. (a) Maximum likelihood trees of RdRp-encoding (nt) sequences of strains GII.P2 (n = 72) and GII.P16 (n = 131). Phylogenetic clusters based on VP1 nt sequences (Fig. 1) are shown here; the P2 cluster is represented in red, the P16 cluster in green, and the P2-P16 mix cluster in yellow. Reemerging GII.P16 strains, including GII.P16-GII.2 and GII.P16-GII.4, are indicated in blue. Each strain is represented by a circle, with different colors according to the associated capsid genotypes; GII.P2-GII.2 and GII.P16-GII.2 strains are indicated by red and green, respectively. The numbers on the ancestral nodes of the clusters represent the node support value calculated by the approximate likelihood-ratio test using PhyML. (b) The root-to-tip divergence plot of RdRp-encoding (nt) sequences of GII.P2 and GII.P16 strains. The x axis indicates isolation year, and the y axis shows the root-to-tip divergence on the maximum likelihood phylogenetic tree. The black lines indicate a linear regression line of the root-to-tip divergence and isolation year. Each strain is colored according to its capsid genotype. The genomic regions used spanned nt 4385 to 5104 relative to Snow Mountain virus (GenBank accession number AY134748). (c) The kernel density plot indicates the posterior estimates of substitution rate (substitutions/site/year) among GII.P2, GII.P16, and GII.P16 without reemerging strains. (d) Structural mapping of the amino acid (aa) substitutions in the reemerging GII.P16-GII.2 strains. Conservative mutations in the RdRp aa sequence from the reemerging GII.P16 strains compared to pre-2016 GII.P16 strains (P16 cluster and P2-P16 mix cluster) are mapped on the structural model of a GII.P4 RdRp (PDB number 4QPX). The molecular model was visualized using Chimera v.1.11. aa substitutions are indicated in blue. Residue 291, a residue that was shown by Bull et al. (46) to alter the incorporation rate of GII.4 polymerases, is indicated in red.

Tables

  • Figures
  • Supplemental Material
  • TABLE 1

    Amino acid substitution observed in the consensus VP1 sequence among phylogenetic clusters

    Strain or clusterNo. of sequencesAmino acid substitution at codon:
    71130303335344354364386400
    Reemerging GII.P16-GII.2 (2016)27AVVI (1 V)aS (1G)aGANE
    P2-P16 mix clusterb (2005–2015) (without reemerging strains)9AVVI (2 V)aSGANE
    P16 clusterb (2008–2014)38SVVVT (1 I, 2 A)aAAND (2 E)a
    P2 clusterb (2007–2010)29AI (4 V)aI (1 V)aISGASE
    P2 clusterb (2004–2006)18AVVVSGS (1P)aNE
    • ↵a Data in parentheses represent minor variants within the cluster.

    • ↵b Phylogenetic clustering as shown in Fig. 1, assigned based on the phylogenetic clustering of the VP1-encoding sequence and the associated polymerase genotype.

  • TABLE 2

    Codon position of VP1 sequences (1,337 nucleotides) with positive selection

    Positively selected siteMethod(s) (P value)Domain
    24MEMEa (0.015)N terminus
    78iFELb (0.015)Shell
    275MEME (0.016)P1
    344MEME (0.039)P2 (surface)
    345iFEL (0.01), MEME (0.020)P2 (surface)
    384MEME (0.01)P2 (surface)
    385MEME (0.024)P2 (surface)
    397iFEL (0.045)P2 (surface)
    • ↵a MEME, mixed-effect model of evolution.

    • ↵b iFEL, internal fixed-effect likelihood methods.

Supplemental Material

  • Figures
  • Tables
  • FIG S1

    Histogram of the pairwise p-distance of VP1 nucleotide (nt) sequences among GII.2 strains (n = 151). The pairwise distance data for three strains from the GII.P21 cluster (including two strains with unknown polymerase genotypes [KC998960 and AY660568] and JQ320072) and three from the GII.P22 cluster (DQ366347, AB279554, and AB279556) are indicated in black. The genomic region used spanned nt 5085 to 6713 of Snow Mountain virus (GenBank accession number AY134748). Download FIG S1, TIF file, 1 MB.

    Copyright © 2017 Tohma et al.

    This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

  • FIG S2

    Evolutionary dynamics of the major capsid protein of GII.2 noroviruses. Evolutionary dynamics of VP1 was analyzed by using full-length VP1-encoding (nt) sequences. (a) Maximum likelihood phylogenetic tree of VP1 nt sequences from GII.2 strains. The numbers on the ancestral nodes of the clusters represent the node support values calculated by the approximate likelihood-ratio test using PhyML. (b) Root-to-tip divergence plot of VP1 nt sequences of GII.2 strains. The x axis indicates isolation year, and the y axis shows the root-to-tip divergence on the maximum likelihood phylogenetic tree. The black line indicates a linear regression line of the root-to-tip divergence and isolation year. A total of 134 sequences are included for both analyses. Each strain is colored according to its polymerase genotype. Strains whose polymerase genotype was not available are indicated in black. The genomic region used in the analyses spanned nt 5085 to 6713 relative to Snow Mountain virus (GenBank accession number AY134748). (c) The kernel density plot indicates the posterior estimates of substitution rates (substitutions/site/year) of VP1 nt sequences. To calculate the substitution rate of the VP1 protein associated with different polymerases, strains were clustered using the VP1 phylogenetic tree. Thus, the P2 cluster is shaded in red, the P16 cluster in green, and the P2-P16 mix in yellow. The overall rate is indicated by black shading. Download FIG S2, TIF file, 1.9 MB.

    Copyright © 2017 Tohma et al.

    This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

  • TABLE S1

    Data set. Download TABLE S1, PDF file, 0.04 MB.

    Copyright © 2017 Tohma et al.

    This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

  • FIG S3

    Amino acid diversification of the major capsid protein of GII.2 noroviruses. The heat map indicates the accumulation of aa mutations in the complete VP1 of GII.2 strains (543 amino acids) over time. The x axis indicates the pairwise aa mutations, and the y axis indicates isolation year. The number of pairwise comparisons is represented by the color gradient in the heat map. Download FIG S3, TIF file, 0.3 MB.

    Copyright © 2017 Tohma et al.

    This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

  • FIG S4

    Maximum likelihood phylogenetic tree of VP1 amino acid sequences of GII.2 strains. Each strain is colored according to its polymerase genotype. Strains whose polymerase genotype was not available are indicated in black. A total of 151 sequences were included in the analysis. The genomic region used spanned nt 5085 to 6419 relative to Snow Mountain virus (GenBank accession number AY134748). The numbers on the ancestral nodes of the clusters represent the node support value calculated by the approximate likelihood-ratio test using PhyML. The color shade indicates the phylogenetic clustering based on the polymerase genotype as indicated in Fig. 1. Download FIG S4, TIF file, 1.1 MB.

    Copyright © 2017 Tohma et al.

    This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

  • TABLE S2

    Codon position of complete VP1 sequences (1,629 nt) with positive selection. Download TABLE S2, PDF file, 0.03 MB.

    Copyright © 2017 Tohma et al.

    This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

  • FIG S5

    Maximum likelihood trees of the RdRp amino acid sequences. Trees were built using GII.P2 (n = 72) (a) and GII.P16 (n = 131) (b) strains. Phylogenetic clusters based on VP1 nt sequences (Fig. 1) are shown here; the P2 cluster is represented in red, the P16 cluster in green, and the P2-P16 mix cluster in yellow. Reemerging GII.P16 strains, including GII.P16-GII.2 and GII.P16-GII.4, are indicated in blue. Each strain is represented by a circle, with different colors according to the associated capsid genotypes; GII.P2-GII.2 and GII.P16-GII.2 strains are indicated by red and green, respectively. The numbers on the ancestral nodes of the clusters represent the node support value calculated by the approximate likelihood-ratio test using PhyML. Download FIG S5, TIF file, 1.5 MB.

    Copyright © 2017 Tohma et al.

    This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

  • TABLE S3

    Amino acid substitution observed in the RdRp sequence of reemerging GII.P16 strains. Download TABLE S3, PDF file, 0.03 MB.

    Copyright © 2017 Tohma et al.

    This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

  • FIG S6

    Maximum clade credibility (MCC) tree of GII.2 strains. VP1-encoding sequences of GII.2 strains (n = 151) (a), RdRp-encoding sequences from GII.P2 strains (n = 72) (b), and GII.P16 strains (n = 131) (c) were used for tree reconstructions. Each strain in the VP1 MCC tree is colored according to its polymerase genotype (a), and each strain in the RdRp MCC trees is colored by its capsid genotype (b and c). Strains whose polymerase genotype was not available are indicated in black. Reemerging GII.P16-GII.2 strains are indicated by blue shading. Arrows indicate the means and 95% highest posterior density intervals of the time of the most recent common ancestor. Phylogenetic clustering is indicated by color as described for Fig. 1. The numbers on the ancestral nodes of the clusters represent the posterior probability values of the nodes. The x axis and branch length represent the year. Estimated divergence time is indicated on the node of common ancestors of phylogenetic clustering. Download FIG S6, TIF file, 2.3 MB.

    Copyright © 2017 Tohma et al.

    This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

  • TABLE S4

    Evolutionary models used in BEAST analyses. Download TABLE S4, PDF file, 0.03 MB.

    Copyright © 2017 Tohma et al.

    This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

PreviousNext
Back to top
Download PDF
Citation Tools
Phylogenetic Analyses Suggest that Factors Other Than the Capsid Protein Play a Role in the Epidemic Potential of GII.2 Norovirus
Kentaro Tohma, Cara J. Lepore, Lauren A. Ford-Siltz, Gabriel I. Parra
mSphere May 2017, 2 (3) e00187-17; DOI: 10.1128/mSphereDirect.00187-17

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSphere article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Phylogenetic Analyses Suggest that Factors Other Than the Capsid Protein Play a Role in the Epidemic Potential of GII.2 Norovirus
(Your Name) has forwarded a page to you from mSphere
(Your Name) thought you would be interested in this article in mSphere.
Share
Phylogenetic Analyses Suggest that Factors Other Than the Capsid Protein Play a Role in the Epidemic Potential of GII.2 Norovirus
Kentaro Tohma, Cara J. Lepore, Lauren A. Ford-Siltz, Gabriel I. Parra
mSphere May 2017, 2 (3) e00187-17; DOI: 10.1128/mSphereDirect.00187-17
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

calicivirus
noroviruses
phylogenetic analysis
transmissible gastroenteritis virus

Related Articles

Cited By...

About

  • About mSphere
  • Board of Editors
  • Policies
  • For Reviewers
  • For the Media
  • Embargo Policy
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSphereJ

@ASMicrobiology

       

 

Website feedback

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

Copyright © 2019 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5042