Swiss public health measures associated with reduced SARS-CoV-2 transmission using genome data.

Genome sequences from evolving infectious pathogens allow quantification of case introductions and local transmission dynamics. We sequenced 11,357 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes from Switzerland in 2020-the sixth largest effort globally. Using a representative subset of these data, we estimated viral introductions to Switzerland and their persistence over the course of 2020. We contrasted these estimates with simple null models representing the absence of certain public health measures. We show that Switzerland's border closures decoupled case introductions from incidence in neighboring countries. Under a simple model, we estimate an 86 to 98% reduction in introductions during Switzerland's strictest border closures. Furthermore, the Swiss 2020 partial lockdown roughly halved the time for sampled introductions to die out. Last, we quantified local transmission dynamics once introductions into Switzerland occurred using a phylodynamic model. We found that transmission slowed 35 to 63% upon outbreak detection in summer 2020 but not in fall. This finding may indicate successful contact tracing over summer before overburdening in fall. The study highlights the added value of genome sequencing data for understanding transmission dynamics.

The Leaf Microbiome of Arabidopsis Displays Reproducible Dynamics and Patterns throughout the Growing Season.

Leaves are primarily responsible for the plant's photosynthetic activity. Thus, changes in the leaf microbiota, which includes deleterious and beneficial microbes, can have far-reaching effects on plant fitness and productivity. Identifying the processes and microorganisms that drive these changes over a plant's lifetime is, therefore, crucial. In this study, we analyzed the temporal dynamics in the leaf microbiome of Arabidopsis thaliana, integrating changes in both composition and microbe-microbe interactions via the study of microbial networks. Field-grown Arabidopsis were used to monitor leaf bacterial, fungal and oomycete communities throughout the plant's natural growing season (extending from November to March) over three consecutive years. Our results revealed the existence of conserved temporal patterns, with microbial communities and networks going through a stabilization phase of decreased diversity and variability at the beginning of the plant's growing season. Despite a high turnover in these communities, we identified 19 "core" taxa persisting on Arabidopsis leaves across time and plant generations. With the hypothesis these microbes could be playing key roles in the structuring of leaf microbial communities, we conducted a time-informed microbial network analysis which showed core taxa are not necessarily highly connected network "hubs," and "hubs" alternate with time. Our study shows that leaf microbial communities exhibit reproducible dynamics and patterns, suggesting the potential of using our understanding of temporal trajectories in microbial community composition to design experiments aimed at driving these communities toward desired states. IMPORTANCE Utilizing plant microbiota to promote plant growth and plant health is key to more environmentally friendly agriculture. A major bottleneck in the engineering of plant-beneficial microbial communities is the low persistence of applied microbes under filed conditions, especially considering plant leaves. Indeed, although many leaf-associated microorganisms have the potential to promote plant growth and protect plants from pathogens, few of them are able to survive and thrive over time. In our study, we could show that leaf microbial communities are very variable at the beginning of the plant growing season but become more and more similar and less variable as the season progresses. We further identify a cohort of 19 "core" microbes, systematically present on plant leaves that would make these microbes exceptional candidates for future agricultural applications.

Determinants of SARS-CoV-2 transmission to guide vaccination strategy in an urban area.

Transmission chains within small urban areas (accommodating ∼30 per cent of the European population) greatly contribute to case burden and economic impact during the ongoing coronavirus pandemic and should be a focus for preventive measures to achieve containment. Here, at very high spatio-temporal resolution, we analysed determinants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission in a European urban area, Basel-City (Switzerland). We combined detailed epidemiological, intra-city mobility and socio-economic data sets with whole-genome sequencing during the first SARS-CoV-2 wave. For this, we succeeded in sequencing 44 per cent of all reported cases from Basel-City and performed phylogenetic clustering and compartmental modelling based on the dominating viral variant (B.1-C15324T; 60 per cent of cases) to identify drivers and patterns of transmission. Based on these results we simulated vaccination scenarios and corresponding healthcare system burden (intensive care unit (ICU) occupancy). Transmissions were driven by socio-economically weaker and highly mobile population groups with mostly cryptic transmissions which lacked genetic and identifiable epidemiological links. Amongst more senior population transmission was clustered. Simulated vaccination scenarios assuming 60-90 per cent transmission reduction and 70-90 per cent reduction of severe cases showed that prioritising mobile, socio-economically weaker populations for vaccination would effectively reduce case numbers. However, long-term ICU occupation would also be effectively reduced if senior population groups were prioritised, provided there were no changes in testing and prevention strategies. Reducing SARS-CoV-2 transmission through vaccination strongly depends on the efficacy of the deployed vaccine. A combined strategy of protecting risk groups by extensive testing coupled with vaccination of the drivers of transmission (i.e. highly mobile groups) would be most effective at reducing the spread of SARS-CoV-2 within an urban area.

External Quality Assessment of SARS-CoV-2 Sequencing: an ESGMD-SSM Pilot Trial across 15 European Laboratories.

This first pilot trial on external quality assessment (EQA) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) whole-genome sequencing, initiated by the European Society of Clinical Microbiology and Infectious Diseases (ESCMID) Study Group for Genomic and Molecular Diagnostics (ESGMD) and the Swiss Society for Microbiology (SSM), aims to build a framework between laboratories in order to improve pathogen surveillance sequencing. Ten samples with various viral loads were sent out to 15 clinical laboratories that had free choice of sequencing methods and bioinformatic analyses. The key aspects on which the individual centers were compared were the identification of (i) single nucleotide polymorphisms (SNPs) and indels, (ii) Pango lineages, and (iii) clusters between samples. The participating laboratories used a wide array of methods and analysis pipelines. Most were able to generate whole genomes for all samples. Genomes were sequenced to various depths (up to a 100-fold difference across centers). There was a very good consensus regarding the majority of reporting criteria, but there were a few discrepancies in lineage and cluster assignments. Additionally, there were inconsistencies in variant calling. The main reasons for discrepancies were missing data, bioinformatic choices, and interpretation of data. The pilot EQA was overall a success. It was able to show the high quality of participating laboratories and provide valuable feedback in cases where problems occurred, thereby improving the sequencing setup of laboratories. A larger follow-up EQA should, however, improve on defining the variables and format of the report. Additionally, contamination and/or minority variants should be a further aspect of assessment.

Global Genomic Analysis of SARS-CoV-2 RNA Dependent RNA Polymerase Evolution and Antiviral Drug Resistance.

