Validation of an automated, end-to-end metagenomic sequencing assay for agnostic detection of respiratory viruses.

BACKGROUND: Current molecular diagnostics are limited in the number and type of detectable pathogens. Metagenomic next generation sequencing (mNGS) is an emerging, and increasingly feasible, pathogen-agnostic diagnostic approach. Translational barriers prohibit the widespread adoption of this technology in clinical laboratories. We validate an end-to-end mNGS assay for detection of respiratory viruses. Our assay is optimized to reduce turnaround time, lower cost-per-sample, increase throughput, and deploy secure and actionable bioinformatic results. METHODS: We validated our assay using residual nasopharyngeal swab specimens from Vancouver General Hospital (n = 359), RT-PCR-positive, or negative for Influenza, SARS-CoV-2, and RSV. We quantified sample stability, assay precision, the effect of background nucleic acid levels, and analytical limits of detection. Diagnostic performance metrics were estimated. RESULTS: We report that our mNGS assay is highly precise, semi-quantitative, with analytical limits of detection ranging from 103-104 copies/mL. Our assay is highly specific (100%) and sensitive (61.9% Overall: 86.8%; RT-PCR Ct < 30). Multiplexing capabilities enable processing of up to 55-specimens simultaneously on an Oxford Nanopore GridION device, with results reported within 12-hours. CONCLUSIONS: This study outlines the diagnostic performance and feasibility of mNGS for respiratory viral diagnostics, infection control, and public health surveillance. We addressed translational barriers to widespread mNGS adoption.

Rapid microdissection of tissue sections via laser ablation.

We demonstrate a method for tissue microdissection using scanning laser ablation that is approximately two orders of magnitude faster than conventional laser capture microdissection. Our novel approach uses scanning laser optics and a slide coating under the tissue that can be excited by the laser to selectively eject regions of tissue for further processing. Tissue was dissected at 0.117 s/mm2 without reduction in yield, sequencing insert size or base quality compared with undissected tissue. From eight cases, 58-416 mm2 of tissue was obtained from one to four slides in 7-48 seconds total dissection time per case. These samples underwent exome sequencing and we found the variant allelic fraction increased in regions enriched for tumour as expected. This suggests that our ablation technique may be useful as a tool in both clinical and research labs.

Development of the quantitative PET prostate phantom (Q3P) for improved quality assurance of 18F-PSMA PET imaging in metastatic prostate cancer.

BACKGROUND: Phantoms are commonly used to evaluate and compare the performance of imaging systems given the known ground truth. Positron emission tomography (PET) scanners are routinely validated using the NEMA image quality phantom, in which lesions are modeled using 10 to 37 mm fillable spheres. The NEMA phantom neglects, however, to model focal (3-10-mm), high-uptake lesions that are increasingly observed in prostate-specific membrane antigen (PSMA) PET images. PSMA-targeting radiopharmaceuticals allow for enhanced detection of metastatic prostate cancers. As such, there is significant need to develop an updated phantom which considers both the quantitative and lesion detectability of this new paradigm in oncological PET imaging. PURPOSE: In this work, we present the Quantitative PET Prostate Phantom (Q3P); a portable and modular phantom that can be used to improve and harmonize imaging protocols for 18F-PSMA PET scans. METHODS: A one-piece cylindrical phantom was designed effectively in two halves, which we call modules. Module 1 was designed to mimic lesions in the presence of background, and Module 2 mimicked very high contrast conditions (i.e., very low background) that can be observed in 18F-PSMA PET scans. Shell-less radioactive spheres (3-16-mm) were cast using epoxy resin mixed with sodium-22 (22Na), a long half-life positron emitter with positron range similar to 18F. To establish realistic lesion contrast, the 22Na spheres were mounted in a cylindrical chamber that can be filled with an 18F background (module 1). Thirteen exchangeable spherical cavity inserts (3-37-mm) were machined in two parts and solvent welded together, and filled with 18F (50 kBq/mL) to model lesions with very high contrast (module 2). Five 2.5-min PET scans were acquired on a 5-ring GE Discovery MI PET/CT scanner (General Electric, USA). Lesions were segmented using 41% of SUVmax fixed thresholding (41% FT) and recovery coefficients (RCs) were computed from 5 noise realizations. RESULTS: The manufactured phantom is portable (5.7 kg) and scan preparation takes less than 40 min. The total 22Na activity is 250 kBq, allowing it to be shipped as an exempt package under International Atomic Energy Agency (IAEA) regulations. Recovery coefficients, computed using PSF modeling and no post-reconstruction smoothing, were 130.3% (16 mm), 147.1% (10 mm), 87.2% (6 mm), and 7.0% (3 mm) for RCmax, which decreased to 91.1% (16 mm), 90.6% (10 mm), 53.2% (6 mm), and 3.6% (3 mm) for RCmean in the 22Na spheres. Comparatively, 18F sphere recovery was 110.7% (17 mm), 123.6% (10 mm), 106.5% (7 mm), and 23.3% (3 mm) for RCmax, which was reduced to 76.7% (17 mm), 77.7% (10 mm), 66.8% (7 mm), and 13.5% (3 mm), for RCmean. CONCLUSIONS: A standardized imaging phantom was developed for lesion quantification assessment in 18F-PSMA PET images. The phantom is configurable, providing users with the opportunity to modify background activity levels or sphere sizes according to clinical demands. Distributed to the community, the Q3P phantom has the potential to enable better assessment of lesion quantification and harmonization of 18F-PSMA PET imaging, which may lead to more robust predictive metrics and better outcome prediction in metastatic prostate cancer.

