Identification of chromatin remodeling genes Arid4a and Arid4b as leukemia suppressor genes

Epigenetic regulation of human FOXP3+ Tregs: from homeostasis maintenance to pathogen defense

Regulatory T cells (Tregs), characterized by the expression of Forkhead Box P3 (FOXP3), constitute a distinct subset of T cells crucial for immune regulation. Tregs can exert direct and indirect control over immune homeostasis by releasing inhibitory factors or differentiating into Th-like Treg (Th-Treg), thereby actively contributing to the prevention and treatment of autoimmune diseases. The epigenetic regulation of FOXP3, encompassing DNA methylation, histone modifications, and post-translational modifications, governs the development and optimal suppressive function of Tregs. In addition, Tregs can also possess the ability to maintain homeostasis in diverse microenvironments through non-suppressive mechanisms. In this review, we primarily focus on elucidating the epigenetic regulation of Tregs as well as their multifaceted roles within diverse physiological contexts while looking forward to potential strategies involving augmentation or suppression of Tregs activity for disease management, particularly in light of the ongoing global COVID-19 pandemic.

Urinary Exosomal miR-17-5p Accelerates Bladder Cancer Invasion by Repressing its Target Gene ARID4B and Regulating the Immune Microenvironment

BACKGROUND

Urothelial bladder cancer (BCa) is a common malignant tumor of the urinary system. It has been identified that exosomal miRNAs contribute to the development of BCa. However, its significance and mechanism in the malignant biological behavior of BCa remain unclear. In this study, the influence of exosomal miRNAs on BCa progression was investigated.

METHODS

High-throughput sequencing was conducted to analyze the microRNA-expression profile in urinary exosomes to screen out the key miRNA of muscle-invasive bladder cancer (MIBC). Then, candidate miRNA expression was verified and validated in urinary exosomes and tissue samples. To address the potential role of the candidate miRNA, we overexpressed and knocked down the candidate miRNA and explored its activity in BCa cell lines. Furthermore, the target gene of the selected miRNA was predicted and validated.

RESULTS

The expression profile of miRNAs revealed increased expression of miR-17-5p in MIBC urinary exosomes, and this was later confirmed in urinary exosomes and tissue samples. Cell function studies revealed that exosomal miR-17-5p significantly promoted the growth and invasion of BCa cells. Bioinformatics and luciferase experiments demonstrated that the ARID4B mRNA 3' UTR might be the binding site for miR-17-5p. Low ARID4B levels were linked to high-grade BCa patients and were associated with a better prognosis.

CONCLUSION

Elevated miR-17-5p contributes to BCa progression by targeting ARID4B and influencing the immune system. Based on these findings, miR-17-5p has the potential to be a new therapeutic target for the treatment of BCa.

Comprehensive Atlas of Alternative Splicing Reveals NSRP1 Promoting Adipogenesis through CCDC18

Alternative splicing (AS) plays a crucial role in regulating gene expression, function, and diversity. However, limited reports exist on the identification and comparison of AS in Eastern and Western pigs. Here, we analyzed 243 transcriptome data from eight tissues, integrating information on transcription factors (TFs), selection signals, splicing factors (SFs), and quantitative trait loci (QTL) to comprehensively study alternative splicing events (ASEs) in pigs. Five ASE types were identified, with Mutually Exclusive Exon (MXE) and Skipped Exon (SE) ASEs being the most prevalent. A significant portion of genes with ASEs (ASGs) showed conservation across all eight tissues (63.21-76.13% per tissue). Differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) exhibited tissue specificity, with blood and adipose tissues having more DASGs. Functional enrichment analysis revealed coDASG_DEGs in adipose were enriched in pathways associated with adipose deposition and immune inflammation, while coDASG_DEGs in blood were enriched in pathways related to immune inflammation and metabolism. Adipose deposition in Eastern pigs might be linked to the down-regulation of immune-inflammation-related pathways and reduced insulin resistance. The TFs, selection signals, and SFs appeared to regulate ASEs. Notably, ARID4A (TF), NSRP1 (SF), ANKRD12, IFT74, KIAA2026, CCDC18, NEXN, PPIG, and ROCK1 genes in adipose tissue showed potential regulatory effects on adipose-deposition traits. NSRP1 could promote adipogenesis by regulating alternative splicing and expression of CCDC18. Conducting an in-depth investigation into AS, this study has successfully identified key marker genes essential for pig genetic breeding and the enhancement of meat quality, which will play important roles in promoting the diversity of pork quality and meeting market demand.

Quantitative association between gene expression and blood cell production of individual hematopoietic stem cells in mice

Individual hematopoietic stem cells (HSCs) produce different amounts of blood cells upon transplantation. Taking advantage of the intercellular variation, we developed an experimental and bioinformatic approach to evaluating the quantitative association between gene expression and blood cell production across individual HSCs. We found that most genes associated with blood production exhibit the association only at some levels of blood production. By mapping gene expression with blood production, we identified four distinct patterns of their quantitative association. Some genes consistently correlate with blood production over a range of levels or across all levels, and these genes are found to regulate lymphoid but not myeloid production. Other genes exhibit one or more clear peaks of association. Genes with overlapping peaks are found to be coexpressed in other tissues and share similar molecular functions and regulatory motifs. By dissecting intercellular variations, our findings revealed four quantitative association patterns that reflect distinct dose-response molecular mechanisms modulating the blood cell production of HSCs.

Heterozygous Knockout of ARID4B Using CRISPR/Cas9 Attenuates Some Aggressive Phenotypes in a Breast Cancer Cell Line

Breast cancer is one of the leading causes of death in women around the world. Over time, many genes and mutations that are associated with the development of this disease have been identified. However, the specific role of many genes has not yet been fully elucidated. Higher ARID4B expression has been identified as a risk factor for diverse cancer types. Silencing experiments also showed that ARID4B is associated with developing cancer-associated characteristics. However, no transcriptomic studies have shown the overall cellular effect of loss of function in breast cancer in humans. This study addresses the impact of loss-of-function mutations in breast cancer MCF-7 cells. Using the CRISPR/Cas9 system, we generated mutations that caused heterozygous truncated proteins, isolating three monoclonal lines carrying insertions and deletions in ARID4B. We observed reduced proliferation and migration in in vitro experiments. In addition, from RNA-seq assays, a differential expression analysis shows known and novel deregulated cancer-associate pathways in mutated cells supporting the impact of ARID4B. For example, we found the AKT-PI3K pathway to be altered at the transcript level but through different genes than those reported for ARID4B. Our transcriptomic results also suggest new insights into the role of ARID4B in aggressiveness by the epithelial-to-mesenchymal transition and TGF-β pathways and in metabolism through cholesterol and mevalonate pathways. We also performed exome sequencing to show that no off-target effects were apparent. In conclusion, the ARID4B gene is associated with some aggressive phenotypes in breast cancer cells.

A research review of experimental animal models with myelodysplastic syndrome

Myelodysplastic syndrome (MDS) consists of a group of hematologic tumors that are derived from the clonal proliferation of hematopoietic stem cells, featuring abnormal hematopoietic cell development and ineffective hematopoiesis. Animal models are an important scientific research platform that has been widely applied in the research of human diseases, especially tumors. Animal models with MDS can simulate characteristic human genetic variations and tumor phenotypes. They also provide a reliable platform for the exploration of the pathogenesis and diagnostic markers of MDS as well as for a drug efficacy evaluation. This paper reviews the research status of three animal models and a new spontaneous mouse model with MDS.

Epigenetic insight into effects of prenatal alcohol exposure on stress axis development: Systematic review with meta-analytic approaches

