A robust association test leveraging unknown genetic interactions: Application to cystic fibrosis lung disease
Bibliographic record
Abstract
Abstract For complex traits such as lung disease in Cystic Fibrosis (CF), Gene x Gene or Gene x Environment interactions can impact disease severity but these remain largely unknown. Unaccounted-for genetic interactions introduce a distributional shift in the quantitative trait across the genotypic groups. Joint location and scale tests, or full distributional differences across genotype groups can account for unknown genetic interactions and increase power for gene identification compared with the conventional association test. Here we propose a new joint location and scale test (JLS), a quantile regression-basd JLS (qJLS), that addresses previous limitations. Specifically, qJLS is free of distributional assumptions, thus applies to non-Gaussian traits; is as powerful as the existing JLS tests under Gaussian traits; and is computationally efficient for genome-wide association studies (GWAS). Our simulation studies, which model unknown genetic interactions, demonstrate that qJLS is robust to skewed and heavy-tailed error distributions and is as powerful as other JLS tests in the literature under normality. Without any unknown genetic interaction, qJLS shows a large increase in power with non-Gaussian traits over conventional association tests and is slightly less powerful under normality. We apply the qJLS method to the Canadian CF Gene Modifier Study (n=1,997) and identified a genome-wide significant variant, rs9513900 on chromosome 13, that had not previously been reported to contribute to CF lung disease. qJLS provides a powerful alternative to conventional genetic association tests, where interactions my contribute to a quantitative trait. Author summary Cystic fibrosis (CF) is a genetic disorder caused by loss-of-function variants in CF transmembrane conductance regulator ( CFTR ) gene, leading to disease in several organs and notably the lungs. Even among those who share identical CF causing variants, their lung disease severity is variable, which is presumed to be caused in part by other genes besides CFTR referred to as modifier genes. Several genome-wide association studies of CF lung disease have identified associated loci but these account for only a small fraction of the total CF lung disease heritability. This may be due to other environmental factors such as infections, smoke exposure, socioeconomic status, treatment of lung diseases or a numerous other unknown or unmeasured factors that may interact with modifier genes. A class of new statistical methods can leverage these unknown interactions to better detect putative genetic loci. We provide a comprehensive simulation study that incorporates unknown interactions and we show that these statistical methods perform better than conventional approaches at identifying contributing genetic loci when the assumptions for these approaches are met. We then develop an approach that is robust to the typical normal assumptions, provide software for implementation and we apply it to the Canadian CF Gene Modifier Study to identify novel variants contributing to CF lung disease.
Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.
How this classification was reachedexpand
Full frame machine prediction
Teacher imitationNot calibrated prevalence, not ground truth. Human validation pending. The Gemma side is a direct model label for every work in the frame, read from the title-only record. The Codex side is a classifier learned from the 10,348 direct Codex labels and calibrated to design-weighted sample rates; fields without enough sample support carry no Codex call. Candidate is the union of the two sides; consensus is their intersection. These outputs are machine_predicted_unvalidated and are not human labels.
Distilled classifier scores by category (both heads)
| Category | Codex | Gemma |
|---|---|---|
| Metaresearch | 0.019 | 0.063 |
| Meta-epidemiology (narrow) | 0.001 | 0.000 |
| Meta-epidemiology (broad) | 0.001 | 0.001 |
| Bibliometrics | 0.001 | 0.002 |
| Science and technology studies | 0.001 | 0.002 |
| Scholarly communication | 0.001 | 0.001 |
| Open science | 0.002 | 0.002 |
| Research integrity | 0.001 | 0.001 |
| Insufficient payload (model declined to judge) | 0.002 | 0.000 |
Machine scores (provisional)
The two teacher heads of the student model, read on this work. A score orders the frame for review; it never asserts a category, and the validation status ships verbatim with every row.
Baseline scores from an immature model (maturity gate not passed, 7 training rounds). Scores rank; they never assert a category.
score_only:v0-immature-baseline · verbatim from the scoring run: score_only means the number may rank works, and no category label ships from itClassification
machine, unvalidatedMachine predicted; a candidate call from one source (direct Gemma or distilled Codex), not a consensus.
How this classification was reached, model by model and score by score, is at the end of the page under "How this classification was reached".