Nuclear adaptor Ldb1 regulates a transcriptional program essential for the maintenance of hematopoietic stem cells

Enhancer looping protein LDB1 modulates MYB expression in T-ALL cell lines in vitro by cooperating with master transcription factors

BACKGROUND

Despite significant progress in the prognosis of pediatric T-cell acute lymphoblastic leukemia (T-ALL) in recent decades, a notable portion of children still confronts challenges such as treatment resistance and recurrence, leading to limited options and a poor prognosis. LIM domain-binding protein 1 (LDB1) has been confirmed to exert a crucial role in various physiological and pathological processes. In our research, we aim to elucidate the underlying function and mechanisms of LDB1 within the background of T-ALL.

METHODS

Employing short hairpin RNA (shRNA) techniques, we delineated the functional impact of LDB1 in T-ALL cell lines. Through the application of RNA-Seq, CUT&Tag, and immunoprecipitation assays, we scrutinized master transcription factors cooperating with LDB1 and identified downstream targets under LDB1 regulation.

RESULTS

LDB1 emerges as a critical transcription factor co-activator in cell lines derived from T-ALL. It primarily collaborates with master transcription factors (ERG, ETV6, IRF1) to cooperatively regulate the transcription of downstream target genes. Both in vitro and in vivo experiments affirm the essential fuction of LDB1 in the proliferation and survival of cell lines derived from T-ALL, with MYB identified as a significant downstream target of LDB1.

CONCLUSIONS

To sum up, our research establishes the pivotal fuction of LDB1 in the tumorigenesis and progression of T-ALL cell lines. Mechanistic insights reveal that LDB1 cooperates with ERG, ETV6, and IRF1 to modulate the expression of downstream effector genes. Furthermore, LDB1 controls MYB through remote enhancer modulation, providing valuable mechanistic insights into its involvement in the progression of T-ALL.

Single-cell RNA sequencing illuminates the ontogeny, conservation and diversification of cartilaginous and bony fish lymphocytes

Elucidating cellular architecture and cell-type evolution across species is central to understanding immune system function and susceptibility to disease. Adaptive immunity is a shared trait of the common ancestor of cartilaginous and bony fishes. However, evolutionary features of lymphocytes in these two jawed vertebrates remain unclear. Here, we present a single-cell RNA sequencing atlas of immune cells from cartilaginous (white-spotted bamboo shark) and bony (zebrafish and Chinese tongue sole) fishes. Cross-species comparisons show that the same cell types across different species exhibit similar transcriptional profiles. In the bamboo shark, we identify a phagocytic B cell population expressing several pattern recognition receptors, as well as a T cell sub-cluster co-expressing both T and B cell markers. In contrast to a division by function in the bony fishes, we show close linkage and poor functional specialization among lymphocytes in the cartilaginous fish. Our cross-species single-cell comparison presents a resource for uncovering the origin and evolution of the gnathostome immune system.

A newly identified gene Ahed plays essential roles in murine haematopoiesis

The development of haematopoiesis involves the coordinated action of numerous genes, some of which are implicated in haematological malignancies. However, the biological function of many genes remains elusive and unknown functional genes are likely to remain to be uncovered. Here, we report a previously uncharacterised gene in haematopoiesis, identified by screening mutant embryonic stem cells. The gene, 'attenuated haematopoietic development (Ahed)', encodes a nuclear protein. Conditional knockout (cKO) of Ahed results in anaemia from embryonic day 14.5 onward, leading to prenatal demise. Transplantation experiments demonstrate the incapacity of Ahed-deficient haematopoietic cells to reconstitute haematopoiesis in vivo. Employing a tamoxifen-inducible cKO model, we further reveal that Ahed deletion impairs the intrinsic capacity of haematopoietic cells in adult mice. Ahed deletion affects various pathways, and published databases present cancer patients with somatic mutations in Ahed. Collectively, our findings underscore the fundamental roles of Ahed in lifelong haematopoiesis, implicating its association with malignancies.

MYB: A Key Transcription Factor in the Hematopoietic System Subject to Many Levels of Control

MYB is a master regulator and pioneer factor highly expressed in hematopoietic progenitor cells (HPCs) where it contributes to the reprogramming processes operating during hematopoietic development. MYB plays a complex role being involved in several lineages of the hematopoietic system. At the molecular level, the MYB gene is subject to intricate regulation at many levels through several enhancer and promoter elements, through transcriptional elongation control, as well as post-transcriptional regulation. The protein is modulated by post-translational modifications (PTMs) such as SUMOylation restricting the expression of its downstream targets. Together with a range of interaction partners, cooperating transcription factors (TFs) and epigenetic regulators, MYB orchestrates a fine-tuned symphony of genes expressed during various stages of haematopoiesis. At the same time, the complex MYB system is vulnerable, being a target for unbalanced control and cancer development.

Lin- PU.1dim GATA-1- defines haematopoietic stem cells with long-term multilineage reconstitution activity

Despite extensive characterization of the state and function of haematopoietic stem cells (HSCs), the use of transcription factors to define the HSC population is still limited. We show here that the HSC population in mouse bone marrow can be defined by the distinct expression levels of Spi1 and Gata1. By using a double fluorescence knock-in mouse model, PGdKI, in which the expression levels of PU.1 and GATA-1 are indicated by the expression of GFP and mCherry, respectively, we uncover that the HSCs with lymphoid and myeloid repopulating activity are specifically enriched in a Lin- PU.1dim GATA-1- (LPG) cell subset. In vivo competitive repopulation assays demonstrate that bone marrow cells gated by LPG exhibit haematopoietic reconstitution activity which is comparable to that of classical Lin- Sca1+ c-kit+ (LSK). The integrated analysis of single-cell RNA sequence data from LPG and LSK-gated cells reveals that a transcriptional network governed by core TFs contributes to regulation of HSC multipotency. These discoveries provide new clues for HSC characterization and functional study.

LMO2 promotes the development of AML through interaction with transcription co-regulator LDB1

One of the characteristics of leukemia is that it contains multiple rearrangements of signal transduction genes and overexpression of non-mutant genes, such as transcription factors. As an important regulator of hematopoietic stem cell development and erythropoiesis, LMO2 is considered an effective carcinogenic driver in T cell lines and a marker of poor prognosis in patients with AML with normal karyotype. LDB1 is a key factor in the transformation of thymocytes into T-ALL induced by LMO2, and enhances the stability of carcinogenic related proteins in leukemia. However, the function and mechanism of LMO2 and LDB1 in AML remains unclear. Herein, the LMO2 gene was knocked down to observe its effects on proliferation, survival, and colony formation of NB4, Kasumi-1 and K562 cell lines. Using mass spectrometry and IP experiments, our results showed the presence of LMO2/LDB1 protein complex in AML cell lines, which is consistent with previous studies. Furthermore, in vitro and in vivo experiments revealed that LDB1 is essential for the proliferation and survival of AML cell lines. Analysis of RNA-seq and ChIP-Seq results showed that LDB1 could regulate apoptosis-related genes, including LMO2. In LDB1-deficient AML cell lines, the overexpression of LMO2 partially compensates for the proliferation inhibition. In summary, our findings revealed that LDB1 played an important role in AML as an oncogene, and emphasize the potential importance of the LMO2/LDB1 complex in clinical treatment of patients with AML.

Cell fate decisions, transcription factors and signaling during early retinal development

The development of the vertebrate eyes is a complex process starting from anterior-posterior and dorso-ventral patterning of the anterior neural tube, resulting in the formation of the eye field. Symmetrical separation of the eye field at the anterior neural plate is followed by two symmetrical evaginations to generate a pair of optic vesicles. Next, reciprocal invagination of the optic vesicles with surface ectoderm-derived lens placodes generates double-layered optic cups. The inner and outer layers of the optic cups develop into the neural retina and retinal pigment epithelium (RPE), respectively. In vitro produced retinal tissues, called retinal organoids, are formed from human pluripotent stem cells, mimicking major steps of retinal differentiation in vivo. This review article summarizes recent progress in our understanding of early eye development, focusing on the formation the eye field, optic vesicles, and early optic cups. Recent single-cell transcriptomic studies are integrated with classical in vivo genetic and functional studies to uncover a range of cellular mechanisms underlying early eye development. The functions of signal transduction pathways and lineage-specific DNA-binding transcription factors are dissected to explain cell-specific regulatory mechanisms underlying cell fate determination during early eye development. The functions of homeodomain (HD) transcription factors Otx2, Pax6, Lhx2, Six3 and Six6, which are required for early eye development, are discussed in detail. Comprehensive understanding of the mechanisms of early eye development provides insight into the molecular and cellular basis of developmental ocular anomalies, such as optic cup coloboma. Lastly, modeling human development and inherited retinal diseases using stem cell-derived retinal organoids generates opportunities to discover novel therapies for retinal diseases.

