AICAR confers prophylactic cardioprotection in doxorubicin-induced heart failure in rats.

Doxorubicin (DOX) is a widely used chemotherapeutic agent that can cause serious cardiotoxic side effects, leading to heart failure (HF). Impaired mitochondrial function is thought to be key factor driving progression into HF. We have previously shown in a rat model of DOX-HF that heart failure with reduced ejection fraction correlates with mitochondrial loss and dysfunction. Adenosine monophosphate-dependent kinase (AMPK) is a cellular energy sensor, regulating mitochondrial biogenesis and energy metabolism, including fatty acid oxidation. We hypothesised that AMPK activation could restore mitochondrial function and therefore be a novel cardioprotective strategy for the prevention of DOX-HF. Consequently, we set out to assess whether 5-aminoimidazole-4-carboxamide 1-ß-D-ribofuranoside (AICAR), an activator of AMPK, could prevent cardiac functional decline in this chronic intravenous rat model of DOX-HF. In line with our hypothesis, AICAR improved cardiac systolic function. AICAR furthermore improved cardiac mitochondrial fatty acid oxidation, independent of mitochondrial number, and in the absence of observable AMPK-activation. In addition, we found that AICAR prevented loss of myocardial mass. RNAseq analysis showed that this may be driven by normalisation of pathways associated with ribosome function and protein synthesis, which are impaired in DOX-treated rat hearts. AICAR furthermore prevented dyslipidemia and excessive body-weight loss in DOX-treated rats, which may contribute to preservation of myocardial mass. Though it is unclear whether AICAR exerted its cardioprotective effect through cardiac or extra-cardiac AMPK-activation or via an AMPK-independent effect, these results show promise for the use of AICAR as a cardioprotective agent in DOX-HF to both preserve cardiac function and mass.

Genomic regulation of Krüppel-like-factor family members by corticosteroid receptors in the rat brain.

Hippocampal mineralocorticoid receptors (MRs) and glucocorticoid receptors (GRs) mediate glucocorticoid hormone (GC) action in the hippocampus. These receptors bind to glucocorticoid responsive elements (GREs) within target genes, eliciting transcriptional effects in response to stress and circadian variation. Until recently, little was known about the genome-wide targets of hippocampal MRs and GRs under physiological conditions. Following on from our genome-wide MR and GR ChIP-seq and Ribo-Zero RNA-seq studies on rat hippocampus, we investigated the Krüppel-like factors (KLFs) as targets of MRs and GRs throughout the brain under circadian variation and after acute stress. In particular, Klf2, Klf9 and Klf15 are known to be stress and/or GC responsive and play a role in neurobiological processes including synaptic plasticity and neuronal differentiation. We found increased binding of MR and GR to GREs within Klf2, Klf9 and Klf15 in the hippocampus, amygdala, prefrontal cortex, and neocortex after acute stress and resulting from circadian variation, which was accompanied by upregulation of corresponding hnRNA and mRNA levels. Adrenalectomy abolished transcriptional upregulation of specific Klf genes. These results show that MRs and GRs regulate Klf gene expression throughout the brain following exposure to acute stress or in response to circadian variation, likely alongside other transcription factors.

miR-3156-5p is downregulated in serum of MEN1 patients and regulates expression of MORF4L2.

Multiple endocrine neoplasia type 1 (MEN1), caused by mutations in the MEN1 gene encoding menin, is an autosomal dominant disorder characterised by the combined occurrence of parathyroid, pituitary and pancreatic neuroendocrine tumours (NETs). Development of these tumours is associated with wide variations in their severity, order and ages (from <5 to >80 years), requiring life-long screening. To improve tumour surveillance and quality of life, better circulating biomarkers, particularly for pancreatic NETs that are associated with higher mortality, are required. We, therefore, examined the expression of circulating miRNA in the serum of MEN1 patients. Initial profiling analysis followed by qRT-PCR validation studies identified miR-3156-5p to be significantly downregulated (-1.3 to 5.8-fold, P < 0.05-0.0005) in nine MEN1 patients, compared to matched unaffected relatives. MEN1 knock-down experiments in BON-1 human pancreatic NET cells resulted in reduced MEN1 (49%, P < 0.05), menin (54%, P < 0.05) and miR-3156-5p expression (20%, P < 0.005), compared to control-treated cells, suggesting that miR-3156-5p downregulation is a consequence of loss of MEN1 expression. In silico analysis identified mortality factor 4-like 2 (MOR4FL2) as a potential target of miR-3156-5p, and in vitro functional studies in BON-1 cells transfected with either miR-3156-5p mimic or inhibitors showed that the miR-3156-5p mimic significantly reduced MORF4L2 protein expression (46%, P < 0.005), while miR-3156-5p inhibitor significantly increased MORF4L2 expression (1.5-fold, P < 0.05), compared to control-treated cells, thereby confirming that miR-3156-5p regulates MORF4L2 expression. Thus, the inverse relationship between miR-3156-5p and MORF4L2 expression represents a potential serum biomarker that could facilitate the detection of NET occurrence in MEN1 patients.

Distinct regulation of hippocampal neuroplasticity and ciliary genes by corticosteroid receptors.

Glucocorticoid hormones (GCs) - acting through hippocampal mineralocorticoid receptors (MRs) and glucocorticoid receptors (GRs) - are critical to physiological regulation and behavioural adaptation. We conducted genome-wide MR and GR ChIP-seq and Ribo-Zero RNA-seq studies on rat hippocampus to elucidate MR- and GR-regulated genes under circadian variation or acute stress. In a subset of genes, these physiological conditions resulted in enhanced MR and/or GR binding to DNA sequences and associated transcriptional changes. Binding of MR at a substantial number of sites however remained unchanged. MR and GR binding occur at overlapping as well as distinct loci. Moreover, although the GC response element (GRE) was the predominant motif, the transcription factor recognition site composition within MR and GR binding peaks show marked differences. Pathway analysis uncovered that MR and GR regulate a substantial number of genes involved in synaptic/neuro-plasticity, cell morphology and development, behavior, and neuropsychiatric disorders. We find that MR, not GR, is the predominant receptor binding to >50 ciliary genes; and that MR function is linked to neuronal differentiation and ciliogenesis in human fetal neuronal progenitor cells. These results show that hippocampal MRs and GRs constitutively and dynamically regulate genomic activities underpinning neuronal plasticity and behavioral adaptation to changing environments.

Decreased ATM Function Causes Delayed DNA Repair and Apoptosis in Common Variable Immunodeficiency Disorders.

PURPOSE: Common variable immunodeficiency disorders (CVID) is characterized by low/absent serum immunoglobulins and susceptibility to bacterial infection. Patients can develop an infections-only phenotype or a complex disease course with inflammatory, autoimmune, and/or malignant complications. We hypothesized that deficient DNA repair mechanisms may be responsible for the antibody deficiency and susceptibility to inflammation and cancer in some patients. METHODS: Germline variants were identified following targeted sequencing of n = 252 genes related to DNA repair in n = 38 patients. NanoString nCounter PlexSet assay measured gene expression in n = 20 CVID patients and n = 7 controls. DNA damage and apoptosis were assessed by flow cytometry in n = 34 CVID patients and n = 11 controls. RESULTS: Targeted sequencing supported enrichment of rare genetic variants in genes related to DNA repair pathways with novel and rare likely pathogenic variants identified and an altered gene expression signature that distinguished patients from controls and complex patients from those with an infections-only phenotype. Consistent with this, flow cytometric analyses of lymphocytes following DNA damage revealed a subset of CVID patients whose immune cells have downregulated ATM, impairing the recruitment of other repair factors, delaying repair and promoting apoptosis. CONCLUSION: These data suggest that germline genetics and altered gene expression predispose a subset of CVID patients to increased sensitivity to DNA damage and reduced DNA repair capacity.

Early detection of doxorubicin-induced cardiotoxicity in rats by its cardiac metabolic signature assessed with hyperpolarized MRI.

Doxorubicin (DOX) is a widely used chemotherapeutic agent that can cause serious cardiotoxic side effects culminating in congestive heart failure (HF). There are currently no clinical imaging techniques or biomarkers available to detect DOX-cardiotoxicity before functional decline. Mitochondrial dysfunction is thought to be a key factor driving functional decline, though real-time metabolic fluxes have never been assessed in DOX-cardiotoxicity. Hyperpolarized magnetic resonance imaging (MRI) can assess real-time metabolic fluxes in vivo. Here we show that cardiac functional decline in a clinically relevant rat-model of DOX-HF is preceded by a change in oxidative mitochondrial carbohydrate metabolism, measured by hyperpolarized MRI. The decreased metabolic fluxes were predominantly due to mitochondrial loss and additional mitochondrial dysfunction, and not, as widely assumed hitherto, to oxidative stress. Since hyperpolarized MRI has been successfully translated into clinical trials this opens up the potential to test cancer patients receiving DOX for early signs of cardiotoxicity.

Author Correction: Nuclear-mitochondrial DNA segments resemble paternally inherited mitochondrial DNA in humans.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Nuclear-mitochondrial DNA segments resemble paternally inherited mitochondrial DNA in humans.

Several strands of evidence question the dogma that human mitochondrial DNA (mtDNA) is inherited exclusively down the maternal line, most recently in three families where several individuals harbored a 'heteroplasmic haplotype' consistent with biparental transmission. Here we report a similar genetic signature in 7 of 11,035 trios, with allelic fractions of 5-25%, implying biparental inheritance of mtDNA in 0.06% of offspring. However, analysing the nuclear whole genome sequence, we observe likely large rare or unique nuclear-mitochondrial DNA segments (mega-NUMTs) transmitted from the father in all 7 families. Independently detecting mega-NUMTs in 0.13% of fathers, we see autosomal transmission of the haplotype. Finally, we show the haplotype allele fraction can be explained by complex concatenated mtDNA-derived sequences rearranged within the nuclear genome. We conclude that rare cryptic mega-NUMTs can resemble paternally mtDNA heteroplasmy, but find no evidence of paternal transmission of mtDNA in humans.

Phenotypic and transcriptomic characterization of canine myeloid-derived suppressor cells.

Myeloid-derived suppressor cells (MDSCs) are key players in immune evasion, tumor progression and metastasis. MDSCs accumulate under various pathological states and fall into two functionally and phenotypically distinct subsets that have been identified in humans and mice: polymorphonuclear (PMN)-MDSCs and monocytic (M)-MDSCs. As dogs are an excellent model for human tumor development and progression, we set out to identify PMN-MDSCs and M-MDSCs in clinical canine oncology patients. Canine hypodense MHC class II-CD5-CD21-CD11b+ cells can be subdivided into polymorphonuclear (CADO48A+CD14-) and monocytic (CADO48A-CD14+) MDSC subsets. The transcriptomic signatures of PMN-MDSCs and M-MDSCs are distinct, and moreover reveal a statistically significant similarity between canine and previously published human PMN-MDSC gene expression patterns. As in humans, peripheral blood frequencies of canine PMN-MDSCs and M-MDSCs are significantly higher in dogs with cancer compared to healthy control dogs (PMN-MDSCs: p < 0.001; M-MDSCs: p < 0.01). By leveraging the power of evolution, we also identified additional conserved genes in PMN-MDSCs of multiple species that may play a role in MDSC function. Our findings therefore validate the dog as a model for studying MDSCs in the context of cancer.

Impact of spliceosome mutations on RNA splicing in myelodysplasia: dysregulated genes/pathways and clinical associations.

SF3B1, SRSF2, and U2AF1 are the most frequently mutated splicing factor genes in the myelodysplastic syndromes (MDS). We have performed a comprehensive and systematic analysis to determine the effect of these commonly mutated splicing factors on pre-mRNA splicing in the bone marrow stem/progenitor cells and in the erythroid and myeloid precursors in splicing factor mutant MDS. Using RNA-seq, we determined the aberrantly spliced genes and dysregulated pathways in CD34+ cells of 84 patients with MDS. Splicing factor mutations result in different alterations in splicing and largely affect different genes, but these converge in common dysregulated pathways and cellular processes, focused on RNA splicing, protein synthesis, and mitochondrial dysfunction, suggesting common mechanisms of action in MDS. Many of these dysregulated pathways and cellular processes can be linked to the known disease pathophysiology associated with splicing factor mutations in MDS, whereas several others have not been previously associated with MDS, such as sirtuin signaling. We identified aberrantly spliced events associated with clinical variables, and isoforms that independently predict survival in MDS and implicate dysregulation of focal adhesion and extracellular exosomes as drivers of poor survival. Aberrantly spliced genes and dysregulated pathways were identified in the MDS-affected lineages in splicing factor mutant MDS. Functional studies demonstrated that knockdown of the mitosis regulators SEPT2 and AKAP8, aberrantly spliced target genes of SF3B1 and SRSF2 mutations, respectively, led to impaired erythroid cell growth and differentiation. This study illuminates the effect of the common spliceosome mutations on the MDS phenotype and provides novel insights into disease pathophysiology.

A practical solution for preserving single cells for RNA sequencing.

The design and implementation of single-cell experiments is often limited by their requirement for fresh starting material. We have adapted a method for histological tissue fixation using dithio-bis(succinimidyl propionate) (DSP), or Lomant's Reagent, to stabilise cell samples for single-cell transcriptomic applications. DSP is a reversible cross-linker of free amine groups that has previously been shown to preserve tissue integrity for histology while maintaining RNA integrity and yield in bulk RNA extractions. Although RNA-seq data from DSP-fixed single cells appears to be prone to characteristic artefacts, such as slightly reduced yield of cDNA and a detectable 3' bias in comparison with fresh cells, cell preservation using DSP does not appear to substantially reduce RNA complexity at the gene level. In addition, there is evidence that instantaneous fixation of cells can reduce inter-cell technical variability. The ability of DSP-fixed cells to retain commonly used dyes, such as propidium iodide, enables the tracking of experimental sub-populations and the recording of cell viability at the point of fixation. Preserving cells using DSP will remove several barriers in the staging of single-cell experiments, including the transport of samples and the scheduling of shared equipment for downstream single-cell isolation and processing.

The nature and nurture of cell heterogeneity: accounting for macrophage gene-environment interactions with single-cell RNA-Seq.

BACKGROUND: Single-cell RNA-Seq can be a valuable and unbiased tool to dissect cellular heterogeneity, despite the transcriptome's limitations in describing higher functional phenotypes and protein events. Perhaps the most important shortfall with transcriptomic 'snapshots' of cell populations is that they risk being descriptive, only cataloging heterogeneity at one point in time, and without microenvironmental context. Studying the genetic ('nature') and environmental ('nurture') modifiers of heterogeneity, and how cell population dynamics unfold over time in response to these modifiers is key when studying highly plastic cells such as macrophages. RESULTS: We introduce the programmable Polaris™ microfluidic lab-on-chip for single-cell sequencing, which performs live-cell imaging while controlling for the culture microenvironment of each cell. Using gene-edited macrophages we demonstrate how previously unappreciated knockout effects of SAMHD1, such as an altered oxidative stress response, have a large paracrine signaling component. Furthermore, we demonstrate single-cell pathway enrichments for cell cycle arrest and APOBEC3G degradation, both associated with the oxidative stress response and altered proteostasis. Interestingly, SAMHD1 and APOBEC3G are both HIV-1 inhibitors ('restriction factors'), with no known co-regulation. CONCLUSION: As single-cell methods continue to mature, so will the ability to move beyond simple 'snapshots' of cell populations towards studying the determinants of population dynamics. By combining single-cell culture, live-cell imaging, and single-cell sequencing, we have demonstrated the ability to study cell phenotypes and microenvironmental influences. It's these microenvironmental components - ignored by standard single-cell workflows - that likely determine how macrophages, for example, react to inflammation and form treatment resistant HIV reservoirs.

Factors influencing success of clinical genome sequencing across a broad spectrum of disorders.

To assess factors influencing the success of whole-genome sequencing for mainstream clinical diagnosis, we sequenced 217 individuals from 156 independent cases or families across a broad spectrum of disorders in whom previous screening had identified no pathogenic variants. We quantified the number of candidate variants identified using different strategies for variant calling, filtering, annotation and prioritization. We found that jointly calling variants across samples, filtering against both local and external databases, deploying multiple annotation tools and using familial transmission above biological plausibility contributed to accuracy. Overall, we identified disease-causing variants in 21% of cases, with the proportion increasing to 34% (23/68) for mendelian disorders and 57% (8/14) in family trios. We also discovered 32 potentially clinically actionable variants in 18 genes unrelated to the referral disorder, although only 4 were ultimately considered reportable. Our results demonstrate the value of genome sequencing for routine clinical diagnosis but also highlight many outstanding challenges.

Clinical whole-genome sequencing in severe early-onset epilepsy reveals new genes and improves molecular diagnosis.

In severe early-onset epilepsy, precise clinical and molecular genetic diagnosis is complex, as many metabolic and electro-physiological processes have been implicated in disease causation. The clinical phenotypes share many features such as complex seizure types and developmental delay. Molecular diagnosis has historically been confined to sequential testing of candidate genes known to be associated with specific sub-phenotypes, but the diagnostic yield of this approach can be low. We conducted whole-genome sequencing (WGS) on six patients with severe early-onset epilepsy who had previously been refractory to molecular diagnosis, and their parents. Four of these patients had a clinical diagnosis of Ohtahara Syndrome (OS) and two patients had severe non-syndromic early-onset epilepsy (NSEOE). In two OS cases, we found de novo non-synonymous mutations in the genes KCNQ2 and SCN2A. In a third OS case, WGS revealed paternal isodisomy for chromosome 9, leading to identification of the causal homozygous missense variant in KCNT1, which produced a substantial increase in potassium channel current. The fourth OS patient had a recessive mutation in PIGQ that led to exon skipping and defective glycophosphatidyl inositol biosynthesis. The two patients with NSEOE had likely pathogenic de novo mutations in CBL and CSNK1G1, respectively. Mutations in these genes were not found among 500 additional individuals with epilepsy. This work reveals two novel genes for OS, KCNT1 and PIGQ. It also uncovers unexpected genetic mechanisms and emphasizes the power of WGS as a clinical tool for making molecular diagnoses, particularly for highly heterogeneous disorders.

Mutations in TCF12, encoding a basic helix-loop-helix partner of TWIST1, are a frequent cause of coronal craniosynostosis.

Craniosynostosis, the premature fusion of the cranial sutures, is a heterogeneous disorder with a prevalence of ∼1 in 2,200 (refs. 1,2). A specific genetic etiology can be identified in ∼21% of cases, including mutations of TWIST1, which encodes a class II basic helix-loop-helix (bHLH) transcription factor, and causes Saethre-Chotzen syndrome, typically associated with coronal synostosis. Using exome sequencing, we identified 38 heterozygous TCF12 mutations in 347 samples from unrelated individuals with craniosynostosis. The mutations predominantly occurred in individuals with coronal synostosis and accounted for 32% and 10% of subjects with bilateral and unilateral pathology, respectively. TCF12 encodes one of three class I E proteins that heterodimerize with class II bHLH proteins such as TWIST1. We show that TCF12 and TWIST1 act synergistically in a transactivation assay and that mice doubly heterozygous for loss-of-function mutations in Tcf12 and Twist1 have severe coronal synostosis. Hence, the dosage of TCF12-TWIST1 heterodimers is critical for normal coronal suture development.

Next-generation sequencing (NGS) as a diagnostic tool for retinal degeneration reveals a much higher detection rate in early-onset disease.

Inherited retinal degeneration (IRD) is a common cause of visual impairment (prevalence ∼1/3500). There is considerable phenotype and genotype heterogeneity, making a specific diagnosis very difficult without molecular testing. We investigated targeted capture combined with next-generation sequencing using Nimblegen 12plex arrays and the Roche 454 sequencing platform to explore its potential for clinical diagnostics in two common types of IRD, retinitis pigmentosa and cone-rod dystrophy. 50 patients (36 unknowns and 14 positive controls) were screened, and pathogenic mutations were identified in 25% of patients in the unknown, with 53% in the early-onset cases. All patients with new mutations detected had an age of onset <21 years and 44% had a family history. Thirty-one percent of mutations detected were novel. A de novo mutation in rhodopsin was identified in one early-onset case without a family history. Bioinformatic pipelines were developed to identify likely pathogenic mutations and stringent criteria were used for assignment of pathogenicity. Analysis of sequencing metrics revealed significant variability in capture efficiency and depth of coverage. We conclude that targeted capture and next-generation sequencing are likely to be very useful in a diagnostic setting, but patients with earlier onset of disease are more likely to benefit from using this strategy. The mutation-detection rate suggests that many patients are likely to have mutations in novel genes.

GREVE: Genomic Recurrent Event ViEwer to assist the identification of patterns across individual cancer samples.

