Complete de novo assembly of Wolbachia endosymbiont of Drosophila willistoni using long-read genome sequencing.

Wolbachia is an obligate intracellular 𝛼-proteobacterium which commonly infects arthropods and filarial nematodes. Different strains of Wolbachia are capable of a wide range of regulatory manipulations in many hosts and modulate host cellular differentiation to influence host reproduction. The genetic basis for the majority of these phenotypes is unknown. The w Wil strain from the neotropical fruit fly, Drosophila willistoni , exhibits a remarkably high affinity for host germline-derived cells relative to the soma. This trait could be leveraged for understanding how Wolbachia influences the host germline and for controlling host populations in the field. To further the use of this strain in biological and biomedical research, we sequenced the genome of the w Wil strain isolated from host cell culture cells. Here, we present the first high quality nanopore assembly of w Wil, the Wolbachia endosymbiont of D. willistoni . Our assembly resulted in a circular genome of 1.27 Mb with a BUSCO completeness score of 99.7%. Consistent with other insect-associated Wolbachia strains, comparative genomic analysis revealed that wWil has a highly mosaic genome relative to the closely related wMel strain from Drosophila melanogaster .

Addressing pandemic-wide systematic errors in the SARS-CoV-2 phylogeny.

The SARS-CoV-2 genome occupies a unique place in infection biology - it is the most highly sequenced genome on earth (making up over 20% of public sequencing datasets) with fine scale information on sampling date and geography, and has been subject to unprecedented intense analysis. As a result, these phylogenetic data are an incredibly valuable resource for science and public health. However, the vast majority of the data was sequenced by tiling amplicons across the full genome, with amplicon schemes that changed over the pandemic as mutations in the viral genome interacted with primer binding sites. In combination with the disparate set of genome assembly workflows and lack of consistent quality control (QC) processes, the current genomes have many systematic errors that have evolved with the virus and amplicon schemes. These errors have significant impacts on the phylogeny, and therefore over the last few years, many thousands of hours of researchers time has been spent in "eyeballing" trees, looking for artefacts, and then patching the tree. Given the huge value of this dataset, we therefore set out to reprocess the complete set of public raw sequence data in a rigorous amplicon-aware manner, and build a cleaner phylogeny. Here we provide a global tree of 3,960,704 samples, built from a consistently assembled set of high quality consensus sequences from all available public data as of March 2023, viewable at https://viridian.taxonium.org. Each genome was constructed using a novel assembly tool called Viridian (https://github.com/iqbal-lab-org/viridian), developed specifically to process amplicon sequence data, eliminating artefactual errors and mask the genome at low quality positions. We provide simulation and empirical validation of the methodology, and quantify the improvement in the phylogeny. Phase 2 of our project will address the fact that the data in the public archives is heavily geographically biased towards the Global North. We therefore have contributed new raw data to ENA/SRA from many countries including Ghana, Thailand, Laos, Sri Lanka, India, Argentina and Singapore. We will incorporate these, along with all public raw data submitted between March 2023 and the current day, into an updated set of assemblies, and phylogeny. We hope the tree, consensus sequences and Viridian will be a valuable resource for researchers.