LMO2 expression is frequent in T-lymphoblastic leukemia and correlates with survival, regardless of T-cell stage

T- lymphoblastic leukemia/lymphoma (T-LL) is an aggressive malignancy of immature T-cells with poor overall survival (OS) and in need of new therapies. LIM-domain only 2 (LMO2) is a critical regulator of hematopoietic cell development that can be overexpressed in T-LL due to chromosomal abnormalities. Deregulated LMO2 expression contributes to T-LL development by inducing block of T-cell differentiation and continuous thymocyte self-renewal. However, LMO2 expression and its biologic significance in T-LL remain largely unknown. We analyzed LMO2 expression in 100 initial and follow-up biopsies of T-LL from 67 patients, including 31 (46%) early precursor T-cell (ETP)-ALL, 26 (39%) cortical and 10 (15%) medullary type. LMO2 expression was present in 50 (74.6%) initial biopsies with an average of 87% positive tumor cells (range 30-100%). LMO2 expression in ETP, medullary and cortical T-LLs was not statistically different. In patients with biopsies after initial therapy, LMO2 expression was stable. LMO2 expression was associated with longer OS (p = 0.048) regardless of T-lymphoblast stage or other clinicopathologic features. These findings indicate that LMO2 is a promising new prognostic marker that could predict patients' outcomes and potentially be targeted for novel chemotherapy, i.e. PARP1/2 inhibitors, which have been shown to enhance chemotherapy sensitivity in LMO2 expressing diffuse large B cell lymphoma (DLBCL) tumors by decreasing DNA repair efficiency.

HIRA, a DiGeorge Syndrome Candidate Gene, Confers Proper Chromatin Accessibility on HSCs and Supports All Stages of Hematopoiesis

HIRA is a histone chaperone that deposits the histone variant H3.3 in transcriptionally active genes. In DiGeorge syndromes, a DNA stretch encompassing HIRA is deleted. The syndromes manifest varied abnormalities, including immunodeficiency and thrombocytopenia. HIRA is essential in mice, as total knockout (KO) results in early embryonic death. However, the role of HIRA in hematopoiesis is poorly understood. We investigate hematopoietic cell-specific Hira deletion in mice and show that it dramatically reduces bone marrow hematopoietic stem cells (HSCs), resulting in anemia, thrombocytopenia, and lymphocytopenia. In contrast, fetal hematopoiesis is normal in Hira-KO mice, although fetal HSCs lack the reconstitution capacity. Transcriptome analysis reveals that HIRA is required for expression of many transcription factors and signaling molecules critical for HSCs. ATAC-seq analysis demonstrates that HIRA establishes HSC-specific DNA accessibility, including the SPIB/PU.1 sites. Together, HIRA provides a chromatin environment essential for HSCs, thereby steering their development and survival.

Ldb1 is required for Lmo2 oncogene-induced thymocyte self-renewal and T-cell acute lymphoblastic leukemia

Prolonged or enhanced expression of the proto-oncogene Lmo2 is associated with a severe form of T-cell acute lymphoblastic leukemia (T-ALL), designated early T-cell precursor ALL, which is characterized by the aberrant self-renewal and subsequent oncogenic transformation of immature thymocytes. It has been suggested that Lmo2 exerts these effects by functioning as component of a multi-subunit transcription complex that includes the ubiquitous adapter Ldb1 along with b-HLH and/or GATA family transcription factors; however, direct experimental evidence for this mechanism is lacking. In this study, we investigated the importance of Ldb1 for Lmo2-induced T-ALL by conditional deletion of Ldb1 in thymocytes in an Lmo2 transgenic mouse model of T-ALL. Our results identify a critical requirement for Ldb1 in Lmo2-induced thymocyte self-renewal and thymocyte radiation resistance and for the transition of preleukemic thymocytes to overt T-ALL. Moreover, Ldb1 was also required for acquisition of the aberrant preleukemic ETP gene expression signature in immature Lmo2 transgenic thymocytes. Co-binding of Ldb1 and Lmo2 was detected at the promoters of key upregulated T-ALL driver genes (Hhex, Lyl1, and Nfe2) in preleukemic Lmo2 transgenic thymocytes, and binding of both Ldb1 and Lmo2 at these sites was reduced following Cre-mediated deletion of Ldb1. Together, these results identify a key role for Ldb1, a nonproto-oncogene, in T-ALL and support a model in which Lmo2-induced T-ALL results from failure to downregulate Ldb1/Lmo2-nucleated transcription complexes which normally function to enforce self-renewal in bone marrow hematopoietic progenitors.

E-cadherin is regulated by GATA-2 and marks the early commitment of mouse hematopoietic progenitors to the basophil and mast cell fates

E-cadherin is a calcium-dependent cell-cell adhesion molecule extensively studied for its involvement in tissue formation, epithelial cell behavior, and suppression of cancer. However, E-cadherin expression in the hematopoietic system has not been fully elucidated. Combining single-cell RNA-sequencing analyses and immunophenotyping, we revealed that progenitors expressing high levels of E-cadherin and contained within the granulocyte-monocyte progenitors (GMPs) fraction have an enriched capacity to differentiate into basophils and mast cells. We detected E-cadherin expression on committed progenitors before the expression of other reported markers of these lineages. We named such progenitors pro-BMPs (pro-basophil and mast cell progenitors). Using RNA sequencing, we observed transcriptional priming of pro-BMPs to the basophil and mast cell lineages. We also showed that GATA-2 directly regulates E-cadherin expression in the basophil and mast cell lineages, thus providing a mechanistic connection between the expression of this cell surface marker and the basophil and mast cell fate specification.

Embryonic erythropoiesis and hemoglobin switching require transcriptional repressor ETO2 to modulate chromatin organization

The underlying mechanism of transcriptional co-repressor ETO2 during early erythropoiesis and hemoglobin switching is unclear. We find that absence of ETO2 in mice interferes with down-regulation of PU.1 and GATA2 in the fetal liver, impeding a key step required for commitment to erythroid maturation. In human β-globin transgenic Eto2 null mice and in human CD34+ erythroid progenitor cells with reduced ETO2, loss of ETO2 results in ineffective silencing of embryonic/fetal globin gene expression, impeding hemoglobin switching during erythroid differentiation. ETO2 occupancy genome-wide occurs virtually exclusively at LDB1-complex binding sites in enhancers and ETO2 loss leads to increased enhancer activity and expression of target genes. ETO2 recruits the NuRD nucleosome remodeling and deacetylation complex to regulate histone acetylation and nucleosome occupancy in the β-globin locus control region and γ-globin gene. Loss of ETO2 elevates LDB1, MED1 and Pol II in the locus and facilitates fetal γ-globin/LCR looping and γ-globin transcription. Absence of the ETO2 hydrophobic heptad repeat region impairs ETO2-NuRD interaction and function in antagonizing γ-globin/LCR looping. Our results reveal a pivotal role for ETO2 in erythropoiesis and globin gene switching through its repressive role in the LDB1 complex, affecting the transcription factor and epigenetic environment and ultimately restructuring chromatin organization.

Selection Signatures Analysis Reveals Genes Associated with High-Altitude Adaptation in Tibetan Goats from Nagqu, Tibet

Simple Summary

In the process of domestication, goats have undergone long-term artificial and natural selection, leading to differences among goat breeds and leaving different selection traces on the genome. However, the genetic components underlying high-altitude adaptation remain largely unknown. Here, we genotyped four goat breeds using the Illumina Caprine 50K single nucleotide polymorphism (SNP) Chip. One highland breed (Tibetan goat) compared with three lowland breeds (Huanghuai goat, Taihang goat and Xinjiang goat) to identify the molecular basis of high-altitude adaptation. So, we investigated selection signatures using the di statistic of four goat breeds and some genes in Tibetan goats related to high-altitude adaptation were identified. In addition, q-PCR validated the gene expression level in Tibetan goats and Huanghuai goats. This information may be valuable for the study of the genetic uniqueness of Tibetan goats and increased understanding of the hypoxic adaptation mechanism of Tibetan goats on the plateau.

Abstract

Tibetan goat is an ancient breed, which inhabits the adverse conditions of the plateaus in China. To investigate the role of selection in shaping its genomes, we genotyped Tibetan goats (Nagqu Prefecture, above 4500 m) and three lowland populations (Xinjiang goats, Taihang goats and Huanghuai goats). The result of PCA, neighbor-joining (N-J) tree and model-based clustering showed that the genetic structure between the Tibetan goat and the three lowland populations has significant difference. As demonstrated by the di statistic, we found that some genes were related to the high-altitude adaptation of Tibetan goats. Functional analysis revealed that these genes were enriched in the VEGF (vascular endothelial growth factor) signaling pathway and melanoma, suggesting that nine genes (FGF2, EGFR, AKT1, PTEN, MITF, ENPEP, SIRT6, KDR, and CDC42) might have important roles in the high-altitude adaptation of Nagqu Tibetan goats. We also found that the LEPR gene was under the strongest selection (di value = 16.70), and it could induce upregulation of the hypoxic ventilatory response. In addition, five genes (LEPR, LDB1, EGFR, NOX4 and FGF2) with high di values were analyzed using q-PCR. Among them, we found that LEPR, LDB1 and FGF2 exhibited higher expression in the lungs of the Tibetan goats; LEPR, EGFR and LDB1 exhibited higher expression in the hearts of the Huanghuai goat. Our results suggest that LEPR, LDB1, EGFR and FGF2 genes may be related to the high-altitude adaptation of the goats. These findings improve our understanding of the selection of the high-altitude adaptability of the Nagqu Tibetan goats and provide new theoretical knowledge for the conservation and utilization of germplasm resources.

The transcriptional corepressor CBFA2T3 inhibits all-trans-retinoic acid-induced myeloid gene expression and differentiation in acute myeloid leukemia

CBFA2/RUNX1 partner transcriptional co-repressor 3 (CBFA2T3, also known as MTG16 or ETO2) is a myeloid translocation gene family protein that functions as a master transcriptional corepressor in hematopoiesis. Recently, it has been shown that CBFA2T3 maintains leukemia stem cell gene expression and promotes relapse in acute myeloid leukemia (AML). However, a role for CBFA2T3 in myeloid differentiation of AML has not been reported. Here, we show that CBFA2T3 represses retinoic acid receptor (RAR) target gene expression and inhibits all-trans-retinoic acid (ATRA)-induced myeloid differentiation of AML cells. ChIP-Seq revealed that CBFA2T3 targets the RARα/RXRα cistrome in U937 AML cells, predominantly at myeloid-specific enhancers associated with terminal differentiation. CRISPR/Cas9-mediated abrogation of CBFA2T3 resulted in spontaneous and ATRA-induced activation of myeloid-specific genes in a manner correlated with myeloid differentiation. Importantly, these effects were reversed by CBFA2T3 re-expression. Mechanistic studies showed that CBFA2T3 inhibits RAR target gene transcription by acting at an early step to regulate histone acetyltransferase recruitment, histone acetylation, and chromatin accessibility at RARα target sites, independently of the downstream, RAR-mediated steps of transcription. Finally, we validated the inhibitory effect of CBFA2T3 on RAR in multiple AML subtypes and patient samples. To our knowledge, this is the first study to show that CBFA2T3 down-regulation is both necessary and sufficient for enhancing ATRA-induced myeloid gene expression and differentiation of AML cells. Our findings suggest that CBFA2T3 can serve as a potential target for enhancing AML responsiveness to ATRA differentiation therapies.

LDB1 Enforces Stability on Direct and Indirect Oncoprotein Partners in Leukemia

The LMO2/LDB1 macromolecular complex is critical in hematopoietic stem and progenitor cell specification and in the development of acute leukemia. This complex is comprised of core subunits of LMO2 and LDB1 as well as single-stranded DNA-binding protein (SSBP) cofactors and DNA-binding basic helix-loop-helix (bHLH) and GATA transcription factors. We analyzed the steady-state abundance and kinetic stability of LMO2 and its partners via Halo protein tagging in conjunction with variant proteins deficient in binding their respective direct protein partners. We discovered a hierarchy of protein stabilities (with half-lives in descending order) as follows: LDB1 > SSBP > LMO2 > TAL1. Importantly, LDB1 is a remarkably stable protein that confers enhanced stability upon direct and indirect partners, thereby nucleating the formation of the multisubunit protein complex. The data imply that free subunits are more rapidly degraded than those incorporated within the LMO2/LDB1 complex. Our studies provided significant insights into LMO2/LDB1 macromolecular protein complex assembly and stability, which has implications for understanding its role in blood cell formation and for therapeutically targeting this complex in human leukemias.

A non-linear reverse-engineering method for inferring genetic regulatory networks

Hematopoiesis is a highly complex developmental process that produces various types of blood cells. This process is regulated by different genetic networks that control the proliferation, differentiation, and maturation of hematopoietic stem cells (HSCs). Although substantial progress has been made for understanding hematopoiesis, the detailed regulatory mechanisms for the fate determination of HSCs are still unraveled. In this study, we propose a novel approach to infer the detailed regulatory mechanisms. This work is designed to develop a mathematical framework that is able to realize nonlinear gene expression dynamics accurately. In particular, we intended to investigate the effect of possible protein heterodimers and/or synergistic effect in genetic regulation. This approach includes the Extended Forward Search Algorithm to infer network structure (top-down approach) and a non-linear mathematical model to infer dynamical property (bottom-up approach). Based on the published experimental data, we study two regulatory networks of 11 genes for regulating the erythrocyte differentiation pathway and the neutrophil differentiation pathway. The proposed algorithm is first applied to predict the network topologies among 11 genes and 55 non-linear terms which may be for heterodimers and/or synergistic effect. Then, the unknown model parameters are estimated by fitting simulations to the expression data of two different differentiation pathways. In addition, the edge deletion test is conducted to remove possible insignificant regulations from the inferred networks. Furthermore, the robustness property of the mathematical model is employed as an additional criterion to choose better network reconstruction results. Our simulation results successfully realized experimental data for two different differentiation pathways, which suggests that the proposed approach is an effective method to infer the topological structure and dynamic property of genetic regulations.

Rotational symmetry of the structured Chip/LDB-SSDP core module of the Wnt enhanceosome

The Chip/LIM-domain binding protein (LDB)-single-stranded DNA-binding protein (SSDP) (ChiLS) complex controls numerous cell-fate decisions in animal cells, by mediating transcription of developmental control genes via remote enhancers. ChiLS is recruited to these enhancers by lineage-specific LIM-domain proteins that bind to its Chip/LDB subunit. ChiLS recently emerged as the core module of the Wnt enhanceosome, a multiprotein complex that primes developmental control genes for timely Wnt responses. ChiLS binds to NPFxD motifs within Pygopus (Pygo) and the Osa/ARID1A subunit of the BAF chromatin remodeling complex, which could synergize with LIM proteins in tethering ChiLS to enhancers. Chip/LDB and SSDP both contain N-terminal dimerization domains that constitute the bulk of their structured cores. Here, we report the crystal structures of these dimerization domains, in part aided by DARPin chaperones. We conducted systematic surface scanning by structure-designed mutations, followed by in vitro and in vivo binding assays, to determine conserved surface residues required for binding between Chip/LDB, SSDP, and Pygo-NPFxD. Based on this, and on the 4:2 (SSDP-Chip/LDB) stoichiometry of ChiLS, we derive a highly constrained structural model for this complex, which adopts a rotationally symmetrical SSDP2-LDB2-SSDP2 architecture. Integrity of ChiLS is essential for Pygo binding, and our mutational analysis places the NPFxD pockets on either side of the Chip/LDB dimer, each flanked by an SSDP dimer. The symmetry and multivalency of ChiLS underpin its function as an enhancer module integrating Wnt signals with lineage-specific factors to operate context-dependent transcriptional switches that are pivotal for normal development and cancer.

The dynamic emergence of GATA1 complexes identified in in vitro embryonic stem cell differentiation and in vivo mouse fetal liver

GATA1 is an essential transcriptional regulator of myeloid hematopoietic differentiation towards red blood cells. During erythroid differentiation, GATA1 forms different complexes with other transcription factors such as LDB1, TAL1, E2A and LMO2 ("the LDB1 complex") or with FOG1. The functions of GATA1 complexes have been studied extensively in definitive erythroid differentiation; however, the temporal and spatial formation of these complexes during erythroid development is unknown. We applied proximity ligation assay (PLA) to detect, localize and quantify individual interactions during embryonic stem cell differentiation and in mouse fetal liver (FL) tissue. We show that GATA1/LDB1 interactions appear before the proerythroblast stage and increase in a subset of the CD71⁺/TER119^- cells to activate the terminal erythroid differentiation program in 12.5 day FL. Using Ldb1 and Gata1 knockdown FL cells, we studied the functional contribution of the GATA1/LDB1 complex during differentiation. This shows that the active LDB1 complex appears quite late at the proerythroblast stage of differentiation and confirms the power of PLA in studying the dynamic interaction of proteins in cell differentiation at the single cell level. We provide dynamic insight into the temporal and spatial formation of the GATA1 and LDB1 transcription factor complexes during hematopoietic development and differentiation.

GATA factor transcriptional activity: Insights from genome-wide binding profiles

The members of the GATA family of transcription factors have homologous zinc fingers and bind to similar sequence motifs. Recent advances in genome-wide technologies and the integration of bioinformatics data have led to a better understanding of how GATA factors regulate gene expression; GATA-factor-induced transcriptional and epigenetic changes have now been analyzed at unprecedented levels of detail. Here, we review the results of genome-wide studies of GATA factor occupancy in human and murine cell lines and primary cells (as determined by chromatin immunoprecipitation sequencing), and then discuss the molecular mechanisms underlying the mediation of transcriptional and epigenetic regulation by GATA factors.

Genome-Wide Maps of Transcription Regulatory Elements and Transcription Enhancers in Development and Disease

Gene expression is regulated by numerous elements including enhancers, insulators, transcription factors, and architectural proteins. Regions of DNA distal to the transcriptional start site, called enhancers, play a central role in the temporal and tissue-specific regulation of gene expression through RNA polymerase II. The identification of enhancers and other cis regulatory elements has largely been possible due to advances in next generation sequencing technologies. Enhancers regulate gene expression through chromatin loops mediated by architectural proteins such as YY1, CTCF, the cohesin complex, and LDB1. Additionally, enhancers can be transcribed to produce noncoding RNAs termed enhancer RNAs that likely participate in transcriptional regulation. The central role of enhancers in regulating gene expression implicates them in both normal physiology but also many disease states. The importance of enhancers is evident by the suggested role of SNPs, duplications, and other alterations of enhancer function in many diseases, ranging from cancer to atherosclerosis to chronic kidney disease. Although much progress has been made in recent years, the field of enhancer biology and our knowledge of the cis regulome remains a work in progress. This review will highlight recent seminal studies which demonstrate the role of enhancers in normal physiology and disease pathogenesis. © 2019 American Physiological Society. Compr Physiol 9:439-455, 2019.

Temporal autoregulation during human PU.1 locus SubTAD formation

Epigenetic control of gene expression occurs within discrete spatial chromosomal units called topologically associating domains (TADs), but the exact spatial requirements of most genes are unknown; this is of particular interest for genes involved in cancer. We therefore applied high-resolution chromosomal conformation capture sequencing to map the three-dimensional (3D) organization of the human locus encoding the key myeloid transcription factor PU.1 in healthy monocytes and acute myeloid leukemia (AML) cells. We identified a dynamic ∼75-kb unit (SubTAD) as the genomic region in which spatial interactions between PU.1 gene regulatory elements occur during myeloid differentiation and are interrupted in AML. Within this SubTAD, proper initiation of the spatial chromosomal interactions requires PU.1 autoregulation and recruitment of the chromatin-adaptor protein LDB1 (LIM domain-binding protein 1). However, once these spatial interactions have occurred, LDB1 stabilizes them independently of PU.1 autoregulation. Thus, our data support that PU.1 autoregulates its expression in a "hit-and-run" manner by initiating stable chromosomal loops that result in a transcriptionally active chromatin architecture.

Integrating Enhancer Mechanisms to Establish a Hierarchical Blood Development Program

Hematopoietic development requires the transcription factor GATA-2, and GATA-2 mutations cause diverse pathologies, including leukemia. GATA-2-regulated enhancers increase Gata2 expression in hematopoietic stem/progenitor cells and control hematopoiesis. The +9.5-kb enhancer activates transcription in endothelium and hematopoietic stem cells (HSCs), and its deletion abrogates HSC generation. The -77-kb enhancer activates transcription in myeloid progenitors, and its deletion impairs differentiation. Since +9.5^-/- embryos are HSC deficient, it was unclear whether the +9.5 functions in progenitors or if GATA-2 expression in progenitors solely requires -77. We further dissected the mechanisms using -77;+9.5 compound heterozygous (CH) mice. The embryonic lethal CH mutation depleted megakaryocyte-erythrocyte progenitors (MEPs). While the +9.5 suffices for HSC generation, the -77 and +9.5 must reside on one allele to induce MEPs. The -77 generated burst-forming unit-erythroid through the induction of GATA-1 and other GATA-2 targets. The enhancer circuits controlled signaling pathways that orchestrate a GATA factor-dependent blood development program.

LMO2 is required for TAL1 DNA binding activity and initiation of definitive haematopoiesis at the haemangioblast stage

LMO2 is a bridging factor within a DNA binding complex and is required for definitive haematopoiesis to occur. The developmental stage of the block in haematopoietic specification is not known. We show that Lmo2^−/− mouse embryonic stem cells differentiated to Flk-1⁺ haemangioblasts, but less efficiently to haemogenic endothelium, which only produced primitive haematopoietic progenitors. Genome-wide approaches indicated that LMO2 is required at the haemangioblast stage to position the TAL1/LMO2/LDB1 complex to regulatory elements that are important for the establishment of the haematopoietic developmental program. In the absence of LMO2, the target site recognition of TAL1 is impaired. The lack of LMO2 resulted in altered gene expression levels already at the haemangioblast stage, with transcription factor genes accounting for ∼15% of affected genes. Comparison of Lmo2^−/− with Tal1^−/− Flk-1⁺ cells further showed that TAL1 was required to initiate or sustain Lmo2 expression.

LDB2 inhibits proliferation and migration in liver cancer cells by abrogating HEY1 expression

Hepatocellular carcinoma (HCC) was one of the most common cancers around the world, has very low 5-year survival rate. However, the mechanism of HCC occurrence and development is largely unknown. LDB2 belongs to the LIM-domain binding family and functions as an adaptor for transcriptional regulation. Here we found that LDB2 is downregulated in HCC samples. LDB2 has the ability to inhibit proliferation and migration of hepatocarcinoma cells. We found that the proliferation and migration abilities in HCC sample cells were impaired after LDB2 overexpression and vice versa. In mechanism, we found that LDB2 can recruit BRD7 to HEY1 promoter and then block its expression. HEY1 whose expression is upregulated in HCC acts as an oncogene. In brief, our research reveals a new regulatory mechanism for hepatocarcinoma cell proliferation and migration.

GATA Factor-Regulated Samd14 Enhancer Confers Red Blood Cell Regeneration and Survival in Severe Anemia

