Scalable Bayesian Divergence Time Estimation With Ratio Transformations.

Divergence time estimation is crucial to provide temporal signals for dating biologically important events from species divergence to viral transmissions in space and time. With the advent of high-throughput sequencing, recent Bayesian phylogenetic studies have analyzed hundreds to thousands of sequences. Such large-scale analyses challenge divergence time reconstruction by requiring inference on highly correlated internal node heights that often become computationally infeasible. To overcome this limitation, we explore a ratio transformation that maps the original $N-1$ internal node heights into a space of one height parameter and $N-2$ ratio parameters. To make the analyses scalable, we develop a collection of linear-time algorithms to compute the gradient and Jacobian-associated terms of the log-likelihood with respect to these ratios. We then apply Hamiltonian Monte Carlo sampling with the ratio transform in a Bayesian framework to learn the divergence times in 4 pathogenic viruses (West Nile virus, rabies virus, Lassa virus, and Ebola virus) and the coralline red algae. Our method both resolves a mixing issue in the West Nile virus example and improves inference efficiency by at least 5-fold for the Lassa and rabies virus examples as well as for the algae example. Our method now also makes it computationally feasible to incorporate mixed-effects molecular clock models for the Ebola virus example, confirms the findings from the original study, and reveals clearer multimodal distributions of the divergence times of some clades of interest.

Emergence of the B.1.214.2 SARS-CoV-2 lineage with an Omicron-like spike insertion and a unique upper airway immune signature.

We investigate the emergence, mutation profile, and dissemination of SARS-CoV-2 lineage B.1.214.2, first identified in Belgium in January 2021. This variant, featuring a 3-amino acid insertion in the spike protein similar to the Omicron variant, was speculated to enhance transmissibility or immune evasion. Initially detected in international travelers, it substantially transmitted in Central Africa, Belgium, Switzerland, and France, peaking in April 2021. Our travel-aware phylogeographic analysis, incorporating travel history, estimated the origin to the Republic of the Congo, with primary European entry through France and Belgium, and multiple smaller introductions during the epidemic. We correlate its spread with human travel patterns and air passenger data. Further, upon reviewing national reports of SARS-CoV-2 outbreaks in Belgian nursing homes, we found this strain caused moderately severe outcomes (8.7% case fatality ratio). A distinct nasopharyngeal immune response was observed in elderly patients, characterized by 80% unique signatures, higher B- and T-cell activation, increased type I IFN signaling, and reduced NK, Th17, and complement system activation, compared to similar outbreaks. This unique immune response may explain the variant's epidemiological behavior and underscores the need for nasal vaccine strategies against emerging variants.

Gene Balance Predicts Transcriptional Responses Immediately Following Ploidy Change in Arabidopsis thaliana.

The gene balance hypothesis postulates that there is selection on gene copy number (gene dosage) to preserve the stoichiometric balance among interacting proteins. This presupposes that gene product abundance is governed by gene dosage and that gene dosage responses are consistent for interacting genes in a dosage-balance-sensitive network or complex. Gene dosage responses, however, have rarely been quantified, and the available data suggest that they are highly variable. We sequenced the transcriptomes of two synthetic autopolyploid accessions of Arabidopsis (Arabidopsis thaliana) and their diploid progenitors, as well as one natural tetraploid and its synthetic diploid produced via haploid induction, to estimate transcriptome size and dosage responses immediately following ploidy change. Similar to what has been observed in previous studies, overall transcriptome size does not exhibit a simple doubling in response to genome doubling, and individual gene dosage responses are highly variable in all three accessions, indicating that expression is not strictly coupled with gene dosage. Nonetheless, putatively dosage balance-sensitive gene groups (Gene Ontology terms, metabolic networks, gene families, and predicted interacting proteins) exhibit smaller and more coordinated dosage responses than do putatively dosage-insensitive gene groups, suggesting that constraints on dosage balance operate immediately following whole-genome duplication and that duplicate gene retention patterns are shaped by selection to preserve dosage balance.

Helicobacter pylori diversification during chronic infection within a single host generates sub-populations with distinct phenotypes.

