Human papillomavirus type 16 (HPV16) is the most common cause of cervical cancer, but most infections are transient with lesions not progressing to cancer. There is a lack of specific biomarkers for early cancer risk stratification. This study aimed to explore the intrahost HPV16 genomic variation in longitudinal samples from HPV16-infected women with different cervical lesion severity (normal, low-grade, and high-grade). The TaME-seq deep sequencing protocol was used to generate whole genome HPV16 sequences of 102 samples collected over time from 40 individuals. Single nucleotide variants (SNVs) and intrahost SNVs (iSNVs) were identified in the viral genomes. A majority of individuals had a unique set of SNVs and these SNVs were stable over time. Overall, the number of iSNVs and APOBEC3-induced iSNVs were significantly lower in high-grade relative to normal and low-grade samples. A significant increase in the number of APOBEC3-induced iSNVs over time was observed for normal samples when compared to high-grade. Our results indicates that the lower incidence of iSNVs and APOBEC3-induced iSNVs in high-grade lesions may have implications for novel biomarkers discoveries, potentially aiding early stratification of HPV-induced cervical precancerous lesions.

BACKGROUND: Rhinovirus (RV) infections can progress from the upper (URT) to lower (LRT) respiratory tract in immunocompromised individuals, causing high rates of fatal pneumonia. Little is known about how RV evolves within hosts during infection.
METHODS: We sequenced RV complete genomes from 12 hematopoietic cell transplant patients with infection for up to 190 days from both URT (nasal wash, NW) and LRT (bronchoalveolar lavage, BAL). Metagenomic and amplicon next-generation sequencing were used to track the emergence and evolution of intrahost single nucleotide variants (iSNVs).
RESULTS: Identical RV intrahost populations in matched NW and BAL specimens indicated no genetic adaptation is required for RV to progress from URT to LRT. Coding iSNVs were 2.3-fold more prevalent in capsid over nonstructural genes. iSNVs modeled were significantly more likely to be found in capsid surface residues, but were not preferentially located in known RV-neutralizing antibody epitopes. Newly emergent, genotype-matched iSNV haplotypes from immunocompromised individuals in 2008-2010 could be detected in Seattle-area community RV sequences in 2020-2021.
CONCLUSIONS: RV infections in immunocompromised hosts can progress from URT to LRT with no specific evolutionary requirement. Capsid proteins carry the highest variability and emergent mutations can be detected in other, including future, RV sequences.

Genetic analysis of intra-host viral populations provides unique insight into pre-emergent mutations that may contribute to the genotype of future variants. Clinical samples positive for SARS-CoV-2 collected in California during the first months of the pandemic were sequenced to define the dynamics of mutation emergence as the virus became established in the state. Deep sequencing of 90 nasopharyngeal samples showed that many mutations associated with the establishment of SARS-CoV-2 globally were present at varying frequencies in a majority of the samples, even those collected as the virus was first detected in the US. A subset of mutations that emerged months later in consensus sequences were detected as subconsensus members of intra-host populations. Spike mutations P681H, H655Y, and V1104L were detected prior to emergence in variant genotypes, mutations were detected at multiple positions within the furin cleavage site, and pre-emergent mutations were identified in the nucleocapsid and the envelope genes. Because many of the samples had a very high depth of coverage, a bioinformatics pipeline, \"Mappgene\", was established that uses both iVar and LoFreq variant calling to enable identification of very low-frequency variants. This enabled detection of a spike protein deletion present in many samples at low frequency and associated with a variant of concern.

UNASSIGNED: Low frequency intrahost single nucleotide variants (iSNVs) of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) have been increasingly recognised as predictive indicators of positive selection. Particularly as growing numbers of SARS-CoV-2 variants of interest (VOI) and concern (VOC) emerge. However, the dynamics of subgenomic RNA (sgRNA) expression and its impact on genomic diversity and infection outcome remain poorly understood. This study aims to investigate and quantify iSNVs and sgRNA expression in single and longitudinally sampled cohorts over the course of mild and severe SARS-CoV-2 infection, benchmarked against an in vitro infection model.
UNASSIGNED: Two clinical cohorts of SARS-CoV-2 positive cases in New South Wales, Australia collected between March 2020 and August 2021 were sequenced. Longitudinal samples from cases hospitalised due to SARS-CoV-2 infection (severe) (n = 16) were analysed and compared with cases that presented with SARS-CoV-2 symptoms but were not hospitalised (mild) (n = 23). SARS-CoV-2 genomic diversity profiles were also examined from daily sampling of culture experiments for three SARS-CoV-2 variants (Lineage A, B.1.351, and B.1.617.2) cultured in VeroE6 C1008 cells (n = 33).
UNASSIGNED: Intrahost single nucleotide variants were detected in 83% (19/23) of the mild cohort cases and 100% (16/16) of the severe cohort cases. SNP profiles remained relatively fixed over time, with an average of 1.66 SNPs gained or lost, and an average of 4.2 and 5.9 low frequency variants per patient were detected in severe and mild infection, respectively. sgRNA was detected in 100% (25/25) of the mild genomes and 92% (24/26) of the severe genomes. Total sgRNA expressed across all genes in the mild cohort was significantly higher than that of the severe cohort. Significantly higher expression levels were detected in the spike and the nucleocapsid genes. There was significantly less sgRNA detected in the culture dilutions than the clinical cohorts.
UNASSIGNED: The positions and frequencies of iSNVs in the severe and mild infection cohorts were dynamic overtime, highlighting the importance of continual monitoring, particularly during community outbreaks where multiple SARS-CoV-2 variants may co-circulate. sgRNA levels can vary across patients and the overall level of sgRNA reads compared to genomic RNA can be less than 1%. The relative contribution of sgRNA to the severity of illness warrants further investigation given the level of variation between genomes. Further monitoring of sgRNAs will improve the understanding of SARS-CoV-2 evolution and the effectiveness of therapeutic and public health containment measures during the pandemic.

