Combinatorial transcription factor (TF) interactions control cellular phenotypes and, therefore, underpin stem cell formation, maintenance, and differentiation. Here, we report the genome-wide binding patterns and combinatorial interactions for ten key regulators of blood stem/progenitor cells (SCL/TAL1, LYL1, LMO2, GATA2, RUNX1, MEIS1, PU.1, ERG, FLI-1, and GFI1B), thus providing the most comprehensive TF data set for any adult stem/progenitor cell type to date. Genome-wide computational analysis of complex binding patterns, followed by functional validation, revealed the following: first, a previously unrecognized combinatorial interaction between a heptad of TFs (SCL, LYL1, LMO2, GATA2, RUNX1, ERG, and FLI-1). Second, we implicate direct protein-protein interactions between four key regulators (RUNX1, GATA2, SCL, and ERG) in stabilizing complex binding to DNA. Third, Runx1(+/-)::Gata2(+/-) compound heterozygous mice are not viable with severe hematopoietic defects at midgestation. Taken together, this study demonstrates the power of genome-wide analysis in generating novel functional insights into the transcriptional control of stem and progenitor cells.

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.

Mean platelet volume (MPV) is increased in myocardial and cerebral infarction and is an independent and strong predictor for postevent morbidity and mortality. We conducted a genome-wide association study (GWAS), the KORA (Kooperative Gesundheitsforschung in der Region Augsburg) F3 500K study, and found MPV to be strongly associated with three common single-nucleotide polymorphisms (SNPs): rs7961894 located within intron 3 of WDR66 on chromosome 12q24.31, rs12485738 upstream of the ARHGEF3 on chromosome 3p13-p21, and rs2138852 located upstream of TAOK1 on chromosome 17q11.2. We replicated all three SNPs in another GWAS from the UK and in two population-based samples from Germany. In a combined analysis including 10,048 subjects, the SNPs had p values of 7.24 x 10(-48) for rs7961894, 3.81 x 10(-27) for rs12485738, and 7.19 x 10(-28) for rs2138852. These three quantitative trait loci together accounted for 4%-5% of the variance in MPV. In-depth sequence analysis of WDR66 in 382 samples from the extremes revealed 20 new variants and a haplotype with three coding SNPs and one SNP at the transcription start site associated with MPV (p = 6.8 x 10(-5)). In addition, expression analysis indicated a direct correlation of WDR66 transcripts and MPV. These findings may not only enhance our understanding of platelet activation and function, but may also provide a focus for several novel research avenues.

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.

Platelets are anuclear cells that are essential for blood clotting. They are produced by large polyploid precursor cells called megakaryocytes. Previous genome-wide association studies in nearly 70,000 individuals indicated that single nucleotide variants (SNVs) in the gene encoding the actin cytoskeletal regulator tropomyosin 4 (TPM4) exert an effect on the count and volume of platelets. Platelet number and volume are independent risk factors for heart attack and stroke. Here, we have identified 2 unrelated families in the BRIDGE Bleeding and Platelet Disorders (BPD) collection who carry a TPM4 variant that causes truncation of the TPM4 protein and segregates with macrothrombocytopenia, a disorder characterized by low platelet count. N-Ethyl-N-nitrosourea-induced (ENU-induced) missense mutations in Tpm4 or targeted inactivation of the Tpm4 locus led to gene dosage-dependent macrothrombocytopenia in mice. All other blood cell counts in Tpm4-deficient mice were normal. Insufficient TPM4 expression in human and mouse megakaryocytes resulted in a defect in the terminal stages of platelet production and had a mild effect on platelet function. Together, our findings demonstrate a nonredundant role for TPM4 in platelet biogenesis in humans and mice and reveal that truncating variants in TPM4 cause a previously undescribed dominant Mendelian platelet disorder.

Raised blood pressure (BP) is a major risk factor for cardiovascular disease. Previous studies have identified 47 distinct genetic variants robustly associated with BP, but collectively these explain only a few percent of the heritability for BP phenotypes. To find additional BP loci, we used a bespoke gene-centric array to genotype an independent discovery sample of 25,118 individuals that combined hypertensive case-control and general population samples. We followed up four SNPs associated with BP at our p < 8.56 × 10(-7) study-specific significance threshold and six suggestively associated SNPs in a further 59,349 individuals. We identified and replicated a SNP at LSP1/TNNT3, a SNP at MTHFR-NPPB independent (r(2) = 0.33) of previous reports, and replicated SNPs at AGT and ATP2B1 reported previously. An analysis of combined discovery and follow-up data identified SNPs significantly associated with BP at p < 8.56 × 10(-7) at four further loci (NPR3, HFE, NOS3, and SOX6). The high number of discoveries made with modest genotyping effort can be attributed to using a large-scale yet targeted genotyping array and to the development of a weighting scheme that maximized power when meta-analyzing results from samples ascertained with extreme phenotypes, in combination with results from nonascertained or population samples. Chromatin immunoprecipitation and transcript expression data highlight potential gene regulatory mechanisms at the MTHFR and NOS3 loci. These results provide candidates for further study to help dissect mechanisms affecting BP and highlight the utility of studying SNPs and samples that are independent of those studied previously even when the sample size is smaller than that in previous studies.

