The Canadian VirusSeq Data Portal & Duotang: open resources for SARS-CoV-2 viral sequences and genomic epidemiology.

The COVID-19 pandemic led to a large global effort to sequence SARS-CoV-2 genomes from patient samples to track viral evolution and inform public health response. Millions of SARS-CoV-2 genome sequences have been deposited in global public repositories. The Canadian COVID-19 Genomics Network (CanCOGeN - VirusSeq), a consortium tasked with coordinating expanded sequencing of SARS-CoV-2 genomes across Canada early in the pandemic, created the Canadian VirusSeq Data Portal, with associated data pipelines and procedures, to support these efforts. The goal of VirusSeq was to allow open access to Canadian SARS-CoV-2 genomic sequences and enhanced, standardized contextual data that were unavailable in other repositories and that meet FAIR standards (Findable, Accessible, Interoperable and Reusable). In addition, the Portal data submission pipeline contains data quality checking procedures and appropriate acknowledgement of data generators that encourages collaboration. From inception to execution, the portal was developed with a conscientious focus on strong data governance principles and practices. Extensive efforts ensured a commitment to Canadian privacy laws, data security standards, and organizational processes. This Portal has been coupled with other resources like Viral AI and was further leveraged by the Coronavirus Variants Rapid Response Network (CoVaRR-Net) to produce a suite of continually updated analytical tools and notebooks. Here we highlight this Portal, including its contextual data not available elsewhere, and the 'Duotang', a web platform that presents key genomic epidemiology and modeling analyses on circulating and emerging SARS-CoV-2 variants in Canada. Duotang presents dynamic changes in variant composition of SARS-CoV-2 in Canada and by province, estimates variant growth, and displays complementary interactive visualizations, with a text overview of the current situation. The VirusSeq Data Portal and Duotang resources, alongside additional analyses and resources computed from the Portal (COVID-MVP, CoVizu), are all open-source and freely available. Together, they provide an updated picture of SARS-CoV-2 evolution to spur scientific discussions, inform public discourse, and support communication with and within public health authorities. They also serve as a framework for other jurisdictions interested in open, collaborative sequence data sharing and analyses.

Zooanthroponotic transmission of SARS-CoV-2 and host-specific viral mutations revealed by genome-wide phylogenetic analysis.

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is a generalist virus, infecting and evolving in numerous mammals, including captive and companion animals, free-ranging wildlife, and humans. Transmission among non-human species poses a risk for the establishment of SARS-CoV-2 reservoirs, makes eradication difficult, and provides the virus with opportunities for new evolutionary trajectories, including the selection of adaptive mutations and the emergence of new variant lineages. Here, we use publicly available viral genome sequences and phylogenetic analysis to systematically investigate the transmission of SARS-CoV-2 between human and non-human species and to identify mutations associated with each species. We found the highest frequency of animal-to-human transmission from mink, compared with lower transmission from other sampled species (cat, dog, and deer). Although inferred transmission events could be limited by sampling biases, our results provide a useful baseline for further studies. Using genome-wide association studies, no single nucleotide variants (SNVs) were significantly associated with cats and dogs, potentially due to small sample sizes. However, we identified three SNVs statistically associated with mink and 26 with deer. Of these SNVs, ~⅔ were plausibly introduced into these animal species from local human populations, while the remaining ~⅓ were more likely derived in animal populations and are thus top candidates for experimental studies of species-specific adaptation. Together, our results highlight the importance of studying animal-associated SARS-CoV-2 mutations to assess their potential impact on human and animal health.

Does exposure to different menstrual products affect the vaginal environment?

The vaginal ecosystem is a key component of women's health. It also represents an ideal system for ecologists to investigate the consequence of perturbations on species diversity and emerging properties between organizational levels. Here, we study how exposure to different types of menstrual products is linked to microbial, immunological, demographic, and behavioural measurements in a cohort of young adult women who reported using more often tampons (n = 107) or menstrual cups (n = 31). We first found that cup users were older and smoked less than tampon users. When analysing health indicators, we detected potential associations between cups use reporting and fungal genital infection. A multivariate analysis confirmed that in our cohort, reporting using cups over tampons was associated with the higher odds ratio to report a fungal genital infection diagnosis by a medical doctor within the last 3 months. We did not detect significant differences between groups in terms of their bacterial vaginal microbiota composition and found marginal differences in the level of expression of 20 cytokines. However, a multivariate analysis of these biological data identified some level of clustering based on the menstrual product type preferred (cups or tampons). These results suggest that exposure to different types of menstrual products could influence menstrual health. Larger studies and studies with a more powered setting are needed to assess the robustness of these associations and identify causal mechanisms.

Concomitant and productive genital infections by HSV-2 and HPV in two young women: A case report.

Human papillomaviruses (HPVs), the most oncogenic virus known to humans, are often associated with Herpes Simplex Virus-2 (HSV-2) infections. The involvement of the latter in cervical cancer is controversial but its long-term infections might modulate the mucosal microenvironment in a way that favors carcinogenesis. We know little about coinfections between HSV-2 and HPVs, and studying the immunological and microbiological dynamics in the early stages of these infections may help identify or rule out potential interactions. We report two cases of concomitant productive, although asymptomatic, HSV-2 and HPV infections in young women (aged 20 and 25). The women were followed up for approximately a year, with clinical visits every two months and weekly self-samples. We performed quantitative analyses of their HSV-2 and HPV viral loads, immunological responses (IgG and IgM antibodies and local cytokines expression profiles), vaginal microbiota composition, as well as demographic and behavior data. We detect interactions between virus loads, immune response, and the vaginal microbiota, which improve our understanding of HSV-2 and HPVs' coinfections and calls for further investigation with larger cohorts.

Population Genomics Approaches for Genetic Characterization of SARS-CoV-2 Lineages.

The genome of the Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2), the pathogen that causes coronavirus disease 2019 (COVID-19), has been sequenced at an unprecedented scale leading to a tremendous amount of viral genome sequencing data. To assist in tracing infection pathways and design preventive strategies, a deep understanding of the viral genetic diversity landscape is needed. We present here a set of genomic surveillance tools from population genetics which can be used to better understand the evolution of this virus in humans. To illustrate the utility of this toolbox, we detail an in depth analysis of the genetic diversity of SARS-CoV-2 in first year of the COVID-19 pandemic. We analyzed 329,854 high-quality consensus sequences published in the GISAID database during the pre-vaccination phase. We demonstrate that, compared to standard phylogenetic approaches, haplotype networks can be computed efficiently on much larger datasets. This approach enables real-time lineage identification, a clear description of the relationship between variants of concern, and efficient detection of recurrent mutations. Furthermore, time series change of Tajima's D by haplotype provides a powerful metric of lineage expansion. Finally, principal component analysis (PCA) highlights key steps in variant emergence and facilitates the visualization of genomic variation in the context of SARS-CoV-2 diversity. The computational framework presented here is simple to implement and insightful for real-time genomic surveillance of SARS-CoV-2 and could be applied to any pathogen that threatens the health of populations of humans and other organisms.

A small number of early introductions seeded widespread transmission of SARS-CoV-2 in Québec, Canada.

BACKGROUND: Québec was the Canadian province most impacted by COVID-19, with 401,462 cases as of September 24th, 2021, and 11,347 deaths due mostly to a very severe first pandemic wave. In April 2020, we assembled the Coronavirus Sequencing in Québec (CoVSeQ) consortium to sequence SARS-CoV-2 genomes in Québec to track viral introduction events and transmission within the province. METHODS: Using genomic epidemiology, we investigated the arrival of SARS-CoV-2 to Québec. We report 2921 high-quality SARS-CoV-2 genomes in the context of > 12,000 publicly available genomes sampled globally over the first pandemic wave (up to June 1st, 2020). By combining phylogenetic and phylodynamic analyses with epidemiological data, we quantify the number of introduction events into Québec, identify their origins, and characterize the spatiotemporal spread of the virus. RESULTS: Conservatively, we estimated approximately 600 independent introduction events, the majority of which happened from spring break until 2 weeks after the Canadian border closed for non-essential travel. Subsequent mass repatriations did not generate large transmission lineages (> 50 sequenced cases), likely due to mandatory quarantine measures in place at the time. Consistent with common spring break and "snowbird" destinations, most of the introductions were inferred to have originated from Europe via the Americas. Once introduced into Québec, viral lineage sizes were overdispersed, with a few lineages giving rise to most infections. Consistent with founder effects, the earliest lineages to arrive tended to spread most successfully. Fewer than 100 viral introductions arrived during spring break, of which 7-12 led to the largest transmission lineages of the first wave (accounting for 52-75% of all sequenced infections). These successful transmission lineages dispersed widely across the province. Transmission lineage size was greatly reduced after March 11th, when a quarantine order for returning travellers was enacted. While this suggests the effectiveness of early public health measures, the biggest transmission lineages had already been ignited prior to this order. CONCLUSIONS: Combined, our results reinforce how, in the absence of tight travel restrictions or quarantine measures, fewer than 100 viral introductions in a week can ensure the establishment of extended transmission chains.

