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.     

A mutant Pfu DNA polymerase designed for advanced uracil-excision DNA engineering

Background  

The combined use of restriction enzymes with PCR has revolutionized molecular cloning, but is inherently restricted by the content of the manipulated DNA sequences. Uracil-excision based cloning is ligase and sequence independent and allows seamless fusion of multiple DNA sequences in simple one-tube reactions, with higher accuracy than overlapping PCR.     

Nonea pisidica (Boraginaceae-Boragineae), a new species from southwest Anatolia and its relationships inferred from karyology and cpDNA sequences

The new species Nonea pisidica (Boraginaceae-Boragineae) is described based on two collections by the authors from the region of Lake Burdur in southwest Anatolia, Turkey. It shows a combination of morphological characteristics concerning habit, flower, and fruit that distinguishes it from N. caspica and N. pallens, the two more similar taxa in Nonea sect. Nonea. Karyological analyses corroborate the separation of the three species, which have different chromosome complements and even base numbers. N. pisidica is characterised by 2n = 30, a complement previously unknown for west Asiatic species of Nonea. The dibasic haploid number x = 15 may have originated through amphidiploid hybridisation between two annual, diploid taxa with x = 7 and x = 8, such as N. lutea and N. caspica. The relationships of the new species were further analysed using trnL_UAA sequences of the chloroplast genome, which were obtained for 16 selected species of Nonea. The resulting phylogeny confirms that it is related to N. caspica and N. lutea, but not to N. pallens, in spite of morphological resemblance. Lack of relationship with the south Mediterranean N. vesicaria, the only other species of Nonea known to have 2n = 30, suggests that amphidiploidy may have played an important role in speciation processes through recurrent and polytopical occurrence.     

Genomic selection for fruit quality traits in apple (Malus×domestica Borkh.)

The genome sequence of apple (Malus×domestica Borkh.) was published more than a year ago, which helped develop an 8K SNP chip to assist in implementing genomic selection (GS). In apple breeding programmes, GS can be used to obtain genomic breeding values (GEBV) for choosing next-generation parents or selections for further testing as potential commercial cultivars at a very early stage. Thus GS has the potential to accelerate breeding efficiency significantly because of decreased generation interval or increased selection intensity. We evaluated the accuracy of GS in a population of 1120 seedlings generated from a factorial mating design of four females and two male parents. All seedlings were genotyped using an Illumina Infinium chip comprising 8,000 single nucleotide polymorphisms (SNPs), and were phenotyped for various fruit quality traits. Random-regression best liner unbiased prediction (RR-BLUP) and the Bayesian LASSO method were used to obtain GEBV, and compared using a cross-validation approach for their accuracy to predict unobserved BLUP-BV. Accuracies were very similar for both methods, varying from 0.70 to 0.90 for various fruit quality traits. The selection response per unit time using GS compared with the traditional BLUP-based selection were very high (>100%) especially for low-heritability traits. Genome-wide average estimated linkage disequilibrium (LD) between adjacent SNPs was 0.32, with a relatively slow decay of LD in the long range (r(2)?=?0.33 and 0.19 at 100 kb and 1,000 kb respectively), contributing to the higher accuracy of GS. Distribution of estimated SNP effects revealed involvement of large effect genes with likely pleiotropic effects. These results demonstrated that genomic selection is a credible alternative to conventional selection for fruit quality traits.     

Engagement of ubiquitination and de-ubiquitination at rostral ventrolateral medulla in experimental brain death

Background

Whereas brain death is a vitally important clinical phenomenon, our contemporary understanding on its underlying cellular mechanisms remains elusive. This study evaluated whether the ubiquitin-proteasome system (UPS) in the rostral ventrolateral medulla (RVLM), a neural substrate that our laboratory identified previously to be intimately related to brain death, is engaged in this fatal process.

Methods

We performed proteomics, Western Blot, real-time PCR, ELISA and pharmacological experiments in conjunction with a clinically relevant experimental endotoxemia model of brain death based on intravenous administration of Escherichia coli lipopolysaccharide in adult male Sprague–Dawley rats.