A variety of antiviral treatments for COVID-19 have been investigated, involving many repurposed drugs. Currently, the SARS-CoV-2 RNA-dependent RNA polymerase (RdRp, encoded by nsp12-nsp7-nsp8) has been targeted by numerous inhibitors, e.g., remdesivir, the only provisionally approved treatment to-date, although the clinical impact of these interventions remains inconclusive. However, the potential emergence of antiviral resistance poses a threat to the efficacy of any successful therapies on a wide scale. Here, we propose a framework to monitor the emergence of antiviral resistance, and as a proof of concept, we address the interaction between RdRp and remdesivir. We show that SARS-CoV-2 RdRp is under purifying selection, that potential escape mutations are rare in circulating lineages, and that those mutations, where present, do not destabilise RdRp. In more than 56,000 viral genomes from 105 countries from the first pandemic wave, we found negative selective pressure affecting nsp12 (Tajima's D = -2.62), with potential antiviral escape mutations in only 0.3% of sequenced genomes. Potential escape mutations included known key residues, such as Nsp12:Val473 and Nsp12:Arg555. Of the potential escape mutations involved globally, in silico structural models found that they were unlikely to be associated with loss of stability in RdRp. No potential escape mutation was found in a local cohort of remdesivir treated patients. Collectively, these findings indicate that RdRp is a suitable drug target, and that remdesivir does not seem to exert high selective pressure. We anticipate our framework to be the starting point of a larger effort for a global monitoring of drug resistance throughout the COVID-19 pandemic.

Obtaining deeper insights into microbiome diversity using a simple method to block host and nontargets in amplicon sequencing.

Profiling diverse microbiomes is revolutionizing our understanding of biological mechanisms and ecologically relevant problems, including metaorganism (host + microbiome) assembly, functions and adaptation. Amplicon sequencing of multiple conserved, phylogenetically informative loci has therefore become an instrumental tool for many researchers. Investigations in many systems are hindered, however, since essential sequencing depth can be lost by amplification of nontarget DNA from hosts or overabundant microorganisms. Here, we introduce "blocking oligos", a low-cost and flexible method using standard oligonucleotides to block amplification of diverse nontargets and software to aid their design. We apply them primarily in leaves, where exceptional challenges with host amplification prevail. A. thaliana-specific blocking oligos applied in eight different target loci reduce undesirable host amplification by up to 90%. To expand applicability, we designed universal 16S and 18S rRNA gene plant blocking oligos for targets that are conserved in diverse plant species and demonstrate that they efficiently block five plant species from five orders spanning monocots and dicots (Bromus erectus, Plantago lanceolata, Lotus corniculatus, Amaranth sp., Arabidopsis thaliana). These can increase alpha diversity discovery without biasing beta diversity patterns and do not compromise microbial load information inherent to plant-derived 16S rRNA gene amplicon sequencing data. Finally, we designed and tested blocking oligos to avoid amplification of 18S rRNA genes of a sporulating oomycete pathogen, demonstrating their effectiveness in applications well beyond plants. Using these tools, we generated a survey of the A. thaliana leaf microbiome based on eight loci targeting bacterial, fungal, oomycete and other eukaryotic microorganisms and discuss complementarity of commonly used amplicon sequencing regions for describing leaf microbiota. This approach has potential to make questions in a variety of study systems more tractable by making amplicon sequencing more targeted, leading to deeper, systems-based insights into microbial discovery. For fast and easy design for blocking oligos for any nontarget DNA in other study systems, we developed a publicly available R package.

SARS-CoV-2 N501Y Introductions and Transmissions in Switzerland from Beginning of October 2020 to February 2021-Implementation of Swiss-Wide Diagnostic Screening and Whole Genome Sequencing.

The rapid spread of the SARS-CoV-2 lineages B.1.1.7 (N501Y.V1) throughout the UK, B.1.351 (N501Y.V2) in South Africa, and P.1 (B.1.1.28.1; N501Y.V3) in Brazil has led to the definition of variants of concern (VoCs) and recommendations for lineage specific surveillance. In Switzerland, during the last weeks of December 2020, we established a nationwide screening protocol across multiple laboratories, focusing first on epidemiological and microbiological definitions. In January 2021, we validated and implemented an N501Y-specific PCR to rapidly screen for VoCs, which are then confirmed using amplicon sequencing or whole genome sequencing (WGS). A total of 13,387 VoCs have been identified since the detection of the first Swiss case in October 2020, with 4194 being B.1.1.7, 172 B.1.351, and 7 P.1. The remaining 9014 cases of VoCs have been described without further lineage specification. Overall, all diagnostic centers reported a rapid increase of the percentage of detected VOCs, with a range of 6 to 46% between 25 to 31 of January 2021 increasing towards 41 to 82% between 22 to 28 of February. A total of 739 N501Y positive genomes were analysed and show a broad range of introduction events to Switzerland. In this paper, we describe the nationwide coordination and implementation process across laboratories, public health institutions, and researchers, the first results of our N501Y-specific variant screening, and the phylogenetic analysis of all available WGS data in Switzerland, that together identified the early introduction events and subsequent community spreading of the VoCs.

SARS-CoV-2 outbreak in a tri-national urban area is dominated by a B.1 lineage variant linked to a mass gathering event.

The first case of SARS-CoV-2 in Basel, Switzerland was detected on February 26th 2020. We present a phylogenetic study to explore viral introduction and evolution during the exponential early phase of the local COVID-19 outbreak from February 26th until March 23rd. We sequenced SARS-CoV-2 naso-oropharyngeal swabs from 746 positive tests that were performed at the University Hospital Basel during the study period. We successfully generated 468 high quality genomes from unique patients and called variants with our COVID-19 Pipeline (COVGAP), and analysed viral genetic diversity using PANGOLIN taxonomic lineages. To identify introduction and dissemination events we incorporated global SARS-CoV-2 genomes and inferred a time-calibrated phylogeny. Epidemiological data from patient questionnaires was used to facilitate the interpretation of phylogenetic observations. The early outbreak in Basel was dominated by lineage B.1 (83·6%), detected first on March 2nd, although the first sample identified belonged to B.1.1. Within B.1, 68·2% of our samples fall within a clade defined by the SNP C15324T ('Basel cluster'), including 157 identical sequences at the root of the 'Basel cluster', some of which we can specifically trace to regional spreading events. We infer the origin of B.1-C15324T to mid-February in our tri-national region. The other genomes map broadly over the global phylogenetic tree, showing several introduction events from and/or dissemination to other regions of the world via travellers. Family transmissions can also be traced in our data. A single lineage variant dominated the outbreak in the Basel area while other lineages, such as the first (B.1.1), did not propagate. A mass gathering event was the predominant initial source of cases, with travel returners and family transmissions to a lesser extent. We highlight the importance of adding specific questions to epidemiological questionnaires, to obtain data on attendance of large gatherings and their locations, as well as travel history, to effectively identify routes of transmissions in up-coming outbreaks. This phylogenetic analysis in concert with epidemiological and contact tracing data, allows connection and interpretation of events, and can inform public health interventions. Trial Registration: ClinicalTrials.gov NCT04351503.

Brief validation of the novel GeneXpert Xpress SARS-CoV-2 PCR assay.

The clinical and epidemiologic management of the SARS-CoV-2 pandemic is critically dependent on molecular assays with short turn-around time. We validated the novel Xpert Xpress SARS-CoV-2 assay using a commercial nucleic acid testing (Roche Cobas 6800). We found an excellent concordance over a range of SARS-CoV-2 loads and across established human coronaviruses.

Early Lotus japonicus root transcriptomic responses to symbiotic and pathogenic fungal exudates.

The objective of this study is to evaluate Lotus japonicus transcriptomic responses to arbuscular mycorrhizal (AM) germinated spore exudates (GSEs), responsible for activating nuclear Ca(2+) spiking in plant root epidermis. A microarray experiment was performed comparing gene expression in Lotus rootlets treated with GSE or water after 24 and 48 h. The transcriptional pattern of selected genes that resulted to be regulated in the array was further evaluated upon different treatments and timings. In particular, Lotus rootlets were treated with: GSE from the pathogenic fungus Colletotrichum trifolii; short chitin oligomers (COs; acknowledged AM fungal signals) and long COs (as activators of pathogenic responses). This experimental set up has revealed that AM GSE generates a strong transcriptomic response in Lotus roots with an extensive defense-related response after 24 h and a subsequent down-regulation after 48 h. A similar subset of defense-related genes resulted to be up-regulated also upon treatment with C. trifolii GSE, although with an opposite trend. Surprisingly, long COs activated both defense-like and symbiosis-related genes. Among the genes regulated in the microarray, promoter-GUS assay showed that LjMATE1 activates in epidermal cells and root hairs.