The impact of DNA input amount and DNA source on the performance of whole-exome sequencing in cancer epidemiology

Validated WGS and WES protocols proved saliva-derived gDNA as an equivalent to blood-derived gDNA for clinical and population genomic analyses

BACKGROUND

Whole exome sequencing (WES) and whole genome sequencing (WGS) have become standard methods in human clinical diagnostics as well as in population genomics (POPGEN). Blood-derived genomic DNA (gDNA) is routinely used in the clinical environment. Conversely, many POPGEN studies and commercial tests benefit from easy saliva sampling. Here, we evaluated the quality of variant call sets and the level of genotype concordance of single nucleotide variants (SNVs) and small insertions and deletions (indels) for WES and WGS using paired blood- and saliva-derived gDNA isolates employing genomic reference-based validated protocols.

METHODS

The genomic reference standard Coriell NA12878 was repeatedly analyzed using optimized WES and WGS protocols, and data calls were compared with the truth dataset published by the Genome in a Bottle Consortium. gDNA was extracted from the paired blood and saliva samples of 10 participants and processed using the same protocols. A comparison of paired blood-saliva call sets was performed in the context of WGS and WES genomic reference-based technical validation results.

RESULTS

The quality pattern of called variants obtained from genomic-reference-based technical replicates correlates with data calls of paired blood-saliva-derived samples in all levels of tested examinations despite a higher rate of non-human contamination found in the saliva samples. The F1 score of 10 blood-to-saliva-derived comparisons ranged between 0.8030-0.9998 for SNVs and between 0.8883-0.9991 for small-indels in the case of the WGS protocol, and between 0.8643-0.999 for SNVs and between 0.7781-1.000 for small-indels in the case of the WES protocol.

CONCLUSION

Saliva may be considered an equivalent material to blood for genetic analysis for both WGS and WES under strict protocol conditions. The accuracy of sequencing metrics and variant-detection accuracy is not affected by choosing saliva as the gDNA source instead of blood but much more significantly by the genomic context, variant types, and the sequencing technology used.

A polygenic score associated with fracture risk in breast cancer patients treated with aromatase inhibitors

Identifying women at high risk of osteoporotic fracture from aromatase inhibitor (AI) therapy for breast cancer is largely based on known risk factors for healthy postmenopausal women, which might not accurately reflect the risk in breast cancer patients post-AI therapy. To determine whether a polygenic score associated with fracture in healthy women is also significant in women treated with AIs for breast cancer, we used data from a prospective observational cohort of 2152 women diagnosed with hormonal receptor positive breast cancer treated with AIs as the initial endocrine therapy and examined a polygenic score of heel quantitative ultrasound speed of sound (gSOS) in relation to incident osteoporotic fracture after AI therapy during a median 6.1 years of follow up after AI initiation. In multivariable models, patients with the second and third highest tertiles (T) versus the lowest tertile of gSOS had significantly lower risk of fracture (T2: adjusted HR = 0.61, 95% CI: 0.46-0.80; T3: adjusted HR = 0.53, 95% CI: 0.40-0.70). The lower risk of fracture in patients with the highest tertile of gSOS remained significant after further adjustment for BMD at the hip (T3: adjusted HR = 0.62, 95% CI: 0.42-0.91). In conclusion, our analysis showed gSOS as a novel genetic predictor for fracture risk independent of BMD among breast cancer patients treated with AIs. Future studies are warranted to evaluate the performance of incorporating gSOS in prediction models for the risk of AI-related fracture in breast cancer patients.

Sample Collection Methods in Upper Gastrointestinal Research

In recent years, significant translational research advances have been made in the upper gastrointestinal (GI) research field. Endoscopic evaluation is a reasonable option for acquiring upper GI tissue for research purposes because it has minimal risk and can be applied to unresectable gastric cancer. The optimal number of biopsy samples and sample storage is crucial and might influence results. Furthermore, the methods for sample acquisition can be applied differently according to the research purpose; however, there have been few reports on methods for sample collection from endoscopic biopsies. In this review, we suggested a protocol for collecting study samples for upper GI research, including microbiome, DNA, RNA, protein, single-cell RNA sequencing, and organoid culture, through a comprehensive literature review. For microbiome analysis, one or two pieces of biopsied material obtained using standard endoscopic forceps may be sufficient. Additionally, 5 mL of gastric fluid and 3-4 mL of saliva is recommended for microbiome analyses. At least one gastric biopsy tissue is necessary for most DNA or RNA analyses, while proteomics analysis may require at least 2-3 biopsy tissues. Single cell-RNA sequencing requires at least 3-5 tissues and additional 1-2 tissues, if possible. For successful organoid culture, multiple sampling is necessary to improve the quality of specimens.