Results

Proteomics, Western blot and enzyme activity analyses demonstrated that polyubiquitination was preserved and de-ubiquitination by ubiquitin C-terminal hydrolase isozyme-L1 (UCH-L1) was sustained, alongside increased monoubiquitin availability or proteasome activity in RVLM over the course of experimental endotoxemia. However, real-time PCR revealed no significant alteration in proteasome subunit alpha type-1, ubiquitin or UCH-L1 at mRNA level. Functionally, whereas microinjection into the bilateral RVLM of proteasome inhibitors (lactacystin or proteasome inhibitor II) potentiated survival, an inhibitor of ubiquitin-recycling (ubiquitin aldehyde) or an UCH-L1 inhibitor exacerbated mortality.

Conclusions

We proposed previously that the progression towards brain death entails a tug-of-war between pro-death and pro-life programs in RVLM. It is conceivable that ubiquitination or de-ubiquitination in RVLM participate in brain death by regulating the degradation of the proteins involved in those programs.     

Ectopic Wnt signal determines the eyeless phenotype of zebrafish masterblind mutant

masterblind (mbl) is a zebrafish mutation characterised by the absence or reduction in size of the telencephalon, optic vesicles and olfactory placodes. We show that inhibition of Gsk3beta in zebrafish embryos either by overexpression of dominant negative dn gsk3beta mRNA or by lithium treatment after the midblastula transition phenocopies mbl. The loss of anterior neural tissue in mbl and lithium-treated embryos is preceded by posteriorization of presumptive anterior neuroectoderm during gastrulation, which is evident from the anterior shift of marker genes Otx2 and Wnt1. Heterozygous mbl embryos showed increased sensitivity to inhibition of GSK3beta by lithium or dn Xgsk3beta that led to the loss of eyes. Overexpression of gsk3beta mRNA rescued eyes and the wild-type fgf8 expression of homozygous mbl embryos. emx1 that delineates the telencephalon is expanded and shifted ventroanteriorly in mbl embryos. In contrast to fgf8, the emx1 expression domain was not restored upon overexpression of gsk3beta mRNA. These experiments place mbl as an antagonist of the Wnt pathway in parallel or upstream of the complex consisting of Axin, APC and Gsk3beta that binds and phosphorylates beta-catenin, thereby destabilising it. mbl maps on LG 3 close to a candidate gene axin1. In mbl we detected a point mutation in the conserved minimal Gsk3beta-binding domain of axin1 leading to a leucine to glutamine substitution at position 399. Overexpression of wild-type axin1 mRNA rescued mbl completely, demonstrating that mutant axin1 is responsible for the mutant phenotype. Overexpression of mutant L399Q axin1 in wild-type embryos resulted in a dose-dependent dominant negative activity as demonstrated by the loss of telencephalon and eyes. We suggest that the function of Axin1/Mbl protein is to antagonise the Wnt signal and in doing so to establish and maintain the most anterior CNS. Our findings provide new insights into the mechanisms by which the Wnt pathway generates anteroposterior polarity of the neural plate.     

Validation of Heart Failure Events in the Antihypertensive and Lipid Lowering Treatment to Prevent Heart Attack Trial (ALLHAT) Participants Assigned to Doxazosin and Chlorthalidone

