Spi-1/PU.1 oncogene accelerates DNA replication fork elongation and promotes genetic instability in the absence of DNA breakage

A multi-omics approach for biomarker discovery in neuroblastoma: a network-based framework

Neuroblastoma (NB) is one of the leading causes of cancer-associated death in children. MYCN amplification is a prominent genetic marker for NB, and its targeting to halt NB progression is difficult to achieve. Therefore, an in-depth understanding of the molecular interactome of NB is needed to improve treatment outcomes. Analysis of NB multi-omics unravels valuable insight into the interplay between MYCN transcriptional and miRNA post-transcriptional modulation. Moreover, it aids in the identification of various miRNAs that participate in NB development and progression. This study proposes an integrated computational framework with three levels of high-throughput NB data (mRNA-seq, miRNA-seq, and methylation array). Similarity Network Fusion (SNF) and ranked SNF methods were utilized to identify essential genes and miRNAs. The specified genes included both miRNA-target genes and transcription factors (TFs). The interactions between TFs and miRNAs and between miRNAs and their target genes were retrieved where a regulatory network was developed. Finally, an interaction network-based analysis was performed to identify candidate biomarkers. The candidate biomarkers were further analyzed for their potential use in prognosis and diagnosis. The candidate biomarkers included three TFs and seven miRNAs. Four biomarkers have been previously studied and tested in NB, while the remaining identified biomarkers have known roles in other types of cancer. Although the specific molecular role is yet to be addressed, most identified biomarkers possess evidence of involvement in NB tumorigenesis. Analyzing cellular interactome to identify potential biomarkers is a promising approach that can contribute to optimizing efficient therapeutic regimens to target NB vulnerabilities.

SPI1 activates TGF-β1/PI3K/Akt signaling through transcriptional upregulation of FKBP12 to support the mesenchymal phenotype of glioma stem cells

Glioma stem cells (GSCs) exhibit diverse molecular subtypes with the mesenchymal (MES) population representing the most malignant variant. The oncogenic potential of Salmonella pathogenicity island 1 (SPI1), an oncogenic transcription factor, has been established across various human malignancies. In this study, we explored the association between the SPI1 pathway and the MES GSC phenotype. Through comprehensive analysis of the Cancer Genome Atlas and Chinese Glioma Genome Atlas glioma databases, along with patient-derived GSC cultures, we analyzed SPI1 expression. Using genetic knockdown and overexpression techniques, we assessed the functional impact of SPI1 on GSC MES marker expression, invasion, proliferation, self-renewal, and sensitivity to radiation in vitro, as well as its influence on tumor formation in vivo. Additionally, we investigated the downstream signaling cascades activated by SPI1. Our findings revealed a positive correlation between elevated SPI1 expression and the MES phenotype, which in turn, correlated with poor survival. SPI1 enhanced GSC MES differentiation, self-renewal, and radioresistance in vitro, promoting tumorigenicity in vivo. Mechanistically, SPI1 augmented the transcriptional activity of both TGF-β1 and FKBP12 while activating the non-canonical PI3K/Akt pathway. Notably, inhibition of TGF-β1/PI3K/Akt signaling partially attenuated SPI1-induced GSC MES differentiation and its associated malignant phenotype. Collectively, our results underscore SPI1's role in activating TGF-β1/PI3K/Akt signaling through transcriptional upregulation of FKBP12, thereby supporting the aggressive MES phenotype of GSCs. Therefore, SPI1 emerges as a potential therapeutic target in glioma treatment.

Accelerated DNA replication fork speed due to loss of R-loops in myelodysplastic syndromes with SF3B1 mutation

Myelodysplastic syndromes (MDS) with mutated SF3B1 gene present features including a favourable outcome distinct from MDS with mutations in other splicing factor genes SRSF2 or U2AF1. Molecular bases of these divergences are poorly understood. Here we find that SF3B1-mutated MDS show reduced R-loop formation predominating in gene bodies associated with intron retention reduction, not found in U2AF1- or SRSF2-mutated MDS. Compared to erythroblasts from SRSF2- or U2AF1-mutated patients, SF3B1-mutated erythroblasts exhibit augmented DNA synthesis, accelerated replication forks, and single-stranded DNA exposure upon differentiation. Importantly, histone deacetylase inhibition using vorinostat restores R-loop formation, slows down DNA replication forks and improves SF3B1-mutated erythroblast differentiation. In conclusion, loss of R-loops with associated DNA replication stress represents a hallmark of SF3B1-mutated MDS ineffective erythropoiesis, which could be used as a therapeutic target.

