Variant Calling Pipeline Research Articles

Case report: A case study of variant calling pipeline selection effect on the molecular diagnostics outcome

Next-generation sequencing technologies have not only defined a breakthrough in medical genetics, but also been able to enter routine clinical practice to determine individual genetic susceptibilities. Modern technological developments are routinely introduced to genetic analysis overtaking the established approaches, potentially raising a number of challenges. To what extent is the advantage of new methodologies in synthetic metrics, such as precision and recall, more important than stability and reproducibility? Could differences in the technical protocol for calling variants be crucial to the diagnosis and, by extension, the patient’s treatment strategy? A regulatory review process may delay the incorporation of potentially beneficial technologies, resulting in missed opportunities to make the right medical decisions. On the other hand, a blind adoption of new technologies based solely on synthetic metrics of precision and recall can lead to incorrect conclusions and adverse outcomes for the specific patient. Here, we use the example of a patient with a WHO-diagnosed desmoplastic/nodular SHH-medulloblastoma to explore how the choice of DNA variant search protocol affects the genetic diagnostics outcome.

A phylogenetics and variant calling pipeline to support SARS-CoV-2 genomic epidemiology in the UK.

In response to the escalating SARS-CoV-2 pandemic, in March 2020 the COVID-19 Genomics UK (COG-UK) consortium was established to enable national-scale genomic surveillance in the UK. By the end of 2020, 49% of all SARS-CoV-2 genome sequences globally had been generated as part of the COG-UK programme, and to date, this system has generated >3 million SARS-CoV-2 genomes. Rapidly and reliably analysing this unprecedented number of genomes was an enormous challenge. To fulfil this need and to inform public health decision-making, we developed a centralized pipeline that performs quality control, alignment, and variant calling and provides the global phylogenetic context of sequences. We present this pipeline and describe how we tailored it as the pandemic progressed to scale with the increasing amounts of data and to provide the most relevant analyses on a daily basis.

StratoMod: predicting sequencing and variant calling errors with interpretable machine learning

Despite the variety in sequencing platforms, mappers, and variant callers, no single pipeline is optimal across the entire human genome. Therefore, developers, clinicians, and researchers need to make tradeoffs when designing pipelines for their application. Currently, assessing such tradeoffs relies on intuition about how a certain pipeline will perform in a given genomic context. We present StratoMod, which addresses this problem using an interpretable machine-learning classifier to predict germline variant calling errors in a data-driven manner. We show StratoMod can precisely predict recall using Hifi or Illumina and leverage StratoMod’s interpretability to measure contributions from difficult-to-map and homopolymer regions for each respective outcome. Furthermore, we use Statomod to assess the effect of mismapping on predicted recall using linear vs. graph-based references, and identify the hard-to-map regions where graph-based methods excelled and by how much. For these we utilize our draft benchmark based on the Q100 HG002 assembly, which contains previously-inaccessible difficult regions. Furthermore, StratoMod presents a new method of predicting clinically relevant variants likely to be missed, which is an improvement over current pipelines which only filter variants likely to be false. We anticipate this being useful for performing precise risk-reward analyses when designing variant calling pipelines.

Fast and accurate DNASeq variant calling workflow composed of LUSH toolkit

BackgroundWhole genome sequencing (WGS) is becoming increasingly prevalent for molecular diagnosis, staging and prognosis because of its declining costs and the ability to detect nearly all genes associated with a patient's disease. The currently widely accepted variant calling pipeline, GATK, is limited in terms of its computational speed and efficiency, which cannot meet the growing analysis needs.ResultsHere, we propose a fast and accurate DNASeq variant calling workflow that is purely composed of tools from LUSH toolkit. The precision and recall measurements indicate that both the LUSH and GATK pipelines exhibit high levels of consistency, with precision and recall rates exceeding 99% on the 30x NA12878 dataset. In terms of processing speed, the LUSH pipeline outperforms the GATK pipeline, completing 30x WGS data analysis in just 1.6 h, which is approximately 17 times faster than GATK. Notably, the LUSH_HC tool completes the processing from BAM to VCF in just 12 min, which is around 76 times faster than GATK.ConclusionThese findings suggest that the LUSH pipeline is a highly promising alternative to the GATK pipeline for WGS data analysis, with the potential to significantly improve bedside analysis of acutely ill patients, large-scale cohort data analysis, and high-throughput variant calling in crop breeding programs. Furthermore, the LUSH pipeline is highly scalable and easily deployable, allowing it to be readily applied to various scenarios such as clinical diagnosis and genomic research.