We conducted a systematic review with meta-analytic elements using publicly available Gene Expression Omnibus (GEO) datasets to determine the role of epigenetic mechanisms in prenatal alcohol exposure (PAE)-induced hypothalamic-pituitary-adrenal (HPA) axis dysfunctions in offspring. Several studies have demonstrated that PAE has long-term consequences on HPA axis functions in offspring. Some studies determined that alcohol-induced epigenetic alterations during fetal development persist in adulthood. However, additional research is needed to understand the major epigenetic events leading to alcohol-induced teratogenesis of the HPA axis. Our network analysis of GEO datasets identified key pathways relevant to alcohol-mediated histone modifications, DNA methylation, and miRNA involvement associated with PAE-induced alterations of the HPA axis. Our analysis indicated that PAE perturbated the epigenetic machinery to activate corticotrophin-releasing hormone, while it suppressed opioid, glucocorticoid receptor, and circadian clock genes. These results help to further our understanding of the epigenetic basis of alcohol's effects on HPA axis development.

Bone marrow confined IL-6 signaling mediates the progression of myelodysplastic syndromes to acute myeloid leukemia

Myelodysplastic syndromes (MDS) are age-related myeloid neoplasms with increased risk of progression to acute myeloid leukemia (AML). The mechanisms of transformation of MDS to AML are poorly understood, especially in relation to the aging microenvironment. We previously established an mDia1/miR-146a double knockout (DKO) mouse model phenocopying MDS. These mice develop age-related pancytopenia with oversecretion of proinflammatory cytokines. Here, we found that most of the DKO mice underwent leukemic transformation at 12–14 months of age. These mice showed myeloblast replacement of fibrotic bone marrow and widespread leukemic infiltration. Strikingly, depletion of IL-6 in these mice largely rescued the leukemic transformation and markedly extended survival. Single-cell RNA sequencing analyses revealed that DKO leukemic mice had increased monocytic blasts that were reduced with IL-6 knockout. We further revealed that the levels of surface and soluble IL-6 receptor (IL-6R) in the bone marrow were significantly increased in high-risk MDS patients. Similarly, IL-6R was also highly expressed in older DKO mice. Blocking of IL-6 signaling significantly ameliorated AML progression in the DKO model and clonogenicity of CD34-positive cells from MDS patients. Our study establishes a mouse model of progression of age-related MDS to AML and indicates the clinical significance of targeting IL-6 signaling in treating high-risk MDS.

A novel heterozygous missense variant of the ARID4A gene identified in Han Chinese families with schizophrenia-diagnosed siblings that interferes with DNA-binding activity

ARID4A plays an important role in regulating gene expression and cell proliferation. ARID4A belongs to the AT-rich interaction domain (ARID)-containing family, and a PWWP domain immediately precedes its ARID region. The molecular mechanism and structural basis of ARID4A are largely unknown. Whole-exome sequencing (WES) revealed that a novel heterozygous missense variant, ARID4A c.1231 C > G (p.His411Asp), was associated with schizophrenia (SCZ) in this study. We determined the crystal structure of the PWWP-ARID tandem at 2.05 Å, revealing an unexpected mode in which ARID4A assembles with its PWWP and ARID from a structural and functional supramodule. Our results further showed that compared with the wild type, the p.His411Asp ARID mutant protein adopts a less compact conformation and exhibits a weaker dsDNA-binding ability. The p.His411Asp mutation decreased the number of cells that were arrested in the G0-G1 phase and caused more cells to progress to the G2-M phase. In addition, the missense mutation promoted the proliferation of HEK293T cells. In conclusion, our data provide evidence that ARID4A p.His411Asp could cause a conformational change in the ARID4A ARID domain, influence the DNA binding function, and subsequently disturb the cell cycle arrest in the G1 phase. ARID4A is likely a susceptibility gene for SCZ; thus, these findings provide new insight into the role of ARID4A in psychiatric disorders.

The Coming of Age of Preclinical Models of MDS

Myelodysplastic syndromes (MDS) are a heterogeneous group of clonal bone-marrow diseases with ineffective hematopoiesis resulting in cytopenias and morphologic dysplasia of hematopoietic cells. MDS carry a wide spectrum of genetic abnormalities, ranging from chromosomal abnormalities such as deletions/additions, to recurrent mutations affecting the spliceosome, epigenetic modifiers, or transcription factors. As opposed to AML, research in MDS has been hindered by the lack of preclinical models that faithfully replicate the complexity of the disease and capture the heterogeneity. The complex molecular landscape of the disease poses a unique challenge when creating transgenic mouse-models. In addition, primary MDS cells are difficult to manipulate ex vivo limiting in vitro studies and resulting in a paucity of cell lines and patient derived xenograft models. In recent years, progress has been made in the development of both transgenic and xenograft murine models advancing our understanding of individual contributors to MDS pathology as well as the complex primary interplay of genetic and microenvironment aberrations. We here present a comprehensive review of these transgenic and xenograft models for MDS and future directions.

Reynolds L, Enquobahrie DA, L. Fitzpatrick A, Kerr KF, J. Budoff M, Lee SI, Siscovick D, D. Kaufman J. Epigenome-wide analysis of long-term air pollution exposure and DNA methylation in monocytes: results from the Multi-Ethnic Study of Atherosclerosis

Air pollution might affect atherosclerosis through DNA methylation changes in cells crucial to atherosclerosis, such as monocytes. We conducted an epigenome-wide study of DNA methylation in CD14+ monocytes and long-term ambient air pollution exposure in adults participating in the Multi-Ethnic Study of Atherosclerosis (MESA). We also assessed the association between differentially methylated signals and cis-gene expression. Using spatiotemporal models, one-year average concentrations of outdoor fine particulate matter (PM_2.5) and oxides of nitrogen (NO_X) were estimated at participants' homes. We assessed DNA methylation and gene expression using Illumina 450k and HumanHT-12 v4 Expression BeadChips, respectively (n = 1,207). We used bump hunting and site-specific approaches to identify differentially methylated signals (false discovery rate of 0.05) and used linear models to assess associations between differentially methylated signals and cis-gene expression. Four differentially methylated regions (DMRs) located on chromosomes 5, 6, 7, and 16 (within or near SDHAP3, ZFP57, HOXA5, and PRM1, respectively) were associated with PM_2.5. The DMRs on chromosomes 5 and 6 also associated with NO_X. The DMR on chromosome 5 had the smallest p-value for both PM_2.5 (p = 1.4×10^-6) and NO_X (p = 7.7×10^-6). Three differentially methylated CpGs were identified for PM_2.5, and cg05926640 (near TOMM20) had the smallest p-value (p = 5.6×10^-8). NO_X significantly associated with cg11756214 within ZNF347 (p = 5.6×10^-8). Several differentially methylated signals were also associated with cis-gene expression. The DMR located on chromosome 7 was associated with the expression of HOXA5, HOXA9, and HOXA10. The DMRs located on chromosomes 5 and 16 were associated with expression of MRPL36 and DEXI, respectively. The CpG cg05926640 was associated with expression of ARID4B, IRF2BP2, and TOMM20. We identified differential DNA methylation in monocytes associated with long-term air pollution exposure. Methylation signals associated with gene expression might help explain how air pollution contributes to cardiovascular disease.

Positive Selection and Enhancer Evolution Shaped Lifespan and Body Mass in Great Apes

