Dataset Information

Speeding up eQTL scans in the BXD population using GPUs.

ABSTRACT: The BXD family of mouse strains are an important reference population for systems biology and genetics that have been fully sequenced and deeply phenotyped. To facilitate interactive use of genotype-phenotype relations using many massive omics data sets for this and other segregating populations, we have developed new algorithms and code that enable near-real-time whole-genome quantitative trait locus (QTL) scans for up to one million traits. By using easily parallelizable operations including matrix multiplication, vectorized operations, and element-wise operations, our method is more than 700 times faster than a R/qtl linear model genome scan using 16 threads. We used parallelization of different CPU threads as well as GPUs. We found that the speed advantage of GPUs is dependent on problem size and shape (the number of cases, number of genotypes, and number of traits). Our approach is ideal for interactive web services, such as GeneNetwork.org that need to display results in real-time. Our implementation is available as the Julia language package LiteQTL at https://github.com/senresearch/LiteQTL.jl.

SUBMITTER: Trotter C

PROVIDER: S-EPMC8664437 | biostudies-literature | 2021 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Speeding up eQTL scans in the BXD population using GPUs.

Trotter Chelsea C Kim Hyeonju H Farage Gregory G Prins Pjotr P Williams Robert W RW Broman Karl W KW Sen Śaunak Ś

G3 (Bethesda, Md.) 20211201 12

The BXD family of mouse strains are an important reference population for systems biology and genetics that have been fully sequenced and deeply phenotyped. To facilitate interactive use of genotype-phenotype relations using many massive omics data sets for this and other segregating populations, we have developed new algorithms and code that enable near-real-time whole-genome quantitative trait locus (QTL) scans for up to one million traits. By using easily parallelizable operations including m ...[more]

PMID: 34499130

Similar Datasets

Project description:BackgroundGene Set Enrichment Analysis (GSEA) is a popular method to reveal significant dependencies between predefined sets of gene symbols and observed phenotypes by evaluating the deviation of gene expression values between cases and controls. An established measure of inter-class deviation, the enrichment score, is usually computed using a weighted running sum statistic over the whole set of gene symbols. Due to the lack of analytic expressions the significance of enrichment scores is determined using a non-parametric estimation of their null distribution by permuting the phenotype labels of the probed patients. Accordingly, GSEA is a time-consuming task due to the large number of required permutations to accurately estimate the nominal p-value - a circumstance that is even more pronounced during multiple hypothesis testing since its estimate is lower-bounded by the inverse number of samples in permutation space.ResultsWe present rapidGSEA - a software suite consisting of two tools for facilitating permutation-based GSEA: cudaGSEA and ompGSEA. cudaGSEA is a CUDA-accelerated tool using fine-grained parallelization schemes on massively parallel architectures while ompGSEA is a coarse-grained multi-threaded tool for multi-core CPUs. Nominal p-value estimation of 4,725 gene sets on a data set consisting of 20,639 unique gene symbols and 200 patients (183 cases + 17 controls) each probing one million permutations takes 19 hours on a Xeon CPU and less than one hour on a GeForce Titan X GPU while the established GSEA tool from the Broad Institute (broadGSEA) takes roughly 13 days.ConclusioncudaGSEA outperforms broadGSEA by around two orders-of-magnitude on a single Tesla K40c or GeForce Titan X GPU. ompGSEA provides around one order-of-magnitude speedup to broadGSEA on a standard Xeon CPU. The rapidGSEA suite is open-source software and can be downloaded at https://github.com/gravitino/cudaGSEA as standalone application or package for the R framework.

Project description:BackgroundEfficient management of drug-resistant tuberculosis relies on fast diagnostics. To accelerate phenotypic drug susceptibility testing [pDST] for Mycobacterium tuberculosis [TB], we introduce TRACeR-TB, a test that infers drug resistance from antibiotic-specific mRNA biomarkers.MethodsTo develop TRACeR-TB, target genes were first identified through RNA sequencing experiments conducted on two drug-exposed, susceptible strains for four antitubercular drugs. Based on these findings, we designed drug-specific multiplex Quantigene panels to quantify mRNA levels of 8-9 biomarkers per drug (class), directly from crude cell lysates. The performance of TRACeR-TB was compared to the widely used Mycobacteria Growth Indicator Tube [MGIT] pDST by subjecting 238 strains with diverse drug resistance profiles to both methods, and aligning results to genotypic data. Furthermore, we explored TRACeR-TB's potential for evaluating molecules that enhance antibiotic efficacy, and investigated its applicability in macrophage models to assess Mtb's intracellular stress responses to drugs.FindingsAntituberculosis drugs trigger distinct transcriptional stress responses in susceptible, but not resistant bacilli, enabling a differentiation of the antibiotic phenotype in only 6 h. Validation on 238 strains showed TRACeR-TB had 100% (95% CI: 93·1-100%) sensitivity and 89·5% (95% CI: 74·7-97·2%) specificity compared to, respectively, 82·3% (95% CI: 69·2%-91·5%) and 94·8% (95% CI: 81·9%-99·4%) for MGIT pDST. TRACeR-TB specificity is likely underestimated due to the inclusion of isolates harbouring uncharacterised mutations. TRACeR-TB demonstrated 100% concordance with MGIT for drugs with reliable MGIT outcomes (moxifloxacin and isoniazid). Additionally, its sensitivity outperformed current rifampicin testing, detecting resistance in all borderline-resistant strains that MGIT missed, and bedaquiline testing. Furthermore, the assay detected the predicted effect of a novel drug booster and the intracellular drug-induced stress in macrophage models, highlighting its potential for drug optimisation.InterpretationTRACeR-TB is a complementary addition to current DSTs and can have a substantial impact on the TB diagnostics field. This tool can also play a vital role in identifying resistance mutations, thereby closing gaps in genotypic knowledge, and contribute to drug discovery and development.FundingInstitut Pasteur, Agence Nationale de la Recherche.

Dataset Information

Speeding up eQTL scans in the BXD population using GPUs.

Publications

Speeding up eQTL scans in the BXD population using GPUs.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets