Intra-host Diversity Research Articles

Organ-specific genome diversity of replication-competent SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is not always confined to the respiratory system, as it impacts people on a broad clinical spectrum from asymptomatic to severe systemic manifestations resulting in death. Further, accumulation of intra-host single nucleotide variants during prolonged SARS-CoV-2 infection may lead to emergence of variants of concern (VOCs). Still, information on virus infectivity and intra-host evolution across organs is sparse. We report a detailed virological analysis of thirteen postmortem coronavirus disease 2019 (COVID-19) cases that provides proof of viremia and presence of replication-competent SARS-CoV-2 in extrapulmonary organs of immunocompromised patients, including heart, kidney, liver, and spleen (NCT04366882). In parallel, we identify organ-specific SARS-CoV-2 genome diversity and mutations of concern N501Y, T1027I, and Y453F, while the patient had died long before reported emergence of VOCs. These mutations appear in multiple organs and replicate in Vero E6 cells, highlighting their infectivity. Finally, we show two stages of fatal disease evolution based on disease duration and viral loads in lungs and plasma. Our results provide insights about the pathogenesis and intra-host evolution of SARS-CoV-2 and show that COVID-19 treatment and hygiene measures need to be tailored to specific needs of immunocompromised patients, even when respiratory symptoms cease.

Complex infections in vivax malaria: the more you look, the more you find

The human malaria parasite Plasmodium vivax commonly causes complex multiclonal infections. Recently, Dia et al. have developed innovative methods for single-cell sequencing (SCS) of P. vivax infections by adapting an approach used previously for Plasmodium falciparum. Their studies provide fascinating new insights into P. vivax intrahost diversity and evolution.

Identification and Characterization of Swine Influenza Virus H1N1 Variants Generated in Vaccinated and Nonvaccinated, Challenged Pigs.

Influenza viruses represent a continuous threat to both animal and human health. The 2009 H1N1 A influenza pandemic highlighted the importance of a swine host in the adaptation of influenza viruses to humans. Nowadays, one of the most extended strategies used to control swine influenza viruses (SIVs) is the trivalent vaccine application, whose formulation contains the most frequently circulating SIV subtypes H1N1, H1N2, and H3N2. These vaccines do not provide full protection against the virus, allowing its replication, evolution, and adaptation. To better understand the main mechanisms that shape viral evolution, here, the SIV intra-host diversity was analyzed in samples collected from both vaccinated and nonvaccinated animals challenged with the H1N1 influenza A virus. Twenty-eight whole SIV genomes were obtained by next-generation sequencing, and differences in nucleotide variants between groups were established. Substitutions were allocated along all influenza genetic segments, while the most relevant nonsynonymous substitutions were allocated in the NS1 protein on samples collected from vaccinated animals, suggesting that SIV is continuously evolving despite vaccine application. Moreover, new viral variants were found in both vaccinated and nonvaccinated pigs, showing relevant substitutions in the HA, NA, and NP proteins, which may increase viral fitness under field conditions.

Twenty years of evolution and diversification of digitaria streak virus in Digitaria setigera.

Within the family Geminiviridae, the emergence of new species results from their high mutation and recombination rates. In this study, we report the variability and evolution of digitaria streak virus (DSV), a mastrevirus isolated in 1986 from the grass Digitaria setigera in an island of the Vanuatu archipelago. Viral DNA of DSV samples was amplified from D. setigera specimens, derived from the naturally infected original plant, which were propagated in different laboratories in France and Italy for more than 20 years. From the consensus sequences, the nucleotide substitution rate was estimated for the period between a sample and the original sequence published in 1987, as well as for the period between samples. In addition, the intra-host genetic complexity and diversity of 8 DSV populations with a total of 165 sequenced haplotypes was characterized. The evolutionary rate of DSV was estimated to be between 1.13 × 10−4 and 9.87 × 10−4 substitutions/site/year, within the ranges observed in other single-stranded DNA viruses and RNA viruses. Bioinformatic analyses revealed high variability and heterogeneity in DSV populations, which confirmed that mutant spectra are continuously generated and are organized as quasispecies. The analysis of polymorphisms revealed nucleotide substitution biases in viral genomes towards deamination and oxidation of single-stranded DNA. The differences in variability in each of the genomic regions reflected a dynamic and modular evolution in the mutant spectra that was not reflected in the consensus sequences. Strikingly, the most variable region of the DSV genome, encoding the movement protein, showed rapid fixation of the mutations in the consensus sequence and a concomitant dN/dS ratio of 6.130, which suggests strong positive selection in this region. Phylogenetic analyses revealed a possible divergence in three genetic lineages from the original Vanuatu DSV isolate.