A microwell platform for high-throughput longitudinal phenotyping and selective retrieval of organoids

Organoids are powerful experimental models for studying the ontogeny and progression of various diseases including cancer. Organoids are conventionally cultured in bulk using an extracellular matrix mimic. However, bulk-cultured organoids physically overlap, making it impossible to track the growth of individual organoids over time in high throughput. Moreover, local spatial variations in bulk matrix properties make it difficult to assess whether observed phenotypic heterogeneity between organoids results from intrinsic cell differences or differences in the microenvironment. Here, we developed a microwell-based method that enables high-throughput quantification of image-based parameters for organoids grown from single cells, which can further be retrieved from their microwells for molecular profiling. Coupled with a deep learning image-processing pipeline, we characterized phenotypic traits including growth rates, cellular movement, and apical-basal polarity in two CRISPR-engineered human gastric organoid models, identifying genomic changes associated with increased growth rate and changes in accessibility and expression correlated with apical-basal polarity. A record of this paper's transparent peer review process is included in the supplemental information.

FANCA deficiency promotes leukaemic progression by allowing the emergence of cells carrying oncogenic driver mutations

Leukaemia is caused by the clonal evolution of a cell that accumulates mutations/genomic rearrangements, allowing unrestrained cell growth. However, recent identification of leukaemic mutations in the blood cells of healthy individuals revealed that additional events are required to expand the mutated clones for overt leukaemia. Here, we assessed the functional consequences of deleting the Fanconi anaemia A (Fanca) gene, which encodes a DNA damage response protein, in Spi1 transgenic mice that develop preleukaemic syndrome. FANCA loss increases SPI1-associated disease penetrance and leukaemic progression without increasing the global mutation load of leukaemic clones. However, a high frequency of leukaemic FANCA-depleted cells display heterozygous activating mutations in known oncogenes, such as Kit or Nras, also identified but at low frequency in FANCA-WT mice with preleukaemic syndrome, indicating that FANCA counteracts the emergence of oncogene mutated leukaemic cells. A unique transcriptional signature is associated with the leukaemic status of FANCA-depleted cells, leading to activation of MDM4, NOTCH and Wnt/β-catenin pathways. We show that NOTCH signalling improves the proliferation capacity of FANCA-deficient leukaemic cells. Collectively, our observations indicate that loss of the FANC pathway, known to control genetic instability, fosters the expansion of leukaemic cells carrying oncogenic mutations rather than mutation formation. FANCA loss may contribute to this leukaemogenic progression by reprogramming transcriptomic landscape of the cells.

Accurate Identification of DNA Replication Origin by Fusing Epigenomics and Chromatin Interaction Information

DNA replication initiation is a complex process involving various genetic and epigenomic signatures. The correct identification of replication origins (ORIs) could provide important clues for the study of a variety of diseases caused by replication. Here, we design a computational approach named iORI-Epi to recognize ORIs by incorporating epigenome-based features, sequence-based features, and 3D genome-based features. The iORI-Epi displays excellent robustness and generalization ability on both training datasets and independent datasets of K562 cell line. Further experiments confirm that iORI-Epi is highly scalable in other cell lines (MCF7 and HCT116). We also analyze and clarify the regulatory role of epigenomic marks, DNA motifs, and chromatin interaction in DNA replication initiation of eukaryotic genomes. Finally, we discuss gene enrichment pathways from the perspective of ORIs in different replication timing states and heuristically dissect the effect of promoters on replication initiation. Our computational methodology is worth extending to ORI identification in other eukaryotic species.

Topoisomerase 1-dependent R-loop deficiency drives accelerated replication and genomic instability

DNA replication is a complex process tightly regulated to ensure faithful genome duplication, and its perturbation leads to DNA damage and genomic instability. Replication stress is commonly associated with slow and stalled replication forks. Recently, accelerated replication has emerged as a non-canonical form of replication stress. However, the molecular basis underlying fork acceleration is largely unknown. Here, we show that mutated HRAS activation leads to increased topoisomerase 1 (TOP1) expression, causing aberrant replication fork acceleration and DNA damage by decreasing RNA-DNA hybrids or R-loops. In these cells, restoration of TOP1 expression or mild replication inhibition rescues the perturbed replication and reduces DNA damage. Furthermore, TOP1 or RNaseH1 overexpression induces accelerated replication and DNA damage, highlighting the importance of TOP1 equilibrium in regulating R-loop homeostasis to ensure faithful DNA replication and genome integrity. Altogether, our results dissect a mechanism of oncogene-induced DNA damage by aberrant replication fork acceleration.

•

Increased TOP1 expression by mutated RAS reduces R loops

•

Low R-loop levels promote accelerated replication and DNA damage

•

TOP1 restoration or mild replication inhibition rescue DNA acceleration and damage

•

High TOP1 expression is associated with replication mutagenesis in cancer

Replication stress triggered by nucleotide pool imbalance drives DNA damage and cGAS-STING pathway activation in NAFLD

Non-alcoholic steatotic liver disease (NAFLD) is the most common cause of chronic liver disease worldwide. NAFLD has a major effect on the intrinsic proliferative properties of hepatocytes. Here, we investigated the mechanisms underlying the activation of DNA damage response during NAFLD. Proliferating mouse NAFLD hepatocytes harbor replication stress (RS) with an alteration of the replication fork's speed and activation of ATR pathway, which is sufficient to cause DNA breaks. Nucleotide pool imbalance occurring during NAFLD is the key driver of RS. Remarkably, DNA lesions drive cGAS/STING pathway activation, a major component of cells' intrinsic immune response. The translational significance of this study was reiterated by showing that lipid overload in proliferating HepaRG was sufficient to induce RS and nucleotide pool imbalance. Moreover, livers from NAFLD patients displayed nucleotide pathway deregulation and cGAS/STING gene alteration. Altogether, our findings shed light on the mechanisms by which damaged NAFLD hepatocytes might promote disease progression.

Comprehensive Analysis of Immune Implications and Prognostic Value of SPI1 in Gastric Cancer

Background

The transcription factor Spi-1 proto-oncogene (SPI1, also known as PU.1) is a key regulator of signal communication in the immune system and is essential for the development of myeloid cells and lymphocytes. However, the potential role of SPI1 in gastric cancer (GC) and the correlations between SPI1 and immune infiltration remain unclear.

Methods

In the present study, multiple databases including ONCOMINE, TIMER, Kaplan–Meier Plotter, and The Cancer Genome Atlas were used to explore the expression levels and prognostic value of SPI1 in GC. cBioPortal was used to explore the possible reasons for the increased expression of SPI1 in GC. The correlations between SPI1 expression and tumor-infiltrating immune cells (TICs) were analyzed using CIBERSORT and TIMER. Gene set enrichment analysis was used to determine the biological function of SPI1 in the development of GC. In addition, a risk signature based on SPI1-related immunomodulators was constructed to accurately evaluate the prognosis of patients with GC. The upregulation of SPI1 expression in GC was further confirmed through immunohistochemistry, western blotting, and real-time quantitative PCR (RT-qPCR) assay.

Results

The expression of SPI1 was increased significantly in GC according to multiple databases, and high expression of SPI1 was related to poor prognosis and progression of GC. The main factor influencing the high expression of SPI1 mRNA in GC may be diploidy, not DNA methylation. Moreover, immunohistochemistry, western blotting, and RT-qPCR assays also confirmed the upregulated expression of SPI1 in GC. CIBERSORT analysis revealed that SPI1 expression was correlated with seven types of TICs (naive B cells, resting memory CD4 T cells, activated memory CD4 T cells, activated natural killer cells, resting natural killer cells, M2 macrophages, and resting dendritic cells). Gene set enrichment analysis indicated that SPI1 might be related to immune activation in GC and participate in cell cycle regulation. In addition, based on SPI1-related immunomodulators, we developed multiple-gene risk prediction signatures and constructed a nomogram that can independently predict the clinical outcome of GC.

Conclusion

The results of the present study suggest that SPI1 has a critical role in determining the prognosis of GC patients and may be a potential immunotherapeutic target.

A Disease-Causing Single Amino Acid Deletion in the Coiled-Coil Domain of RAD50 Impairs MRE11 Complex Functions in Yeast and Humans

The MRE11-RAD50-NBS1 complex plays a central role in response to DNA double-strand breaks. Here, we identify a patient with bone marrow failure and developmental defects caused by biallelic RAD50 mutations. One of the mutations creates a null allele, whereas the other (RAD50E1035Δ) leads to the loss of a single residue in the heptad repeats within the RAD50 coiled-coil domain. This mutation represents a human RAD50 separation-of-function mutation that impairs DNA repair, DNA replication, and DNA end resection without affecting ATM-dependent DNA damage response. Purified recombinant proteins indicate that RAD50E1035Δ impairs MRE11 nuclease activity. The corresponding mutation in Saccharomyces cerevisiae causes severe thermosensitive defects in both DNA repair and Tel1ATM-dependent signaling. These findings demonstrate that a minor heptad break in the RAD50 coiled coil suffices to impede MRE11 complex functions in human and yeast. Furthermore, these results emphasize the importance of the RAD50 coiled coil to regulate MRE11-dependent DNA end resection in humans.

Co-option of the lineage-specific LAVA retrotransposon in the gibbon genome

Co-option of transposable elements (TEs) to become part of existing or new enhancers is an important mechanism for evolution of gene regulation. However, contributions of lineage-specific TE insertions to recent regulatory adaptations remain poorly understood. Gibbons present a suitable model to study these contributions as they have evolved a lineage-specific TE called LAVA (LINE-AluSz-VNTR-Alu _LIKE), which is still active in the gibbon genome. The LAVA retrotransposon is thought to have played a role in the emergence of the highly rearranged structure of the gibbon genome by disrupting transcription of cell cycle genes. In this study, we investigated whether LAVA may have also contributed to the evolution of gene regulation by adopting enhancer function. We characterized fixed and polymorphic LAVA insertions across multiple gibbons and found 96 LAVA elements overlapping enhancer chromatin states. Moreover, LAVA was enriched in multiple transcription factor binding motifs, was bound by an important transcription factor (PU.1), and was associated with higher levels of gene expression in cis We found gibbon-specific signatures of purifying/positive selection at 27 LAVA insertions. Two of these insertions were fixed in the gibbon lineage and overlapped with enhancer chromatin states, representing putative co-opted LAVA enhancers. These putative enhancers were located within genes encoding SETD2 and RAD9A, two proteins that facilitate accurate repair of DNA double-strand breaks and prevent chromosomal rearrangement mutations. Co-option of LAVA in these genes may have influenced regulation of processes that preserve genome integrity. Our findings highlight the importance of considering lineage-specific TEs in studying evolution of gene regulatory elements.

Mechanisms of Oncogene-Induced Replication Stress: Jigsaw Falling into Place

Oncogene activation disturbs cellular processes and accommodates a complex landscape of changes in the genome that contribute to genomic instability, which accelerates mutation rates and promotes tumorigenesis. Part of this cellular turmoil involves deregulation of physiologic DNA replication, widely described as replication stress. Oncogene-induced replication stress is an early driver of genomic instability and is attributed to a plethora of factors, most notably aberrant origin firing, replication-transcription collisions, reactive oxygen species, and defective nucleotide metabolism.Significance: Replication stress is a fundamental step and an early driver of tumorigenesis and has been associated with many activated oncogenes. Deciphering the mechanisms that contribute to the replication stress response may provide new avenues for targeted cancer treatment. In this review, we discuss the latest findings on the DNA replication stress response and examine the various mechanisms through which activated oncogenes induce replication stress. Cancer Discov; 8(5); 537-55. ©2018 AACR.

The Spi1/PU.1 transcription factor accelerates replication fork progression by increasing PP1 phosphatase in leukemia

Oncogenes trigger replicative stress that can lead to genetic instability, which participates in cancer progression. Thus, determining how cells cope with replicative stress can help our understanding of oncogenesis and lead to the identification of new antitumor treatment targets. We previously showed that constitutive overexpression of the oncogenic transcription factor Spi1/PU.1 leads to pre-leukemic cells that have a shortened S phase duration with an increased replication fork speed and increased mutability in the absence of DNA breaks. Here, we demonstrate that the S phase checkpoint protein CHK1 is maintained in a low phosphorylation state in Spi1/PU.1-overexpressing cells and provide evidence that this is not due to negative control of its primary kinase ATR. Notably, we found that the expression of the CHK1 phosphatase PP1α is increased in Spi1/PU.1-overexpressing cells. By exogenously modulating its activity, we demonstrate that PP1α is required to maintain CHK1 in a dephosphorylated state and, more importantly, that it is responsible for the accelerated replication fork progression in Spi1/PU.1-overexpressing cells. These results identify a novel pathway by which an oncogene influences replication in the absence of DNA damage.

Senescence is a Spi1-induced anti-proliferative mechanism in primary hematopoietic cells

Transcriptional deregulation caused by epigenetic or genetic alterations is a major cause of leukemic transformation. The Spi1/PU.1 transcription factor is a key regulator of many steps of hematopoiesis, and limits self-renewal of hematopoietic stem cells. The deregulation of its expression or activity contributes to leukemia, in which Spi1 can be either an oncogene or a tumor suppressor. Herein we explored whether cellular senescence, an anti-tumoral pathway that restrains cell proliferation, is a mechanism by which Spi1 limits hematopoietic cell expansion, and thus prevents the development of leukemia. We show that Spi1 overexpression triggers cellular senescence both in primary fibroblasts and hematopoietic cells. Erythroid and myeloid lineages are both prone to Spi1-induced senescence. In hematopoietic cells, Spi1-induced senescence requires its DNA-binding activity and a functional p38MAPK14 pathway but is independent of a DNA-damage response. In contrast, in fibroblasts, Spi1-induced senescence is triggered by a DNA-damage response. Importantly, using our well-established Spi1 transgenic leukemia mouse model, we demonstrate that Spi1 overexpression also induces senescence in erythroid progenitors of the bone marrow in vivo before the onset of the pre-leukemic phase of erythroleukemia. Remarkably, the senescence response is lost during the progression of the disease and erythroid blasts do not display a higher expression of Dec1 and CDKN1A, two of the induced senescence markers in young animals. These results bring indirect evidence that leukemia develops from cells which have bypassed Spi1-induced senescence. Overall, our results reveal senescence as a Spi1-induced anti-proliferative mechanism that may be a safeguard against the development of acute myeloid leukemia.

The pan-cancer pathological regulatory landscape

Dysregulation of the normal gene expression program is the cause of a broad range of diseases, including cancer. Detecting the specific perturbed regulators that have an effect on the generation and the development of the disease is crucial for understanding the disease mechanism and for taking decisions on efficient preventive and curative therapies. Moreover, detecting such perturbations at the patient level is even more important from the perspective of personalized medicine. We applied the Transcription Factor Target Enrichment Analysis, a method that detects the activity of transcription factors based on the quantification of the collective transcriptional activation of their targets, to a large collection of 5607 cancer samples covering eleven cancer types. We produced for the first time a comprehensive catalogue of altered transcription factor activities in cancer, a considerable number of them significantly associated to patient’s survival. Moreover, we described several interesting TFs whose activity do not change substantially in the cancer with respect to the normal tissue but ultimately play an important role in patient prognostic determination, which suggest they might be promising therapeutic targets. An additional advantage of this method is that it allows obtaining personalized TF activity estimations for individual patients.

Cdk5 promotes DNA replication stress checkpoint activation through RPA-32 phosphorylation, and impacts on metastasis free survival in breast cancer patients

Cyclin dependent kinase 5 (Cdk5) is a determinant of PARP inhibitor and ionizing radiation (IR) sensitivity. Here we show that Cdk5-depleted (Cdk5-shRNA) HeLa cells show higher sensitivity to S-phase irradiation, chronic hydroxyurea exposure, and 5-fluorouracil and 6-thioguanine treatment, with hydroxyurea and IR sensitivity also seen in Cdk5-depleted U2OS cells. As Cdk5 is not directly implicated in DNA strand break repair we investigated in detail its proposed role in the intra-S checkpoint activation. While Cdk5-shRNA HeLa cells showed altered basal S-phase dynamics with slower replication velocity and fewer active origins per DNA megabase, checkpoint activation was impaired after a hydroxyurea block. Cdk5 depletion was associated with reduced priming phosphorylations of RPA32 serines 29 and 33 and SMC1-Serine 966 phosphorylation, lower levels of RPA serine 4 and 8 phosphorylation and DNA damage measured using the alkaline Comet assay, gamma-H2AX signal intensity, RPA and Rad51 foci, and sister chromatid exchanges resulting in impaired intra-S checkpoint activation and subsequently higher numbers of chromatin bridges. In vitro kinase assays coupled with mass spectrometry demonstrated that Cdk5 can carry out the RPA32 priming phosphorylations on serines 23, 29, and 33 necessary for this checkpoint activation. In addition we found an association between lower Cdk5 levels and longer metastasis free survival in breast cancer patients and survival in Cdk5-depleted breast tumor cells after treatment with IR and a PARP inhibitor. Taken together, these results show that Cdk5 is necessary for basal replication and replication stress checkpoint activation and highlight clinical opportunities to enhance tumor cell killing.

Direct non transcriptional role of NF-Y in DNA replication

NF-Y is a heterotrimeric transcription factor, which plays a pioneer role in the transcriptional control of promoters containing the CCAAT-box, among which genes involved in cell cycle regulation, apoptosis and DNA damage response. The knock-down of the sequence-specific subunit NF-YA triggers defects in S-phase progression, which lead to apoptotic cell death. Here, we report that NF-Y has a critical function in DNA replication progression, independent from its transcriptional activity. NF-YA colocalizes with early DNA replication factories, its depletion affects the loading of replisome proteins to DNA, among which Cdc45, and delays the passage from early to middle-late S phase. Molecular combing experiments are consistent with a role for NF-Y in the control of fork progression. Finally, we unambiguously demonstrate a direct non-transcriptional role of NF-Y in the overall efficiency of DNA replication, specifically in the DNA elongation process, using a Xenopus cell-free system. Our findings broaden the activity of NF-Y on a DNA metabolism other than transcription, supporting the existence of specific TFs required for proper and efficient DNA replication.

Analysis of Genomic Sequence Motifs for Deciphering Transcription Factor Binding and Transcriptional Regulation in Eukaryotic Cells

Eukaryotic genomes contain a variety of structured patterns: repetitive elements, binding sites of DNA and RNA associated proteins, splice sites, and so on. Often, these structured patterns can be formalized as motifs and described using a proper mathematical model such as position weight matrix and IUPAC consensus. Two key tasks are typically carried out for motifs in the context of the analysis of genomic sequences. These are: identification in a set of DNA regions of over-represented motifs from a particular motif database, and de novo discovery of over-represented motifs. Here we describe existing methodology to perform these two tasks for motifs characterizing transcription factor binding. When applied to the output of ChIP-seq and ChIP-exo experiments, or to promoter regions of co-modulated genes, motif analysis techniques allow for the prediction of transcription factor binding events and enable identification of transcriptional regulators and co-regulators. The usefulness of motif analysis is further exemplified in this review by how motif discovery improves peak calling in ChIP-seq and ChIP-exo experiments and, when coupled with information on gene expression, allows insights into physical mechanisms of transcriptional modulation.

Direct Visualization of DNA Replication Dynamics in Zebrafish Cells

Spatiotemporal regulation of DNA replication in the S-phase nucleus has been extensively studied in mammalian cells because it is tightly coupled with the regulation of other nuclear processes such as transcription. However, little is known about the replication dynamics in nonmammalian cells. Here, we analyzed the DNA replication processes of zebrafish (Danio rerio) cells through the direct visualization of replicating DNA in the nucleus and on DNA fiber molecules isolated from the nucleus. We found that zebrafish chromosomal DNA at the nuclear interior was replicated first, followed by replication of DNA at the nuclear periphery, which is reminiscent of the spatiotemporal regulation of mammalian DNA replication. However, the relative duration of interior DNA replication in zebrafish cells was longer compared to mammalian cells, possibly reflecting zebrafish-specific genomic organization. The rate of replication fork progression and ori-to-ori distance measured by the DNA combing technique were ∼ 1.4 kb/min and 100 kb, respectively, which are comparable to those in mammalian cells. To our knowledge, this is a first report that measures replication dynamics in zebrafish cells.

TFmiR: a web server for constructing and analyzing disease-specific transcription factor and miRNA co-regulatory networks

TFmiR is a freely available web server for deep and integrative analysis of combinatorial regulatory interactions between transcription factors, microRNAs and target genes that are involved in disease pathogenesis. Since the inner workings of cells rely on the correct functioning of an enormously complex system of activating and repressing interactions that can be perturbed in many ways, TFmiR helps to better elucidate cellular mechanisms at the molecular level from a network perspective. The provided topological and functional analyses promote TFmiR as a reliable systems biology tool for researchers across the life science communities. TFmiR web server is accessible through the following URL: http://service.bioinformatik.uni-saarland.de/tfmir.

Algorithms for network-based identification of differential regulators from transcriptome data: a systematic evaluation

Identification of differential regulators is critical to understand the dynamics of cellular systems and molecular mechanisms of diseases. Several computational algorithms have recently been developed for this purpose by using transcriptome and network data. However, it remains largely unclear which algorithm performs better under a specific condition. Such knowledge is important for both appropriate application and future enhancement of these algorithms. Here, we systematically evaluated seven main algorithms (TED, TDD, TFactS, RIF1, RIF2, dCSA_t2t, and dCSA_r2t), using both simulated and real datasets. In our simulation evaluation, we artificially inactivated either a single regulator or multiple regulators and examined how well each algorithm detected known gold standard regulators. We found that all these algorithms could effectively discern signals arising from regulatory network differences, indicating the validity of our simulation schema. Among the seven tested algorithms, TED and TFactS were placed first and second when both discrimination accuracy and robustness against data variation were considered. When applied to two independent lung cancer datasets, both TED and TFactS replicated a substantial fraction of their respective differential regulators. Since TED and TFactS rely on two distinct features of transcriptome data, namely differential co-expression and differential expression, both may be applied as mutual references during practical application.

Epigenetic silencing of Bim transcription by Spi-1/PU.1 promotes apoptosis resistance in leukaemia

Deregulation of transcriptional networks contributes to haematopoietic malignancies. The transcription factor Spi-1/PU.1 is a master regulator of haematopoiesis and its alteration leads to leukaemia. Spi-1 overexpression inhibits differentiation and promotes resistance to apoptosis in erythroleukaemia. Here, we show that Spi-1 inhibits mitochondrial apoptosis in vitro and in vivo through the transcriptional repression of Bim, a proapoptotic factor. BIM interacts with MCL-1 that behaves as a major player in the survival of the preleukaemic cells. The repression of BIM expression reduces the amount of BIM-MCL-1 complexes, thus increasing the fraction of potentially active antiapoptotic MCL-1. We then demonstrate that Spi-1 represses Bim transcription by binding to the Bim promoter and by promoting the trimethylation of histone 3 on lysine 27 (H3K27me3, a repressive histone mark) on the Bim promoter. The PRC2 repressive complex of Polycomb is directly responsible for the deposit of H3K27me3 mark at the Bim promoter. SUZ12 and the histone methyltransferase EZH2, two PRC2 subunits bind to the Bim promoter at the same location than H3K27me3, distinct of the Spi-1 DNA binding site. As Spi-1 interacts with SUZ12 and EZH2, these results indicate that Spi-1 modulates the activity of PRC2 without directly recruiting the complex to the site of its activity on the chromatin. Our results identify a new mechanism whereby Spi-1 represses transcription and provide mechanistic insights on the antiapoptotic function of a transcription factor mediated by the epigenetic control of gene expression.

Extracting coordinated patterns of DNA methylation and gene expression in ovarian cancer

OBJECTIVE

DNA methylation, a regulator of gene expression, plays an important role in diverse biological processes including developmental process, carcinogenesis and aging. In particular, aberrant DNA methylation has been largely observed in several types of cancers. Currently, it is important to extract disease-specific gene sets associated with the regulation of DNA methylation.

MATERIALS AND METHODS

Here we propose a novel approach to find the minimum regulatory units of genes, co-methylated and co-expressed gene pairs (MEGP) that are highly correlated gene pairs between DNA methylation and gene expression showing the co-regulatory relationship. To evaluate whether our method is applicable to extract disease-associated genes, we applied our method to a large-scale dataset from the Cancer Genome Atlas extracting significantly associated MEGP and analyzed their functional correlation.

RESULTS

We observed that many MEGP physically interacted with each other and showed high semantic similarity with gene ontology terms. Furthermore, we performed gene set enrichment tests to identify how they are correlated in a complex biological process. Our MEGP were highly enriched in the biological pathway associated with ovarian cancers.

CONCLUSIONS

Our approach is useful for discovering coordinated epigenetic markers associated with specific diseases.

Primary microcephaly, impaired DNA replication, and genomic instability caused by compound heterozygous ATR mutations

Ataxia telangiectasia-mutated (ATM) and ataxia telangiectasia and Rad3-related (ATR) kinases are two key regulators of DNA-damage responses (DDR) that are mainly activated in response to DNA double-strand breaks and single-stranded DNA damages, respectively. Seckel syndrome, a rare genetic disorder characterized by a microcephaly and a markedly reduced body size, has been associated with defective ATR-dependent DNA damage signaling. However, the only human genetic ATR defect reported so far is a hypomorphic splicing mutation identified in five related individuals with Seckel syndrome. Here, we report the first case of primary microcephaly with compound heterozygous mutations in ATR: a 540 kb genomic deletion on one allele and a missense mutation leading to splice dysregulation on the other, which ultimately lead to a sharp decrease in ATR expression. DNA combing technology revealed a profound spontaneous alteration of several DNA replication parameters in patient's cells and FISH analyses highlighted the genomic instability caused by ATR deficiency. Collectively, our results emphasize the crucial role for ATR in the control of DNA replication, and reinforce the complementary and nonredundant contributions of ATM and ATR in human cells to face DNA damages and warrant genome integrity.

Spi-1/PU.1 activates transcription through clustered DNA occupancy in erythroleukemia