Within primates, the great apes are outliers both in terms of body size and lifespan, since they include the largest and longest-lived species in the order. Yet, the molecular bases underlying such features are poorly understood. Here, we leveraged an integrated approach to investigate multiple sources of molecular variation across primates, focusing on over 10,000 genes, including approximately 1,500 previously associated with lifespan, and additional approximately 9,000 for which an association with longevity has never been suggested. We analyzed dN/dS rates, positive selection, gene expression (RNA-seq), and gene regulation (ChIP-seq). By analyzing the correlation between dN/dS, maximum lifespan, and body mass, we identified 276 genes whose rate of evolution positively correlates with maximum lifespan in primates. Further, we identified five genes, important for tumor suppression, adaptive immunity, metastasis, and inflammation, under positive selection exclusively in the great ape lineage. RNA-seq data, generated from the liver of six species representing all the primate lineages, revealed that 8% of approximately 1,500 genes previously associated with longevity are differentially expressed in apes relative to other primates. Importantly, by integrating RNA-seq with ChIP-seq for H3K27ac (which marks active enhancers), we show that the differentially expressed longevity genes are significantly more likely than expected to be located near a novel "ape-specific" enhancer. Moreover, these particular ape-specific enhancers are enriched for young transposable elements, and specifically SINE-Vntr-Alus. In summary, we demonstrate that multiple evolutionary forces have contributed to the evolution of lifespan and body size in primates.

Differentiation of fetal hematopoietic stem cells requires ARID4B to restrict autocrine KITLG/KIT-Src signaling

Balance between the hematopoietic stem cell (HSC) duality to either possess self-renewal capacity or differentiate into multipotency progenitors (MPPs) is crucial for maintaining homeostasis of the hematopoietic stem/progenitor cell (HSPC) compartment. To retain the HSC self-renewal activity, KIT, a receptor tyrosine kinase, in HSCs is activated by its cognate ligand KITLG originating from niche cells. Here, we show that AT-rich interaction domain 4B (ARID4B) interferes with KITLG/KIT signaling, consequently allowing HSC differentiation. Conditional Arid4b knockout in mouse hematopoietic cells blocks fetal HSC differentiation, preventing hematopoiesis. Mechanistically, ARID4B-deficient HSCs self-express KITLG and overexpress KIT. As to downstream pathways of KITLG/KIT signaling, inhibition of Src family kinases rescues the HSC differentiation defect elicited by ARID4B loss. In summary, the intrinsic ARID4B-KITLG/KIT-Src axis is an HSPC regulatory program that enables the differentiation state, while KIT stimulation by KITLG from niche cells preserves the HSPC undifferentiated pool.

Hematopoietic stem cells (HSCs) at the top of the hematopoietic hierarchy are able to self-renew and differentiate to mature blood cells. Young et al. report that an HSC self-control mechanism established by ARID4B ensures HSC differentiation. ARID4B-deficient HSCs produce KITLG to stimulate KIT, leading to blockage of HSC differentiation and eventual hematopoietic failure.

Mouse Models of CMML

Chronic myelomonocytic leukemia (CMML) is a rare and challenging type of myeloproliferative neoplasm. Poor prognosis and high mortality, associated predominantly with progression to secondary acute myeloid leukemia (sAML), is still an unsolved problem. Despite a growing body of knowledge about the molecular repertoire of this disease, at present, the prognostic significance of CMML-associated mutations is controversial. The absence of available CMML cell lines and the small number of patients with CMML make pre-clinical testing and clinical trials complicated. Currently, specific therapy for CMML has not been approved; most of the currently available therapeutic approaches are based on myelodysplastic syndrome (MDS) and other myeloproliferative neoplasm (MNP) studies. In this regard, the development of the robust CMML animal models is currently the focus of interest. This review describes important studies concerning animal models of CMML, examples of methodological approaches, and the obtained hematologic phenotypes.

Structural Insight into Chromatin Recognition by Multiple Domains of the Tumor Suppressor RBBP1

Retinoblastoma-binding protein 1 (RBBP1) is involved in gene regulation, epigenetic regulation, and disease processes. RBBP1 contains five domains with DNA-binding or histone-binding activities, but how RBBP1 specifically recognizes chromatin is still unknown. An AT-rich interaction domain (ARID) in RBBP1 was proposed to be the key region for DNA-binding and gene suppression. Here, we first determined the solution structure of a tandem PWWP-ARID domain mutant of RBBP1 after deletion of a long flexible acidic loop L12 in the ARID domain. NMR titration results indicated that the ARID domain interacts with DNA with no GC- or AT-rich preference. Surprisingly, we found that the loop L12 binds to the DNA-binding region of the ARID domain as a DNA mimic and inhibits DNA binding. The loop L12 can also bind weakly to the Tudor and chromobarrel domains of RBBP1, but binds more strongly to the DNA-binding region of the histone H2A-H2B heterodimer. Furthermore, both the loop L12 and DNA can enhance the binding of the chromobarrel domain to H3K4me3 and H4K20me3. Based on these results, we propose a model of chromatin recognition by RBBP1, which highlights the unexpected multiple key roles of the disordered acidic loop L12 in the specific binding of RBBP1 to chromatin.

ARID4B is critical for mouse embryonic stem cell differentiation towards mesoderm and endoderm, linking epigenetics to pluripotency exit

Distinct cell types emerge from embryonic stem cells through a precise and coordinated execution of gene expression programs during lineage commitment. This is established by the action of lineage specific transcription factors along with chromatin complexes. Numerous studies have focused on epigenetic factors that affect embryonic stem cells (ESC) self-renewal and pluripotency. However, the contribution of chromatin to lineage decisions at the exit from pluripotency has not been as extensively studied. Using a pooled epigenetic shRNA screen strategy, we identified chromatin-related factors critical for differentiation toward mesodermal and endodermal lineages. Here we reveal a critical role for the chromatin protein, ARID4B. Arid4b-deficient mESCs are similar to WT mESCs in the expression of pluripotency factors and their self-renewal. However, ARID4B loss results in defects in up-regulation of the meso/endodermal gene expression program. It was previously shown that Arid4b resides in a complex with SIN3A and HDACS 1 and 2. We identified a physical and functional interaction of ARID4B with HDAC1 rather than HDAC2, suggesting functionally distinct Sin3a subcomplexes might regulate cell fate decisions Finally, we observed that ARID4B deficiency leads to increased H3K27me3 and a reduced H3K27Ac level in key developmental gene loci, whereas a subset of genomic regions gain H3K27Ac marks. Our results demonstrate that epigenetic control through ARID4B plays a key role in the execution of lineage-specific gene expression programs at pluripotency exit.

ARID4B Knockdown Suppresses PI3K/AKT Signaling and Induces Apoptosis in Human Glioma Cells

PURPOSE

Glioblastoma multiforme is a highly malignant primary brain cancer with a poor prognosis. We recently reported that ARID4B could potentially serve as a biomarker associated with poor survival in glioma patients. However, the function of ARID4B in human gliomas remains unclear. The aim of this study is to investigate the molecular cell biology role of ARID4B in human glioma cells.

MATERIALS AND METHODS

Gene Expression Omnibus (GEO) and Human Protein Atlas (HPA) datasets were analyzed for the expression of ARID4B in WHO pathological grading, overall survival and immunohistochemical staining. Using quantitative RT-PCR and Western blotting, those findings were confirmed in normal brain tissue and glioma cell lines. ARID4B knockdown was conducted via lentivirus-based transfection of small hairpin RNA in human glioma cells to investigate cell proliferation, cell cycle, and apoptosis.

RESULTS

