Spruce giga-genomes: structurally similar yet distinctive with differentially expanding gene families and rapidly evolving genes.

Spruces (Picea spp.) are coniferous trees widespread in boreal and mountainous forests of the northern hemisphere, with large economic significance and enormous contributions to global carbon sequestration. Spruces harbor very large genomes with high repetitiveness, hampering their comparative analysis. Here, we present and compare the genomes of four different North American spruces: the genome assemblies for Engelmann spruce (Picea engelmannii) and Sitka spruce (Picea sitchensis) together with improved and more contiguous genome assemblies for white spruce (Picea glauca) and for a naturally occurring introgress of these three species known as interior spruce (P. engelmannii × glauca × sitchensis). The genomes were structurally similar, and a large part of scaffolds could be anchored to a genetic map. The composition of the interior spruce genome indicated asymmetric contributions from the three ancestral genomes. Phylogenetic analysis of the nuclear and organelle genomes revealed a topology indicative of ancient reticulation. Different patterns of expansion of gene families among genomes were observed and related with presumed diversifying ecological adaptations. We identified rapidly evolving genes that harbored high rates of non-synonymous polymorphisms relative to synonymous ones, indicative of positive selection and its hitchhiking effects. These gene sets were mostly distinct between the genomes of ecologically contrasted species, and signatures of convergent balancing selection were detected. Stress and stimulus response was identified as the most frequent function assigned to expanding gene families and rapidly evolving genes. These two aspects of genomic evolution were complementary in their contribution to divergent evolution of presumed adaptive nature. These more contiguous spruce giga-genome sequences should strengthen our understanding of conifer genome structure and evolution, as their comparison offers clues into the genetic basis of adaptation and ecology of conifers at the genomic level. They will also provide tools to better monitor natural genetic diversity and improve the management of conifer forests. The genomes of four closely related North American spruces indicate that their high similarity at the morphological level is paralleled by the high conservation of their physical genome structure. Yet, the evidence of divergent evolution is apparent in their rapidly evolving genomes, supported by differential expansion of key gene families and large sets of genes under positive selection, largely in relation to stimulus and environmental stress response.

BioThings SDK: a toolkit for building high-performance data APIs in biomedical research.