Helicobacter pylori chronically infects the stomach of approximately half of the world's population. Manifestation of clinical diseases associated with H. pylori infection, including cancer, is driven by strain properties and host responses; and as chronic infection persists, both are subject to change. Previous studies have documented frequent and extensive within-host bacterial genetic variation. To define how within-host diversity contributes to phenotypes related to H. pylori pathogenesis, this project leverages a collection of 39 clinical isolates acquired prospectively from a single subject at two time points and from multiple gastric sites. During the six years separating collection of these isolates, this individual, initially harboring a duodenal ulcer, progressed to gastric atrophy and concomitant loss of acid secretion. Whole genome sequence analysis identified 1,767 unique single nucleotide polymorphisms (SNPs) across isolates and a nucleotide substitution rate of 1.3x10-4 substitutions/site/year. Gene ontology analysis identified cell envelope genes among the genes with excess accumulation of nonsynonymous SNPs (nSNPs). A maximum likelihood tree based on genetic similarity clusters isolates from each time point separately. Within time points, there is segregation of subgroups with phenotypic differences in bacterial morphology, ability to induce inflammatory cytokines, and mouse colonization. Higher inflammatory cytokine induction in recent isolates maps to shared polymorphisms in the Cag PAI protein, CagY, while rod morphology in a subgroup of recent isolates mapped to eight mutations in three distinct helical cell shape determining (csd) genes. The presence of subgroups with unique genetic and phenotypic properties suggest complex selective forces and multiple niches within the stomach during chronic infection.

Nextstrain: real-time tracking of pathogen evolution.

Summary: Understanding the spread and evolution of pathogens is important for effective public health measures and surveillance. Nextstrain consists of a database of viral genomes, a bioinformatics pipeline for phylodynamics analysis, and an interactive visualization platform. Together these present a real-time view into the evolution and spread of a range of viral pathogens of high public health importance. The visualization integrates sequence data with other data types such as geographic information, serology, or host species. Nextstrain compiles our current understanding into a single accessible location, open to health professionals, epidemiologists, virologists and the public alike. Availability and implementation: All code (predominantly JavaScript and Python) is freely available from github.com/nextstrain and the web-application is available at nextstrain.org.

Genomic epidemiology supports multiple introductions and cryptic transmission of Zika virus in Colombia.

BACKGROUND: Colombia was the second most affected country during the American Zika virus (ZIKV) epidemic, with over 109,000 reported cases. Despite the scale of the outbreak, limited genomic sequence data were available from Colombia. We sought to sequence additional samples and use genomic epidemiology to describe ZIKV dynamics in Colombia. METHODS: We sequenced ZIKV genomes directly from clinical diagnostic specimens and infected Aedes aegypti samples selected to cover the temporal and geographic breadth of the Colombian outbreak. We performed phylogeographic analysis of these genomes, along with other publicly-available ZIKV genomes from the Americas, to estimate the frequency and timing of ZIKV introductions to Colombia. RESULTS: We attempted PCR amplification on 184 samples; 19 samples amplified sufficiently to perform sequencing. Of these, 8 samples yielded sequences with at least 50% coverage. Our phylogeographic reconstruction indicates two separate introductions of ZIKV to Colombia, one of which was previously unrecognized. We find that ZIKV was first introduced to Colombia in February 2015 (95%CI: Jan 2015 - Apr 2015), corresponding to 5 to 8 months of cryptic ZIKV transmission prior to confirmation in September 2015. Despite the presence of multiple introductions, we find that the majority of Colombian ZIKV diversity descends from a single introduction. We find evidence for movement of ZIKV from Colombia into bordering countries, including Peru, Ecuador, Panama, and Venezuela. CONCLUSIONS: Similarly to genomic epidemiological studies of ZIKV dynamics in other countries, we find that ZIKV circulated cryptically in Colombia. More accurately dating when ZIKV was circulating refines our definition of the population at risk. Additionally, our finding that the majority of ZIKV transmission within Colombia was attributable to transmission between individuals, rather than repeated travel-related importations, indicates that improved detection and control might have succeeded in limiting the scale of the outbreak within Colombia.