The cost of whole-genome bisulfite sequencing (WGBS) remains a bottleneck for many studies and it is therefore imperative to extract as much information as possible from a given dataset. This is particularly important because even at the recommend 30X coverage for reference methylomes, up to 50% of high-resolution features such as differentially methylated positions (DMPs) cannot be called with current methods as determined by saturation analysis. To address this limitation, we have developed a tool that dynamically segments WGBS methylomes into blocks of comethylation (COMETs) from which lost information can be recovered in the form of differentially methylated COMETs (DMCs). Using this tool, we demonstrate recovery of ∼30% of the lost DMP information content as DMCs even at very low (5X) coverage. This constitutes twice the amount that can be recovered using an existing method based on differentially methylated regions (DMRs). In addition, we explored the relationship between COMETs and haplotypes in lymphoblastoid cell lines of African and European origin. Using best fit analysis, we show COMETs to be correlated in a population-specific manner, suggesting that this type of dynamic segmentation may be useful for integrated (epi)genome-wide association studies in the future.

Obesity is associated with vitamin D deficiency, and both are areas of active public health concern. We explored the causality and direction of the relationship between body mass index (BMI) and 25-hydroxyvitamin D [25(OH)D] using genetic markers as instrumental variables (IVs) in bi-directional Mendelian randomization (MR) analysis.

The burden of cardiovascular disease (CVD) cannot be fully addressed by therapy targeting known pathophysiological pathways. Even with stringent control of all risk factors CVD events are only diminished by half. A number of additional pathways probably play a role in the development of CVD and might serve as novel therapeutic targets. Genome wide expression studies represent a powerful tool to identify such novel pathways. We compared the expression profiles in monocytes from twenty two young male patients with premature familial CAD with those from controls matched for age, sex and smoking status, without a family history of CVD. Since all patients were on statins and aspirin treatment, potentially affecting the expression of genes in monocytes, twelve controls were subsequently treated with simvastatin and aspirin for 6 and 2 weeks, respectively. By whole genome expression arrays six genes were identified to have differential expression in the monocytes of patients versus controls; ABCA1, ABCG1 and RGS1 were downregulated in patients, whereas ADRB2, FOLR3 and GSTM1 were upregulated. Differential expression of all genes, apart from GSTM1, was confirmed by qPCR. Aspirin and statins altered gene expression of ABCG1 and ADBR2. All finding were validated in a second group of twenty four patients and controls. Differential expression of ABCA1, RSG1 and ADBR2 was replicated. In conclusion, we identified these 3 genes to be expressed differently in CAD cases which might play a role in the pathogenesis of atherosclerotic vascular disease.

We aimed to assess whether pri-miRNA SNPs (miSNPs) could influence monocyte gene expression, either through marginal association or by interacting with polymorphisms located in 3'UTR regions (3utrSNPs). We then conducted a genome-wide search for marginal miSNPs effects and pairwise miSNPs × 3utrSNPs interactions in a sample of 1,467 individuals for which genome-wide monocyte expression and genotype data were available. Statistical associations that survived multiple testing correction were tested for replication in an independent sample of 758 individuals with both monocyte gene expression and genotype data. In both studies, the hsa-mir-1279 rs1463335 was found to modulate in cis the expression of LYZ and in trans the expression of CNTN6, CTRC, COPZ2, KRT9, LRRFIP1, NOD1, PCDHA6, ST5 and TRAF3IP2 genes, supporting the role of hsa-mir-1279 as a regulator of several genes in monocytes. In addition, we identified two robust miSNPs × 3utrSNPs interactions, one involving HLA-DPB1 rs1042448 and hsa-mir-219-1 rs107822, the second the H1F0 rs1894644 and hsa-mir-659 rs5750504, modulating the expression of the associated genes.As some of the aforementioned genes have previously been reported to reside at disease-associated loci, our findings provide novel arguments supporting the hypothesis that the genetic variability of miRNAs could also contribute to the susceptibility to human diseases.

