Abstract:Replicable signals from different yet conceptually related studies provide stronger scientific evidence and more powerful inference. We introduce STAREG, a statistical method for replicability analysis of high throughput experiments, and apply it to analyze spatial transcriptomic studies. STAREG uses summary statistics from multiple studies of high throughput experiments and models the the joint distribution of p -values accounting for the heterogeneity of different studies. It effectively controls the false discovery rate (FDR) and has higher power by information borrowing. Moreover, it provides different rankings of important genes. With the EM algorithm in combination with pool-adjacent-violator-algorithm (PAVA), STAREG is scalable to datasets with millions of genes without any tuning parameters. Analyzing two pairs of spatially resolved transcriptomic datasets, we are able to make biological discoveries that otherwise cannot be obtained by using existing methods. Irreplicable research wastes time, money, and/or resources. Approximately $28 billion is estimated to be spent on preclinical research that cannot be replicated every year in the United States alone. Possible causes of irreplicable research may include experimental design, laboratory practices, and data analysis. We focus on data analysis. The past two decades have witnessed the expansion and increased availability of genomic data from high-throughput experiments. Due to privacy concerns or logistic reasons, raw data can be difficult to access but summary data such as p -values are readily available. We introduce STAREG, which jointly analyzes p -values from multiple genomic datasets that target the same scientific question with different populations or different technologies. This allows us to have more convincing and robust findings. STAREG is computationally scalable with solid statistical analysis. Moreover, it is versatile, platform-independent, and only requires p -values as input. By analyzing data sets from spatially resolved transcriptomic studies, we make biological discoveries that otherwise cannot be obtained with existing methods.

A powerful approach to identify replicable variants in genome-wide association studies

Replicability analysis for genome-wide association studies

Estimating Reproducibility in Genome-Wide Association Studies

Technical reproducibility of genotyping SNP arrays used in genome-wide association studies.

Robust Reference Powered Association Test of Genome-Wide Association Studies

Replication of Known and Identification of Novel Associations in Biobank-Scale Datasets: A Survey Using UK Biobank and FinnGen

Validating GWAS Findings through Reverse Engineering of Contingency Tables

A robust statistical method for Genome-wide association analysis of human copy number variation

STAREG: Statistical replicability analysis of high throughput experiments with applications to spatial transcriptomic studies

Assessing reproducibility of inherited variants detected with short-read whole genome sequencing

Replicability analysis of high dimensional data accounting for dependence

Genome-wide interaction-based association analysis identified multiple new susceptibility Loci for common diseases.

Nonreplication in Genetic Studies of Complex Diseases—Lessons Learned from Studies of Osteoporosis and Tentative Remedies

Beyond guilty by association at scale: searching for causal variants on the basis of genome-wide summary statistics

JUMP: replicability analysis of high-throughput experiments with applications to spatial transcriptomic studies

Pitfall of genome-wide association studies: Sources of inconsistency in genotypes and their effects

Extracting replicable associations across multiple studies: algorithms for controlling the false discovery rate

Statistical analysis for genome-wide association study

A Novel Efficient Algorithm for Common Variants Genotyping from Low-Coverage Sequencing Data

A nonparametric test for association with multiple loci in the retrospective case-control study.

Powerful gene-based testing by integrating long-range chromatin interactions and knockoff genotypes