Assessment of the latest NGS enrichment capture methods in clinical context

Recurrence Prediction by Circulating Tumor DNA in the Patient with Colorectal Liver Metastases After Hepatectomy: A Prospective Biomarker Study

BACKGROUND

The recurrence rate after hepatic resection of colorectal liver metastases (CRLM) remains high. This study aimed to investigate postoperative circulating tumor DNA (ctDNA) based on ultra-deep next-generation sequencing (NGS) to predict patient recurrence and survival.

METHODS

Using the high-throughput NGS method tagged with a dual-indexed unique molecular identifier, named the CRLM-specific 25-gene panel (J25), this study sequenced ctDNA in peripheral blood samples collected from 134 CRLM patients who underwent hepatectomy after postoperative day 6.

RESULTS

Of 134 samples, 42 (31.3%) were shown to be ctDNA-positive, and 37 resulted in recurrence. Kaplan-Meier survival analysis showed that disease-free survival (DFS) in the ctDNA-positive subgroup was significantly shorter than in the ctDNA-negative subgroup (hazard ratio [HR], 2.96; 95% confidence interval [CI], 1.91-4.6; p < 0.05). When the 42 ctDNA-positive samples were further divided by the median of the mean allele frequency (AF, 0.1034%), the subgroup with higher AFs showed a significantly shorter DFS than the subgroup with lower AFs (HR, 1.98; 95% CI, 1.02-3.85; p < 0.05). The ctDNA-positive patients who received adjuvant chemotherapy longer than 2 months showed a significantly longer DFS than those who received treatment for 2 months or less (HR, 0.377; 95% CI, 0.189-0.751; p < 0.05). Uni- and multivariable Cox regression indicated two factors independently correlated with prognosis: ctDNA positivity and no preoperative chemotherapy.

CONCLUSION

The study demonstrated that ctDNA status 6 days postoperatively could sensitively and accurately predict recurrence for patients with CRLM using the J25 panel.

Whole genome sequencing of Mycobacterium bovis directly from clinical tissue samples without culture

Advancement in next generation sequencing offers the possibility of routine use of whole genome sequencing (WGS) for Mycobacterium bovis (M. bovis) genomes in clinical reference laboratories. To date, the M. bovis genome could only be sequenced if the mycobacteria were cultured from tissue. This requirement for culture has been due to the overwhelmingly large amount of host DNA present when DNA is prepared directly from a granuloma. To overcome this formidable hurdle, we evaluated the usefulness of an RNA-based targeted enrichment method to sequence M. bovis DNA directly from tissue samples without culture. Initial spiking experiments for method development were established by spiking DNA extracted from tissue samples with serially diluted M. bovis BCG DNA at the following concentration range: 0.1 ng/μl to 0.1 pg/μl (10^-1 to 10^-4). Library preparation, hybridization and enrichment was performed using SureSelect custom capture library RNA baits and the SureSelect XT HS2 target enrichment system for Illumina paired-end sequencing. The method validation was then assessed using direct WGS of M. bovis DNA extracted from tissue samples from naturally (n = 6) and experimentally (n = 6) infected animals with variable Ct values. Direct WGS of spiked DNA samples achieved 99.1% mean genome coverage (mean depth of coverage: 108×) and 98.8% mean genome coverage (mean depth of coverage: 26.4×) for tissue samples spiked with BCG DNA at 10^-1 (mean Ct value: 20.3) and 10^-2 (mean Ct value: 23.4), respectively. The M. bovis genome from the experimentally and naturally infected tissue samples was successfully sequenced with a mean genome coverage of 99.56% and depth of genome coverage ranging from 9.2× to 72.1×. The spoligoyping and M. bovis group assignment derived from sequencing DNA directly from the infected tissue samples matched that of the cultured isolates from the same sample. Our results show that direct sequencing of M. bovis DNA from tissue samples has the potential to provide accurate sequencing of M. bovis genomes significantly faster than WGS from cultures in research and diagnostic settings.

NGS method by library enrichment for rapid pestivirus purity testing in biologics

NGS sequencing was evaluated to understand its added value for animal health vaccine candidates. We have previously established the proof of concept for its application in purity testing on several Master Seeds. Here we evaluate the NGS method after enrichment to detect pestiviruses. To achieve this, we conducted a spiking study using 6 viruses, consisting of 3 pestiviruses and 3 other RNA-viruses at different concentrations into cell suspension. A deep Illumina random sequencing of all nucleic acids (DNA and RNA) was performed. The bioinformatics analysis including both assembly into contigs and annotation were processed using viral public databases for the spiked viruses' identification. Here we present the results of spiking experiments for the simultaneous spike of 6 viruses at 100-10 and 1 TCID₅₀/ml. Using Illumina sequencing, the 3 pestiviruses were all detected at the highest concentration, and even at the lowest one such as 1 TCID₅₀/ml for CSFV. Regarding the other viruses, they were not detected at all. Overall, the study showed consistent results for specific detection of pestiviruses with an increase of sensitivity after enrichment. The sensitivity of NGS evaluated by virus spiking experiments of cells demonstrated that NGS method is a valuable and sensitive tool for specific agent detection required in purity testing during vaccine development. This NGS method should be considered as an alternative tool of current purity testing for the prospective testing of biological products.

High-throughput method for the hybridisation-based targeted enrichment of long genomic fragments for PacBio third-generation sequencing

Hybridisation-based targeted enrichment is a widely used and well-established technique in high-throughput second-generation short-read sequencing. Despite the high potential to genetically resolve highly repetitive and variable genomic sequences by, for example PacBio third-generation sequencing, targeted enrichment for long fragments has not yet established the same high-throughput due to currently existing complex workflows and technological dependencies. We here describe a scalable targeted enrichment protocol for fragment sizes of >7 kb. For demonstration purposes we developed a custom blood group panel of challenging loci. Test results achieved > 65% on-target rate, good coverage (142.7×) and sufficient coverage evenness for both non-paralogous and paralogous targets, and sufficient non-duplicate read counts (83.5%) per sample for a highly multiplexed enrichment pool of 16 samples. We genotyped the blood groups of nine patients employing highly accurate phased assemblies at an allelic resolution that match reference blood group allele calls determined by SNP array and NGS genotyping. Seven Genome-in-a-Bottle reference samples achieved high recall (96%) and precision (99%) rates. Mendelian error rates were 0.04% and 0.13% for the included Ashkenazim and Han Chinese trios, respectively. In summary, we provide a protocol and first example for accurate targeted long-read sequencing that can be used in a high-throughput fashion.

Next Generation and Other Sequencing Technologies in Diagnostic Microbiology and Infectious Diseases

Next-generation sequencing (NGS) technologies have become increasingly available for use in the clinical microbiology diagnostic environment. There are three main applications of these technologies in the clinical microbiology laboratory: whole genome sequencing (WGS), targeted metagenomics sequencing and shotgun metagenomics sequencing. These applications are being utilized for initial identification of pathogenic organisms, the detection of antimicrobial resistance mechanisms and for epidemiologic tracking of organisms within and outside hospital systems. In this review, we analyze these three applications and provide a comprehensive summary of how these applications are currently being used in public health, basic research, and clinical microbiology laboratory environments. In the public health arena, WGS is being used to identify and epidemiologically track food borne outbreaks and disease surveillance. In clinical hospital systems, WGS is used to identify multi-drug-resistant nosocomial infections and track the transmission of these organisms. In addition, we examine how metagenomics sequencing approaches (targeted and shotgun) are being used to circumvent the traditional and biased microbiology culture methods to identify potential pathogens directly from specimens. We also expand on the important factors to consider when implementing these technologies, and what is possible for these technologies in infectious disease diagnosis in the next 5 years.

