Abstract
GWASs have successfully identified numerous genetic variants linked to complex traits, but traditional univariate approaches often fail to capture shared genetic architecture across multiple phenotypes. As the scale of genomic data continues to increase, the demand for more efficient multi-phenotype analysis methods has become particularly critical. In addition, the issue of population structure must also be properly addressed to ensure robust and unbiased results. Multivariate methods for multi-phenotype analysis, such as GAMMA, address this by combining linear mixed models with multivariate distance matrix regression to account for population structure; however, since these methods utilize computationally intensive models, developing efficient implementations is essential for practical analysis. Although GAMMA is a well-designed and effective tool, its original implementation relies on multiple programming environments and requires frequent data exchanges between components. These factors increase computational burden and complicate installation and execution for users unfamiliar with programming, making practical applications, particularly for high-dimensional datasets, challenging. Here, we present GAMMA-RAY, a high-performance C++ implementation that streamlines the computational pipeline, leverages parallel processing, and employs efficient matrix operations to achieve substantial reductions in runtime and memory usage. GAMMA-RAY provides both a user-friendly web-based interface for non-programmers and a standalone version for secure local execution. We further applied GAMMA-RAY to a yeast dataset and identified putative trans-eQTLs, in which several variants overlapped with previously reported cis- and trans-eQTLs. In addition, functional enrichment analysis revealed that the associated trans-eGenes are enriched, a conclusion consistently supported by biological annotation resources and underscoring the biological significance of these results.