Performance characterization of PCR-free whole genome sequencing for clinical diagnosis

To evaluate the performance of polymerase chain reaction (PCR)-free whole genome sequencing (WGS) for clinical diagnosis, and thereby revealing how experimental parameters affect variant detection.Five NA12878 samples were sequenced using MGISEQ-2000. NA12878 samples underwent WGS with differing deoxyribonucleic acid (DNA) input and library preparation protocol (PCR-based vs PCR-free protocols for library preparation). The depth of coverage and genotype quality of each sample were compared. The performance of each sample was measured for sensitivity, coverage of depth and breadth of coverage of disease-related genes, and copy number variants. We also developed a systematic WGS pipeline (PCR-free) for the analysis of 11 clinical cases.In general, NA12878-2 (PCR-free WGS) showed better depth of coverage and genotype quality distribution than NA12878-1 (PCR-based WGS). With a mean depth of ∼40×, the sensitivity of homozygous and heterozygous single nucleotide polymorphisms (SNPs) of NA12878-2 showed higher sensitivity (>99.77% and >99.82%) than NA12878-1, and positive predictive value exceeded 99.98% and 99.07%. The sensitivity and positive predictive value of homozygous and heterozygous indels for NA12878-2 (PCR-free WGS) showed great improvement than NA128878-1. The breadths of coverage for disease-related genes and copy number variants are slightly better for samples with PCR-free library preparation protocol than the sample with PCR-based library preparation protocol. DNA input also influences the performance of variant detection in samples with PCR-free WGS. All the 19 previously confirmed variants in 11 clinical cases were successfully detected by our WGS pipeline (PCR free).Different experimental parameters may affect variant detection for clinical WGS. Clinical scientists should know the range of sensitivity of variants for different methods of WGS, which would be useful when interpreting and delivering clinical reports.

Whole exome sequencing studies in epilepsy: A deep analysis of the published literature

To evaluate the quality of whole-exome sequencing (WES) reporting in the epilepsy literature. We aimed to assess the quality of reporting of WES in epilepsy. We compared studies based on journal type and if outcome reporting biases exist. We used a self-constructed benchmark to quantitatively analyze studies. We included 451 publications. Reporting was heterogeneous with poor reporting of (1) ACMG guideline application 13% and (2) Human Phenotype Ontology (HPO) numbers in 3% of studies, 3) VUS in 19%. Predictors of reporting included journal type and journal impact factor. Date of publication and publication type were not predictors of poor reporting. Pairwise comparisons of genetics versus neurology journals using relative risks yielded significant differences in reporting of ACMG guideline application (RR 1.88, 95% CI 1.04-3.38); HPO numbers (RR 8.62, 95% CI 1.08-63.37) and deposition of findings to ClinVar (RR 2.50, 95% CI 1.03-6.1). Reporting of WES literature is heterogeneous in quality, and poor reporting hinders collaboration and accession of data into large databases like OMIM and OrphaNet. This study highlights reporting bias in this area and, formal structural guidelines like the CONSORT guidelines used in the reporting of clinical trials are needed to address the issue.

Direct RNA Nanopore Sequencing of SARS-CoV-2 Extracted from Critical Material from Swabs

In consideration of the increasing prevalence of COVID-19 cases in several countries and the resulting demand for unbiased sequencing approaches, we performed a direct RNA sequencing (direct RNA seq.) experiment using critical oropharyngeal swab samples collected from Italian patients infected with SARS-CoV-2 from the Palermo region in Sicily. Here, we identified the sequences SARS-CoV-2 directly in RNA extracted from critical samples using the Oxford Nanopore MinION technology without prior cDNA retrotranscription. Using an appropriate bioinformatics pipeline, we could identify mutations in the nucleocapsid (N) gene, which have been reported previously in studies conducted in other countries. In conclusion, to the best of our knowledge, the technique used in this study has not been used for SARS-CoV-2 detection previously owing to the difficulties in the extraction of RNA of sufficient quantity and quality from routine oropharyngeal swabs. Despite these limitations, this approach provides the advantages of true native RNA sequencing and does not include amplification steps that could introduce systematic errors. This study can provide novel information relevant to the current strategies adopted in SARS-CoV-2 next-generation sequencing.

A benchmarking study of SARS-CoV-2 whole-genome sequencing protocols using COVID-19 patient samples

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an emerging new type of coronavirus that is responsible for the COVID-19 pandemic and the unprecedented global health emergency. Whole-genome sequencing (WGS) of SARS-CoV-2 plays a critical role in understanding the disease. Performance variation exists across SARS-CoV-2 viral WGS technologies, but there is currently no benchmarking study comparing different WGS sequencing protocols. We compared seven different SARS-CoV-2 WGS library protocols using RNA from patient nasopharyngeal swab samples under two storage conditions with low and high viral inputs. We found large differences in mappability and genome coverage, and variations in sensitivity, reproducibility, and precision of single-nucleotide variant calling across different protocols. For certain amplicon-based protocols, an appropriate primer trimming step is critical for accurate single-nucleotide variant calling. We ranked the performance of protocols based on six different metrics. Our findings offer guidance in choosing appropriate WGS protocols to characterize SARS-CoV-2 and its evolution.

Evaluation of saliva as a source of accurate whole-genome and microbiome sequencing data

This study sets out to establish the suitability of saliva-based whole-genome sequencing (WGS) through a comparison against blood-based WGS. To fully appraise the observed differences, we developed a novel technique of pseudo-replication. We also investigated the potential of characterizing individual salivary microbiomes from non-human DNA fragments found in saliva. We observed that the majority of discordant genotype calls between blood and saliva fell into known regions of the human genome that are typically sequenced with low confidence; and could be identified by quality control measures. Pseudo-replication demonstrated that the levels of discordance between blood- and saliva-derived WGS data were entirely similar to what one would expect between technical replicates if an individual's blood or saliva had been sequenced twice. Finally, we successfully sequenced salivary microbiomes in parallel to human genomes as demonstrated by a comparison against the Human Microbiome Project.

Characterizing sensitivity and coverage of clinical WGS as a diagnostic test for genetic disorders

BACKGROUND

Due to its reduced cost and incomparable advantages, WGS is likely to lead to changes in clinical diagnosis of rare and undiagnosed diseases. However, the sensitivity and breadth of coverage of clinical WGS as a diagnostic test for genetic disorders has not been fully evaluated.

METHODS

Here, the performance of WGS in NA12878, the YH cell line, and the Chinese trios were measured by assessing their sensitivity, PPV, depth and breadth of coverage using MGISEQ-2000. We also compared the performance of WES and WGS using NA12878. The sensitivity and PPV were tested using the family-based trio design for the Chinese trios. We further developed a systematic WGS pipeline for the analysis of 8 clinical cases.

RESULTS

In general, the sensitivity and PPV for SNV/indel detection increased with mean depth and reached a plateau at an ~ 40X mean depth using down-sampling samples of NA12878. With a mean depth of 40X, the sensitivity of homozygous and heterozygous SNPs of NA12878 was > 99.25% and > 99.50%, respectively, and the PPV was 99.97% and 98.96%. Homozygous and heterozygous indels showed lower sensitivity and PPV. The sensitivity and PPV were still not 100% even with a mean depth of ~ 150X. We also observed a substantial variation in the sensitivity of CNV detection across different tools, especially in CNVs with a size less than 1 kb. In general, the breadth of coverage for disease-associated genes and CNVs increased with mean depth. The sensitivity and coverage of WGS (~ 40X) was better than WES (~ 120X). Among the Chinese trios with an ~ 40X mean depth, the sensitivity among offspring was > 99.48% and > 96.36% for SNP and indel detection, and the PPVs were 99.86% and 97.93%. All 12 previously validated variants in the 8 clinical cases were successfully detected using our WGS pipeline.

CONCLUSIONS

The current standard of a mean depth of 40X may be sufficient for SNV/indel detection and identification of most CNVs. It would be advisable for clinical scientists to determine the range of sensitivity and PPV for different classes of variants for a particular WGS pipeline, which would be useful when interpreting and delivering clinical reports.