Human cytomegalovirus multiple-strain infections and viral population diversity in haematopoietic stem cell transplant recipients analysed by high-throughput sequencing

Human cytomegalovirus (HCMV) is an important opportunistic pathogen in allogeneic haematopoietic stem cell transplant (HSCT) recipients. High-throughput sequencing of target-enriched libraries was performed to characterise the diversity of HCMV strains present in this high-risk group. Forty-four HCMV-DNA-positive plasma specimens (median viral input load 321 IU per library) collected at defined time points from 23 HSCT recipients within 80 days of transplantation were sequenced. The genotype distribution for 12 hypervariable HCMV genes and the number of HCMV strains present (i.e. single- vs. multiple-strain infection) were determined for 29 samples from 16 recipients. Multiple-strain infection was observed in seven of these 16 recipients, and five of these seven recipients had the donor (D)/recipient (R) HCMV-serostatus combination D + R + . A very broad range of genotypes was detected, with an intrahost composition that was generally stable over time. Multiple-strain infection was not associated with particular virological or clinical features, such as altered levels or duration of antigenaemia, development of acute graft-versus-host disease or increased mortality. In conclusion, despite relatively low viral plasma loads, a high frequency of multiple-strain HCMV infection and a high strain complexity were demonstrated in systematically collected clinical samples from this cohort early after HSCT. However, robust evaluation of the pathogenic role of intrahost viral diversity and multiple-strain infection will require studies enrolling larger numbers of recipients.

Genomic Epidemiology of SARS-CoV-2 in Pakistan

COVID-19 has swept globally and Pakistan is no exception. To investigate the initial introductions and transmissions of the SARS-CoV-2 in Pakistan, we performed the largest genomic epidemiology study of COVID-19 in Pakistan and generated 150 complete SARS-CoV-2 genome sequences from samples collected from March 16 to June 1, 2020. We identified a total of 347 mutated positions, 31 of which were over-represented in Pakistan. Meanwhile, we found over 1000 intra-host single-nucleotide variants (iSNVs). Several of them occurred concurrently, indicating possible interactions among them or coevolution. Some of the high-frequency iSNVs in Pakistan were not observed in the global population, suggesting strong purifying selections. The genomic epidemiology revealed five distinctive spreading clusters. The largest cluster consisted of 74 viruses which were derived from different geographic locations of Pakistan and formed a deep hierarchical structure, indicating an extensive and persistent nation-wide transmission of the virus that was probably attributed to a signature mutation (G8371T in ORF1ab) of this cluster. Furthermore, 28 putative international introductions were identified, several of which are consistent with the epidemiological investigations. In all, this study has inferred the possible pathways of introductions and transmissions of SARS-CoV-2 in Pakistan, which could aid ongoing and future viral surveillance and COVID-19 control.

Intra-host evolution during SARS-CoV-2 prolonged infection.