Performance comparison of Agilent new SureSelect All Exon v8 probes with v7 probes for exome sequencing

Exome sequencing is becoming a routine in health care, because it increases the chance of pinpointing the genetic cause of an individual patient's condition and thus making an accurate diagnosis. It is important for facilities providing genetic services to keep track of changes in the technology of exome capture in order to maximize throughput while reducing cost per sample. In this study, we focused on comparing the newly released exome probe set Agilent SureSelect Human All Exon v8 and the previous probe set v7. In preparation for higher throughput of exome sequencing using the DNBSEQ-G400, we evaluated target design, coverage statistics, and variants across these two different exome capture products. Although the target size of the v8 design has not changed much compared to the v7 design (35.24 Mb vs 35.8 Mb), the v8 probe design allows you to call more of SNVs (+ 3.06%) and indels (+ 8.49%) with the same number of raw reads per sample on the common target regions (34.84 Mb). Our results suggest that the new Agilent v8 probe set for exome sequencing yields better data quality than the current Agilent v7 set.

System analysis of the sequencing quality of human whole exome samples on BGI NGS platform

Human exome sequencing is a classical method used in most medical genetic applications. The leaders in the field are the manufacturers of enrichment kits based on hybridization of cRNA or cDNA biotinylated probes specific for a genomic region of interest. Recently, the platforms manufactured by the Chinese company MGI Tech have become widespread in Europe and Asia. The reliability and quality of the obtained data are already beyond any doubt. However, only a few kits compatible with these sequencers can be used for such specific tasks as exome sequencing. We developed our own solution for library pre-capture pooling and exome enrichment with Agilent probes. In this work, using a set of the standard benchmark samples from the Platinum Genome collection, we demonstrate that the qualitative and quantitative parameters of our protocol which we called "RSMU_exome" exceed those of the MGI Tech kit. Our protocol allows for identifying more SNV and indels, generates fewer PCR duplicates, enables pooling of more samples in a single enrichment procedure, and requires less raw data to obtain results comparable with the MGI Tech's protocol. The cost of our protocol is also lower than that of MGI Tech's solution.

Consistent count region-copy number variation (CCR-CNV): an expandable and robust tool for clinical diagnosis of copy number variation at the exon level using next-generation sequencing data

PURPOSE

Despite the importance of exonic copy number variations (CNVs) in human genetic diseases, reliable next-generation sequencing-based methods for detecting them are unavailable. We developed an expandable and robust exonic CNV detection tool called consistent count region (CCR)-CNV.

METHODS

In total, about 1000 samples of the truth set were used for validating CCR-CNV. We compared CCR-CNV performance with 2 well-known CNV tools. Finally, to overcome the limitations of CCR-CNV, we devised a combined approach.

RESULTS

The mean sensitivity and specificity of CCR-CNV alone were above 95%, which was superior to that of other CNV tools, such as DECoN and Atlas-CNV. However, low covered region and positive predictive value and high false discovery rate act as obstacles to its use in clinical settings. The combined approach showed much improved performance than CCR-CNV alone.

CONCLUSION

In this study, we present a novel diagnostic tool that allows the identification of exonic CNVs with high confidence using various reagents and clinical next-generation sequencing platforms. We validated this method using the largest multiple ligation-dependent probe amplification-confirmed data set, including sufficient copy normal control data. The approach, combined with existing CNV tools, allows the implementation of CCR-CNV in clinical settings.

Dealing with Pseudogenes in Molecular Diagnostics in the Next Generation Sequencing Era

Presence of pseudogenes is a dreadful issue in next generation sequencing (NGS), because their contamination can interfere with the detection of variants in the genuine gene and generate false positive and false negative variants.In this chapter we focus on issues related to the application of NGS strategies for analysis of genes with pseudogenes in a clinical setting. The degree to which a pseudogene impacts the ability to accurately detect and map variants in its parent gene depends on the degree of similarity (homology) with the parent gene itself. Hereby, target enrichment and mapping strategies are crucial factors to avoid "contaminating" pseudogene sequences. For target enrichment, we describe advantages and disadvantages of PCR- and capture-based strategies. For mapping strategies, we discuss crucial parameters that need to be considered to accurately distinguish sequences of functional genes from pseudogenic sequences. Finally, we discuss some examples of genes associated with Mendelian disorders, for which interesting NGS approaches are described to avoid interference with pseudogene sequences.

Next-Generation Sequencing in Newborn Screening: A Review of Current State

Newborn screening was first introduced at the beginning of the 1960s with the successful implementation of the first phenylketonuria screening programs. Early expansion of the included disorders was slow because each additional disorder screened required a separate test. Subsequently, the technological advancements of biochemical methodology enabled the scaling-up of newborn screening, most notably with the implementation of tandem mass spectrometry. In recent years, we have witnessed a remarkable progression of high-throughput sequencing technologies, which has resulted in a continuous decrease of both cost and time required for genetic analysis. This has enabled more widespread use of the massive multiparallel sequencing. Genomic sequencing is now frequently used in clinical applications, and its implementation in newborn screening has been intensively advocated. The expansion of newborn screening has raised many clinical, ethical, legal, psychological, sociological, and technological concerns over time. This review provides an overview of the current state of next-generation sequencing regarding newborn screening including current recommendations and potential challenges for the use of such technologies in newborn screening.

Pathologists at the Leading Edge of Optimizing the Tumor Tissue Journey for Diagnostic Accuracy and Molecular Testing

OBJECTIVES

Tumor biomarker analyses accompanying immuno-oncology therapies are coupled with a tumor tissue journey aiming to guide tissue procurement and allow for accurate diagnosis and delivery of test results. The engagement of pathologists in the tumor tissue journey is essential because they are able to link the preanalytic requirements of this process with pathologic evaluation and clinical information, ultimately influencing treatment decisions for patients with cancer. The aim of this review is to provide suggestions on how cancer diagnosis and the delivery of molecular test results may be optimized, based on the needs and available resources of institutions, by placing the tumor tissue journey under the leadership of pathologists.

METHODS

Literature searches on PubMed and personal experience provided the necessary material to satisfy the objectives of this review.

RESULTS

Pathologists are usually involved across many steps of the tumor tissue journey and have the requisite knowledge to ensure its efficiency.

CONCLUSIONS

The expansion of oncology diagnostic testing emphasizes the need for pathologists to acquire a leadership role in the multidisciplinary effort to optimize the accuracy, completeness, and delivery of diagnoses guiding personalized treatments.

Ancient DNA from the Asiatic Wild Dog (Cuon alpinus) from Europe

The Asiatic wild dog (Cuon alpinus), restricted today largely to South and Southeast Asia, was widespread throughout Eurasia and even reached North America during the Pleistocene. Like many other species, it suffered from a huge range loss towards the end of the Pleistocene and went extinct in most of its former distribution. The fossil record of the dhole is scattered and the identification of fossils can be complicated by an overlap in size and a high morphological similarity between dholes and other canid species. We generated almost complete mitochondrial genomes for six putative dhole fossils from Europe. By using three lines of evidence, i.e., the number of reads mapping to various canid mitochondrial genomes, the evaluation and quantification of the mapping evenness along the reference genomes and phylogenetic analysis, we were able to identify two out of six samples as dhole, whereas four samples represent wolf fossils. This highlights the contribution genetic data can make when trying to identify the species affiliation of fossil specimens. The ancient dhole sequences are highly divergent when compared to modern dhole sequences, but the scarcity of dhole data for comparison impedes a more extensive analysis.

Urine as a high-quality source of host genomic DNA from wild populations

The ability to generate genomic data from wild animal populations has the potential to give unprecedented insight into the population history and dynamics of species in their natural habitats. However, for many species, it is impossible legally, ethically or logistically to obtain tissue samples of quality sufficient for genomic analyses. In this study we evaluate the success of multiple sources of genetic material (faeces, urine, dentin and dental calculus) and several capture methods (shotgun, whole-genome, exome) in generating genome-scale data in wild eastern chimpanzees (Pan troglodytes schweinfurthii) from Gombe National Park, Tanzania. We found that urine harbours significantly more host DNA than other sources, leading to broader and deeper coverage across the genome. Urine also exhibited a lower rate of allelic dropout. We found exome sequencing to be far more successful than both shotgun sequencing and whole-genome capture at generating usable data from low-quality samples such as faeces and dental calculus. These results highlight urine as a promising and untapped source of DNA that can be noninvasively collected from wild populations of many species.

Novel CRISPR-based sequence specific enrichment methods for target loci and single base mutations

The programmable sequence specificity of CRISPR has found uses in gene editing and diagnostics. This manuscript describes an additional application of CRISPR through a family of novel DNA enrichment technologies. CAMP (CRISPR Associated Multiplexed PCR) and cCAMP (chimeric CRISPR Associated Multiplexed PCR) utilize the sequence specificity of the Cas9/sgRNA complex to target loci for the ligation of a universal adapter that is used for subsequent amplification. cTRACE (chimeric Targeting Rare Alleles with CRISPR-based Enrichment) also applies this method to use Cas9/sgRNA to target loci for the addition of universal adapters, however it has an additional selection for specific mutations through the use of an allele-specific primer. These three methods can produce multiplex PCR that significantly reduces the optimization required for every target. The methods are also not specific to any downstream analytical platform. We additionally will present a mutation specific enrichment technology that is non-amplification based and leaves the DNA in its native state: TRACE (Targeting Rare Alleles with CRISPR-based Enrichment). TRACE utilizes the Cas9/sgRNA complex to sterically protect the ends of targeted sequences from exonuclease activity which digests both the normal variant as well as any off-target sequences.

Detection of Structural Variants in Circulating Cell-Free DNA from Sarcoma Patients Using Next Generation Sequencing

Circulating tumour DNA (ctDNA) analysis using next generation sequencing (NGS) is being implemented in clinical practice for treatment stratification and disease monitoring. However, using ctDNA to detect structural variants, a common occurrence in sarcoma, can be challenging. Here, we use a sarcoma-specific targeted NGS panel to identify translocations and copy number variants in a cohort of 12 tissue specimens and matched circulating cell-free DNA (cfDNA) from soft tissue sarcoma patients, including alveolar rhabdomyosarcoma (n = 2), Ewing's Sarcoma (n = 2), synovial sarcoma (n = 2), extraskeletal myxoid chondrosarcoma (n = 1), clear cell sarcoma (n = 1), undifferentiated round cell sarcoma (n = 1), myxoid liposarcoma (n = 1), alveolar soft part cell sarcoma (n = 1) and dedifferentiated liposarcoma (n = 1). Structural variants were detected in 11/12 (91.6%) and 6/12 (50%) of tissue and plasma samples, respectively. Structural variants were detected in cfDNA at variant allele frequencies >0.2% with an average sequencing depth of 1026×. The results from this cohort show clinical potential for using NGS in ctDNA to aid in the diagnosis and clinical monitoring of sarcomas and warrant additional studies in larger cohorts.

Helpful Criteria When Implementing NGS Panels in Childhood Lymphoblastic Leukemia

The development of Next-Generation Sequencing (NGS) has provided useful diagnostic, prognostic, and therapeutic strategies for individualized management of B-cell precursor acute lymphoblastic leukemia (BCP-ALL) patients. Consequently, NGS is rapidly being established in clinical practice. However, the technology’s complexity, bioinformatics analysis, and the different available options difficult a broad consensus between different laboratories in its daily routine introduction. This collaborative study among Spanish centers was aimed to assess the feasibility, pros, and cons of our customized panel and other commercial alternatives of NGS-targeted approaches. The custom panel was tested in three different sequencing centers. We used the same samples to assess other commercial panels (Oncomine^TM Childhood Cancer Research Assay; Archer^®FusionPlex^® ALL, and Human Comprehensive Cancer Panel GeneRead Panel v2^®). Overall, the panels showed a good performance in different centers and platforms, but each NGS approach presented some issues, as well as pros and cons. Moreover, a previous consensus on the analysis and reporting following international guidelines would be preferable to improve the concordance in results among centers. Our study shows the challenges posed by NGS methodology and the need to consider several aspects of the chosen NGS-targeted approach and reach a consensus before implementing it in daily practice.

A 4.6 Mb Inversion Leading to PCDH15-LINC00844 and BICC1-PCDH15 Fusion Transcripts as a New Pathogenic Mechanism Implicated in Usher Syndrome Type 1

Usher type 1 syndrome is a rare autosomal recessive disorder involving congenital severe-to-profound hearing loss, development of vision impairment in the first decade, and severe balance difficulties. The PCDH15 gene, one of the five genes implicated in this disease, is involved in 8–20% of cases. In this study, we aimed to identify and characterize the two causal variants in a French patient with typical Usher syndrome clinical features. Massively parallel sequencing-based gene panel and screening for large rearrangements were used, which detected a single multi-exon deletion in the PCDH15 gene. As the second pathogenic event was likely localized in the unscreened regions of the gene, PCDH15 transcripts from cultured nasal cells were analyzed and revealed a loss of junction between exon 13 and exon 14. This aberration could be explained by the identification of two fusion transcripts, PCDH15-LINC00844 and BICC1-PCDH15, originating from a 4.6 Mb inversion. This complex chromosomal rearrangement could not be detected by our diagnostic approach but was instead characterized by long-read sequencing, which offers the possibility of detecting balanced structural variants (SVs). This finding extends our knowledge of the mutational spectrum of the PCDH15 gene with the first ever identification of a large causal paracentric inversion of chromosome 10 and illustrates the utility of screening balanced SVs in an exhaustive molecular diagnostic approach.

Shedding light on dark genes: enhanced targeted resequencing by optimizing the combination of enrichment technology and DNA fragment length

The exome contains many obscure regions difficult to explore with current short-read sequencing methods. Repetitious genomic regions prevent the unique alignment of reads, which is essential for the identification of clinically-relevant genetic variants. Long-read technologies attempt to resolve multiple-mapping regions, but they still produce many sequencing errors. Thus, a new approach is required to enlighten the obscure regions of the genome and rescue variants that would be otherwise neglected. This work aims to improve the alignment of multiple-mapping reads through the extension of the standard DNA fragment size. As Illumina can sequence fragments up to 550 bp, we tested different DNA fragment lengths using four major commercial WES platforms and found that longer DNA fragments achieved a higher genotypability. This metric, which indicates base calling calculated by combining depth of coverage with the confidence of read alignment, increased from hundreds to thousands of genes, including several associated with clinical phenotypes. While depth of coverage has been considered crucial for the assessment of WES performance, we demonstrated that genotypability has a greater impact in revealing obscure regions, with ~1% increase in variant calling in respect to shorter DNA fragments. Results confirmed that this approach enlightened many regions previously not explored.

