Cross-oncopanel study reveals high sensitivity and accuracy with overall analytical performance depending on genomic regions.

BACKGROUND: Targeted sequencing using oncopanels requires comprehensive assessments of accuracy and detection sensitivity to ensure analytical validity. By employing reference materials characterized by the U.S. Food and Drug Administration-led SEquence Quality Control project phase2 (SEQC2) effort, we perform a cross-platform multi-lab evaluation of eight Pan-Cancer panels to assess best practices for oncopanel sequencing. RESULTS: All panels demonstrate high sensitivity across targeted high-confidence coding regions and variant types for the variants previously verified to have variant allele frequency (VAF) in the 5-20% range. Sensitivity is reduced by utilizing VAF thresholds due to inherent variability in VAF measurements. Enforcing a VAF threshold for reporting has a positive impact on reducing false positive calls. Importantly, the false positive rate is found to be significantly higher outside the high-confidence coding regions, resulting in lower reproducibility. Thus, region restriction and VAF thresholds lead to low relative technical variability in estimating promising biomarkers and tumor mutational burden. CONCLUSION: This comprehensive study provides actionable guidelines for oncopanel sequencing and clear evidence that supports a simplified approach to assess the analytical performance of oncopanels. It will facilitate the rapid implementation, validation, and quality control of oncopanels in clinical use.

A verified genomic reference sample for assessing performance of cancer panels detecting small variants of low allele frequency.

BACKGROUND: Oncopanel genomic testing, which identifies important somatic variants, is increasingly common in medical practice and especially in clinical trials. Currently, there is a paucity of reliable genomic reference samples having a suitably large number of pre-identified variants for properly assessing oncopanel assay analytical quality and performance. The FDA-led Sequencing and Quality Control Phase 2 (SEQC2) consortium analyze ten diverse cancer cell lines individually and their pool, termed Sample A, to develop a reference sample with suitably large numbers of coding positions with known (variant) positives and negatives for properly evaluating oncopanel analytical performance. RESULTS: In reference Sample A, we identify more than 40,000 variants down to 1% allele frequency with more than 25,000 variants having less than 20% allele frequency with 1653 variants in COSMIC-related genes. This is 5-100× more than existing commercially available samples. We also identify an unprecedented number of negative positions in coding regions, allowing statistical rigor in assessing limit-of-detection, sensitivity, and precision. Over 300 loci are randomly selected and independently verified via droplet digital PCR with 100% concordance. Agilent normal reference Sample B can be admixed with Sample A to create new samples with a similar number of known variants at much lower allele frequency than what exists in Sample A natively, including known variants having allele frequency of 0.02%, a range suitable for assessing liquid biopsy panels. CONCLUSION: These new reference samples and their admixtures provide superior capability for performing oncopanel quality control, analytical accuracy, and validation for small to large oncopanels and liquid biopsy assays.

Evaluating the analytical validity of circulating tumor DNA sequencing assays for precision oncology.

Circulating tumor DNA (ctDNA) sequencing is being rapidly adopted in precision oncology, but the accuracy, sensitivity and reproducibility of ctDNA assays is poorly understood. Here we report the findings of a multi-site, cross-platform evaluation of the analytical performance of five industry-leading ctDNA assays. We evaluated each stage of the ctDNA sequencing workflow with simulations, synthetic DNA spike-in experiments and proficiency testing on standardized, cell-line-derived reference samples. Above 0.5% variant allele frequency, ctDNA mutations were detected with high sensitivity, precision and reproducibility by all five assays, whereas, below this limit, detection became unreliable and varied widely between assays, especially when input material was limited. Missed mutations (false negatives) were more common than erroneous candidates (false positives), indicating that the reliable sampling of rare ctDNA fragments is the key challenge for ctDNA assays. This comprehensive evaluation of the analytical performance of ctDNA assays serves to inform best practice guidelines and provides a resource for precision oncology.

De novo assembly of the olive fruit fly (Bactrocera oleae) genome with linked-reads and long-read technologies minimizes gaps and provides exceptional Y chromosome assembly.