Long-term infection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) represents a challenge to virus dispersion and the control of coronavirus disease 2019 (COVID-19) pandemic. The reason why some people have prolonged infection and how the virus persists for so long are still not fully understood. Recent studies suggested that the accumulation of intra-host single nucleotide variants (iSNVs) over the course of the infection might play an important role in persistence as well as emergence of mutations of concern. For this reason, we aimed to investigate the intra-host evolution of SARS-CoV-2 during prolonged infection. Thirty-three patients who remained reverse transcription polymerase chain reaction (RT-PCR) positive in the nasopharynx for on average 18 days from the symptoms onset were included in this study. Whole-genome sequences were obtained for each patient at two different time points. Phylogenetic, populational, and computational analyses of viral sequences were consistent with prolonged infection without evidence of coinfection in our cohort. We observed an elevated within-host genomic diversity at the second time point samples positively correlated with cycle threshold (Ct) values (lower viral load). Direct transmission was also confirmed in a small cluster of healthcare professionals that shared the same workplace by the presence of common iSNVs. A differential accumulation of missense variants between the time points was detected targeting crucial structural and non-structural proteins such as Spike and helicase. Interestingly, longitudinal acquisition of iSNVs in Spike protein coincided in many cases with SARS-CoV-2 reactive and predicted T cell epitopes. We observed a distinguishing pattern of mutations over the course of the infection mainly driven by increasing A→U and decreasing G→A signatures. G→A mutations may be associated with RNA-editing enzyme activities; therefore, the mutational profiles observed in our analysis were suggestive of innate immune mechanisms of the host cell defense. Therefore, we unveiled a dynamic and complex landscape of host and pathogen interaction during prolonged infection of SARS-CoV-2, suggesting that the host’s innate immunity shapes the increase of intra-host diversity. Our findings may also shed light on possible mechanisms underlying the emergence and spread of new variants resistant to the host immune response as recently observed in COVID-19 pandemic.

Assessment of two-pool multiplex long-amplicon nanopore sequencing of SARS-CoV-2.

Genomic surveillance of severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) plays an important role in COVID‐19 pandemic control and elimination efforts, especially by elucidating its global transmission network and illustrating its viral evolution. The deployment of multiplex PCR assays that target SARS‐CoV‐2 followed by either massively parallel or nanopore sequencing is a widely‐used strategy to obtain genome sequences from primary samples. However, multiplex PCR‐based sequencing carries an inherent bias of sequencing depth among different amplicons, which may cause uneven coverage. Here we developed a two‐pool, long‐amplicon 36‐plex PCR primer panel with ~1000‐bp amplicon lengths for full‐genome sequencing of SARS‐CoV‐2. We validated the panel by assessing nasopharyngeal swab samples with a <30 quantitative reverse transcription PCR cycle threshold value and found that ≥90% of viral genomes could be covered with high sequencing depths (≥20% mean depth). In comparison, the widely‐used ARTIC panel yielded 79%–88% high‐depth genome regions. We estimated that ~5 Mbp nanopore sequencing data may ensure a >95% viral genome coverage with a ≥10‐fold depth and may generate reliable genomes at consensus sequence levels. Nanopore sequencing yielded false‐positive variations with frequencies of supporting reads <0.8, and the sequencing errors mostly occurred on the 5′ or 3′ ends of reads. Thus, nanopore sequencing could not elucidate intra‐host viral diversity.

Molecular Detection and Characterization of Blastocystis sp. and Enterocytozoon bieneusi in Cattle in Northern Spain.

