Allelic mapping bias in RNA-sequencing is not a major confounder in eQTL studies

Background

RNA sequencing (RNA-seq) is the current gold-standard method to quantify gene expression for expression quantitative trait locus (eQTL) studies. However, a potential caveat in these studies is that RNA-seq reads carrying the non-reference allele of variant loci can have lower probability to map correctly to the reference genome, which could bias gene quantifications and cause false positive eQTL associations. In this study, we analyze the effect of this allelic mapping bias in eQTL discovery.

Results

We simulate RNA-seq read mapping over 9.5 M common SNPs and indels, with 15.6% of variants showing biased mapping rate for reference versus non-reference reads. However, removing potentially biased RNA-seq reads from an eQTL dataset of 185 individuals has a very small effect on gene and exon quantifications and eQTL discovery. We detect only a handful of likely false positive eQTLs, and overall eQTL SNPs show no significant enrichment for high mapping bias.

Conclusion

Our results suggest that RNA-seq quantifications are generally robust against allelic mapping bias, and that this does not have a severe effect on eQTL discovery. Nevertheless, we provide our catalog of putatively biased loci to allow better controlling for mapping bias to obtain more accurate results in future RNA-seq studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0467-2) contains supplementary material, which is available to authorized users.     

Experimental validation of methods for differential gene expression analysis and sample pooling in RNA-seq

Background

Massively parallel cDNA sequencing (RNA-seq) experiments are gradually superseding microarrays in quantitative gene expression profiling. However, many biologists are uncertain about the choice of differentially expressed gene (DEG) analysis methods and the validity of cost-saving sample pooling strategies for their RNA-seq experiments. Hence, we performed experimental validation of DEGs identified by Cuffdiff2, edgeR, DESeq2 and Two-stage Poisson Model (TSPM) in a RNA-seq experiment involving mice amygdalae micro-punches, using high-throughput qPCR on independent biological replicate samples. Moreover, we sequenced RNA-pools and compared their results with sequencing corresponding individual RNA samples.

Results

False-positivity rate of Cuffdiff2 and false-negativity rates of DESeq2 and TSPM were high. Among the four investigated DEG analysis methods, sensitivity and specificity of edgeR was relatively high. We documented the pooling bias and that the DEGs identified in pooled samples suffered low positive predictive values.

Conclusions

Our results highlighted the need for combined use of more sensitive DEG analysis methods and high-throughput validation of identified DEGs in future RNA-seq experiments. They indicated limited utility of sample pooling strategies for RNA-seq in similar setups and supported increasing the number of biological replicate samples.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1767-y) contains supplementary material, which is available to authorized users.     

DupChecker: a bioconductor package for checking high-throughput genomic data redundancy in meta-analysis

Background

Meta-analysis has become a popular approach for high-throughput genomic data analysis because it often can significantly increase power to detect biological signals or patterns in datasets. However, when using public-available databases for meta-analysis, duplication of samples is an often encountered problem, especially for gene expression data. Not removing duplicates could lead false positive finding, misleading clustering pattern or model over-fitting issue, etc in the subsequent data analysis.

Results

We developed a Bioconductor package Dupchecker that efficiently identifies duplicated samples by generating MD5 fingerprints for raw data. A real data example was demonstrated to show the usage and output of the package.

Conclusions

Researchers may not pay enough attention to checking and removing duplicated samples, and then data contamination could make the results or conclusions from meta-analysis questionable. We suggest applying DupChecker to examine all gene expression data sets before any data analysis step.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-323) contains supplementary material, which is available to authorized users.     

Analysis of Genes Expression of Spodoptera exigua Larvae upon AcMNPV Infection

Background

The impact of Autographa californica multicapsid nucleopolyhedrovirus (AcMNPV) infection on host gene expression in Spodoptera exigua 4th instar larvae was investigated through the use of 454 sequencing-based RNA-seq of cDNA libraries developed from insects challenged with active AcMNPV or heat-inactivated AcMNPV.

Methodology/Principal Findings

By comparing the two cDNA libraries, we show that 201 host genes are significantly up-regulated and 234 genes are significantly down-regulated by active AcMNPV infection. Down-regulated host genes included genes encoding antimicrobial peptides, namely three gloverin isoforms and an attacin, indicating that the viral infection actively repressed the expression of a portion of the host immune gene repertoire. Another interesting group of down-regulated host genes included genes encoding two juvenile hormone binding proteins and a hexamerin, all of which are involved in juvenile hormone regulation. The expression of these genes was enhanced by the topical application of Juvenile Hormone III (JHIII) in the insects challenged with heat-inactivated AcMNPV. However, infection with the active virus strongly suppresses the expression of these three genes, regardless of the absence or presence of JHIII.

Conclusions/Significance

Using RNA-seq, we have identified groups of immune-regulated and juvenile hormone-regulated genes that are suppressed by infection with active AcMNPV. This information and further studies on the regulation of host gene expression by AcMNPV will provide the tools needed to enhance the utility of the virus as an effective protein expression system and as an insecticide.     

The RNA-binding protein hnRNPLL induces a T cell alternative splicing program delineated by differential intron retention in polyadenylated RNA