Sphingolipids are fundamental to membrane trafficking, apoptosis, and cell differentiation and proliferation. KDSR or 3-keto-dihydrosphingosine reductase is an essential enzyme for de novo sphingolipid synthesis, and pathogenic mutations in KDSR result in the severe skin disorder erythrokeratodermia variabilis et progressiva-4 Four of the eight reported cases also had thrombocytopenia but the underlying mechanism has remained unexplored. Here we expand upon the phenotypic spectrum of KDSR deficiency with studies in two siblings with novel compound heterozygous variants associated with thrombocytopenia, anemia, and minimal skin involvement. We report a novel phenotype of progressive juvenile myelofibrosis in the propositus, with spontaneous recovery of anemia and thrombocytopenia in the first decade of life. Examination of bone marrow biopsies showed megakaryocyte hyperproliferation and dysplasia. Megakaryocytes obtained by culture of CD34+ stem cells confirmed hyperproliferation and showed reduced proplatelet formation. The effect of KDSR insufficiency on the sphingolipid profile was unknown, and was explored in vivo and in vitro by a broad metabolomics screen that indicated activation of an in vivo compensatory pathway that leads to normalization of downstream metabolites such as ceramide. Differentiation of propositus-derived induced pluripotent stem cells to megakaryocytes followed by expression of functional KDSR showed correction of the aberrant cellular and biochemical phenotypes, corroborating the critical role of KDSR in proplatelet formation. Finally, Kdsr depletion in zebrafish recapitulated the thrombocytopenia and showed biochemical changes similar to those observed in the affected siblings. These studies support an important role for sphingolipids as regulators of cytoskeletal organization during megakaryopoiesis and proplatelet formation.

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.

Blood cells play essential roles in human health, underpinning physiological processes such as immunity, oxygen transport, and clotting, which when perturbed cause a significant global health burden. Here we integrate data from UK Biobank and a large-scale international collaborative effort, including data for 563,085 European ancestry participants, and discover 5,106 new genetic variants independently associated with 29 blood cell phenotypes covering a range of variation impacting hematopoiesis. We holistically characterize the genetic architecture of hematopoiesis, assess the relevance of the omnigenic model to blood cell phenotypes, delineate relevant hematopoietic cell states influenced by regulatory genetic variants and gene networks, identify novel splice-altering variants mediating the associations, and assess the polygenic prediction potential for blood traits and clinical disorders at the interface of complex and Mendelian genetics. These results show the power of large-scale blood cell trait GWAS to interrogate clinically meaningful variants across a wide allelic spectrum of human variation.

Restless legs syndrome is a prevalent chronic neurological disorder with potentially severe mental and physical health consequences. Clearer understanding of the underlying pathophysiology is needed to improve treatment options. We did a meta-analysis of genome-wide association studies (GWASs) to identify potential molecular targets.

Linking non-coding genetic variants associated with the risk of diseases or disease-relevant traits to target genes is a crucial step to realize GWAS potential in the introduction of precision medicine. Here we set out to determine the mechanisms underpinning variant association with platelet quantitative traits using cell type-matched epigenomic data and promoter long-range interactions. We identify potential regulatory functions for 423 of 565 (75%) non-coding variants associated with platelet traits and we demonstrate, through ex vivo and proof of principle genome editing validation, that variants in super enhancers play an important role in controlling archetypical platelet functions.

Iron is essential for many biological functions and iron deficiency and overload have major health implications. We performed a meta-analysis of three genome-wide association studies from Iceland, the UK and Denmark of blood levels of ferritin (N = 246,139), total iron binding capacity (N = 135,430), iron (N = 163,511) and transferrin saturation (N = 131,471). We found 62 independent sequence variants associating with iron homeostasis parameters at 56 loci, including 46 novel loci. Variants at DUOX2, F5, SLC11A2 and TMPRSS6 associate with iron deficiency anemia, while variants at TF, HFE, TFR2 and TMPRSS6 associate with iron overload. A HBS1L-MYB intergenic region variant associates both with increased risk of iron overload and reduced risk of iron deficiency anemia. The DUOX2 missense variant is present in 14% of the population, associates with all iron homeostasis biomarkers, and increases the risk of iron deficiency anemia by 29%. The associations implicate proteins contributing to the main physiological processes involved in iron homeostasis: iron sensing and storage, inflammation, absorption of iron from the gut, iron recycling, erythropoiesis and bleeding/menstruation.