SUMMARY: To meet the increased need of making biomedical resources more accessible and reusable, Web Application Programming Interfaces (APIs) or web services have become a common way to disseminate knowledge sources. The BioThings APIs are a collection of high-performance, scalable, annotation as a service APIs that automate the integration of biological annotations from disparate data sources. This collection of APIs currently includes MyGene.info, MyVariant.info and MyChem.info for integrating annotations on genes, variants and chemical compounds, respectively. These APIs are used by both individual researchers and application developers to simplify the process of annotation retrieval and identifier mapping. Here, we describe the BioThings Software Development Kit (SDK), a generalizable and reusable toolkit for integrating data from multiple disparate data sources and creating high-performance APIs. This toolkit allows users to easily create their own BioThings APIs for any data type of interest to them, as well as keep APIs up-to-date with their underlying data sources. AVAILABILITY AND IMPLEMENTATION: The BioThings SDK is built in Python and released via PyPI (https://pypi.org/project/biothings/). Its source code is hosted at its github repository (https://github.com/biothings/biothings.api). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Comparison of Video Versus Direct Laryngoscopy: A Prospective Prehospital Air Medical Services Study.

OBJECTIVE: Videolaryngoscopy (VL) in the prehospital setting remains controversial, with conflicting data on its utility. We compared C-MAC VL (Karl Storz, Tuttlingen, Germany) versus direct laryngoscopy (DL) in the prehospital setting, recording the grade of the glottic view, first pass success (FPS), overall success, and equipment functionality. METHODS: We conducted a prospective observational study with a convenience sample of 49 adult patients who were intubated by flight crew nurses and paramedics using the C-MAC videolaryngoscope from April to November 2013. We compared Cormack-Lehane (CL) grades of view for DL and VL, intubation success rates, and equipment functionality. RESULTS: CL grades 1 or 2 were obtained with 24 patients (49%) with DL and 45 patients (92%) with VL. Of the 25 patients (51%) who had a CL grade 3 or 4 view on DL, 22 of those patients (88%) converted to a CL grade 1 or 2 with VL (P < .001). There was an overall success rate of 96% and an FPS rate of 71%. The C-MAC videolaryngoscope was functional during intubation 100% of the time. CONCLUSION: VL improved glottic visualization compared with DL. The FPS and overall intubation success rates were similar to other published prehospital studies using VL. The C-MAC provided reliable, high-quality video despite demanding prehospital conditions.

Native bacterial communities and Listeria monocytogenes survival in soils collected from the Lower Mainland of British Columbia, Canada.

Soil is an important reservoir for Listeria monocytogenes, a foodborne pathogen implicated in numerous produce-related outbreaks. Our objectives were to (i) compare the survival of L. monocytogenes among three soils, (ii) compare the native bacterial communities across these soils, and (iii) investigate relationships between L. monocytogenes survival, native bacterial communities, and soil properties. Listeria spp. populations were monitored on PALCAM agar in three soils inoculated with L. monocytogenes (∼5 × 106 CFU/g): conventionally farmed (CS), grassland transitioning to conventionally farmed (TS), and uncultivated grassland (GS). Bacterial diversity of the soils was analyzed using 16S rRNA targeted amplicon sequencing. A 2 log reduction of Listeria spp. was observed in all soils within 10 days, but at a significantly lower rate in GS (Fisher's least significant difference test; p < 0.05). Survival correlated with increased moisture and a neutral pH. GS showed the highest microbial diversity. Acidobacteria was the dominant phylum differentiating CS and TS from GS, and was negatively correlated with pH, carbon, nitrogen, and moisture. High moisture content and neutral pH are likely to increase the ability of L. monocytogenes to persist in soil. This study confirmed that native bacterial communities and short-term survival of L. monocytogenes varies across soils.

Improved white spruce (Picea glauca) genome assemblies and annotation of large gene families of conifer terpenoid and phenolic defense metabolism.

White spruce (Picea glauca), a gymnosperm tree, has been established as one of the models for conifer genomics. We describe the draft genome assemblies of two white spruce genotypes, PG29 and WS77111, innovative tools for the assembly of very large genomes, and the conifer genomics resources developed in this process. The two white spruce genotypes originate from distant geographic regions of western (PG29) and eastern (WS77111) North America, and represent elite trees in two Canadian tree-breeding programs. We present an update (V3 and V4) for a previously reported PG29 V2 draft genome assembly and introduce a second white spruce genome assembly for genotype WS77111. Assemblies of the PG29 and WS77111 genomes confirm the reconstructed white spruce genome size in the 20 Gbp range, and show broad synteny. Using the PG29 V3 assembly and additional white spruce genomics and transcriptomics resources, we performed MAKER-P annotation and meticulous expert annotation of very large gene families of conifer defense metabolism, the terpene synthases and cytochrome P450s. We also comprehensively annotated the white spruce mevalonate, methylerythritol phosphate and phenylpropanoid pathways. These analyses highlighted the large extent of gene and pseudogene duplications in a conifer genome, in particular for genes of secondary (i.e. specialized) metabolism, and the potential for gain and loss of function for defense and adaptation.

Mutational and structural analysis of diffuse large B-cell lymphoma using whole-genome sequencing.

Diffuse large B-cell lymphoma (DLBCL) is a genetically heterogeneous cancer composed of at least 2 molecular subtypes that differ in gene expression and distribution of mutations. Recently, application of genome/exome sequencing and RNA-seq to DLBCL has revealed numerous genes that are recurrent targets of somatic point mutation in this disease. Here we provide a whole-genome-sequencing-based perspective of DLBCL mutational complexity by characterizing 40 de novo DLBCL cases and 13 DLBCL cell lines and combining these data with DNA copy number analysis and RNA-seq from an extended cohort of 96 cases. Our analysis identified widespread genomic rearrangements including evidence for chromothripsis as well as the presence of known and novel fusion transcripts. We uncovered new gene targets of recurrent somatic point mutations and genes that are targeted by focal somatic deletions in this disease. We highlight the recurrence of germinal center B-cell-restricted mutations affecting genes that encode the S1P receptor and 2 small GTPases (GNA13 and GNAI2) that together converge on regulation of B-cell homing. We further analyzed our data to approximate the relative temporal order in which some recurrent mutations were acquired and demonstrate that ongoing acquisition of mutations and intratumoral clonal heterogeneity are common features of DLBCL. This study further improves our understanding of the processes and pathways involved in lymphomagenesis, and some of the pathways mutated here may indicate new avenues for therapeutic intervention.

Mutational evolution in a lobular breast tumour profiled at single nucleotide resolution.

Recent advances in next generation sequencing have made it possible to precisely characterize all somatic coding mutations that occur during the development and progression of individual cancers. Here we used these approaches to sequence the genomes (>43-fold coverage) and transcriptomes of an oestrogen-receptor-alpha-positive metastatic lobular breast cancer at depth. We found 32 somatic non-synonymous coding mutations present in the metastasis, and measured the frequency of these somatic mutations in DNA from the primary tumour of the same patient, which arose 9 years earlier. Five of the 32 mutations (in ABCB11, HAUS3, SLC24A4, SNX4 and PALB2) were prevalent in the DNA of the primary tumour removed at diagnosis 9 years earlier, six (in KIF1C, USP28, MYH8, MORC1, KIAA1468 and RNASEH2A) were present at lower frequencies (1-13%), 19 were not detected in the primary tumour, and two were undetermined. The combined analysis of genome and transcriptome data revealed two new RNA-editing events that recode the amino acid sequence of SRP9 and COG3. Taken together, our data show that single nucleotide mutational heterogeneity can be a property of low or intermediate grade primary breast cancers and that significant evolution can occur with disease progression.

Panorama: a targeted proteomics knowledge base.

Panorama is a web application for storing, sharing, analyzing, and reusing targeted assays created and refined with Skyline,1 an increasingly popular Windows client software tool for targeted proteomics experiments. Panorama allows laboratories to store and organize curated results contained in Skyline documents with fine-grained permissions, which facilitates distributed collaboration and secure sharing of published and unpublished data via a web-browser interface. It is fully integrated with the Skyline workflow and supports publishing a document directly to a Panorama server from the Skyline user interface. Panorama captures the complete Skyline document information content in a relational database schema. Curated results published to Panorama can be aggregated and exported as chromatogram libraries. These libraries can be used in Skyline to pick optimal targets in new experiments and to validate peak identification of target peptides. Panorama is open-source and freely available. It is distributed as part of LabKey Server,2 an open source biomedical research data management system. Laboratories and organizations can set up Panorama locally by downloading and installing the software on their own servers. They can also request freely hosted projects on https://panoramaweb.org , a Panorama server maintained by the Department of Genome Sciences at the University of Washington.

Assembling the 20 Gb white spruce (Picea glauca) genome from whole-genome shotgun sequencing data.

UNLABELLED: White spruce (Picea glauca) is a dominant conifer of the boreal forests of North America, and providing genomics resources for this commercially valuable tree will help improve forest management and conservation efforts. Sequencing and assembling the large and highly repetitive spruce genome though pushes the boundaries of the current technology. Here, we describe a whole-genome shotgun sequencing strategy using two Illumina sequencing platforms and an assembly approach using the ABySS software. We report a 20.8 giga base pairs draft genome in 4.9 million scaffolds, with a scaffold N50 of 20,356 bp. We demonstrate how recent improvements in the sequencing technology, especially increasing read lengths and paired end reads from longer fragments have a major impact on the assembly contiguity. We also note that scalable bioinformatics tools are instrumental in providing rapid draft assemblies. AVAILABILITY: The Picea glauca genome sequencing and assembly data are available through NCBI (Accession#: ALWZ0100000000 PID: PRJNA83435). http://www.ncbi.nlm.nih.gov/bioproject/83435.

Genome and transcriptome analyses of the mountain pine beetle-fungal symbiont Grosmannia clavigera, a lodgepole pine pathogen.

In western North America, the current outbreak of the mountain pine beetle (MPB) and its microbial associates has destroyed wide areas of lodgepole pine forest, including more than 16 million hectares in British Columbia. Grosmannia clavigera (Gc), a critical component of the outbreak, is a symbiont of the MPB and a pathogen of pine trees. To better understand the interactions between Gc, MPB, and lodgepole pine hosts, we sequenced the ∼30-Mb Gc genome and assembled it into 18 supercontigs. We predict 8,314 protein-coding genes, and support the gene models with proteome, expressed sequence tag, and RNA-seq data. We establish that Gc is heterothallic, and report evidence for repeat-induced point mutation. We report insights, from genome and transcriptome analyses, into how Gc tolerates conifer-defense chemicals, including oleoresin terpenoids, as they colonize a host tree. RNA-seq data indicate that terpenoids induce a substantial antimicrobial stress in Gc, and suggest that the fungus may detoxify these chemicals by using them as a carbon source. Terpenoid treatment strongly activated a ∼100-kb region of the Gc genome that contains a set of genes that may be important for detoxification of these host-defense chemicals. This work is a major step toward understanding the biological interactions between the tripartite MPB/fungus/forest system.

An open-label feasibility trial of transdermal cannabidiol for hand osteoarthritis.

Hand osteoarthritis (OA) is an irreversible degenerative condition causing chronic pain and impaired functionality. Existing treatment options are often inadequate. Cannabidiol (CBD) has demonstrated analgesic and anti-inflammatory effects in preclinical models of arthritis. In this open-label feasibility trial, participants with symptomatically active hand OA applied a novel transdermal CBD gel (4% w/w) three times a day for four weeks to their most painful hand. Changes in daily self-reported pain scores were measured on a 0-10 Numeric Pain Rating Scale (NPRS). Hand functionality was determined via daily grip strength measures using a Bluetooth equipped squeeze ball and self-report questionnaire. Quality of life (QoL) ratings around sleep, anxiety, stiffness and fatigue were also measured. All self-report measures and grip strength data were gathered via smartphone application. Urinalysis was conducted at trial end to determine systemic absorption of CBD. Eighteen participants were consented and 15 completed the trial. Pain ratings were significantly reduced over time from pre-treatment baseline including current pain (- 1.91 ± 0.35, p < 0.0001), average pain (- 1.92 ± 0.35, p < 0.0001) and maximum pain (- 1.97 ± 0.34, p < 0.0001) (data represent mean reduction on a 0-10 NPRS scale ± standard error of the mean (SEM)). A significant increase in grip strength in the treated hand (p < 0.0001) was observed although self-reported functionality did not improve. There were significant (p < 0.005) improvements in three QoL measures: fatigue, stiffness and anxiety. CBD and its metabolites were detected at low concentrations in all urine samples. Measured reductions in pain and increases in grip strength seen during treatment reverted back towards baseline during the washout phase. In summary, pain, grip strength and QoL measures, using smartphone technology, was shown to improve over time following transdermal CBD application suggesting feasibility of this intervention in relieving osteoarthritic hand pain. Proof of efficacy, however, requires further confirmation in a placebo-controlled randomised trial.Trial registration: ANZCTR public trials registry (ACTRN12621001512819, 05/11/2021).

Updated genome assembly and annotation of Paenibacillus larvae, the agent of American foulbrood disease of honey bees.

BACKGROUND: As scientists continue to pursue various 'omics-based research, there is a need for high quality data for the most fundamental 'omics of all: genomics. The bacterium Paenibacillus larvae is the causative agent of the honey bee disease American foulbrood. If untreated, it can lead to the demise of an entire hive; the highly social nature of bees also leads to easy disease spread, between both individuals and colonies. Biologists have studied this organism since the early 1900s, and a century later, the molecular mechanism of infection remains elusive. Transcriptomics and proteomics, because of their ability to analyze multiple genes and proteins in a high-throughput manner, may be very helpful to its study. However, the power of these methodologies is severely limited without a complete genome; we undertake to address that deficiency here. RESULTS: We used the Illumina GAIIx platform and conventional Sanger sequencing to generate a 182-fold sequence coverage of the P. larvae genome, and assembled the data using ABySS into a total of 388 contigs spanning 4.5 Mbp. Comparative genomics analysis against fully-sequenced soil bacteria P. JDR2 and P. vortex showed that regions of poor conservation may contain putative virulence factors. We used GLIMMER to predict 3568 gene models, and named them based on homology revealed by BLAST searches; proteases, hemolytic factors, toxins, and antibiotic resistance enzymes were identified in this way. Finally, mass spectrometry was used to provide experimental evidence that at least 35% of the genes are expressed at the protein level. CONCLUSIONS: This update on the genome of P. larvae and annotation represents an immense advancement from what we had previously known about this species. We provide here a reliable resource that can be used to elucidate the mechanism of infection, and by extension, more effective methods to control and cure this widespread honey bee disease.

Smokefree streets: a pilot study of methods to inform policy.

INTRODUCTION: Smokefree street policies are relatively rare, and little has been published on the methods for establishing an evidence base to inform such policy making. We aimed to (a) pilot methods for such data collection in New Zealand, a country where local governments are actively pursuing outdoor smokefree policies and (a) to provide data on smoking behavior, attitudes toward smokefree policies, and levels of smoke exposure on streets in Wellington. METHODS: Three methods were piloted: (a) systematic observation of smoking behavior by observers walking a standard route of major streets, the "Golden Mile" (GM) in Wellington (n = 42 observation runs); (b) measurement of fine particulate levels (PM(2.5)) along this route and with purposeful sampling in selected settings; and (c) an attitudinal survey of pedestrians along sections of this route. RESULTS: Each of the 3 methods proved to be feasible in this urban setting. A total of 932 smokers were observed during 21 hr of observation, an average of 7 smokers every 10 min of walking. Air monitoring indicated fine particulate exposure. Levels of (mean) PM(2.5) were 1.5 times higher during periods when smoking was observed than when they were not (9.3 vs. 6.3 µg/m(3), p = .002). Dose-response patterns were observed for smoking proximity and for smoker numbers. Surveying pedestrians (n = 220) with a brief questionnaire achieved an 81% response rate and was able to identify variation in support for a smokefree GM by different groups (overall support was 55.9%, 95% CI = 49.3%-62.4%). Reasons for support were also identified, for example, perceived health hazards, at 34.1%, was the main reason. CONCLUSIONS: These methods can provide information that may contribute to the smokefree streets policymaking process and may also be relevant to informing other smokefree outdoor policies.

Base excision repair deficiency signatures implicate germline and somatic MUTYH aberrations in pancreatic ductal adenocarcinoma and breast cancer oncogenesis.

We report a case of early-onset pancreatic ductal adenocarcinoma in a patient harboring biallelic MUTYH germline mutations, whose tumor featured somatic mutational signatures consistent with defective MUTYH-mediated base excision repair and the associated driver KRAS transversion mutation p.Gly12Cys. Analysis of an additional 730 advanced cancer cases (N = 731) was undertaken to determine whether the mutational signatures were also present in tumors from germline MUTYH heterozygote carriers or if instead the signatures were only seen in those with biallelic loss of function. We identified two patients with breast cancer each carrying a pathogenic germline MUTYH variant with a somatic MUTYH copy loss leading to the germline variant being homozygous in the tumor and demonstrating the same somatic signatures. Our results suggest that monoallelic inactivation of MUTYH is not sufficient for C:G>A:T transversion signatures previously linked to MUTYH deficiency to arise (N = 9), but that biallelic complete loss of MUTYH function can cause such signatures to arise even in tumors not classically seen in MUTYH-associated polyposis (N = 3). Although defective MUTYH is not the only determinant of these signatures, MUTYH germline variants may be present in a subset of patients with tumors demonstrating elevated somatic signatures possibly suggestive of MUTYH deficiency (e.g., COSMIC Signature 18, SigProfiler SBS18/SBS36, SignatureAnalyzer SBS18/SBS36).

Therapeutic Implication of Genomic Landscape of Adult Metastatic Sarcoma.

PURPOSE: This study investigated therapeutic potential of integrated genome and transcriptome profiling of metastatic sarcoma, a rare but extremely heterogeneous group of aggressive mesenchymal malignancies with few systemic therapeutic options. METHODS: Forty-three adult patients with advanced or metastatic non-GI stromal tumor sarcomas of various histology subtypes who were enrolled in the Personalized OncoGenomics program at BC Cancer were included in this study. Fresh tumor tissues along with blood samples underwent whole-genome and transcriptome sequencing. RESULTS: The most frequent genomic alterations in this cohort are large-scale structural variation and somatic copy number variation. Outlier RNA expression as well as somatic copy number variations, structural variations, and small mutations together suggest the presence of one or more potential therapeutic targets in the majority of patients in our cohort. Point mutations or deletions in known targetable cancer genes are rare; for example, tuberous sclerosis complex 2 provides a rationale for targeting the mammalian target of rapamycin pathway, resulting in a few patients with exceptional clinical benefit from everolimus. In addition, we observed recurrent 17p11-12 amplifications, which seem to be a sarcoma-specific event. This may suggest that this region harbors an oncogene(s) that is significant for sarcoma tumorigenesis. Furthermore, some sarcoma tumors carrying a distinct mutational signature suggestive of homologous recombination deficiency seem to demonstrate sensitivity to double-strand DNA-damaging agents. CONCLUSION: Integrated large-scale genomic analysis may provide insights into potential therapeutic targets as well as novel biologic features of metastatic sarcomas that could fuel future experimental and clinical research and help design biomarker-driven basket clinical trials for novel therapeutic strategies.

Whole-genome and transcriptome profiling of a metastatic thyroid-like follicular renal cell carcinoma.

Thyroid-like follicular renal cell carcinoma (TLFRCC) is a rare cancer with few reports of metastatic disease. Little is known regarding genomic characteristics and therapeutic targets. We present the clinical, pathologic, genomic, and transcriptomic analyses of a case of a 27-yr-old male with TLFRCC who presented initially with bone metastases of unknown primary. Genomic DNA from peripheral blood and metastatic tumor samples were sequenced. A transcriptome of 280 million sequence reads was generated from the same tumor sample. Tumor somatic expression profiles were analyzed to detect aberrant expression. Genomic and transcriptomic data sets were integrated to reveal dysregulation in pathways and identify potential therapeutic targets. Integrative genomic analysis with The Cancer Genome Atlas (TCGA) data set revealed the following outliers in gene expression profiles: CDK6 (81st percentile), MYC (99th percentile), AR (100th percentile), PDGFRA and PDGFRB (99th and 100th percentiles, respectively), and MAP2K2 (86th percentile). The patient received first-line sunitinib to target PDGFRA and PDGFRB and had stable disease for >6 mo, followed by nivolumab upon progression. To the authors' knowledge, this is the first reported case of comprehensive somatic genomic analyses in a patient with metastatic TLFRCC. Somatic analyses provided molecular confirmation of the primary site of cancer and potential therapeutic strategies in a rare disease with little evidence of efficacy on systemic therapy.

Genome-Enhanced Detection and Identification (GEDI) of plant pathogens.

Plant diseases caused by fungi and Oomycetes represent worldwide threats to crops and forest ecosystems. Effective prevention and appropriate management of emerging diseases rely on rapid detection and identification of the causal pathogens. The increase in genomic resources makes it possible to generate novel genome-enhanced DNA detection assays that can exploit whole genomes to discover candidate genes for pathogen detection. A pipeline was developed to identify genome regions that discriminate taxa or groups of taxa and can be converted into PCR assays. The modular pipeline is comprised of four components: (1) selection and genome sequencing of phylogenetically related taxa, (2) identification of clusters of orthologous genes, (3) elimination of false positives by filtering, and (4) assay design. This pipeline was applied to some of the most important plant pathogens across three broad taxonomic groups: Phytophthoras (Stramenopiles, Oomycota), Dothideomycetes (Fungi, Ascomycota) and Pucciniales (Fungi, Basidiomycota). Comparison of 73 fungal and Oomycete genomes led the discovery of 5,939 gene clusters that were unique to the targeted taxa and an additional 535 that were common at higher taxonomic levels. Approximately 28% of the 299 tested were converted into qPCR assays that met our set of specificity criteria. This work demonstrates that a genome-wide approach can efficiently identify multiple taxon-specific genome regions that can be converted into highly specific PCR assays. The possibility to easily obtain multiple alternative regions to design highly specific qPCR assays should be of great help in tackling challenging cases for which higher taxon-resolution is needed.

Genome sequences of six Phytophthora species threatening forest ecosystems.

The Phytophthora genus comprises of some of the most destructive plant pathogens and attack a wide range of hosts including economically valuable tree species, both angiosperm and gymnosperm. Many known species of Phytophthora are invasive and have been introduced through nursery and agricultural trade. As part of a larger project aimed at utilizing genomic data for forest disease diagnostics, pathogen detection and monitoring (The TAIGA project: Tree Aggressors Identification using Genomic Approaches; http://taigaforesthealth.com/), we sequenced the genomes of six important Phytophthora species that are important invasive pathogens of trees and a serious threat to the international trade of forest products. This genomic data was used to develop highly sensitive and specific detection assays and for genome comparisons and to make evolutionary inferences and will be useful to the broader plant and tree health community. These WGS data have been deposited in the International Nucleotide Sequence Database Collaboration (DDBJ/ENA/GenBank) under the accession numbers AUPN01000000, AUVH01000000, AUWJ02000000, AUUF02000000, AWVV02000000 and AWVW02000000.