In the present study, our analysis of GEO datasets showed that ARID4B mRNA expression is higher in WHO grade IV tumors (n = 81) than in non-tumor control tissue (n = 23, P <0.0001). ARID4B knockdown suppressed glioma cell proliferation and induced G1 phase arrest via the PI3K/AKT pathway. It also increased expression of HDAC1, leading to higher acetyl-p53 and acetyl-H3 levels and reduced glioma cell migration and invasion. These effects were mediated via downregulation of AKT pathway components, including p-mTOR, p-PI3K and p-AKT. ARID4B knockdown also led to downregulation of Cyclin D1, which increased apoptosis in human glioma cells.

CONCLUSION

These findings that ARID4B expression correlates positively with WHO pathologic grading in glioma. ARID4B knockdown suppresses PI3K/AKT signaling and induces apoptosis in human glioma cells. These results suggests that ARID4B acts as an oncogene in human gliomas.

Structural basis for the DNA-binding activity of human ARID4B Tudor domain

Human ARID4A and ARID4B are homologous proteins that are important in controlling gene expression and epigenetic regulation but have distinct functions. Previous studies have shown that the N-terminal domain of ARID4A is an unusual interdigitated double Tudor domain with DNA-binding activity. However, how the Tudor domain of ARID4B differs from that of ARID4A remains unknown. Here, we found that the ARID4B Tudor domain has significantly weaker DNA affinity than the ARID4A Tudor domain despite sharing more than 80% sequence identity. Structure determination and DNA titration analysis indicated that the ARID4B Tudor domain is also an interdigitated double Tudor domain with a DNA-binding surface similar to ARID4A. We identified a residue close to the DNA-binding site of the Tudor domain that differs between ARID4A and ARID4B. The Leu50 in ARID4A is Glu50 in ARID4B, and the latter forms salt bridges with two lysine residues at the DNA-binding surface. This causes a decrease in the strength of positive charge, thus reducing DNA-binding affinity while significantly increasing protein stability. We also found that a C-terminal extension region enhances the DNA-binding affinity of the ARID4B Tudor domain. This C-terminal extension is disordered and contains a positively charged RGR motif, providing an additional DNA-binding site. Finally, sequence and phylogenetic analyses indicated that the residue differences and the presence of the RGR extension region are conserved. These results provide new insight into the functional differences between ARID4A and ARID4B proteins, as well as elucidating the function of the disordered regions in these proteins.

System Genetics Including Causal Inference Identify Immune Targets for Coronary Artery Disease and the Lifespan

BACKGROUND

Randomized clinical trials indicate that the immune response plays a significant role in coronary artery disease (CAD), a disorder impacting the lifespan potential. However, the identification of targets critical to the immune response in atheroma is still hampered by a lack of solid inference.

METHODS

Herein, we implemented a system genetics approach to identify causally associated immune targets implicated in atheroma. We leveraged genome-wide association studies to perform mapping and Mendelian randomization to assess causal associations between gene expression in blood cells with CAD and the lifespan. Expressed genes (eGenes) were prioritized in network and in single-cell expression derived from plaque immune cells.

RESULTS

Among 840 CAD-associated blood eGenes, 37 were predicted causally associated with CAD and 6 were also associated with the parental lifespan in Mendelian randomization. In multivariable Mendelian randomization, the impact of eGenes on the lifespan potential was mediated by the CAD risk. Predicted causal eGenes were central in network. FLT1 and CCR5 were identified as targets of approved drugs, whereas 22 eGenes were deemed tractable for the development of small molecules and antibodies. Analyses of plaque immune single-cell expression identified predicted causal eGenes enriched in macrophages (GPX1, C4orf3) and involved in ligand-receptor interactions (CCR5).

CONCLUSIONS

We identified 37 blood eGenes predicted causally associated with CAD. The predicted expression for 6 eGenes impacted the lifespan potential through the risk of CAD. Prioritization based on network, annotations, and single-cell expression identified targets deemed tractable for the development of drugs and for drug repurposing.

Germline IKAROS dimerization haploinsufficiency causes hematologic cytopenias and malignancies

IKAROS is a transcription factor forming homo- and heterodimers and regulating lymphocyte development and function. Germline mutations affecting the IKAROS N-terminal DNA binding domain, acting in a haploinsufficient or dominant-negative manner, cause immunodeficiency. Herein, we describe 4 germline heterozygous IKAROS variants affecting its C-terminal dimerization domain, via haploinsufficiency, in 4 unrelated families. Index patients presented with hematologic disease consisting of cytopenias (thrombocytopenia, anemia, neutropenia)/Evans syndrome and malignancies (T-cell acute lymphoblastic leukemia, Burkitt lymphoma). These dimerization defective mutants disrupt homo- and heterodimerization in a complete or partial manner, but they do not affect the wild-type allele function. Moreover, they alter key mechanisms of IKAROS gene regulation, including sumoylation, protein stability, and the recruitment of the nucleosome remodeling and deacetylase complex; none affected in N-terminal DNA binding defects. These C-terminal dimerization mutations are largely associated with hematologic disorders, display dimerization haploinsufficiency and incomplete clinical penetrance, and differ from previously reported allelic variants in their mechanism of action. Dimerization mutants contribute to the growing spectrum of IKAROS-associated diseases displaying a genotype-phenotype correlation.

Circ-SMARCA5 suppresses colorectal cancer progression via downregulating miR-39-3p and upregulating ARID4B

BACKGROUND

Circular RNAs are crucial in tumorigenesis. However, little is known about their functions in colorectal cancer (CRC). Circ-SMARCA5 was found to be an oncogene or tumor suppresser in different types of cancers, but its exact role in CRC remains unknown. Here, we aim to identify the role of circ-SMARCA5 in CRC development.

METHODS

Circ-SMARCA5 expression was determined by qRT-PCR. CRC cell proliferation, migration, and invasion were detected by CCK-8, wound healing, and Transwell assays, respectively. Bioinformatics analysis was performed to predict target genes. The interaction of microRNA (miR) with circ-SMARCA5 or target genes was detected using luciferase reporter assay. Xenograft model was established to determine the effect of circ-SMARCA5 on CRC tumor growth in vivo.

RESULTS

Circ-SMARCA5 expression was dramatically decreased in CRC cell lines and tissues. Circ-SMARCA5 overexpression inhibited CRC cell proliferation, migration and invasion. MiR-93-3p was predicted as a target of circ-SMARCA5 and its overexpression attenuated the anti-tumor effect of circ-SMARCA5 on CRC cells. Furthermore, we predicted AT-rich interaction domain 4B (ARID4B) as the target of miR-39-3p. Functional analysis showed that circ-SMARCA5 upregulated ARID4B expression via miR-39-3p. Additionally, in vivo studies demonstrated that circ-SMARCA5 suppressed CRC tumor progression.

CONCLUSION

Circ-SMARCA5 functions as a tumor suppressor by upregulating ARID4B expression via sponging miR-39-3p, and thereby inhibited CRC progression.

A regulatory role for CHD2 in myelopoiesis

The transcriptional program that dictates haematopoietic cell fate and differentiation requires an epigenetic regulatory and memory function, provided by a network of epigenetic factors that regulate DNA methylation, post-translational histone modifications and chromatin structure. Disturbed epigenetic regulation causes perturbations in the blood cell differentiation program that results in various types of haematopoietic disorders. Thus, accurate epigenetic regulation is essential for functional haematopoiesis. In this study, we used a CRISPR-Cas9 screening approach to identify new epigenetic regulators in myeloid differentiation. We designed a Chromatin-UMI CRISPR guide library targeting 1092 epigenetic regulators. Phorbol 12-myristate 13-acetate (PMA) treatment of the chronic myeloid leukaemia cell line K-562 was used as a megakaryocytic myeloid differentiation model. Both previously described developmental epigenetic regulators and novel factors were identified in our screen. In this study, we validated and characterized a role for the chromatin remodeller CHD2 in myeloid proliferation and megakaryocytic differentiation.