Methodological differences can affect sequencing depth with a possible impact on the accuracy of genetic diagnosis

For a better interpretation of variants, evidence-based databases, such as ClinVar, compile data on the presumed relationships between variants and phenotypes. In this study, we aimed to analyze the pattern of sequencing depth in variants from whole-exome sequencing data in the 1000 Genomes project phase 3, focusing on the variants present in the ClinVar database that were predicted to affect protein-coding regions. We demonstrate that the distribution of the sequencing depth varies across different sequencing centers (pair-wise comparison, p < 0.001). Most importantly, we found that the distribution pattern of sequencing depth is specific to each facility, making it possible to correctly assign 96.9% of the samples to their sequencing center. Thus, indicating the presence of a systematic bias, related to the methods used in the different facilities, which generates significant variations in breadth and depth in whole-exome sequencing data in clinically relevant regions. Our results show that methodological differences, leading to significant heterogeneity in sequencing depth, may potentially influence the accuracy of genetic diagnosis. Furthermore, our findings highlight how it is still challenging to integrate results from different sequencing centers, which may also have an impact on genomic research.

Diagnostic Yield of Next-Generation Sequencing in Patients With Chronic Kidney Disease of Unknown Etiology

Advances in next-generation sequencing (NGS) techniques, including whole exome sequencing, have facilitated cost-effective sequencing of large regions of the genome, enabling the implementation of NGS in clinical practice. Chronic kidney disease (CKD) is a major contributor to global burden of disease and is associated with an increased risk of morbidity and mortality. CKD can be caused by a wide variety of primary renal disorders. In about one in five CKD patients, no primary renal disease diagnosis can be established. Moreover, recent studies indicate that the clinical diagnosis may be incorrect in a substantial number of patients. Both the absence of a diagnosis or an incorrect diagnosis can have therapeutic implications. Genetic testing might increase the diagnostic accuracy in patients with CKD, especially in patients with unknown etiology. The diagnostic utility of NGS has been shown mainly in pediatric CKD cohorts, while emerging data suggest that genetic testing can also be a valuable diagnostic tool in adults with CKD. In addition to its implications for unexplained CKD, NGS can contribute to the diagnostic process in kidney diseases with an atypical presentation, where it may lead to reclassification of the primary renal disease diagnosis. So far, only a few studies have reported on the diagnostic yield of NGS-based techniques in patients with unexplained CKD. Here, we will discuss the potential diagnostic role of gene panels and whole exome sequencing in pediatric and adult patients with unexplained and atypical CKD.

Detection and identification of transgenic events by next generation sequencing combined with enrichment technologies

Next generation sequencing (NGS) is a promising tool for analysing the quality and safety of food and feed products. The detection and identification of genetically modified organisms (GMOs) is complex, as the diversity of transgenic events and types of structural elements introduced in plants continue to increase. In this paper, we show how a strategy that combines enrichment technologies with NGS can be used to detect a large panel of structural elements and partially or completely reconstruct the new sequence inserted into the plant genome in a single analysis, even at low GMO percentages. The strategy of enriching sequences of interest makes the approach applicable even to mixed products, which was not possible before due to insufficient coverage of the different genomes present. This approach is also the first step towards a more complete characterisation of agrifood products in a single analysis.

Clinical sequencing: From raw data to diagnosis with lifetime value

High-throughput sequencing (HTS) has revolutionized genetics by enabling the detection of sequence variants at hitherto unprecedented large scale. Despite these advances, however, there are still remaining challenges in the complete coverage of targeted regions (genes, exome or genome) as well as in HTS data analysis and interpretation. Moreover, it is easy to get overwhelmed by the plethora of available methods and tools for HTS. Here, we review the step-by-step process from the generation of sequence data to molecular diagnosis of Mendelian diseases. Highlighting advantages and limitations, this review addresses the current state of (1) HTS technologies, considering targeted, whole-exome, and whole-genome sequencing on short- and long-read platforms; (2) read alignment, variant calling and interpretation; as well as (3) regulatory issues related to genetic counseling, reimbursement, and data storage.

Development and Application of Duplex Sequencing Strategy for Cell-Free DNA-Based Longitudinal Monitoring of Stage IV Colorectal Cancer

Potential applications of cell-free DNA (cfDNA)-based molecular profiling have used in patients with diverse malignant tumors. However, capturing all cfDNA that originates from tumor cells and identifying true variants present in this minute fraction remain challenges to the widespread application of cfDNA-based liquid biopsies in the clinical setting. In this study, we evaluate a systematic approach and identify key components of wet bench and bioinformatics strategies to address these challenges. We found that concentration of enrichment oligonucleotides, elements of the library preparation, and the structure of adaptors are critical for achieving high enrichment of target regions, retaining variant allele frequencies accurately throughout all involved steps of library preparation, and obtaining high variant coverage. We developed a dual molecular barcode-integrated error elimination strategy to remove sequencing artifacts and a background error correction strategy to distinguish true variants from abundant false-positive variants. We further describe a clinical application of this cfDNA-based duplex sequencing approach that can be used to monitor disease progression in patients with stage IV colorectal cancer. The findings also suggest that cfDNA-based molecular testing observations are highly concordant with observations obtained by traditional imaging methods. Overall, the findings presented in this study have potential implications for early detection of cancer, identification of minimal residual disease, and evaluation of therapeutic responses in patients with cancer.

Preanalytic Variables and Tissue Stewardship for Reliable Next-Generation Sequencing (NGS) Clinical Analysis

An enduring goal of personalized medicine in cancer is the ability to identify patients who are likely to respond to specific therapies. Our growing understanding of the biology and molecular signatures of individual tumor types has facilitated the identification of predictive biomarkers and has led to an increasing number of diagnostic tests to be performed, often as serial and distinct assays on limited tumor specimens. The biomarker diagnostics field has been revolutionized by next-generation sequencing (NGS), which provides a comprehensive overview of the genomic profile of a tumor. Many preanalytic variables can influence the accuracy and reliability of NGS results. Standardization of preanalytic variables is, however, complicated by the plethora of specimen acquisition and processing methods. Variables across the tissue journey, including specimen acquisition, specimen fixation, and sectioning, as well as postfixation processing, such as nucleic acid extraction, library preparation, and choice of sequencing methods, are critical for the reliability of NGS analysis; thus, standardization would be beneficial. In this article, each step in the tissue journey is outlined, with specific focus on preanalytic variables that can influence NGS results. Practical considerations for standardization of these variables are provided to facilitate accurate, reliable, and reproducible NGS-based molecular characterization of tumors, ultimately informing diagnosis and guiding treatment.

MammaPrint and BluePrint Molecular Diagnostics Using Targeted RNA Next-Generation Sequencing Technology

Next-generation DNA sequencing is rapidly becoming an indispensable tool for genome-directed cancer diagnostics, but next-generation RNA sequencing (RNA-seq) is currently not standardly used in clinical diagnostics for expression assessment. However, multigene RNA diagnostic assays are used increasingly in the routine diagnosis of early-stage breast cancer. Two of the most widely used tests are currently available only as a central laboratory service, which limits their clinical use. We evaluated the use of RNA-seq as a decentralized method to perform such tests. The MammaPrint and BluePrint RNA-seq tests were found to be equivalent to the clinically validated microarray tests. The RNA-seq tests were highly reproducible when performed in different locations and were stable over time. The MammaPrint RNA-seq test was clinically validated. Our data demonstrate that RNA-seq can be used as a decentralized platform, yielding results substantially equivalent to results derived from the predicate diagnostic device.