An enhancer with amalgamated E-box and GATA motifs (+9.5) controls expression of the regulator of hematopoiesis GATA-2. While similar GATA-2-occupied elements are common in the genome, occupancy does not predict function, and GATA-2-dependent genetic networks are incompletely defined. A "+9.5-like" element resides in an intron of Samd14 (Samd14-Enh) encoding a sterile alpha motif (SAM) domain protein. Deletion of Samd14-Enh in mice strongly decreased Samd14 expression in bone marrow and spleen. Although steady-state hematopoiesis was normal, Samd14-Enh^-/- mice died in response to severe anemia. Samd14-Enh stimulated stem cell factor/c-Kit signaling, which promotes erythrocyte regeneration. Anemia activated Samd14-Enh by inducing enhancer components and enhancer chromatin accessibility. Thus, a GATA-2/anemia-regulated enhancer controls expression of an SAM domain protein that confers survival in anemia. We propose that Samd14-Enh and an ensemble of anemia-responsive enhancers are essential for erythrocyte regeneration in stress erythropoiesis, a vital process in pathologies, including β-thalassemia, myelodysplastic syndrome, and viral infection.

GATA3 Abundance Is a Critical Determinant of T Cell Receptor β Allelic Exclusion

Allelic exclusion describes the essential immunological process by which feedback repression of sequential DNA rearrangements ensures that only one autosome expresses a functional T or B cell receptor. In wild-type mammals, approximately 60% of cells have recombined the DNA of one T cell receptor β (TCRβ) V-to-DJ-joined allele in a functional configuration, while the second allele has recombined only the DJ sequences; the other 40% of cells have recombined the V to the DJ segments on both alleles, with only one of the two alleles predicting a functional TCRβ protein. Here we report that the transgenic overexpression of GATA3 leads predominantly to biallelic TCRβ gene (Tcrb) recombination. We also found that wild-type immature thymocytes can be separated into distinct populations based on intracellular GATA3 expression and that GATA3^LO cells had almost exclusively recombined only one Tcrb locus (that predicted a functional receptor sequence), while GATA3^HI cells had uniformly recombined both Tcrb alleles (one predicting a functional and the other predicting a nonfunctional rearrangement). These data show that GATA3 abundance regulates the recombination propensity at the Tcrb locus and provide new mechanistic insight into the historic immunological conundrum for how Tcrb allelic exclusion is mediated.

LIM domain-binding 1 maintains the terminally differentiated state of pancreatic β cells

The recognition of β cell dedifferentiation in type 2 diabetes raises the translational relevance of mechanisms that direct and maintain β cell identity. LIM domain-binding protein 1 (LDB1) nucleates multimeric transcriptional complexes and establishes promoter-enhancer looping, thereby directing fate assignment and maturation of progenitor populations. Many terminally differentiated endocrine cell types, however, remain enriched for LDB1, but its role is unknown. Here, we have demonstrated a requirement for LDB1 in maintaining the terminally differentiated status of pancreatic β cells. Inducible ablation of LDB1 in mature β cells impaired insulin secretion and glucose homeostasis. Transcriptomic analysis of LDB1-depleted β cells revealed the collapse of the terminally differentiated gene program, indicated by a loss of β cell identity genes and induction of the endocrine progenitor factor neurogenin 3 (NEUROG3). Lineage tracing confirmed that LDB1-depleted, insulin-negative β cells express NEUROG3 but do not adopt alternate endocrine cell fates. In primary mouse islets, LDB1 and its LIM homeodomain-binding partner islet 1 (ISL1) were coenriched at chromatin sites occupied by pancreatic and duodenal homeobox 1 (PDX1), NK6 homeobox 1 (NKX6.1), forkhead box A2 (FOXA2), and NK2 homeobox 2 (NKX2.2) - factors that co-occupy active enhancers in 3D chromatin domains in human islets. Indeed, LDB1 was enriched at active enhancers in human islets. Thus, LDB1 maintains the terminally differentiated state of β cells and is a component of active enhancers in both murine and human islets.

LIM-Only Protein 4 (LMO4) and LIM Domain Binding Protein 1 (LDB1) Promote Growth and Metastasis of Human Head and Neck Cancer (LMO4 and LDB1 in Head and Neck Cancer)

Squamous cell carcinoma of the head and neck (HNSCC) accounts for more than 300,000 deaths worldwide per year as a consequence of tumor cell invasion of adjacent structures or metastasis. LIM-only protein 4 (LMO4) and LIM-domain binding protein 1 (LDB1), two directly interacting transcriptional adaptors that have important roles in normal epithelial cell differentiation, have been associated with increased metastasis, decreased differentiation, and shortened survival in carcinoma of the breast. Here, we implicate two LDB1-binding proteins, single-stranded binding protein 2 (SSBP2) and 3 (SSBP3), in controlling LMO4 and LDB1 protein abundance in HNSCC and in regulating specific tumor cell functions in this disease. First, we found that the relative abundance of LMO4, LDB1, and the two SSBPs correlated very significantly in a panel of human HNSCC cell lines. Second, expression of these proteins in tumor primaries and lymph nodes involved by metastasis were concordant in 3 of 3 sets of tissue. Third, using a Matrigel invasion and organotypic reconstruct assay, CRISPR/Cas9-mediated deletion of LDB1 in the VU-SCC-1729 cell line, which is highly invasive of basement membrane and cellular monolayers, reduced tumor cell invasiveness and migration, as well as proliferation on tissue culture plastic. Finally, inactivation of the LDB1 gene in these cells decreased growth and vascularization of xenografted human tumor cells in vivo. These data show that LMO4, LDB1, and SSBP2 and/or SSBP3 regulate metastasis, proliferation, and angiogenesis in HNSCC and provide the first evidence that SSBPs control LMO4 and LDB1 protein abundance in a cancer context.

Cadmium modulates hematopoietic stem and progenitor cells and skews toward myelopoiesis in mice

The heavy metal cadmium (Cd) is known to modulate immunity and cause osteoporosis. However, how Cd influences on hematopoiesis remain largely unknown. Herein, we show that wild-type C57BL/6 (B6) mice exposed to Cd for 3months had expanded bone marrow (BM) populations of long-term hematopoietic stem cells (LT-HSCs), common myeloid progenitors (CMPs) and granulocyte-macrophage progenitors (GMPs), while having reduced populations of multipotent progenitors (MPPs) and common lymphoid progenitors (CLPs). A competitive mixed BM transplantation assay indicates that BM from Cd-treated mice had impaired LT-HSC ability to differentiate into mature cells. In accordance with increased myeloid progenitors and decreased lymphoid progenitors, the BM and spleens of Cd-treated mice had more monocytes and/or neutrophils and fewer B cells and T cells. Cd impaired the ability of the non-hematopoietic system to support LT-HSCs, in that lethally irradiated Cd-treated recipients transplanted with normal BM cells had reduced LT-HSCs after the hematopoietic system was fully reconstituted. This is consistent with reduced osteoblasts, a known critical component for HSC niche, observed in Cd-treated mice. Conversely, lethally irradiated control recipients transplanted with BM cells from Cd-treated mice had normal LT-HSC reconstitution. Furthermore, both control mice and Cd-treated mice that received Alendronate, a clinical drug used for treating osteoporosis, had BM increases of LT-HSCs. Thus, the results suggest Cd increase of LT-HSCs is due to effects on HSCs and not on osteoblasts, although, Cd causes osteoblast reduction and impaired niche function for maintaining HSCs. Furthermore, Cd skews HSCs toward myelopoiesis.

The stage-dependent roles of Ldb1 and functional redundancy with Ldb2 in mammalian retinogenesis

The Lim domain-binding proteins are key co-factor proteins that assemble with LIM domains of the LMO/LIM-HD family to form functional complexes that regulate cell proliferation and differentiation. Using conditional mutagenesis and comparative phenotypic analysis, we analyze the function of Ldb1 and Ldb2 in mouse retinal development, and demonstrate overlapping and specific functions of both proteins. Ldb1 interacts with Lhx2 in the embryonic retina and both Ldb1 and Ldb2 play a key role in maintaining the pool of retinal progenitor cells. This is accomplished by controlling the expression of the Vsx2 and Rax, and components of the Notch and Hedgehog signaling pathways. Furthermore, the Ldb1/Ldb2-mediated complex is essential for generation of early-born photoreceptors through the regulation of Rax and Crx. Finally, we demonstrate functional redundancy between Ldb1 and Ldb2. Ldb1 can fully compensate the loss of Ldb2 during all phases of retinal development, whereas Ldb2 alone is sufficient to sustain activity of Lhx2 in both early- and late-stage RPCs and in Müller glia. By contrast, loss of Ldb1 disrupts activity of the LIM domain factors in neuronal precursors. An intricate regulatory network exists that is mediated by Ldb1 and Ldb2, and promotes RPC proliferation and multipotency; it also controls specification of mammalian retina cells.