A network view of microRNA and gene interactions in different pathological stages of colon cancer

BACKGROUND

Colon cancer is one of the common cancers in human. Although the number of annual cases has decreased drastically, prognostic screening and translational methods can be improved. Hence, it is critical to understand the molecular mechanisms of disease progression and prognosis.

RESULTS

In this study, we develop a new strategy for integrating microRNA and gene expression profiles together with clinical information toward understanding the regulation of colon cancer. Particularly, we use this approach to identify microRNA and gene expression networks that are specific to certain pathological stages. To demonstrate the application of our method, we apply this approach to identify microRNA and gene interactions that are specific to pathological stages of colon cancer in The Cancer Genome Atlas (TCGA) datasets.

CONCLUSIONS

Our results show that there are significant differences in network connections between miRNAs and genes in different pathological stages of colon cancer. These findings point to a hypothesis that these networks signify different roles of microRNA and gene regulation in the pathogenesis and tumorigenesis of colon cancer.

Proteomic mapping by rapamycin-dependent targeting of APEX2 identifies binding partners of VAPB at the inner nuclear membrane

Vesicle-associated membrane protein-associated protein B (VAPB) is a tail-anchored protein that is present at several contact sites of the endoplasmic reticulum (ER). We now show by immunoelectron microscopy that VAPB also localizes to the inner nuclear membrane (INM). Using a modified enhanced ascorbate peroxidase 2 (APEX2) approach with rapamycin-dependent targeting of the peroxidase to a protein of interest, we searched for proteins that are in close proximity to VAPB, particularly at the INM. In combination with stable isotope labeling with amino acids in cell culture (SILAC), we confirmed many well-known interaction partners at the level of the ER with a clear distinction between specific and nonspecific hits. Furthermore, we identified emerin, TMEM43, and ELYS as potential interaction partners of VAPB at the INM and the nuclear pore complex, respectively.

Resonance assignments for the tandem PWWP-ARID domains of human RBBP1

Retinoblastoma-binding protein 1 (RBBP1), also known as AT-rich interaction domain 4A (ARID4A), is a tumour suppressor involved in the regulation of the epigenetic programming in leukemia and Prader-Willi/Angelman syndromes. The ARID domain of RBBP1 binds to DNA non-specifically and has gene suppression activity. However, no structural data has been obtained for the human RBBP1 ARID domain so far. Here we report the near-complete ¹H, ¹³C, ¹⁵N backbone and side-chain NMR assignment of a 27 kDa tandem PWWP-ARID domain construct that spans residues 171-414 with the removal of a short disordered region between the two domains. The predicted secondary structure based on the assigned chemical shifts is consistent with the structures of the isolated PWWP domain of human RBBP1 previously solved and the homologous ARID domains of other proteins.

Determination of genetic aberrations and novel transcripts involved in the pathogenesis of oligodendroglioma using array comparative genomic hybridization and next generation sequencing

The aim of the present study was to determine the genetic aberrations and novel transcripts, particularly the fusion transcripts, involved in the pathogenesis of low-grade and anaplastic oligodendroglioma. In the present study, tissue samples were obtained from patients with oligodendroglioma and additionally from archived tissue samples from the Brain Tumor Tissue Bank of the Brain Tumor Foundation of Canada. Six samples were obtained, three of which were low-grade oligodendroglioma and the other three anaplastic oligodendroglioma. DNA and RNA were extracted from each tissue sample. The resulting genomic DNA was then hybridized using the Agilent CytoSure 4×180K oligonucleotide array. Human reference DNA and samples were labeled using Cy3 cytidine 5′-triphosphate (CTP) and Cy5 CTP, respectively, while human Cot-1 DNA was used to reduce non-specific binding. Microarray-based comparative genomic hybridization data was then analyzed for genetic aberrations using the Agilent Cytosure Interpret software v3.4.2. The total RNA isolated from each sample was mixed with oligo dT magnetic beads to enrich for poly(A) mRNA. cDNAs were then synthesized and subjected to end-repair, poly(A) addition and connected using sequencing adapters using the Illumina TruSeq RNA Sample Preparation kit. The fragments were then purified and selected as templates for polymerase chain reaction amplification. The final library was constructed with fragments between 350–450 base pairs and sequenced using deep transcriptome sequencing on an Illumina HiSeq 2500 sequencer. The array comparative genomic hybridization revealed numerous amplifications and deletions on several chromosomes in all samples. However, the most interesting result was from the next generation sequencing, where one anaplastic oligodendroglioma sample was demonstrated to have five novel fusion genes that may potentially serve a critical role in tumor pathogenesis and progression.

Downregulation of ARID4A and ARID4B promote tumor progression and directly regulated by microRNA-30d in patient with prostate cancer

AT-rich interaction domain 4A (ARID4A) and AT-rich interaction domain 4B (ARID4B), which are both the AT-rich interaction domain (ARID) family, have been reported to be oncogene or tumor suppressor gene in various human malignances, but there is no involvement about their functions in prostate cancer (PCa). Our previous study has reported that microRNA-30d (miR-30d) expression can predicted poor clinical prognosis in PCa, however, the underlying mechanisms of miR-30d have not been fully described. The aim of our study is to investigate the expression relevance between miR-30d and ARID4A or ARID4B, and examine the clinical significance and biological function of ARID4A and AIRD4B in PCa. In this study, both ARID4A and ARID4B were identified as the target genes of miR-30d. In addition, the mRNA expression of miR-30d in PCa tissues were significantly negative correlated with ARID4A (Pearson correlation coefficient = -0.313, P = 0.001) and ARID4B (Pearson correlation coefficient = -0.349, P < 0.001), while there was a positive correlation between ARID4A and ARID4B (Pearson correlation coefficient = 0.865, P < 0.001). Moreover, both ARID4A and ARID4B were significantly downregulated in PCa tissues with high Gleason scores (P = 0.005, P = 0.033), PSA failure (P = 0.012, P = 0.05) and short biochemical recurrent-free survival (P = 0.033, P = 0.031). Furthermore, the knockout expression of ARID4A and ARID4B promoted PCa cell proliferation, migration and invasion in vitro. In conclusion, our results indicated that ARID4A and ARID4B may serve as tumor suppressor in PCa progression, suggesting that they might be the potential therapeutic targets in prostate cancer.

A fast exact functional test for directional association and cancer biology applications

Directional association measured by functional dependency can answer important questions on relationships between variables, for example, in discovery of molecular interactions in biological systems. However, when one has no prior information about the functional form of a directional association, there is not a widely established statistical procedure to detect such an association. To address this issue, here we introduce an exact functional test for directional association by examining the strength of functional dependency. It is effective in promoting functional patterns by reducing statistical power on non-functional patterns. We designed an algorithm to carry out the test using a fast branch-and-bound strategy, which achieved a substantial speedup over brute-force enumeration. On data from an epidemiological study of liver cancer, the test identified the hepatitis status of a subject as the most influential risk factor among others for the cancer phenotype. On human lung cancer transcriptome data, the test selected 1049 transcription start sites of putative noncoding RNAs directionally associated with lung cancers, stronger than 95% of 589 curated cancer genes. These predictions include non-monotonic interaction patterns, to which other routine tests were insensitive. Complementing symmetric (non-directional) association methods such as Fisher's exact test, the exact functional test is a unique exact statistical test for evaluating evidence for causal relationships.