Contemporary and historical human migration patterns shape hepatitis B virus diversity.

Infection by hepatitis B virus (HBV) is responsible for approximately 296 million chronic cases of hepatitis B, and roughly 880,000 deaths annually. The global burden of HBV is distributed unevenly, largely owing to the heterogeneous geographic distribution of its subtypes, each of which demonstrates different severity and responsiveness to antiviral therapy. It is therefore crucial to the global public health response to HBV that the spatiotemporal spread of each genotype is well characterized. In this study, we describe a collection of 133 newly sequenced HBV strains from recent African immigrants upon their arrival in Belgium. We incorporate these sequences-all of which we determine to come from genotypes A, D, and E-into a large-scale phylogeographic study with genomes sampled across the globe. We focus on investigating the spatio-temporal processes shaping the evolutionary history of the three genotypes we observe. We incorporate several recently published ancient HBV genomes for genotypes A and D to aid our analysis. We show that different spatio-temporal processes underlie the A, D, and E genotypes with the former two having originated in southeastern Asia, after which they spread across the world. The HBV E genotype is estimated to have originated in Africa, after which it spread to Europe and the Americas. Our results highlight the use of phylogeographic reconstruction as a tool to understand the recent spatiotemporal dynamics of HBV, and highlight the importance of supporting vulnerable populations in accordance with the needs presented by specific HBV genotypes.

Two Separate Clusters of SARS-CoV-2 Delta Variant Infections in a Group of 41 Students Travelling from India: An Illustration of the Need for Rigorous Testing and Quarantine.

We report two clusters of SARS-CoV-2 B.1.617.2 (Delta variant) infections in a group of 41 Indian nursing students who travelled from New Delhi, India, to Belgium via Paris, France. All students tested negative before departure and had a second negative antigen test upon arrival in Paris. Upon arrival in Belgium, the students were quarantined in eight different houses. Four houses remained COVID-free during the 24 days of follow-up, while all 27 residents of the other four houses developed an infection during quarantine, including the four residents who were fully vaccinated and the two residents who were partially vaccinated. Genome sequencing revealed two distinct clusters affecting one and three houses, respectively. In this group of students, vaccination status did not seem to prevent infection nor decrease the viral load. No severe symptoms were reported. Extensive contact tracing and 3 months of nationwide genomic surveillance confirmed that these outbreaks were successfully contained and did not contribute to secondary community transmission in Belgium. These clusters highlight the importance of repeated testing and quarantine measures among travelers coming from countries experiencing a surge of infections, as all infections were detected 6 days or more after arrival.

Genome-wide diversity of Zika virus: Exploring spatio-temporal dynamics to guide a new nomenclature proposal.

The Zika virus (ZIKV) disease caused a public health emergency of international concern that started in February 2016. The overall number of ZIKV-related cases increased until November 2016, after which it declined sharply. While the evaluation of the potential risk and impact of future arbovirus epidemics remains challenging, intensified surveillance efforts along with a scale-up of ZIKV whole-genome sequencing provide an opportunity to understand the patterns of genetic diversity, evolution, and spread of ZIKV. However, a classification system that reflects the true extent of ZIKV genetic variation is lacking. Our objective was to characterize ZIKV genetic diversity and phylodynamics, identify genomic footprints of differentiation patterns, and propose a dynamic classification system that reflects its divergence levels. We analysed a curated dataset of 762 publicly available sequences spanning the full-length coding region of ZIKV from across its geographical span and collected between 1947 and 2021. The definition of genetic groups was based on comprehensive evolutionary dynamics analyses, which included recombination and phylogenetic analyses, within- and between-group pairwise genetic distances comparison, detection of selective pressure, and clustering analyses. Evidence for potential recombination events was detected in a few sequences. However, we argue that these events are likely due to sequencing errors as proposed in previous studies. There was evidence of strong purifying selection, widespread across the genome, as also detected for other arboviruses. A total of 50 sites showed evidence of positive selection, and for a few of these sites, there was amino acid (AA) differentiation between genetic clusters. Two main genetic clusters were defined, ZA and ZB, which correspond to the already characterized 'African' and 'Asian' genotypes, respectively. Within ZB, two subgroups, ZB.1 and ZB.2, represent the Asiatic and the American (and Oceania) lineages, respectively. ZB.1 is further subdivided into ZB.1.0 (a basal Malaysia sequence sampled in the 1960s and a recent Indian sequence), ZB.1.1 (South-Eastern Asia, Southern Asia, and Micronesia sequences), and ZB.1.2 (very similar sequences from the outbreak in Singapore). ZB.2 is subdivided into ZB.2.0 (basal American sequences and the sequences from French Polynesia, the putative origin of South America introduction), ZB.2.1 (Central America), and ZB.2.2 (Caribbean and North America). This classification system does not use geographical references and is flexible to accommodate potential future lineages. It will be a helpful tool for studies that involve analyses of ZIKV genomic variation and its association with pathogenicity and serve as a starting point for the public health surveillance and response to on-going and future epidemics and to outbreaks that lead to the emergence of new variants.

