Biostatistical Aspects of Whole Genome Sequencing Studies: Preprocessing and Quality Control.

Rapid advances in high-throughput DNA sequencing technologies have enabled large-scale whole genome sequencing (WGS) studies. Before performing association analysis between phenotypes and genotypes, preprocessing and quality control (QC) of the raw sequence data need to be performed. Because many biostatisticians have not been working with WGS data so far, we first sketch Illumina's short-read sequencing technology. Second, we explain the general preprocessing pipeline for WGS studies. Third, we provide an overview of important QC metrics, which are applied to WGS data: on the raw data, after mapping and alignment, after variant calling, and after multisample variant calling. Fourth, we illustrate the QC with the data from the GENEtic SequencIng Study Hamburg-Davos (GENESIS-HD), a study involving more than 9000 human whole genomes. All samples were sequenced on an Illumina NovaSeq 6000 with an average coverage of 35× using a PCR-free protocol. For QC, one genome in a bottle (GIAB) trio was sequenced in four replicates, and one GIAB sample was successfully sequenced 70 times in different runs. Fifth, we provide empirical data on the compression of raw data using the DRAGEN original read archive (ORA). The most important quality metrics in the application were genetic similarity, sample cross-contamination, deviations from the expected Het/Hom ratio, relatedness, and coverage. The compression ratio of the raw files using DRAGEN ORA was 5.6:1, and compression time was linear by genome coverage. In summary, the preprocessing, joint calling, and QC of large WGS studies are feasible within a reasonable time, and efficient QC procedures are readily available.

Cost-effectiveness of a telemonitoring programme in patients with cardiovascular diseases compared with standard of care.

OBJECTIVES: The main aim of this work was to analyse the cost-effectiveness of an integrated care concept (NICC) that combines telemonitoring with the support of a care centre in addition to guideline therapy for patients. Secondary aims were to compare health utility and health-related quality of life (QoL) between NICC and standard of care (SoC). METHODS: The randomised controlled CardioCare MV Trial compared NICC and SoC in patients from Mecklenburg-West Pomerania (Germany) with atrial fibrillation, heart failure or treatment-resistant hypertension. QoL was measured using the EQ-5D-5L at baseline, 6 months and 1 year follow-up. Quality-adjusted life years (QALYs), EQ5D utility scores, Visual Analogue Scale (VAS) Scores and VAS adjusted life years (VAS-AL) were calculated. Cost data were obtained from health insurance companies, and the payer perspective was taken in health economic analyses. Quantile regression was used with adjustments for stratification variables. RESULTS: The net benefit of NICC (QALY) was 0.031 (95% CI 0.012 to 0.050; p=0.001) in this trial involving 957 patients. EQ5D Index values, VAS-ALs and VAS were larger for NICC compared with SoC at 1 year follow-up (all p≤0.004). Direct cost per patient and year were €323 (CI €157 to €489) lower in the NICC group. When 2000 patients are served by the care centre, NICC is cost-effective if one is willing to pay €10 652 per QALY per year. CONCLUSION: NICC was associated with higher QoL and health utility. The programme is cost-effective if one is willing to pay approximately €11 000 per QALY per year.

Comparison of calling pipelines for whole genome sequencing: an empirical study demonstrating the importance of mapping and alignment.

Rapid advances in high-throughput DNA sequencing technologies have enabled the conduct of whole genome sequencing (WGS) studies, and several bioinformatics pipelines have become available. The aim of this study was the comparison of 6 WGS data pre-processing pipelines, involving two mapping and alignment approaches (GATK utilizing BWA-MEM2 2.2.1, and DRAGEN 3.8.4) and three variant calling pipelines (GATK 4.2.4.1, DRAGEN 3.8.4 and DeepVariant 1.1.0). We sequenced one genome in a bottle (GIAB) sample 70 times in different runs, and one GIAB trio in triplicate. The truth set of the GIABs was used for comparison, and performance was assessed by computation time, F1 score, precision, and recall. In the mapping and alignment step, the DRAGEN pipeline was faster than the GATK with BWA-MEM2 pipeline. DRAGEN showed systematically higher F1 score, precision, and recall values than GATK for single nucleotide variations (SNVs) and Indels in simple-to-map, complex-to-map, coding and non-coding regions. In the variant calling step, DRAGEN was fastest. In terms of accuracy, DRAGEN and DeepVariant performed similarly and both superior to GATK, with slight advantages for DRAGEN for Indels and for DeepVariant for SNVs. The DRAGEN pipeline showed the lowest Mendelian inheritance error fraction for the GIAB trios. Mapping and alignment played a key role in variant calling of WGS, with the DRAGEN outperforming GATK.