Background

Retention of a subset of introns in spliced polyadenylated mRNA is emerging as a frequent, unexplained finding from RNA deep sequencing in mammalian cells.

Results

Here we analyze intron retention in T lymphocytes by deep sequencing polyadenylated RNA. We show a developmentally regulated RNA-binding protein, hnRNPLL, induces retention of specific introns by sequencing RNA from T cells with an inactivating Hnrpll mutation and from B lymphocytes that physiologically downregulate Hnrpll during their differentiation. In Ptprc mRNA encoding the tyrosine phosphatase CD45, hnRNPLL induces selective retention of introns flanking exons 4 to 6; these correspond to the cassette exons containing hnRNPLL binding sites that are skipped in cells with normal, but not mutant or low, hnRNPLL. We identify similar patterns of hnRNPLL-induced differential intron retention flanking alternative exons in 14 other genes, representing novel elements of the hnRNPLL-induced splicing program in T cells. Retroviral expression of a normally spliced cDNA for one of these targets, Senp2, partially corrects the survival defect of Hnrpll-mutant T cells. We find that integrating a number of computational methods to detect genes with differentially retained introns provides a strategy to enrich for alternatively spliced exons in mammalian RNA-seq data, when complemented by RNA-seq analysis of purified cells with experimentally perturbed RNA-binding proteins.

Conclusions

Our findings demonstrate that intron retention in mRNA is induced by specific RNA-binding proteins and suggest a biological significance for this process in marking exons that are poised for alternative splicing.     

A universal protocol to generate consensus level genome sequences for foot-and-mouth disease virus and other positive-sense polyadenylated RNA viruses using the Illumina MiSeq

Background

Next-Generation Sequencing (NGS) is revolutionizing molecular epidemiology by providing new approaches to undertake whole genome sequencing (WGS) in diagnostic settings for a variety of human and veterinary pathogens. Previous sequencing protocols have been subject to biases such as those encountered during PCR amplification and cell culture, or are restricted by the need for large quantities of starting material. We describe here a simple and robust methodology for the generation of whole genome sequences on the Illumina MiSeq. This protocol is specific for foot-and-mouth disease virus (FMDV) or other polyadenylated RNA viruses and circumvents both the use of PCR and the requirement for large amounts of initial template.

Results

The protocol was successfully validated using five FMDV positive clinical samples from the 2001 epidemic in the United Kingdom, as well as a panel of representative viruses from all seven serotypes. In addition, this protocol was successfully used to recover 94% of an FMDV genome that had previously been identified as cell culture negative. Genome sequences from three other non-FMDV polyadenylated RNA viruses (EMCV, ERAV, VESV) were also obtained with minor protocol amendments. We calculated that a minimum coverage depth of 22 reads was required to produce an accurate consensus sequence for FMDV O. This was achieved in 5 FMDV/O/UKG isolates and the type O FMDV from the serotype panel with the exception of the 5′ genomic termini and area immediately flanking the poly(C) region.

Conclusions

We have developed a universal WGS method for FMDV and other polyadenylated RNA viruses. This method works successfully from a limited quantity of starting material and eliminates the requirement for genome-specific PCR amplification. This protocol has the potential to generate consensus-level sequences within a routine high-throughput diagnostic environment.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-828) contains supplementary material, which is available to authorized users.     

Impact of ozone exposure on the response to glucocorticoid in a mouse model of asthma: involvements of p38 MAPK and MKP-1

Background

Molecular mechanisms involved in the oxidative stress induced glucocorticoids insensitivity remain elusive. The mitogen-activated protein kinase phosphatase (MKP) 1 mediates a part of glucocorticoids action and can be modified by exogenous oxidants. Whether oxidant ozone (O₃) can affect the function of MKP-1 and hence blunt the response to corticotherapy is not clear.

Methods

Here we employed a murine model of asthma established with ovalbumin (OVA) sensitization and challenge to evaluate the influence of O₃ on the inhibitory effect of dexamethasone on AHR and airway inflammation, and by administration of SB239063, a selective p38 MAPK inhibitor, to explore the underlying involvements of the activation of p38 MAPK and the expression of MKP-1.

Results

Ozone exposure not only aggravated the pulmonary inflammation and AHR, but also decreased the inhibitory effects of dexamethasone, accompanied by the elevated oxidative stress, airway neutrophilia, enhanced phosphorylation of p38 MAPK, and upregulated expression of IL-17. Administration of SB239063 caused significant inhibition of the p38 MAPK phosphorylation, alleviation of the airway neutrophilia, and decrement of the ozone-induced IL-17 expression, and partly restored the ozone-impaired effects of dexamethasone. Ozone exposure not only decreased the protein expression of MKP-1, but also diminished the dexamethasone-mediated induction process of MKP-1 mRNA and protein expression.

Conclusions

The glucocorticoids insensitivity elicited by ozone exposure on current asthma model may involve the enhanced phosphorylation of p38 MAPK and disturbed expression of MKP-1.

Electronic supplementary material

The online version of this article (doi:10.1186/s12931-014-0126-x) contains supplementary material, which is available to authorized users.