Two Years of Genomic Surveillance in Belgium during the SARS-CoV-2 Pandemic to Attain Country-Wide Coverage and Monitor the Introduction and Spread of Emerging Variants.

An adequate SARS-CoV-2 genomic surveillance strategy has proven to be essential for countries to obtain a thorough understanding of the variants and lineages being imported and successfully established within their borders. During 2020, genomic surveillance in Belgium was not structurally implemented but performed by individual research laboratories that had to acquire the necessary funds themselves to perform this important task. At the start of 2021, a nationwide genomic surveillance consortium was established in Belgium to markedly increase the country's genomic sequencing efforts (both in terms of intensity and representativeness), to perform quality control among participating laboratories, and to enable coordination and collaboration of research projects and publications. We here discuss the genomic surveillance efforts in Belgium before and after the establishment of its genomic sequencing consortium, provide an overview of the specifics of the consortium, and explore more details regarding the scientific studies that have been published as a result of the increased number of Belgian SARS-CoV-2 genomes that have become available.

Regional effect on the molecular clock rate of protein evolution in Eutherian and Metatherian genomes.

BACKGROUND: Different types of proteins diverge at vastly different rates. Moreover, the same type of protein has been observed to evolve with different rates in different phylogenetic lineages. In the present study we measured the rates of protein evolution in Eutheria (placental mammals) and Metatheria (marsupials) on a genome-wide basis and we propose that the gene position in the genome landscape has an important influence on the rate of protein divergence. RESULTS: We analyzed a protein-encoding gene set (n = 15,727) common to 16 mammals (12 Eutheria and 4 Metatheria). Using sliding windows that averaged regional effects of protein divergence we constructed landscapes in which strong and lineage-specific regional effects were seen on the molecular clock rate of protein divergence. Within each lineage, the relatively high rates were preferentially found in subtelomeric chromosomal regions. Such regions were observed to contain important and well-studied loci for fetal growth, uterine function and the generation of diversity in the adaptive repertoire of immunoglobulins. CONCLUSIONS: A genome landscape approach visualizes lineage-specific regional differences between Eutherian and Metatherian rates of protein evolution. This phenomenon of chromosomal position is a new element that explains at least part of the lineage-specific effects and differences between proteins on the molecular clock rates.

Exploiting genomic surveillance to map the spatio-temporal dispersal of SARS-CoV-2 spike mutations in Belgium across 2020.

At the end of 2020, several new variants of SARS-CoV-2-designated variants of concern-were detected and quickly suspected to be associated with a higher transmissibility and possible escape of vaccine-induced immunity. In Belgium, this discovery has motivated the initiation of a more ambitious genomic surveillance program, which is drastically increasing the number of SARS-CoV-2 genomes to analyse for monitoring the circulation of viral lineages and variants of concern. In order to efficiently analyse the massive collection of genomic data that are the result of such increased sequencing efforts, streamlined analytical strategies are crucial. In this study, we illustrate how to efficiently map the spatio-temporal dispersal of target mutations at a regional level. As a proof of concept, we focus on the Belgian province of Liège that has been consistently sampled throughout 2020, but was also one of the main epicenters of the second European epidemic wave. Specifically, we employ a recently developed phylogeographic workflow to infer the regional dispersal history of viral lineages associated with three specific mutations on the spike protein (S98F, A222V and S477N) and to quantify their relative importance through time. Our analytical pipeline enables analysing large data sets and has the potential to be quickly applied and updated to track target mutations in space and time throughout the course of an epidemic.

