Dataset Information

Identification of genetic variants associated with alternative splicing using sQTLseekeR.

ABSTRACT: Identification of genetic variants affecting splicing in RNA sequencing population studies is still in its infancy. Splicing phenotype is more complex than gene expression and ought to be treated as a multivariate phenotype to be recapitulated completely. Here we represent the splicing pattern of a gene as the distribution of the relative abundances of a gene's alternative transcript isoforms. We develop a statistical framework that uses a distance-based approach to compute the variability of splicing ratios across observations, and a non-parametric analogue to multivariate analysis of variance. We implement this approach in the R package sQTLseekeR and use it to analyze RNA-Seq data from the Geuvadis project in 465 individuals. We identify hundreds of single nucleotide polymorphisms (SNPs) as splicing QTLs (sQTLs), including some falling in genome-wide association study SNPs. By developing the appropriate metrics, we show that sQTLseekeR compares favorably with existing methods that rely on univariate approaches, predicting variants that behave as expected from mutations affecting splicing.

SUBMITTER: Monlong J

PROVIDER: S-EPMC4143934 | biostudies-literature | 2014 Aug

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Identification of genetic variants associated with alternative splicing using sQTLseekeR.

Monlong Jean J Calvo Miquel M Ferreira Pedro G PG Guigó Roderic R

Nature communications 20140820

Identification of genetic variants affecting splicing in RNA sequencing population studies is still in its infancy. Splicing phenotype is more complex than gene expression and ought to be treated as a multivariate phenotype to be recapitulated completely. Here we represent the splicing pattern of a gene as the distribution of the relative abundances of a gene's alternative transcript isoforms. We develop a statistical framework that uses a distance-based approach to compute the variability of sp ...[more]

PMID: 25140736

Similar Datasets

Project description:ObjectiveAccumulating evidence has shown that aberrant alternative splicing events are closely associated with the onset and development of cancer. However, whether genetic variants-associated alternative splicing is linked to risk of endometrial cancer remains largely uncertain.MethodsWe identified single nucleotide polymorphisms (SNPs) locates in the splicing number trait locus (sQTL) of endometrial cancer using the CancerSplicing QTL database. In parallel with bioinformatics analysis, we conducted a case-control study comprising 2,000 cases and 2,013 controls to assess the association between identified SNP which possesses mRNA splicing function and endometrial cancer susceptibility. Furthermore, we used the Kaplan-Meier Plotter, The Human Protein Atlas, SPNR, and Spliceman2 databases for sQTL and differential gene expression analyses to identify the genetic variant which most potentially influence the risk of endometrial cancer through alternative splicing to reveal the potential mechanism by which candidate SNPs regulate the risk of endometrial cancer.ResultsThe results indicated that SNP rs7128029 A<G was significantly associated with an increased risk of endometrial cancer (odds ratio=1.384; 95% confidence interval=1.038-1.964). Moreover, the carcinogenic effect of SNP rs7128029 A<G was consistently revealed by propensity matching analysis, an additive model, and a dominant model. Importantly, sQTL analysis showed that SNP rs7128029 could affect the transcriptional modification of PSMD13 via regulating the exon skipping of PSMD13. When the rs7128029 allele was mutated from A to G, the expression of exon 2 of PSMD13 was markedly lower (p<0.001). Furthermore, compared with participants who had higher PSMD13 expression, those who had lower PSMD13 expression had shorter survival times.ConclusionThese findings suggest that SNP rs7128029-mediated alternative splicing events in PSMD13 are associated with endometrial cancer risk and may be a potential early screening biomarker for endometrial cancer-susceptible populations.

Project description:Paratuberculosis is a major endemic disease caused by Mycobacterium avium subspecies paratuberculosis (MAP) infection and leads to huge economic loss in the dairy sector worldwide. Alternative splicing (AS) events, playing indispensable regulatory roles in many protein functions and biological pathways, are shown to be associated with complex traits and diseases. In this study, by integrating the RNA sequencing (RNA-seq) data of 24 samples from three tissues (peripheral blood, jejunum and salivary gland) of Holstein cows, we obtained 2,706,541,696 uniquely mapped reads in total that represented 12,870 expressed genes, and detected 4285 differentially expressed genes (DEGs) between MAP-infected and healthy cows (p < 0.05). Of them, 92 differentially expressed splicing factors (DESFs) were included. Further, 119, 150 and 68 differential alternative splicing (DAS) events between MAP-infected and healthy cows were identified in peripheral blood, jejunum and salivary glands, respectively. Of note, six DAS events were highly and significantly correlated with the DESFs (R2 > 0.9; p < 0.01), and their corresponding genes (COPI coat complex subunit gamma 2gene (COPG2), kinesin family member 2C gene (KIF2C), exocyst complex component 7 (EXOC7), Rab9 effector protein with kelch motifs gene (RABEPK), deoxyribonuclease 1 gene (DNASE1) and early endosome antigen 1gene (EEA1)) were significantly enriched in immune response such as vesicle-mediated transport, regulation of acute inflammatory response and tuberculosis through gene ontology (GO) and KEGG analysis. KS test showed that the DAS events in the EXOC7 and KIF2C genes indeed displayed differences between MAP-infected cows and healthy cows. The DAS in EXOC7 might produce a new protein sequence with lack of 23 amino acids, and the DAS in KIF2C induced a stop codon of premature occurrence and resulted in a lack of functional domain. In summary, this study identified the DAS events and corresponding genes related to MAP-infection base on the RNA-seq data from multiple tissues of Holstein cows, providing novel insights into the regulatory mechanisms underpinning paratuberculosis in dairy cattle.

Dataset Information

Identification of genetic variants associated with alternative splicing using sQTLseekeR.

Publications

Identification of genetic variants associated with alternative splicing using sQTLseekeR.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets