Evaluation of variant calling algorithms for wastewater-based epidemiology using mixed populations of SARS-CoV-2 variants in synthetic and wastewater samples

Wastewater-based epidemiology has been used extensively throughout the COVID-19 (coronavirus disease 19) pandemic to detect and monitor the spread and prevalence of SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) and its variants. It has proven an excellent, complementary tool to clinical sequencing, supporting the insights gained and helping to make informed public-health decisions. Consequently, many groups globally have developed bioinformatics pipelines to analyse sequencing data from wastewater. Accurate calling of mutations is critical in this process and in the assignment of circulating variants; yet, to date, the performance of variant-calling algorithms in wastewater samples has not been investigated. To address this, we compared the performance of six variant callers (VarScan, iVar, GATK, FreeBayes, LoFreq and BCFtools), used widely in bioinformatics pipelines, on 19 synthetic samples with known ratios of three different SARS-CoV-2 variants of concern (VOCs) (Alpha, Beta and Delta), as well as 13 wastewater samples collected in London between the 15th and 18th December 2021. We used the fundamental parameters of recall (sensitivity) and precision (specificity) to confirm the presence of mutational profiles defining specific variants across the six variant callers. Our results show that BCFtools, FreeBayes and VarScan found the expected variants with higher precision and recall than GATK or iVar, although the latter identified more expected defining mutations than other callers. LoFreq gave the least reliable results due to the high number of false-positive mutations detected, resulting in lower precision. Similar results were obtained for both the synthetic and wastewater samples.

City-wide wastewater genomic surveillance through the successive emergence of SARS-CoV-2 Alpha and Delta variants

Genomic surveillance of SARS-CoV-2 has provided a critical evidence base for public health decisions throughout the pandemic. Sequencing data from clinical cases has helped to understand disease transmission and the spread of novel variants. Genomic wastewater surveillance can offer important, complementary information by providing frequency estimates of all variants circulating in a population without sampling biases. Here we show that genomic SARS-CoV-2 wastewater surveillance can detect fine-scale differences within urban centres, specifically within the city of Liverpool, UK, during the emergence of Alpha and Delta variants between November 2020 and June 2021. Furthermore, wastewater and clinical sequencing match well in the estimated timing of new variant rises and the first detection of a new variant in a given area may occur in either clinical or wastewater samples. The study's main limitation was sample quality when infection prevalence was low in spring 2021, resulting in a lower resolution of the rise of the Delta variant compared to the rise of the Alpha variant in the previous winter. The correspondence between wastewater and clinical variant frequencies demonstrates the reliability of wastewater surveillance. However, discrepancies in the first detection of the Alpha variant between the two approaches highlight that wastewater monitoring can also capture missing information, possibly resulting from asymptomatic cases or communities less engaged with testing programmes, as found by a simultaneous surge testing effort across the city.

Microbiome function predicts amphibian chytridiomycosis disease dynamics

BACKGROUND

The fungal pathogen Batrachochytrium dendrobatidis (Bd) threatens amphibian biodiversity and ecosystem stability worldwide. Amphibian skin microbial community structure has been linked to the clinical outcome of Bd infections, yet its overall functional importance is poorly understood.

METHODS

Microbiome taxonomic and functional profiles were assessed using high-throughput bacterial 16S rRNA and fungal ITS2 gene sequencing, bacterial shotgun metagenomics and skin mucosal metabolomics. We sampled 56 wild midwife toads (Alytes obstetricans) from montane populations exhibiting Bd epizootic or enzootic disease dynamics. In addition, to assess whether disease-specific microbiome profiles were linked to microbe-mediated protection or Bd-induced perturbation, we performed a laboratory Bd challenge experiment whereby 40 young adult A. obstetricans were exposed to Bd or a control sham infection. We measured temporal changes in the microbiome as well as functional profiles of Bd-exposed and control animals at peak infection.

RESULTS

Microbiome community structure and function differed in wild populations based on infection history and in experimental control versus Bd-exposed animals. Bd exposure in the laboratory resulted in dynamic changes in microbiome community structure and functional differences, with infection clearance in all but one infected animal. Sphingobacterium, Stenotrophomonas and an unclassified Commamonadaceae were associated with wild epizootic dynamics and also had reduced abundance in laboratory Bd-exposed animals that cleared infection, indicating a negative association with Bd resistance. This was further supported by microbe-metabolite integration which identified functionally relevant taxa driving disease outcome, of which Sphingobacterium and Bd were most influential in wild epizootic dynamics. The strong correlation between microbial taxonomic community composition and skin metabolome in the laboratory and field is inconsistent with microbial functional redundancy, indicating that differences in microbial taxonomy drive functional variation. Shotgun metagenomic analyses support these findings, with similar disease-associated patterns in beta diversity. Analysis of differentially abundant bacterial genes and pathways indicated that bacterial environmental sensing and Bd resource competition are likely to be important in driving infection outcomes.