Clinical metagenomics

Clinical metagenomic next-generation sequencing (mNGS), the comprehensive analysis of microbial and host genetic material (DNA and RNA) in samples from patients, is rapidly moving from research to clinical laboratories. This emerging approach is changing how physicians diagnose and treat infectious disease, with applications spanning a wide range of areas, including antimicrobial resistance, the microbiome, human host gene expression (transcriptomics) and oncology. Here, we focus on the challenges of implementing mNGS in the clinical laboratory and address potential solutions for maximizing its impact on patient care and public health.

The Hidden Genomics of Chlamydia trachomatis

The application of whole-genome sequencing has moved us on from sequencing single genomes to defining unravelling population structures in different niches, and at the -species, -serotype or even -genus level, and in local, national and global settings. This has been instrumental in cataloguing and revealing a huge a range of diversity in this bacterium, when at first we thought there was little. Genomics has challenged assumptions, added insight, as well as confusion and glimpses of truths. What is clear is that at a time when we start to realise the extent and nature of the diversity contained within a genus or a species like this, the huge depth of knowledge communities have developed, through cell biology, as well as the new found molecular approaches will be more precious than ever to link genotype to phenotype. Here we detail the technological developments and insights we have seen during the relatively short time since we began to see the hidden genome of Chlamydia trachomatis.

Rational "Error Elimination" Approach to Evaluating Molecular Barcoded Next-Generation Sequencing Data Identifies Low-Frequency Mutations in Hematologic Malignancies

The emergence of highly sensitive molecular diagnostic approaches, such as droplet digital PCR, has allowed the accurate identification of low-frequency variant alleles in clinical specimens; however, the multiplex capabilities of droplet digital PCR for variant detection are inadequate. The incorporation of molecular barcodes or unique IDs into next-generation sequencing libraries through PCR has enabled the detection of low-frequency variant alleles across multiple genomic regions. However, rational library preparation and sequencing data analytic strategies that integrate molecular barcodes have rarely been applied to clinical settings. In this study, we evaluated the parameters that are crucial in the use of molecular barcodes in next-generation sequencing for genotyping clinical specimens from patients with hematologic malignancies. The uniform incorporation of molecular barcodes into DNA templates through PCR was found to be crucial, and the extent of uniformity was governed by multiple interdependent variables. An error elimination strategy was developed for removing sequencing background errors by using molecular barcode sequence information as an alternative to the conventional error correction approach. This approach was successfully used to identify mutations with frequencies as low as 0.15%, and the clonal heterogeneity of hematologic malignancies was revealed. These findings have implications for elucidating heterogeneity and temporal and spatial clonal evolution, evaluating response to therapy, and monitoring relapse in patients with hematologic malignancies.

A novel CRISPR/Cas9 associated technology for sequence-specific nucleic acid enrichment

Massively parallel sequencing technologies have made it possible to generate large quantities of sequence data. However, as research-associated information is transferred into clinical practice, cost and throughput constraints generally require sequence-specific targeted analyses. Therefore, sample enrichment methods have been developed to meet the needs of clinical sequencing applications. However, current amplification and hybrid capture enrichment methods are limited in the contiguous length of sequences for which they are able to enrich. PCR based amplification also loses methylation data and other native DNA features. We have developed a novel technology (Negative Enrichment) where we demonstrate targeting long (>10 kb) genomic regions of interest. We use the specificity of CRISPR-Cas9 single guide RNA (Cas9/sgRNA) complexes to define 5' and 3' termini of sequence-specific loci in genomic DNA, targeting 10 to 36 kb regions. The complexes were found to provide protection from exonucleases, by protecting the targeted sequences from degradation, resulting in enriched, double-strand, non-amplified target sequences suitable for next-generation sequencing library preparation or other downstream analyses.

Understanding the basis of Ehlers-Danlos syndrome in the era of the next-generation sequencing

Ehlers-Danlos syndrome (EDS) is a clinically and genetically heterogeneous group of heritable connective tissue disorders (HCTDs) defined by joint laxity, skin alterations, and joint hypermobility. The latest EDS classification recognized 13 subtypes in which the clinical and genetic phenotypes are often overlapping, making the diagnosis rather difficult and strengthening the importance of the molecular diagnostic confirmation. New genetic techniques such as next-generation sequencing (NGS) gave the opportunity to identify the genetic bases of unresolved EDS types and support clinical counseling. To date, the molecular defects have been identified in 19 genes, mainly in those encoding collagen, its modifying enzymes or other constituents of the extracellular matrix (ECM). In this review we summarize the contribution of NGS technologies to the current knowledge of the genetic background in different EDS subtypes.

Low-cost cross-taxon enrichment of mitochondrial DNA using in-house synthesised RNA probes

Hybridization capture with in-solution oligonucleotide probes has quickly become the preferred method for enriching specific DNA loci from degraded or ancient samples prior to high-throughput sequencing (HTS). Several companies synthesize sets of probes for in-solution hybridization capture, but these commercial reagents are usually expensive. Methods for economical in-house probe synthesis have been described, but they do not directly address one of the major advantages of commercially synthesised probes: that probe sequences matching many species can be synthesised in parallel and pooled. The ability to make “phylogenetically diverse” probes increases the cost-effectiveness of commercial probe sets, as they can be used across multiple projects (or for projects involving multiple species). However, it is labour-intensive to replicate this with in-house methods, as template molecules must first be generated for each species of interest. While it has been observed that probes can be used to enrich for phylogenetically distant targets, the ability of this effect to compensate for the lack of phylogenetically diverse probes in in-house synthesised probe sets has not been tested. In this study, we present a refined protocol for in-house RNA probe synthesis and evaluated the ability of probes generated using this method from a single species to successfully enrich for the target locus in phylogenetically distant species. We demonstrated that probes synthesized using long-range PCR products from a placental mammal mitochondrion (Bison spp.) could be used to enrich for mitochondrial DNA in birds and marsupials (but not plants). Importantly, our results were obtained for approximately a third of the cost of similar commercially available reagents.

Molecular Analysis for Characterizing Transgenic Events

To develop a commercial trait product, a large number of transgenic events are often produced to obtain the event with desired level of expression. It is crucial to develop efficient and sensitive molecular characterization methods to advance events with stable transgene expression, free of vector backbone sequences and without major changes to the native genome caused by transgene insertion. Here, we discuss a variety of analytical tools, including quantitative PCR (qPCR), Southern blot analysis, and various sequencing technologies, which have been widely used to determine the insert copy number, presence/absence of vector backbone sequences, integrity of the T-DNA, and genomic location of the T-DNA insertion. Moreover, since the discovery of RNA interference in 1998 (Fire et al., Nature 391:806-811, 1998), RNAi has emerged as another powerful tool in in the development of a new transgenic trait for insect control. RNAi creates a double-stranded RNA duplex as the active molecule which forms a strong secondary structure, resulting in challenges for detection. In addition to molecular analysis at the DNA level, this chapter describes detection methods of the active molecules (i.e., double-stranded RNA) for RNAi-based traits.

Rapid and highly-specific generation of targeted DNA sequencing libraries enabled by linking capture probes with universal primers