BACKGROUND: The olive fruit fly, Bactrocera oleae, is the most important pest in the olive fruit agribusiness industry. This is because female flies lay their eggs in the unripe fruits and upon hatching the larvae feed on the fruits thus destroying them. The lack of a high-quality genome and other genomic and transcriptomic data has hindered progress in understanding the fly's biology and proposing alternative control methods to pesticide use. RESULTS: Genomic DNA was sequenced from male and female Demokritos strain flies, maintained in the laboratory for over 45 years. We used short-, mate-pair-, and long-read sequencing technologies to generate a combined male-female genome assembly (GenBank accession GCA_001188975.2). Genomic DNA sequencing from male insects using 10x Genomics linked-reads technology followed by mate-pair and long-read scaffolding and gap-closing generated a highly contiguous 489 Mb genome with a scaffold N50 of 4.69 Mb and L50 of 30 scaffolds (GenBank accession GCA_001188975.4). RNA-seq data generated from 12 tissues and/or developmental stages allowed for genome annotation. Short reads from both males and females and the chromosome quotient method enabled identification of Y-chromosome scaffolds which were extensively validated by PCR. CONCLUSIONS: The high-quality genome generated represents a critical tool in olive fruit fly research. We provide an extensive RNA-seq data set, and genome annotation, critical towards gaining an insight into the biology of the olive fruit fly. In addition, elucidation of Y-chromosome sequences will advance our understanding of the Y-chromosome's organization, function and evolution and is poised to provide avenues for sterile insect technique approaches.

Unique, dual-indexed sequencing adapters with UMIs effectively eliminate index cross-talk and significantly improve sensitivity of massively parallel sequencing.

BACKGROUND: Sample index cross-talk can result in false positive calls when massively parallel sequencing (MPS) is used for sensitive applications such as low-frequency somatic variant discovery, ancient DNA investigations, microbial detection in human samples, or circulating cell-free tumor DNA (ctDNA) variant detection. Therefore, the limit-of-detection of an MPS assay is directly related to the degree of index cross-talk. RESULTS: Cross-talk rates up to 0.29% were observed when using standard, combinatorial adapters, resulting in 110,180 (0.1% cross-talk rate) or 1,121,074 (0.29% cross-talk rate) misassigned reads per lane in non-patterned and patterned Illumina flow cells, respectively. Here, we demonstrate that using unique, dual-matched indexed adapters dramatically reduces index cross-talk to ≤1 misassigned reads per flow cell lane. While the current study was performed using dual-matched indices, using unique, dual-unrelated indices would also be an effective alternative. CONCLUSIONS: For sensitive downstream analyses, the use of combinatorial indices for multiplexed hybrid capture and sequencing is inappropriate, as it results in an unacceptable number of misassigned reads. Cross-talk can be virtually eliminated using dual-matched indexed adapters. These results suggest that use of such adapters is critical to reduce false positive rates in assays that aim to identify low allele frequency events, and strongly indicate that dual-matched adapters be implemented for all sensitive MPS applications.

Extensive sequencing of seven human genomes to characterize benchmark reference materials.

The Genome in a Bottle Consortium, hosted by the National Institute of Standards and Technology (NIST) is creating reference materials and data for human genome sequencing, as well as methods for genome comparison and benchmarking. Here, we describe a large, diverse set of sequencing data for seven human genomes; five are current or candidate NIST Reference Materials. The pilot genome, NA12878, has been released as NIST RM 8398. We also describe data from two Personal Genome Project trios, one of Ashkenazim Jewish ancestry and one of Chinese ancestry. The data come from 12 technologies: BioNano Genomics, Complete Genomics paired-end and LFR, Ion Proton exome, Oxford Nanopore, Pacific Biosciences, SOLiD, 10X Genomics GemCode WGS, and Illumina exome and WGS paired-end, mate-pair, and synthetic long reads. Cell lines, DNA, and data from these individuals are publicly available. Therefore, we expect these data to be useful for revealing novel information about the human genome and improving sequencing technologies, SNP, indel, and structural variant calling, and de novo assembly.

A hybrid approach for de novo human genome sequence assembly and phasing.

Despite tremendous progress in genome sequencing, the basic goal of producing a phased (haplotype-resolved) genome sequence with end-to-end contiguity for each chromosome at reasonable cost and effort is still unrealized. In this study, we describe an approach to performing de novo genome assembly and experimental phasing by integrating the data from Illumina short-read sequencing, 10X Genomics linked-read sequencing, and BioNano Genomics genome mapping to yield a high-quality, phased, de novo assembled human genome.

Health and population effects of rare gene knockouts in adult humans with related parents.

Examining complete gene knockouts within a viable organism can inform on gene function. We sequenced the exomes of 3222 British adults of Pakistani heritage with high parental relatedness, discovering 1111 rare-variant homozygous genotypes with predicted loss of function (knockouts) in 781 genes. We observed 13.7% fewer homozygous knockout genotypes than we expected, implying an average load of 1.6 recessive-lethal-equivalent loss-of-function (LOF) variants per adult. When genetic data were linked to the individuals' lifelong health records, we observed no significant relationship between gene knockouts and clinical consultation or prescription rate. In this data set, we identified a healthy PRDM9-knockout mother and performed phased genome sequencing on her, her child, and control individuals. Our results show that meiotic recombination sites are localized away from PRDM9-dependent hotspots. Thus, natural LOF variants inform on essential genetic loci and demonstrate PRDM9 redundancy in humans.

Haplotyping germline and cancer genomes with high-throughput linked-read sequencing.

Haplotyping of human chromosomes is a prerequisite for cataloguing the full repertoire of genetic variation. We present a microfluidics-based, linked-read sequencing technology that can phase and haplotype germline and cancer genomes using nanograms of input DNA. This high-throughput platform prepares barcoded libraries for short-read sequencing and computationally reconstructs long-range haplotype and structural variant information. We generate haplotype blocks in a nuclear trio that are concordant with expected inheritance patterns and phase a set of structural variants. We also resolve the structure of the EML4-ALK gene fusion in the NCI-H2228 cancer cell line using phased exome sequencing. Finally, we assign genetic aberrations to specific megabase-scale haplotypes generated from whole-genome sequencing of a primary colorectal adenocarcinoma. This approach resolves haplotype information using up to 100 times less genomic DNA than some methods and enables the accurate detection of structural variants.

TMTC1 and TMTC2 are novel endoplasmic reticulum tetratricopeptide repeat-containing adapter proteins involved in calcium homeostasis.

The endoplasmic reticulum (ER) is organized in part by adapter proteins that nucleate the formation of large protein complexes. Tetratricopeptide repeats (TPR) are well studied protein structural motifs that support intermolecular protein-protein interactions. TMTC1 and TMTC2 were identified by an in silico search as TPR-containing proteins possessing N-terminal ER targeting signal sequences and multiple hydrophobic segments, suggestive of polytopic membrane proteins that are targeted to the secretory pathway. A variety of cell biological and biochemical assays was employed to demonstrate that TMTC1 and TMTC2 are both ER resident integral membrane proteins with multiple clusters of TPR domains oriented within the ER lumen. Proteomic analysis followed by co-immunoprecipitation verification found that both proteins associated with the ER calcium uptake pump SERCA2B, and TMTC2 also bound to the carbohydrate-binding chaperone calnexin. Live cell calcium measurements revealed that overexpression of either TMTC1 or TMTC2 caused a reduction of calcium released from the ER following stimulation, whereas the knockdown of TMTC1 or TMTC2 increased the stimulated calcium released. Together, these results implicate TMTC1 and TMTC2 as ER proteins involved in ER calcium homeostasis.

Viroporins customize host cells for efficient viral propagation.

Viruses are intracellular parasites that must access the host cell machinery to propagate. Viruses hijack the host cell machinery to help with entry, replication, packaging, and release of progeny to infect new cells. To carry out these diverse functions, viruses often transform the cellular environment using viroporins, a growing class of viral-encoded membrane proteins that promote viral proliferation. Viroporins modify the integrity of host membranes, thereby stimulating the maturation of viral infection, and are critical for virus production and dissemination. Significant advances in molecular and cell biological approaches have helped to uncover some of the roles that viroporins serve in the various stages of the viral life cycle. In this study, the ability of viroporins to modify the cellular environment will be discussed, with particular emphasis on their role in the stepwise progression of the viral life cycle.

SV40 late protein VP4 forms toroidal pores to disrupt membranes for viral release.

Nonenveloped viruses are generally released from the cell by the timely lysis of host cell membranes. SV40 has been used as a model virus for the study of the lytic nonenveloped virus life cycle. The expression of SV40 VP4 at later times during infection is concomitant with cell lysis. To investigate the role of VP4 in viral release and its mechanism of action, VP4 was expressed and purified from bacteria as a fusion protein for use in membrane disruption assays. Purified VP4 perforated membranes as demonstrated by the release of fluorescent markers encapsulated within large unilamellar vesicles or liposomes. Dynamic light scattering results revealed that VP4 treatment did not cause membrane lysis or change the size of the liposomes. Liposomes encapsulated with 4,4-difluoro-5,7-dimethyl-4-bora-3a,4a-diaza-3-indacene-labeled streptavidin were used to show that VP4 formed stable pores in membranes. These VP4 pores had an inner diameter of 1-5 nm. Asymmetrical liposomes containing pyrene-labeled lipids in the outer monolayer were employed to monitor transbilayer lipid diffusion. Consistent with VP4 forming toroidal pore structures in membranes, VP4 induced transbilayer lipid diffusion or lipid flip-flop. Altogether, these studies support a central role for VP4 acting as a viroporin in the disruption of cellular membranes to trigger SV40 viral release by forming toroidal pores that unite the outer and inner leaflets of membrane bilayers.

The viroporin activity of the minor structural proteins VP2 and VP3 is required for SV40 propagation.

For nonenveloped viruses such as Simian Virus 40, the mechanism used to translocate viral components across membranes is poorly understood. Previous results indicated that the minor structural proteins, VP2 and VP3, might act as membrane proteins during infection. Here, purified VP2 and VP3 were found to form pores in host cell membranes. To identify possible membrane domains, individual hydrophobic domains from VP2 and VP3 were expressed in a model protein and tested for their ability to integrate into membranes. Several domains from the late proteins supported endoplasmic reticulum membrane insertion as transmembrane domains. Mutations in VP2 and VP3 were engineered that inhibited membrane insertion and pore formation. When these mutations were introduced into the viral genome, viral propagation was inhibited. This comprehensive approach revealed that the viroporin activity of VP2 and VP3 was inhibited by targeted disruptions of individual hydrophobic domains and the loss of membrane disruption activity impaired viral infection.

The Simian virus 40 late viral protein VP4 disrupts the nuclear envelope for viral release.

Simian virus 40 (SV40) appears to initiate cell lysis by expressing the late viral protein VP4 at the end of infection to aid in virus dissemination. To investigate the contribution of VP4 to cell lysis, VP4 was expressed in mammalian cells where it was predominantly observed along the nuclear periphery. The integrity of the nuclear envelope was compromised in these cells, resulting in the mislocalization of a soluble nuclear marker. Using assays that involved the cellular expression of VP4 or the treatment of cells with purified VP4, we found that the central hydrophobic domain and a proximal C-terminal nuclear localization signal of VP4 were required for (i) cytolysis associated with prolonged expression; (ii) nuclear envelope accumulation; and (iii) disruption of the nuclear, red blood cell, or host cell membranes. Furthermore, a conserved proline within the hydrophobic domain was required for membrane perforation, suggesting that this residue was crucial for VP4 cytolytic activity. These results indicate that VP4 forms pores in the nuclear membrane leading to lysis and virus release.

The SV40 late protein VP4 is a viroporin that forms pores to disrupt membranes for viral release.

Nonenveloped viruses are generally released by the timely lysis of the host cell by a poorly understood process. For the nonenveloped virus SV40, virions assemble in the nucleus and then must be released from the host cell without being encapsulated by cellular membranes. This process appears to involve the well-controlled insertion of viral proteins into host cellular membranes rendering them permeable to large molecules. VP4 is a newly identified SV40 gene product that is expressed at late times during the viral life cycle that corresponds to the time of cell lysis. To investigate the role of this late expressed protein in viral release, water-soluble VP4 was expressed and purified as a GST fusion protein from bacteria. Purified VP4 was found to efficiently bind biological membranes and support their disruption. VP4 perforated membranes by directly interacting with the membrane bilayer as demonstrated by flotation assays and the release of fluorescent markers encapsulated into large unilamellar vesicles or liposomes. The central hydrophobic domain of VP4 was essential for membrane binding and disruption. VP4 displayed a preference for membranes comprised of lipids that replicated the composition of the plasma membranes over that of nuclear membranes. Phosphatidylethanolamine, a lipid found at high levels in bacterial membranes, was inhibitory against the membrane perforation activity of VP4. The disruption of membranes by VP4 involved the formation of pores of ∼3 nm inner diameter in mammalian cells including permissive SV40 host cells. Altogether, these results support a central role of VP4 acting as a viroporin in the perforation of cellular membranes to trigger SV40 viral release.