A high-throughput pipeline for DNA/RNA/small RNA purification from tissue samples for sequencing.

High-throughput total nucleic acid (TNA) purification methods based on solid-phase reversible immobilization (SPRI) beads produce TNA suitable for both genomic and transcriptomic applications. Even so, small RNA species, including miRNA, bind weakly to SPRI beads under standard TNA purification conditions, necessitating a separate workflow using column-based methods that are difficult to automate. Here, an SPRI-based high-throughput TNA purification protocol that recovers DNA, RNA and small RNA, called GSC-modified RLT+ Aline bead-based protocol (GRAB-ALL), which incorporates modifications to enhance small RNA recovery is presented. GRAB-ALL was benchmarked against existing nucleic acid purification workflows and GRAB-ALL efficiently purifies TNA, including small RNA, for next-generation sequencing applications in a plate-based format suitable for automated high-throughput sample preparation.

Design of a double acting pneumatic cartilage loading device for magnetic resonance imaging.

Studies of osteoarthritis initiation and progression that measure strain in cartilage require physiological loading levels. Many studies use magnetic resonance (MR) imaging, which necessitates a MR-compatible loading device. In this study, the design and validation of a new device, the cartilage compressive actuator (CCA), is presented. The CCA is designed for high-field (e.g., 9.4 T) small-bore MR scanners, and meets a number of design criteria. These criteria include capability for testing bone-cartilage samples, MR compatibility, constant load and incremental strain application, a water-tight specimen chamber, remote control, and real time displacement feedback. The mechanical components in the final design include an actuating piston, a connecting chamber, and a sealed specimen chamber. An electro-pneumatic system applies compression, and an optical Fibre-Bragg grating (FBG) sensor provides live displacement feedback. A logarithmic relationship was observed between force exerted by the CCA and pressure (R2 = 0.99), with a peak output force of 653 ± 2 N. The relationship between FBG sensor wavelength and displacement was linear when calibrated both outside (R2 = 0.99) and inside (R2 = 0.98) the MR scanner. Average slope was similar between the two validation tests, with a slope of -4.2 nm/mm observed inside the MR scanner and -4.3 to -4.5 nm/mm observed outside the MR scanner. This device meets all design criteria and represents an improvement over published designs. Future work should incorporate a closed feedback loop to allow for cyclical loading of specimens.

Extraction-free clinical detection of SARS-CoV-2 virus from saline gargle samples using Hamilton STARlet liquid handler.