Abstract C102: Polyethnic-1000: Advancing cancer genomics by studying ethnically diverse, underserved patient populations in New York

Abstract Prompted by the considerable biases toward European ancestry cases in large cancer genomics repositories, the New York Genome Center initiated an effort to diversify the public datasets that are forming the foundations of current cancer research, in order to optimize the unbiased applications of modern genomically-driven discoveries to a larger population. Our effort was focused on 8 cancer types (breast, prostate, bladder, lung, pancreatic, colon and endometrial cancer, as well as multiple myeloma) and generated Illumina tumor-normal whole-genome sequencing and tumor transcriptome sequencing. The majority of the samples were retrospectively obtained in clinics and hospitals of the New York region and the analyses are conducted with collaborators in these institutions. Our dataset is currently composed of 851 tumor-normal WGS and 471 RNASeq, and is progressively being distributed to the scientific community in partnership with ICGC-ARGO. Overall, 75% of the cohort is predominantly of non-European ancestry, as estimated by genetic ancestry analysis using the 1000Genomes continental references.We analyzed the samples with NYGC’s somatic variant calling pipeline and will present results from the 8 cancer types related to germline and somatic variants (single-nucleotide variants, indels, copy-number alterations, complex structural variants) as well as mutational signatures, neoantigen prediction, fusion genes and expression profiles, in relation to genetic ancestry estimation and clinical features of the patients. We anticipate that the dataset will significantly diversify the genetic ancestry of tumor profiles in public databases, and serve to form hypotheses regarding the contribution of ancestry on cancer disparities. The experience of the Polyethnic-1000 project informs our contribution to the Cancer Grand Challenge SAMBAI, which focuses on breast, prostate and pancreatic cancers in the African Diaspora. SAMBAI will further investigate the causes of cancer disparity by combining genomic data with social determinants of health and exposomics information, deeper immune profiling and new genomics technologies such as long reads sequencing and the use of a pangenome reference. Citation Format: Nicolas Robine, Lara Winterkorn, Tim Chu, Zoe Goldstein, Will Hooper, Ali Oku, Heather M. Geiger, Melissa B. Davis, Onyinye Balogun. Polyethnic-1000: Advancing cancer genomics by studying ethnically diverse, underserved patient populations in New York [abstract]. In: Proceedings of the 17th AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2024 Sep 21-24; Los Angeles, CA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2024;33(9 Suppl):Abstract nr C102.

198 Discovery of key functional SNP markers associated with feed efficiency in beef cattle

Abstract Beef production contributes to approximately 2.4% of the total Greenhouse Gas emissions in Canada. Improving beef cattle feed efficiency (FE) may lead to improved energy partitioning and resource use, thereby improving the sustainability of the beef industry. Our objective was to identify differentially expressed key regulatory genes and associated novel functional SNP markers linked to FE in Canadian beef cattle using RNA-Sequencing (RNA-Seq) of rumen tissue collected from 48 beef cattle [n = 16 Angus, n = 16 Charolais, n = 16 Kinsella (Composite Hybrid including Angus, Charolais, Galloway, Hereford, Holstein, Brown Swiss, and Simmental)] selected for extreme FE phenotypes. In total, 11 key regulatory genes (MYH1, MYL2, MYLPF, TNNC2, EIF4B, RHOD, TCEANC, CKM, ENSBTAG00000040518, SERPINB2, and USP43) were significantly differentially expressed (DE) between extreme Residual Feed Intake (RFI) groups (low-RFI n = 8, high-RFI n = 8 per breed) using CLC Genomics Workbench (FDR < 0.05; |FC| >2). Using an optimized RNA-Seq variant calling pipeline using STAR and BCFtools, a total of [total (unique to low-RFI, unique to high-RFI)] 75 (36, 39), 78 (42, 36), and 53 (35, 18) uniquely fixed functional SNPs were located within coding regions of these 11 functional candidate genes, in low- and high-RFI animals in the Angus, Charolais, and Kinsella breeds, respectively. Considering all functional SNPs uniquely identified in low- or high-RFI groups for all breed comparisons, the majority of SNPs were identified in MYH1, EIF4B, and SERPINB2 genes, which function together in metabolically demanding biological processes (P < 0.05) related to muscle contraction, muscle system, and muscle filament function processes. Additionally, EIF4B and SERPINB2 genes were found to have a role in signaling pathways that coordinate cell growth and immune function, respectively. Using RNA-Seq to identify key regulatory genes and associated functional SNPs linked to FE may uncover important genetic markers that influence the regulation of FE in beef cattle.

NCBench: providing an open, reproducible, transparent, adaptable, and continuous benchmark approach for DNA-sequencing-based variant calling

We present the results of the human genomic small variant calling benchmarking initiative of the German Research Foundation (DFG) funded Next Generation Sequencing Competence Network (NGS-CN) and the German Human Genome-Phenome Archive (GHGA). In this effort, we developed NCBench, a continuous benchmarking platform for the evaluation of small genomic variant callsets in terms of recall, precision, and false positive/negative error patterns. NCBench is implemented as a continuously re-evaluated open-source repository. We show that it is possible to entirely rely on public free infrastructure (Github, Github Actions, Zenodo) in combination with established open-source tools. NCBench is agnostic of the used dataset and can evaluate an arbitrary number of given callsets, while reporting the results in a visual and interactive way. We used NCBench to evaluate over 40 callsets generated by various variant calling pipelines available in the participating groups that were run on three exome datasets from different enrichment kits and at different coverages. While all pipelines achieve high overall quality, subtle systematic differences between callers and datasets exist and are made apparent by NCBench.These insights are useful to improve existing pipelines and develop new workflows. NCBench is meant to be open for the contribution of any given callset. Most importantly, for authors, it will enable the omission of repeated re-implementation of paper-specific variant calling benchmarks for the publication of new tools or pipelines, while readers will benefit from being able to (continuously) observe the performance of tools and pipelines at the time of reading instead of at the time of writing.

Variant calling in polyploids for population and quantitative genetics.

Advancements in genome assembly and sequencing technology have made whole genome sequence (WGS) data and reference genomes accessible to study polyploid species. Compared to popular reduced-representation sequencing approaches, the genome-wide coverage and greater marker density provided by WGS data can greatly improve our understanding of polyploid species and polyploid biology. However, biological features that make polyploid species interesting also pose challenges in read mapping, variant identification, and genotype estimation. Accounting for characteristics in variant calling like allelic dosage uncertainty, homology between subgenomes, and variance in chromosome inheritance mode can reduce errors. Here, I discuss the challenges of variant calling in polyploid WGS data and discuss where potential solutions can be integrated into a standard variant calling pipeline.

Emergence of the novel sixth Candida auris Clade VI in Bangladesh.

Candida auris, initially identified in 2009, has rapidly become a critical concern due to its antifungal resistance and significant mortality rates in healthcare-associated outbreaks. To date, whole-genome sequencing (WGS) has identified five unique clades of C. auris, with some strains displaying resistance to all primary antifungal drug classes. In this study, we presented the first WGS analysis of C. auris from Bangladesh, describing its origins, transmission dynamics, and antifungal susceptibility testing (AFST) profile. Ten C. auris isolates collected from hospital settings in Bangladesh were initially identified by CHROMagar Candida Plus, followed by VITEK2 system, and later sequenced using Illumina NextSeq 550 system. Reference-based phylogenetic analysis and variant calling pipelines were used to classify the isolates in different clades. All isolates aligned ~90% with the Clade I C. auris B11205 reference genome. Of the 10 isolates, 8 were clustered with Clade I isolates, highlighting a South Asian lineage prevalent in Bangladesh. Remarkably, the remaining two isolates formed a distinct cluster, exhibiting >42,447 single-nucleotide polymorphism differences compared to their closest Clade IV counterparts. This significant variation corroborates the emergence of a sixth clade (Clade VI) of C. auris in Bangladesh, with potential for international transmission. AFST results showed that 80% of the C. auris isolates were resistant to fluconazole and voriconazole, whereas Clade VI isolates were susceptible to azoles, echinocandins, and pyrimidine analogue. Genomic sequencing revealed ERG11_Y132F mutation conferring azole resistance while FCY1_S70R mutation found inconsequential in describing 5-flucytosine resistance. Our study underscores the pressing need for comprehensive genomic surveillance in Bangladesh to better understand the emergence, transmission dynamics, and resistance profiles of C. auris infections. Unveiling the discovery of a sixth clade (Clade VI) accentuates the indispensable role of advanced sequencing methodologies.IMPORTANCECandida auris is a nosocomial fungal pathogen that is commonly misidentified as other Candida species. Since its emergence in 2009, this multidrug-resistant fungus has become one of the five urgent antimicrobial threats by 2019. Whole-genome sequencing (WGS) has proven to be the most accurate identification technique of C. auris which also played a crucial role in the initial discovery of this pathogen. WGS analysis of C. auris has revealed five distinct clades where isolates of each clade differ among themselves based on pathogenicity, colonization, infection mechanism, as well as other phenotypic characteristics. In Bangladesh, C. auris was first reported in 2019 from clinical samples of a large hospital in Dhaka city. To understand the origin, transmission dynamics, and antifungal-resistance profile of C. auris isolates circulating in Bangladesh, we conducted a WGS-based surveillance study on two of the largest hospital settings in Dhaka, Bangladesh.

