Within-host evolutionary dynamics of seasonal and pandemic human influenza A viruses in young children
- Department of Medical Microbiology & Infection Prevention, Amsterdam University Medical Center, Netherlands
- National Institute of Hygiene and Epidemiology, Viet Nam
- Ha Nam Centre for Disease Control, Viet Nam
- Children's Hospital 2, Viet Nam
- Children's Hospital 1, Viet Nam
- Vietnam National Children's Hospital, Viet Nam
- Siriraj Hospital, Mahidol University, Thailand
- Hospital for Tropical Diseases, Viet Nam
- National Hospital for Tropical Diseases, Viet Nam
- Oxford University Clinical Research Unit, Viet Nam
- Radboud Medical Centre, Radboud University, Netherlands
- Nuffield Department of Medicine, University of Oxford, United Kingdom
- Peter Doherty Institute for Infection and Immunity, University of Melbourne, Australia
- WHO Collaborating Centre for Reference and Research on Influenza, Australia
- Centre for Infectious Disease Control, National Institute for Public Health and the Environment, Netherlands
Figures
Figure 1 with 5 supplements
Genetic diversity of within-host influenza A virus populations.
Box plots summarizing the number of intra-host single-nucleotide variants (iSNVs; median, interquartile range [IQR], and whiskers extending within median ±1.5 × IQR) identified in samples with …
Figure 1—figure supplement 1
Sequence coverage across all influenza gene segments and samples.
Black line plots the mean coverage for a sliding window of 50 base pairs (stepsize = 25 base pairs). The interquartile range is shaded in dark pink while the full range is denoted in light pink. (A) …
Figure 1—figure supplement 2
Frequencies of nucleotide variants found in A/H3N2 viral reads sequenced from overlapping amplicons.
Each circle represents a nucleotide variant site (with frequency estimated between 0.02 and 0.98) found in reads attributed to at least two different amplicons (at least 100× coverage for each …
Figure 1—figure supplement 3
Maximum-likelihood phylogeny of putative majority (consensus) and minority whole-genome sequences (by concatenating all eight gene segments) of A/H3N2 virus samples.
Tip names are given in the format: ‘Patient ID_Days since symptom onset_putative consensus or minority sequence.’ The tree is rooted to the A/Brisbane/10/2007 virus (H3N2_Bris07; EPI_ISL_103644). …
Figure 1—figure supplement 4
Maximum-likelihood phylogeny of putative majority (consensus) and minority whole-genome sequences (by concatenating all eight gene segments) of H1N1pdm09 virus samples.
Tip names are given in the format: ‘Patient ID_Days since symptom onset_putative consensus or minority sequence.’ The tree is rooted to the A/California/04/2009 virus (H1N1pdm09_Cali09; …
Figure 1—figure supplement 5
Pearson’s correlation between the first day of oseltamivir treatment administered to patients and the last day on which viral samples with cycle threshold (CT) values ≤35 were collected.
Time is measured by number of days since symptom onset. Each point represents a patient included in this study who was treated with oseltamivir (Supplementary file 4).
Figure 2 with 4 supplements
Box plots (median, interquartile range [IQR], and whiskers extending within median ±1.5 × IQR) summarizing the empirical within-host evolutionary rates of (A) seasonal A/H3N2 viruses and (B) pandemic A/H1N1pdm09 viruses.
Top panel shows the evolutionary rate of individual gene segments over all timepoints () while the bottom panel depicts the genome-wide evolutionary rate () for each day since symptom onset. All …
Figure 2—figure supplement 1
Box plots (median, interquartile range [IQR], and whiskers extending within median ±1.5 × IQR) summarizing the empirical within-host evolutionary rates () of different H3N2 viral gene segments.
All rates are stratified by substitution type (synonymous – blue; nonsynonymous – red; stop codon – gray). Wilcoxon signed-rank tests were performed to assess if the paired synonymous and …
Figure 2—figure supplement 2
Box plots (median, interquartile range [IQR], and whiskers extending within median ±1.5 × IQR) summarizing the empirical within-host evolutionary rates () of different H1N1pdm09 viral gene segments.
All rates are stratified by substitution type (synonymous – blue; nonsynonymous – red; stop codon – gray). Wilcoxon signed-rank test was not performed here due to the low number of samples collected …
Figure 2—figure supplement 3
Linear regression of within-host synonymous and nonsynonymous evolutionary rates of within-host A/H3N2 virus samples.
Each plotted line is the linearly regressed line to the evolutionary rates computed for each A/H3N2-infected individual. Based on our findings, we expect that synonymous rates correlate negatively …
Figure 2—figure supplement 4
Box plots (median, interquartile range [IQR], and whiskers extending within median ±1.5 × IQR) summarizing the empirical daily within-host evolutionary rates of seasonal A/H3N2 viruses.
Variants that could potentially be PCR artifacts were removed (i.e., those found under the 75th percentile [6.3%] of frequency range of variants located in overlapping amplicons but were only …
Figure 3 with 2 supplements
Histogram of the mean number of minority intra-host single-nucleotide variants (iSNVs) identified per sample across all.
(A) A/H3N2 and (B) A/H1N1pdm09 virus specimens, sorted by frequency bins of 5% and substitution type (synonymous – blue; nonsynonymous – red; stop-codon – gray).
Figure 3—figure supplement 1
Histogram of the mean number of minority intra-host single-nucleotide variants (iSNVs) identified per sample across all H3N2 viral gene segments across all samples sorted by frequency bins of 5% and substitution type (synonymous – blue; nonsynonymous – red; gray – stop codon).
Figure 3—figure supplement 2
Histogram of the mean number of minority intra-host single-nucleotide variants (iSNVs) identified across all H1N1pdm09 viral gene segments across all samples sorted by frequency bins of 5% and substitution type (synonymous – blue; nonsynonymous – red; stop-codon – gray).
Figure 4 with 9 supplements
Intra-host single-nucleotide variants in within-host IAV populations.
(A) Breakdown of intra-host single-nucleotide variants (iSNVs) identified in seasonal A/H3N2 virus samples. The top panels plot the nucleotide positions where iSNVs were found in at least two …
Figure 4—figure supplement 1
Plots of intra-host hemagglutinin (HA) amino acid variants in A/H3N2-infected individuals.
Top panel shows the number of subjects where nonsynonymous variants were found in the respective protein site. Different canonical antigenic sites of the HA protein are colored (HA numbering based …
Figure 4—figure supplement 2
Plots of intra-host neuraminidase (NA) amino acid variants in A/H3N2-infected individuals.
Top panel shows the number of subjects where nonsynonymous variants were found in the respective protein site. Bottom panel plots selected as well as parallel amino acid mutations found in multiple …
Figure 4—figure supplement 3
Plots of intra-host nucleoprotein (NP) amino acid variants in A/H3N2-infected individuals.
Top panel shows the number of subjects where nonsynonymous variants were found in the respective protein site. Bottom panel plots selected as well as parallel amino acid mutations found in multiple …
Figure 4—figure supplement 4
Plots of M2 protein intra-host amino acid variants in A/H3N2-infected individuals.
Top panel shows the number of subjects where nonsynonymous variants were found in the respective protein site. Bottom panel plots selected as well as parallel amino acid mutations found in multiple …
Figure 4—figure supplement 5
Plots of hemagglutinin (HA) intra-host amino acid variants in A/H1N1pdm09-infected individuals.
Top panel shows the number of subjects where nonsynonymous variants were found in the respective protein site. Different canonical antigenic sites of the HA protein are colored (HA numbering based …
Figure 4—figure supplement 6
Plots of neuraminidase (NA) intra-host amino acid variants in A/H1N1pdm09-infected individuals.
Top panel shows the number of subjects where nonsynonymous variants were found in the respective protein site. Bottom panel plots selected as well as parallel amino acid mutations found in multiple …
Figure 4—figure supplement 7
Plots of M2 protein intra-host amino acid variants in A/H1N1pdm09-infected individuals.
Top panel shows the number of subjects where nonsynonymous variants were found in the respective protein site. Bottom panel plots selected as well as parallel amino acid mutations found in multiple …
Figure 4—figure supplement 8
Frequency distributions of intra-host single-nucleotide variants (iSNVs) below the 2% variant calling threshold found in nucleotide positions NP-1150 and M-917 that encode for amino acid sites NP-384 and M2-77, respectively.
All A/H3N2 virus samples collected from all patients with site coverage above the 100× are included. The distributions were compared to that of neighboring sites, ±10 nucleotide positions adjacent …
Figure 4—figure supplement 9
Plots of within-host recurring A/H3N2 amino acid variants NP-G384R and M2-R77* based on variant calls and frequencies after remapping sample reads to their respective sample consensus sequence.
Figure 5
The trimeric and monomeric crystal structures of nucleoprotein (PDB: 3ZDP) (Chenavas et al., 2013) of influenza A viruses.
Amino acid sites with potentially linked epistatic amino acid variants as tabulated in Table 1 are separately colored, with their corresponding positions annotated on the monomeric structure.
Figure 6 with 1 supplement
Evolutionary rates computed from forward-time Monte Carlo within-host simulations for different mean deleterious effects () of nonsynonymous mutations.
We assumed that synonymous mutations are neutral for all simulations. (A) Neutral expectation where all nonsynonymous mutations are neutral (). We tested our hypotheses where the majority of …
Figure 6—figure supplement 1
Number of virions () against replicative generation () based on a target cell-limited within-host model.
Blue line with markers denotes the population size computed from the model. When virions, we assumed that remained constant at 107 (pink dashed line) to reduce computational costs of …
Appendix 1—figure 1
Genetic diversity of within-host influenza A virus populations as estimated by nucleotide diversity statistic.
Box plots summarizing the statistic (intra-host single-nucleotide variants [iSNVs]; median, interquartile range [IQR], and whiskers extending within median ±1.5 × IQR) computed for samples with …
Appendix 1—figure 2
A/H1N1pdm09 virus household transmission pairs.
(A) Schematic of A/H1N1pdm09 virus household transmission pairs identified by epidemiological linkage and plotted based on timing of sample collection (left panel) Intra-host single-nucleotide …
Tables
Table 1
FoldX stability predictions of likely linked nonsynonymous minority variants found in A/H3N2 nucleoprotein.
The mean and standard deviation (SD) values reported are based on the results of five distinct simulations. Variants with mean kcal/mol are deemed to be stabilizing while destabilizing mutants …
| (kcal/mol) | ||
|---|---|---|
| Variants | Mean | SD |
| G384R | −3.84 | 0.06 |
| M426I | 2.61 | 0.05 |
| G384R,M426I | −0.42 | 0.06 |
| G102R | 4.87 | 0.00 |
| G384R,G102R | 0.76 | 0.09 |
| A493T | 11.96 | 0.30 |
| G384R,A493T | 5.56 | 0.19 |
| V197I | −3.11 | 0.02 |
| S353Y | −1.97 | 0.68 |
| V197I,S353Y | −4.48 | 0.14 |
Table 2
Parameter values used in the within-host model.
| Parameter | Meaning | Value (units) | Source |
|---|---|---|---|
| - | Number of hours per replicative generation | 6 hr | Assumption |
| Average number of virions produced by an infected cell | 100 virions | Frensing et al., 2016 | |
| Initial target cell population size | 4 × 108 virions | Hadjichrysanthou et al., 2016 | |
| Per capita decay rate | Two per-generation | Assumption | |
| Within-host basic reproduction number | 5 | Hadjichrysanthou et al., 2016 | |
| µ | Per-site, per-generation mutation rate | 3 × 10−5 per-site, per-generation | McCrone et al., 2020 |
Additional files
-
Supplementary file 1
Mean number of nonsynonymous (NS), synonymous (S), and stop codon (Stop) variants per sample for each gene segment as well as the corresponding NS/S ratio.
- https://cdn.elifesciences.org/articles/68917/elife-68917-supp1-v2.xlsx
-
Supplementary file 2
Potentially linked nonsynonymous variants in within-host A/H1N1pdm09 and A/H3N2 virus samples.
Sample names are given in the format of 'Patient ID_Days since symptom onset.' Both linkage disequilibrium () and the normalized measures are tabulated alongside the inferred maximum-likelihood haplotype frequencies ( and are the haplotype frequencies with variant i or ii only while is the frequency of haplotypes encoding both variants).
- https://cdn.elifesciences.org/articles/68917/elife-68917-supp2-v2.xlsx
-
Supplementary file 3
A/H3N2 segment-specific primers.
- https://cdn.elifesciences.org/articles/68917/elife-68917-supp3-v2.xlsx
-
Supplementary file 4
Patients metadata (provided as an Excel file).
- https://cdn.elifesciences.org/articles/68917/elife-68917-supp4-v2.xlsx
-
Supplementary file 5
Acknowledgment table of reference sequences downloaded from GISAID.
- https://cdn.elifesciences.org/articles/68917/elife-68917-supp5-v2.xlsx
-
Transparent reporting form
- https://cdn.elifesciences.org/articles/68917/elife-68917-transrepform-v2.docx