Acute leukemias are characterized by deregulation of transcriptional networks that control the lineage specificity of gene expression. The aberrant overexpression of the Spi-1/PU.1 transcription factor leads to erythroleukemia. To determine how Spi-1 mechanistically influences the transcriptional program, we combined a ChIP-seq analysis with transcriptional profiling in cells from an erythroleukemic mouse model. We show that Spi-1 displays a selective DNA-binding that does not often cause transcriptional modulation. We report that Spi-1 controls transcriptional activation and repression partially through distinct Spi-1 recruitment to chromatin. We revealed several parameters impacting on Spi-1-mediated transcriptional activation. Gene activation is facilitated by Spi-1 occupancy close to transcriptional starting site of genes devoid of CGIs. Moreover, in those regions Spi-1 acts by binding to multiple motifs tightly clustered and with similar orientation. Finally, in contrast to the myeloid and lymphoid B cells in which Spi-1 exerts a physiological activity, in the erythroleukemic cells, lineage-specific cooperating factors do not play a prevalent role in Spi-1-mediated transcriptional activation. Thus, our work describes a new mechanism of gene activation through clustered site occupancy of Spi-1 particularly relevant in regard to the strong expression of Spi-1 in the erythroleukemic cells.

DNA replication fading as proliferating cells advance in their commitment to terminal differentiation

Terminal differentiation is the process by which cycling cells stop proliferating to start new specific functions. It involves dramatic changes in chromatin organization as well as gene expression. In the present report we used cell flow cytometry and genome wide DNA combing to investigate DNA replication during murine erythroleukemia-induced terminal cell differentiation. The results obtained indicated that the rate of replication fork movement slows down and the inter-origin distance becomes shorter during the precommitment and commitment periods before cells stop proliferating and accumulate in G1. We propose this is a general feature caused by the progressive heterochromatinization that characterizes terminal cell differentiation.

Evidence for sequential and increasing activation of replication origins along replication timing gradients in the human genome

Genome-wide replication timing studies have suggested that mammalian chromosomes consist of megabase-scale domains of coordinated origin firing separated by large originless transition regions. Here, we report a quantitative genome-wide analysis of DNA replication kinetics in several human cell types that contradicts this view. DNA combing in HeLa cells sorted into four temporal compartments of S phase shows that replication origins are spaced at 40 kb intervals and fire as small clusters whose synchrony increases during S phase and that replication fork velocity (mean 0.7 kb/min, maximum 2.0 kb/min) remains constant and narrowly distributed through S phase. However, multi-scale analysis of a genome-wide replication timing profile shows a broad distribution of replication timing gradients with practically no regions larger than 100 kb replicating at less than 2 kb/min. Therefore, HeLa cells lack large regions of unidirectional fork progression. Temporal transition regions are replicated by sequential activation of origins at a rate that increases during S phase and replication timing gradients are set by the delay and the spacing between successive origin firings rather than by the velocity of single forks. Activation of internal origins in a specific temporal transition region is directly demonstrated by DNA combing of the IGH locus in HeLa cells. Analysis of published origin maps in HeLa cells and published replication timing and DNA combing data in several other cell types corroborate these findings, with the interesting exception of embryonic stem cells where regions of unidirectional fork progression seem more abundant. These results can be explained if origins fire independently of each other but under the control of long-range chromatin structure, or if replication forks progressing from early origins stimulate initiation in nearby unreplicated DNA. These findings shed a new light on the replication timing program of mammalian genomes and provide a general model for their replication kinetics.

Eukaryotic chromosomes replicate from multiple replication origins that fire at different times in S phase. The mechanisms that specify origin position and firing time and coordinate origins to ensure complete genome duplication are unclear. Previous studies proposed either that origins are arranged in temporally coordinated groups or fire independently of each other in a stochastic manner. Here, we have performed a quantitative analysis of human genome replication kinetics using a combination of DNA combing, which reveals local patterns of origin firing and replication fork progression on single DNA molecules, and massive sequencing of newly replicated DNA, which reveals the population-averaged replication timing profile of the entire genome. We show that origins are activated synchronously in large regions of uniform replication timing but more gradually in temporal transition regions and that the rate of origin firing increases as replication progresses. Large regions of unidirectional fork progression are abundant in embryonic stem cells but rare in differentiated cells. We propose a model in which replication forks progressing from early origins stimulate initiation in nearby unreplicated DNA in a manner that explains the shape of the replication timing profile. These results provide a fundamental insight into the temporal regulation of mammalian genome replication.