Whole-genome sequencing data of Kazakh individuals

Objectives

Kazakhstan is a Central Asian crossroad of European and Asian populations situated along the way of the Great Silk Way. The territory of Kazakhstan has historically been inhabited by nomadic tribes and today is the multi-ethnic country with the dominant Kazakh ethnic group. We sequenced and analyzed the whole-genomes of five ethnic healthy Kazakh individuals with high coverage using next-generation sequencing platform. This whole-genome sequence data of healthy Kazakh individuals can be a valuable reference for biomedical studies investigating disease associations and population-wide genomic studies of ethnically diverse Central Asian region.

Data description

Blood samples have been collected from five ethnic healthy Kazakh individuals living in Kazakhstan. The genomic DNA was extracted from blood and sequenced. Sequencing was performed on Illumina HiSeq2000 next-generation sequencing platform. We sequenced and analyzed the whole-genomes of ethnic Kazakh individuals with the coverage ranging from 26 to 32X. Ranging from 98.85 to 99.58% base pairs were totally mapped and aligned on the human reference genome GRCh37 hg19. Het/Hom and Ts/Tv ratios for each whole genome ranged from 1.35 to 1.49 and from 2.07 to 2.08, respectively. Sequencing data are available in the National Center for Biotechnology Information SRA database under the accession number PRJNA374772.

Schizophrenic Psychosis Symptoms in a Background of Mild-To-Moderate Carnitine Palmitoyltransferase II Deficiency: A Case Report

Schizophrenia is a multifaceted mental illness characterized by cognitive and neurobehavioral abnormalities. Carnitine palmitoyltransferase II (CPT II) deficiency is a metabolic disorder resulting in impaired transport of long-chain fatty acids from the cytosol to the mitochondrial inner membrane, where fatty acid β-oxidation takes place. Here, we present an interesting clinical case of an adolescent male that presented with psychosis and a history of mild-to-moderate CPT II deficiency. To identify germline genetic variation that may contribute to the phenotypes observed, we performed whole-exome sequencing on DNA from the proband, unaffected fraternal twin, and biological parents. The proband was identified to be homozygous for the p.Val368Ile and heterozygous for the p.Met647Val variant in CPT2. Each of these variants are benign on their own; however, their combined effect is unclear. Further, variation was identified in the dopamine β-hydroxylase (DBH) gene (c.339+2T>C), which may contribute to decreased activity of DBH; however, based on the patient’s presentation, severe DBH deficiency is unlikely. In conclusion, the variants identified in this study do not clearly explain the observed patient phenotypes, indicating that the complex phenotypes are likely caused by an interplay of genetic and environmental factors that warrant further investigation.

Quality of whole genome sequencing from blood versus saliva derived DNA in cardiac patients

Background

Whole-genome sequencing (WGS) is becoming an increasingly important tool for detecting genomic variation. Blood derived DNA is the current standard for WGS for research or clinical purposes but may not always be feasible to acquire. The usability of DNA from saliva for WGS is not known. We compared the quality of WGS between blood versus saliva derived DNA.

Methods

WGS was performed in DNA from 531 blood and 502 saliva samples (including 5 paired samples) from participants enrolled in a heart disease biorepository. We compared the proportion of sequencing reads that mapped to non-human sources (microbiome), the sequencing coverage, and the yield and concordance of single nucleotide variant (SNV) and copy number variant (CNV) calls between blood and saliva genomes.

Results

Of 531 blood and 502 saliva samples, 46% saliva DNA failed quality control (QC) requirements for WGS compared to 6% QC failure for blood DNA. An average of 10.7% WGS reads in the saliva samples mapped to the human oral microbiome compared to 0.09% WGS reads in blood samples. However, these reads were readily excluded by excluding reads that did not map to the human reference genome. Sequencing coverage met or exceeded the target sequencing depth of 30x in all the blood samples and 4 of the 5 saliva samples; the fifth saliva sample had an average sequencing depth of 22.6x. Over 95% of SNVs identified in saliva were concordant with those identified in blood across the genome, within all gene coding regions, and within cardiovascular disease-related gene coding regions. Rare SNVs, defined as those with a minor allele frequency of less than 1% in the Genome Aggregation Database, had a lower concordance of 90% between blood and saliva genomes. CNVs had only 76% concordance between blood and saliva samples.

