Analysis of allele-specific expression in mouse liver by RNA-Seq: a comparison with Cis-eQTL identified using genetic linkage.
Ontology highlight
ABSTRACT:
SUBMITTER: Lagarrigue S
PROVIDER: S-EPMC3813844 | biostudies-literature | 2013 Nov
REPOSITORIES: biostudies-literature
Publications
Analysis of allele-specific expression in mouse liver by RNA-Seq: a comparison with Cis-eQTL identified using genetic linkage.
Lagarrigue Sandrine S Martin Lisa L Hormozdiari Farhad F Roux Pierre-François PF Pan Calvin C van Nas Atila A Demeure Olivier O Cantor Rita R Ghazalpour Anatole A Eskin Eleazar E Lusis Aldons J AJ
Genetics 20130911 3
We report an analysis of allele-specific expression (ASE) and parent-of-origin expression in adult mouse liver using next generation sequencing (RNA-Seq) of reciprocal crosses of heterozygous F1 mice from the parental strains C57BL/6J and DBA/2J. We found a 60% overlap between genes exhibiting ASE and putative cis-acting expression quantitative trait loci (cis-eQTL) identified in an intercross between the same strains. We discuss the various biological and technical factors that contribute to th ...[more]
Similar Datasets
Project description:Studies of expression quantitative trait loci (eQTLs) offer insight into the molecular mechanisms of loci that were found to be associated with complex diseases and the mechanisms can be classified into cis- and trans-acting regulation. At present, high-throughput RNA sequencing (RNA-seq) is rapidly replacing expression microarrays to assess gene expression abundance. Unlike microarrays that only measure the total expression of each gene, RNA-seq also provides information on allele-specific expression (ASE), which can be used to distinguish cis-eQTLs from trans-eQTLs and, more importantly, enhance cis-eQTL mapping. However, assessing the cis-effect of a candidate eQTL on a gene requires knowledge of the haplotypes connecting the candidate eQTL and the gene, which cannot be inferred with certainty. The existing two-stage approach that first phases the candidate eQTL against the gene and then treats the inferred phase as observed in the association analysis tends to attenuate the estimated cis-effect and reduce the power for detecting a cis-eQTL. In this article, we provide a maximum-likelihood framework for cis-eQTL mapping with RNA-seq data. Our approach integrates the inference of haplotypes and the association analysis into a single stage, and is thus unbiased and statistically powerful. We also develop a pipeline for performing a comprehensive scan of all local eQTLs for all genes in the genome by controlling for false discovery rate, and implement the methods in a computationally efficient software program. The advantages of the proposed methods over the existing ones are demonstrated through realistic simulation studies and an application to empirical breast cancer data from The Cancer Genome Atlas project.
| S-EPMC4642818 | biostudies-literature
Project description:As RNA-seq is replacing gene expression microarrays to assess genome-wide transcription abundance, gene expression Quantitative Trait Locus (eQTL) studies using RNA-seq have emerged. RNA-seq delivers two novel features that are important for eQTL studies. First, it provides information on allele-specific expression (ASE), which is not available from gene expression microarrays. Second, it generates unprecedentedly rich data to study RNA-isoform expression. In this paper, we review current methods for eQTL mapping using ASE and discuss some future directions. We also review existing works that use RNA-seq data to study RNA-isoform expression and we discuss the gaps between these works and isoform-specific eQTL mapping.
| S-EPMC3650863 | biostudies-literature
Project description:RNA-Sequencing (RNA-Seq) experiments have been optimized for library preparation, mapping, and gene expression estimation. These methods, however, have revealed weaknesses in the next stages of analysis of differential expression, with results sensitive to systematic sample stratification or, in more extreme cases, to outliers. Further, a method to assess normalization and adjustment measures imposed on the data is lacking.To address these issues, we utilize previously published eQTLs as a novel gold standard at the center of a framework that integrates DNA genotypes and RNA-Seq data to optimize analysis and aid in the understanding of genetic variation and gene expression. After detecting sample contamination and sequencing outliers in RNA-Seq data, a set of previously published brain eQTLs was used to determine if sample outlier removal was appropriate. Improved replication of known eQTLs supported removal of these samples in downstream analyses. eQTL replication was further employed to assess normalization methods, covariate inclusion, and gene annotation. This method was validated in an independent RNA-Seq blood data set from the GTEx project and a tissue-appropriate set of eQTLs. eQTL replication in both data sets highlights the necessity of accounting for unknown covariates in RNA-Seq data analysis.As each RNA-Seq experiment is unique with its own experiment-specific limitations, we offer an easily-implementable method that uses the replication of known eQTLs to guide each step in one's data analysis pipeline. In the two data sets presented herein, we highlight not only the necessity of careful outlier detection but also the need to account for unknown covariates in RNA-Seq experiments.
| S-EPMC3890578 | biostudies-literature
Project description:BackgroundThe simplest definition of cis-eQTLs versus trans, refers to genetic variants that affect expression in an allele specific manner, with implications on underlying mechanism. Yet, due to technical limitations of expression microarrays, the vast majority of eQTL studies performed in the last decade used a genomic distance based definition as a surrogate for cis, therefore exploring local rather than cis-eQTLs.ResultsIn this study we use RNAseq to explore allele specific expression (ASE) in adipose tissue of male and female F1 mice, produced from reciprocal crosses of C57BL/6J and DBA/2J strains. Comparison of the identified cis-eQTLs, to local-eQTLs, that were obtained from adipose tissue expression in two previous population based studies in our laboratory, yields poor overlap between the two mapping approaches, while both local-eQTL studies show highly concordant results. Specifically, local-eQTL studies show ~60% overlap between themselves, while only 15-20% of local-eQTLs are identified as cis by ASE, and less than 50% of ASE genes are recovered in local-eQTL studies. Utilizing recently published ENCODE data, we also find that ASE genes show significant bias for SNPs prevalence in DNase I hypersensitive sites that is ASE direction specific.ConclusionsWe suggest a new approach to analysis of allele specific expression that is more sensitive and accurate than the commonly used fisher or chi-square statistics. Our analysis indicates that technical differences between the cis and local-eQTL approaches, such as differences in genomic background or sex specificity, account for relatively small fraction of the discrepancy. Therefore, we suggest that the differences between two eQTL mapping approaches may facilitate sorting of SNP-eQTL interactions into true cis and trans, and that a considerable portion of local-eQTL may actually represent trans interactions.
| S-EPMC4089026 | biostudies-literature
Project description:RNA-seq may replace gene expression microarrays in the near future. Using RNA-seq, the expression of a gene can be estimated using the total number of sequence reads mapped to that gene, known as the total read count (TReC). Traditional expression quantitative trait locus (eQTL) mapping methods, such as linear regression, can be applied to TReC measurements after they are properly normalized. In this article, we show that eQTL mapping, by directly modeling TReC using discrete distributions, has higher statistical power than the two-step approach: data normalization followed by linear regression. In addition, RNA-seq provides information on allele-specific expression (ASE) that is not available from microarrays. By combining the information from TReC and ASE, we can computationally distinguish cis- and trans-eQTL and further improve the power of cis-eQTL mapping. Both simulation and real data studies confirm the improved power of our new methods. We also discuss the design issues of RNA-seq experiments. Specifically, we show that by combining TReC and ASE measurements, it is possible to minimize cost and retain the statistical power of cis-eQTL mapping by reducing sample size while increasing the number of sequence reads per sample. In addition to RNA-seq data, our method can also be employed to study the genetic basis of other types of sequencing data, such as chromatin immunoprecipitation followed by DNA sequencing data. In this article, we focus on eQTL mapping of a single gene using the association-based method. However, our method establishes a statistical framework for future developments of eQTL mapping methods using RNA-seq data (e.g., linkage-based eQTL mapping), and the joint study of multiple genetic markers and/or multiple genes.
| S-EPMC3218220 | biostudies-literature
Project description:BACKGROUND:Genomic imprinting, resulting in parent-of-origin specific gene expression, plays a critical role in mammalian development. Here, we apply allele-specific RNA-seq on isogenic B6D2F1 mice to assay imprinted genes in tissues from early embryonic tissues between E3.5 and E7.25 and in pluripotent cell lines to evaluate maintenance of imprinted gene expression. For the cell lines, we include embryonic stem cells (ESCs) and epiblast stem cells (EpiSCs) derived from fertilized embryos and from embryos obtained after nuclear transfer (NT) or parthenogenetic activation (PGA). RESULTS:As homozygous genomic regions of PGA-derived cells are not compatible with allele-specific RNA-seq, we developed an RNA-seq-based genotyping strategy allowing identification of informative heterozygous regions. Global analysis shows that proper imprinted gene expression as observed in embryonic tissues is largely lost in the ESC lines included in this study, which mainly consisted of female ESCs. Differentiation of ESC lines to embryoid bodies or NPCs does not restore monoallelic expression of imprinted genes, neither did reprogramming of the serum-cultured ESCs to the pluripotent ground state by the use of 2 kinase inhibitors. Fertilized EpiSC and EpiSC-NT lines largely maintain imprinted gene expression, as did EpiSC-PGA lines that show known paternally expressed genes being silent and known maternally expressed genes consistently showing doubled expression. Notably, two EpiSC-NT lines show aberrant silencing of Rian and Meg3, two critically imprinted genes in mouse iPSCs. With respect to female EpiSC, most of the lines displayed completely skewed X inactivation suggesting a (near) clonal origin. CONCLUSIONS:Altogether, our analysis provides a comprehensive overview of imprinted gene expression in pluripotency and provides a benchmark to allow identification of cell lines that faithfully maintain imprinted gene expression and therefore retain full developmental potential.
| S-EPMC6376749 | biostudies-literature
Project description:The differentiation and maturation trajectories of fetal liver stem/progenitor cells (LSPCs) are not fully understood at single-cell resolution, and a priori knowledge of limited biomarkers could restrict trajectory tracking.We employed marker-free single-cell RNA-Seq to characterize comprehensive transcriptional profiles of 507 cells randomly selected from seven stages between embryonic day 11.5 and postnatal day 2.5 during mouse liver development, and also 52 Epcam-positive cholangiocytes from postnatal day 3.25 mouse livers. LSPCs in developing mouse livers were identified via marker-free transcriptomic profiling. Single-cell resolution dynamic developmental trajectories of LSPCs exhibited contiguous but discrete genetic control through transcription factors and signaling pathways. The gene expression profiles of cholangiocytes were more close to that of embryonic day 11.5 rather than other later staged LSPCs, cuing the fate decision stage of LSPCs. Our marker-free approach also allows systematic assessment and prediction of isolation biomarkers for LSPCs.Our data provide not only a valuable resource but also novel insights into the fate decision and transcriptional control of self-renewal, differentiation and maturation of LSPCs.
| S-EPMC5715535 | biostudies-literature
Project description:Inhibitory GABAergic interneurons create different brain activity patterns that correlate with behavioural states. In this characterizing study, we used single-cell RNA-Seq to analyse anatomically- and electrophysiologically identified hippocampal oriens-lacunosum moleculare (OLM) interneurons. OLMs express somatostatin (Sst), generate feedback inhibition and play important roles in theta oscillations and fear encoding. Although an anatomically- and biophysically homogenous population, OLMs presumably comprise of two functionally distinct types with different developmental origins, inferred from the expression pattern of serotonin type-3a (5-HT3a, or Htr3a) receptor subunit and 5-HT excitability in a set of OLMs. To broadly characterize OLM cells, we used the Sst-Cre and the BAC transgenic Htr3a-Cre mouse lines and separately analysed SstCre-OLM and Htr3aCre-OLM types. We found a surprisingly consistent expression of Npy in OLMs, which was previously not associated with the identity of this type. Our analyses furthermore revealed uniform expression of developmental origin-related genes, including transcription factors and neurexin isoforms, without providing support for the current view that OLMs may originate from multiple neurogenic zones. Together, we found that OLMs constitute a highly homogenous transcriptomic population. Finally, our results revealed surprisingly infrequent expression of Htr3a in only ~10% of OLMs and an apparently specific expression of the 5-HT3b subunit-coding gene Htr3b in Htr3aCre-OLMs, but not in SstCre-OLMs. However, additional in situ hybridization experiments suggested that the differential expression of Htr3b may represent an unexpected consequence arising from the design of the Htr3a-Cre BAC transgenic line.
| S-EPMC6973274 | biostudies-literature
Project description:BACKGROUND:The chitinase-like protein YKL-40 plays a major role in inhibiting the inflammasome. Deregulation of inflammasome activation is emerging as a key modulator of pathologic airway inflammation in patients with asthma. We determined whether cis-expression quantitative trait loci (eQTLs) of the gene that encodes YKL-40, chitinase 3-like 1 (CHI3L1), are involved in the onset of asthma or in specific asthma phenotypes. METHODS:This case-control study, which was conducted at the University of Tsukuba, Japan, included a total of 2709 adults from the Tsukuba genome-wide association study (GWAS) cohort (734 healthy volunteers and 237 asthma patients), the Tsukuba replication cohort (375 healthy adult volunteers and 381 adult asthma patients), and the Hokkaido replication cohort (554 healthy adult volunteers and 428 adult asthma patients). Among 34 cis-eQTLs in CHI3L1 in the lung, rs946261 was associated with adult asthma in these Japanese cohorts. The genetic impact of rs946261 on asthma was also examined according to the age at onset and adult asthma clusters. RESULTS:In the Tsukuba GWAS cohort, the C allele at rs946261 was significantly associated with reduced expression of CHI3L1 mRNA in the lung and with development of asthma (odds ratio (OR) 1.27; P?=?0.036). The association was also observed following analysis of the three Japanese cohorts (OR 1.16; P?=?0.013). A stronger association was found with late-onset asthma that developed at 41?years of age or later (OR 1.24; 95% confidence interval (CI) 1.07-1.45; P?=?0.0058) and with a specific asthma phenotype characterized by late onset, less atopy, and mild airflow obstruction (OR 1.29; 95% CI 1.03-1.61; P?=?0.027). CONCLUSIONS:The genotype consisting of the cis-eQTL allele that reduces expression of CHI3L1 was specifically associated with late-onset adult asthma. Given the important role of YKL-40 in many pathophysiological processes, including cell growth, migration, chemotaxis, reorganization, and tissue remodeling, it may be involved in an important pathogenic role in the establishment of inflammation and remodeling in asthmatic airways. Our findings may indicate the presence of a specific endotype related to exaggerated activation of YKL-40 in the pathogenesis of late-onset adult asthma.
| S-EPMC6444873 | biostudies-literature
Project description:Expression Quantitative Trait Loci (eQTL) analysis enables characterisation of functional genetic variation influencing expression levels of individual genes. In outbread populations, including humans, eQTLs are commonly analysed using the conventional linear model, adjusting for relevant covariates, assuming an allelic dosage model and a Gaussian error term. However, gene expression data generally have noise that induces heavy-tailed errors relative to the Gaussian distribution and often include atypical observations, or outliers. Such departures from modelling assumptions can lead to an increased rate of type II errors (false negatives), and to some extent also type I errors (false positives). Careful model checking can reduce the risk of type-I errors but often not type II errors, since it is generally too time-consuming to carefully check all models with a non-significant effect in large-scale and genome-wide studies. Here we propose the application of a robust linear model for eQTL analysis to reduce adverse effects of deviations from the assumption of Gaussian residuals. We present results from a simulation study as well as results from the analysis of real eQTL data sets. Our findings suggest that in many situations robust models have the potential to provide more reliable eQTL results compared to conventional linear models, particularly in respect to reducing type II errors due to non-Gaussian noise. Post-genomic data, such as that generated in genome-wide eQTL studies, are often noisy and frequently contain atypical observations. Robust statistical models have the potential to provide more reliable results and increased statistical power under non-Gaussian conditions. The results presented here suggest that robust models should be considered routinely alongside other commonly used methodologies for eQTL analysis.
| S-EPMC4436354 | biostudies-literature