A HaemAtlas: characterizing gene expression in differentiated human blood cells

Measuring cell identity in noisy biological systems

Global gene expression measurements are increasingly obtained as a function of cell type, spatial position within a tissue and other biologically meaningful coordinates. Such data should enable quantitative analysis of the cell-type specificity of gene expression, but such analyses can often be confounded by the presence of noise. We introduce a specificity measure Spec that quantifies the information in a gene's complete expression profile regarding any given cell type, and an uncertainty measure dSpec, which measures the effect of noise on specificity. Using global gene expression data from the mouse brain, plant root and human white blood cells, we show that Spec identifies genes with variable expression levels that are nonetheless highly specific of particular cell types. When samples from different individuals are used, dSpec measures genes’ transcriptional plasticity in each cell type. Our approach is broadly applicable to mapped gene expression measurements in stem cell biology, developmental biology, cancer biology and biomarker identification. As an example of such applications, we show that Spec identifies a new class of biomarkers, which exhibit variable expression without compromising specificity. The approach provides a unifying theoretical framework for quantifying specificity in the presence of noise, which is widely applicable across diverse biological systems.

Exome sequencing identifies NBEAL2 as the causative gene for gray platelet syndrome

Gray platelet syndrome (GPS) is a predominantly recessive platelet disorder that is characterized by mild thrombocytopenia with large platelets and a paucity of α-granules; these abnormalities cause mostly moderate but in rare cases severe bleeding. We sequenced the exomes of four unrelated individuals and identified NBEAL2 as the causative gene; it has no previously known function but is a member of a gene family that is involved in granule development. Silencing of nbeal2 in zebrafish abrogated thrombocyte formation.

Protein kinase C mediates platelet secretion and thrombus formation through protein kinase D2

Platelets are highly specialized blood cells critically involved in hemostasis and thrombosis. Members of the protein kinase C (PKC) family have established roles in regulating platelet function and thrombosis, but the molecular mechanisms are not clearly understood. In particular, the conventional PKC isoform, PKCα, is a major regulator of platelet granule secretion, but the molecular pathway from PKCα to secretion is not defined. Protein kinase D (PKD) is a family of 3 kinases activated by PKC, which may represent a step in the PKC signaling pathway to secretion. In the present study, we show that PKD2 is the sole PKD member regulated downstream of PKC in platelets, and that the conventional, but not novel, PKC isoforms provide the upstream signal. Platelets from a gene knock-in mouse in which 2 key phosphorylation sites in PKD2 have been mutated (Ser707Ala/Ser711Ala) show a significant reduction in agonist-induced dense granule secretion, but not in α-granule secretion. This deficiency in dense granule release was responsible for a reduced platelet aggregation and a marked reduction in thrombus formation. Our results show that in the molecular pathway to secretion, PKD2 is a key component of the PKC-mediated pathway to platelet activation and thrombus formation through its selective regulation of dense granule secretion.

The genetics of common variation affecting platelet development, function and pharmaceutical targeting

Common variant effects on human platelet function and response to anti-platelet treatment have traditionally been studied using candidate gene approaches involving a limited number of variants and genes. These studies have often been undertaken in clinically defined cohorts. More recently, studies have applied genome-wide scans in larger population samples than prior candidate studies, in some cases scanning relatively healthy individuals. These studies demonstrate synergy with some prior candidate gene findings (e.g., GP6, ADRA2A) but also uncover novel loci involved in platelet function. Here, I summarise findings on common genetic variation influencing platelet development, function and therapeutics. Taken together, candidate gene and genome-wide studies begin to account for common variation in platelet function and provide information that may ultimately be useful in pharmacogenetic applications in the clinic. More than 50 loci have been identified with consistent associations with platelet phenotypes in ≥ 2 populations. Several variants are under further study in clinical trials relating to anti-platelet therapies. In order to have useful clinical applications, variants must have large effects on a modifiable outcome. Regardless of clinical applications, studies of common genetic influences, even of small effect, offer additional insights into platelet biology including the importance of intracellular signalling and novel receptors. Understanding of common platelet-related genetics remains behind parallel fields (e.g., lipids, blood pressure) due to challenges in phenotype ascertainment. Further work is necessary to discover and characterise loci for platelet function, and to assess whether these loci contribute to disease aetiologies or response to therapeutics.

CD248+ stromal cells are associated with progressive chronic kidney disease

Stromal fibroblasts are the primary cells of the kidney that produce fibrotic matrix. CD248 is a stromal marker expressed on fibroblasts and pericytes within the human kidney. Here, we tested whether CD248 expression in the kidney colocalizes with fibrosis and if it is associated with known determinants of chronic kidney disease (CKD). CD248 expression was located and quantified in situ by immunohistochemistry in kidney biopsies from 93 patients with IgA nephropathy and compared with 22 archived biopsies encompassing normal kidney tissue as control. In normal kidney tissue, CD248 was expressed by resident pericytes, stromal fibroblasts, and was upregulated in human CKD. The expression was linked to known determinants of renal progression. This relationship was maintained in a multivariate analysis with CD248 expression linked to renal survival. CD248 was expressed by a population of α-smooth muscle actin (SMA)(+) myofibroblasts and α-SMA(-) stromal cells but not expressed on CD45(+) leukocytes. Thus, CD248 defines a subset of stromal cells, including but not limited to some myofibroblasts, linked to albuminuria and tubulointerstitial damage during tissue remodeling in CKD.

The stromal cell antigen CD248 (endosialin) is expressed on naive CD8+ human T cells and regulates proliferation

CD248 (endosialin) is a transmembrane glycoprotein that is dynamically expressed on pericytes and fibroblasts during tissue development, tumour neovascularization and inflammation. Its role in tissue remodelling is associated with increased stromal cell proliferation and migration. We show that CD248 is also uniquely expressed by human, but not mouse (C57BL/6), CD8(+) naive T cells. CD248 is found only on CD8(+) CCR7(+) CD11a(low) naive T cells and on CD8 single-positive T cells in the thymus. Transfection of the CD248 negative T-cell line MOLT-4 with CD248 cDNA surprisingly reduced cell proliferation. Knock-down of CD248 on naive CD8 T cells increased cell proliferation. These data demonstrate opposing functions for CD248 on haematopoietic (CD8(+)) versus stromal cells and suggests that CD248 helps to maintain naive CD8(+) human T cells in a quiescent state.

Regulating thrombus growth and stability to achieve an optimal response to injury

An optimal platelet response to injury can be defined as one in which blood loss is restrained and haemostasis is achieved without the penalty of further tissue damage caused by unwarranted vascular occlusion. This brief review considers some of the ways in which thrombus growth and stability can be regulated so that an optimal platelet response can be achieved in vivo. Three related topics are considered. The first focuses on intracellular mechanisms that regulate the early events of platelet activation downstream of G protein coupled receptors for agonists such as thrombin, thromboxane A(2) and ADP. The second considers the ways in which signalling events that are dependent on stable contacts between platelets can influence the state of platelet activation and thus affect thrombus growth and stability. The third focuses on the changes that are experienced by platelets as they move from their normal environment in freely-flowing plasma to a very different environment within the growing haemostatic plug, an environment in which the narrowing gaps and junctions between platelets not only facilitate communication, but also increasingly limit both the penetration of plasma and the exodus of platelet-derived bioactive molecules.

BCL2/BCL-X(L) inhibition induces apoptosis, disrupts cellular calcium homeostasis, and prevents platelet activation

Apoptosis in megakaryocytes results in the formation of platelets. The role of apoptotic pathways in platelet turnover and in the apoptotic-like changes seen after platelet activation is poorly understood. ABT-263 (Navitoclax), a specific inhibitor of antiapoptotic BCL2 proteins, which is currently being evaluated in clinical trials for the treatment of leukemia and other malignancies, induces a dose-limiting thrombocytopenia. In this study, the relationship between BCL2/BCL-X(L) inhibition, apoptosis, and platelet activation was investigated. Exposure to ABT-263 induced apoptosis but repressed platelet activation by physiologic agonists. Notably, ABT-263 induced an immediate calcium response in platelets and the depletion of intracellular calcium stores, indicating that on BCL2/BCL-X(L) inhibition platelet activation is abrogated because of a diminished calcium signaling. By comparing the effects of ABT-263 and its analog ABT-737 on platelets and leukemia cells from the same donor, we show, for the first time, that these BCL2/BCL-X(L) inhibitors do not offer any selective toxicity but induce apoptosis at similar concentrations in leukemia cells and platelets. However, reticulated platelets are less sensitive to apoptosis, supporting the hypothesis that treatment with ABT-263 induces a selective loss of older platelets and providing an explanation for the transient thrombocytopenia observed on ABT-263 treatment.

Cell subset prediction for blood genomic studies

BACKGROUND

Genome-wide transcriptional profiling of patient blood samples offers a powerful tool to investigate underlying disease mechanisms and personalized treatment decisions. Most studies are based on analysis of total peripheral blood mononuclear cells (PBMCs), a mixed population. In this case, accuracy is inherently limited since cell subset-specific differential expression of gene signatures will be diluted by RNA from other cells. While using specific PBMC subsets for transcriptional profiling would improve our ability to extract knowledge from these data, it is rarely obvious which cell subset(s) will be the most informative.

RESULTS

We have developed a computational method (Subset Prediction from Enrichment Correlation, SPEC) to predict the cellular source for a pre-defined list of genes (i.e. a gene signature) using only data from total PBMCs. SPEC does not rely on the occurrence of cell subset-specific genes in the signature, but rather takes advantage of correlations with subset-specific genes across a set of samples. Validation using multiple experimental datasets demonstrates that SPEC can accurately identify the source of a gene signature as myeloid or lymphoid, as well as differentiate between B cells, T cells, NK cells and monocytes. Using SPEC, we predict that myeloid cells are the source of the interferon-therapy response gene signature associated with HCV patients who are non-responsive to standard therapy.

CONCLUSIONS

SPEC is a powerful technique for blood genomic studies. It can help identify specific cell subsets that are important for understanding disease and therapy response. SPEC is widely applicable since only gene expression profiles from total PBMCs are required, and thus it can easily be used to mine the massive amount of existing microarray or RNA-seq data.

Maps of open chromatin guide the functional follow-up of genome-wide association signals: application to hematological traits

Turning genetic discoveries identified in genome-wide association (GWA) studies into biological mechanisms is an important challenge in human genetics. Many GWA signals map outside exons, suggesting that the associated variants may lie within regulatory regions. We applied the formaldehyde-assisted isolation of regulatory elements (FAIRE) method in a megakaryocytic and an erythroblastoid cell line to map active regulatory elements at known loci associated with hematological quantitative traits, coronary artery disease, and myocardial infarction. We showed that the two cell types exhibit distinct patterns of open chromatin and that cell-specific open chromatin can guide the finding of functional variants. We identified an open chromatin region at chromosome 7q22.3 in megakaryocytes but not erythroblasts, which harbors the common non-coding sequence variant rs342293 known to be associated with platelet volume and function. Resequencing of this open chromatin region in 643 individuals provided strong evidence that rs342293 is the only putative causative variant in this region. We demonstrated that the C- and G-alleles differentially bind the transcription factor EVI1 affecting PIK3CG gene expression in platelets and macrophages. A protein-protein interaction network including up- and down-regulated genes in Pik3cg knockout mice indicated that PIK3CG is associated with gene pathways with an established role in platelet membrane biogenesis and thrombus formation. Thus, rs342293 is the functional common variant at this locus; to the best of our knowledge this is the first such variant to be elucidated among the known platelet quantitative trait loci (QTLs). Our data suggested a molecular mechanism by which a non-coding GWA index SNP modulates platelet phenotype.

Arginine transport in human monocytic leukemia THP-1 cells during macrophage differentiation

