Background

As human genome analysis technologies advance, research to identify target genes that cause diseases has been actively pursued. Genome-Wide Association Studies (GWAS) have succeeded in discovering numerous disease-related loci. However, because most variants are located in non-coding regions that do not encode proteins, it has been difficult to identify the causal genes that directly cause disease.

To overcome these limitations, Transcriptome-Wide Association Studies (TWAS) have become a key tool in statistical genetics. TWAS statistically models the correlation between gene expression levels and disease phenotypes by integrating expression Quantitative Trait Loci (eQTL) data. Researchers have used this to estimate genes that increase or decrease disease risk in specific tissues or cells.

Existing analysis pipelines possess a structural vulnerability: they fail to completely exclude the phenomenon of Linkage Disequilibrium (LD), where adjacent genetic variants are inherited together. This is because co-regulation occurs, where the expression profiles of surrounding genes change together due to linkage with nearby marker variants rather than the causal variant itself. Concerns have been steadily raised that false positive signals might be mistaken for disease causes or that the direction of causality might be reversed.

Key Findings

Researchers mathematically modeled and quantified the confounding effect of LD-driven co-regulation of expression profiles on TWAS results across cell types. They precisely disentangled the genetic structure shared by target genes harboring causal variants and nearby non-causal genes, applying it to both simulations and real genomic data.

Analysis revealed that co-regulation entangled with nearby tagging variants goes beyond simple statistical noise, fundamentally altering the direction of gene-disease associations. A reversal phenomenon was systematically observed, in which genes that actually suppress disease occurrence were incorrectly estimated as risk factors promoting disease, obscured by strong linkage disequilibrium (LD) signals at adjacent loci.

As expression data across tissues and cell types became more complexly entangled, the probability of such signal interference significantly increased. Particularly in complex polygenic diseases where multiple genes are clustered in a narrow genomic region, it was numerically proven that a non-negligible proportion of false correlation coefficients derived from co-regulation disturbance pass through existing correction models. The researchers clearly demonstrated, using statistical metrics, the limitation of existing TWAS algorithms that fail to incorporate both cell-specific transcriptional regulatory networks and local LD structures.

While previous methodologies focused primarily on expanding sample sizes or improving the predictive power of machine learning models, this study differentiates its approach by presenting a quantitative metric to separately measure the multi-pleiotropic signal entanglement caused by associated variants.

Significance and Outlook

This analysis challenges the prevailing dogma in transcriptome research regarding methods for prioritizing disease candidate genes. In the early stages of discovering numerous drug candidates, there have been many cases where development failed in later phases because the directionality suggested by TWAS was blindly trusted, leading to the development of therapeutics that inhibit or promote the functions of target genes. If the direction of disease association was reversed due to statistical illusion, massive resources and time would inevitably be wasted in the preclinical verification process.

The field of genomic epidemiology now faces the challenge of introducing a standard validation framework that verifies local variant structures alongside, rather than accepting transcriptome-wide association study (TWAS) results as the sole truth. It is essential to include a stage of physically validating statistical correlations by integrating high-resolution single-cell transcriptome data with functional genomics experiments based on CRISPR gene editing.

At the same time, existing multi-ancestry genomic reference panels remain limited in their ability to fully capture the complexity of tissue-specific eQTLs and differences in linkage disequilibrium (LD) structures across ancestries. Follow-up research must focus on significantly expanding data from non-European populations and developing multivariate statistical correction algorithms that integrate the transcriptional variability of the cellular microenvironment.