A robust biostatistical method leverages informative but uncertainly determined qPCR data for biomarker detection, early diagnosis, and treatment

• Published in: PloS one Vol. 17; no. 1; p. e0263070
• Main Authors: Zhuang, Wei, Camacho, Luísa, Silva, Camila S., Thomson, Michael, Snyder, Kevin
• Format: Journal Article
• Language: English
• Published: United States Public Library of Science 31.01.2022; Public Library of Science (PLoS)
• Subjects: Biological materials; Biology and life sciences; Biomarkers; Carbon monoxide; Data analysis; Decision making; Diagnosis; Early Diagnosis; Electric power distribution; Gene expression; Genetic markers; Genetic screening; Health aspects; Humans; Identification and classification; Management; Medical diagnosis; Methods; MicroRNAs; Models, Theoretical; Normal distribution; Nucleic acids; Physical Sciences; Polymerase chain reaction; Real-Time Polymerase Chain Reaction; Research and Analysis Methods; Robust statistics; Robustness; Samples; Skewed distributions; Software; Statistical analysis; Statistical methods; Statistics; Toxicology; United States; United States > US; Arkansas; Maryland
• ISSN: 1932-6203; 1932-6203
• DOI: 10.1371/journal.pone.0263070

As a common medium-throughput technique, qPCR (quantitative real-time polymerase chain reaction) is widely used to measure levels of nucleic acids. In addition to accurate and complete data, experimenters have unavoidably observed some incomplete and uncertainly determined qPCR data because of intrinsically low overall amounts of biological materials, such as nucleic acids present in biofluids. When there are samples with uncertainly determined qPCR data, some investigators apply the statistical complete-case method by excluding the subset of samples with uncertainly determined data from analysis (CO), while others simply choose not to analyze (CNA) these datasets altogether. To include as many observations as possible in analysis for interesting differential changes between groups, some investigators set incomplete observations equal to the maximum quality qPCR cycle (MC), such as 32 and 40. Although straightforward, these methods may decrease the sample size, skew the data distribution, and compromise statistical power and research reproducibility across replicate qPCR studies. To overcome the shortcomings of the existing, commonly-used qPCR data analysis methods and to join the efforts in advancing statistical analysis in rigorous preclinical research, we propose a robust nonparametric statistical cycle-to-threshold method (CTOT) to analyze incomplete qPCR data for two-group comparisons. CTOT incorporates important characteristics of qPCR data and time-to-event statistical methodology, resulting in a novel analytical method for qPCR data that is built around good quality data from all subjects, certainly determined or not. Considering the benchmark full data (BFD), we compared the abilities of CTOT, CO, MC, and CNA statistical methods to detect interesting differential changes between groups with informative but uncertainly determined qPCR data. Our simulations and applications show that CTOT improves the power of detecting and confirming differential changes in many situations over the three commonly used methods without excess type I errors. The robust nonparametric statistical method of CTOT helps leverage qPCR technology and increase the power to detect differential changes that may assist decision making with respect to biomarker detection and early diagnosis, with the goal of improving the management of patient healthcare.