SUMMARY: GREVE has been developed to assist with the identification of recurrent genomic aberrations across cancer samples. The exact characterization of such aberrations remains a challenge despite the availability of increasing amount of data, from SNParray to next-generation sequencing. Furthermore, genomic aberrations in cancer are especially difficult to handle because they are, by nature, unique to the patients. However, their recurrence in specific regions of the genome has been shown to reflect their relevance in the development of tumors. GREVE makes use of previously characterized events to identify such regions and focus any further analysis. AVAILABILITY: GREVE is available through a web interface and open-source application (http://www.well.ox.ac.uk/GREVE).

Coexpression network analysis in abdominal and gluteal adipose tissue reveals regulatory genetic loci for metabolic syndrome and related phenotypes.

Metabolic Syndrome (MetS) is highly prevalent and has considerable public health impact, but its underlying genetic factors remain elusive. To identify gene networks involved in MetS, we conducted whole-genome expression and genotype profiling on abdominal (ABD) and gluteal (GLU) adipose tissue, and whole blood (WB), from 29 MetS cases and 44 controls. Co-expression network analysis for each tissue independently identified nine, six, and zero MetS-associated modules of coexpressed genes in ABD, GLU, and WB, respectively. Of 8,992 probesets expressed in ABD or GLU, 685 (7.6%) were expressed in ABD and 51 (0.6%) in GLU only. Differential eigengene network analysis of 8,256 shared probesets detected 22 shared modules with high preservation across adipose depots (D(ABD-GLU) = 0.89), seven of which were associated with MetS (FDR P<0.01). The strongest associated module, significantly enriched for immune response-related processes, contained 94/620 (15%) genes with inter-depot differences. In an independent cohort of 145/141 twins with ABD and WB longitudinal expression data, median variability in ABD due to familiality was greater for MetS-associated versus un-associated modules (ABD: 0.48 versus 0.18, P = 0.08; GLU: 0.54 versus 0.20, P = 7.8×10(-4)). Cis-eQTL analysis of probesets associated with MetS (FDR P<0.01) and/or inter-depot differences (FDR P<0.01) provided evidence for 32 eQTLs. Corresponding eSNPs were tested for association with MetS-related phenotypes in two GWAS of >100,000 individuals; rs10282458, affecting expression of RARRES2 (encoding chemerin), was associated with body mass index (BMI) (P = 6.0×10(-4)); and rs2395185, affecting inter-depot differences of HLA-DRB1 expression, was associated with high-density lipoprotein (P = 8.7×10(-4)) and BMI-adjusted waist-to-hip ratio (P = 2.4×10(-4)). Since many genes and their interactions influence complex traits such as MetS, integrated analysis of genotypes and coexpression networks across multiple tissues relevant to clinical traits is an efficient strategy to identify novel associations.

A fine-scale chimpanzee genetic map from population sequencing.

To study the evolution of recombination rates in apes, we developed methodology to construct a fine-scale genetic map from high-throughput sequence data from 10 Western chimpanzees, Pan troglodytes verus. Compared to the human genetic map, broad-scale recombination rates tend to be conserved, but with exceptions, particularly in regions of chromosomal rearrangements and around the site of ancestral fusion in human chromosome 2. At fine scales, chimpanzee recombination is dominated by hotspots, which show no overlap with those of humans even though rates are similarly elevated around CpG islands and decreased within genes. The hotspot-specifying protein PRDM9 shows extensive variation among Western chimpanzees, and there is little evidence that any sequence motifs are enriched in hotspots. The contrasting locations of hotspots provide a natural experiment, which demonstrates the impact of recombination on base composition.

The use of genome-wide eQTL associations in lymphoblastoid cell lines to identify novel genetic pathways involved in complex traits.

The integrated analysis of genotypic and expression data for association with complex traits could identify novel genetic pathways involved in complex traits. We profiled 19,573 expression probes in Epstein-Barr virus-transformed lymphoblastoid cell lines (LCLs) from 299 twins and correlated these with 44 quantitative traits (QTs). For 939 expressed probes correlating with more than one QT, we investigated the presence of eQTL associations in three datasets of 57 CEU HapMap founders and 86 unrelated twins. Genome-wide association analysis of these probes with 2.2 m SNPs revealed 131 potential eQTLs (1,989 eQTL SNPs) overlapping between the HapMap datasets, five of which were in cis (58 eQTL SNPs). We then tested 535 SNPs tagging the eQTL SNPs, for association with the relevant QT in 2,905 twins. We identified nine potential SNP-QT associations (P<0.01) but none significantly replicated in five large consortia of 1,097-16,129 subjects. We also failed to replicate previous reported eQTL associations with body mass index, plasma low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides levels derived from lymphocytes, adipose and liver tissue. Our results and additional power calculations suggest that proponents may have been overoptimistic in the power of LCLs in eQTL approaches to elucidate regulatory genetic effects on complex traits using the small datasets generated to date. Nevertheless, larger tissue-specific expression data sets relevant to specific traits are becoming available, and should enable the adoption of similar integrated analyses in the near future.