CONCLUSIONS

Bd infection drives altered microbiome taxonomic and functional profiles across laboratory and field environments. Our application of multi-omics analyses in experimental and field settings robustly predicts Bd disease dynamics and identifies novel candidate biomarkers of infection. Video Abstract.

Generation and transmission of interlineage recombinants in the SARS-CoV-2 pandemic

We present evidence for multiple independent origins of recombinant SARS-CoV-2 viruses sampled from late 2020 and early 2021 in the United Kingdom. Their genomes carry single-nucleotide polymorphisms and deletions that are characteristic of the B.1.1.7 variant of concern but lack the full complement of lineage-defining mutations. Instead, the remainder of their genomes share contiguous genetic variation with non-B.1.1.7 viruses circulating in the same geographic area at the same time as the recombinants. In four instances, there was evidence for onward transmission of a recombinant-origin virus, including one transmission cluster of 45 sequenced cases over the course of 2 months. The inferred genomic locations of recombination breakpoints suggest that every community-transmitted recombinant virus inherited its spike region from a B.1.1.7 parental virus, consistent with a transmission advantage for B.1.1.7's set of mutations.

Post-epizootic microbiome associations across communities of neotropical amphibians

Microbiome-pathogen interactions are increasingly recognized as an important element of host immunity. While these host-level interactions will have consequences for community disease dynamics, the factors which influence host microbiomes at larger scales are poorly understood. We here describe landscape-scale pathogen-microbiome associations within the context of post-epizootic amphibian chytridiomycosis, a disease caused by the panzootic chytrid fungus Batrachochytrium dendrobatidis. We undertook a survey of Neotropical amphibians across altitudinal gradients in Ecuador ~30 years following the observed amphibian declines and collected skin swab-samples which were metabarcoded using both fungal (ITS-2) and bacterial (r16S) amplicons. The data revealed marked variation in patterns of both B. dendrobatidis infection and microbiome structure that are associated with host life history. Stream breeding amphibians were most likely to be infected with B. dendrobatidis. This increased probability of infection was further associated with increased abundance and diversity of non-Batrachochytrium chytrid fungi in the skin and environmental microbiome. We also show that increased alpha diversity and the relative abundance of fungi are lower in the skin microbiome of adult stream amphibians compared to adult pond-breeding amphibians, an association not seen for bacteria. Finally, stream tadpoles exhibit lower proportions of predicted protective microbial taxa than pond tadpoles, suggesting reduced biotic resistance. Our analyses show that host breeding ecology strongly shapes pathogen-microbiome associations at a landscape scale, a trait that may influence resilience in the face of emerging infectious diseases.

Compensatory growth renders Tcf7l1a dispensable for eye formation despite its requirement in eye field specification

The vertebrate eye originates from the eye field, a domain of cells specified by a small number of transcription factors. In this study, we show that Tcf7l1a is one such transcription factor that acts cell-autonomously to specify the eye field in zebrafish. Despite the much-reduced eye field in tcf7l1a mutants, these fish develop normal eyes revealing a striking ability of the eye to recover from a severe early phenotype. This robustness is not mediated through genetic compensation at neural plate stage; instead, the smaller optic vesicle of tcf7l1a mutants shows delayed neurogenesis and continues to grow until it achieves approximately normal size. Although the developing eye is robust to the lack of Tcf7l1a function, it is sensitised to the effects of additional mutations. In support of this, a forward genetic screen identified mutations in hesx1, cct5 and gdf6a, which give synthetically enhanced eye specification or growth phenotypes when in combination with the tcf7l1a mutation.

Left and right eyes develop independently, yet they consistently grow to roughly the same size in humans and other creatures. How they do this remains a mystery, though scientists have learned that both eyes originate from a single group of cells in the developing nervous system called the eye field. As development progresses, the eye field splits in two, and buds into the two separate compartments from which each eye forms. As the eyes grow, the cells in each compartment specialize, or ‘differentiate’, to make working left and right eyes. Scientists often study eye development in zebrafish embryos because it is easy to see each step in the process.

Now, Young at al. show that zebrafish with a mutation that causes the eye field to be half its normal size go on to form normal-sized eyes. Somehow these developing embryos overcome this deleterious mutation. It turns out that the eyes of zebrafish with this mutation grow for a longer period of time than typical zebrafish eyes. This change allows the mutant fish’s eyes to catch up and reach normal size. When Young et al. removed some cells from one of the forming eyes of normal zebrafish embryos they found that same thing happened. The smaller eye developed for a longer time and delayed its differentiation until both eyes were the same size. Conversely, when eyes developed from a larger than normal eye field, growth stopped prematurely and differentiation began early preventing the eyes from ending up oversized.