As part of the COVID-19 pandemic, clinical laboratories have been faced with massive increases in testing, resulting in sample collection systems, reagent, and staff shortages. We utilized self-collected saline gargle samples to optimize high throughput SARS-CoV-2 multiplex polymerase chain reaction (PCR) testing in order to minimize cost and technologist time. This was achieved through elimination of nucleic acid extraction and automation of sample handling on a widely available robotic liquid handler, Hamilton STARlet. A customized barcode scanning script for reading the sample ID by the Hamilton STARlet's software system was developed to allow primary tube sampling. Use of pre-frozen SARS-CoV-2 assay reaction mixtures reduced assay setup time. In both validation and live testing, the assay produced no false positive or false negative results. Of the 1060 samples tested during validation, 3.6% (39/1060) of samples required retesting as they were either single gene positive, had internal control failure or liquid aspiration error. Although the overall turnaround time was only slightly faster in the automated workflow (185 min vs 200 min), there was a 76% reduction in hands-on time, potentially reducing staff fatigue and burnout. This described process from sample self-collection to automated direct PCR testing significantly reduces the total burden on healthcare systems in terms of human resources and reagent requirements.

Automated Library Construction and Analysis for High-Throughput Nanopore Sequencing of SARS-CoV-2.

BACKGROUND: To support the implementation of high-throughput pipelines suitable for SARS-CoV-2 sequencing and analysis in a clinical laboratory, we developed an automated sample preparation and analysis workflow. METHODS: We used the established ARTIC protocol with approximately 400 bp amplicons sequenced on Oxford Nanopore's MinION. Sequences were analyzed using Nextclade, assigning both a clade and quality score to each sample. RESULTS: A total of 2179 samples on twenty-five 96-well plates were sequenced. Plates of purified RNA were processed within 12 h, sequencing required up to 24 h, and analysis of each pooled plate required 1 h. The use of samples with known threshold cycle (Ct) values enabled normalization, acted as a quality control check, and revealed a strong correlation between sample Ct values and successful analysis, with 85% of samples with Ct < 30 achieving a "good" Nextclade score. Less abundant samples responded to enrichment with the fraction of Ct > 30 samples achieving a "good" classification rising by 60% after addition of a post-ARTIC PCR normalization. Serial dilutions of 3 variant of concern samples, diluted from approximately Ct = 16 to approximately Ct = 50, demonstrated successful sequencing to Ct = 37. The sample set contained a median of 24 mutations per sample and a total of 1281 unique mutations with reduced sequence read coverage noted in some regions of some samples. A total of 10 separate strains were observed in the sample set, including 3 variants of concern prevalent in British Columbia in the spring of 2021. CONCLUSIONS: We demonstrated a robust automated sequencing pipeline that takes advantage of input Ct values to improve reliability.

Construction of Strand-seq libraries in open nanoliter arrays.

Single-cell Strand-seq generates directional genomic information to study DNA repair, assemble genomes, and map structural variation onto chromosome-length haplotypes. We report a nanoliter-volume, one-pot (OP) Strand-seq library preparation protocol in which reagents are added cumulatively, DNA purification steps are avoided, and enzymes are inactivated with a thermolabile protease. OP-Strand-seq libraries capture 10%-25% of the genome from a single-cell with reduced costs and increased throughput.

Optimization of magnetic bead-based nucleic acid extraction for SARS-CoV-2 testing using readily available reagents.

The COVID-19 pandemic has highlighted the need for generic reagents and flexible systems in diagnostic testing. Magnetic bead-based nucleic acid extraction protocols using 96-well plates on open liquid handlers are readily amenable to meet this need. Here, one such approach is rigorously optimized to minimize cross-well contamination while maintaining sensitivity.

Reduction in Doses to Organs at Risk and Normal Tissue During Breast Radiation Therapy With a Carbon-Fiber Adjustable Reusable Accessory.