With deep sequencing of virus genomes within the hosts, intrahost single nucleotide variations (iSNVs) have been used for analyses of virus genome variation and evolution, which is indicated to correlate with viral pathogenesis and disease severity. Little is known about the features of iSNVs among DNA viruses. We performed the epidemiological and laboratory investigation of one outbreak of adenovirus. The whole genomes of viruses in both original oral swabs and cell-cultured virus isolates were deeply sequenced. We identified 737 iSNVs in the viral genomes sequenced from original samples and 46 viral iSNVs in cell-cultured isolates, with 33 iSNVs shared by original samples and cultured isolates. Meanwhile, we found these 33 iSNVs were shared by different patients, among which, three hot spot areas 6367-6401, 9213-9247, and 10 584-10 606 within the functional genes of the adenovirus genome were found. Notably, the substitution rates of iSNVs were closely correlated with the clinical and immune indicators of the patients. Especially a positive correlation to neutrophils was found, indicating a predictable biomarker of iSNV dynamics. Our findings demonstrated the neutrophil-correlated dynamic evolution features of the iSNVs within adenoviruses, which indicates a virus-host interaction during human infection of a DNA virus.

The novel coronavirus (SARS-CoV-2) has become a pandemic and is threatening human health globally. Here, we report nine newly evolved SARS-CoV-2 single nucleotide polymorphism (SNP) alleles those underwent a rapid increase (seven cases) or decrease (two cases) in their frequency for 30-80% in the initial four months, which are further confirmed by intrahost single nucleotide variation analysis using raw sequence data including 8,217 samples. The nine SNPs are mostly (8/9) located in the coding region and are mainly (6/9) nonsynonymous substitutions. The nine SNPs show a complete linkage in SNP pairs and belong to three different linkage groups, named LG_1 to LG_3. Analyses in population genetics show signatures of adaptive selection toward the mutants in LG_1, but no signal of selection for LG_2. Population genetic analysis results on LG_3 show geological differentiation. Analyses on geographic COVID-19 cases and published clinical data provide evidence that the mutants in LG_1 and LG_3 benefit virus replication and those in LG_1 have a positive correlation with the disease severity in COVID-19-infected patients. The mutants in LG_2 show a bias toward mildness of the disease based on available public clinical data. Our findings may be instructive for epidemiological surveys and disease control of COVID-19 in the future.

COVID-19 (coronavirus disease 2019) has caused a major epidemic worldwide; however, much is yet to be known about the epidemiology and evolution of the virus partly due to the scarcity of full-length SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) genomes reported. One reason is that the challenges underneath sequencing SARS-CoV-2 directly from clinical samples have not been completely tackled, i.e., sequencing samples with low viral load often results in insufficient viral reads for analyses.
We applied a novel multiplex PCR amplicon (amplicon)-based and hybrid capture (capture)-based sequencing, as well as ultra-high-throughput metatranscriptomic (meta) sequencing in retrieving complete genomes, inter-individual and intra-individual variations of SARS-CoV-2 from serials dilutions of a cultured isolate, and eight clinical samples covering a range of sample types and viral loads. We also examined and compared the sensitivity, accuracy, and other characteristics of these approaches in a comprehensive manner.
We demonstrated that both amplicon and capture methods efficiently enriched SARS-CoV-2 content from clinical samples, while the enrichment efficiency of amplicon outran that of capture in more challenging samples. We found that capture was not as accurate as meta and amplicon in identifying between-sample variations, whereas amplicon method was not as accurate as the other two in investigating within-sample variations, suggesting amplicon sequencing was not suitable for studying virus-host interactions and viral transmission that heavily rely on intra-host dynamics. We illustrated that meta uncovered rich genetic information in the clinical samples besides SARS-CoV-2, providing references for clinical diagnostics and therapeutics. Taken all factors above and cost-effectiveness into consideration, we proposed guidance for how to choose sequencing strategy for SARS-CoV-2 under different situations.
This is, to the best of our knowledge, the first work systematically investigating inter- and intra-individual variations of SARS-CoV-2 using amplicon- and capture-based whole-genome sequencing, as well as the first comparative study among multiple approaches. Our work offers practical solutions for genome sequencing and analyses of SARS-CoV-2 and other emerging viruses.