Negligible effects of read trimming on the accuracy of germline short variant calling in the human genome

Background Next generation sequencing (NGS) has become a standard tool in the molecular diagnostics of Mendelian disease, and the precision of such diagnostics is greatly affected by the accuracy of variant calling from sequencing data. Recently, we have comprehensively evaluated the performance of multiple variant calling pipelines. However, no systematic analysis of the effects of read trimming on variant discovery with modern variant calling software has yet been performed. Methods In this work, we systematically evaluated the effects of adapters on the performance of 8 variant calling and filtering methods using 14 standard reference Genome-in-a-Bottle (GIAB) samples. Variant calls were compared to the ground truth variant sets, and the effect of adapter trimming with different tools was assessed using major performance metrics (precision, recall, and F1 score). Results We show that adapter trimming has no effect on the accuracy of the best-performing variant callers (e.g., DeepVariant) on whole-genome sequencing (WGS) data. For whole-exome sequencing (WES) datasets subtle improvement of accuracy was observed in some of the samples. In high-coverage WES data (~200x mean coverage), adapter removal allowed for discovery of 2-4 additional true positive variants in only two out of seven datasets tested. Moreover, this effect was not dependent on the median insert size and proportion of adapter sequences in reads. Surprisingly, the effect of trimming on variant calling was reversed when moderate coverage (~80-100x) WES data was used. Finally, we show that some of the recently developed machine learning-based variant callers demonstrate greater dependence on the presence of adapters in reads. Conclusions Taken together, our results indicate that adapter removal is unnecessary when calling germline variants, but suggest that preprocessing methods should be carefully chosen when developing and using machine learning-based variant analysis methods.

Scalable and efficient DNA sequencing analysis on different compute infrastructures aiding variant discovery.

DNA variation analysis has become indispensable in many aspects of modern biomedicine, most prominently in the comparison of normal and tumor samples. Thousands of samples are collected in local sequencing efforts and public databases requiring highly scalable, portable, and automated workflows for streamlined processing. Here, we present nf-core/sarek 3, a well-established, comprehensive variant calling and annotation pipeline for germline and somatic samples. It is suitable for any genome with a known reference. We present a full rewrite of the original pipeline showing a significant reduction of storage requirements by using the CRAM format and runtime by increasing intra-sample parallelization. Both are leading to a 70% cost reduction in commercial clouds enabling users to do large-scale and cross-platform data analysis while keeping costs and CO2 emissions low. The code is available at https://nf-co.re/sarek.

Abstract 4926: Advancements in somatic variant calling from UG100 whole genome and whole exome sequencing data