BTR: training asynchronous Boolean models using single-cell expression data

Background

Rapid technological innovation for the generation of single-cell genomics data presents new challenges and opportunities for bioinformatics analysis. One such area lies in the development of new ways to train gene regulatory networks. The use of single-cell expression profiling technique allows the profiling of the expression states of hundreds of cells, but these expression states are typically noisier due to the presence of technical artefacts such as drop-outs. While many algorithms exist to infer a gene regulatory network, very few of them are able to harness the extra expression states present in single-cell expression data without getting adversely affected by the substantial technical noise present.

Results

Here we introduce BTR, an algorithm for training asynchronous Boolean models with single-cell expression data using a novel Boolean state space scoring function. BTR is capable of refining existing Boolean models and reconstructing new Boolean models by improving the match between model prediction and expression data. We demonstrate that the Boolean scoring function performed favourably against the BIC scoring function for Bayesian networks. In addition, we show that BTR outperforms many other network inference algorithms in both bulk and single-cell synthetic expression data. Lastly, we introduce two case studies, in which we use BTR to improve published Boolean models in order to generate potentially new biological insights.

Conclusions

BTR provides a novel way to refine or reconstruct Boolean models using single-cell expression data. Boolean model is particularly useful for network reconstruction using single-cell data because it is more robust to the effect of drop-outs. In addition, BTR does not assume any relationship in the expression states among cells, it is useful for reconstructing a gene regulatory network with as few assumptions as possible. Given the simplicity of Boolean models and the rapid adoption of single-cell genomics by biologists, BTR has the potential to make an impact across many fields of biomedical research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1235-y) contains supplementary material, which is available to authorized users.

Transcription factor defects causing platelet disorders

Recent years have seen increasing recognition of a subgroup of inherited platelet function disorders which are due to defects in transcription factors that are required to regulate megakaryopoiesis and platelet production. Thus, germline mutations in the genes encoding the haematopoietic transcription factors RUNX1, GATA-1, FLI1, GFI1b and ETV6 have been associated with both quantitative and qualitative platelet abnormalities, and variable bleeding symptoms in the affected patients. Some of the transcription factor defects are also associated with an increased predisposition to haematologic malignancies (RUNX1, ETV6), abnormal erythropoiesis (GATA-1, GFI1b, ETV6) and immune dysfunction (FLI1). The persistence of MYH10 expression in platelets is a surrogate marker for FLI1 and RUNX1 defects. Characterisation of the transcription factor defects that give rise to platelet function disorders, and of the genes that are differentially regulated as a result, are yielding insights into the roles of these genes in platelet formation and function.

Role of the RB-Interacting Proteins in Stem Cell Biology

Human retinoblastoma gene RB1 is the first tumor suppressor gene (TSG) isolated by positional cloning in 1986. RB is extensively studied for its ability to regulate cell cycle by binding to E2F1 and inhibiting the transcriptional activity of the latter. In human embryonic stem cells (ESCs), only a minute trace of RB is found in complex with E2F1. Increased activity of RB triggers differentiation, cell cycle arrest, and cell death. On the other hand, inactivation of the entire RB family (RB1, RBL1, and RBL2) in human ESC induces G2/M arrest and cell death. These observations indicate that both loss and overactivity of RB could be lethal for the stemness of cells. A question arises why inactive RB is required for the survival and stemness of cells? To shed some light on this question, we analyzed the RB-binding proteins. In this review we have focused on 27 RB-binding partners that may have potential roles in different aspects of stem cell biology.

Ldb1 Is Essential for the Development of Isthmic Organizer and Midbrain Dopaminergic Neurons

LIM domain-binding protein 1 (Ldb1) is a nuclear cofactor that interacts with LIM homeodomain proteins to form multiprotein complexes that are important for transcription regulation. Ldb1 has been shown to play essential roles in various processes during mouse embryogenesis. To determine the role of Ldb1 in mid- and hindbrain development, we have generated a conditional mutant with a specific deletion of the Ldb1 in the Engrailed-1-expressing region of the developing mid- and hindbrain. Our study showed that the deletion impaired the expression of signaling molecules, such as fibroblast growth factor 8 (FGF8) and Wnt1, in the isthmic organizer and the expression of Shh in the ventral midbrain. The midbrain and the cerebellum were severely reduced in size, and the midbrain dopaminergic (mDA) neurons were missing in the mutant. These defects are identical to the phenotype that has been observed previously in mice with a deletion of the LIM homeodomain gene Lmx1b. Our results thus provide genetic evidence supporting that Ldb1 and Lmx1b function cooperatively to regulate mid- and hindbrain development. In addition, we found that mouse embryonic stem cells lacking Ldb1 failed to generate several types of differentiated neurons, including the mDA neurons, serotonergic neurons, cholinergic neurons, and olfactory bulb neurons, indicating an essential cell-autonomous role for Ldb1 in the development of these neurons.

Cofactors of LIM Domains Associate with Estrogen Receptor α to Regulate the Expression of Noncoding RNA H19 and Corneal Epithelial Progenitor Cell Function

Cofactors of LIM domain proteins, CLIM1 and CLIM2, are widely expressed transcriptional cofactors that are recruited to gene regulatory regions by DNA-binding proteins, including LIM domain transcription factors. In the cornea, epithelium-specific expression of a dominant negative (DN) CLIM under the keratin 14 (K14) promoter causes blistering, wounding, inflammation, epithelial hyperplasia, and neovascularization followed by epithelial thinning and subsequent epidermal-like differentiation of the corneal epithelium. The defects in corneal epithelial differentiation and cell fate determination suggest that CLIM may regulate corneal progenitor cells and the transition to differentiation. Consistent with this notion, the K14-DN-Clim corneal epithelium first exhibits increased proliferation followed by fewer progenitor cells with decreased proliferative potential. In vivo ChIP-sequencing experiments with corneal epithelium show that CLIM binds to and regulates numerous genes involved in cell adhesion and proliferation, including limbally enriched genes. Intriguingly, CLIM associates primarily with non-LIM homeodomain motifs in corneal epithelial cells, including that of estrogen receptor α. Among CLIM targets is the noncoding RNA H19 whose deregulation is associated with Silver-Russell and Beckwith-Wiedemann syndromes. We demonstrate here that H19 negatively regulates corneal epithelial proliferation. In addition to cell cycle regulators, H19 affects the expression of multiple cell adhesion genes. CLIM interacts with estrogen receptor α at the H19 locus, potentially explaining the higher expression of H19 in female than male corneas. Together, our results demonstrate an important role for CLIM in regulating the proliferative potential of corneal epithelial progenitors and identify CLIM downstream target H19 as a regulator of corneal epithelial proliferation and adhesion.

GATA Factor-G-Protein-Coupled Receptor Circuit Suppresses Hematopoiesis

Hematopoietic stem cells (HSCs) originate from hemogenic endothelium within the aorta-gonad-mesonephros (AGM) region of the mammalian embryo. The relationship between genetic circuits controlling stem cell genesis and multi-potency is not understood. A Gata2 cis element (+9.5) enhances Gata2 expression in the AGM and induces the endothelial to HSC transition. We demonstrated that GATA-2 rescued hematopoiesis in +9.5^−/− AGMs. As G-protein-coupled receptors (GPCRs) are the most common targets for FDA-approved drugs, we analyzed the GPCR gene ensemble to identify GATA-2-regulated GPCRs. Of the 20 GATA-2-activated GPCR genes, four were GATA-1-activated, and only Gpr65 expression resembled Gata2. Contrasting with the paradigm in which GATA-2-activated genes promote hematopoietic stem and progenitor cell genesis/function, our mouse and zebrafish studies indicated that GPR65 suppressed hematopoiesis. GPR65 established repressive chromatin at the +9.5 site, restricted occupancy by the activator Scl/TAL1, and repressed Gata2 transcription. Thus, a Gata2 cis element creates a GATA-2-GPCR circuit that limits positive regulators that promote hematopoiesis.

•

GATA-2 rescues +9.5^−/− AGM hematopoietic activity

•

GATA-2 upregulates Gpr65, which encodes a negative regulator of hematopoiesis

•

GPR65 suppresses hematopoiesis by repressing Gata2 expression

•

GPR65 represses Gata2 expression by increasing H4K20me1, restricting Scl/TAL1 occupancy

Chromatin looping as a target for altering erythroid gene expression

The β-hemoglobinopathies are the most common monogenic disorders in humans, with symptoms arising after birth when the fetal γ-globin genes are silenced and the adult β-globin gene is activated. There is a growing appreciation that genome organization and the folding of chromosomes are key determinants of gene transcription. Underlying this function is the activity of transcriptional enhancers that increase the transcription of target genes over long linear distances. To accomplish this, enhancers engage in close physical contact with target promoters through chromosome folding or looping that is orchestrated by protein complexes that bind to both sites and stabilize their interaction. We find that enhancer activity can be redirected with concomitant changes in gene transcription. Both targeting the β-globin locus control region (LCR) to the γ-globin gene in adult erythroid cells by tethering and epigenetic unmasking of a silenced γ-globin gene lead to increased frequency of LCR/γ-globin contacts and reduced LCR/β-globin contacts. The outcome of these manipulations is robust, pancellular γ-globin transcription activation with a concomitant reduction in β-globin transcription. These examples show that chromosome looping may be considered a therapeutic target for gene activation in β-thalassemia and sickle cell disease.

3off2: A network reconstruction algorithm based on 2-point and 3-point information statistics

Background

The reconstruction of reliable graphical models from observational data is important in bioinformatics and other computational fields applying network reconstruction methods to large, yet finite datasets. The main network reconstruction approaches are either based on Bayesian scores, which enable the ranking of alternative Bayesian networks, or rely on the identification of structural independencies, which correspond to missing edges in the underlying network. Bayesian inference methods typically require heuristic search strategies, such as hill-climbing algorithms, to sample the super-exponential space of possible networks. By contrast, constraint-based methods, such as the PC and IC algorithms, are expected to run in polynomial time on sparse underlying graphs, provided that a correct list of conditional independencies is available. Yet, in practice, conditional independencies need to be ascertained from the available observational data, based on adjustable statistical significance levels, and are not robust to sampling noise from finite datasets.

Results

We propose a more robust approach to reconstruct graphical models from finite datasets. It combines constraint-based and Bayesian approaches to infer structural independencies based on the ranking of their most likely contributing nodes. In a nutshell, this local optimization scheme and corresponding 3off2 algorithm iteratively “take off” the most likely conditional 3-point information from the 2-point (mutual) information between each pair of nodes. Conditional independencies are thus derived by progressively collecting the most significant indirect contributions to all pairwise mutual information. The resulting network skeleton is then partially directed by orienting and propagating edge directions, based on the sign and magnitude of the conditional 3-point information of unshielded triples. The approach is shown to outperform both constraint-based and Bayesian inference methods on a range of benchmark networks. The 3off2 approach is then applied to the reconstruction of the hematopoiesis regulation network based on recent single cell expression data and is found to retrieve more experimentally ascertained regulations between transcription factors than with other available methods.

Conclusions

The novel information-theoretic approach and corresponding 3off2 algorithm combine constraint-based and Bayesian inference methods to reliably reconstruct graphical models, despite inherent sampling noise in finite datasets. In particular, experimentally verified interactions as well as novel predicted regulations are established on the hematopoiesis regulatory networks based on single cell expression data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0856-x) contains supplementary material, which is available to authorized users.

The LMO2 oncogene regulates DNA replication in hematopoietic cells

Oncogenic transcription factors are commonly activated in acute leukemias and subvert normal gene expression networks to reprogram hematopoietic progenitors into preleukemic stem cells, as exemplified by LIM-only 2 (LMO2) in T-cell acute lymphoblastic leukemia (T-ALL). Whether or not these oncoproteins interfere with other DNA-dependent processes is largely unexplored. Here, we show that LMO2 is recruited to DNA replication origins by interaction with three essential replication enzymes: DNA polymerase delta (POLD1), DNA primase (PRIM1), and minichromosome 6 (MCM6). Furthermore, tethering LMO2 to synthetic DNA sequences is sufficient to transform these sequences into origins of replication. We next addressed the importance of LMO2 in erythroid and thymocyte development, two lineages in which cell cycle and differentiation are tightly coordinated. Lowering LMO2 levels in erythroid progenitors delays G1-S progression and arrests erythropoietin-dependent cell growth while favoring terminal differentiation. Conversely, ectopic expression in thymocytes induces DNA replication and drives these cells into cell cycle, causing differentiation blockade. Our results define a novel role for LMO2 in directly promoting DNA synthesis and G1-S progression.

The Hematopoietic Stem and Progenitor Cell Cistrome: GATA Factor-Dependent cis-Regulatory Mechanisms

Transcriptional regulators mediate the genesis and function of the hematopoietic system by binding complex ensembles of cis-regulatory elements to establish genetic networks. While thousands to millions of any given cis-element resides in a genome, how transcriptional regulators select these sites and how site attributes dictate functional output is not well understood. An instructive system to address this problem involves the GATA family of transcription factors that control vital developmental and physiological processes and are linked to multiple human pathologies. Although GATA factors bind DNA motifs harboring the sequence GATA, only a very small subset of these abundant motifs are occupied in genomes. Mechanistic studies revealed a unique configuration of a GATA factor-regulated cis-element consisting of an E-box and a downstream GATA motif separated by a short DNA spacer. GATA-1- or GATA-2-containing multiprotein complexes at these composite elements control transcription of genes critical for hematopoietic stem cell emergence in the mammalian embryo, hematopoietic progenitor cell regulation, and erythroid cell maturation. Other constituents of the complex include the basic helix-loop-loop transcription factor Scl/TAL1, its heterodimeric partner E2A, and the Lim domain proteins LMO2 and LDB1. This chapter reviews the structure/function of E-box-GATA composite cis-elements, which collectively constitute an important sector of the hematopoietic stem and progenitor cell cistrome.

LMO2 Oncoprotein Stability in T-Cell Leukemia Requires Direct LDB1 Binding

LMO2 is a component of multisubunit DNA-binding transcription factor complexes that regulate gene expression in hematopoietic stem and progenitor cell development. Enforced expression of LMO2 causes leukemia by inducing hematopoietic stem cell-like features in T-cell progenitor cells, but the biochemical mechanisms of LMO2 function have not been fully elucidated. In this study, we systematically dissected the LMO2/LDB1-binding interface to investigate the role of this interaction in T-cell leukemia. Alanine scanning mutagenesis of the LIM interaction domain of LDB1 revealed a discrete motif, R(320)LITR, required for LMO2 binding. Most strikingly, coexpression of full-length, wild-type LDB1 increased LMO2 steady-state abundance, whereas coexpression of mutant proteins deficient in LMO2 binding compromised LMO2 stability. These mutant LDB1 proteins also exerted dominant negative effects on growth and transcription in diverse leukemic cell lines. Mass spectrometric analysis of LDB1 binding partners in leukemic lines supports the notion that LMO2/LDB1 function in leukemia occurs in the context of multisubunit complexes, which also protect the LMO2 oncoprotein from degradation. Collectively, these data suggest that the assembly of LMO2 into complexes, via direct LDB1 interaction, is a potential molecular target that could be exploited in LMO2-driven leukemias resistant to existing chemotherapy regimens.

A Systems Approach Identifies Essential FOXO3 Functions at Key Steps of Terminal Erythropoiesis

Circulating red blood cells (RBCs) are essential for tissue oxygenation and homeostasis. Defective terminal erythropoiesis contributes to decreased generation of RBCs in many disorders. Specifically, ineffective nuclear expulsion (enucleation) during terminal maturation is an obstacle to therapeutic RBC production in vitro. To obtain mechanistic insights into terminal erythropoiesis we focused on FOXO3, a transcription factor implicated in erythroid disorders. Using an integrated computational and experimental systems biology approach, we show that FOXO3 is essential for the correct temporal gene expression during terminal erythropoiesis. We demonstrate that the FOXO3-dependent genetic network has critical physiological functions at key steps of terminal erythropoiesis including enucleation and mitochondrial clearance processes. FOXO3 loss deregulated transcription of genes implicated in cell polarity, nucleosome assembly and DNA packaging-related processes and compromised erythroid enucleation. Using high-resolution confocal microscopy and imaging flow cytometry we show that cell polarization is impaired leading to multilobulated Foxo3^-/- erythroblasts defective in nuclear expulsion. Ectopic FOXO3 expression rescued Foxo3^-/- erythroblast enucleation-related gene transcription, enucleation defects and terminal maturation. Remarkably, FOXO3 ectopic expression increased wild type erythroblast maturation and enucleation suggesting that enhancing FOXO3 activity may improve RBCs production. Altogether these studies uncover FOXO3 as a novel regulator of erythroblast enucleation and terminal maturation suggesting FOXO3 modulation might be therapeutic in disorders with defective erythroid maturation.

Red blood cells (RBCs) are highly specialized cells that transport oxygen throughout the body and are essential for survival. However, RBCs have a limited lifespan and need to be replenished continuously by stem cells in the bone marrow. Mammalian RBCs are unique in that in order to fully mature they exclude their nucleus and other organelles. Mechanisms involved in these processes are not well understood at the molecular level. Defects in any of the these processes may lead to red blood cell defects, a decreased capacity to transport oxygen and/or a block in red blood cell production in vitro. Therefore, understanding how these processes are regulated at the molecular level can lead to promising new therapies for red blood cell defects and improved methods of generating red blood cells in a dish. Here, using an integrated computational and experimental biology approach, we found that the nuclear factor FOXO3 is a crucial regulator of red blood cell production by coordinating the expression of many of the genes specific for terminal maturation of red blood cells. Furthermore we found that FOXO3 can even increase the production of normal red blood cells in culture raising the possibility that enhancing FOXO3 may have a therapeutic use. Our studies identify FOXO3 as a novel regulator of RBC enucleation and terminal erythropoiesis.

An ancient Pygo-dependent Wnt enhanceosome integrated by Chip/LDB-SSDP

TCF/LEF factors are ancient context-dependent enhancer-binding proteins that are activated by β-catenin following Wnt signaling. They control embryonic development and adult stem cell compartments, and their dysregulation often causes cancer. β-catenin-dependent transcription relies on the NPF motif of Pygo proteins. Here, we use a proteomics approach to discover the Chip/LDB-SSDP (ChiLS) complex as the ligand specifically binding to NPF. ChiLS also recognizes NPF motifs in other nuclear factors including Runt/RUNX2 and Drosophila ARID1, and binds to Groucho/TLE. Studies of Wnt-responsive dTCF enhancers in the Drosophila embryonic midgut indicate how these factors interact to form the Wnt enhanceosome, primed for Wnt responses by Pygo. Together with previous evidence, our study indicates that ChiLS confers context-dependence on TCF/LEF by integrating multiple inputs from lineage and signal-responsive factors, including enhanceosome switch-off by Notch. Its pivotal function in embryos and stem cells explain why its integrity is crucial in the avoidance of cancer.

A novel TGFβ modulator that uncouples R-Smad/I-Smad-mediated negative feedback from R-Smad/ligand-driven positive feedback

As some of the most widely utilised intercellular signalling molecules, transforming growth factor β (TGFβ) superfamily members play critical roles in normal development and become disrupted in human disease. Establishing appropriate levels of TGFβ signalling involves positive and negative feedback, which are coupled and driven by the same signal transduction components (R-Smad transcription factor complexes), but whether and how the regulation of the two can be distinguished are unknown. Genome-wide comparison of published ChIP-seq datasets suggests that LIM domain binding proteins (Ldbs) co-localise with R-Smads at a substantial subset of R-Smad target genes including the locus of inhibitory Smad7 (I-Smad7), which mediates negative feedback for TGFβ signalling. We present evidence suggesting that zebrafish Ldb2a binds and directly activates the I-Smad7 gene, whereas it binds and represses the ligand gene, Squint (Sqt), which drives positive feedback. Thus, the fine tuning of TGFβ signalling derives from positive and negative control by Ldb2a. Expression of ldb2a is itself activated by TGFβ signals, suggesting potential feed-forward loops that might delay the negative input of Ldb2a to the positive feedback, as well as the positive input of Ldb2a to the negative feedback. In this way, precise gene expression control by Ldb2a enables an initial build-up of signalling via a fully active positive feedback in the absence of buffering by the negative feedback. In Ldb2a-deficient zebrafish embryos, homeostasis of TGFβ signalling is perturbed and signalling is stably enhanced, giving rise to excess mesoderm and endoderm, an effect that can be rescued by reducing signalling by the TGFβ family members, Nodal and BMP. Thus, Ldb2a is critical to the homeostatic control of TGFβ signalling and thereby embryonic patterning.

Early stages of induction of anterior head ectodermal properties in Xenopus embryos are mediated by transcriptional cofactor ldb1

BACKGROUND

Specific molecules involved in early inductive signaling from anterior neural tissue to the placodal ectoderm to establish a lens-forming bias, as well as their regulatory factors, remain largely unknown. In this study, we sought to identify and characterize these molecules.

RESULTS

Using an expression cloning strategy to isolate genes with lens-inducing activity, we identified the transcriptional cofactor ldb1. This, together with evidence for its nuclear dependence, suggests its role as a regulatory factor, not a direct signaling molecule. We propose that ldb1 mediates induction of early lens genes in our functional assay by transcriptional activation of lens-inducing signals. Gain-of-function assays demonstrate that the inductive activity of the anterior neural plate on head ectodermal structures can be augmented by ldb1. Loss-of-function assays show that knockdown of ldb1 leads to decreased expression of early lens and retinal markers and subsequently to defects in eye development.

CONCLUSIONS

The functional cloning, expression pattern, overexpression, and knockdown data show that an ldb1-regulated mechanism acts as an early signal for Xenopus lens induction.

Requirement for ssbp2 in hematopoietic stem cell maintenance and stress response

Transcriptional mechanisms governing hematopoietic stem cell (HSC) quiescence, self-renewal, and differentiation are not fully understood. Sequence-specific ssDNA-binding protein 2 (SSBP2) is a candidate acute myelogenous leukemia (AML) suppressor gene located at chromosome 5q14. SSBP2 binds the transcriptional adaptor protein Lim domain-binding protein 1 (LDB1) and enhances LDB1 stability to regulate gene expression. Notably, Ldb1 is essential for HSC specification during early development and maintenance in adults. We previously reported shortened lifespan and greater susceptibility to B cell lymphomas and carcinomas in Ssbp2(-/-) mice. However, whether Ssbp2 plays a regulatory role in normal HSC function and leukemogenesis is unknown. In this study, we provide several lines of evidence to demonstrate a requirement for Ssbp2 in the function and transcriptional program of hematopoietic stem and progenitor cells (HSPCs) in vivo. We found that hematopoietic tissues were hypoplastic in Ssbp2(-/-) mice, and the frequency of lymphoid-primed multipotent progenitor cells in bone marrow was reduced. Other significant features of these mice were delayed recovery from 5-fluorouracil treatment and diminished multilineage reconstitution in lethally irradiated bone marrow recipients. Dramatic reduction of Notch1 transcripts and increased expression of transcripts encoding the transcription factor E2a and its downstream target Cdkn1a also distinguished Ssbp2(-/-) HSPCs from wild-type HSPCs. Finally, a tendency toward coordinated expression of SSBP2 and the AML suppressor NOTCH1 in a subset of the Cancer Genome Atlas AML cases suggested a role for SSBP2 in AML pathogenesis. Collectively, our results uncovered a critical regulatory function for SSBP2 in HSPC gene expression and function.