L-arginine metabolism in myeloid cells plays a central role in the processes of macrophage activation and in the regulation of immune responses. In this study, we investigated arginine transport activity and the expression of the related transporter genes during the differentiation of monocytes to macrophages. We show here that the induction of THP-1 monocyte differentiation by PMA markedly increases the expression of SLC7A7 mRNA and of y(+)LAT1 protein and consequently, the activity of system y(+)L-mediated arginine transport. Conversely, the activity of system y(+) decreases during macrophage differentiation as a result of a reduction in CAT1 protein expression. The PMA-induced, macrophage-differentiated phenotype and the increased activity of system y(+)L through the induction of SLC7A7 gene are mediated by the specific activation of PKCβ. SLC7A7 gene silencing causes a significant reduction of system y(+)L activity and a subsequent, marked increase of arginine and lysine cell content, thus suggesting that in macrophagic cells, system y(+)L activity is mainly directed outwardly. Differentiating agents other than PMA, i.e., VD3 and ATRA, are equally effective in the stimulation of system y(+)L transport activity through the increased expression of SLC7A7 mRNA and y(+)LAT1 protein. Moreover, we found that also during differentiation of human monocytes from peripheral blood SLC7A7 mRNA and system y(+)L activity are increased. These findings point to SLC7A7 gene as a marker of macrophage differentiation.

Genome-wide analysis of simultaneous GATA1/2, RUNX1, FLI1, and SCL binding in megakaryocytes identifies hematopoietic regulators

Hematopoietic differentiation critically depends on combinations of transcriptional regulators controlling the development of individual lineages. Here, we report the genome-wide binding sites for the five key hematopoietic transcription factors--GATA1, GATA2, RUNX1, FLI1, and TAL1/SCL--in primary human megakaryocytes. Statistical analysis of the 17,263 regions bound by at least one factor demonstrated that simultaneous binding by all five factors was the most enriched pattern and often occurred near known hematopoietic regulators. Eight genes not previously appreciated to function in hematopoiesis that were bound by all five factors were shown to be essential for thrombocyte and/or erythroid development in zebrafish. Moreover, one of these genes encoding the PDZK1IP1 protein shared transcriptional enhancer elements with the blood stem cell regulator TAL1/SCL. Multifactor ChIP-Seq analysis in primary human cells coupled with a high-throughput in vivo perturbation screen therefore offers a powerful strategy to identify essential regulators of complex mammalian differentiation processes.

Platelets at work in primary hemostasis

When platelet numbers are low or when their function is disabled, the risk of bleeding is high, which on the one hand indicates that in normal life vascular damage is a rather common event and that hence the role of platelets in maintaining a normal hemostasis is a continuously ongoing physiological process. Upon vascular injury, platelets instantly adhere to the exposed extracellular matrix resulting in platelet activation and aggregation to form a hemostatic plug. This self-amplifying mechanism nevertheless requires a tight control to prevent uncontrolled platelet aggregate formation that eventually would occlude the vessel. Therefore endothelial cells produce inhibitory compounds such as prostacyclin and nitric oxide that limit the growth of the platelet thrombus to the damaged area. With this review, we intend to give an integrated survey of the platelet response to vascular injury in normal hemostasis.

Global gene expression analysis of human erythroid progenitors

Understanding the pattern of gene expression during erythropoiesis is crucial for a synthesis of erythroid developmental biology. Here, we isolated 4 distinct populations at successive erythropoietin-dependent stages of erythropoiesis, including the terminal, pyknotic stage. The transcriptome was determined using Affymetrix arrays. First, we demonstrated the importance of using defined cell populations to identify lineage and temporally specific patterns of gene expression. Cells sorted by surface expression profile not only express significantly fewer genes than unsorted cells but also demonstrate significantly greater differences in the expression levels of particular genes between stages than unsorted cells. Second, using standard software, we identified more than 1000 transcripts not previously observed to be differentially expressed during erythroid maturation, 13 of which are highly significantly terminally regulated, including RFXAP and SMARCA4. Third, using matched filtering, we identified 12 transcripts not previously reported to be continuously up-regulated in maturing human primary erythroblasts. Finally, using transcription factor binding site analysis, we identified potential transcription factors that may regulate gene expression during terminal erythropoiesis. Our stringent lists of differentially regulated and continuously expressed transcripts containing many genes with undiscovered functions in erythroblasts are a resource for future functional studies of erythropoiesis. Our Human Erythroid Maturation database is available at https://cellline.molbiol.ox.ac.uk/eryth/index.html. [corrected].

The role of meis1 in primitive and definitive hematopoiesis during zebrafish development

BACKGROUND

The Meis1 protein represents an important cofactor for Hox and Pbx1 and is implicated in human and murine leukemias. Though much is known about the role of meis1 in leukemogenesis, its function in normal hematopoiesis remains largely unclear. Here we characterized the role of the proto-oncogene, meis1, during zebrafish primitive and definitive hematopoiesis.

DESIGN AND METHODS

Zebrafish embryos were stained with o-dianisidine to detect hemoglobin-containing cells and Sudan black to quantify neutrophils. The numbers of other cells (scl-, gata1- and alas2-positive cells) were also quantified by measuring the corresponding stained areas of the embryos. We used anti-Meis1 antibody and whole mount immunohistochemistry to determine the pattern of expression of Meis1 during zebrafish development and then analyzed the functional role of Meis1 by knocking-down the meis1 gene.

RESULTS

Using antisense morpholino oligomers to interrupt meis1 expression we found that, although primitive macrophage development could occur unhampered, posterior erythroid differentiation required meis1, and its absence resulted in a severe decrease in the number of mature erythrocytes. Furthermore a picture emerged that meis1 exerts important effects on later stages of erythrocyte maturation and that these effects are independent of gata1, but under the control of scl. In addition, meis1 morpholino knock-down led to dramatic single arteriovenous tube formation. We also found that knock-down of pbx1 resulted in a phenotype that was strikingly similar to that of meis1 knock-down zebrafish.

CONCLUSIONS

These results imply that meis1, jointly with pbx1, regulates primitive hematopoiesis as well as vascular development.

MicroRNA-15a and -16-1 act via MYB to elevate fetal hemoglobin expression in human trisomy 13

Many human aneuploidy syndromes have unique phenotypic consequences, but in most instances it is unclear whether these phenotypes are attributable to alterations in the dosage of specific genes. In human trisomy 13, there is delayed switching and persistence of fetal hemoglobin (HbF) and elevation of embryonic hemoglobin in newborns. Using partial trisomy cases, we mapped this trait to chromosomal band 13q14; by examining the genes in this region, two microRNAs, miR-15a and -16-1, appear as top candidates for the elevated HbF levels. Indeed, increased expression of these microRNAs in primary human erythroid progenitor cells results in elevated fetal and embryonic hemoglobin gene expression. Moreover, we show that a direct target of these microRNAs, MYB, plays an important role in silencing the fetal and embryonic hemoglobin genes. Thus we demonstrate how the developmental regulation of a clinically important human trait can be better understood through the genetic and functional study of aneuploidy syndromes and suggest that miR-15a, -16-1, and MYB may be important therapeutic targets to increase HbF levels in patients with sickle cell disease and β-thalassemia.

Comparative transcriptional and phenotypic peripheral blood analysis of kidney recipients under cyclosporin A or sirolimus monotherapy

Due to its low level of nephrotoxicity and capacity to harness tolerogenic pathways, sirolimus (SRL) has been proposed as an alternative to calcineurin inhibitors in transplantation. The exact mechanisms underlying its unique immunosuppressive profile in humans, however, are still not well understood. In the current study, we aimed to depict the in vivo effects of SRL in comparison with cyclosporin A (CSA) by employing gene expression profiling and multiparameter flow cytometry on blood cells collected from stable kidney recipients under immunosuppressant monotherapy. SRL recipients displayed an increased frequency of CD4 + CD25highFoxp3 + T cells. However, this was accompanied by an increased number of effector memory T cells and by enrichment in NFkB-related pro-inflammatory expression pathways and monocyte and NK cell lineage-specific transcripts. Furthermore, measurement of a transcriptional signature characteristic of operationally tolerant kidney recipients failed to detect differences between SRL and CSA-treated recipients. In conclusion, we show here that the blood transcriptional profile induced by SRL monotherapy in vivo does not resemble that of operationally tolerant recipients and is dominated by innate immune cells and NFkB-related pro-inflammatory events. These data provide novel insights on the complex effects of SLR on the immune system in clinical transplantation.

Transcription profiling in human platelets reveals LRRFIP1 as a novel protein regulating platelet function

Within the healthy population, there is substantial, heritable, and interindividual variability in the platelet response. We explored whether a proportion of this variability could be accounted for by interindividual variation in gene expression. Through a correlative analysis of genome-wide platelet RNA expression data from 37 subjects representing the normal range of platelet responsiveness within a cohort of 500 subjects, we identified 63 genes in which transcript levels correlated with variation in the platelet response to adenosine diphosphate and/or the collagen-mimetic peptide, cross-linked collagen-related peptide. Many of these encode proteins with no reported function in platelets. An association study of 6 of the 63 genes in 4235 cases and 6379 controls showed a putative association with myocardial infarction for COMMD7 (COMM domain-containing protein 7) and a major deviation from the null hypo thesis for LRRFIP1 [leucine-rich repeat (in FLII) interacting protein 1]. Morpholino-based silencing in Danio rerio identified a modest role for commd7 and a significant effect for lrrfip1 as positive regulators of thrombus formation. Proteomic analysis of human platelet LRRFIP1-interacting proteins indicated that LRRFIP1 functions as a component of the platelet cytoskeleton, where it interacts with the actin-remodeling proteins Flightless-1 and Drebrin. Taken together, these data reveal novel proteins regulating the platelet response.

Detection of human hematopoietic stem cell engraftment in the livers of adult immunodeficient mice by an optimized flow cytometric method

Immunodeficient NOD.Cg-Prkdc(scid) Il2rg(tm1Wjl)/SzJ (NSG) mice are a valuable resource to study human hematopoietic stem cells. Prolonged multilineage hematopoiesis indicates stem cell engraftment and generally is measured by flow cytometry. In this study, we took advantage of the multi-parameter detection afforded by modern flow cytometers to optimize detection of human hematopoiesis in NSG mice. Antigens widely expressed by mouse or human cells were evaluated as markers to distinguish mixtures of these cells to optimize and test the limits of chimerism detection. The bone marrow, spleen, and liver of NSG mice transplanted with human hematopoietic cells were analyzed for evidence of engraftment.Mouse bone marrow cells were best marked for exclusion by staining with a combination of CD45, TER-119, and anti-H-2K(d) monoclonal antibodies, whereas live human cells were most accurately identified by elimination of cell doublets and positive staining for CD59. Human stem cells (CD34(++)CD133(+)CD38(low)) and progenitors were detected in the bone marrow and liver, but not in the spleen. An unusual pattern of myeloid antigen expression was detected in the bone marrow and CD3(+)CD4(+)CD8(+) T-cells were detected in the spleen. We concluded that multicolor flow cytometric analysis that clearly distinguishes mouse and human cells offers accurate detection of human chimerism in NSG mice. Human hematopoiesis can be detected in the bone marrow and liver of NSG mice with T-lymphopoiesis, possibly occurring in the spleen.

Platelet physiology and antiplatelet agents

Apart from the central beneficial role platelets play in hemostasis, they are also involved in atherothrombotic diseases. Here, we review the current knowledge of platelet intracellular signal transduction pathways involved in platelet adhesion, activation, amplification of the activation signal and aggregation, as well as pathways limiting platelet aggregation. A thorough understanding of these pathways allows explanation of the mechanism of action of existing antiplatelet agents, but also helps to identify targets for novel drug development.

SCLD: a Stem Cell Lineage Database for the annotation of cell types and developmental lineages