Conclusions

High quality saliva samples that meet stringent QC criteria can be used for WGS when blood-derived DNA is not available or is not suitable. Saliva DNA provides an acceptable yield of SNV calls but has a lower yield for CNV calls compared to blood DNA.

Recommendations on Collecting and Storing Samples for Genetic Studies in Hearing and Tinnitus Research

Objectives:

Research on the genetic basis of tinnitus is still in its first steps. A group of scientists dedicated to tinnitus genetics within European Tinnitus Network (TINNET) network recognize that further progress requires multicenter collaborative efforts for defining contributing genes. The purpose of the present work is to provide instructions regarding collection, processing, storage, and shipment of samples intended for genetic studies in auditory research.

Design:

One part of the recommendations has a general character; another part is of particular importance for auditory healthcare practitioners such as otolaryngology physicians, audiologists, and general practitioners.

Results:

We provide a set of instructions and various options for obtaining samples. We give advice regarding sample processing, storage, and shipment and define the minimal and essential clinical information that should accompany the samples collected for genetic processing.

Conclusions:

These recommendations offer a basis to standardize and optimize collaborations between geneticists and healthcare practitioners specialized in tinnitus and hearing disorders.

Comparison of SureSelect and Nextera Exome Capture Performance in Single-Cell Sequencing

BACKGROUND

Advances in single-cell sequencing provide unprecedented opportunities for clinical examination of circulating tumor cells, cancer stem cells, and other rare cells responsible for disease progression and drug resistance. On the genomic level, single-cell whole exome sequencing (scWES) started to gain popularity with its unique potentials in characterizing mutational landscapes at a single-cell level. Currently, there is little known about the performance of different exome capture kits in scWES. Nextera rapid capture (NXT; Illumina, Inc.) has been the only exome capture kit recommended for scWES by Fluidigm C1, a widely accessed system in single-cell preparation.

RESULTS

In this study, we compared the performance of NXT following Fluidigm's protocol with Agilent SureSelectXT Target Enrichment System (AGL), another exome capture kit widely used for bulk sequencing. We created DNA libraries of 192 single cells isolated from spheres grown from a melanoma specimen using Fluidigm C1. Twelve high-yield cells were selected to perform dual-exome capture and sequencing using AGL and NXT in parallel. After mapping and coverage analysis, AGL outperformed NXT in coverage uniformity, mapping rates of reads, exome capture rates, and low PCR duplicate rates. For germline variant calling, AGL achieved better performance in overlap with known variants in dbSNP and transition-transversion ratios. Using calls from high coverage bulk sequencing from blood DNA as the golden standard, AGL-based scWES demonstrated high positive predictive values, and medium to high sensitivity. Lastly, we evaluated somatic mutation calling by comparing single-cell data with the matched blood sequence as control. On average, 300 mutations were identified in each cell. In 10 of 12 cells, higher numbers of mutations were identified using AGL than NXT, probably caused by coverage depth. When mutations are adequately covered in both AGL and NXT data, the two methods showed very high concordance (93-100% per cell).

CONCLUSIONS

Our results suggest that AGL can also be used for scWES when there is sufficient DNA, and it yields better data quality than the current Fluidigm's protocol using NXT.

Measuring coverage and accuracy of whole-exome sequencing in clinical context

PURPOSE

To evaluate the coverage and accuracy of whole-exome sequencing (WES) across vendors.

METHODS

Blood samples from three trios underwent WES at three vendors. Relative performance of the three WES services was measured for breadth and depth of coverage. The false-negative rates (FNRs) were estimated using the segregation pattern within each trio.

RESULTS

Mean depth of coverage for all genes was 189.0, 124.9, and 38.3 for the three vendor services. Fifty-five of the American College of Medical Genetics and Genomics 56 genes, but only 56 of 63 pharmacogenes, were 100% covered at 10 × in at least one of the nine individuals for all vendors; however, there was substantial interindividual variability. For the two vendors with mean depth of coverage >120 ×, analytic positive predictive values (aPPVs) exceeded 99.1% for single-nucleotide variants and homozygous indels, and sensitivities were 98.9-99.9%; however, heterozygous indels showed lower accuracy and sensitivity. Among the trios, FNRs in the offspring were 0.07-0.62% at well-covered variants concordantly called in both parents.

