Optimizing expression quantitative trait locus mapping workflows for single-cell studies.
Ontology highlight
ABSTRACT: Background
Single-cell RNA sequencing (scRNA-seq) has enabled the unbiased, high-throughput quantification of gene expression specific to cell types and states. With the cost of scRNA-seq decreasing and techniques for sample multiplexing improving, population-scale scRNA-seq, and thus single-cell expression quantitative trait locus (sc-eQTL) mapping, is increasingly feasible. Mapping of sc-eQTL provides additional resolution to study the regulatory role of common genetic variants on gene expression across a plethora of cell types and states and promises to improve our understanding of genetic regulation across tissues in both health and disease.Results
While previously established methods for bulk eQTL mapping can, in principle, be applied to sc-eQTL mapping, there are a number of open questions about how best to process scRNA-seq data and adapt bulk methods to optimize sc-eQTL mapping. Here, we evaluate the role of different normalization and aggregation strategies, covariate adjustment techniques, and multiple testing correction methods to establish best practice guidelines. We use both real and simulated datasets across single-cell technologies to systematically assess the impact of these different statistical approaches.Conclusion
We provide recommendations for future single-cell eQTL studies that can yield up to twice as many eQTL discoveries as default approaches ported from bulk studies.
SUBMITTER: Cuomo ASE
PROVIDER: S-EPMC8223300 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Coordinated regulation of gene expression levels across a series of experimental conditions provides valuable information about the functions of correlated transcripts. The consideration of gene expression correlation over a time or tissue dimension has proved valuable in predicting gene function. Here, we consider correlations over a genetic dimension. In addition to identifying coregulated genes, the genetic dimension also supplies us with information about the genomic locations of putative regulatory loci. We calculated correlations among approximately 45,000 expression traits derived from 60 individuals in an F2 sample segregating for obesity and diabetes. By combining the correlation results with linkage mapping information, we were able to identify regulatory networks, make functional predictions for uncharacterized genes, and characterize novel members of known pathways. We found evidence of coordinate regulation of 174 G protein-coupled receptor protein signaling pathway expression traits. Of the 174 traits, 50 had their major LOD peak within 10 cM of a locus on Chromosome 2, and 81 others had a secondary peak in this region. We also characterized a Riken cDNA clone that showed strong correlation with stearoyl-CoA desaturase 1 expression. Experimental validation confirmed that this clone is involved in the regulation of lipid metabolism. We conclude that trait correlation combined with linkage mapping can reveal regulatory networks that would otherwise be missed if we studied only mRNA traits with statistically significant linkages in this small cross. The combined analysis is more sensitive compared with linkage mapping alone.
| S-EPMC1331977 | biostudies-literature
Project description:Genome-wide association study (GWAS) can be used to identify genes that increase the risk of psychiatric diseases. However, much of the disease heritability is still unexplained, suggesting that there are genes to be discovered. Functional annotation of the genetic variants may increase the power of GWAS to identify disease genes, by providing prior information that can be used in Bayesian analysis or in reducing the number of tests. Expression quantitative trait loci (eQTLs) are genomic loci that regulate gene expression. Genetic mapping of eQTLs can help reveal novel functional effects of thousands of single nucleotide polymorphisms (SNPs). The present review mainly focused on the current knowledge on brain eQTL mapping, and discussed some major methodological issues and their possible solutions. The frequently ignored problems of batch effects, covariates, and multiple testing were emphasized, since they can lead to false positives and false negatives. The future application of eQTL data in GWAS analysis was also discussed.
| S-EPMC3074249 | biostudies-other
Project description:Expression quantitative trait loci (eQTL) can provide a link between disease susceptibility variants discovered by genetic association studies and biology. To date, eQTL mapping studies have been primarily conducted in healthy individuals from population-based cohorts. Genetic effects have been known to be context-specific and vary with changing environmental stimuli. We conducted a transcriptome- and genome-wide eQTL mapping study in a cohort of patients with idiopathic or heritable pulmonary arterial hypertension (PAH) using RNA sequencing (RNAseq) data from whole blood. We sought confirmation from three published population-based eQTL studies, including the GTEx Project, and followed up potentially novel eQTL not observed in the general population. In total, we identified 2314 eQTL of which 90% were cis-acting and 75% were confirmed by at least one of the published studies. While we observed a higher GWAS trait colocalization rate among confirmed eQTL, colocalisation rate of novel eQTL reported for lung-related phenotypes was twice as high as that of confirmed eQTL. Functional enrichment analysis of genes with novel eQTL in PAH highlighted immune-related processes, a suspected contributor to PAH. These potentially novel eQTL specific to or active in PAH could be useful in understanding genetic risk factors for other diseases that share common mechanisms with PAH.
| S-EPMC7690609 | biostudies-literature
Project description:Coordinated regulation of gene expression levels across a series of experimental conditions provides valuable information about the functions of correlated transcripts. To map gene regulatory pathways, we used microarray-derived gene expression measurements in 60 individuals of an F2 sample segregating for diabetes. We performed correlation analysis among ~40,000 expression traits. By combining correlation among expression traits and linkage mapping information, we were able to identify regulatory networks, make functional predictions to uncharacterized genes, and characterize novel members of known pathways. Using 36 seed traits, we found evidence of coordinate regulation of 160 G-protein coupled receptor (GPCR) pathway expression traits. Of the 160 traits, 50 had their major LOD peak within 8 cM of a locus on chromosome 2, and 81 others had a secondary peak in this region. A previously uncharacterized Riken cDNA clone, which showed strong correlation with stearoyl CoA desaturase 1 expression, was experimentally validated to be responsive to conditions that regulate lipid metabolism. Using linkage mapping, we identified multiple genes whose expression is under the control of transcription regulatory loci. Trait-correlation combined with linkage mapping can reveal regulatory networks that would otherwise be missed if we only studied mRNA traits with statistically significant linkages in this small cross. The combined analysis is more sensitive compared with linkage mapping only. References: Kendziorski C., M. Chen, M. Yuan, H. Lan, and A.D. Attie. Statistical Methods for Expression Quantitative Trait Loci (eQTL) Mapping. Biometrics, to appear, 2005. Lan H, Chen M, Flowers JB, Yandell BS, Stapleton DS, et al. (2006) Combined Expression Trait Correlations and Expression Quantitative Trait Locus Mapping. PLoS Genet 2(1): e6. Keywords: Genetics of gene expression
2005-09-21 | GSE3330 | GEO
Project description:Power calculation prior to a genetic experiment can help investigators choose the optimal sample size to detect a quantitative trait locus (QTL). Without the guidance of power analysis, an experiment may be underpowered or overpowered. Either way will result in wasted resource. QTL mapping and genome-wide association studies (GWAS) are often conducted using a linear mixed model (LMM) with controls of population structure and polygenic background using markers of the whole genome. Power analysis for such a mixed model is often conducted via Monte Carlo simulations. In this study, we derived a non-centrality parameter for the Wald test statistic for association, which allows analytical power analysis. We show that large samples are not necessary to detect a biologically meaningful QTL, say explaining 5% of the phenotypic variance. Several R functions are provided so that users can perform power analysis to determine the minimum sample size required to detect a given QTL with a certain statistical power or calculate the statistical power with given sample size and known values of other population parameters.
| S-EPMC6781134 | biostudies-literature
Project description:Genetic mapping studies in the mouse and other model organisms are used to search for genes underlying complex phenotypes. Traditional genetic mapping studies that employ single-generation crosses have poor mapping resolution and limit discovery to loci that are polymorphic between the two parental strains. Multiparent outbreeding populations address these shortcomings by increasing the density of recombination events and introducing allelic variants from multiple founder strains. However, multiparent crosses present new analytical challenges and require specialized software to take full advantage of these benefits. Each animal in an outbreeding population is genetically unique and must be genotyped using a high-density marker set; regression models for mapping must accommodate multiple founder alleles, and complex breeding designs give rise to polygenic covariance among related animals that must be accounted for in mapping analysis. The Diversity Outbred (DO) mice combine the genetic diversity of eight founder strains in a multigenerational breeding design that has been maintained for >16 generations. The large population size and randomized mating ensure the long-term genetic stability of this population. We present a complete analytical pipeline for genetic mapping in DO mice, including algorithms for probabilistic reconstruction of founder haplotypes from genotyping array intensity data, and mapping methods that accommodate multiple founder haplotypes and account for relatedness among animals. Power analysis suggests that studies with as few as 200 DO mice can detect loci with large effects, but loci that account for <5% of trait variance may require a sample size of up to 1000 animals. The methods described here are implemented in the freely available R package DOQTL.
| S-EPMC4169154 | biostudies-literature
Project description:Time to flowering and maturity in soybean is controlled by loci E1 to E5, and E7 to E9. These loci were assigned to molecular linkage groups (MLGs) except for E5. This study was conducted to map the E5 locus using F2 populations expected to segregate for E5. F2 populations were subjected to quantitative trait locus (QTL) analysis for days to flowering (DF) and maturity (DM). In Harosoy-E5 × Clark-e2 population, QTLs for DF and DM were found at a similar position with E2. In Harosoy × Clark-e2E5 population, QTLs for DF and DM were found in MLG D1a and B1, respectively. In Harosoy-E5Dt2 × Clark-e2 population, a QTL for DF was found in MLG B1. Thus, results from these populations were not fully consistent, and no candidate QTL for E5 was found. In Harosoy × PI 80837 population, from which E5 was originally identified, QTLs corresponding to E1 and E3 were found, but none for E5 existed. Harosoy and PI 80837 had the e2-ns allele whereas Harosoy-E5 had the E2-dl allele. The E2-dl allele of Harosoy-E5 may have been generated by outcrossing and may be responsible for the lateness of Harosoy-E5. We conclude that a unique E5 gene may not exist.
| S-EPMC4902463 | biostudies-literature
Project description:Fertility traits in humans are heritable, however, little is known about the genes that influence reproductive outcomes or the genetic variants that contribute to differences in these traits between individuals, particularly women. To address this gap in knowledge, we performed an unbiased genome-wide expression quantitative trait locus (eQTL) mapping study to identify common regulatory (expression) single nucleotide polymorphisms (eSNPs) in mid-secretory endometrium. We identified 423 cis-eQTLs for 132 genes that were significant at a false discovery rate (FDR) of 1%. After pruning for strong LD (r2 >0.95), we tested for associations between eSNPs and fecundability (the ability to get pregnant), measured as the length of the interval to pregnancy, in 117 women. Two eSNPs were associated with fecundability at a FDR of 5%; both were in the HLA region and were eQTLs for the TAP2 gene (P = 1.3x10-4) and the HLA-F gene (P = 4.0x10-4), respectively. The effects of these SNPs on fecundability were replicated in an independent sample. The two eSNPs reside within or near regulatory elements in decidualized human endometrial stromal cells. Our study integrating eQTL mapping in a primary tissue with association studies of a related phenotype revealed novel genes and associated alleles with independent effects on fecundability, and identified a central role for two HLA region genes in human implantation success.
| S-EPMC4957750 | biostudies-literature
Project description:In multiple quantitative trait locus (QTL) mapping, a high-dimensional sparse regression model is usually employed to account for possible multiple linked QTLs. The QTL model may include closely linked and thus highly correlated genetic markers, especially when high-density marker maps are used in QTL mapping because of the advancement in sequencing technology. Although existing algorithms, such as Lasso, empirical Bayesian Lasso (EBlasso) and elastic net (EN) are available to infer such QTL models, more powerful methods are highly desirable to detect more QTLs in the presence of correlated QTLs. We developed a novel empirical Bayesian EN (EBEN) algorithm for multiple QTL mapping that inherits the efficiency of our previously developed EBlasso algorithm. Simulation results demonstrated that EBEN provided higher power of detection and almost the same false discovery rate compared with EN and EBlasso. Particularly, EBEN can identify correlated QTLs that the other two algorithms may fail to identify. When analyzing a real dataset, EBEN detected more effects than EN and EBlasso. EBEN provides a useful tool for inferring high-dimensional sparse model in multiple QTL mapping and other applications. An R software package 'EBEN' implementing the EBEN algorithm is available on the Comprehensive R Archive Network (CRAN).
| S-EPMC4815589 | biostudies-literature
Project description:BACKGROUND: Etiology of complex disorders, such as cataract and neurodegenerative diseases including age-related macular degeneration (AMD), remains poorly understood due to the paucity of animal models, fully replicating the human disease. Previously, two quantitative trait loci (QTLs) associated with early cataract, AMD-like retinopathy, and some behavioral aberrations in senescence-accelerated OXYS rats were uncovered on chromosome 1 in a cross between OXYS and WAG rats. To confirm the findings, we generated interval-specific congenic strains, WAG/OXYS-1.1 and WAG/OXYS-1.2, carrying OXYS-derived loci of chromosome 1 in the WAG strain. Both congenic strains displayed early cataract and retinopathy but differed clinically from OXYS rats. Here we applied a high-throughput RNA sequencing (RNA-Seq) strategy to facilitate nomination of the candidate genes and functional pathways that may be responsible for these differences and can contribute to the development of the senescence-accelerated phenotype of OXYS rats. RESULTS: First, the size and map position of QTL-derived congenic segments were determined by comparative analysis of coding single-nucleotide polymorphisms (SNPs), which were identified for OXYS, WAG, and congenic retinal RNAs after sequencing. The transferred locus was not what we expected in WAG/OXYS-1.1 rats. In rat retina, 15442 genes were expressed. Coherent sets of differentially expressed genes were identified when we compared RNA-Seq retinal profiles of 20-day-old WAG/OXYS-1.1, WAG/OXYS-1.2, and OXYS rats. The genes most different in the average expression level between the congenic strains included those generally associated with the Wnt, integrin, and TGF-? signaling pathways, widely involved in neurodegenerative processes. Several candidate genes (including Arhgap33, Cebpg, Gtf3c1, Snurf, Tnfaip3, Yme1l1, Cbs, Car9 and Fn1) were found to be either polymorphic in the congenic loci or differentially expressed between the strains. These genes may contribute to the development of cataract and retinopathy. CONCLUSIONS: This study is the first RNA-Seq analysis of the rat retinal transcriptome generated with 40 mln sequencing read depth. The integration of QTL and transcriptomic analyses in our study forms the basis of future research into the relationship between the candidate genes within the congenic regions and specific changes in the retinal transcriptome as possible causal mechanisms that underlie age-associated disorders.
| S-EPMC4303943 | biostudies-literature