Outbreak investigation of SARS-CoV-2 transmission in an emergency childcare centre.

OBJECTIVE: To characterize SARS-CoV-2 transmission following a COVID-19 outbreak in an emergency childcare centre (ECCC) in April 2020 in Quebec, Canada. METHODS: The study population consisted of all the children and employees who attended the ECCC as well as household contacts of the confirmed COVID-19 cases. Of the 120 individuals in the study, five cases were confirmed by epidemiological link and 25 were identified as COVID-19 by RT-PCR among which 19 were analyzed by viral whole genome sequencing. Descriptive epidemiology, social network visualization, and phylogenetic analysis were used to characterize viral transmission. RESULTS: Phylogenetic analysis identified two separate introductions of distinct lineages of SARS-CoV-2 and estimated an average effective reproductive number of Re = 1.9 (range 0.9-4.9) with a mean doubling time of 3.2 days (range 2.1-5.2). The first and most prevalent lineage was introduced by two asymptomatic children who were likely infected by their parent, a confirmed COVID-19 case working in a long-term care centre. Among infected household adults, attack rates were significantly higher in mothers than in fathers (risk ratio = 4.5; 95% CI 1.1-18.7). The extent of transmission makes it one of the largest documented outbreaks in a daycare in Canada. CONCLUSION: The analyses carried out showed the probable origin and direction of the transmission of the infection (adult-child, child-adult, and child-child), thus highlighting how asymptomatic children can efficiently transmit SARS-CoV-2.

RéSUMé: OBJECTIF: Caractériser la transmission du SRAS-CoV-2 à la suite d'une éclosion de COVID-19 dans un service de garde d'urgence en milieu scolaire (SGUMS) en avril 2020 au Québec, Canada. MéTHODES: La population à l'étude était composée de tous les enfants et employés ayant fréquenté le SGUMS ainsi que les contacts familiaux des cas confirmés de COVID-19. Sur les 120 personnes à l'étude, cinq cas ont été confirmés par lien épidémiologique et 25 par RT-PCR. Parmi ces derniers, 19 ont été analysés par séquençage viral du génome entier. La caractérisation de la transmission a été réalisée à l'aide d'analyses descriptives et phylogénétiques ainsi que de la visualisation de réseaux sociaux. RéSULTATS: L'analyse phylogénétique a identifié deux introductions de lignées distinctes du SRAS-CoV-2 et un taux de reproduction net Re = 1,9 (étendue 0,9­4,9) avec un temps moyen de doublement de 3,2 jours (étendue 2,1­5,2). La première lignée, et la plus répandue, a été introduite par deux enfants asymptomatiques qui ont probablement été infectés par leur parent, un travailleur de la santé atteint de COVID-19. Dans les noyaux familiaux, les taux d'attaque étaient significativement plus élevés chez les mères que chez les pères (rapport de risque = 4,5 ; IC à 95 % 1,1­18,7). L'ampleur de la transmission en fait de celle-ci la plus importante éclosion documentée dans un service de garde au Canada. CONCLUSION: Cette étude a permis de déterminer l'origine et la direction probables de la transmission de l'infection (adulte-enfant, enfant-adulte et enfant-enfant) et démontrer que les enfants asymptomatiques peuvent transmettre le SRAS-CoV-2.

Cytokine response following perturbation of the cervicovaginal milieu during HPV genital infection.

Human papillomaviruses (HPVs) are oncogenic viruses causing most cervical cancers. Highly prevalent in young, sexually active women, only a minority of HPV infections persist. To better characterize the immuno-modulatory impact of early HPV infections, we measured changes in a panel of 20 cytokines in cervicovaginal samples collected from young women who were tested for HPV and self-reported for genital inflammation and infection symptoms. Multi-factor statistical analyses revealed that increased IL-1Alpha and IL-12/IL-23p40 concentrations were associated with HPV infection and that macrophage inflammatory proteins were associated in particular with high-risk HPV infections. ClinicalTrials.gov identifier NCT02946346.

HPV cervical infections and serological status in vaccinated and unvaccinated women.

Understanding genital infections by Human papillomaviruses (HPVs) remains a major public health issue, especially in countries where vaccine uptake is low. We investigate HPV prevalence and antibody status in 150 women (ages 18 to 25) in Montpellier, France. At inclusion and one month later, cervical swabs, blood samples and questionnaires (for demographics and behavioural variables) were collected. Oncogenic, non-vaccine genotypes HPV51, HPV66, HPV53, and HPV52 were the most frequently detected viral genotypes overall. Vaccination status, which was well-balanced in the cohort, showed the strongest (protective) effect against HPV infections, with an associated odds ratio for alphapapillomavirus detection of 0.45 (95% confidence interval: [0.22;0.58]). We also identified significant effects of age, number of partners, body mass index, and contraception status on HPV detection and on coinfections. Type-specific IgG serological status was also largely explained by the vaccination status. IgM seropositivity was best explained by HPV detection at inclusion only. Finally, we identify a strong significant effect of vaccination on genotype prevalence, with a striking under-representation of HPV51 in vaccinated women. Variations in HPV prevalence correlate with key demographic and behavioural variables. The cross-protective effect of the vaccine against HPV51 merits further investigation.

Does diversity beget diversity in microbiomes?

Microbes are embedded in complex communities where they engage in a wide array of intra- and inter-specific interactions. The extent to which these interactions drive or impede microbiome diversity is not well understood. Historically, two contrasting hypotheses have been suggested to explain how species interactions could influence diversity. 'Ecological Controls' (EC) predicts a negative relationship, where the evolution or migration of novel types is constrained as niches become filled. In contrast, 'Diversity Begets Diversity' (DBD) predicts a positive relationship, with existing diversity promoting the accumulation of further diversity via niche construction and other interactions. Using high-throughput amplicon sequencing data from the Earth Microbiome Project, we provide evidence that DBD is strongest in low-diversity biomes, but weaker in more diverse biomes, consistent with biotic interactions initially favouring the accumulation of diversity (as predicted by DBD). However, as niches become increasingly filled, diversity hits a plateau (as predicted by EC).

Detecting within-host interactions from genotype combination prevalence data.

Parasite genetic diversity can provide information on disease transmission dynamics but most mathematical and statistical frameworks ignore the exact combinations of genotypes in infections. We introduce and validate a new method that combines explicit epidemiological modelling of coinfections and regression-Approximate Bayesian Computing (ABC) to detect within-host interactions. Using a susceptible-infected-susceptible (SIS) model, we show that, if sufficiently strong, within-host parasite interactions can be detected from epidemiological data. We also show that, in this simple setting, this detection is robust even in the face of some level of host heterogeneity in behaviour. These simulations results offer promising applications to analyse large datasets of multiple infection prevalence data, such as those collected for genital infections by Human Papillomaviruses (HPVs).

Natural history, dynamics, and ecology of human papillomaviruses in genital infections of young women: protocol of the PAPCLEAR cohort study.

INTRODUCTION: Human papillomaviruses (HPVs) are responsible for one-third of all cancers caused by infections. Most HPV studies focus on chronic infections and cancers, and we know little about the early stages of the infection. Our main objective is to better understand the course and natural history of cervical HPV infections in healthy, unvaccinated and vaccinated, young women, by characterising the dynamics of various infection-related populations (virus, epithelial cells, vaginal microbiota and immune effectors). Another objective is to analyse HPV diversity within hosts, and in the study population, in relation to co-factors (lifestyle characteristics, vaccination status, vaginal microbiota, human genetics). METHODS AND ANALYSIS: The PAPCLEAR study is a single center longitudinal study following 150 women, aged 18-25 years, for up to 2 years. Visits occur every 2 or 4 months (depending on HPV status) during which several variables are measured, such as behaviours (via questionnaires), vaginal pH, HPV presence and viral load (via qPCR), local concentrations of cytokines (via MesoScale Discovery technology) and immune cells (via flow cytometry). Additional analyses are outsourced, such as titration of circulating anti-HPV antibodies, vaginal microbiota sequencing (16S and ITS1 loci) and human genotyping. To increase the statistical power of the epidemiological arm of the study, an additional 150 women are screened cross-sectionally. Finally, to maximise the resolution of the time series, participants are asked to perform weekly self-samples at home. Statistical analyses will involve classical tools in epidemiology, genomics and virus kinetics, and will be performed or coordinated by the Centre National de la Recherche Scientifique (CNRS) in Montpellier. ETHICS AND DISSEMINATION: This study has been approved by the Comité de Protection des Personnes Sud Méditerranée I (reference number 2016-A00712-49); by the Comité Consultatif sur le Traitement de l'Information en matière de Recherche dans le domaine de la Santé (reference number 16.504); by the Commission Nationale Informatique et Libertés (reference number MMS/ABD/AR1612278, decision number DR-2016-488) and by the Agence Nationale de Sécurité du Médicament et des Produits de Santé (reference 20160072000007). Results will be published in preprint servers, peer-reviewed journals and disseminated through conferences. TRIAL REGISTRATION NUMBER: NCT02946346; Pre-results.

Modelling the evolution of viral oncogenesis.

Most human oncogenic viruses share several characteristics, such as being DNA viruses, having long (co)evolutionary histories with their hosts and causing either latent or chronic infections. They can reach high prevalences while causing relatively low case mortality, which makes them quite fit according to virulence evolution theory. After analysing the life histories of DNA oncoviruses, we use a mathematical modelling approach to investigate how the virus life cycle may generate selective pressures favouring or acting against oncogenesis at the within-host or at the between-host level. In particular, we focus on two oncoprotein activities, namely extending cell life expectancy and increasing cell proliferation rate. These have immediate benefits (increasing viral population size) but can be associated with fitness costs at the epidemiological level (increasing recovery rate or risk of cancer) thus creating evolutionary trade-offs. We interpret the results of our nested model in light of the biological features and identify future perspectives for modelling oncovirus dynamics and evolution. This article is part of the theme issue 'Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses'.

Epithelial stratification shapes infection dynamics.

Infections of stratified epithelia contribute to a large group of common diseases, such as dermatological conditions and sexually transmitted diseases. To investigate how epithelial structure affects infection dynamics, we develop a general ecology-inspired model for stratified epithelia. Our model allows us to simulate infections, explore new hypotheses and estimate parameters that are difficult to measure with tissue cell cultures. We focus on two contrasting pathogens: Chlamydia trachomatis and Human papillomaviruses (HPV). Using cervicovaginal parameter estimates, we find that key infection symptoms can be explained by differential interactions with the layers, while clearance and pathogen burden appear to be bottom-up processes. Cell protective responses to infections (e.g. mucus trapping) generally lowered pathogen load but there were specific effects based on infection strategies. Our modeling approach opens new perspectives for 3D tissue culture experimental systems of infections and, more generally, for developing and testing hypotheses related to infections of stratified epithelia.

Why Human Papillomavirus Acute Infections Matter.

Most infections by human papillomaviruses (HPVs) are `acute', that is non-persistent. Yet, for HPVs, as for many other oncoviruses, there is a striking gap between our detailed understanding of chronic infections and our limited data on the early stages of infection. Here we argue that studying HPV acute infections is necessary and timely. Focusing on early interactions will help explain why certain infections are cleared while others become chronic or latent. From a molecular perspective, descriptions of immune effectors and pro-inflammatory pathways during the initial stages of infections have the potential to lead to novel treatments or to improved handling algorithms. From a dynamical perspective, adopting concepts from spatial ecology, such as meta-populations or meta-communities, can help explain why HPV acute infections sometimes last for years. Furthermore, cervical cancer screening and vaccines impose novel iatrogenic pressures on HPVs, implying that anticipating any viral evolutionary response remains essential. Finally, hints at the associations between HPV acute infections and fertility deserve further investigation given their high, worldwide prevalence. Overall, understanding asymptomatic and benign infections may be instrumental in reducing HPV virulence.

Potential Impact of Sexual Transmission on Ebola Virus Epidemiology: Sierra Leone as a Case Study.

BACKGROUND: Sexual transmission of Ebola virus disease (EVD) 6 months after onset of symptoms has been recently documented, and Ebola virus RNA has been detected in semen of survivors up to 9 months after onset of symptoms. As countries affected by the 2013-2015 epidemic in West Africa, by far the largest to date, are declared free of Ebola virus disease (EVD), it remains unclear what threat is posed by rare sexual transmission events that could arise from survivors. METHODOLOGY/PRINCIPAL FINDINGS: We devised a compartmental mathematical model that includes sexual transmission from convalescent survivors: a SEICR (susceptible-exposed-infectious-convalescent-recovered) transmission model. We fitted the model to weekly incidence of EVD cases from the 2014-2015 epidemic in Sierra Leone. Sensitivity analyses and Monte Carlo simulations showed that a 0.1% per sex act transmission probability and a 3-month convalescent period (the two key unknown parameters of sexual transmission) create very few additional cases, but would extend the epidemic by 83 days [95% CI: 68-98 days] (p < 0.0001) on average. Strikingly, a 6-month convalescent period extended the average epidemic by 540 days (95% CI: 508-572 days), doubling the current length, despite an insignificant rise in the number of new cases generated. CONCLUSIONS/SIGNIFICANCE: Our results show that reductions in the per sex act transmission probability via abstinence and condom use should reduce the number of sporadic sexual transmission events, but will not significantly reduce the epidemic size and may only minimally shorten the length of time the public health community must maintain response preparedness. While the number of infectious survivors is expected to greatly decline over the coming months, our results show that transmission events may still be expected for quite some time as each event results in a new potential cluster of non-sexual transmission. Precise measurement of the convalescent period is thus important for planning ongoing surveillance efforts.

Could the human papillomavirus vaccines drive virulence evolution?

The human papillomavirus (HPV) vaccines hold great promise for preventing several cancers caused by HPV infections. Yet little attention has been given to whether HPV could respond evolutionarily to the new selection pressures imposed on it by the novel immunity response created by the vaccine. Here, we present and theoretically validate a mechanism by which the vaccine alters the transmission-recovery trade-off that constrains HPV's virulence such that higher oncogene expression is favoured. With a high oncogene expression strategy, the virus is able to increase its viral load and infected cell population before clearance by the vaccine, thus improving its chances of transmission. This new rapid cell-proliferation strategy is able to circulate between hosts with medium to high turnover rates of sexual partners. We also discuss the importance of better quantifying the duration of challenge infections and the degree to which a vaccinated host can shed virus. The generality of the models presented here suggests a wider applicability of this mechanism, and thus highlights the need to investigate viral oncogenicity from an evolutionary perspective.

Quantifying the epidemic spread of Ebola virus (EBOV) in Sierra Leone using phylodynamics.

Measuring epidemic parameters early in an outbreak is essential to inform control efforts. Using the viral genome sequence and collection date from 78 infections in the 2014 Ebola virus outbreak in Sierra Leone, we estimate key epidemiological parameters such as infectious period duration (approximately 71 hours) and date of the first case in Sierra Leone (approximately April 25th). We also estimate the effective reproduction number, Re, (approximately 1.26), which is the number of secondary infections effectively caused by an infected individual and accounts for public health control measures. This study illustrates that phylodynamics methods, applied during the initial phase of an outbreak on fewer and more easily attainable data, can yield similar estimates to count-based epidemiological studies.

Clinical applications of pathogen phylogenies.

Innovative sequencing techniques now allow the routine access of whole genomes of pathogens, generating vast amounts of data. Phylogenetic trees are a common method for synthesizing this information. Unfortunately, these trees are often seen only as a visual support to guide researchers, thus neglecting the value of employing phylogenetic techniques to perform hypothesis testing on clinical questions. These include investigating how a pathogen spreads within a patient, or whether the infection severity (often measured by virus load) is controlled by viral genetics. Advances in methodology mean the time is ripe for combining phylogenies with clinical data to better understand and fight infectious diseases.