Overexpression of ARID4B predicts poor survival in patients with hepatocellular carcinoma

AT-rich interaction domain 4B (ARID4B), which belongs to the ARID family, is heavily involved in cell growth and differentiation and is closely associated with many types of tumors. However, the role of this protein in hepatocellular carcinoma (HCC) remains unknown. In this study, we used data from The Cancer Genome Atlas and Gene Expression Omnibus to analyze ARID4B expression in HCC. We subjected 15 pairs of fresh-frozen tissue samples to quantitative real-time polymerase chain reaction and Western blotting analyses to investigate ARID4B expression. We also subjected 157 formalin-fixed, paraffin-embedded HCC tissue samples to immunohistochemical analysis to detect ARID4B expression and to determine the clinical significance of ARID4B expression in HCC. The bioinformatics analysis, quantitative real-time polymerase chain reaction, and Western blotting results showed that ARID4B was highly expressed in HCC tissues compared with adjacent normal liver tissues. High ARID4B expression was strongly correlated with tumor number (P = .02), vascular invasion (P = .004), Edmondson-Steiner grades (P = .000), and tumor-node-metastasis stages (P = .001). Moreover, Kaplan-Meier and Cox proportional-hazards analyses indicated that high ARID4B expression was significantly associated with poor survival in patients with HCC and that ARID4B was an independent prognostic factor for overall survival and disease-free survival in patients with HCC. In conclusion, our results suggest that ARID4B acts as an oncogene in HCC and can therefore serve as a biomarker for the prognoses of patients with HCC.

CALR mutational status identifies different disease subtypes of essential thrombocythemia showing distinct expression profiles

Polycythemia vera (PV) and essential thrombocythemia (ET) are Philadelphia-negative myeloproliferative neoplasms (MPNs) characterized by erythrocytosis and thrombocytosis, respectively. Approximately 95% of PV and 50–70% of ET patients harbor the V617F mutation in the exon 14 of JAK2 gene, while about 20–30% of ET patients carry CALRins5 or CALRdel52 mutations. These ET CALR-mutated subjects show higher platelet count and lower thrombotic risk compared to JAK2-mutated patients. Here, we showed that CALR-mutated and JAK2V617F-positive CD34+ cells display different gene and miRNA expression profiles. Indeed, we highlighted several pathways differentially activated between JAK2V617F- and CALR-mutated progenitors, i.e., mTOR, MAPK/PI3K, and MYC pathways. Furthermore, we unveiled that the expression of several genes involved in DNA repair, chromatin remodeling, splicing, and chromatid cohesion are decreased in CALR-mutated cells. According to the low risk of thrombosis in CALR-mutated patients, we also found the downregulation of several genes involved in thrombin signaling and platelet activation. As a whole, these data support the model that CALR-mutated ET could be considered as a distinct disease entity from JAK2V617F-positive MPNs and may provide the molecular basis supporting the different clinical features of these patients.

miR-494-3p overexpression promotes megakaryocytopoiesis in primary myelofibrosis hematopoietic stem/progenitor cells by targeting SOCS6

Primary myelofibrosis (PMF) is a chronic Philadelphia-negative myeloproliferative neoplasm characterized by hematopoietic stem cell-derived clonal myeloproliferation, involving especially the megakaryocyte lineage. To better characterize how the altered expression of microRNAs might contribute to PMF pathogenesis, we have previously performed the integrative analysis of gene and microRNA expression profiles of PMF hematopoietic stem/progenitor cells (HSPCs), which allowed us to identify miR-494-3p as the upregulated microRNA predicted to target the highest number of downregulated mRNAs.To elucidate the role of miR-494-3p in hematopoietic differentiation, in the present study we demonstrated that miR-494-3p enforced expression in normal HSPCs promotes megakaryocytopoiesis. Gene expression profiling upon miR-494-3p overexpression allowed the identification of genes commonly downregulated both after microRNA overexpression and in PMF CD34+ cells. Among them, suppressor of cytokine signaling 6 (SOCS6) was confirmed to be a miR-494-3p target by luciferase assay. Western blot analysis showed reduced level of SOCS6 protein as well as STAT3 activation in miR-494-3p overexpressing cells. Furthermore, transient inhibition of SOCS6 expression in HSPCs demonstrated that SOCS6 silencing stimulates megakaryocytopoiesis, mimicking the phenotypic effects observed upon miR-494-3p overexpression. Finally, to disclose the contribution of miR-494-3p upregulation to PMF pathogenesis, we performed inhibition experiments in PMF HSPCs, which showed that miR-494-3p silencing led to SOCS6 upregulation and impaired megakaryocyte differentiation.Taken together, our results describe for the first time the role of miR-494-3p during normal HSPC differentiation and suggest that its increased expression, and the subsequent downregulation of its target SOCS6, might contribute to the megakaryocyte hyperplasia commonly observed in PMF patients.

Appraising the relevance of DNA copy number loss and gain in prostate cancer using whole genome DNA sequence data

A variety of models have been proposed to explain regions of recurrent somatic copy number alteration (SCNA) in human cancer. Our study employs Whole Genome DNA Sequence (WGS) data from tumor samples (n = 103) to comprehensively assess the role of the Knudson two hit genetic model in SCNA generation in prostate cancer. 64 recurrent regions of loss and gain were detected, of which 28 were novel, including regions of loss with more than 15% frequency at Chr4p15.2-p15.1 (15.53%), Chr6q27 (16.50%) and Chr18q12.3 (17.48%). Comprehensive mutation screens of genes, lincRNA encoding sequences, control regions and conserved domains within SCNAs demonstrated that a two-hit genetic model was supported in only a minor proportion of recurrent SCNA losses examined (15/40). We found that recurrent breakpoints and regions of inversion often occur within Knudson model SCNAs, leading to the identification of ZNF292 as a target gene for the deletion at 6q14.3-q15 and NKX3.1 as a two-hit target at 8p21.3-p21.2. The importance of alterations of lincRNA sequences was illustrated by the identification of a novel mutational hotspot at the KCCAT42, FENDRR, CAT1886 and STCAT2 loci at the 16q23.1-q24.3 loss. Our data confirm that the burden of SCNAs is predictive of biochemical recurrence, define nine individual regions that are associated with relapse, and highlight the possible importance of ion channel and G-protein coupled-receptor (GPCR) pathways in cancer development. We concluded that a two-hit genetic model accounts for about one third of SCNA indicating that mechanisms, such haploinsufficiency and epigenetic inactivation, account for the remaining SCNA losses.

Turning the tide in myelodysplastic/myeloproliferative neoplasms

Myelodysplastic syndromes/myeloproliferative neoplasms (MDS/MPN) are aggressive myeloid malignancies recognized as a distinct category owing to their unique combination of dysplastic and proliferative features. Although current classification schemes still emphasize morphology and exclusionary criteria, disease-defining somatic mutations and/or germline predisposition alleles are increasingly incorporated into diagnostic algorithms. The developing picture suggests that phenotypes are driven mostly by epigenetic mechanisms that reflect a complex interplay between genotype, physiological processes such as ageing and interactions between malignant haematopoietic cells and the stromal microenvironment of the bone marrow. Despite the rapid accumulation of genetic knowledge, therapies have remained nonspecific and largely inefficient. In this Review, we discuss the pathogenesis of MDS/MPN, focusing on the relationship between genotype and phenotype and the molecular underpinnings of epigenetic dysregulation. Starting with the limitations of current therapies, we also explore how the available mechanistic data may be harnessed to inform strategies to develop rational and more effective treatments, and which gaps in our knowledge need to be filled to translate biological understanding into clinical progress.

miR-376c promotes carcinogenesis and serves as a plasma marker for gastric carcinoma

Gastric carcinoma is highly prevalent throughout the world. Understanding the pathogenesis of this disease will benefit diagnosis and resolution. Studies show that miRNAs are involved in the tumorigenesis of gastric carcinoma. An initial screening followed by subsequent validation identified that miR-376c is up-regulated in gastric carcinoma tissue and the plasma of patients with the disease. In addition, the urinary level of miR-376c is also significantly increased in gastric carcinoma patients. The plasma miR-376c level was validated as a biomarker for gastric carcinoma, including early stage tumors. The induction of miR-376c was found to enrich the proliferation, migration and anchorage-independent growth of carcinoma cells and, furthermore, the repression of the expression of endogenous miR-376c was able to reduce such oncogenic phenotypes. ARID4A gene is a direct target of miR-376c. Knockdown of endogenous ARID4A increased the oncogenicity of carcinoma cells, while ARID4A was found to be drastically down-regulated in tumor tissue. Thus, expression levels of miR-376c and ARID4A mRNA tended to be opposing in tumor tissue. Our results demonstrate that miR-376c functions by suppressing ARID4A expression, which in turn enhances the oncogenicity of gastric carcinoma cells. It seems likely that the level of miR-376c in plasma and urine could act as invaluable markers for the detection of gastric carcinoma.

Differential Sensitivity of Target Genes to Translational Repression by miR-17~92

MicroRNAs (miRNAs) are thought to exert their functions by modulating the expression of hundreds of target genes and each to a small degree, but it remains unclear how small changes in hundreds of target genes are translated into the specific function of a miRNA. Here, we conducted an integrated analysis of transcriptome and translatome of primary B cells from mutant mice expressing miR-17~92 at three different levels to address this issue. We found that target genes exhibit differential sensitivity to miRNA suppression and that only a small fraction of target genes are actually suppressed by a given concentration of miRNA under physiological conditions. Transgenic expression and deletion of the same miRNA gene regulate largely distinct sets of target genes. miR-17~92 controls target gene expression mainly through translational repression and 5’UTR plays an important role in regulating target gene sensitivity to miRNA suppression. These findings provide molecular insights into a model in which miRNAs exert their specific functions through a small number of key target genes.

MicroRNAs (miRNAs) are small RNAs encoded by our genome. Each miRNA binds hundreds of target mRNAs and performs specific functions. It is thought that miRNAs exert their function by reducing the expression of all these target genes and each to a small degree. However, these target genes often have very diverse functions. It has been unclear how small changes in hundreds of target genes with diverse functions are translated into the specific function of a miRNA. Here we take advantage of recent technical advances to globally examine the mRNA and protein levels of 868 target genes regulated by miR-17~92, the first oncogenic miRNA, in mutant mice with transgenic overexpression or deletion of this miRNA gene. We show that miR-17~92 regulates target gene expression mainly at the protein level, with little effect on mRNA. Surprisingly, only a small fraction of target genes respond to miR-17~92 expression changes. Further studies show that the sensitivity of target genes to miR-17~92 is determined by a non-coding region of target mRNA. Our findings demonstrate that not every target gene is equal, and suggest that the function of a miRNA is mediated by a small number of key target genes.

ARID4B is a good biomarker to predict tumour behaviour and decide WHO grades in gliomas and meningiomas

Aims

Although ARID4B is known to promote tumour metastasis in breast cancer and inhibit transformation and progression in leukaemia, the possible effect of ARID4B on primary brain tumours (PBTs) is not well characterised. We tested the hypothesis that expression of ARID4B correlates with WHO grade and survival in patients with PBTs.

Methods

Western blot analysis was performed on protein lysates prepared from normal brain tissue and glioma cell lines (U87MG, LN229, GBM8401 and U118MG). Subsequently, immunohistochemical analysis of ARID4B was performed on 2 tissue microarrays, including 12 normal brain tissues, 63 meningiomas with different subtypes, 232 gliomas of various grades and degrees of differentiation, 8 central neurocytomas and 4 chordomas. The ARID4B immunostaining score was calculated by multiplying the intensity score by the percentage of tumour cells expressing ARID4B.

Results

In vitro, ARID4B protein expression was increased in some glioma cell lines. In addition, the average ARID4B immunostaining score was 38.03, 79.09, 129.76 and 119.32, respectively, in gliomas of WHO grade I, II, III and IV. Higher ARID4B immunostaining score was significantly correlated with more advanced WHO grade of gliomas (p=7.4×10–6) and meningiomas. Finally, higher ARID4B expression tended to the shorter survival rates, but did not reach statistical significance.

Conclusions

ARID4B overexpression presented in most of PBTs, rather than non-neoplastic brain tissue, and correlated with WHO grades in meningiomas and gliomas. Therefore, ARID4B is a satisfactory biomarker to highlight tumour component and predict tumour behaviour in primary brain neoplasms.

Screening of biomarkers for prediction of response to and prognosis after chemotherapy for breast cancers

OBJECTIVE

To screen the biomarkers having the ability to predict prognosis after chemotherapy for breast cancers.

METHODS

Three microarray data of breast cancer patients undergoing chemotherapy were collected from Gene Expression Omnibus database. After preprocessing, data in GSE41112 were analyzed using significance analysis of microarrays to screen the differentially expressed genes (DEGs). The DEGs were further analyzed by Differentially Coexpressed Genes and Links to construct a function module, the prognosis efficacy of which was verified by the other two datasets (GSE22226 and GSE58644) using Kaplan-Meier plots. The involved genes in function module were subjected to a univariate Cox regression analysis to confirm whether the expression of each prognostic gene was associated with survival.

RESULTS

A total of 511 DEGs between breast cancer patients who received chemotherapy or not were obtained, consisting of 421 upregulated and 90 downregulated genes. Using the Differentially Coexpressed Genes and Links package, 1,244 differentially coexpressed genes (DCGs) were identified, among which 36 DCGs were regulated by the transcription factor complex NFY (NFYA, NFYB, NFYC). These 39 genes constructed a gene module to classify the samples in GSE22226 and GSE58644 into three subtypes and these subtypes exhibited significantly different survival rates. Furthermore, several genes of the 39 DCGs were shown to be significantly associated with good (such as CDC20) and poor (such as ARID4A) prognoses following chemotherapy.

CONCLUSION

Our present study provided a serial of biomarkers for predicting the prognosis of chemotherapy or targets for development of alternative treatment (ie, CDC20 and ARID4A) in breast cancer patients.

Insights into myelodysplastic syndromes from current preclinical models

In recent years, there has been significant progress made in our understanding of the molecular genetics of myelodysplastic syndromes (MDS). Using massively parallel sequencing techniques, recurring mutations are identified in up to 80% of MDS cases, including many with a normal karyotype. The differential role of some of these mutations in the initiation and progression of MDS is starting to be elucidated. Engineering candidate genes in mice to model MDS has contributed to recent insights into this complex disease. In this review, we examine currently available mouse models, with detailed discussion of selected models. Finally, we highlight some advances made in our understanding of MDS biology, and conclude with discussions of questions that remain unanswered.

Pediatric T-cell lymphoblastic leukemia evolves into relapse by clonal selection, acquisition of mutations and promoter hypomethylation

Relapsed precursor T-cell acute lymphoblastic leukemia is characterized by resistance against chemotherapy and is frequently fatal. We aimed at understanding the molecular mechanisms resulting in relapse of T-cell acute lymphoblastic leukemia and analyzed 13 patients at first diagnosis, remission and relapse by whole exome sequencing, targeted ultra-deep sequencing, multiplex ligation dependent probe amplification and DNA methylation array. Compared to primary T-cell acute lymphoblastic leukemia, in relapse the number of single nucleotide variants and small insertions and deletions approximately doubled from 11.5 to 26. Targeted ultra-deep sequencing sensitively detected subclones that were selected for in relapse. The mutational pattern defined two types of relapses. While both are characterized by selection of subclones and acquisition of novel mutations, 'type 1' relapse derives from the primary leukemia whereas 'type 2' relapse originates from a common pre-leukemic ancestor. Relapse-specific changes included activation of the nucleotidase NT5C2 resulting in resistance to chemotherapy and mutations of epigenetic modulators, exemplified by SUZ12, WHSC1 and SMARCA4. While mutations present in primary leukemia and in relapse were enriched for known drivers of leukemia, relapse-specific changes revealed an association with general cancer-promoting mechanisms. This study thus identifies mechanisms that drive progression of pediatric T-cell acute lymphoblastic leukemia to relapse and may explain the characteristic treatment resistance of this condition.

Androgen Receptor Coactivator ARID4B Is Required for the Function of Sertoli Cells in Spermatogenesis

Defects in spermatogenesis, a process that produces spermatozoa inside seminiferous tubules of the testis, result in male infertility. Spermatogenic progression is highly dependent on a microenvironment provided by Sertoli cells, the only somatic cells and epithelium of seminiferous tubules. However, genes that regulate such an important activity of Sertoli cells are poorly understood. Here, we found that AT-rich interactive domain 4B (ARID4B), is essential for the function of Sertoli cells to regulate spermatogenesis. Specifically, we generated Sertoli cell-specific Arid4b knockout (Arid4bSCKO) mice, and showed that the Arid4bSCKO male mice were completely infertile with impaired testis development and significantly reduced testis size. Importantly, severe structural defects accompanied by loss of germ cells and Sertoli cell-only phenotype were found in many seminiferous tubules of the Arid4bSCKO testes. In addition, maturation of Sertoli cells was significantly delayed in the Arid4bSCKO mice, associated with delayed onset of spermatogenesis. Spermatogenic progression was also defective, showing an arrest at the round spermatid stage in the Arid4bSCKO testes. Interestingly, we showed that ARID4B functions as a "coactivator" of androgen receptor and is required for optimal transcriptional activation of reproductive homeobox 5, an androgen receptor target gene specifically expressed in Sertoli cells and critical for spermatogenesis. Together, our study identified ARID4B to be a key regulator of Sertoli cell function important for male germ cell development.

Blood and immune cell engineering: Cytoskeletal contractility and nuclear rheology impact cell lineage and localization: Biophysical regulation of hematopoietic differentiation and trafficking

Clinical success with human hematopoietic stem cell (HSC) transplantation establishes a paradigm for regenerative therapies with other types of stem cells. However, it remains generally challenging to therapeutically treat tissues after engineering of stem cells in vitro. Recent studies suggest that stem and progenitor cells sense physical features of their niches. Here, we review biophysical contributions to lineage decisions, maturation, and trafficking of blood and immune cells. Polarized cellular contractility and nuclear rheology are separately shown to be functional markers of a hematopoietic hierarchy that predict the ability of a lineage to traffic in and out of the bone marrow niche. These biophysical determinants are regulated by a set of structural molecules, including cytoplasmic myosin-II and nuclear lamins, which themselves are modulated by a diverse range of transcriptional and post-translational mechanisms. Small molecules that target these mechanobiological circuits, along with novel bioengineering methods, could prove broadly useful in programming blood and immune cells for therapies ranging from blood transfusions to immune attack of tumors.

Randomness and preserved patterns in cancer network

Breast cancer has been reported to account for the maximum cases among all female cancers till date. In order to gain a deeper insight into the complexities of the disease, we analyze the breast cancer network and its normal counterpart at the proteomic level. While the short range correlations in the eigenvalues exhibiting universality provide an evidence towards the importance of random connections in the underlying networks, the long range correlations along with the localization properties reveal insightful structural patterns involving functionally important proteins. The analysis provides a benchmark for designing drugs which can target a subgraph instead of individual proteins.

Recent advances in the ARID family: focusing on roles in human cancer

The human AT-rich interaction domain (ARID) family contains seven subfamilies and 15 members characterized by having an ARID. Members of the ARID family have the ability to regulate transcription and are involved in cell differentiation and proliferation. Accumulating evidence suggests that ARID family members are involved in cancer-related signaling pathways, highly mutated or differentially expressed in tumor tissues, and act as predictive factors for cancer prognosis or therapeutic outcome. Here we review the molecular biology and clinical studies concerned with the role played by the ARID family in cancer. This may contribute to our understanding of the initiation and progression of cancer from a novel point of view, as well as providing potential targets for cancer therapy.

Exome sequencing reveals frequent inactivating mutations in ARID1A, ARID1B, ARID2 and ARID4A in microsatellite unstable colorectal cancer

ARID1A has been identified as a novel tumor suppressor gene in ovarian cancer and subsequently in various other tumor types. ARID1A belongs to the ARID domain containing gene family, which comprises of 15 genes involved, for example, in transcriptional regulation, proliferation and chromatin remodeling. In this study, we used exome sequencing data to analyze the mutation frequency of all the ARID domain containing genes in 25 microsatellite unstable (MSI) colorectal cancers (CRCs) as a first systematic effort to characterize the mutation pattern of the whole ARID gene family. Genes which fulfilled the selection criteria in this discovery set (mutations in at least 4/25 [16%] samples, including at least one nonsense or splice site mutation) were chosen for further analysis in an independent validation set of 21 MSI CRCs. We found that in addition to ARID1A, which was mutated in 39% of the tumors (18/46), also ARID1B (13%, 6/46), ARID2 (13%, 6/46) and ARID4A (20%, 9/46) were frequently mutated. In all these genes, the mutations were distributed along the entire length of the gene, thus distinguishing them from typical MSI target genes previously described. Our results indicate that in addition to ARID1A, other members of the ARID gene family may play a role in MSI CRC.

Retinoblastoma-binding protein 1 has an interdigitated double Tudor domain with DNA binding activity

Retinoblastoma-binding protein 1 (RBBP1) is a tumor and leukemia suppressor that binds both methylated histone tails and DNA. Our previous studies indicated that RBBP1 possesses a Tudor domain, which cannot bind histone marks. In order to clarify the function of the Tudor domain, the solution structure of the RBBP1 Tudor domain was determined by NMR and is presented here. Although the proteins are unrelated, the RBBP1 Tudor domain forms an interdigitated double Tudor structure similar to the Tudor domain of JMJD2A, which is an epigenetic mark reader. This indicates the functional diversity of Tudor domains. The RBBP1 Tudor domain structure has a significant area of positively charged surface, which reveals a capability of the RBBP1 Tudor domain to bind nucleic acids. NMR titration and isothermal titration calorimetry experiments indicate that the RBBP1 Tudor domain binds both double- and single-stranded DNA with an affinity of 10-100 μM; no apparent DNA sequence specificity was detected. The DNA binding mode and key interaction residues were analyzed in detail based on a model structure of the Tudor domain-dsDNA complex, built by HADDOCK docking using the NMR data. Electrostatic interactions mediate the binding of the Tudor domain with DNA, which is consistent with NMR experiments performed at high salt concentration. The DNA-binding residues are conserved in Tudor domains of the RBBP1 protein family, resulting in conservation of the DNA-binding function in the RBBP1 Tudor domains. Our results provide further insights into the structure and function of RBBP1.