Abstract Somatic variant calling involves the identification of genomic alterations that occur in somatic cells, requiring deep coverage to enable high sensitivity for low-frequency variants. Characterizing somatic variants across the entire genome therefore benefits from novel cost-efficient sequencing platforms, such as UG100. Here, we present optimization of variant calling tools for short and structural variants on WGS and WES data from UG100. For calling short variants, we optimized DeepVariant (DV) for somatic calling using data from matched tumor-normal sample pairs, improving both variant calling accuracy and pipeline running time (up to 10-fold). We defined the task of somatic variant calling as deciding if the pileup image containing reads from the tumor and normal samples represents a true somatic variant (vs a germline variant or artifact). The challenge of robust variant calling using deep learning models is exacerbated in somatic calling, where sequencing depth and coverage variability are typically high. Our optimized DV overcomes these challenges by several data sampling strategies. First, allele-frequency preserving down-sampling reduces randomness of read sub-sampling in high coverage regions. Second, alternative allele prioritization samples alt-allele supporting reads first allowing to call variants at very high coverage loci without sacrificing sensitivity and computational efficiency. Finally, a Panel-of-Normals based on targeted WES data provides an additional improvement of precision for this assay type. We used these strategies to train two models, one for tumor characterization using WGS (T/N coverage: 40x-150x/40x-100x), and one for deep WES (T/N coverage: >500x/>120x). We called variants on simulated tumors using the WGS model. For VAF>10% the model showed SNV recall >98% and indel recall >95% with false-positive rate of 0.2/Mb. For VAF range of 5-10%, indel recall was 67% and SNV recall was 86%. To demonstrate the utility of our somatic variant calling, we applied the models to call somatic variants from well characterized cancer cell lines: COLO829, HCC1395 and HCC1143. Results showed F1>90% for variants with VAF>10%. The WES model was used to reliably call variants at VAF>5% on simulated tumors with average SNV recall of 99% with precision >99% and indel recall >86% with precision >94%. To analyze structural and copy-number variations, we optimized the assembly engine of GRIDSS to enable fast calling of structural variations and demonstrate that Control-FREEC can be used to call copy number variants. SV calling on COLO829/COLO829BL achieved sensitivity >95%. In conclusion, our research highlights the utility of UG100 within the field of oncology, demonstrating its capacity for comprehensive and precise somatic variant detection, both on WGS and WES data. Citation Format: Doron Shem-Tov, Maya Levy, Gil Hornung, Ilya Soifer, Hila Benjamin, Ariel Jaimovich, Adam Blattler, William Brandler, Robert Sugar, Isaac Kinde, Omer Barad, Doron Lipson. Advancements in somatic variant calling from UG100 whole genome and whole exome sequencing data [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 4926.

ArCH: improving the performance of clonal hematopoiesis variant calling and interpretation.

The acquisition of somatic mutations in hematopoietic stem and progenitor stem cells with resultant clonal expansion, termed clonal hematopoiesis (CH), is associated with increased risk of hematologic malignancies and other adverse outcomes. CH is generally present at low allelic fractions, but clonal expansion and acquisition of additional mutations leads to hematologic cancers in a small proportion of individuals. With high depth and high sensitivity sequencing, CH can be detected in most adults and its clonal trajectory mapped over time. However, accurate CH variant calling is challenging due to the difficulty in distinguishing low frequency CH mutations from sequencing artifacts. The lack of well-validated bioinformatic pipelines for CH calling may contribute to lack of reproducibility in studies of CH. Here, we developed ArCH, an Artifact filtering Clonal Hematopoiesis variant calling pipeline for detecting single nucleotide variants and short insertions/deletions by combining the output of four variant calling tools and filtering based on variant characteristics and sequencing error rate estimation. ArCH is an end-to-end cloud-based pipeline optimized to accept a variety of inputs with customizable parameters adaptable to multiple sequencing technologies, research questions, and datasets. Using deep targeted sequencing data generated from six acute myeloid leukemia patient tumor: normal dilutions, 31 blood samples with orthogonal validation, and 26 blood samples with technical replicates, we show that ArCH improves the sensitivity and positive predictive value of CH variant detection at low allele frequencies compared to standard application of commonly used variant calling approaches. The code for this workflow is available at: https://github.com/kbolton-lab/ArCH.

Clinical application of tumour-in-normal contamination assessment from whole genome sequencing

The unexpected contamination of normal samples with tumour cells reduces variant detection sensitivity, compromising downstream analyses in canonical tumour-normal analyses. Leveraging whole-genome sequencing data available at Genomics England, we develop a tool for normal sample contamination assessment, which we validate in silico and against minimal residual disease testing. From a systematic review of 771\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$771$$\\end{document} patients with haematological malignancies and sarcomas, we find contamination across a range of cancer clinical indications and DNA sources, with highest prevalence in saliva samples from acute myeloid leukaemia patients, and sorted CD3+ T-cells from myeloproliferative neoplasms. Further exploration reveals 108 hotspot mutations in genes associated with haematological cancers at risk of being subtracted by standard variant calling pipelines. Our work highlights the importance of contamination assessment for accurate somatic variants detection in research and clinical settings, especially with large-scale sequencing projects being utilised to deliver accurate data from which to make clinical decisions for patient care.

Systematic comparison of variant calling pipelines of target genome sequencing cross multiple next-generation sequencers.

Targeted genomic sequencing (TS) greatly benefits precision oncology by rapidly detecting genetic variations with better accuracy and sensitivity owing to its high sequencing depth. Multiple sequencing platforms and variant calling tools are available for TS, making it excruciating for researchers to choose. Therefore, benchmarking study across different platforms and pipelines available for TS is imperative. In this study, we performed a TS of Reference OncoSpan FFPE (HD832) sample enriched by TSO500 panel using four commercially available sequencers, and analyzed the output 50 datasets using five commonly-used bioinformatics pipelines. We systematically investigated the sequencing quality and variant detection sensitivity, expecting to provide optimal recommendations for future research. Four sequencing platforms returned highly concordant results in terms of base quality (Q20 > 94%), sequencing coverage (>97%) and depth (>2000×). Benchmarking revealed good concordance of variant calling across different platforms and pipelines, among which, FASTASeq 300 platform showed the highest sensitivity (100%) and precision (100%) in high-confidence variants calling when analyzed by SNVer and VarScan 2 algorithms. Furthermore, this sequencer demonstrated the shortest sequencing time (∼21h) at the sequencing mode PE150. Through the intersection of 50 datasets generated in this study, we recommended a novel set of variant genes outside the truth set published by HD832, expecting to replenish HD832 for future research on tumor variant diagnosis. Besides, we applied these five tools to another panel (TargetSeq One) for Twist cfDNA Pan-cancer Reference Standard, comprehensive consideration of SNP and InDel sensitivity, SNVer and VarScan 2 performed best among them. Furthermore, SNVer and VarScan 2 also performed best for six cancer cell lines samples regarding SNP and InDel sensitivity. Considering the dissimilarity of variant calls across different pipelines for datasets from the same platform, we recommended an integration of multiple tools to improve variant calling sensitivity and accuracy for the cancer genome. Illumina and GeneMind technologies can be used independently or together by public health laboratories performing tumor TS. SNVer and VarScan 2 perform better regarding variant detection sensitivity for three typical tumor samples. Our study provides a standardized target sequencing resource to benchmark new bioinformatics protocols and sequencing platforms.

Incorporating CNV analysis improves the yield of exome sequencing for rare monogenic disorders-an important consideration for resource-constrained settings.

Exome sequencing (ES) is a recommended first-tier diagnostic test for many rare monogenic diseases. It allows for the detection of both single-nucleotide variants (SNVs) and copy number variants (CNVs) in coding exonic regions of the genome in a single test, and this dual analysis is a valuable approach, especially in limited resource settings. Single-nucleotide variants are well studied; however, the incorporation of copy number variant analysis tools into variant calling pipelines has not been implemented yet as a routine diagnostic test, and chromosomal microarray is still more widely used to detect copy number variants. Research shows that combined single and copy number variant analysis can lead to a diagnostic yield of up to 58%, increasing the yield with as much as 18% from the single-nucleotide variant only pipeline. Importantly, this is achieved with the consideration of computational costs only, without incurring any additional sequencing costs. This mini review provides an overview of copy number variant analysis from exome data and what the current recommendations are for this type of analysis. We also present an overview on rare monogenic disease research standard practices in resource-limited settings. We present evidence that integrating copy number variant detection tools into a standard exome sequencing analysis pipeline improves diagnostic yield and should be considered a significantly beneficial addition, with relatively low-cost implications. Routine implementation in underrepresented populations and limited resource settings will promote generation and sharing of CNV datasets and provide momentum to build core centers for this niche within genomic medicine.

A Fast, Reproducible, High-throughput Variant Calling Workflow for Population Genomics.

The increasing availability of genomic resequencing data sets and high-quality reference genomes across the tree of life present exciting opportunities for comparative population genomic studies. However, substantial challenges prevent the simple reuse of data across different studies and species, arising from variability in variant calling pipelines, data quality, and the need for computationally intensive reanalysis. Here, we present snpArcher, a flexible and highly efficient workflow designed for the analysis of genomic resequencing data in nonmodel organisms. snpArcher provides a standardized variant calling pipeline and includes modules for variant quality control, data visualization, variant filtering, and other downstream analyses. Implemented in Snakemake, snpArcher is user-friendly, reproducible, and designed to be compatible with high-performance computing clusters and cloud environments. To demonstrate the flexibility of this pipeline, we applied snpArcher to 26 public resequencing data sets from nonmammalian vertebrates. These variant data sets are hosted publicly to enable future comparative population genomic analyses. With its extensibility and the availability of public data sets, snpArcher will contribute to a broader understanding of genetic variation across species by facilitating the rapid use and reuse of large genomic data sets.

Vcfdist: accurately benchmarking phased small variant calls in human genomes

Accurately benchmarking small variant calling accuracy is critical for the continued improvement of human whole genome sequencing. In this work, we show that current variant calling evaluations are biased towards certain variant representations and may misrepresent the relative performance of different variant calling pipelines. We propose solutions, first exploring the affine gap parameter design space for complex variant representation and suggesting a standard. Next, we present our tool vcfdist and demonstrate the importance of enforcing local phasing for evaluation accuracy. We then introduce the notion of partial credit for mostly-correct calls and present an algorithm for clustering dependent variants. Lastly, we motivate using alignment distance metrics to supplement precision-recall curves for understanding variant calling performance. We evaluate the performance of 64 phased Truth Challenge V2 submissions and show that vcfdist improves measured insertion and deletion performance consistency across variant representations from R2 = 0.97243 for baseline vcfeval to 0.99996 for vcfdist.

Benchmarking and improving the performance of variant-calling pipelines with RecallME.

The steady increment of Whole Genome/Exome sequencing and the development of novel Next Generation Sequencing-based gene panels requires continuous testing and validation of variant calling (VC) pipelines and the detection of sequencing-related issues to be maintained up-to-date and feasible for the clinical settings. State of the art tools are reliable when used to compute standard performance metrics. However, the need for an automated software to discriminate between bioinformatic and sequencing issues and to optimize VC parameters remains unmet. The aim of the current work is to present RecallME, a bioinformatic suite that tracks down difficult-to-detect variants as insertions and deletions in highly repetitive regions, thus providing the maximum reachable recall for both single nucleotide variants and small insertion and deletions and to precisely guide the user in the pipeline optimization process. Source code is freely available under MIT license at https://github.com/mazzalab-ieo/recallme. RecallME web application is available at https://translational-oncology-lab.shinyapps.io/recallme/. To use RecallME, users must obtain a license for ANNOVAR by themselves.

Direct Comparative Analysis of a Pharmacogenomics Panel with PacBio Hifi® Long-Read and Illumina Short-Read Sequencing.

Pharmacogenetics (PGx) aims to determine genetic signatures that can be used in clinical settings to individualize treatment for each patient, including anti-cancer drugs, anti-psychotics, and painkillers. Taken together, a better understanding of the impacts of genetic variants on the corresponding protein function or expression permits the prediction of the pharmacological response: responders, non-responders, and those with adverse drug reactions (ADRs). This work provides a comparison between innovative long-read sequencing (LRS) and short-read sequencing (SRS) techniques. The gene panel captured using PacBio HiFi® sequencing was tested on thirteen clinical samples on GENTYANE's platform. SRS, using a comprehensive pharmacogenetics panel, was performed in routine settings at the Civil Hospitals of Lyon. We focused on complex regions analysis, including copy number variations (CNVs), structural variants, repeated regions, and phasing-haplotyping for three key pharmacogenes: CYP2D6, UGT1A1, and NAT2. Variants and the corresponding expected star (*) alleles were reported. Although only 38.4% concordance was found for haplotype determination and 61.5% for diplotype, this did not affect the metabolism scoring. A better accuracy of LRS was obtained for the detection of the CYP2D6*5 haplotype in the presence of the duplicated wild-type CYP2D6*2 form. A total concordance was performed for UGT1A1 TA repeat detection. Direct phasing using the LRS approach allowed us to correct certain NAT2 profiles. Combining an optimized variant-calling pipeline and with direct phasing analysis, LRS is a robust technique for PGx analysis that can minimize the risk of mis-haplotyping.