BACKGROUND: The Antihypertensive and Lipid Lowering Treatment to Prevent Heart Attack Trial (ALLHAT) is a randomized, double-blind, active-controlled trial designed to compare the rate of coronary heart disease events in high-risk hypertensive participants initially randomized to a diuretic (chlorthalidone) versus each of three alternative antihypertensive drugs: alpha-adrenergic blocker (doxazosin), ACE-inhibitor (lisinopril), and calcium-channel blocker (amlodipine). Combined cardiovascular disease risk was significantly increased in the doxazosin arm compared to the chlorthalidone arm (RR 1.25; 95% CI, 1.17-1.33; P <.001), with a doubling of heart failure (fatal, hospitalized, or non-hospitalized but treated) (RR 2.04; 95% CI, 1.79-2.32; P <.001). Questions about heart failure diagnostic criteria led to steps to validate these events further. METHODS AND RESULTS: Baseline characteristics (age, race, sex, blood pressure) did not differ significantly between treatment groups (P <.05) for participants with heart failure events. Post-event pharmacologic management was similar in both groups and generally conformed to accepted heart failure therapy. Central review of a small sample of cases showed high adherence to ALLHAT heart failure criteria. Of 105 participants with quantitative ejection fraction measurements provided, (67% by echocardiogram, 31% by catheterization), 29/46 (63%) from the chlorthalidone group and 41/59 (70%) from the doxazosin group were at or below 40%. Two-year heart failure case-fatalities (22% and 19% in the doxazosin and chlorthalidone groups, respectively) were as expected and did not differ significantly (RR 0.96; 95% CI, 0.67-1.38; P = 0.83). CONCLUSION: Results of the validation process supported findings of increased heart failure in the ALLHAT doxazosin treatment arm compared to the chlorthalidone treatment arm.     

Markov chain Monte Carlo for mapping a quantitative trait locus in outbred populations

A Bayesian approach is presented for mapping a quantitative trait locus (QTL) using the 'Fernando and Grossman' multivariate Normal approximation to QTL inheritance. For this model, a Bayesian implementation that includes QTL position is problematic because standard Markov chain Monte Carlo (MCMC) algorithms do not mix, i.e. the QTL position gets stuck in one marker interval. This is because of the dependence of the covariance structure for the QTL effects on the adjacent markers and may be typical of the 'Fernando and Grossman' model. A relatively new MCMC technique, simulated tempering, allows mixing and so makes possible inferences about QTL position based on marginal posterior probabilities. The model was implemented for estimating variance ratios and QTL position using a continuous grid of allowed positions and was applied to simulated data of a standard granddaughter design. The results showed a smooth mixing of QTL position after implementation of the simulated tempering sampler. In this implementation, map distance between QTL and its flanking markers was artificially stretched to reduce the dependence of markers and covariance. The method generalizes easily to more complicated applications and can ultimately contribute to QTL mapping in complex, heterogeneous, human, animal or plant populations.     

Defunct brain stem cardiovascular regulation underlies cardiovascular collapse associated with methamphetamine intoxication

Background  

Intoxication from the psychostimulant methamphetamine (METH) because of cardiovascular collapse is a common cause of death within the abuse population. For obvious reasons, the heart has been taken as the primary target for this METH-induced toxicity. The demonstration that failure of brain stem cardiovascular regulation, rather than the heart, holds the key to cardiovascular collapse induced by the pesticide mevinphos implicates another potential underlying mechanism. The present study evaluated the hypothesis that METH effects acute cardiovascular depression by dampening the functional integrity of baroreflex via an action on brain stem nuclei that are associated with this homeostatic mechanism.     

应用ARIMA模型预测宝安区某街道其它感染性腹泻发病率的探讨

目的:探讨应用ARIMA模型预测宝安区某街道其它感染性腹泻发病率的可行性。方法:应用SPSSl3.0软件对2005年～2009年宝安区某街道其它感染性腹泻逐月发病率进行ARIMA模型建模拟合,用所得到的模型对2010年各月发病率进行预测,并评价其预测效果。结果:宝安区桌街道其它感染性腹泻发病率每年11月为发病高峰,ARIMA(0,1,1)(0,1,0)12模型是其拟合的最佳模型,其预测结果和实际值绝对误差的绝对值最大为930.47,最小为1.96,平均值214.83,平均相对误差百分比39.04%。结论:模型虽然起到一定的预测效果,但预测精度仍存在误差,可通过积累新的周期数据对ARIMA模型进行修正和重新拟合,也可尝试新的预测方法或其他模型,才能加强和保证预测的精度。