CONCLUSION

The current standard of 120 × coverage for clinical WES may be insufficient for consistent breadth of coverage across the exome. Ordering clinicians and researchers would benefit from vendors' reports that estimate sensitivity and aPPV, including depth of coverage across the exome.

A framework for the estimation of the proportion of true discoveries in single nucleotide variant detection studies for human data

Any single nucleotide variant detection study could benefit from a fast and cheap method of measuring the quality of variant call list. It is advantageous to be able to see how the call list quality is affected by different variant filtering thresholds and other adjustments to the study parameters. Here we look into a possibility of estimating the proportion of true positives in a single nucleotide variant call list for human data. Using whole-exome and whole-genome gold standard data sets for training, we focus on building a generic model that only relies on information available from any variant caller. We assess and compare the performance of different candidate models based on their practical accuracy. We find that the generic model delivers decent accuracy most of the time. Further, we conclude that its performance could be improved substantially by leveraging the variant quality metrics that are specific to each variant calling tool.

Pitfalls of improperly procured adjacent non-neoplastic tissue for somatic mutation analysis using next-generation sequencing

Background

The rapid adoption of next-generation sequencing provides an efficient system for detecting somatic alterations in neoplasms. The detection of such alterations requires a matched non-neoplastic sample for adequate filtering of non-somatic events such as germline polymorphisms. Non-neoplastic tissue adjacent to the excised neoplasm is often used for this purpose as it is simultaneously collected and generally contains the same tissue type as the neoplasm. Following NGS analysis, we and others have frequently observed low-level somatic mutations in these non-neoplastic tissues, which may impose additional challenges to somatic mutation detection as it complicates germline variant filtering.

Methods

We hypothesized that the low-level somatic mutation observed in non-neoplastic tissues may be entirely or partially caused by inadvertent contamination by neoplastic cells during the surgical pathology gross assessment or tissue procurement process. To test this hypothesis, we applied a systematic protocol designed to collect multiple grossly non-neoplastic tissues using different methods surrounding each single neoplasm. The procedure was applied in two breast cancer lumpectomy specimens. In each case, all samples were first sequenced by whole-exome sequencing to identify somatic mutations in the neoplasm and determine their presence in the adjacent non-neoplastic tissues. We then generated ultra-deep coverage using targeted sequencing to assess the levels of contamination in non-neoplastic tissue samples collected under different conditions.

Results

Contamination levels in non-neoplastic tissues ranged up to 3.5 and 20.9 % respectively in the two cases tested, with consistent pattern correlated with the manner of grossing and procurement. By carefully controlling the conditions of various steps during this process, we were able to eliminate any detectable contamination in both patients.

Conclusion

The results demonstrated that the process of tissue procurement contributes to the level of contamination in non-neoplastic tissue, and contamination can be reduced to below detectable levels by using a carefully designed collection method. A standard protocol dedicated for acquiring adjacent non-neoplastic tissue that minimizes neoplasm contamination should be implemented for all somatic mutation detection studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-016-0226-1) contains supplementary material, which is available to authorized users.

The minimal amount of starting DNA for Agilent's hybrid capture-based targeted massively parallel sequencing

Targeted capture massively parallel sequencing is increasingly being used in clinical settings, and as costs continue to decline, use of this technology may become routine in health care. However, a limited amount of tissue has often been a challenge in meeting quality requirements. To offer a practical guideline for the minimum amount of input DNA for targeted sequencing, we optimized and evaluated the performance of targeted sequencing depending on the input DNA amount. First, using various amounts of input DNA, we compared commercially available library construction kits and selected Agilent’s SureSelect-XT and KAPA Biosystems’ Hyper Prep kits as the kits most compatible with targeted deep sequencing using Agilent’s SureSelect custom capture. Then, we optimized the adapter ligation conditions of the Hyper Prep kit to improve library construction efficiency and adapted multiplexed hybrid selection to reduce the cost of sequencing. In this study, we systematically evaluated the performance of the optimized protocol depending on the amount of input DNA, ranging from 6.25 to 200 ng, suggesting the minimal input DNA amounts based on coverage depths required for specific applications.