PURPOSE: This pilot study (ClinicalTrials.gov NCT04543851) investigates a novel breast positioning device using a low density, high tensile carbon-fiber cradle to support the breast, remove the inframammary fold, and reduce dose to organs at risk for whole breast radiation therapy in the supine position. METHODS AND MATERIALS: Thirty patients with inframammary folds ≥1 cm or lateral ptosis in supine treatment position were planned with standard positioning and with a carbon-fiber Adjustable Reusable Accessory (CARA) breast support. Twenty patients received whole breast with or without regional nodal irradiation with 42.5 Gy in 16 fractions or 50 Gy in 25 fractions using CARA. Median body mass index was 32 in this study. RESULTS: CARA removed all inframammary folds and reduced V20Gyipsilateral lung, V105%breast, and V50% body, without compromising target coverage. Median (range) V20Gyipsilateral lung for whole breast radiation therapy was 12.3% (1.4%-28.7%) with standard of care versus 10.9% (1.2%-17.3%) with CARA (Wilcoxon P = .005). Median V105% breast was 8.0% (0.0%-29%) with standard of care versus 4.0% (0.0%-23%) with CARA (P = .006) and median V50% body was 3056 mL (1476-5285 mL) versus 2780 mL (1415-5123 mL) with CARA (P = .001). CARA was compatible with deep inspiration breath hold and achieved median V25Gyheart = 0.1% (range 0%-1.9%) for all patients with left breast cancer. Skin reactions with CARA were consistent with historical data and daily variation in treatment setup was consistent with standard supine positioning. CONCLUSIONS: CARA can reduce V105%breast, lung and normal tissue dose, and remove the inframammary fold for breast patients with large or pendulous breasts and high body mass index treated in the supine position, without compromising target coverage. CARA will undergo further study in a randomized controlled trial.

A Scalable Strand-Specific Protocol Enabling Full-Length Total RNA Sequencing From Single Cells.

RNA sequencing (RNAseq) has been widely used to generate bulk gene expression measurements collected from pools of cells. Only relatively recently have single-cell RNAseq (scRNAseq) methods provided opportunities for gene expression analyses at the single-cell level, allowing researchers to study heterogeneous mixtures of cells at unprecedented resolution. Tumors tend to be composed of heterogeneous cellular mixtures and are frequently the subjects of such analyses. Extensive method developments have led to several protocols for scRNAseq but, owing to the small amounts of RNA in single cells, technical constraints have required compromises. For example, the majority of scRNAseq methods are limited to sequencing only the 3' or 5' termini of transcripts. Other protocols that facilitate full-length transcript profiling tend to capture only polyadenylated mRNAs and are generally limited to processing only 96 cells at a time. Here, we address these limitations and present a novel protocol that allows for the high-throughput sequencing of full-length, total RNA at single-cell resolution. We demonstrate that our method produced strand-specific sequencing data for both polyadenylated and non-polyadenylated transcripts, enabled the profiling of transcript regions beyond only transcript termini, and yielded data rich enough to allow identification of cell types from heterogeneous biological samples.

Whole-slide laser microdissection for tumour enrichment.

The practical application of genome-scale technologies to precision oncology research requires flexible tissue processing strategies that can be used to differentially select both tumour and normal cell populations from formalin-fixed, paraffin-embedded tissues. As tumour sequencing scales towards clinical implementation, practical difficulties in scheduling and obtaining fresh tissue biopsies at scale, including blood samples as surrogates for matched 'normal' DNA, have focused attention on the use of formalin-preserved clinical samples collected routinely for diagnostic purposes. In practice, such samples often contain both tumour and normal cells which, if correctly partitioned, could be used to profile both tumour and normal genomes, thus identifying somatic alterations. Here we report a semi-automated method for laser microdissecting entire slide-mounted tissue sections to enrich for cells of interest with sufficient yield for whole genome and transcriptome sequencing. Using this method, we demonstrated enrichment of tumour material from mixed tumour-normal samples by up to 67%. Leveraging new methods that allow for the extraction of high-quality nucleic acids from small amounts of formalin-fixed tissues, we further showed that the method was successful in yielding sequence data of sufficient quality for use in BC Cancer's Personalized OncoGenomics (POG) program. © 2020 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Complete Mitochondrial Genome of a Gymnosperm, Sitka Spruce (Picea sitchensis), Indicates a Complex Physical Structure.

Plant mitochondrial genomes vary widely in size. Although many plant mitochondrial genomes have been sequenced and assembled, the vast majority are of angiosperms, and few are of gymnosperms. Most plant mitochondrial genomes are smaller than a megabase, with a few notable exceptions. We have sequenced and assembled the complete 5.5-Mb mitochondrial genome of Sitka spruce (Picea sitchensis), to date, one of the largest mitochondrial genomes of a gymnosperm. We sequenced the whole genome using Oxford Nanopore MinION, and then identified contigs of mitochondrial origin assembled from these long reads based on sequence homology to the white spruce mitochondrial genome. The assembly graph shows a multipartite genome structure, composed of one smaller 168-kb circular segment of DNA, and a larger 5.4-Mb single component with a branching structure. The assembly graph gives insight into a putative complex physical genome structure, and its branching points may represent active sites of recombination.

Evaluation of protocols for rRNA depletion-based RNA sequencing of nanogram inputs of mammalian total RNA.

Next generation RNA-sequencing (RNA-seq) is a flexible approach that can be applied to a range of applications including global quantification of transcript expression, the characterization of RNA structure such as splicing patterns and profiling of expressed mutations. Many RNA-seq protocols require up to microgram levels of total RNA input amounts to generate high quality data, and thus remain impractical for the limited starting material amounts typically obtained from rare cell populations, such as those from early developmental stages or from laser micro-dissected clinical samples. Here, we present an assessment of the contemporary ribosomal RNA depletion-based protocols, and identify those that are suitable for inputs as low as 1-10 ng of intact total RNA and 100-500 ng of partially degraded RNA from formalin-fixed paraffin-embedded tissues.

Clonal Decomposition and DNA Replication States Defined by Scaled Single-Cell Genome Sequencing.

Accurate measurement of clonal genotypes, mutational processes, and replication states from individual tumor-cell genomes will facilitate improved understanding of tumor evolution. We have developed DLP+, a scalable single-cell whole-genome sequencing platform implemented using commodity instruments, image-based object recognition, and open source computational methods. Using DLP+, we have generated a resource of 51,926 single-cell genomes and matched cell images from diverse cell types including cell lines, xenografts, and diagnostic samples with limited material. From this resource we have defined variation in mitotic mis-segregation rates across tissue types and genotypes. Analysis of matched genomic and image measurements revealed correlations between cellular morphology and genome ploidy states. Aggregation of cells sharing copy number profiles allowed for calculation of single-nucleotide resolution clonal genotypes and inference of clonal phylogenies and avoided the limitations of bulk deconvolution. Finally, joint analysis over the above features defined clone-specific chromosomal aneuploidy in polyclonal populations.

An in vitro study for the dosimetric and radiobiological validation of respiratory gating in conventional and hypofractionated radiotherapy of the lung: effect of dose, dose rate, and breathing pattern.

Stereotactic body radiotherapy (SBRT) of the lung has become a standard of care for early-stage inoperable non-small cell lung cancer (NSCLC). A common strategy to manage respiratory motion is gating, which inevitably results in an increase in treatment time, especially in irregularly-breathing patients. Flattening-filter free (FFF) beams allow for delivery of the treatment at a higher dose rate, therefore counteracting the lengthened treatment time due to frequent interruption of the beam during gated radiotherapy. In this study, we perform our in vitro evaluation of the dosimetric and radiobiological effect of gated lung SBRT with simultaneous integrated boost (SIB) using both flattened and FFF beams. A moving thorax-shaped phantom with inserts and applicators was used for simulation, planning, gated treatment delivery measurements and in vitro tests. The effects of gating window, dose rate, and breathing pattern were evaluated. Planned doses represented a typical conventional fractionation, 200 cGy per fraction with SIB to 240 cGy, flattened beam only, and SBRT, 800 cGy with SIB to 900 cGy, flattened and FFF beams. Ideal, as well as regular and irregular patient-specific breathing patterns with and without gating were used. A survival assay for lung adenocarcinoma A549 cell line was performed. Delivered dose was within 6% for locations planned to receive 200 and 800 cGy and within 4% for SIB locations. Time between first beam-on and last beam-off varied from approximately 1.5 min for conventional fractionation, 200/240 cGy, to 10.5 min for gated SBRT, 800/900 cGy doses, flattened beam and irregular breathing motion pattern. With FFF beams dose delivery time was shorter by a factor of 2-3, depending on the gating window and breathing pattern. We have found that, for the most part, survival depended on dose and not on dose rate, gating window, or breathing regularity.

Rise of the Machines: Advances in Deep Learning for Cancer Diagnosis.

Deep learning refers to a set of computer models that have recently been used to make unprecedented progress in the way computers extract information from images. These algorithms have been applied to tasks in numerous medical specialties, most extensively radiology and pathology, and in some cases have attained performance comparable to human experts. Furthermore, it is possible that deep learning could be used to extract data from medical images that would not be apparent by human analysis and could be used to inform on molecular status, prognosis, or treatment sensitivity. In this review, we outline the current developments and state-of-the-art in applying deep learning for cancer diagnosis, and discuss the challenges in adapting the technology for widespread clinical deployment.

A high-throughput protocol for isolating cell-free circulating tumor DNA from peripheral blood.

The analysis of cell-free circulating tumor DNA (ctDNA) is potentially a less invasive, more dynamic assessment of cancer progression and treatment response than characterizing solid tumor biopsies. Standard isolation methods require separation of plasma by centrifugation, a time-consuming step that complicates automation. To address these limitations, we present an automatable magnetic bead-based ctDNA isolation method that eliminates centrifugation to purify ctDNA directly from peripheral blood (PB). To develop and test our method, ctDNA from cancer patients was purified from PB and plasma. We found that allelic fractions of somatic single-nucleotide variants from target gene capture libraries were comparable, indicating that the PB ctDNA purification method may be a suitable replacement for the plasma-based protocols currently in use.

Imaging-Based 3-Dimensional Printing for Improved Maxillofacial Presurgical Planning: A Single Center Case Series.

PURPOSE: 3-D printing is an increasingly widespread technology that allows physical models to be constructed based on cross-sectional medical imaging data. We sought to develop a pipeline for production of 3-dimensional (3-D) models for presurgical planning and assess the value of these models for surgeons and patients. METHODS: In this institutional review board-approved, single-center case series, participating surgeons identified cases for 3-D model printing, and after obtaining patient consent, a 3-D model was produced for each of the 7 participating patients based on preoperative cross-sectional imaging. Each model was given to the surgeon to use during the surgical consent discussion and preoperative planning. Patients and surgeons completed questionnaires evaluating the quality and usefulness of the models. RESULTS: The 3-D models improved surgeon confidence in their operative approach, influencing the choice of operative approach in the majority of cases. Patients and surgeons reported that the model improved patient comprehension of the surgery during the consent discussion, including risks and benefits of the surgery. Model production time was as little as 4 days, and the average per-model cost was $350. CONCLUSIONS: 3-D printed models are useful presurgical tools from both surgeon and patient perspectives. Development of local hospital-based 3-D printing capabilities enables model production with rapid turnaround and modest cost, representing a value-added service for radiologists to offer their surgical colleagues.

MAVIS: merging, annotation, validation, and illustration of structural variants.

Summary: Reliably identifying genomic rearrangements and interpreting their impact is a key step in understanding their role in human cancers and inherited genetic diseases. Many short read algorithmic approaches exist but all have appreciable false negative rates. A common approach is to evaluate the union of multiple tools increasing sensitivity, followed by filtering to retain specificity. Here we describe an application framework for the rapid generation of structural variant consensus, unique in its ability to visualize the genetic impact and context as well as process both genome and transcriptome data. Availability and implementation: http://mavis.bcgsc.ca. Supplementary information: Supplementary data are available at Bioinformatics online.