Some enteric parasites causing zoonotic diseases in livestock have been poorly studied or even neglected. This is the case in stramenopile Blastocystis sp. and the microsporidia Enterocytozoon bieneusi in Spain. This transversal molecular epidemiological survey aims to estimate the prevalence and molecular diversity of Blastocystis sp. and E. bieneusi in cattle faecal samples (n = 336) in the province of Álava, Northern Spain. Initial detection of Blastocystis and E. bieneusi was carried out by polymerase chain reaction (PCR) and Sanger sequencing of the small subunit (ssu) rRNA gene and internal transcribed spacer (ITS) region, respectively. Intra-host Blastocystis subtype diversity was further investigated by next generation amplicon sequencing (NGS) of the ssu rRNA gene in those samples that tested positive by conventional PCR. Amplicons compatible with Blastocystis sp. and E. bieneusi were observed in 32.1% (108/336, 95% CI: 27.2–37.4%) and 0.6% (2/336, 95% CI: 0.0–1.4%) of the cattle faecal samples examined, respectively. Sanger sequencing produced ambiguous/unreadable sequence data for most of the Blastocystis isolates sequenced. NGS allowed the identification of 10 Blastocystis subtypes including ST1, ST3, ST5, ST10, ST14, ST21, ST23, ST24, ST25, and ST26. All Blastocystis-positive isolates involved mixed infections of 2–8 STs in a total of 31 different combinations. The two E. bieneusi sequences were confirmed as potentially zoonotic genotype BEB4. Our data demonstrate that Blastocystis mixed subtype infections are extremely frequent in cattle in the study area. NGS was particularly suited to discern underrepresented subtypes or mixed subtype infections that were undetectable or unreadable by Sanger sequencing. The presence of zoonotic Blastocystis ST1, ST3, and ST5, and E. bieneusi BEB4 suggest cross-species transmission and a potential risk of human infection/colonization.

In Silico Analysis of Dengue Virus Serotype 2 Mutations Detected at the Intrahost Level in Patients with Different Clinical Outcomes.

ABSTRACTIntrahost genetic diversity is thought to facilitate arbovirus adaptation to changing environments and hosts, and it may also be linked to viral pathogenesis. Intending to shed light on the viral determinants for severe dengue pathogenesis, we previously analyzed the DENV-2 intrahost genetic diversity in 68 patients clinically classified as dengue fever (n = 31), dengue with warning signs (n = 19), and severe dengue (n = 18), performing viral whole-genome deep sequencing from clinical samples with an amplicon-free approach. From it, we identified a set of 141 relevant mutations distributed throughout the viral genome that deserved further attention. Therefore, we employed molecular modeling to recreate three-dimensional models of the viral proteins and secondary RNA structures to map the mutations and assess their potential effects. Results showed that, in general lines, disruptive variants were identified primarily among dengue fever cases. In contrast, potential immune-escape variants were associated mainly with warning signs and severe cases, in line with the latter’s longer intrahost evolution times. Furthermore, several mutations were located on protein-surface regions, with no associated function. They could represent sites of further investigation, as the interaction of viral and host proteins is critical for both host immunomodulation and virus hijacking of the cellular machinery. The present analysis provides new information about the implications of the intrahost genetic diversity of DENV-2, contributing to the knowledge about the viral factors possibly involved in its pathogenesis within the human host. Strengthening our results with functional studies could allow many of these variants to be considered in the design of therapeutic or prophylactic compounds and the improvement of diagnostic assays.IMPORTANCE Previous evidence showed that intrahost genetic diversity in arboviruses may be linked to viral pathogenesis and that one or a few amino acid replacements within a single protein are enough to modify a biological feature of an RNA virus. To assess dengue virus serotype 2 determinants potentially involved in pathogenesis, we previously analyzed the intrahost genetic diversity of the virus in patients with different clinical outcomes and identified a set of 141 mutations that deserved further study. Thus, through a molecular modeling approach, we showed that disruptive variants were identified primarily among cases with mild dengue fever, while potential immune-escape variants were mainly associated with cases of greater severity. We believe that some of the variants pointed out in this study were attractive enough to be potentially considered in future intelligent designs of therapeutic or prophylactic compounds or the improvement of diagnostic tools. The present analysis provides new information about DENV-2 viral factors possibly involved in its pathogenesis within the human host.

Evolutionary Comparisons of Chelonid Alphaherpesvirus 5 (ChHV5) Genomes from Fibropapillomatosis-Afflicted Green (Chelonia mydas), Olive Ridley (Lepidochelys olivacea) and Kemp's Ridley (Lepidochelys kempii) Sea Turtles.

Simple SummaryOur research aims to unravel uncertainties relating to the genetic and viral causes of the debilitating sea turtle disease fibropapillomatosis, which affects all seven species of sea turtle. This disease is likely caused by an alphaherpesvirus (ChHV5) and an environmental trigger (e.g., pollution). Fibropapillomatosis is characterised by multiple benign tumours which grow on the skin, eyes and internal organs, and is becoming a threat to sea turtle conservation globally. ChHV5 research is crucial to better provide effective management and conservation of turtles from this disease. This study aimed to compare ChHV5 genomes between geographic regions and sea turtle species and observe how this virus has evolved and changed. ChHV5 genomes harboured differences within and between geographic regions (88–2793 single nucleotide polymorphisms (SNPs) per sequenced genome). Multiple ChHV5 genes were also found to be under varying selective pressures. Phylogenomic and phylogenetic analyses revealed grouping of the virus, mostly by geography rather than by species, and found differences in ChHV5 genomes between tumours from the same individual. This study pioneers the phylogenomic approach to ChHV5 research. This study provides the most comprehensive picture to-date of whole-genome inter-species ChHV5 diversity and provides important baseline ChHV5 genomic data for future comparisons.The spreading global sea turtle fibropapillomatosis (FP) epizootic is threatening some of Earth’s ancient reptiles, adding to the plethora of threats faced by these keystone species. Understanding this neoplastic disease and its likely aetiological pathogen, chelonid alphaherpesvirus 5 (ChHV5), is crucial to understand how the disease impacts sea turtle populations and species and the future trajectory of disease incidence. We generated 20 ChHV5 genomes, from three sea turtle species, to better understand the viral variant diversity and gene evolution of this oncogenic virus. We revealed previously underappreciated genetic diversity within this virus (with an average of 2035 single nucleotide polymorphisms (SNPs), 1.54% of the ChHV5 genome) and identified genes under the strongest evolutionary pressure. Furthermore, we investigated the phylogeny of ChHV5 at both genome and gene level, confirming the propensity of the virus to be interspecific, with related variants able to infect multiple sea turtle species. Finally, we revealed unexpected intra-host diversity, with up to 0.15% of the viral genome varying between ChHV5 genomes isolated from different tumours concurrently arising within the same individual. These findings offer important insights into ChHV5 biology and provide genomic resources for this oncogenic virus.

Difference of Intrahost Dynamics of the Second Human Pegivirus and Hepatitis C Virus in HPgV-2/HCV-Coinfected Patients.

BackgroundThe second human pegivirus (HPgV-2) and hepatitis C virus (HCV) belong to the Flaviviridae family and share some common genome features. However, the two viruses exhibit significantly different genetic diversity. The comparison of intrahost dynamics of HPgV-2 and HCV that mainly reflect virus-host interactions is needed to elucidate their intrahost difference of genetic diversity and the possible mechanisms.MethodsIntrahost single nucleotide variations (iSNVs) were identified by means of next-generation sequencing from both cross-sectional and longitudinal samples from HPgV-2- and HCV-coinfected patients. The levels of human cytokines were quantified in the patient before and after HCV elimination by the treatment of direct-acting antivirals (DAA).ResultsUnlike HCV, the viral sequences of HPgV-2 are highly conserved among HPgV-2-infected patients. However, iSNV analysis confirmed the intrahost variation or quasispecies of HPgV-2. Almost all iSNVs of HPgV-2 did not accumulate or transmit within host over time, which may explain the highly conserved HPgV-2 consensus sequence. Intrahost variation of HPgV-2 mainly causes nucleotide transition in particular at the 3rd codon position and synonymous substitutions, indicating purifying or negative selection posed by host immune system. Cytokine data further indicate that HPgV-2 infection alone may not efficiently stimulate innate immune responses since proinflammatory cytokine expression dramatically decreased with elimination of HCV.ConclusionThis study provided new insights into the intrahost genomic variations and evolutionary dynamics of HPgV-2 as well as the impact of host immune selection and virus polymerase on virus evolution. The different genetic diversity of HPgV-2 and HCV makes HPgV-2 a potential new model to investigate RNA virus diversity and the mechanism of viral polymerase in modulating virus replication.