Trait-associated genetic variants affect complex phenotypes primarily via regulatory mechanisms on the transcriptome. To investigate the genetics of gene expression, we performed cis- and trans-expression quantitative trait locus (eQTL) analyses using blood-derived expression from 31,684 individuals through the eQTLGen Consortium. We detected cis-eQTL for 88% of genes, and these were replicable in numerous tissues. Distal trans-eQTL (detected for 37% of 10,317 trait-associated variants tested) showed lower replication rates, partially due to low replication power and confounding by cell type composition. However, replication analyses in single-cell RNA-seq data prioritized intracellular trans-eQTL. Trans-eQTL exerted their effects via several mechanisms, primarily through regulation by transcription factors. Expression of 13% of the genes correlated with polygenic scores for 1,263 phenotypes, pinpointing potential drivers for those traits. In summary, this work represents a large eQTL resource, and its results serve as a starting point for in-depth interpretation of complex phenotypes.

Cardiometabolic diseases are frequently polygenic in architecture, comprising a large number of risk alleles with small effects spread across the genome1-3. Polygenic scores (PGS) aggregate these into a metric representing an individual's genetic predisposition to disease. PGS have shown promise for early risk prediction4-7 and there is an open question as to whether PGS can also be used to understand disease biology8. Here, we demonstrate that cardiometabolic disease PGS can be used to elucidate the proteins underlying disease pathogenesis. In 3,087 healthy individuals, we found that PGS for coronary artery disease, type 2 diabetes, chronic kidney disease and ischaemic stroke are associated with the levels of 49 plasma proteins. Associations were polygenic in architecture, largely independent of cis and trans protein quantitative trait loci and present for proteins without quantitative trait loci. Over a follow-up of 7.7 years, 28 of these proteins associated with future myocardial infarction or type 2 diabetes events, 16 of which were mediators between polygenic risk and incident disease. Twelve of these were druggable targets with therapeutic potential. Our results demonstrate the potential for PGS to uncover causal disease biology and targets with therapeutic potential, including those that may be missed by approaches utilizing information at a single locus.

The implementation of high-throughput sequencing (HTS) technologies in research and diagnostic laboratories has linked many new genes to rare bleeding, thrombotic, and platelet disorders (BTPD), and revealed multiple genetic variants linked to those disorders, many of them being of uncertain pathogenicity when considering the accepted evidence (variant consequence, frequency in control datasets, number of reported patients, prediction models, and functional assays). The sequencing effort has also resulted in resources for gathering disease-causing variants associated with specific genes, but for BTPD, such well-curated databases exist only for a few genes. On the other hand, submissions by individuals or diagnostic laboratories to the variant database ClinVar are hampered by the lack of a submission process tailored to capture the specific features of hemostatic diseases. As we move toward the implementation of HTS in the diagnosis of BTPD, the Scientific and Standardization Committee for Genetics in Thrombosis and Haemostasis has developed and tested a REDCap-based interface, aimed at the community, to submit curated genetic variants for diagnostic-grade BTPD genes. Here, we describe the use of the interface and the initial submission of 821 variants from 30 different centers covering 14 countries. This open-access variant resource will be shared with the community to improve variant classification and regular bulk data transfer to ClinVar.

Each day, about 1012 erythrocytes and platelets are released into the bloodstream. This substantial output from hematopoietic stem cells is tightly regulated by transcriptional and epigenetic factors. Whether and how circular RNAs (circRNAs) contribute to the differentiation and/or identity of hematopoietic cells is to date not known. We recently reported that erythrocytes and platelets contain the highest levels and numbers of circRNAs among hematopoietic cells. Here, we provide the first detailed analysis of circRNA expression during erythroid and megakaryoid differentiation. CircRNA expression not only significantly increased upon enucleation, but also had limited overlap between progenitor cells and mature cells, suggesting that circRNA expression stems from regulated processes rather than resulting from mere accumulation. To study circRNA function in hematopoiesis, we first compared the expression levels of circRNAs with the translation efficiency of their mRNA counterpart. We found that only one out of 2531 (0.04%) circRNAs associated with mRNA-translation regulation. Furthermore, irrespective of thousands of identified putative open reading frames, deep ribosome-footprinting sequencing, and mass spectrometry analysis provided little evidence for translation of endogenously expressed circRNAs. In conclusion, circRNAs alter their expression profile during terminal hematopoietic differentiation, yet their contribution to regulate cellular processes remains enigmatic.