Stem cell biology has experienced explosive growth over the past decade as researchers attempt to generate therapeutically relevant cell types in the laboratory. Recapitulation of endogenous developmental trajectories is a dominant paradigm in the design of directed differentiation protocols, and attempts to guide stem cell differentiation are often based explicitly on knowledge of in vivo development. Therefore, when designing protocols, stem cell biologists rely heavily upon information including (i) cell type-specific gene expression profiles, (ii) anatomical and developmental relationships between cells and tissues and (iii) signals important for progression from progenitors to target cell types. Here, we present the Stem Cell Lineage Database (SCLD) (http://scld.mcb.uconn.edu) that aims to unify this information into a single resource where users can easily store and access information about cell type gene expression, cell lineage maps and stem cell differentiation protocols for both human and mouse stem cells and endogenous developmental lineages. By establishing the SCLD, we provide scientists with a centralized location to organize access and share data, dispute and resolve contentious relationships between cell types and within lineages, uncover discriminating cell type marker panels and design directed differentiation protocols.

Deciphering the human platelet sheddome

Activated platelets shed surface proteins, potentially modifying platelet function as well as providing a source of bioactive fragments. Previous studies have identified several constituents of the platelet sheddome, but the full extent of shedding is unknown. Here we have taken a global approach, analyzing protein fragments in the supernate of activated platelets using mass spectroscopy and looking for proteins originating from platelet membranes. After removing plasma proteins and microparticles, 1048 proteins were identified, including 69 membrane proteins. Nearly all of the membrane proteins had been detected previously, but only 10 had been shown to be shed in platelets. The remaining 59 are candidates subject to confirmation. Based on spectral counts, protein representation in the sheddome varies considerably. As proof of principle, we validated one of the less frequently detected proteins, semaphorin 7A, which had not previously been identified in platelets. Surface expression, cleavage, and shedding of semaphorin 7A were demonstrated, as was its association with α-granules. Finally, cleavage of semaphorin 7A and 12 other proteins was substantially reduced by an inhibitor of ADAM17, a known sheddase. These results define a subset of membrane proteins as sheddome candidates, forming the basis for further studies examining the impact of ectodomain shedding on platelet function.

T1DBase: update 2011, organization and presentation of large-scale data sets for type 1 diabetes research

T1DBase (http://www.t1dbase.org) is web platform, which supports the type 1 diabetes (T1D) community. It integrates genetic, genomic and expression data relevant to T1D research across mouse, rat and human and presents this to the user as a set of web pages and tools. This update describes the incorporation of new data sets, tools and curation efforts as well as a new website design to simplify site use. New data sets include curated summary data from four genome-wide association studies relevant to T1D, HaemAtlas—a data set and tool to query gene expression levels in haematopoietic cells and a manually curated table of human T1D susceptibility loci, incorporating genetic overlap with other related diseases. These developments will continue to support T1D research and allow easy access to large and complex T1D relevant data sets.

Multicenter validation of the diagnostic accuracy of a blood-based gene expression test for assessing obstructive coronary artery disease in nondiabetic patients

BACKGROUND

Diagnosing obstructive coronary artery disease (CAD) in at-risk patients can be challenging and typically requires both noninvasive imaging methods and coronary angiography, the gold standard. Previous studies have suggested that peripheral blood gene expression can indicate the presence of CAD.

OBJECTIVE

To validate a previously developed 23-gene, expression-based classification test for diagnosis of obstructive CAD in nondiabetic patients.

DESIGN

Multicenter prospective trial with blood samples obtained before coronary angiography. (ClinicalTrials.gov registration number: NCT00500617) SETTING: 39 centers in the United States.

PATIENTS

An independent validation cohort of 526 nondiabetic patients with a clinical indication for coronary angiography.

MEASUREMENTS

Receiver-operating characteristic (ROC) analysis of classifier score measured by real-time polymerase chain reaction, additivity to clinical factors, and reclassification of patient disease likelihood versus disease status defined by quantitative coronary angiography. Obstructive CAD was defined as 50% or greater stenosis in 1 or more major coronary arteries by quantitative coronary angiography.

RESULTS

The area under the ROC curve (AUC) was 0.70 ± 0.02 (P < 0.001); the test added to clinical variables (Diamond-Forrester method) (AUC, 0.72 with the test vs. 0.66 without; P = 0.003) and added somewhat to an expanded clinical model (AUC, 0.745 with the test vs. 0.732 without; P = 0.089). The test improved net reclassification over both the Diamond-Forrester method and the expanded clinical model (P < 0.001). At a score threshold that corresponded to a 20% likelihood of obstructive CAD (14.75), the sensitivity and specificity were 85% and 43% (yielding a negative predictive value of 83% and a positive predictive value of 46%), with 33% of patient scores below this threshold.

LIMITATION

Patients with chronic inflammatory disorders, elevated levels of leukocytes or cardiac protein markers, or diabetes were excluded.

CONCLUSION

A noninvasive whole-blood test based on gene expression and demographic characteristics may be useful for assessing obstructive CAD in nondiabetic patients without known CAD.

PRIMARY FUNDING SOURCE

CardioDx.

Complex exon-intron marking by histone modifications is not determined solely by nucleosome distribution

It has recently been shown that nucleosome distribution, histone modifications and RNA polymerase II (Pol II) occupancy show preferential association with exons (“exon-intron marking”), linking chromatin structure and function to co-transcriptional splicing in a variety of eukaryotes. Previous ChIP-sequencing studies suggested that these marking patterns reflect the nucleosomal landscape. By analyzing ChIP-chip datasets across the human genome in three cell types, we have found that this marking system is far more complex than previously observed. We show here that a range of histone modifications and Pol II are preferentially associated with exons. However, there is noticeable cell-type specificity in the degree of exon marking by histone modifications and, surprisingly, this is also reflected in some histone modifications patterns showing biases towards introns. Exon-intron marking is laid down in the absence of transcription on silent genes, with some marking biases changing or becoming reversed for genes expressed at different levels. Furthermore, the relationship of this marking system with splicing is not simple, with only some histone modifications reflecting exon usage/inclusion, while others mirror patterns of exon exclusion. By examining nucleosomal distributions in all three cell types, we demonstrate that these histone modification patterns cannot solely be accounted for by differences in nucleosome levels between exons and introns. In addition, because of inherent differences between ChIP-chip array and ChIP-sequencing approaches, these platforms report different nucleosome distribution patterns across the human genome. Our findings confound existing views and point to active cellular mechanisms which dynamically regulate histone modification levels and account for exon-intron marking. We believe that these histone modification patterns provide links between chromatin accessibility, Pol II movement and co-transcriptional splicing.

Genome-wide meta-analyses identifies seven loci associated with platelet aggregation in response to agonists

Platelet function mediates both beneficial and harmful effects on human health, but few genes are known to contribute to variability in this process. We tested association of 2.5 million SNPs with platelet aggregation responses to three agonists (ADP, epinephrine and collagen) in two cohorts of European ancestry (N Collapse

Key Words Collapse

MESH Headings

• Adenosine Diphosphate/pharmacology
• Adrenergic alpha-Agonists/pharmacology
• Cohort Studies
• Collagen/pharmacology
• Epinephrine/pharmacology
• Female
• Genome-Wide Association Study
• Genotype
• Humans
• Male
• Meta-Analysis as Topic
• Phenotype
• Platelet Aggregation/drug effects
• Platelet Aggregation/genetics
• Polymorphism, Single Nucleotide
• Prospective Studies

Collapse

Grants

• R01 HL087698 NHLBI NIH HHS
• N02 HL064278 NHLBI NIH HHS
• N01 HC025195 NHLBI NIH HHS
• M01-RR000052 NCRR NIH HHS
• N01 HC025195 NHLBI NIH HHS
• M01 RR000052 NCRR NIH HHS
• U01 HL72518 NHLBI NIH HHS
• U01 HL072518 NHLBI NIH HHS
• Z99 HL999999 Intramural NIH HHS
• N02-HL-6-4278 NHLBI NIH HHS
• N01-HC-25195 NHLBI NIH HHS
• R01 HL048157 NHLBI NIH HHS
• M01 RR005280 NCRR NIH HHS

Collapse

Affiliation(s)

Andrew D Johnson National Heart, Lung, and Blood Institute's The Framingham Heart Study, Framingham, Massachusetts, USA.
Lisa R Yanek
Ming-Huei Chen
Nauder Faraday
Martin G Larson
Geoffrey Tofler
Shiow J Lin
Aldi T Kraja
Michael A Province
Qiong Yang
Diane M Becker
Christopher J O'Donnell
Lewis C Becker

Collapse

Collaborators Collapse

Gfi1-cells and circuits: unraveling transcriptional networks of development and disease

PURPOSE OF REVIEW

The review will integrate current knowledge of transcriptional circuits whose dysregulation leads to autoimmunity, neutropenia and leukemia.

RECENT FINDINGS

Growth factor independent-1 (Gfi1) is a transcriptional repressor with essential roles in controlling hematopoietic stem cell biology, myeloid and lymphoid differentiation and lymphocyte effector functions. Recent work has suggested that Gfi1 competes or collaborates with other transcription factors to modulate transcription programs and lineage decisions.

SUMMARY

Gfi1 is central to several transcriptional circuits whose dysregulation leads to abnormal or malignant hematopoiesis. These functional relationships are conserved from Drosophila development. Such conserved pathways represent central oncogenic or 'gatekeeper' pathways that are pivotal to understanding the process of cellular transformation, and illustrate key targets for clinical intervention.

Phorbol ester-evoked Ca2+ signaling in human platelets is via autocrine activation of P(2X1) receptors, not a novel non-capacitative Ca2+ entry

SUMMARY BACKGROUND

Platelets are reported to possess a protein kinase C (PKC)-dependent non-capacitative Ca(2+)entry (NCCE) pathway. The phorbol ester, phorbol, 12-myristate, 13-acetate (PMA) has been suggested to stimulate platelet NCCE. Recently we demonstrated important roles in store-operated Ca(2+)entry in human platelets for Na(+)/Ca(2+) exchangers (NCXs) and autocrine signaling between platelets after dense granule secretion. As PMA evokes dense granule secretion, we have investigated the role of NCXs and autocrine signaling in PMA-evoked Ca(2+)entry.

OBJECTIVES

To investigate the roles of NCXs and dense granule secretion in PMA-evoked Ca(2+)signaling in human platelets.

METHODS

Fura-2- or sodium-binding benzofuran isophthalate (SBFI)-loaded platelets were used to monitor cytosolic Ca(2+)or Na(+) concentrations. Dense granule secretion was monitored as ATP release using luciferin-luciferase.

RESULTS

The NCX inhibitors KB-R7943 or SN-6, and removal of extracellular Na(+), significantly reduced PMA-evoked Ca(2+)entry. PMA-evoked dense granule secretion was almost abolished by pretreatment with the PKC inhibitor Ro-31-8220 and significantly slowed by KB-R7943. The P(2X1) antagonists Ro-0437626 or MRS-2159, or desensitization of P(2X1) receptors by prior treatment with alpha,beta-Methylene-ATP or omitting apyrase from the medium, reduced PMA-evoked Ca(2+)entry. Ro-0437626 or chelation of extracellular Ca(2+) slowed but did not abolish PMA-evoked ATP release, indicating that PMA-evoked dense granule secretion does not require P(2X1) receptor activation but is accelerated by P(2X1)-mediated Ca(2+)entry. The presence of NCX3 in human platelets was demonstrated by Western blotting.

CONCLUSION

PMA-evoked Ca(2+)entry results from an NCX3-dependent dense granule secretion and subsequent P(2X1) receptor activation by secreted ATP, rather than activation of a novel NCCE pathway.

Polycomb repressor complex 2 regulates HOXA9 and HOXA10, activating ID2 in NK/T-cell lines

BACKGROUND

NK- and T-cells are closely related lymphocytes, originating from the same early progenitor cells during hematopoiesis. In these differentiation processes deregulation of developmental genes may contribute to leukemogenesis. Here, we compared expression profiles of NK- and T-cell lines for identification of aberrantly expressed genes in T-cell acute lymphoblastic leukemia (T-ALL) which physiologically regulate the differentiation program of the NK-cell lineage.

RESULTS

This analysis showed high expression levels of HOXA9, HOXA10 and ID2 in NK-cell lines in addition to T-cell line LOUCY, suggesting leukemic deregulation therein. Overexpression experiments, chromatin immuno-precipitation and promoter analysis demonstrated that HOXA9 and HOXA10 directly activated expression of ID2. Concomitantly elevated expression levels of HOXA9 and HOXA10 together with ID2 in cell lines containing MLL translocations confirmed this form of regulation in both ALL and acute myeloid leukemia. Overexpression of HOXA9, HOXA10 or ID2 resulted in repressed expression of apoptosis factor BIM. Furthermore, profiling data of genes coding for chromatin regulators of homeobox genes, including components of polycomb repressor complex 2 (PRC2), indicated lacking expression of EZH2 in LOUCY and exclusive expression of HOP in NK-cell lines. Subsequent treatment of T-cell lines JURKAT and LOUCY with DZNep, an inhibitor of EZH2/PRC2, resulted in elevated and unchanged HOXA9/10 expression levels, respectively. Moreover, siRNA-mediated knockdown of EZH2 in JURKAT enhanced HOXA10 expression, confirming HOXA10-repression by EZH2. Additionally, profiling data and overexpression analysis indicated that reduced expression of E2F cofactor TFDP1 contributed to the lack of EZH2 in LOUCY. Forced expression of HOP in JURKAT cells resulted in reduced HOXA10 and ID2 expression levels, suggesting enhancement of PRC2 repression.

CONCLUSIONS

Our results show that major differentiation factors of the NK-cell lineage, including HOXA9, HOXA10 and ID2, were (de)regulated via PRC2 which therefore contributes to T-cell leukemogenesis.

Towards a more functional concept of causality in cancer research

Advances in molecular technologies challenge the different concepts of causality in biology, epidemiology and multistage mathematical models. The lack of integration of the different aspects of causality into a common framework could postpone our attempts to build a human causal model of carcinogenesis. We present here some aspects of differences in methodology, terminology and traditions between the scientific disciplines and propose a research strategy using functional analyses of the transcriptome and epigenetics to illuminate causality in complex biological systems. Overcoming the challenges of biological material collection suitable for such analyses into a prospective design, this could give unique opportunities for verification of mechanistic information from basic biological research in a human model system. The ultimate goal is to obtain a dynamic causal description of the different carcinogenesis stages. The success of this novel approach depends on the biological relationship between the gene expression of the somatic driver mutations or co-expressed genes in tumours and the gene expressions mirrored in peripheral blood along the different stages of carcinogenesis. The use of gene expression profiles and epigenetics could produce a functional concept of causality to explain the human multistage carcinogenic process.

Biomarker discovery in heterogeneous tissue samples -taking the in-silico deconfounding approach

Background

For heterogeneous tissues, such as blood, measurements of gene expression are confounded by relative proportions of cell types involved. Conclusions have to rely on estimation of gene expression signals for homogeneous cell populations, e.g. by applying micro-dissection, fluorescence activated cell sorting, or in-silico deconfounding. We studied feasibility and validity of a non-negative matrix decomposition algorithm using experimental gene expression data for blood and sorted cells from the same donor samples. Our objective was to optimize the algorithm regarding detection of differentially expressed genes and to enable its use for classification in the difficult scenario of reversely regulated genes. This would be of importance for the identification of candidate biomarkers in heterogeneous tissues.

Results

Experimental data and simulation studies involving noise parameters estimated from these data revealed that for valid detection of differential gene expression, quantile normalization and use of non-log data are optimal. We demonstrate the feasibility of predicting proportions of constituting cell types from gene expression data of single samples, as a prerequisite for a deconfounding-based classification approach.

Classification cross-validation errors with and without using deconfounding results are reported as well as sample-size dependencies. Implementation of the algorithm, simulation and analysis scripts are available.

Conclusions

The deconfounding algorithm without decorrelation using quantile normalization on non-log data is proposed for biomarkers that are difficult to detect, and for cases where confounding by varying proportions of cell types is the suspected reason. In this case, a deconfounding ranking approach can be used as a powerful alternative to, or complement of, other statistical learning approaches to define candidate biomarkers for molecular diagnosis and prediction in biomedicine, in realistically noisy conditions and with moderate sample sizes.

Platelets release novel thiol isomerase enzymes which are recruited to the cell surface following activation

The thiol isomerase enzymes protein disulphide isomerase (PDI) and endoplasmic reticulum protein 5 (ERp5) are released by resting and activated platelets. These re-associate with the cell surface where they modulate a range of platelet responses including adhesion, secretion and aggregation. Recent studies suggest the existence of yet uncharacterised platelet thiol isomerase proteins. This study aimed to identify which other thiol isomerase enzymes are present in human platelets. Through the use of immunoblotting, flow cytometry, cell-surface biotinylation and gene array analysis, we report the presence of five additional thiol isomerases in human and mouse platelets and megakaryocytes, namely; ERp57, ERp72, ERp44, ERp29 and TMX3. ERp72, ERp57, ERp44 and ERp29 are released by platelets and relocate to the cell surface following platelet activation. The transmembrane thiol isomerase TMX3 was also detected on the platelet surface but does not increase following activation. Extracellular PDI is also implicated in the regulation of coagulation by the modulation of tissue factor activity. ERp57 was identified within platelet-derived microparticle fractions, suggesting that ERp57 may also be involved in the regulation of coagulation as well as platelet function. These data collectively implicate the expanding family of platelet-surface thiol isomerases in the regulation of haemostasis.

Canonical Wnt signaling negatively regulates platelet function

Wnts regulate important intracellular signaling events, and dysregulation of the Wnt pathway has been linked to human disease. Here, we uncover numerous Wnt canonical effectors in human platelets where Wnts, their receptors, and downstream signaling components have not been previously described. We demonstrate that the Wnt3a ligand inhibits platelet adhesion, activation, dense granule secretion, and aggregation. Wnt3a also altered platelet shape change and inhibited the activation of the small GTPase RhoA. In addition, we found the Wnt-beta-catenin signaling pathway to be functional in platelets. Finally, disruption of the Wnt Frizzled 6 receptor in the mouse resulted in a hyperactivatory platelet phenotype and a reduced sensitivity to Wnt3a. Taken together our studies reveal a novel functional role for Wnt signaling in regulating anucleate platelet function and may provide a tractable target for future antiplatelet therapy.

Multiple loci influence erythrocyte phenotypes in the CHARGE Consortium

Measurements of erythrocytes within the blood are important clinical traits and can indicate various hematological disorders. We report here genome-wide association studies (GWAS) for six erythrocyte traits, including hemoglobin concentration (Hb), hematocrit (Hct), mean corpuscular volume (MCV), mean corpuscular hemoglobin (MCH), mean corpuscular hemoglobin concentration (MCHC) and red blood cell count (RBC). We performed an initial GWAS in cohorts of the CHARGE Consortium totaling 24,167 individuals of European ancestry and replication in additional independent cohorts of the HaemGen Consortium totaling 9,456 individuals. We identified 23 loci significantly associated with these traits in a meta-analysis of the discovery and replication cohorts (combined P values ranging from 5 x 10(-8) to 7 x 10(-86)). Our findings include loci previously associated with these traits (HBS1L-MYB, HFE, TMPRSS6, TFR2, SPTA1) as well as new associations (EPO, TFRC, SH2B3 and 15 other loci). This study has identified new determinants of erythrocyte traits, offering insight into common variants underlying variation in erythrocyte measures.

A functional genomics approach reveals novel quantitative trait loci associated with platelet signaling pathways

Platelet response to activation varies widely between individuals but shows interindividual consistency and strong heritability. The genetic basis of this variation has not been properly explored. We therefore systematically measured the effect on function of sequence variation in 97 candidate genes in the collagen and adenosine-diphosphate (ADP) signaling pathways. Resequencing of the genes in 48 European DNA samples nearly doubled the number of known single nucleotide polymorphisms (SNPs) and informed the selection of 1327 SNPs for genotyping in 500 healthy Northern European subjects with known platelet responses to collagen-related peptide (CRP-XL) and ADP. This identified 17 novel associations with platelet function (P < .005) accounting for approximately 46% of the variation in response. Further investigations with platelets of known genotype explored the mechanisms behind some of the associations. SNPs in PEAR1 associated with increased platelet response to CRP-XL and increased PEAR1 protein expression after platelet degranulation. The minor allele of a 3' untranslated region (UTR) SNP (rs2769668) in VAV3 was associated with higher protein expression (P = .03) and increased P-selectin exposure after ADP activation (P = .004). Furthermore the minor allele of the intronic SNP rs17786144 in ITPR1 modified Ca(2+) levels after activation with ADP (P < .004). These data provide novel insights into key hubs within platelet signaling networks.

The discovery of genes implicated in myocardial infarction

The era of Genome-Wide Association Studies (GWAS) commenced in 2007 with the study of the Wellcome Trust Case Control Consortium (WTCCC) which for the first time ever showed that risk loci can be identified by scanning the complete genome for sequence variation in large numbers of cases of disease and healthy controls. We and others have expanded on this effort and successfully identified the first 11 risk loci for myocardial infarction (MI) and coronary artery disease (CAD). Studies on quantitative traits provide an alternative approach to identify MI/CAD risk loci. This review captures the early successes in the emerging field of disease genomics.

Functional genomics in zebrafish permits rapid characterization of novel platelet membrane proteins

In this study, we demonstrate the suitability of the vertebrate Danio rerio (zebrafish) for functional screening of novel platelet genes in vivo by reverse genetics. Comparative transcript analysis of platelets and their precursor cell, the megakaryocyte, together with nucleated blood cell elements, endothelial cells, and erythroblasts, identified novel platelet membrane proteins with hitherto unknown roles in thrombus formation. We determined the phenotype induced by antisense morpholino oligonucleotide (MO)-based knockdown of 5 of these genes in a laser-induced arterial thrombosis model. To validate the model, the genes for platelet glycoprotein (GP) IIb and the coagulation protein factor VIII were targeted. MO-injected fish showed normal thrombus initiation but severely impaired thrombus growth, consistent with the mouse knockout phenotypes, and concomitant knockdown of both resulted in spontaneous bleeding. Knockdown of 4 of the 5 novel platelet proteins altered arterial thrombosis, as demonstrated by modified kinetics of thrombus initiation and/or development. We identified a putative role for BAMBI and LRRC32 in promotion and DCBLD2 and ESAM in inhibition of thrombus formation. We conclude that phenotypic analysis of MO-injected zebrafish is a fast and powerful method for initial screening of novel platelet proteins for function in thrombosis.

Identification of Tspan9 as a novel platelet tetraspanin and the collagen receptor GPVI as a component of tetraspanin microdomains

Platelets are essential for wound healing and inflammatory processes, but can also play a deleterious role by causing heart attack and stroke. Normal platelet activation is dependent on tetraspanins, a superfamily of glycoproteins that function as 'organisers' of cell membranes by recruiting other receptors and signalling proteins into tetraspanin-enriched microdomains. However, our understanding of how tetraspanin microdomains regulate platelets is hindered by the fact that only four of the 33 mammalian tetraspanins have been identified in platelets. This is because of a lack of antibodies to most tetraspanins and difficulties in measuring mRNA, due to low levels in this anucleate cell. To identify potentially platelet-expressed tetraspanins, mRNA was measured in their nucleated progenitor cell, the megakaryocyte, using serial analysis of gene expression and DNA microarrays. Amongst 19 tetraspanins identified in megakaryocytes, Tspan9, a previously uncharacterized tetraspanin, was relatively specific to these cells. Through generating the first Tspan9 antibodies, Tspan9 expression was found to be tightly regulated in platelets. The relative levels of CD9, CD151, Tspan9 and CD63 were 100, 14, 6 and 2 respectively. Since CD9 was expressed at 49000 cell surface copies per platelet, this suggested a copy number of 2800 Tspan9 molecules. Finally, Tspan9 was shown to be a component of tetraspanin microdomains that included the collagen receptor GPVI (glycoprotein VI) and integrin alpha6beta1, but not the von Willebrand receptor GPIbalpha or the integrins alphaIIbbeta3 or alpha2beta1. These findings suggest a role for Tspan9 in regulating platelet function in concert with other platelet tetraspanins and their associated proteins.