Variant Analysis of SARS-CoV-2 Genomes from Belgian Military Personnel Engaged in Overseas Missions and Operations.

More than a year after the first identification of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) as the causative agent of the 2019 coronavirus disease (COVID-19) in China, the emergence and spread of genomic variants of this virus through travel raise concerns regarding the introduction of lineages in previously unaffected regions, requiring adequate containment strategies. Concomitantly, such introductions fuel worries about a possible increase in transmissibility and disease severity, as well as a possible decrease in vaccine efficacy. Military personnel are frequently deployed on missions around the world. As part of a COVID-19 risk mitigation strategy, Belgian Armed Forces that engaged in missions and operations abroad were screened (7683 RT-qPCR tests), pre- and post-mission, for the presence of SARS-CoV-2, including the identification of viral lineages. Nine distinct viral genotypes were identified in soldiers returning from operations in Niger, the Democratic Republic of the Congo, Afghanistan, and Mali. The SARS-CoV-2 variants belonged to major clades 19B, 20A, and 20B (Nextstrain nomenclature), and included "variant of interest" B.1.525, "variant under monitoring" A.27, as well as lineages B.1.214, B.1, B.1.1.254, and A (pangolin nomenclature), some of which are internationally monitored due to the specific mutations they harbor. Through contact tracing and phylogenetic analysis, we show that isolation and testing policies implemented by the Belgian military command appear to have been successful in containing the influx and transmission of these distinct SARS-CoV-2 variants into military and civilian populations.

Emergence and spread of SARS-CoV-2 lineage B.1.620 with variant of concern-like mutations and deletions.

Distinct SARS-CoV-2 lineages, discovered through various genomic surveillance initiatives, have emerged during the pandemic following unprecedented reductions in worldwide human mobility. We here describe a SARS-CoV-2 lineage - designated B.1.620 - discovered in Lithuania and carrying many mutations and deletions in the spike protein shared with widespread variants of concern (VOCs), including E484K, S477N and deletions HV69Δ, Y144Δ, and LLA241/243Δ. As well as documenting the suite of mutations this lineage carries, we also describe its potential to be resistant to neutralising antibodies, accompanying travel histories for a subset of European cases, evidence of local B.1.620 transmission in Europe with a focus on Lithuania, and significance of its prevalence in Central Africa owing to recent genome sequencing efforts there. We make a case for its likely Central African origin using advanced phylogeographic inference methodologies incorporating recorded travel histories of infected travellers.

Evolution and rapid spread of a reassortant A(H3N2) virus that predominated the 2017-2018 influenza season.

The 2017-2018 North American influenza season caused more hospitalizations and deaths than any year since the 2009 H1N1 pandemic. The majority of recorded influenza infections were caused by A(H3N2) viruses, with most of the virus's North American diversity falling into the A2 clade. Within A2, we observe a subclade which we call A2/re that rose to comprise almost 70 per cent of A(H3N2) viruses circulating in North America by early 2018. Unlike most fast-growing clades, however, A2/re contains no amino acid substitutions in the hemagglutinin (HA) segment. Moreover, hemagglutination inhibition assays did not suggest substantial antigenic differences between A2/re viruses and viruses sampled during the 2016-2017 season. Rather, we observe that the A2/re clade was the result of a reassortment event that occurred in late 2016 or early 2017 and involved the combination of the HA and PB1 segments of an A2 virus with neuraminidase (NA) and other segments a virus from the clade A1b. The success of this clade shows the need for antigenic analysis that targets NA in addition to HA. Our results illustrate the potential for non-HA drivers of viral success and necessitate the need for more thorough tracking of full viral genomes to better understand the dynamics of influenza epidemics.