Evidence map›Paper›PMID 28249565›Full record

ArticleBMC bioinformatics2017

Statistical method to compare massive parallel sequencing pipelines.

M H Elsensohn, N Leblay, S Dimassi, A Campan-Fournier, A Labalme, F Roucher-Boulez, D Sanlaville, G Lesca, C Bardel, P Roy

Abstract read
In one paragraph

Article in BMC bioinformatics, 2017. The graph could read no effect estimate from its abstract, so it casts no vote on the map. Cited by 2 papers.

0numbers the graph read from it
0cells of the map it votes in
2citing papers in PubMed
–field-weighted citation impact
1 · What the graph read from it

What it found

Each row is one number read from the abstract, on the scale the paper reported it, with its interval. Left of the dashed line favours the treatment, right favours the comparator. Under each row is the sentence it came from. New to these charts? A ten-minute tutorial.

The abstract states no effect estimate the extractor could read, or names no intervention and outcome on the map, so this paper lights no cell and moves no belief. It is still indexed, cited and linked below.

2 · The registry

The trial behind it

Trials whose registry record cites this paper, or whose number appears in the abstract. A trial that started after this paper was published is citing it as background, not reporting it.

Neither the registry nor the abstract names a trial number. If this is a trial report, that itself is worth knowing.

3 · Its place in the literature

Who cites it

2 citing papers in PubMed.

  1. Article
  2. Review
4 · The record

Corrections and comments

PubMed lists nothing against this paper. Absence here is not a guarantee, only a check that was made.

5 · Who and what money

Authors and funding

10 authors.

M H ElsensohnService de Biostatistique-Bioinformatique, Hospices Civils de Lyon, 162 avenue Lacassagne, F-69003, Lyon, France. mad-helenie.elsensohn@chu-lyon.fr.
N LeblayService de Biostatistique-Bioinformatique, Hospices Civils de Lyon, 162 avenue Lacassagne, F-69003, Lyon, France.
S DimassiUniversité de Lyon, Lyon, France.
A Campan-FournierUniversité de Lyon, Lyon, France.
A LabalmeService de Génétique, Hospices Civils de Lyon, Lyon, France.
F Roucher-BoulezUniversité de Lyon, Lyon, France.
D SanlavilleUniversité de Lyon, Lyon, France.
G LescaUniversité de Lyon, Lyon, France.
C BardelService de Biostatistique-Bioinformatique, Hospices Civils de Lyon, 162 avenue Lacassagne, F-69003, Lyon, France.
P RoyService de Biostatistique-Bioinformatique, Hospices Civils de Lyon, 162 avenue Lacassagne, F-69003, Lyon, France.

Funding

No grant is acknowledged in the PubMed record.

6 · The paper itself

Abstract

backgroundToday, sequencing is frequently carried out by Massive Parallel Sequencing (MPS) that cuts drastically sequencing time and expenses. Nevertheless, Sanger sequencing remains the main validation method to confirm the presence of variants. The analysis of MPS data involves the development of several bioinformatic tools, academic or commercial. We present here a statistical method to compare MPS pipelines and test it in a comparison between an academic (BWA-GATK) and a commercial pipeline (TMAP-NextGENe®), with and without reference to a gold standard (here, Sanger sequencing), on a panel of 41 genes in 43 epileptic patients. This method used the number of variants to fit log-linear models for pairwise agreements between pipelines. To assess the heterogeneity of the margins and the odds ratios of agreement, four log-linear models were used: a full model, a homogeneous-margin model, a model with single odds ratio for all patients, and a model with single intercept. Then a log-linear mixed model was fitted considering the biological variability as a random effect.

resultsAmong the 390,339 base-pairs sequenced, TMAP-NextGENe® and BWA-GATK found, on average, 2253.49 and 1857.14 variants (single nucleotide variants and indels), respectively. Against the gold standard, the pipelines had similar sensitivities (63.47% vs. 63.42%) and close but significantly different specificities (99.57% vs. 99.65%; p < 0.001). Same-trend results were obtained when only single nucleotide variants were considered (99.98% specificity and 76.81% sensitivity for both pipelines).

conclusionsThe method allows thus pipeline comparison and selection. It is generalizable to all types of MPS data and all pipelines.

Indexed as

Models, StatisticalComputational BiologyEpilepsyHigh-Throughput Nucleotide SequencingHumansINDEL MutationPolymorphism, Single NucleotideSequence Analysis, DNAMassive parallel sequencingNext-generation sequencingPipeline comparisonSensitivitySpecificityStatistical methods

Identifiers

PMID28249565
PMCPMC5333416

What OpenQuestion holds

Textmetadata
LicenceCC BY
Read underepoch 390

Registered trials

None linked

Read under generation 80e0d062 · epoch 390. Bibliography from PubMed, PubMed Central and OpenAlex; grants from NIH RePORTER; trial links from ClinicalTrials.gov; estimates, votes and beliefs from the OpenQuestion graph.