2020 Hindsight: Should evolutionary virologists have expected the unexpected during a pandemic?

2020 Hindsight: Should evolutionary virologists have expected the unexpected during a pandemic?

Characterization of Acute HCV Infection and Transmission Networks in People Who Currently Inject Drugs in Catalonia: Usefulness of Dried Blood Spots.

Accurate identification of recent HCV infections is critical for tracing the extent and mechanisms of ongoing transmission. We aimed to validate dried blood spot (DBS) samples for the assessment of Hepatitis C virus (HCV) genetic diversity and to determine epidemiological parameters including incidence, determinants of acute infection, and phylogenetic clustering in people who inject drugs (PWID). HCV nonstructural protein 5B next-generation sequencing was performed from plasma and/or DBS in 220 viremic PWID from the HepCdetect II study. No significant differences were found in consensus sequences or Shannon entropy (SE) intrahost diversity estimate between paired plasma/DBS specimens. SE values were used to identify acute infections with 93.3% sensitivity (95% CI, 0.81-1.06) and 95.0% specificity (95% CI, 0.88-1.02) in a set of well-defined controls. An acute HCV infection (either primary infection or reinfection) was detected in 13.5% of viremic participants and was associated with age ≤30years (OR,8.09), injecting less than daily (OR,4.35), ≤5years of injected drug use (OR,3.43), sharing cocaine snorting straws (OR,2.89), and being unaware of their HCV status (OR,3.62). Annualized HCV incidence was estimated between 31 and 59/100 person-years. On phylogenetic analysis, 46.8% of viremic cases were part of a transmission pair or cluster; age ≤30years (OR,6.16), acute infection (OR,5.73), and infection with subtype 1a (OR,4.78) were independently associated with this condition. The results obtained from plasma and DBS characterize PWID with acute infection and those involved in ongoing HCV transmission and allow estimating incidence from cross-sectional data. This information is critical for the design and assessment of targeted harm reduction programs and test-and-treat interventions and to facilitate monitoring of HCV elimination in this key population.

Evidence of novel genetic variants of Grapevine rupestris stem pitting-associated virus and intra-host diversity in Indian grapevine cultivars

Evidence of novel genetic variants of Grapevine rupestris stem pitting-associated virus and intra-host diversity in Indian grapevine cultivars

Genetic diversity and Multilocus Sequence Typing Analysis of Bartonella henselae in domestic cats from Southeastern Brazil

Bartonella henselae is the causative agent for the infectious disease Cat Scratch Disease (CSD), which can be fatal. Domestic and wild felines are known to be its main mammal reservoirs. The present study aimed to investigate the occurrence and genetic diversity of Bartonella spp. in cats sampled in São Paulo (SP) and Minas Gerais (MG) States, Southeastern Brazil. Based on a quantitative real-time PCR (qPCR) assay, a Bartonella sp. nuoG gene fragment was detected in 39.9% (122/306) of the blood samples (46/151 cats of SP; 76/155 cats of MG). The blood samples were submitted to a pre-enrichment culture technique that allowed the detection of 12 additional positive samples, which showed to be negative in the qPCR using DNA blood samples as templates. Furthermore, five B. henselae isolates were obtained from qPCR-negative samples for both blood and pre-enrichment culture. Seven out of 24 Ctenocephalides felis fleas were positive for Bartonella spp. in the qPCR assay; 4/7 positive fleas were collected from Bartonella-negative cats. Twenty-three rpoB B. henselae cloned sequences were obtained from nine cats’ blood samples, showing the occurrence of 13 different genotypes. Median-joining network and SplitsTree distance analysis showed that the obtained sequences represented distinct B. henselae genotypes when compared to those previously deposited in GenBank. Intra-host diversity was found, since different rpoB genotypes of B. henselae were detected in individual single cats. Bartonella henselae isolates showed two allelic profiles (ST37 in cats from MG state and ST9 in SP state) by MLST (Multilocus Sequence Typing) based on sequencing of eight molecular markers. The present study is the first molecular report of Bartonella sp. in cats from Minas Gerais State. In summary, this body of work showed the occurrence of different B. henselae rpoB genotypes at an intra-reservoir host level. Based on qPCR from blood samples and pre-enrichment liquid culture and isolation, occurrence of 33.1% (50/151) and 56.8% (88/155) for Bartonella sp. was found in cats from SP and MG states, respectively. Two different allelic profiles of B. henselae were found in cats from the states of São Paulo (ST9) and Minas Gerais (ST37), suggesting a clonal evolution of Bartonellae in a certain geographical region.

