Dataset Information

Comparison of quantitative trait loci methods: Total expression and allelic imbalance method in brain RNA-seq.

ABSTRACT:

Background

Of the 108 Schizophrenia (SZ) risk-loci discovered through genome-wide association studies (GWAS), 96 are not altering the sequence of any protein. Evidence linking non-coding risk-SNPs and genes may be established using expression quantitative trait loci (eQTL). However, other approaches such allelic expression quantitative trait loci (aeQTL) also may be of use.

Methods

We applied both the eQTL and aeQTL analysis to a biobank of deeply sequenced RNA from 680 dorso-lateral pre-frontal cortex (DLPFC) samples. For each of 340 genes proximal to the SZ risk-SNPs, we asked how much SNP-genotype affected total expression (eQTL), as well as how much the expression ratio between the two alleles differed from 1:1 as a consequence of the risk-SNP genotype (aeQTL).

Results

We analyzed overlap with comparable eQTL-findings: 16 of the 30 risk-SNPs known to have gene-level eQTL also had gene-level aeQTL effects. 6 of 21 risk-SNPs with known splice-eQTL had exon-aeQTL effects. 12 novel potential risk genes were identified with the aeQTL approach, while 55 tested SNP-pairs were found as eQTL but not aeQTL. Of the tested 108 loci we could find at least one gene to be associated with 21 of the risk-SNPs using gene-level aeQTL, and with an additional 18 risk-SNPs using exon-level aeQTL.

Conclusion

Our results suggest that the aeQTL strategy complements the eQTL approach to susceptibility gene identification.

SUBMITTER: Gadin JR

PROVIDER: S-EPMC6576752 | biostudies-literature | 2019

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Comparison of quantitative trait loci methods: Total expression and allelic imbalance method in brain RNA-seq.

Gådin Jesper R JR Buil Alfonso A Colantuoni Carlo C Jaffe Andrew E AE Nielsen Jacob J Shin Joo-Heon JH Hyde Thomas M TM Kleinman Joel E JE Plath Niels N Eriksson Per P Brunak Søren S Didriksen Michael M Weinberger Daniel R DR Folkersen Lasse L

PloS one 20190617 6

<h4>Background</h4>Of the 108 Schizophrenia (SZ) risk-loci discovered through genome-wide association studies (GWAS), 96 are not altering the sequence of any protein. Evidence linking non-coding risk-SNPs and genes may be established using expression quantitative trait loci (eQTL). However, other approaches such allelic expression quantitative trait loci (aeQTL) also may be of use.<h4>Methods</h4>We applied both the eQTL and aeQTL analysis to a biobank of deeply sequenced RNA from 680 dorso-late ...[more]

PMID: 31206532

Similar Datasets

Project description:BackgroundOne method of identifying cis regulatory differences is to analyze allele-specific expression (ASE) and identify cases of allelic imbalance (AI). RNA-seq is the most common way to measure ASE and a binomial test is often applied to determine statistical significance of AI. This implicitly assumes that there is no bias in estimation of AI. However, bias has been found to result from multiple factors including: genome ambiguity, reference quality, the mapping algorithm, and biases in the sequencing process. Two alternative approaches have been developed to handle bias: adjusting for bias using a statistical model and filtering regions of the genome suspected of harboring bias. Existing statistical models which account for bias rely on information from DNA controls, which can be cost prohibitive for large intraspecific studies. In contrast, data filtering is inexpensive and straightforward, but necessarily involves sacrificing a portion of the data.ResultsHere we propose a flexible Bayesian model for analysis of AI, which accounts for bias and can be implemented without DNA controls. In lieu of DNA controls, this Poisson-Gamma (PG) model uses an estimate of bias from simulations. The proposed model always has a lower type I error rate compared to the binomial test. Consistent with prior studies, bias dramatically affects the type I error rate. All of the tested models are sensitive to misspecification of bias. The closer the estimate of bias is to the true underlying bias, the lower the type I error rate. Correct estimates of bias result in a level alpha test.ConclusionsTo improve the assessment of AI, some forms of systematic error (e.g., map bias) can be identified using simulation. The resulting estimates of bias can be used to correct for bias in the PG model, without data filtering. Other sources of bias (e.g., unidentified variant calls) can be easily captured by DNA controls, but are missed by common filtering approaches. Consequently, as variant identification improves, the need for DNA controls will be reduced. Filtering does not significantly improve performance and is not recommended, as information is sacrificed without a measurable gain. The PG model developed here performs well when bias is known, or slightly misspecified. The model is flexible and can accommodate differences in experimental design and bias estimation.

Project description:Genetic dissection of the S rat genome has provided strong evidence for the presence of two interacting blood pressure (BP) quantitative trait loci (QTLs), termed QTL1 and QTL2, on rat chromosome 5. However, the identities of the underlying interacting genetic factors remain unknown. Further experiments targeted to identify the interacting genetic factors by the substitution mapping approach alone are difficult because of the interdependency of natural recombinations to occur at the two QTLs. We hypothesized that the interacting genetic factors underlying these two QTLs may interact at the level of gene transcription and thereby represent expression QTLs (eQTLs). To detect these interacting eQTLs, a custom QTL chip containing the annotated genes within QTL1 and QTL2 was developed and used to conduct a transcriptional profiling study of S and two congenic strains that retain either one or both the QTLs. The results uncovered an interaction between two transcription factors, DMRTA2 and NFIA. Further, the âbiological signatureâ elicited by these two transcription factors was differential between the congenic strain that retained LEW alleles at both QTL1 and 2 compared to the congenic strain that retained LEW alleles at QTL1 alone. A network of transcription factors potentially affecting BP could be traced, lending support to our hypothesis. Pairs of Cy5 and Cy3 labeled targets were co-hybridized onto either a custom long oligonucleotide microarray for the interrogation of 231 genes encompassed by QTL1 and QTL2, or a TIGR rat cDNA array consisting of 26,401 probe elements representing 20,465 unique non-QTL genes. A âflip-dyeâ or âbalanced blockâ design was used as the experimental method of choice to account for potential dye-bias labeling effects. Six âbalanced blockâ normalized files are submitted for the long oligonucleotide array interrogating the hearts from S versus S.LEW(5)x6x9 animals, fourteen âflip-dyeâ hybridizations are submitted for the cDNA array interrogating the hearts from S versus S.LEW(5)x6x9 animals, and twelve âflip-dyeâ hybridizations are submitted for the cDNA array interrogating the hearts from S versus S.LEW(5)x6x11 animals. GPR and MEV files cannot be located.

Dataset Information

Comparison of quantitative trait loci methods: Total expression and allelic imbalance method in brain RNA-seq.

Background

Methods

Results

Conclusion

Publications

Comparison of quantitative trait loci methods: Total expression and allelic imbalance method in brain RNA-seq.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets