Dysregulated COMT Expression in Fragile X Syndrome

Transcriptional and proteomics analyses in human fragile X syndrome (FXS) neurons identified markedly reduced expression of COMT, a key enzyme involved in the metabolism of catecholamines, including dopamine, epinephrine and norepinephrine. FXS is the most common genetic cause of intellectual disability and autism spectrum disorders. COMT encodes for catechol-o-methyltransferase and its association with neuropsychiatric disorders and cognitive function has been extensively studied. We observed a significantly reduced level of COMT in in FXS human neural progenitors and neurons, as well as hippocampal neurons from Fmr1 null mice. We show that deficits in COMT were associated with an altered response in an assay of dopaminergic activity in Fmr1 null mice. These findings demonstrate that loss of FMRP downregulates COMT expression and affects dopamine signaling in FXS, and supports the notion that targeting catecholamine metabolism may be useful in regulating certain neuropsychiatric aspects of FXS.

Inhibition of miR-199a-3p in a murine hypertrophic cardiomyopathy (HCM) model attenuates fibrotic remodeling

Background

Hypertrophic cardiomyopathy (HCM) is an autosomal dominant genetic disorder, characterized by cardiomyocyte hypertrophy, cardiomyocyte disarray and fibrosis, which has a prevalence of ∼1: 200-500 and predisposes individuals to heart failure and sudden death. The mechanisms through which diverse HCM-causing mutations cause cardiac dysfunction remain mostly unknown and their identification may reveal new therapeutic avenues. MicroRNAs (miRNAs) have emerged as critical regulators of gene expression and disease phenotype in various pathologies. We explored whether miRNAs could play a role in HCM pathogenesis and offer potential therapeutic targets.

Methods and results

Using high-throughput miRNA expression profiling and qPCR analysis in two distinct mouse models of HCM, we found that miR-199a-3p expression levels are upregulated in mutant mice compared to age- and treatment-matched wild-type mice. We also found that miR-199a-3p expression is enriched in cardiac non-myocytes compared to cardiomyocytes. When we expressed miR-199a-3p mimic in cultured murine primary cardiac fibroblasts and analyzed the conditioned media by proteomics, we found that several extracellular matrix (ECM) proteins (e.g., TSP2, FBLN3, COL11A1, LYOX) were differentially secreted (data are available via ProteomeXchange with identifier PXD042904). We confirmed our proteomics findings by qPCR analysis of selected mRNAs and demonstrated that miR-199a-3p mimic expression in cardiac fibroblasts drives upregulation of ECM gene expression, including Tsp2, Fbln3, Pcoc1, Col1a1 and Col3a1. To examine the role of miR-199a-3p in vivo, we inhibited its function using lock-nucleic acid (LNA)-based inhibitors (antimiR-199a-3p) in an HCM mouse model. Our results revealed that progression of cardiac fibrosis is attenuated when miR-199a-3p function is inhibited in mild-to-moderate HCM. Finally, guided by computational target prediction algorithms, we identified mRNAs Cd151 and Itga3 as direct targets of miR-199a-3p and have shown that miR-199a-3p mimic expression negatively regulates AKT activation in cardiac fibroblasts.

Conclusions

Altogether, our results suggest that miR-199a-3p may contribute to cardiac fibrosis in HCM through its actions in cardiac fibroblasts. Thus, inhibition of miR-199a-3p in mild-to-moderate HCM may offer therapeutic benefit in combination with complementary approaches that target the primary defect in cardiac myocytes.

Widespread dysregulation of mRNA splicing implicates RNA processing in the development and progression of Huntington's disease

BACKGROUND

In Huntington's disease (HD), a CAG repeat expansion mutation in the Huntingtin (HTT) gene drives a gain-of-function toxicity that disrupts mRNA processing. Although dysregulation of gene splicing has been shown in human HD post-mortem brain tissue, post-mortem analyses are likely confounded by cell type composition changes in late-stage HD, limiting the ability to identify dysregulation related to early pathogenesis.

METHODS

To investigate gene splicing changes in early HD, we performed alternative splicing analyses coupled with a proteogenomics approach to identify early CAG length-associated splicing changes in an established isogenic HD cell model.

FINDINGS

We report widespread neuronal differentiation stage- and CAG length-dependent splicing changes, and find an enrichment of RNA processing, neuronal function, and epigenetic modification-related genes with mutant HTT-associated splicing. When integrated with a proteomics dataset, we identified several of these differential splicing events at the protein level. By comparing with human post-mortem and mouse model data, we identified common patterns of altered splicing from embryonic stem cells through to post-mortem striatal tissue.

INTERPRETATION

We show that widespread splicing dysregulation in HD occurs in an early cell model of neuronal development. Importantly, we observe HD-associated splicing changes in our HD cell model that were also identified in human HD striatum and mouse model HD striatum, suggesting that splicing-associated pathogenesis possibly occurs early in neuronal development and persists to later stages of disease. Together, our results highlight splicing dysregulation in HD which may lead to disrupted neuronal function and neuropathology.

FUNDING

This research is supported by the Lee Kong Chian School of Medicine, Nanyang Technological University Singapore Nanyang Assistant Professorship Start-Up Grant, the Singapore Ministry of Education under its Singapore Ministry of Education Academic Research Fund Tier 1 (RG23/22), the BC Children's Hospital Research Institute Investigator Grant Award (IGAP), and a Scholar Award from the Michael Smith Health Research BC.

Mitochondrial control of microglial phagocytosis by the translocator protein and hexokinase 2 in Alzheimer's disease

Microglial phagocytosis is an energetically demanding process that plays a critical role in the removal of toxic protein aggregates in Alzheimer's disease (AD). Recent evidence indicates that a switch in energy production from mitochondrial respiration to glycolysis disrupts this important protective microglial function and may provide therapeutic targets for AD. Here, we demonstrate that the translocator protein (TSPO) and a member of its mitochondrial complex, hexokinase-2 (HK), play critical roles in microglial respiratory-glycolytic metabolism and phagocytosis. Pharmacological and genetic loss-of-function experiments showed that TSPO is critical for microglial respiratory metabolism and energy supply for phagocytosis, and its expression is enriched in phagocytic microglia of AD mice. Meanwhile, HK controlled glycolytic metabolism and phagocytosis via mitochondrial binding or displacement. In cultured microglia, TSPO deletion impaired mitochondrial respiration and increased mitochondrial recruitment of HK, inducing a switch to glycolysis and reducing phagocytosis. To determine the functional significance of mitochondrial HK recruitment, we developed an optogenetic tool for reversible control of HK localization. Displacement of mitochondrial HK inhibited glycolysis and improved phagocytosis in TSPO-knockout microglia. Mitochondrial HK recruitment also coordinated the inflammatory switch to glycolysis that occurs in response to lipopolysaccharide in normal microglia. Interestingly, cytosolic HK increased phagocytosis independent of its metabolic activity, indicating an immune signaling function. Alzheimer's beta amyloid drastically stimulated mitochondrial HK recruitment in cultured microglia, which may contribute to microglial dysfunction in AD. Thus, targeting mitochondrial HK may offer an immunotherapeutic approach to promote phagocytic microglial function in AD.

A high-resolution map of human RNA translation

Translated small open reading frames (smORFs) can have important regulatory roles and encode microproteins, yet their genome-wide identification has been challenging. We determined the ribosome locations across six primary human cell types and five tissues and detected 7,767 smORFs with translational profiles matching those of known proteins. The human genome was found to contain highly cell-type- and tissue-specific smORFs and a subset that encodes highly conserved amino acid sequences. Changes in the translational efficiency of upstream-encoded smORFs (uORFs) and the corresponding main ORFs predominantly occur in the same direction. Integration with 456 mass-spectrometry datasets confirms the presence of 603 small peptides at the protein level in humans and provides insights into the subcellular localization of these small proteins. This study provides a comprehensive atlas of high-confidence translated smORFs derived from primary human cells and tissues in order to provide a more complete understanding of the translated human genome.

Ermin deficiency leads to compromised myelin, inflammatory milieu, and susceptibility to demyelinating insult

Ermin is an actin-binding protein found almost exclusively in the central nervous system (CNS) as a component of myelin sheaths. Although Ermin has been predicted to play a role in the formation and stability of myelin sheaths, this has not been directly examined in vivo. Here, we show that Ermin is essential for myelin sheath integrity and normal saltatory conduction. Loss of Ermin in mice caused de-compacted and fragmented myelin sheaths and led to slower conduction along with progressive neurological deficits. RNA sequencing of the corpus callosum, the largest white matter structure in the CNS, pointed to inflammatory activation in aged Ermin-deficient mice, which was corroborated by increased levels of microgliosis and astrogliosis. The inflammatory milieu and myelin abnormalities were further associated with increased susceptibility to immune-mediated demyelination insult in Ermin knockout mice. Supporting a possible role of Ermin deficiency in inflammatory white matter disorders, a rare inactivating mutation in the ERMN gene was identified in multiple sclerosis patients. Our findings demonstrate a critical role for Ermin in maintaining myelin integrity. Given its near-exclusive expression in myelinating oligodendrocytes, Ermin deficiency represents a compelling "inside-out" model of inflammatory dysmyelination and may offer a new paradigm for the development of myelin stability-targeted therapies.

Mitchell-Riley syndrome iPSCs exhibit reduced pancreatic endoderm differentiation due to a mutation in RFX6

Mitchell-Riley syndrome (MRS) is caused by recessive mutations in the regulatory factor X6 gene (RFX6) and is characterised by pancreatic hypoplasia and neonatal diabetes. To determine why individuals with MRS specifically lack pancreatic endocrine cells, we micro-CT imaged a 12-week-old foetus homozygous for the nonsense mutation RFX6 c.1129C>T, which revealed loss of the pancreas body and tail. From this foetus, we derived iPSCs and show that differentiation of these cells in vitro proceeds normally until generation of pancreatic endoderm, which is significantly reduced. We additionally generated an RFX6HA reporter allele by gene targeting in wild-type H9 cells to precisely define RFX6 expression and in parallel performed in situ hybridisation for RFX6 in the dorsal pancreatic bud of a Carnegie stage 14 human embryo. Both in vitro and in vivo, we find that RFX6 specifically labels a subset of PDX1-expressing pancreatic endoderm. In summary, RFX6 is essential for efficient differentiation of pancreatic endoderm, and its absence in individuals with MRS specifically impairs formation of endocrine cells of the pancreas head and tail.

TLR4 signals through islet macrophages to alter cytokine secretion during diabetes

Toll-like receptors (TLRs), particularly TLR4, may act as immune sensors for metabolic stress signals such as lipids and link tissue metabolic changes to innate immunity. TLR signalling is not only tissue-dependent but also cell-type dependent and recent studies suggest that TLRs are not restricted to innate immune cells alone. Pancreatic islets, a hub of metabolic hormones and cytokines, respond to TLR signalling. However, the source of TLR signalling within the islet remain poorly understood. Uncovering the specific cell source and its role in mediating TLR signalling, especially within type 2 diabetes (T2D) islet will yield new targets to tackle islet inflammation, hormone secretion dysregulation and ultimately diabetes. In the present study, we immuno-characterised TLRs linked to pancreatic islets in both healthy and obese diabetic mice. We found that while TLRs1-4 and TLR9 were expressed in mouse islets, these TLRs did not co-localise with insulin-producing β-cells. β-Cells from obese diabetic mice were also devoid of these TLRs. While TLR immunoreactivity in obese mice islets increased, this was driven mostly by increased islet endothelial cell and islet macrophage presence. Analysis of human islet single-cell RNA-seq databases revealed that macrophages were an important source of islet TLRs. However, only TLR4 and TLR8 showed variation and cell-type specificity in their expression patterns. Cell depletion experiments in isolated mouse islets showed that TLR4 signalled through macrophages to alter islet cytokine secretome. Together, these studies suggest that islet macrophages are a dominant source of TLR4-mediated signalling in both healthy and diabetic islets.

Integrated systems-genetic analyses reveal a network target for delaying glioma progression

Objective

To identify a convergent, multitarget proliferation characteristic for astrocytoma transformation that could be targeted for therapy discovery.

Methods

Using an integrated functional genomics approach, we prioritized networks associated with astrocytoma progression using the following criteria: differential co‐expression between grade II and grade III IDH1‐mutated and 1p/19q euploid astrocytomas, preferential enrichment for genetic risk to cancer, association with patient survival and sample‐level genomic features. Drugs targeting the identified multitarget network characteristic for astrocytoma transformation were computationally predicted using drug transcriptional perturbation data and validated using primary human astrocytoma cells.

Results

A single network, M2, consisting of 177 genes, was associated with glioma progression on the basis of the above criteria. Functionally, M2 encoded physically interacting proteins regulating cell cycle processes and analysis of genome‐wide gene‐regulatory interactions using mutual information and DNA–protein interactions revealed the known regulators of cell cycle processes FoxM1, B‐Myb, and E2F2 as key regulators of M2. These results suggest functional disruption of M2 via gene mutation or altered expression as a convergent pathway regulating astrocytoma transformation. By considering M2 as a multitarget drug target regulating astrocytoma transformation, we identified several drugs that are predicted to restore M2 expression in anaplastic astrocytoma toward its low‐grade profile and of these, we validated the known antiproliferative drug resveratrol as down‐regulating multiple nodes of M2 including at nanomolar concentrations achievable in human cerebrospinal fluid by oral dosing.

Interpretation

Our results identify M2 as a multitarget network characteristic for astrocytoma progression and encourage M2‐based drug screening to identify new compounds for preventing glioma transformation.

Widespread Translational Control of Fibrosis in the Human Heart by RNA-Binding Proteins

Supplemental Digital Content is available in the text.

Fibrosis is a common pathology in many cardiac disorders and is driven by the activation of resident fibroblasts. The global posttranscriptional mechanisms underlying fibroblast-to-myofibroblast conversion in the heart have not been explored.

Ethnicity-Specific Skeletal Muscle Transcriptional Signatures and Their Relevance to Insulin Resistance in Singapore

CONTEXT

Insulin resistance (IR) and obesity differ among ethnic groups in Singapore, with the Malays more obese yet less IR than Asian-Indians. However, the molecular basis underlying these differences is not clear.

OBJECTIVE

As the skeletal muscle (SM) is metabolically relevant to IR, we investigated molecular pathways in SM that are associated with ethnic differences in IR, obesity, and related traits.

DESIGN, SETTING, AND MAIN OUTCOME MEASURES

We integrated transcriptomic, genomic, and phenotypic analyses in 156 healthy subjects representing three major ethnicities in the Singapore Adult Metabolism Study.

PATIENTS

This study contains Chinese (n = 63), Malay (n = 51), and Asian-Indian (n = 42) men, aged 21 to 40 years, without systemic diseases.

RESULTS

We found remarkable diversity in the SM transcriptome among the three ethnicities, with >8000 differentially expressed genes (40% of all genes expressed in SM). Comparison with blood transcriptome from a separate Singaporean cohort showed that >95% of SM expression differences among ethnicities were unique to SM. We identified a network of 46 genes that were specifically downregulated in Malays, suggesting dysregulation of components of cellular respiration in SM of Malay individuals. We also report 28 differentially expressed gene clusters, four of which were also enriched for genes that were found in genome-wide association studies of metabolic traits and disease and correlated with variation in IR, obesity, and related traits.

CONCLUSION

We identified extensive gene-expression changes in SM among the three Singaporean ethnicities and report specific genes and molecular pathways that might underpin and explain the differences in IR among these ethnic groups.

Inhibition of profibrotic microRNA-21 affects platelets and their releasate

Fibrosis is a major contributor to organ disease for which no specific therapy is available. MicroRNA-21 (miR-21) has been implicated in the fibrogenetic response, and inhibitors of miR-21 are currently undergoing clinical trials. Here, we explore how miR-21 inhibition may attenuate fibrosis using a proteomics approach. Transfection of miR-21 mimic or inhibitor in murine cardiac fibroblasts revealed limited effects on extracellular matrix (ECM) protein secretion. Similarly, miR-21–null mouse hearts showed an unaltered ECM composition. Thus, we searched for additional explanations as to how miR-21 might regulate fibrosis. In plasma samples from the community-based Bruneck Study, we found a marked correlation of miR-21 levels with several platelet-derived profibrotic factors, including TGF-β1. Pharmacological miR-21 inhibition with an antagomiR reduced the platelet release of TGF-β1 in mice. Mechanistically, Wiskott-Aldrich syndrome protein, a negative regulator of platelet TGF-β1 secretion, was identified as a direct target of miR-21. miR-21–null mice had lower platelet and leukocyte counts compared with littermate controls but higher megakaryocyte numbers in the bone marrow. Thus, to our knowledge this study reports a previously unrecognized effect of miR-21 inhibition on platelets. The effect of antagomiR-21 treatment on platelet TGF-β1 release, in particular, may contribute to the antifibrotic effects of miR-21 inhibitors.

MicroRNA-21 inhibition may convey its therapeutic benefits in fibrosis through its action in bone marrow cells rather than targeting fibroblasts directly.

Genome-wide analysis of differential RNA editing in epilepsy

The recoding of genetic information through RNA editing contributes to proteomic diversity, but the extent and significance of RNA editing in disease is poorly understood. In particular, few studies have investigated the relationship between RNA editing and disease at a genome-wide level. Here, we developed a framework for the genome-wide detection of RNA sites that are differentially edited in disease. Using RNA-sequencing data from 100 hippocampi from mice with epilepsy (pilocarpine–temporal lobe epilepsy model) and 100 healthy control hippocampi, we identified 256 RNA sites (overlapping with 87 genes) that were significantly differentially edited between epileptic cases and controls. The degree of differential RNA editing in epileptic mice correlated with frequency of seizures, and the set of genes differentially RNA-edited between case and control mice were enriched for functional terms highly relevant to epilepsy, including “neuron projection” and “seizures.” Genes with differential RNA editing were preferentially enriched for genes with a genetic association to epilepsy. Indeed, we found that they are significantly enriched for genes that harbor nonsynonymous de novo mutations in patients with epileptic encephalopathy and for common susceptibility variants associated with generalized epilepsy. These analyses reveal a functional convergence between genes that are differentially RNA-edited in acquired symptomatic epilepsy and those that contribute risk for genetic epilepsy. Taken together, our results suggest a potential role for RNA editing in the epileptic hippocampus in the occurrence and severity of epileptic seizures.

Extracellular matrix proteomics identifies molecular signature of symptomatic carotid plaques

BACKGROUND. The identification of patients with high-risk atherosclerotic plaques prior to the manifestation of clinical events remains challenging. Recent findings question histology- and imaging-based definitions of the “vulnerable plaque,” necessitating an improved approach for predicting onset of symptoms.

METHODS. We performed a proteomics comparison of the vascular extracellular matrix and associated molecules in human carotid endarterectomy specimens from 6 symptomatic versus 6 asymptomatic patients to identify a protein signature for high-risk atherosclerotic plaques. Proteomics data were integrated with gene expression profiling of 121 carotid endarterectomies and an analysis of protein secretion by lipid-loaded human vascular smooth muscle cells. Finally, epidemiological validation of candidate biomarkers was performed in two community-based studies.

RESULTS. Proteomics and at least one of the other two approaches identified a molecular signature of plaques from symptomatic patients that comprised matrix metalloproteinase 9, chitinase 3-like-1, S100 calcium binding protein A8 (S100A8), S100A9, cathepsin B, fibronectin, and galectin-3-binding protein. Biomarker candidates measured in 685 subjects in the Bruneck study were associated with progression to advanced atherosclerosis and incidence of cardiovascular disease over a 10-year follow-up period. A 4-biomarker signature (matrix metalloproteinase 9, S100A8/S100A9, cathepsin D, and galectin-3-binding protein) improved risk prediction and was successfully replicated in an independent cohort, the SAPHIR study.

CONCLUSION. The identified 4-biomarker signature may improve risk prediction and diagnostics for the management of cardiovascular disease. Further, our study highlights the strength of tissue-based proteomics for biomarker discovery.

FUNDING. UK: British Heart Foundation (BHF); King’s BHF Center; and the National Institute for Health Research Biomedical Research Center based at Guy’s and St Thomas’ NHS Foundation Trust and King’s College London in partnership with King’s College Hospital. Austria: Federal Ministry for Transport, Innovation and Technology (BMVIT); Federal Ministry of Science, Research and Economy (BMWFW); Wirtschaftsagentur Wien; and Standortagentur Tirol.

Rare and common epilepsies converge on a shared gene regulatory network providing opportunities for novel antiepileptic drug discovery

BACKGROUND

The relationship between monogenic and polygenic forms of epilepsy is poorly understood and the extent to which the genetic and acquired epilepsies share common pathways is unclear. Here, we use an integrated systems-level analysis of brain gene expression data to identify molecular networks disrupted in epilepsy.

RESULTS

We identified a co-expression network of 320 genes (M30), which is significantly enriched for non-synonymous de novo mutations ascertained from patients with monogenic epilepsy and for common variants associated with polygenic epilepsy. The genes in the M30 network are expressed widely in the human brain under tight developmental control and encode physically interacting proteins involved in synaptic processes. The most highly connected proteins within the M30 network were preferentially disrupted by deleterious de novo mutations for monogenic epilepsy, in line with the centrality-lethality hypothesis. Analysis of M30 expression revealed consistent downregulation in the epileptic brain in heterogeneous forms of epilepsy including human temporal lobe epilepsy, a mouse model of acquired temporal lobe epilepsy, and a mouse model of monogenic Dravet (SCN1A) disease. These results suggest functional disruption of M30 via gene mutation or altered expression as a convergent mechanism regulating susceptibility to epilepsy broadly. Using the large collection of drug-induced gene expression data from Connectivity Map, several drugs were predicted to preferentially restore the downregulation of M30 in epilepsy toward health, most notably valproic acid, whose effect on M30 expression was replicated in neurons.

CONCLUSIONS

Taken together, our results suggest targeting the expression of M30 as a potential new therapeutic strategy in epilepsy.

"Young at heart": Regenerative potential linked to immature cardiac phenotypes

The adult human myocardium is incapable of regeneration; yet, the zebrafish (Danio rerio) can regenerate damaged myocardium. Similar to the zebrafish heart, hearts of neonatal, but not adult mice are capable of myocardial regeneration. We performed a proteomics analysis of adult zebrafish hearts and compared their protein expression profile to hearts from neonatal and adult mice. Using difference in-gel electrophoresis (DIGE), there was little overlap between the proteome from adult mouse (> 8 weeks old) and adult zebrafish (18 months old) hearts. Similarly, there was a significant degree of mismatch between the protein expression in neonatal and adult mouse hearts. Enrichment analysis of the selected proteins revealed over-expression of DNA synthesis-related proteins in the cardiac proteome of the adult zebrafish heart similar to neonatal and 4 days old mice, whereas in hearts of adult mice there was a mitochondria-related predominance in protein expression. Importantly, we noted pronounced differences in the myofilament composition: the adult zebrafish heart lacks many of the myofilament proteins of differentiated adult cardiomyocytes such as the ventricular isoforms of myosin light chains and nebulette. Instead, troponin I and myozenin 1 were expressed as skeletal isoforms rather than cardiac isoforms. The relative immaturity of the adult zebrafish heart was further supported by cardiac microRNA data. Our assessment of zebrafish and mammalian hearts challenges the assertions on the translational potential of cardiac regeneration in the zebrafish model. The immature myofilament composition of the fish heart may explain why adult mouse and human cardiomyocytes lack this endogenous repair mechanism.

•

Proteomics reveals minimal overlap between adult mouse and adult zebrafish hearts.

•

Gene expression analysis confirms profound differences in myofilament composition.

•

The adult zebrafish heart is more similar to a newborn mouse heart.

•

The relative immaturity is further supported by cardiac microRNA data.

Angiogenic microRNAs Linked to Incidence and Progression of Diabetic Retinopathy in Type 1 Diabetes

Circulating microRNAs (miRNAs) have emerged as novel biomarkers of diabetes. The current study focuses on the role of circulating miRNAs in patients with type 1 diabetes and their association with diabetic retinopathy. A total of 29 miRNAs were quantified in serum samples (n = 300) using a nested case-control study design in two prospective cohorts of the DIabetic REtinopathy Candesartan Trial (DIRECT): PROTECT-1 and PREVENT-1. The PREVENT-1 trial included patients without retinopathy at baseline; the PROTECT-1 trial included patients with nonproliferative retinopathy at baseline. Two miRNAs previously implicated in angiogenesis, miR-27b and miR-320a, were associated with incidence and with progression of retinopathy: the odds ratio per SD higher miR-27b was 0.57 (95% CI 0.40, 0.82; P = 0.002) in PREVENT-1, 0.78 (0.57, 1.07; P = 0.124) in PROTECT-1, and 0.67 (0.50, 0.92; P = 0.012) combined. The respective odds ratios for higher miR-320a were 1.57 (1.07, 2.31; P = 0.020), 1.43 (1.05, 1.94; P = 0.021), and 1.48 (1.17, 1.88; P = 0.001). Proteomics analyses in endothelial cells returned the antiangiogenic protein thrombospondin-1 as a common target of both miRNAs. Our study identifies two angiogenic miRNAs, miR-320a and miR-27b, as potential biomarkers for diabetic retinopathy.

Association of MicroRNAs and YRNAs With Platelet Function

RATIONALE

Platelets shed microRNAs (miRNAs). Plasma miRNAs change on platelet inhibition. It is unclear whether plasma miRNA levels correlate with platelet function.

OBJECTIVE

To link small RNAs to platelet reactivity.

METHODS AND RESULTS

Next-generation sequencing of small RNAs in plasma revealed 2 peaks at 22 to 23 and 32 to 33 nucleotides corresponding to miRNAs and YRNAs, respectively. Among YRNAs, predominantly, fragments of RNY4 and RNY5 were detected. Plasma miRNAs and YRNAs were measured in 125 patients with a history of acute coronary syndrome who had undergone detailed assessment of platelet function 30 days after the acute event. Using quantitative real-time polymerase chain reactions, 92 miRNAs were assessed in patients with acute coronary syndrome on different antiplatelet therapies. Key platelet-related miRNAs and YRNAs were correlated with platelet function tests. MiR-223 (rp=0.28; n=121; P=0.002), miR-126 (rp=0.22; n=121; P=0.016), and other abundant platelet miRNAs and YRNAs showed significant positive correlations with the vasodilator-stimulated phosphoprotein phosphorylation assay. YRNAs, miR-126, and miR-223 were also among the small RNAs showing the greatest dependency on platelets and strongly correlated with plasma levels of P-selectin, platelet factor 4, and platelet basic protein in the population-based Bruneck study (n=669). A single-nucleotide polymorphism that facilitates processing of pri-miR-126 to mature miR-126 accounted for a rise in circulating platelet activation markers. Inhibition of miR-126 in mice reduced platelet aggregation. MiR-126 directly and indirectly affects ADAM9 and P2Y12 receptor expression.

CONCLUSIONS

Levels of platelet-related plasma miRNAs and YRNAs correlate with platelet function tests in patients with acute coronary syndrome and platelet activation markers in the general population. Alterations in miR-126 affect platelet reactivity.

Comparative analysis of statistical methods used for detecting differential expression in label-free mass spectrometry proteomics

Label-free LC-MS/MS proteomics has proven itself to be a powerful method for evaluating protein identification and quantification from complex samples. For comparative proteomics, several methods have been used to detect the differential expression of proteins from such data. We have assessed seven methods used across the literature for detecting differential expression from spectral count quantification: Student's t-test, significance analysis of microarrays (SAM), normalised spectral abundance factor (NSAF), normalised spectral abundance factor-power law global error model (NSAF-PLGEM), spectral index (SpI), DESeq and QSpec. We used 2000 simulated datasets as well as publicly available data from a proteomic standards study to assess the ability of these methods to detect differential expression in varying effect sizes and proportions of differentially expressed proteins. At two false discovery rate (FDR) levels, we find that several of the methods detect differential expression within the data with reasonable precision, others detect differential expression at the expense of low precision, and finally, others which fail to identify any differentially expressed proteins. The inability of these seven methods to fully capture the differential landscape, even at the largest effect size, illustrates some of the limitations of the existing technologies and the statistical methodologies.

SIGNIFICANCE

In label-free mass spectrometry experiments, protein identification and quantification have always been important, but there is now a growing focus on comparative proteomics. Detecting differential expression in protein levels can inform on important biological mechanisms and provide direction for further study. Given the high cost and labour intensive nature of validation experiments, statistical methods are important for prioritising proteins of interest. Here, we have performed a comparative analysis to investigate the statistical methodologies for detecting differential expression and provide a reference for future experimental designs. This article is part of a Special Issue entitled: Computational Proteomics.

MicroRNA profiles in hippocampal granule cells and plasma of rats with pilocarpine-induced epilepsy--comparison with human epileptic samples

The identification of biomarkers of the transformation of normal to epileptic tissue would help to stratify patients at risk of epilepsy following brain injury, and inform new treatment strategies. MicroRNAs (miRNAs) are an attractive option in this direction. In this study, miRNA microarrays were performed on laser-microdissected hippocampal granule cell layer (GCL) and on plasma, at different time points in the development of pilocarpine-induced epilepsy in the rat: latency, first spontaneous seizure and chronic epileptic phase. Sixty-three miRNAs were differentially expressed in the GCL when considering all time points. Three main clusters were identified that separated the control and chronic phase groups from the latency group and from the first spontaneous seizure group. MiRNAs from rats in the chronic phase were compared to those obtained from the laser-microdissected GCL of epileptic patients, identifying several miRNAs (miR-21-5p, miR-23a-5p, miR-146a-5p and miR-181c-5p) that were up-regulated in both human and rat epileptic tissue. Analysis of plasma samples revealed different levels between control and pilocarpine-treated animals for 27 miRNAs. Two main clusters were identified that segregated controls from all other groups. Those miRNAs that are altered in plasma before the first spontaneous seizure, like miR-9a-3p, may be proposed as putative biomarkers of epileptogenesis.

Long-term therapeutic silencing of miR-33 increases circulating triglyceride levels and hepatic lipid accumulation in mice

Plasma high-density lipoprotein (HDL) levels show a strong inverse correlation with atherosclerotic vascular disease. Previous studies have demonstrated that antagonism of miR-33 in vivo increases circulating HDL and reverse cholesterol transport (RCT), thereby reducing the progression and enhancing the regression of atherosclerosis. While the efficacy of short-term anti-miR-33 treatment has been previously studied, the long-term effect of miR-33 antagonism in vivo remains to be elucidated. Here, we show that long-term therapeutic silencing of miR-33 increases circulating triglyceride (TG) levels and lipid accumulation in the liver. These adverse effects were only found when mice were fed a high-fat diet (HFD). Mechanistically, we demonstrate that chronic inhibition of miR-33 increases the expression of genes involved in fatty acid synthesis such as acetyl-CoA carboxylase (ACC) and fatty acid synthase (FAS) in the livers of mice treated with miR-33 antisense oligonucleotides. We also report that anti-miR-33 therapy enhances the expression of nuclear transcription Y subunit gamma (NFYC), a transcriptional regulator required for DNA binding and full transcriptional activation of SREBP-responsive genes, including ACC and FAS. Taken together, these results suggest that persistent inhibition of miR-33 when mice are fed a high-fat diet (HFD) might cause deleterious effects such as moderate hepatic steatosis and hypertriglyceridemia. These unexpected findings highlight the importance of assessing the effect of chronic inhibition of miR-33 in non-human primates before we can translate this therapy to humans.

Systems genetics identifies Sestrin 3 as a regulator of a proconvulsant gene network in human epileptic hippocampus

Gene-regulatory network analysis is a powerful approach to elucidate the molecular processes and pathways underlying complex disease. Here we employ systems genetics approaches to characterize the genetic regulation of pathophysiological pathways in human temporal lobe epilepsy (TLE). Using surgically acquired hippocampi from 129 TLE patients, we identify a gene-regulatory network genetically associated with epilepsy that contains a specialized, highly expressed transcriptional module encoding proconvulsive cytokines and Toll-like receptor signalling genes. RNA sequencing analysis in a mouse model of TLE using 100 epileptic and 100 control hippocampi shows the proconvulsive module is preserved across-species, specific to the epileptic hippocampus and upregulated in chronic epilepsy. In the TLE patients, we map the trans-acting genetic control of this proconvulsive module to Sestrin 3 (SESN3), and demonstrate that SESN3 positively regulates the module in macrophages, microglia and neurons. Morpholino-mediated Sesn3 knockdown in zebrafish confirms the regulation of the transcriptional module, and attenuates chemically induced behavioural seizures in vivo.

Genetic analysis of the cardiac methylome at single nucleotide resolution in a model of human cardiovascular disease

Epigenetic marks such as cytosine methylation are important determinants of cellular and whole-body phenotypes. However, the extent of, and reasons for inter-individual differences in cytosine methylation, and their association with phenotypic variation are poorly characterised. Here we present the first genome-wide study of cytosine methylation at single-nucleotide resolution in an animal model of human disease. We used whole-genome bisulfite sequencing in the spontaneously hypertensive rat (SHR), a model of cardiovascular disease, and the Brown Norway (BN) control strain, to define the genetic architecture of cytosine methylation in the mammalian heart and to test for association between methylation and pathophysiological phenotypes. Analysis of 10.6 million CpG dinucleotides identified 77,088 CpGs that were differentially methylated between the strains. In F1 hybrids we found 38,152 CpGs showing allele-specific methylation and 145 regions with parent-of-origin effects on methylation. Cis-linkage explained almost 60% of inter-strain variation in methylation at a subset of loci tested for linkage in a panel of recombinant inbred (RI) strains. Methylation analysis in isolated cardiomyocytes showed that in the majority of cases methylation differences in cardiomyocytes and non-cardiomyocytes were strain-dependent, confirming a strong genetic component for cytosine methylation. We observed preferential nucleotide usage associated with increased and decreased methylation that is remarkably conserved across species, suggesting a common mechanism for germline control of inter-individual variation in CpG methylation. In the RI strain panel, we found significant correlation of CpG methylation and levels of serum chromogranin B (CgB), a proposed biomarker of heart failure, which is evidence for a link between germline DNA sequence variation, CpG methylation differences and pathophysiological phenotypes in the SHR strain. Together, these results will stimulate further investigation of the molecular basis of locally regulated variation in CpG methylation and provide a starting point for understanding the relationship between the genetic control of CpG methylation and disease phenotypes.

Epigenetic marks provide information that is not encoded in the primary DNA sequence itself but in modifications of genomic DNA and of the associated proteins. Methylation of genomic DNA at cytosine residues is an important epigenetic modification that is associated with developmental processes, carcinogenesis and other diseases. Genome-wide extent of, and reasons for inter-individual differences in cytosine methylation, and their association with phenotypic variation are poorly characterised. To address these questions we have determined and compared the genome-wide methylation patterns in heart tissue of two inbred rat strains, the spontaneously hypertensive rat, an animal model of human disease and a control rat strain. Comparison of methylation differences between genetically identical animals from the same strain and differences between animals from different strains allowed us to quantify association of epigenetic and genetic differences. We show that differences in an individual's germline DNA sequence are important determinants of the variability in methylation between individuals. Comparison with previous reports implicates common mechanisms for regulation of cytosine methylation that are highly conserved across species. Finally, we find correlation between a proposed blood biomarker for heart failure and variation in DNA methylation, suggesting a link between germline DNA sequence variation, methylation and a disease-related phenotype.

Role of miR-195 in aortic aneurysmal disease

RATIONALE

Abdominal aortic aneurysms constitute a degenerative process in the aortic wall. Both the miR-29 and miR-15 families have been implicated in regulating the vascular extracellular matrix.

OBJECTIVE

Our aim was to assess the effect of the miR-15 family on aortic aneurysm development.

METHODS AND RESULTS

Among the miR-15 family members, miR-195 was differentially expressed in aortas of apolipoprotein E-deficient mice on angiotensin II infusion. Proteomics analysis of the secretome of murine aortic smooth muscle cells, after miR-195 manipulation, revealed that miR-195 targets a cadre of extracellular matrix proteins, including collagens, proteoglycans, elastin, and proteins associated with elastic microfibrils, albeit miR-29b showed a stronger effect, particularly in regulating collagens. Systemic and local administration of cholesterol-conjugated antagomiRs revealed better inhibition of miR-195 compared with miR-29b in the uninjured aorta. However, in apolipoprotein E-deficient mice receiving angiotensin II, silencing of miR-29b, but not miR-195, led to an attenuation of aortic dilation. Higher aortic elastin expression was accompanied by an increase of matrix metalloproteinases 2 and 9 in mice treated with antagomiR-195. In human plasma, an inverse correlation of miR-195 was observed with the presence of abdominal aortic aneurysms and aortic diameter.

CONCLUSIONS

We provide the first evidence that miR-195 may contribute to the pathogenesis of aortic aneurysmal disease. Although inhibition of miR-29b proved more effective in preventing aneurysm formation in a preclinical model, miR-195 represents a potent regulator of the aortic extracellular matrix. Notably, plasma levels of miR-195 were reduced in patients with abdominal aortic aneurysms suggesting that microRNAs might serve as a noninvasive biomarker of abdominal aortic aneurysms.

Redox state of pentraxin 3 as a novel biomarker for resolution of inflammation and survival in sepsis

In an endotoxaemic mouse model of sepsis, a tissue-based proteomics approach for biomarker discovery identified long pentraxin 3 (PTX3) as the lead candidate for inflamed myocardium. When the redox-sensitive oligomerization state of PTX3 was further investigated, PTX3 accumulated as an octamer as a result of disulfide-bond formation in heart, kidney, and lung—common organ dysfunctions seen in patients with sepsis. Oligomeric moieties of PTX3 were also detectable in circulation. The oligomerization state of PTX3 was quantified over the first 11 days in critically ill adult patients with sepsis. On admission day, there was no difference in the oligomerization state of PTX3 between survivors and non-survivors. From day 2 onward, the conversion of octameric to monomeric PTX3 was consistently associated with a greater survival after 28 days of follow-up. For example, by day 2 post-admission, octameric PTX3 was barely detectable in survivors, but it still constituted more than half of the total PTX3 in non-survivors (p < 0.001). Monomeric PTX3 was inversely associated with cardiac damage markers NT-proBNP and high-sensitivity troponin I and T. Relative to the conventional measurements of total PTX3 or NT-proBNP, the oligomerization of PTX3 was a superior predictor of disease outcome.

Lipidomics Profiling and Risk of Cardiovascular Disease in the Prospective Population-Based Bruneck Study

Background—

The bulk of cardiovascular disease risk is not explained by traditional risk factors. Recent advances in mass spectrometry allow the identification and quantification of hundreds of lipid species. Molecular lipid profiling by mass spectrometry may improve cardiovascular risk prediction.

Methods and Results—

Lipids were extracted from 685 plasma samples of the prospective population-based Bruneck Study (baseline evaluation in 2000). One hundred thirty-five lipid species from 8 different lipid classes were profiled by shotgun lipidomics with the use of a triple-quadrupole mass spectrometer. Levels of individual species of cholesterol esters (CEs), lysophosphatidylcholines, phosphatidylcholines, phosphatidylethanolamines (PEs), sphingomyelins, and triacylglycerols (TAGs) were associated with cardiovascular disease over a 10-year observation period (2000–2010, 90 incident events). Among the lipid species with the strongest predictive value were TAGs and CEs with a low carbon number and double-bond content, including TAG(54:2) and CE(16:1), as well as PE(36:5) ( P =5.1×10 −7 , 2.2×10 −4 , and 2.5×10 −3 , respectively). Consideration of these 3 lipid species on top of traditional risk factors resulted in improved risk discrimination and classification for cardiovascular disease (cross-validated ΔC index, 0.0210 [95% confidence interval, 0.0010-0.0422]; integrated discrimination improvement, 0.0212 [95% confidence interval, 0.0031-0.0406]; and continuous net reclassification index, 0.398 [95% confidence interval, 0.175-0.619]). A similar shift in the plasma fatty acid composition was associated with cardiovascular disease in the UK Twin Registry (n=1453, 45 cases).

Conclusions—

This study applied mass spectrometry-based lipidomics profiling to population-based cohorts and identified molecular lipid signatures for cardiovascular disease. Molecular lipid species constitute promising new biomarkers that outperform the conventional biochemical measurements of lipid classes currently used in clinics.

Extracellular matrix secretion by cardiac fibroblasts: role of microRNA-29b and microRNA-30c

RATIONALE

MicroRNAs (miRNAs), in particular miR-29b and miR-30c, have been implicated as important regulators of cardiac fibrosis.

OBJECTIVE

To perform a proteomics comparison of miRNA effects on extracellular matrix secretion by cardiac fibroblasts.

METHODS AND RESULTS

Mouse cardiac fibroblasts were transfected with pre-/anti-miR of miR-29b and miR-30c, and their conditioned medium was analyzed by mass spectrometry. miR-29b targeted a cadre of proteins involved in fibrosis, including multiple collagens, matrix metalloproteinases, and leukemia inhibitory factor, insulin-like growth factor 1, and pentraxin 3, 3 predicted targets of miR-29b. miR-29b also attenuated the cardiac fibroblast response to transforming growth factor-β. In contrast, miR-30c had little effect on extracellular matrix production but opposite effects regarding leukemia inhibitory factor and insulin-like growth factor 1. Both miRNAs indirectly affected cardiac myocytes. On transfection with pre-miR-29b, the conditioned medium of cardiac fibroblasts lost its ability to support adhesion of rat ventricular myocytes and led to a significant reduction of cardiac myocyte proteins (α-actinin, cardiac myosin-binding protein C, and cardiac troponin I). Similarly, cardiomyocytes derived from mouse embryonic stem cells atrophied under pre-miR-29 conditioned medium, whereas pre-miR-30c conditioned medium had a prohypertrophic effect. Levels of miR-29a, miR-29c, and miR-30c, but not miR-29b, were significantly reduced in a mouse model of pathological but not physiological hypertrophy. Treatment with antagomiRs to miR-29b induced excess fibrosis after aortic constriction without overt deterioration in cardiac function.

CONCLUSIONS

Our proteomic analysis revealed novel molecular targets of miRNAs that are linked to a fibrogenic cardiac phenotype. Such comprehensive screening methods are essential to define the concerted actions of miRNAs in cardiovascular disease.

Proteomic Identification of Matrix Metalloproteinase Substrates in the Human Vasculature

Background—

Matrix metalloproteinases (MMPs) play a key role in cardiovascular disease, in particular aneurysm formation and plaque rupture. Surprisingly, little is known about MMP substrates in the vasculature.

Methods and Results—

We used a proteomics approach to identify vascular substrates for 3 MMPs, 1 of each of the 3 major classes of MMPs: Human arteries were incubated with MMP-3 (a member of stromelysins), MMP-9 (considered a gelatinase), and MMP-14 (considered a member of the collagenases and of the membrane-bound MMPs). Candidate substrates were identified by mass spectrometry based on increased release from the arterial tissue on digestion, spectral evidence for proteolytic degradation after gel separation, and identification of nontryptic cleavage sites. Using this approach, novel candidates were identified, including extracellular matrix proteins associated with the basement membrane, elastic fibers (emilin-1), and other extracellular proteins (periostin, tenascin-X). Seventy-four nontryptic cleavage sites were detected, many of which were shared among different MMPs. The proteomics findings were validated by immunoblotting and by digesting recombinant/purified proteins with exogenous MMPs. As proof-of-principle, results were related to in vivo pathology by searching for corresponding degradation products in human aortic tissue with different levels of endogenous MMP-9.

Conclusions—

The application of proteomics to identify MMP targets is a new frontier in cardiovascular research. Our current classification of MMPs based on few substrates is an oversimplification of a complex area of biology. This study provides a more comprehensive assessment of potential MMP substrates in the vasculature and represents a valuable resource for future investigations.

Circulating microRNAs as novel biomarkers for platelet activation

RATIONALE

MicroRNA (miRNA) biomarkers are attracting considerable interest. Effects of medication, however, have not been investigated thus far.

OBJECTIVE

To analyze changes in plasma miRNAs in response to antiplatelet therapy.

METHODS AND RESULTS

Profiling for 377 miRNAs was performed in platelets, platelet microparticles, platelet-rich plasma, platelet-poor plasma, and serum. Platelet-rich plasma showed markedly higher levels of miRNAs than serum and platelet-poor plasma. Few abundant platelet miRNAs, such as miR-24, miR-197, miR-191, and miR-223, were also increased in serum compared with platelet-poor plasma. In contrast, antiplatelet therapy significantly reduced miRNA levels. Using custom-made quantitative real-time polymerase chain reaction plates, 92 miRNAs were assessed in a dose-escalation study in healthy volunteers at 4 different time points: at baseline without therapy, at 1 week with 10 mg prasugrel, at 2 weeks with 10 mg prasugrel plus 75 mg aspirin, and at 3 weeks with 10 mg prasugrel plus 300 mg aspirin. Findings in healthy volunteers were confirmed by individual TaqMan quantitative real-time polymerase chain reaction assays (n=9). Validation was performed in an independent cohort of patients with symptomatic atherosclerosis (n=33), who received low-dose aspirin at baseline. Plasma levels of platelet miRNAs, such as miR-223, miR-191, and others, that is, miR-126 and miR-150, decreased on further platelet inhibition.

CONCLUSIONS

Our study demonstrated a substantial platelet contribution to the circulating miRNA pool and identified miRNAs responsive to antiplatelet therapy. It also highlights that antiplatelet therapy and preparation of blood samples could be confounding factors in case-control studies relating plasma miRNAs to cardiovascular disease.

Effects of perhexiline-induced fuel switch on the cardiac proteome and metabolome

Perhexiline is a potent anti-anginal drug used for treatment of refractory angina and other forms of heart disease. It provides an oxygen sparing effect in the myocardium by creating a switch from fatty acid to glucose metabolism through partial inhibition of carnitine palmitoyltransferase 1 and 2. However, the precise molecular mechanisms underlying the cardioprotective effects elicited by perhexiline are not fully understood. The present study employed a combined proteomics, metabolomics and computational approach to characterise changes in murine hearts upon treatment with perhexiline. According to results based on difference in-gel electrophoresis, the most profound change in the cardiac proteome related to the activation of the pyruvate dehydrogenase complex. Metabolomic analysis by high-resolution nuclear magnetic resonance spectroscopy showed lower levels of total creatine and taurine in hearts of perhexiline-treated mice. Creatine and taurine levels were also significantly correlated in a cross-correlation analysis of all metabolites. Computational modelling suggested that far from inducing a simple shift from fatty acid to glucose oxidation, perhexiline may cause complex rebalancing of carbon and nucleotide phosphate fluxes, fuelled by increased lactate and amino acid uptake, to increase metabolic flexibility and to maintain cardiac output. This article is part of a Special Issue entitled "Focus on Cardiac Metabolism".

Proteomics: from single molecules to biological pathways

The conventional reductionist approach to cardiovascular research investigates individual candidate factors or linear signalling pathways but ignores more complex interactions in biological systems. The advent of molecular profiling technologies that focus on a global characterization of whole complements allows an exploration of the interconnectivity of pathways during pathophysiologically relevant processes, but has brought about the issue of statistical analysis and data integration. Proteins identified by differential expression as well as those in protein–protein interaction networks identified through experiments and through computational modelling techniques can be used as an initial starting point for functional analyses. In combination with other ‘-omics’ technologies, such as transcriptomics and metabolomics, proteomics explores different aspects of disease, and the different pillars of observations facilitate the data integration in disease-specific networks. Ultimately, a systems biology approach may advance our understanding of cardiovascular disease processes at a ‘biological pathway’ instead of a ‘single molecule’ level and accelerate progress towards disease-modifying interventions.

Systems-level approaches reveal conservation of trans-regulated genes in the rat and genetic determinants of blood pressure in humans

AIMS

Human genome-wide association studies (GWAS) of hypertension identified only few susceptibility loci with large effect that were replicated across populations. The vast majority of genes detected by GWAS has small effect and the regulatory mechanisms through which these genetic variants cause disease remain mostly unclear. Here, we used comparative genomics between human and an established rat model of hypertension to explore the transcriptional mechanisms mediating the effect of genes identified in 15 hypertension GWAS.

METHODS AND RESULTS

Time series analysis of radiotelemetric blood pressure (BP) was performed to assess 11 parameters of BP variation in recombinant inbred strains derived from the spontaneously hypertensive rat. BP data were integrated with ∼27 000 expression quantative trait loci (eQTLs) mapped across seven tissues, detecting >8000 significant associations between eQTL genes and BP variation in the rat. We then compiled a large catalogue of human genes from GWAS of hypertension and identified a subset of 2292 rat-human orthologous genes. Expression levels for 795 (34%) of these genes correlated with BP variation across rat tissues: 51 genes were cis-regulated, whereas 459 were trans-regulated and enriched for 'calcium signalling pathway' (P = 9.6 × 10(-6)) and 'ion channel' genes (P = 3.5 × 10(-7)), which are important determinants of hypertension. We identified 158 clusters of trans-eQTLs, annotated the underlying 'master regulator' genes and found significant over-representation in the human hypertension gene set (enrichment P = 5 × 10(-4)).

CONCLUSION

We showed extensive conservation of trans-regulated genes and their master regulators between rat and human hypertension. These findings reveal that small-effect genes associated with hypertension by human GWAS are likely to exert their action through coordinate regulation of pathogenic pathways.

Determining signalling nodes for apoptosis by a genetic high-throughput screen

BACKGROUND

With the ever-increasing information emerging from the various sequencing and gene annotation projects, there is an urgent need to elucidate the cellular functions of the newly discovered genes. The genetically regulated cell suicide of apoptosis is especially suitable for such endeavours as it is governed by a vast number of factors.

METHODOLOGY/PRINCIPAL FINDINGS

We have set up a high-throughput screen in 96-well microtiter plates for genes that induce apoptosis upon their individual transfection into human cells. Upon screening approximately 100,000 cDNA clones we determined 74 genes that initiate this cellular suicide programme. A thorough bioinformatics analysis of these genes revealed that 91% are novel apoptosis regulators. Careful sequence analysis and functional annotation showed that the apoptosis factors exhibit a distinct functional distribution that distinguishes the cell death process from other signalling pathways. While only a minority of classic signal transducers were determined, a substantial number of the genes fall into the transporter- and enzyme-category. The apoptosis factors are distributed throughout all cellular organelles and many signalling circuits, but one distinct signalling pathway connects at least some of the isolated genes. Comparisons with microarray data suggest that several genes are dysregulated in specific types of cancers and degenerative diseases.

CONCLUSIONS/SIGNIFICANCE

Many unknown genes for cell death were revealed through our screen, supporting the enormous complexity of cell death regulation. Our results will serve as a repository for other researchers working with genomics data related to apoptosis or for those seeking to reveal novel signalling pathways for cell suicide.

ESS++: a C++ objected-oriented algorithm for Bayesian stochastic search model exploration

SUMMARY

ESS++ is a C++ implementation of a fully Bayesian variable selection approach for single and multiple response linear regression. ESS++ works well both when the number of observations is larger than the number of predictors and in the 'large p, small n' case. In the current version, ESS++ can handle several hundred observations, thousands of predictors and a few responses simultaneously. The core engine of ESS++ for the selection of relevant predictors is based on Evolutionary Monte Carlo. Our implementation is open source, allowing community-based alterations and improvements.

AVAILABILITY

C++ source code and documentation including compilation instructions are available under GNU licence at http://bgx.org.uk/software/ESS.html.

New insights into the genetic control of gene expression using a Bayesian multi-tissue approach

The majority of expression quantitative trait locus (eQTL) studies have been carried out in single tissues or cell types, using methods that ignore information shared across tissues. Although global analysis of RNA expression in multiple tissues is now feasible, few integrated statistical frameworks for joint analysis of gene expression across tissues combined with simultaneous analysis of multiple genetic variants have been developed to date. Here, we propose Sparse Bayesian Regression models for mapping eQTLs within individual tissues and simultaneously across tissues. Testing these on a set of 2,000 genes in four tissues, we demonstrate that our methods are more powerful than traditional approaches in revealing the true complexity of the eQTL landscape at the systems-level. Highlighting the power of our method, we identified a two-eQTL model (cis/trans) for the Hopx gene that was experimentally validated and was not detected by conventional approaches. We showed common genetic regulation of gene expression across four tissues for ∼27% of transcripts, providing >5 fold increase in eQTLs detection when compared with single tissue analyses at 5% FDR level. These findings provide a new opportunity to uncover complex genetic regulatory mechanisms controlling global gene expression while the generality of our modelling approach makes it adaptable to other model systems and humans, with broad application to analysis of multiple intermediate and whole-body phenotypes.

Integrated analysis of genome-wide genetic polymorphisms and gene expression profiles from different tissues or cell types has been highly successful in identifying genes modulating complex phenotypes in animal models and humans. However, an important limitation of the current approaches consists in their sole application to individual tissues, thus ignoring information shared across different tissues. To uncover complex genetic regulatory mechanisms controlling gene expression at the whole organism's level, it is essential to develop appropriate analytical methods for the analysis of genome-wide genetic polymorphisms and gene expression profiles simultaneously in multiple tissues. This paper presents a novel, fully integrated Bayesian approach for mapping the genetic components of gene expression within and across multiple tissues. In addition to increased power and enhanced mapping resolution when compared with traditional approaches, our model directly provides information on potential systemic effects on transcriptional profiles and co-existing local (cis) and distant (trans) genetic control of gene expression. We also discuss the possibility to extend our approach for the analysis of different phenotypes, and other study designs, thus providing an integrated computational tool to explore the genetic control underlying transcriptional regulation at the systems-level, beyond the single tissue resolution.

Differences in DBA/1J and DBA/2J reveal lipid QTL genes

Recent advances in mouse genomics have revealed considerable variation in the form of single-nucleotide polymorphisms (SNPs) among common inbred strains. This has made it possible to characterize closely related strains and to identify genes that differ; such genes may be causal for quantitative phenotypes. The mouse strains DBA/1J and DBA/2J differ by just 5.6% at the SNP level. These strains exhibit differences in a number of metabolic and lipid phenotypes, such as plasma levels of triglycerides (TGs) and HDL. A cross between these strains revealed multiple quantitative trait loci (QTLs) in 294 progeny. We identified significant TG QTLs on chromosomes (Chrs) 1, 2, 3, 4, 8, 9, 10, 11, 12, 13, 14, 16, and 19, and significant HDL QTLs on Chrs 3, 9, and 16. Some QTLs mapped to chromosomes with limited variability between the two strains, thus facilitating the identification of candidate genes. We suggest that Tshr is the QTL gene for Chr 12 TG and HDL levels and that Ihh may account for the TG QTL on Chr 1. This cross highlights the advantage of crossing closely related strains for subsequent identification of QTL genes.