HPV16 and HPV18 type-specific APOBEC3 and integration profiles in different diagnostic categories of cervical samples

Human papillomavirus (HPV) 16 and 18 are the most predominant types in cervical cancer. Only a small fraction of HPV infections progress to cancer, indicating that additional factors and genomic events contribute to the carcinogenesis, such as minor nucleotide variation caused by APOBEC3 and chromosomal integration.We analysed intra-host minor nucleotide variants (MNVs) and integration in HPV16 and HPV18 positive cervical samples with different morphology. Samples were sequenced using an HPV whole genome sequencing protocol TaME-seq. A total of 80 HPV16 and 51 HPV18 positive samples passed the sequencing depth criteria of 300× reads, showing the following distribution: non-progressive disease (HPV16 n = 21, HPV18 n = 12); cervical intraepithelial neoplasia (CIN) grade 2 (HPV16 n = 27, HPV18 n = 9); CIN3/adenocarcinoma in situ (AIS) (HPV16 n = 27, HPV18 n = 30); cervical cancer (HPV16 n = 5).Similar numbers of MNVs in HPV16 and HPV18 samples were observed for most viral genes, with the exception of HPV18 E4 with higher numbers across clinical categories. APOBEC3 signatures were observed in HPV16 lesions, while similar mutation patterns were not detected for HPV18. The proportion of samples with integration was 13% for HPV16 and 59% for HPV18 positive samples, with a noticeable portion located within or close to cancer-related genes.

Precision Tracing of Household Dengue Spread Using Inter- and Intra-Host Viral Variation Data, Kamphaeng Phet, Thailand.

Dengue control approaches are best informed by granular spatial epidemiology of these viruses, yet reconstruction of inter- and intra-household transmissions is limited when analyzing case count, serologic, or genomic consensus sequence data. To determine viral spread on a finer spatial scale, we extended phylogenomic discrete trait analyses to reconstructions of house-to-house transmissions within a prospective cluster study in Kamphaeng Phet, Thailand. For additional resolution and transmission confirmation, we mapped dengue intra-host single nucleotide variants on the taxa of these time-scaled phylogenies. This approach confirmed 19 household transmissions and revealed that dengue disperses an average of 70 m per day between households in these communities. We describe an evolutionary biology framework for the resolution of dengue transmissions that cannot be differentiated based on epidemiologic and consensus genome data alone. This framework can be used as a public health tool to inform control approaches and enable precise tracing of dengue transmissions.

INSaFLU-TELE-Vir: an open web-based bioinformatics suite for influenza and SARS-CoV-2 genome-based surveillance

A new era of virus surveillance is emerging based on the real-time monitoring of virus evolution at whole-genome scale (World Health Organization 2021). Although national and international health authorities have strongly recommended this technological transition, especially for influenza and SARS-CoV-2 (World Health Organization 2021, Revez et al. 2017), the implementation of genomic surveillance can be particularly challenging due to the lack of bioinformatics infrastructures and/or expertise to process and interpret next-generation sequencing (NGS) data (Oakeson et al. 2017). We developed and implemented INSaFLU-TELE-Vir platform (https://insaflu.insa.pt/) (Borges et al. 2018), which is an influenza- and SARS-CoV-2-oriented bioinformatics free web-based suite that handles primary NGS data (reads) towards the automatic generation of the main “genetic requests'' for effective and timely laboratory surveillance. By handling NGS data collected from any amplicon-based schema (making it applicable for other pathogens), INSaFLU-TELE-Vir enables any laboratory to perform multi-step and intensive bioinformatics analyses in a user-oriented manner without requiring advanced training. INSaFLU-TELE-Vir handles NGS data collected from distinct sequencing technologies (Illumina, Ion Torrent and Oxford Nanopore Technologies), with the possibility of constructing comparative analyses using different technologies. It gives access to user-restricted sample databases and project management, being a transparent and flexible tool specifically designed to automatically update project outputs as more samples are uploaded. Data integration is thus cumulative and scalable, fitting the need for both routine surveillance and outbreak investigation activities. The bioinformatics pipeline consists of six core steps: read quality analysis and improvement, human betacoronaviruses (including SARS-CoV-2 Pango lineages) and influenza type/subtype classification, mutation detection and consensus generation, coverage analysis, alignment/phylogeny, intra-host minor variant detection (and automatic detection of putative mixed infections). read quality analysis and improvement, human betacoronaviruses (including SARS-CoV-2 Pango lineages) and influenza type/subtype classification, mutation detection and consensus generation, coverage analysis, alignment/phylogeny, intra-host minor variant detection (and automatic detection of putative mixed infections). The multiple outputs are provided in nomenclature-stable and standardized formats that can be visualized and explored in situ or through multiple compatible downstream applications for fine-tuned data analysis. Novel features are being implemented into the INSaFLU-TELE-Vir bioinformatics toolkit as part of the OHEJP TELE-Vir (https://onehealthejp.eu/jrp-tele-vir/) project, including rapid detection of selected genotype-phenotype associations, and enhanced geotemporal data visualization. All the code is available in github (https://github.com/INSaFLU) with the possibility of a local docker installation (https://github.com/INSaFLU/docker). A detailed documentation and tutorial is also available (https://insaflu.readthedocs.io/en/latest/). In summary, INSaFLU supplies public health laboratories and researchers with an open and user-friendly framework, potentiating a strengthened and timely multi-country genome-based virus surveillance.

Deep amplicon sequencing highlights low intra-host genetic variability of Echinococcus multilocularis and high prevalence of the European-type haplotypes in coyotes and red foxes in Alberta, Canada.

Echinococcus multilocularis (Em) is a zoonotic parasite considered a global emergent pathogen. Recent findings indicate that the parasite is expanding its range in North America and that European-type haplotypes are circulating in western Canada. However, genetic analyses are usually conducted only on a few parasites out of thousands of individuals within each definitive host, likely underestimating the prevalence of less common haplotypes. Moreover, mixed infections with several mtDNA haplotypes in the same host have been reported, but their relative abundance within the host was never estimated. We aimed to 1) estimate the frequency of co-infections of different Em haplotypes in coyotes (Canis latrans) and red foxes (Vulpes vulpes) from western Canada and their relative abundance within the definitive hosts, 2) detect less prevalent haplotypes by sampling a larger proportion of the parasite subpopulation per host, and 3) investigate differences in the distribution of Em haplotypes in these main definitive hosts; foxes and coyotes. We extracted DNA from ~10% of the worm subpopulation per host (20 foxes and 47 coyotes) and used deep amplicon sequencing (NGS technology) on four loci, targeting the most polymorphic regions from the mitochondrial genes cox1 (814 bp), nad1 (344 bp), and cob (387 bp). We detected the presence of mixed infections with multiple Em haplotypes and with different Echinococcus species including Em and E. granulosus s.l. genotypes G8/G10, low intraspecific diversity of Em, and a higher abundance of the European-type haplotypes in both hosts. Our results suggest a population expansion of the European over the North American strain in Alberta and a limited distribution of some European-type haplotypes. Our findings indicate that deep amplicon sequencing represents a valuable tool to characterize Em in multiple hosts, to assess the current distribution and possible origins of the European strain in North America. The potential use of next-generation sequencing technologies is particularly important to understand the patterns of geographic expansion of this parasite.