Synthesis of 53 tissue and cell line expression QTL datasets reveals master eQTLs

Background

Gene expression genetic studies in human tissues and cells identify cis- and trans-acting expression quantitative trait loci (eQTLs). These eQTLs provide insights into regulatory mechanisms underlying disease risk. However, few studies systematically characterized eQTL results across cell and tissues types. We synthesized eQTL results from >50 datasets, including new primary data from human brain, peripheral plaque and kidney samples, in order to discover features of human eQTLs.

Results

We find a substantial number of robust cis-eQTLs and far fewer trans-eQTLs consistent across tissues. Analysis of 45 full human GWAS scans indicates eQTLs are enriched overall, and above nSNPs, among positive statistical signals in genetic mapping studies, and account for a significant fraction of the strongest human trait effects. Expression QTLs are enriched for gene centricity, higher population allele frequencies, in housekeeping genes, and for coincidence with regulatory features, though there is little evidence of 5′ or 3′ positional bias. Several regulatory categories are not enriched including microRNAs and their predicted binding sites and long, intergenic non-coding RNAs. Among the most tissue-ubiquitous cis-eQTLs, there is enrichment for genes involved in xenobiotic metabolism and mitochondrial function, suggesting these eQTLs may have adaptive origins. Several strong eQTLs (CDK5RAP2, NBPFs) coincide with regions of reported human lineage selection. The intersection of new kidney and plaque eQTLs with related GWAS suggest possible gene prioritization. For example, butyrophilins are now linked to arterial pathogenesis via multiple genetic and expression studies. Expression QTL and GWAS results are made available as a community resource through the NHLBI GRASP database [http://apps.nhlbi.nih.gov/grasp/].

Conclusions

Expression QTLs inform the interpretation of human trait variability, and may account for a greater fraction of phenotypic variability than protein-coding variants. The synthesis of available tissue eQTL data highlights many strong cis-eQTLs that may have important biologic roles and could serve as positive controls in future studies. Our results indicate some strong tissue-ubiquitous eQTLs may have adaptive origins in humans. Efforts to expand the genetic, splicing and tissue coverage of known eQTLs will provide further insights into human gene regulation.

Electronic supplementary material

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

Segregation of regulatory polymorphisms with effects on the gluteus medius transcriptome in a purebred pig population

Background

The main goal of the present study was to analyse the genetic architecture of mRNA expression in muscle, a tissue with an outmost economic importance for pig breeders. Previous studies have used F₂ crosses to detect porcine expression QTL (eQTL), so they contributed with data that mostly represents the between-breed component of eQTL variation. Herewith, we have analysed eQTL segregation in an outbred Duroc population using two groups of animals with divergent fatness profiles. This approach is particularly suitable to analyse the within-breed component of eQTL variation, with a special emphasis on loci involved in lipid metabolism.

Methodology/Principal Findings

GeneChip Porcine Genome arrays (Affymetrix) were used to determine the mRNA expression levels of gluteus medius samples from 105 Duroc barrows. A whole-genome eQTL scan was carried out with a panel of 116 microsatellites. Results allowed us to detect 613 genome-wide significant eQTL unevenly distributed across the pig genome. A clear predominance of trans- over cis-eQTL, was observed. Moreover, 11 trans-regulatory hotspots affecting the expression levels of four to 16 genes were identified. A Gene Ontology study showed that regulatory polymorphisms affected the expression of muscle development and lipid metabolism genes. A number of positional concordances between eQTL and lipid trait QTL were also found, whereas limited evidence of a linear relationship between muscle fat deposition and mRNA levels of eQTL regulated genes was obtained.

Conclusions/Significance

Our data provide substantial evidence that there is a remarkable amount of within-breed genetic variation affecting muscle mRNA expression. Most of this variation acts in trans and influences biological processes related with muscle development, lipid deposition and energy balance. The identification of the underlying causal mutations and the ascertainment of their effects on phenotypes would allow gaining a fundamental perspective about how complex traits are built at the molecular level.     

Identification of microRNAs associated with allergic airway disease using a genetically diverse mouse population

Background

Allergic airway diseases (AADs) such as asthma are characterized in part by granulocytic airway inflammation. The gene regulatory networks that govern granulocyte recruitment are poorly understood, but evidence is accruing that microRNAs (miRNAs) play an important role. To identify miRNAs that may underlie AADs, we used two complementary approaches that leveraged the genotypic and phenotypic diversity of the Collaborative Cross (CC) mouse population. In the first approach, we sought to identify miRNA expression quantitative trait loci (eQTL) that overlap QTL for AAD-related phenotypes. Specifically, CC founder strains and incipient lines of the CC were sensitized and challenged with house dust mite allergen followed by measurement of granulocyte recruitment to the lung. Total lung RNA was isolated and miRNA was measured using arrays for CC founders and qRT-PCR for incipient CC lines.

Results

Among CC founders, 92 miRNAs were differentially expressed. We measured the expression of 40 of the most highly expressed of these 92 miRNAs in the incipient lines of the CC and identified 18 eQTL corresponding to 14 different miRNAs. Surprisingly, half of these eQTL were distal to the corresponding miRNAs, and even on different chromosomes. One of the largest-effect local miRNA eQTL was for miR-342-3p, for which we identified putative causal variants by bioinformatic analysis of the effects of single nucleotide polymorphisms on RNA structure. None of the miRNA eQTL co-localized with QTL for eosinophil or neutrophil recruitment. In the second approach, we constructed putative miRNA/mRNA regulatory networks and identified three miRNAs (miR-497, miR-351 and miR-31) as candidate master regulators of genes associated with neutrophil recruitment. Analysis of a dataset from human keratinocytes transfected with a miR-31 inhibitor revealed two target genes in common with miR-31 targets correlated with neutrophils, namely Oxsr1 and Nsf.

Conclusions

miRNA expression in the allergically inflamed murine lung is regulated by genetic loci that are smaller in effect size compared to mRNA eQTL and often act in trans. Thus our results indicate that the genetic architecture of miRNA expression is different from mRNA expression. We identified three miRNAs, miR-497, miR-351 and miR-31, that are candidate master regulators of genes associated with neutrophil recruitment. Because miR-31 is expressed in airway epithelia and is predicted to target genes with known links to neutrophilic inflammation, we suggest that miR-31 is a potentially novel regulator of airway inflammation.

Electronic supplementary material

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

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.     

Comparative genomic analysis of clinical and environmental strains provides insight into the pathogenicity and evolution of Vibrio parahaemolyticus

Background

Vibrio parahaemolyticus is a Gram-negative halophilic bacterium. Infections with the bacterium could become systemic and can be life-threatening to immunocompromised individuals. Genome sequences of a few clinical isolates of V. parahaemolyticus are currently available, but the genome dynamics across the species and virulence potential of environmental strains on a genome-scale have not been described before.

Results

Here we present genome sequences of four V. parahaemolyticus clinical strains from stool samples of patients and five environmental strains in Hong Kong. Phylogenomics analysis based on single nucleotide polymorphisms revealed a clear distinction between the clinical and environmental isolates. A new gene cluster belonging to the biofilm associated proteins of V. parahaemolyticus was found in clincial strains. In addition, a novel small genomic island frequently found among clinical isolates was reported. A few environmental strains were found harboring virulence genes and prophage elements, indicating their virulence potential. A unique biphenyl degradation pathway was also reported. A database for V. parahaemolyticus (http://kwanlab.bio.cuhk.edu.hk/vp) was constructed here as a platform to access and analyze genome sequences and annotations of the bacterium.

Conclusions

We have performed a comparative genomics analysis of clinical and environmental strains of V. parahaemolyticus. Our analyses could facilitate understanding of the phylogenetic diversity and niche adaptation of this bacterium.

Electronic supplementary material

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

Allele-specific expression and eQTL analysis in mouse adipose tissue

Background

The 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.

Results

In 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.

Conclusions

We 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.

Electronic supplementary material

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

Genetic Determinants of On-Aspirin Platelet Reactivity: Focus on the Influence of PEAR1

Background

Platelet aggregation during aspirin treatment displays considerable inter-individual variability. A genetic etiology likely exists, but it remains unclear to what extent genetic polymorphisms determine platelet aggregation in aspirin-treated individuals.

Aim

To identify platelet-related single nucleotide polymorphisms (SNPs) influencing platelet aggregation during aspirin treatment. Furthermore, we explored to what extent changes in cyclooxygenase-1 activity and platelet activation may explain such influence.

Methods

We included 985 Danish patients with stable coronary artery disease treated with aspirin 75 mg/day mono antiplatelet therapy. Patients were genotyped for 16 common SNPs in platelet-related genes using standard PCR-based methods (TaqMan). Platelet aggregation was evaluated by whole blood platelet aggregometry employing Multiplate Analyzer (agonists: arachidonic acid and collagen) and VerifyNow Aspirin. Serum thromboxane B₂ was measured to confirm aspirin adherence and was used as a marker of cyclooxygenase-1 activity. Soluble P-selectin was used as marker of platelet activation. Platelet aggregation, cyclooxygenase-1 activity, and platelet activation were compared across genotypes in adjusted analyses.

Results

The A-allele of the rs12041331 SNP in the platelet endothelial aggregation receptor-1 (PEAR1) gene was associated with reduced platelet aggregation and increased platelet activation, but not with cyclooxygenase-1 activity. Platelet aggregation was unaffected by the other SNPs analyzed.

Conclusion

A common genetic variant in PEAR1 (rs12041331) reproducibly influenced platelet aggregation in aspirin-treated patients with coronary artery disease. The exact biological mechanism remains elusive, but the effect of this polymorphism may be related to changes in platelet activation. Furthermore, 14 SNPs previously suggested to influence aspirin efficacy were not associated with on-aspirin platelet aggregation.

Clinical Trial Registration

ClinicalTrials.gov {"type":"clinical-trial","attrs":{"text":"NCT01383304","term_id":"NCT01383304"}}NCT01383304     

Using the One-Lung Method to Link p38 to Pro-Inflammatory Gene Expression during Overventilation in C57BL/6 and BALB/c Mice

Introduction

The mechanisms of ventilator-induced lung injury (VILI), including the role of MAP kinases, are frequently studied in different mouse strains. A useful model for such studies is the isolated perfused mouse lung. As a further development we present the one-lung method that permits to continue perfusion and ventilation of the right lung after removal of the left lung. This method was used to compare the effect of high pressure ventilation (HPV) on pro-inflammatory signaling events in two widely used mouse strains (C57BL/6, BALB/c) and to further define the role of p38 in VILI.

Methods

Lungs were perfused and ventilated for 30 min under control conditions before they were randomized to low (8 cm H₂O) or high (25 cm H₂O) pressure ventilation (HPV) for 210 min, with the left lung being removed after 180 min. In the left lung we measured the phosphorylation of p38, JNK, ERK and Akt kinase, and in the right lung gene expression and protein concentrations of Il1b, Il6, Tnf, Cxcl1, Cxcl2, and Areg.

Results

Lung mechanics and kinase activation were similar in both mouse strains. HPV increased all genes (except Tnf in BALB/c) and all mediators in both strains. The gene expression of mRNA for Il1b, Il6, Cxcl1 and Cxcl2 was higher in BALB/c mice. Backward regression of the kinase data at t = 180 min with the gene and protein expression data at t = 240 min suggested that p38 controls HPV-induced gene expression, but not protein production. This hypothesis was confirmed in experiments with the p38-kinase inhibitor SB203580.

Conclusions

The one-lung method is useful for mechanistic studies in the lungs. While C57BL/6 show diminished pro-inflammatory responses during HPV, lung mechanics and mechanotransduction processes appear to be similar in both mouse strains. Finally, the one-lung method allowed us to link p38 to gene expression during VILI.     

MEXPRESS: visualizing expression,DNA methylation and clinical TCGA data

Background

In recent years, increasing amounts of genomic and clinical cancer data have become publically available through large-scale collaborative projects such as The Cancer Genome Atlas (TCGA). However, as long as these datasets are difficult to access and interpret, they are essentially useless for a major part of the research community and their scientific potential will not be fully realized. To address these issues we developed MEXPRESS, a straightforward and easy-to-use web tool for the integration and visualization of the expression, DNA methylation and clinical TCGA data on a single-gene level (http://mexpress.be).

Results

In comparison to existing tools, MEXPRESS allows researchers to quickly visualize and interpret the different TCGA datasets and their relationships for a single gene, as demonstrated for GSTP1 in prostate adenocarcinoma. We also used MEXPRESS to reveal the differences in the DNA methylation status of the PAM50 marker gene MLPH between the breast cancer subtypes and how these differences were linked to the expression of MPLH.

Conclusions

We have created a user-friendly tool for the visualization and interpretation of TCGA data, offering clinical researchers a simple way to evaluate the TCGA data for their genes or candidate biomarkers of interest.

Electronic supplementary material

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

Phase I clinical trial of systemically administered TUSC2(FUS1)-nanoparticles mediating functional gene transfer in humans

Background

Tumor suppressor gene TUSC2/FUS1 (TUSC2) is frequently inactivated early in lung cancer development. TUSC2 mediates apoptosis in cancer cells but not normal cells by upregulation of the intrinsic apoptotic pathway. No drug strategies currently exist targeting loss-of–function genetic abnormalities. We report the first in-human systemic gene therapy clinical trial of tumor suppressor gene TUSC2.

Methods

Patients with recurrent and/or metastatic lung cancer previously treated with platinum-based chemotherapy were treated with escalating doses of intravenous N-[1-(2,3-dioleoyloxy)propyl]-N,N,N-trimethylammonium chloride (DOTAP):cholesterol nanoparticles encapsulating a TUSC2 expression plasmid (DOTAP:chol-TUSC2) every 3 weeks.

Results

Thirty-one patients were treated at 6 dose levels (range 0.01 to 0.09 milligrams per kilogram). The MTD was determined to be 0.06 mg/kg. Five patients achieved stable disease (2.6–10.8 months, including 2 minor responses). One patient had a metabolic response on positron emission tomography (PET) imaging. RT-PCR analysis detected TUSC2 plasmid expression in 7 of 8 post-treatment tumor specimens but not in pretreatment specimens and peripheral blood lymphocyte controls. Proximity ligation assay, performed on paired biopsies from 3 patients, demonstrated low background TUSC2 protein staining in pretreatment tissues compared with intense (10–25 fold increase) TUSC2 protein staining in post-treatment tissues. RT-PCR gene expression profiling analysis of apoptotic pathway genes in two patients with high post-treatment levels of TUSC2 mRNA and protein showed significant post-treatment changes in the intrinsic apoptotic pathway. Twenty-nine genes of the 82 tested in the apoptosis array were identified by Igenuity Pathway Analysis to be significantly altered post-treatment in both patients (Pearson correlation coefficient 0.519; p<0.01).

Conclusions

DOTAP:chol-TUSC2 can be safely administered intravenously in lung cancer patients and results in uptake of the gene by human primary and metastatic tumors, transgene and gene product expression, specific alterations in TUSC2-regulated pathways, and anti-tumor effects (to our knowledge for the first time for systemic DOTAP:cholesterol nanoparticle gene therapy).

Trial Registration

ClinicalTrials.gov {"type":"clinical-trial","attrs":{"text":"NCT00059605","term_id":"NCT00059605"}}NCT00059605