Targeted Next Generation Sequencing (NGS) is being adopted increasingly broadly in many research, commercial and clinical settings. Currently used target capture methods, however, typically require complex and lengthy (sometimes multi-day) workflows that complicates their use in certain applications. In addition, small panels for high sequencing depth applications such as liquid biopsy typically have low on-target rates, resulting in unnecessarily high sequencing cost. We have developed a novel targeted sequencing library preparation method, named Linked Target Capture (LTC), which replaces typical multi-day target capture workflows with a single-day, combined ‘target-capture-PCR’ workflow. This approach uses physically linked capture probes and PCR primers and is expected to work with panel sizes from 100 bp to >10 Mbp. It reduces the time and complexity of the capture workflow, eliminates long hybridization and wash steps and enables rapid library construction and target capture. High on-target read fractions are achievable due to repeated sequence selection in the target-capture-PCR step, thus lowering sequencing cost. We have demonstrated this technology on sample types including cell-free DNA (cfDNA) and formalin-fixed, paraffin-embedded (FFPE) derived DNA, capturing a 35-gene pan-cancer panel, and therein detecting single nucleotide variants, copy number variants, insertions, deletions and gene fusions. With the integration of unique molecular identifiers (UMIs), variants as low as 0.25% abundance were detected, limited by input mass and sequencing depth. Additionally, sequencing libraries were prepared in less than eight hours from extracted DNA to loaded sequencer, demonstrating that LTC holds promise as a broadly applicable tool for rapid, cost-effective and high performance targeted sequencing.

Next Generation Molecular Diagnosis of Hereditary Spastic Paraplegias: An Italian Cross-Sectional Study

Hereditary spastic paraplegia (HSP) refers to a group of genetically heterogeneous neurodegenerative motor neuron disorders characterized by progressive age-dependent loss of corticospinal motor tract function, lower limb spasticity, and weakness. Recent clinical use of next generation sequencing (NGS) methodologies suggests that they facilitate the diagnostic approach to HSP, but the power of NGS as a first-tier diagnostic procedure is unclear. The larger-than-expected genetic heterogeneity-there are over 80 potential disease-associated genes-and frequent overlap with other clinical conditions affecting the motor system make a molecular diagnosis in HSP cumbersome and time consuming. In a single-center, cross-sectional study, spanning 4 years, 239 subjects with a clinical diagnosis of HSP underwent molecular screening of a large set of genes, using two different customized NGS panels. The latest version of our targeted sequencing panel (SpastiSure3.0) comprises 118 genes known to be associated with HSP. Using an in-house validated bioinformatics pipeline and several in silico tools to predict mutation pathogenicity, we obtained a positive diagnostic yield of 29% (70/239), whereas variants of unknown significance (VUS) were found in 86 patients (36%), and 83 cases remained unsolved. This study is among the largest screenings of consecutive HSP index cases enrolled in real-life clinical-diagnostic settings. Its results corroborate NGS as a modern, first-step procedure for molecular diagnosis of HSP. It also disclosed a significant number of new mutations in ultra-rare genes, expanding the clinical spectrum, and genetic landscape of HSP, at least in Italy.

UK circulating strains of human parainfluenza 3: an amplicon based next generation sequencing method and phylogenetic analysis

Background: Human parainfluenza viruses type 3 (HPIV3) are a prominent cause of respiratory infection with a significant impact in both pediatric and transplant patient cohorts.  Currently there is a paucity of whole genome sequence data that would allow for detailed epidemiological and phylogenetic analysis of circulating strains in the UK. Although it is known that HPIV3 peaks annually in the UK, to date there are no whole genome sequences of HPIV3 UK strains available.  Methods: Clinical strains were obtained from HPIV3 positive respiratory patient samples collected between 2011 and 2015.  These were then amplified using an amplicon based method, sequenced on the Illumina platform and assembled using a new robust bioinformatics pipeline. Phylogenetic analysis was carried out in the context of other epidemiological studies and whole genome sequence data currently available with stringent exclusion of significantly culture-adapted strains of HPIV3. Results: In the current paper we have presented twenty full genome sequences of UK circulating strains of HPIV3 and a detailed phylogenetic analysis thereof.  We have analysed the variability along the HPIV3 genome and identified a short hypervariable region in the non-coding segment between the M (matrix) and F (fusion) genes. The epidemiological classifications obtained by using this region and whole genome data were then compared and found to be identical. Conclusions: The majority of HPIV3 strains were observed at different geographical locations and with a wide temporal spread, reflecting the global distribution of HPIV3. Consistent with previous data, a particular subcluster or strain was not identified as specific to the UK, suggesting that a number of genetically diverse strains circulate at any one time. A small hypervariable region in the HPIV3 genome was identified and it was shown that, in the absence of full genome data, this region could be used for epidemiological surveillance of HPIV3.

Detecting Biothreat Agents: From Current Diagnostics to Developing Sensor Technologies

Although a fundamental understanding of the pathogenicity of most biothreat agents has been elucidated and available treatments have increased substantially over the past decades, they still represent a significant public health threat in this age of (bio)terrorism, indiscriminate warfare, pollution, climate change, unchecked population growth, and globalization. The key step to almost all prevention, protection, prophylaxis, post-exposure treatment, and mitigation of any bioagent is early detection. Here, we review available methods for detecting bioagents including pathogenic bacteria and viruses along with their toxins. An introduction placing this subject in the historical context of previous naturally occurring outbreaks and efforts to weaponize selected agents is first provided along with definitions and relevant considerations. An overview of the detection technologies that find use in this endeavor along with how they provide data or transduce signal within a sensing configuration follows. Current "gold" standards for biothreat detection/diagnostics along with a listing of relevant FDA approved in vitro diagnostic devices is then discussed to provide an overview of the current state of the art. Given the 2014 outbreak of Ebola virus in Western Africa and the recent 2016 spread of Zika virus in the Americas, discussion of what constitutes a public health emergency and how new in vitro diagnostic devices are authorized for emergency use in the U.S. are also included. The majority of the Review is then subdivided around the sensing of bacterial, viral, and toxin biothreats with each including an overview of the major agents in that class, a detailed cross-section of different sensing methods in development based on assay format or analytical technique, and some discussion of related microfluidic lab-on-a-chip/point-of-care devices. Finally, an outlook is given on how this field will develop from the perspective of the biosensing technology itself and the new emerging threats they may face.

Modified SureSelectQXT Target Enrichment Protocol for Illumina Multiplexed Sequencing of FFPE Samples

Background

Personalised medicine is nowadays a major objective in oncology. Molecular characterization of tumours through NGS offers the possibility to find possible therapeutic targets in a time- and cost-effective way. However, the low quality and complexity of FFPE DNA samples bring a series of disadvantages for massive parallel sequencing techniques compared to high-quality DNA samples (from blood cells, cell cultures, etc.).

Results

We performed several experiments to understand the behaviour of FFPE DNA samples during the construction of SureSelect^QXT libraries. First, we designed a quality checkpoint for FFPE DNA samples based on the quantification of their amplification capability (qcPCR). We observed that FFPE DNA samples can be classified according to DIN value and qcPCR concentration into unusable, or low-quality (LQ) and good-quality (GQ) DNA. For GQ samples, we increased the amount of input DNA to 150 ng and the digestion time to 30 min, whereas for LQ samples, we used 50 ng of DNA as input but we decreased the digestion time to 1 min. In all cases, we increased the cycles of the pre-hyb PCR to 10 but decreased the cycles of the post-hyb PCR to 8. In addition, we confirmed that using half of the volume of reagents can be beneficial. Finally, in order to obtain better results, we designed a decision flow-chart to achieve a seeding concentration of 12–14 pM for MiSeq Reagent Kit v2.

Conclusions

Our experiments allowed us to unveil the behaviour of low-quality FFPE DNA samples during the construction of SureSelect^QXT libraries. Sequencing results showed that, using our modified SureSelect^QXT protocol, the final percentage of usable reads for low-quality samples was increased more than three times allowing to reach median depth/million reads values of 76.35. This value is equivalent to ~ 0.9 and ~ 0.7 of the values obtained for good-quality FFPE and high-quality DNA respectively.

Electronic supplementary material

The online version of this article (10.1186/s12575-018-0084-7) contains supplementary material, which is available to authorized users.

A comprehensive targeted next-generation sequencing panel for genetic diagnosis of patients with suspected inherited thrombocytopenia

BACKGROUND

Inherited thrombocytopenias (ITs) are a heterogeneous group of disorders characterized by low platelet counts and often disproportionate bleeding with over 30 genes currently implicated. Previously the UK-GAPP study using whole exome sequencing (WES) identified a pathogenic variant in 19 of 47 (40%) patients of which 71% had variants in genes known to cause IT.

AIMS

To employ a targeted next-generation sequencing platform to improve efficiency of diagnostic testing and reduce overall costs.

METHODS

We have developed an IT-specific gene panel as a pre-screen for patients prior to WES using the Agilent SureSelectQXT transposon-based enrichment system.

RESULTS

Thirty-one patients were analyzed using the panel-based sequencing, of which; 10% (3/31) were identified with a classified pathogenic variant, 16% (5/31) were identified with a likely pathogenic variant, 51% (16/31) were identified with variants of unknown significance, and 23% (7/31) were identified with either no variant or a benign variant.

DISCUSSION AND CONCLUSION

Although requiring further clarification of the impact of the genetic variations, the application of an IT-specific next generation sequencing panel is an viable method of pre-screening patients for variants in known IT-causing genes prior to WES. With an added benefit of distinguishing IT from idiopathic thrombocytopenic purpura (ITP) and the potential to identify variants in genes known to have a predisposition to hematological malignancies, it could become a critical step in improving patient clinical management.

UK circulating strains of human parainfluenza 3: an amplicon based next generation sequencing method and phylogenetic analysis

Background: Human parainfluenza viruses type 3 (HPIV3) are a prominent cause of respiratory infection with a significant impact in both pediatric and transplant patient cohorts.  Currently there is a paucity of whole genome sequence data that would allow for detailed epidemiological and phylogenetic analysis of circulating strains in the UK. Although it is known that HPIV3 peaks annually in the UK, to date there are no whole genome sequences of HPIV3 UK strains available.  Methods: Clinical strains were obtained from HPIV3 positive respiratory patient samples collected between 2011 and 2015.  These were then amplified using an amplicon based method, sequenced on the Illumina platform and assembled using a new robust bioinformatics pipeline. Phylogenetic analysis was carried out in the context of other epidemiological studies and whole genome sequence data currently available with stringent exclusion of significantly culture-adapted strains of HPIV3. Results: In the current paper we have presented twenty full genome sequences of UK circulating strains of HPIV3 and a detailed phylogenetic analysis thereof.  We have analysed the variability along the HPIV3 genome and identified a short hypervariable region in the non-coding segment between the M (matrix) and F (fusion) genes. The epidemiological classifications obtained by using this region and whole genome data were then compared and found to be identical. Conclusions: The majority of HPIV3 strains were observed at different geographical locations and with a wide temporal spread, reflecting the global distribution of HPIV3. Consistent with previous data, a particular subcluster or strain was not identified as specific to the UK, suggesting that a number of genetically diverse strains circulate at any one time. A small hypervariable region in the HPIV3 genome was identified and it was shown that, in the absence of full genome data, this region could be used for epidemiological surveillance of HPIV3.

Design and validation of a universal influenza virus enrichment probe set and its utility in deep sequence analysis of primary cloacal swab surveillance samples of wild birds

Influenza virus infections in humans and animals are major public health concerns. In the current study, a set of universal influenza enrichment probes was developed to increase the sensitivity of sequence-based virus detection and characterization for all influenza viruses. This universal influenza enrichment probe set contains 46,953 120nt RNA biotin-labeled probes designed based on all available influenza viral sequences and it can be used to enrich for influenza sequences without prior knowledge of type or subtype. Marked enrichment was demonstrated in influenza A/H1N1, influenza B, and H1-to-H16 hemagglutinin plasmids spiked into human DNA and in cultured influenza A/H2N1 virus. Furthermore, enrichment effects and mixed influenza A virus infections were revealed in wild bird cloacal swab samples. Therefore, this universal influenza virus enrichment probe system can capture and enrich influenza viral sequences selectively and effectively in different samples, especially ones with degraded RNA or containing low amount of influenza RNA.

Targeted Sequencing of Respiratory Viruses in Clinical Specimens for Pathogen Identification and Genome-Wide Analysis

A large number of viruses can individually and concurrently cause various respiratory illnesses. Metagenomic sequencing using next-generation sequencing (NGS) technology is capable of identifying a variety of pathogens. Here, we describe a method using a large panel of oligo probes to enrich sequence targets of 34 respiratory DNA and RNA viruses that reduces non-viral reads in NGS data and achieves high performance of sequencing-based pathogen identification. The approach can be applied to total nucleic acids purified from respiratory swabs stored in viral transport medium. Illumina TruSeq RNA Access Library procedure is used in targeted sequencing of respiratory viruses. The samples are subjected to RNA fragmentation, random reverse transcription, random PCR amplification, and ligation with barcoded library adaptors. The libraries are pooled and subjected to two rounds of enrichments by using a large panel of oligos designed to capture whole genomes of 34 respiratory viruses. The enriched libraries are amplified and sequenced using Illumina MiSeq sequencing system and reagents. This method can achieve viral detection sensitivity comparable with molecular assay and obtain partial to complete genome sequences for each virus to allow accurate genotyping and variant analysis.

CoverView: a sequence quality evaluation tool for next generation sequencing data

Quality assurance and quality control are essential for robust next generation sequencing (NGS). Here we present CoverView, a fast, flexible, user-friendly quality evaluation tool for NGS data. CoverView processes mapped sequencing reads and user-specified regions to report depth of coverage, base and mapping quality metrics with increasing levels of detail from a chromosome-level summary to per-base profiles. CoverView can flag regions that do not fulfil user-specified quality requirements, allowing suboptimal data to be systematically and automatically presented for review. It also provides an interactive graphical user interface (GUI) that can be opened in a web browser and allows intuitive exploration of results.

We have integrated CoverView into our accredited clinical cancer predisposition gene testing laboratory that uses the TruSight Cancer Panel (TSCP). CoverView has been invaluable for optimisation and quality control of our testing pipeline, providing transparent, consistent quality metric information and automatic flagging of regions that fall below quality thresholds. We demonstrate this utility with TSCP data from the Genome in a Bottle reference sample, which CoverView analysed in 13 seconds.

CoverView uses data routinely generated by NGS pipelines, reads standard input formats, and rapidly creates easy-to-parse output text (.txt) files that are customised by a simple configuration file. CoverView can therefore be easily integrated into any NGS pipeline. CoverView and detailed documentation for its use are freely available at github.com/RahmanTeamDevelopment/CoverView/releases and www.icr.ac.uk/CoverView

Homology-based enzymatic DNA fragment assembly-based illumina sequencing library preparation

The current Illumina HiSeq and MiSeq platforms can generate paired-end reads of up to 2 x 250 bp and 2 x 300 bp in length, respectively. These read lengths may be substantially longer than genomic regions of interest when a DNA sequencing library is prepared through a target enrichment-based approach. A sequencing library preparation method has been developed based on the homology-based enzymatic DNA fragment assembly scheme to allow processing of multiple PCR products within a single read. Target sequences were amplified using locus-specific PCR primers with 8 bp tags, and using the tags, homology-based enzymatic DNA assembly was performed with DNA polymerase, T7 exonuclease and T4 DNA ligase. Short PCR amplicons can hence be assembled into a single molecule, along with sequencing adapters specific to the Illumina platforms. As a proof-of-concept experiment, short PCR amplicons (57-66 bp in length) derived from genomic DNA templates of field pea and containing variable nucleotide locations were assembled and sequenced on the MiSeq platform. The results were validated with other genotyping methods. When 5 PCR amplicons were assembled, 4.3 targeted sequences (single-nucleotide polymorphisms) on average were successfully identified within each read. The utility of this for sequencing of short fragments has consequently been demonstrated.

Computational Pipeline for the PGV-001 Neoantigen Vaccine Trial

This paper describes the sequencing protocol and computational pipeline for the PGV-001 personalized vaccine trial. PGV-001 is a therapeutic peptide vaccine targeting neoantigens identified from patient tumor samples. Peptides are selected by a computational pipeline that identifies mutations from tumor/normal exome sequencing and ranks mutant sequences by a combination of predicted Class I MHC affinity and abundance estimated from tumor RNA. The personalized genomic vaccine (PGV) pipeline is modular and consists of independently usable tools and software libraries. We hope that the functionality of these tools may extend beyond the specifics of the PGV-001 trial and enable other research groups in their own neoantigen investigations.

A robust targeted sequencing approach for low input and variable quality DNA from clinical samples

Next-generation deep sequencing of gene panels is being adopted as a diagnostic test to identify actionable mutations in cancer patient samples. However, clinical samples, such as formalin-fixed, paraffin-embedded specimens, frequently provide low quantities of degraded, poor quality DNA. To overcome these issues, many sequencing assays rely on extensive PCR amplification leading to an accumulation of bias and artifacts. Thus, there is a need for a targeted sequencing assay that performs well with DNA of low quality and quantity without relying on extensive PCR amplification. We evaluate the performance of a targeted sequencing assay based on Oligonucleotide Selective Sequencing, which permits the enrichment of genes and regions of interest and the identification of sequence variants from low amounts of damaged DNA. This assay utilizes a repair process adapted to clinical FFPE samples, followed by adaptor ligation to single stranded DNA and a primer-based capture technique. Our approach generates sequence libraries of high fidelity with reduced reliance on extensive PCR amplification—this facilitates the accurate assessment of copy number alterations in addition to delivering accurate single nucleotide variant and insertion/deletion detection. We apply this method to capture and sequence the exons of a panel of 130 cancer-related genes, from which we obtain high read coverage uniformity across the targeted regions at starting input DNA amounts as low as 10 ng per sample. We demonstrate the performance using a series of reference DNA samples, and by identifying sequence variants in DNA from matched clinical samples originating from different tissue types.

A new DNA sequencing technology enables comprehensive genetic analyses of poor-quality tumor samples. Hanlee Ji from Stanford University in California, USA, together with colleagues from a company he cofounded called TOMA Biosciences, tested the performance of a targeted sequencing assay known as oligonucleotide-selective sequencing (OS-Seq). They used the “in-solution” version of OS-Seq, which involves a pre-processing step to remove any damaged DNA and then sequences target regions of the genome to look for duplications, insertions or deletions of DNA segments. Using archival specimens (which often contain low quantities of degraded DNA) from patients with lung and colorectal cancer, the researchers showed they could detect sequence variants in a panel of 130 cancer-related genes. The findings suggest the OS-Seq assay could help inform treatment decisions for cancer patients, even with clinical specimens of low quality.

Enrichment of megabase-sized DNA molecules for single-molecule optical mapping and next-generation sequencing

Next-generation sequencing (NGS) has caused a revolution, yet left a gap: long-range genetic information from native, non-amplified DNA fragments is unavailable. It might be obtained by optical mapping of megabase-sized DNA molecules. Frequently only a specific genomic region is of interest, so here we introduce a method for selection and enrichment of megabase-sized DNA molecules intended for single-molecule optical mapping: DNA from a human cell line is digested by the NotI rare-cutting enzyme and size-selected by pulsed-field gel electrophoresis. For demonstration, more than 600 sub-megabase- to megabase-sized DNA molecules were recovered from the gel and analysed by denaturation-renaturation optical mapping. Size-selected molecules from the same gel were sequenced by NGS. The optically mapped molecules and the NGS reads showed enrichment from regions defined by NotI restriction sites. We demonstrate that the unannotated genome can be characterized in a locus-specific manner via molecules partially overlapping with the annotated genome. The method is a promising tool for investigation of structural variants in enriched human genomic regions for both research and diagnostic purposes. Our enrichment method could potentially work with other genomes or target specified regions by applying other genomic editing tools, such as the CRISPR/Cas9 system.

High-throughput annotation of full-length long noncoding RNAs with capture long-read sequencing

Accurate annotation of genes and their transcripts is a foundation of genomics, but currently no annotation technique combines throughput and accuracy. As a result, reference gene collections remain incomplete-many gene models are fragmentary, and thousands more remain uncataloged, particularly for long noncoding RNAs (lncRNAs). To accelerate lncRNA annotation, the GENCODE consortium has developed RNA Capture Long Seq (CLS), which combines targeted RNA capture with third-generation long-read sequencing. Here we present an experimental reannotation of the GENCODE intergenic lncRNA populations in matched human and mouse tissues that resulted in novel transcript models for 3,574 and 561 gene loci, respectively. CLS approximately doubled the annotated complexity of targeted loci, outperforming existing short-read techniques. Full-length transcript models produced by CLS enabled us to definitively characterize the genomic features of lncRNAs, including promoter and gene structure, and protein-coding potential. Thus, CLS removes a long-standing bottleneck in transcriptome annotation and generates manual-quality full-length transcript models at high-throughput scales.

Combined genetic approaches yield a 48% diagnostic rate in a large cohort of French hearing-impaired patients

Hearing loss is the most common sensory disorder and because of its high genetic heterogeneity, implementation of Massively Parallel Sequencing (MPS) in diagnostic laboratories is greatly improving the possibilities of offering optimal care to patients. We present the results of a two-year period of molecular diagnosis that included 207 French families referred for non-syndromic hearing loss. Our multi-step strategy involved (i) DFNB1 locus analysis, (ii) MPS of 74 genes, and (iii) additional approaches including Copy Number Variations, in silico analyses, minigene studies coupled when appropriate with complete gene sequencing, and a specific assay for STRC. This comprehensive screening yielded an overall diagnostic rate of 48%, equally distributed between DFNB1 (24%) and the other genes (24%). Pathogenic genotypes were identified in 19 different genes, with a high prevalence of GJB2, STRC, MYO15A, OTOF, TMC1, MYO7A and USH2A. Involvement of an Usher gene was reported in 16% of the genotyped cohort. Four de novo variants were identified. This study highlights the need to develop several molecular approaches for efficient molecular diagnosis of hearing loss, as this is crucial for genetic counselling, audiological rehabilitation and the detection of syndromic forms.