The 2022 multi-country monkeypox outbreak concurrent with the ongoing COVID-19 pandemic has further highlighted the need for genomic surveillance and pathogen whole genome sequencing. While metagenomic sequencing approaches have been used to sequence many of the early human monkeypox virus infections, these methods are resource intensive and require samples with high viral DNA concentrations. Given the atypical clinical presentation of cases associated with the current outbreak and uncertainty regarding viral load across both the course of infection and anatomical body sites, there is an urgent need for a more sensitive and broadly applicable sequencing approach. Amplicon-based sequencing (PrimalSeq) was initially developed for sequencing of Zika virus, and later adapted as the main sequencing approach for SARS-CoV-2. Here, we used PrimalScheme to develop a primer scheme for human monkeypox virus that can be used with many sequencing and bioinformatics pipelines implemented during the COVID-19 pandemic. We sequenced clinical samples that tested presumptive positive for monkeypox virus with amplicon-based and metagenomic sequencing approaches. Upon comparison, we found notably higher genome coverage across the virus genome, particularly in higher PCR cycle threshold (lower DNA titer) samples, with minimal amplicon drop-outs, in using the amplicon-based sequencing approach. By sending out primer pool aliquots to laboratories across the United States and internationally, we validated the primer scheme in 12 public health laboratories with their established Illumina or Oxford Nanopore Technologies sequencing workflows and with different sample types across a range of Ct values. Our findings suggest that amplicon-based sequencing increases the success rate across a wider range of viral DNA concentrations, with the PCR Ct value threshold at which laboratories were able to achieve 80% genome coverage at 10X ranging between Ct 25-33. Therefore, it increases the number of samples where near-complete genomes can be generated and it provides a cost-effective and widely applicable alternative to metagenomics for continued human monkeypox virus genomic surveillance. Importantly, we show that the human monkeypox virus primer scheme can be used with currently implemented amplicon-based SARS-CoV-2 sequencing workflows, with minimal change to the protocol.

Genomic surveillance of monkeypox virus (MPXV) during the 2022 outbreak has been mainly focused on single nucleotide polymorphism (SNP) changes. DNA viruses, including MPXV, have a lower SNP mutation rate than RNA viruses due to higher fidelity replication machinery. We identified a large genomic rearrangement in a MPXV sequence from a 2022 case in the state of Minnesota (MN), USA, from an abnormal, uneven MPXV read mapping coverage profile in whole-genome sequencing (WGS) data. We further screened WGS data of 206 U.S. MPXV samples and found seven (3.4 percent) sequenced genomes contained similar abnormal read coverage profiles that suggested putative large deletions or genomic rearrangements. Here, we present three MPXV genomes containing deletions ranging from 2.3 to 15 kb and four genomes containing more complex rearrangements. Five genomic changes were each only seen in one sample, but two sequences from linked cases shared an identical 2.3 kb deletion in the 3’ terminal region. All samples were positive using VAC1 and Clade II (formerly West African)-specific MPXV diagnostic tests; however, large deletions and genomic rearrangements like the ones reported here have the potential to result in viruses in which the target of a PCR diagnostic test is deleted. The emergence of genomic rearrangements during the outbreak may have public health implications and highlight the importance of continued genomic surveillance.

The 2022 multicountry mpox outbreak concurrent with the ongoing Coronavirus Disease 2019 (COVID-19) pandemic further highlighted the need for genomic surveillance and rapid pathogen whole-genome sequencing. While metagenomic sequencing approaches have been used to sequence many of the early mpox infections, these methods are resource intensive and require samples with high viral DNA concentrations. Given the atypical clinical presentation of cases associated with the outbreak and uncertainty regarding viral load across both the course of infection and anatomical body sites, there was an urgent need for a more sensitive and broadly applicable sequencing approach. Highly multiplexed amplicon-based sequencing (PrimalSeq) was initially developed for sequencing of Zika virus, and later adapted as the main sequencing approach for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). Here, we used PrimalScheme to develop a primer scheme for human monkeypox virus that can be used with many sequencing and bioinformatics pipelines implemented in public health laboratories during the COVID-19 pandemic. We sequenced clinical specimens that tested presumptively positive for human monkeypox virus with amplicon-based and metagenomic sequencing approaches. We found notably higher genome coverage across the virus genome, with minimal amplicon drop-outs, in using the amplicon-based sequencing approach, particularly in higher PCR cycle threshold (Ct) (lower DNA titer) samples. Further testing demonstrated that Ct value correlated with the number of sequencing reads and influenced the percent genome coverage. To maximize genome coverage when resources are limited, we recommend selecting samples with a PCR Ct below 31 Ct and generating 1 million sequencing reads per sample. To support national and international public health genomic surveillance efforts, we sent out primer pool aliquots to 10 laboratories across the United States, United Kingdom, Brazil, and Portugal. These public health laboratories successfully implemented the human monkeypox virus primer scheme in various amplicon sequencing workflows and with different sample types across a range of Ct values. Thus, we show that amplicon-based sequencing can provide a rapidly deployable, cost-effective, and flexible approach to pathogen whole-genome sequencing in response to newly emerging pathogens. Importantly, through the implementation of our primer scheme into existing SARS-CoV-2 workflows and across a range of sample types and sequencing platforms, we further demonstrate the potential of this approach for rapid outbreak response.

Lentinula is a broadly distributed group of fungi that contains the cultivated shiitake mushroom, L. edodes . We sequenced 24 genomes representing eight described species and several unnamed lineages of Lentinula from 15 countries on four continents. Lentinula comprises four major clades that arose in the Oligocene, three in the Americas and one in Asia–Australasia. To expand sampling of shiitake mushrooms, we assembled 60 genomes of L. edodes from China that were previously published as raw Illumina reads and added them to our dataset. Lentinula edodes sensu lato (s. lat.) contains three lineages that may warrant recognition as species, one including a single isolate from Nepal that is the sister group to the rest of L. edodes s. lat., a second with 20 cultivars and 12 wild isolates from China, Japan, Korea, and the Russian Far East, and a third with 28 wild isolates from China, Thailand, and Vietnam. Two additional lineages in China have arisen by hybridization among the second and third groups. Genes encoding cysteine sulfoxide lyase ( lecsl ) and γ-glutamyl transpeptidase ( leggt ), which are implicated in biosynthesis of the organosulfur flavor compound lenthionine, have diversified in Lentinula . Paralogs of both genes that are unique to Lentinula ( lecsl 3 and leggt 5b) are coordinately up-regulated in fruiting bodies of L. edodes . The pangenome of L. edodes s. lat. contains 20,308 groups of orthologous genes, but only 6,438 orthogroups (32%) are shared among all strains, whereas 3,444 orthogroups (17%) are found only in wild populations, which should be targeted for conservation.