Though the fish were able to overcome the effects of one mutation to develop normal-sized eyes, adding a second mutation that affected eye development led to unusual sized eyes or absence of eyes. Together the experiments identify genes and mechanisms essential for the formation and size of the eyes. Given that the processes underlying eye formation are very similar in many animals, this new information should help scientists to better understand eye abnormalities in humans.

Development and worldwide use of non-lethal, and minimal population-level impact, protocols for the isolation of amphibian chytrid fungi

Parasitic chytrid fungi have emerged as a significant threat to amphibian species worldwide, necessitating the development of techniques to isolate these pathogens into culture for research purposes. However, early methods of isolating chytrids from their hosts relied on killing amphibians. We modified a pre-existing protocol for isolating chytrids from infected animals to use toe clips and biopsies from toe webbing rather than euthanizing hosts, and distributed the protocol to researchers as part of the BiodivERsA project RACE; here called the RML protocol. In tandem, we developed a lethal procedure for isolating chytrids from tadpole mouthparts. Reviewing a database of use a decade after their inception, we find that these methods have been applied across 5 continents, 23 countries and in 62 amphibian species. Isolation of chytrids by the non-lethal RML protocol occured in 18% of attempts with 207 fungal isolates and three species of chytrid being recovered. Isolation of chytrids from tadpoles occured in 43% of attempts with 334 fungal isolates of one species (Batrachochytrium dendrobatidis) being recovered. Together, these methods have resulted in a significant reduction and refinement of our use of threatened amphibian species and have improved our ability to work with this group of emerging pathogens.

Recent Asian origin of chytrid fungi causing global amphibian declines

Globalized infectious diseases are causing species declines worldwide, but their source often remains elusive. We used whole-genome sequencing to solve the spatiotemporal origins of the most devastating panzootic to date, caused by the fungus Batrachochytrium dendrobatidis, a proximate driver of global amphibian declines. We traced the source of B. dendrobatidis to the Korean peninsula, where one lineage, BdASIA-1, exhibits the genetic hallmarks of an ancestral population that seeded the panzootic. We date the emergence of this pathogen to the early 20th century, coinciding with the global expansion of commercial trade in amphibians, and we show that intercontinental transmission is ongoing. Our findings point to East Asia as a geographic hotspot for B. dendrobatidis biodiversity and the original source of these lineages that now parasitize amphibians worldwide.

Antagonism between Gdf6a and retinoic acid pathways controls timing of retinal neurogenesis and growth of the eye in zebrafish

Maintaining neurogenesis in growing tissues requires a tight balance between progenitor cell proliferation and differentiation. In the zebrafish retina, neuronal differentiation proceeds in two stages with embryonic retinal progenitor cells (RPCs) of the central retina accounting for the first rounds of differentiation, and stem cells from the ciliary marginal zone (CMZ) being responsible for late neurogenesis and growth of the eye. In this study, we analyse two mutants with small eyes that display defects during both early and late phases of retinal neurogenesis. These mutants carry lesions in gdf6a, a gene encoding a BMP family member previously implicated in dorsoventral patterning of the eye. We show that gdf6a mutant eyes exhibit expanded retinoic acid (RA) signalling and demonstrate that exogenous activation of this pathway in wild-type eyes inhibits retinal growth, generating small eyes with a reduced CMZ and fewer proliferating progenitors, similar to gdf6a mutants. We provide evidence that RA regulates the timing of RPC differentiation by promoting cell cycle exit. Furthermore, reducing RA signalling in gdf6a mutants re-establishes appropriate timing of embryonic retinal neurogenesis and restores putative stem and progenitor cell populations in the CMZ. Together, our results support a model in which dorsally expressed gdf6a limits RA pathway activity to control the transition from proliferation to differentiation in the growing eye.

Summary: In the vertebrate eye, dorsally expressed Gdf6a limits RA pathway activity to control the transition from proliferation to differentiation, thereby regulating eye size.

In vivo Wnt signaling tracing through a transgenic biosensor fish reveals novel activity domains

The creation of molecular tools able to unravel in vivo spatiotemporal activation of specific cell signaling events during cell migration, differentiation and morphogenesis is of great relevance to developmental cell biology. Here, we describe the generation, validation and applications of two transgenic reporter lines for Wnt/β-catenin signaling, named TCFsiam, and show that they are reliable and sensitive Wnt biosensors for in vivo studies. We demonstrate that these lines sensitively detect Wnt/β-catenin pathway activity in several cellular contexts, from sensory organs to cardiac valve patterning. We provide evidence that Wnt/β-catenin activity is involved in the formation and maintenance of the zebrafish CNS blood vessel network, on which sox10 neural crest-derived cells migrate and proliferate. We finally show that these transgenic lines allow for screening of Wnt signaling modifying compounds, tissue regeneration assessment as well as evaluation of potential Wnt/β-catenin genetic modulators.