Distal Promoter Research Articles

Enhancer-Insulator Pairing Reveals Heterogeneous Dynamics in Long-Distance 3D Gene Regulation

Cells regulate fates and complex body plans using spatiotemporal signaling cascades that alter gene expression. Short DNA sequences, known as enhancers (50–1500 base pairs), help coordinate these cascades by attracting regulatory proteins that enhance the transcription by binding to distal gene promoters. In humans, there are hundreds of thousands of enhancers dispersed across the genome, which poses a challenging coordination task to prevent unintended gene activation. To mitigate this problem, the genome contains insulator elements that block enhancer-promoter interactions. However, there is an open problem with how the insulation works, especially as enhancer-insulator pairs may be separated by millions of base pairs. Based on recent empirical data from Hi-C experiments, this paper proposes a new mechanism that challenges the common paradigm that rests on specific insulator-insulator interactions. Instead, this paper introduces a stochastic looping model where insulators bind weakly to chromatin rather than other insulators. After calibrating the model to experimental data, we use simulations to study the broad distribution of hitting times between an enhancer and a promoter when insulators are present. We find parameter regimes with large differences between average and most probable hitting times. This makes it difficult to assign a typical timescale and hints at highly defocused regulation times. We also map our computational model onto a resetting problem that allows us to derive several analytical results. Besides offering new insights into enhancer-insulator interactions, our paper advances the understanding of gene regulatory networks and causal connections between genome folding and gene activation. Published by the American Physical Society 2024

RUNX1 isoforms regulate RUNX1 and target genes differentially in platelets-megakaryocytes: association with clinical cardiovascular events

BackgroundHematopoietic transcription factor RUNX1 is expressed from proximal P2 and distal P1 promoters to yield isoforms RUNX1 B and C, respectively. The roles of these isoforms in RUNX1 autoregulation and downstream gene regulation in megakaryocytes and platelets are unknown. ObjectivesTo understand the regulation of RUNX1 and its target genes by RUNX1 isoforms. MethodsWe performed studies on RUNX1 isoforms in megakaryocytic human erythroleukemia (HEL) cells and HeLa cells (lack endogenous RUNX1), in platelets from 85 healthy volunteers administered aspirin or ticagrelor, and on the association of RUNX1 target genes with acute events in 587 patients with cardiovascular disease (CVD). ResultsIn chromatin immunoprecipitation and luciferase promoter assays, RUNX1 isoforms B and C bound and regulated P1 and P2 promoters. In HeLa cells, RUNX1B decreased and RUNX1C increased P1 and P2 activities, respectively. In HEL cells, RUNX1B overexpression decreased RUNX1C and RUNX1A expression; RUNX1C increased RUNX1B and RUNX1A. RUNX1B and RUNX1C regulated target genes (MYL9, F13A1, PCTP, PDE5A, and others) differentially in HEL cells. In platelets, RUNX1B transcripts (by RNA sequencing) correlated negatively with RUNX1C and RUNX1A; RUNX1C correlated positively with RUNX1A. RUNX1B correlated positively with F13A1, PCTP, PDE5A, RAB1B, and others, and negatively with MYL9. In our previous studies, RUNX1C transcripts in whole blood were protective against acute events in CVD patients. We found that higher expression of RUNX1 targets F13A1 and RAB31 associated with acute events. ConclusionRUNX1 isoforms B and C autoregulate RUNX1 and regulate downstream genes in a differential manner, and this is associated with acute events in CVD.

The xenobiotic transporter ABCC4/MRP4 promotes epithelial mesenchymal transition in pancreatic cancer.

The xenobiotic transporter ABCC4/MRP4 is highly expressed in pancreatic ductal adenocarcinoma (PDAC) and correlates with a more aggressive phenotype and metastatic propensity. Here, we show that ABCC4 promotes epithelial-mesenchymal transition (EMT) in PDAC, a hallmark process involving the acquisition of mesenchymal traits by epithelial cells, enhanced cell motility, and chemoresistance. Modulation of ABCC4 levels in PANC-1 and BxPC-3 cell lines resulted in the dysregulation of genes present in the EMT signature. Bioinformatic analysis on several cohorts including tumor samples, primary patient-derived cultured cells, patient-derived xenografts, and cell lines, revealed a positive correlation between ABCC4 expression and EMT markers. We also characterized the ABCC4 cistrome and identified four candidate clusters in the distal promoter and intron one that showed differential binding of pro-epithelial FOXA1 and pro-mesenchymal GATA2 transcription factors in low ABCC4-expressing HPAF-II and high ABCC4-expressing PANC-1 xenografts. HPAF-II xenografts showed exclusive binding of FOXA1, and PANC-1 xenografts exclusive binding of GATA2, at ABCC4 clusters, consistent with their low and high EMT phenotype respectively. Our results underscore ABCC4/MRP4 as a valuable prognostic marker and a potential therapeutic target to treat PDAC subtypes with prominent EMT features, such as the basal-like/squamous subtype, characterized by worse prognosis and no effective therapies.

Long range regulation of transcription scales with genomic distance in a gene specific manner.

While critical for tuning the timing and level of transcription, enhancer communication with distal promoters is not well understood. Here we bypass the need for sequence-specific transcription factors and recruit activators directly using CARGO-VPR, an approach for targeting dCas9-VPR using a multiplexed array of RNA guides. We show that this approach achieves effective activator recruitment to arbitrary genomic sites, even those inaccessible by single dCas9. We utilize CARGO-VPR across the Prdm8-Fgf5 locus in mESCs, where neither gene is expressed. We demonstrate that while activator recruitment to any tested region results in transcriptional induction of at least one gene, the expression level strongly depends on the genomic distance between the promoter and activator recruitment site. However, the expression-distance relationship for each gene scales distinctly in a manner not attributable to differences in 3D contact frequency, promoter DNA sequence or presence of the repressive chromatin marks at the locus.

The Distal Promoter of the B438L Gene of African Swine Fever Virus Is Responsible for the Transcription of the Alternatively Spliced B169L.

The B169L protein (pB169L) of African swine fever virus (ASFV) is a structural protein with an unidentified function during the virus replication. The sequences of the B169L gene and the downstream B438L gene are separated by short intergenic regions. However, the regulatory mode of the gene transcription remains unknown. Here, we identified two distinct promoter regions and two transcription start sites (TSSs) located upstream of the open reading frame (ORF) of B438L. Using the promoter reporter system, we demonstrated that the cis activity of the ORF proximal promoter exhibited significantly higher levels compared with that of the distal promoter located in the B169L gene. Furthermore, transfection with the plasmids with two different promoters for B438L could initiate the transcription and expression of the B438L gene in HEK293T cells, and the cis activity of the ORF proximal promoter also displayed higher activities compared with the distal promoter. Interestingly, the B438L distal promoter also initiated the transcription of the alternatively spliced B169L mRNA (B169L mRNA2) encoding a truncated pB169L (tpB169L) (amino acids 92-169), and the gene transcription efficiency was increased upon mutation of the initiation codon located upstream of the alternatively spliced B169L gene. Taken together, we demonstrated that the distal promoter of B438L gene initiates the transcription of both the B438L mRNA and B169L mRNA2. Comprehensive analysis of the transcriptional regulatory mode of the B438L gene is beneficial for the understanding of the association of B438L protein and pB169L and the construction of the gene-deleted ASFV.

RUNX1 Isoforms Regulate RUNX1 and Target-Genes Differentially in Platelets-Megakaryocytes: Association with Clinical Cardiovascular Events.

Hematopoietic transcription factor RUNX1 is expressed from proximal P2 and distal P1 promoter to yield isoforms RUNX1 B and C, respectively. The roles of these isoforms in RUNX1 autoregulation and downstream-gene regulation in megakaryocytes and platelets are unknown. To understand the regulation of RUNX1 and its target genes by RUNX1 isoforms. We performed studies on RUNX1 isoforms in megakaryocytic HEL cells and HeLa cells (lack endogenous RUNX1), in platelets from 85 healthy volunteers administered aspirin or ticagrelor, and on the association of RUNX1 target genes with acute events in 587 patients with cardiovascular disease (CVD). In chromatin immunoprecipitation and luciferase promoter assays, RUNX1 isoforms B and C bound and regulated P1 and P2 promoters. In HeLa cells RUNX1B decreased and RUNX1C increased P1 and P2 activities, respectively. In HEL cells, RUNX1B overexpression decreased RUNX1C and RUNX1A expression; RUNX1C increased RUNX1B and RUNX1A. RUNX1B and RUNX1C regulated target genes (MYL9, F13A1, PCTP, PDE5A and others) differentially in HEL cells. In platelets RUNX1B transcripts (by RNAseq) correlated negatively with RUNX1C and RUNX1A; RUNX1C correlated positively with RUNX1A. RUNX1B correlated positively with F13A1, PCTP, PDE5A, RAB1B, and others, and negatively with MYL9. In our previous studies, RUNX1C transcripts in whole blood were protective against acute events in CVD patients. We found that higher expression of RUNX1 targets F13A1 and RAB31 associated with acute events. RUNX1 isoforms B and C autoregulate RUNX1 and regulate downstream genes in a differential manner and this associates with acute events in CVD.

ID3 is a novel target gene of p53 and modulates lung cancer cell metastasis

The tumor suppressor p53 prevents cancer development by regulating dozens of target genes with diverse biological functions. Although numerous p53 target genes have been identified to date, the dynamics and function of the regulatory network centered on p53 have not yet been fully elucidated. We herein identified inhibitor of DNA-binding/differentiation-3 (ID3) as a direct p53 target gene. p53 bound the distal promoter of ID3 and positively regulated its transcription. ID3 expression was significantly decreased in clinical lung cancer tissues, and was closely associated with overall survival outcomes in these patients. Functionally, ID3 deficiency promoted the metastatic ability of lung cancer cells through its effects on the transcriptional regulation of CDH1. Furthermore, the ectopic expression of ID3 in p53-knockdown cells restored E-cadherin expression. Collectively, the present results demonstrate that ID3 plays a tumor-suppressive role as a downstream effector of p53 and impedes lung cancer cell metastasis by regulating E-cadherin expression.

Transcriptional regulatory mechanism of NR2F2 and ZNF423 in avian preadipocyte differentiation

In the poultry industry, excessive abdominal fat deposition is not conducive to meat quality. Therefore, selection for optimal fat content levels in poultry has become a major breeding goal. We previously constructed NR2F2 overexpression (NR2F2OE) and knockout (NR2F2Δ/Δ/83-125aa) cell lines using Piggybac and CRISPR/Cas9 techniques, and confirmed that the transcription factor NR2F2 can significantly inhibit the differentiation of avian preadipocytes. In this study, we identified a downstream gene ZNF423 regulated by NR2F2, which is also involved in regulating avian fat deposition.First, we performed transcriptome analysis of the NR2F2-edited lines, which has been proven to be an inhibitor of avian fat deposition in our previous studies. Our findings revealed that NR2F2 affects a series of candidate regulators related to adipogenesis. Among these, we focused on ZNF423, which was significantly down-regulated in the NR2F2OE cell line and up-regulated in the NR2F2Δ/Δ/83-125aa cell line. Next, dual luciferase reporter assay results showed that the DNA-binding domain (DBDΔ72-143aa) of transcription factor NR2F2 may negatively affect the expression of downstream target gene ZNF423 by binding to its distal promoter region (-2356 to -2346). Moreover, we constructed a function analytical model and found that overexpression of ZNF423 significantly facilitated the differentiation of adipocytes in immortalized chicken preadipocytes (ICP1). Consistent with these findings, global transcriptome analysis of the ZNF423-overexpressed cell line (ZNF423OE) further demonstrated that the process of adipogenesis was significantly enriched. These results indicate that ZNF423 is a positive regulator of avian adipocyte differentiation. Overexpression of ZNF423 in the NR2F2OE cell line compensated for the inhibition of fat deposition phenotype, further suggesting that ZNF423 is a downstream target gene of NR2F2.These findings uncover a novel function of ZNF423 in avian adipocyte differentiation and analyzed the transcriptional regulation by its upstream transcription factor NR2F2. Additionally, we identified a list of functional candidate genes, providing important insights for further research on the mechanism of avian fat deposition.

Sustained MYB Activity Is Necessary for Oncogenic Transcription in KMT2A-Rearranged Acute Lymphoblastic Leukemia through Enhancer-Promoter Interactions and Epigenetic Modifications at Enhancers

Aberrant enhancer activity is a hallmark of many cancers including KMT2A-rearranged acute lymphoblastic leukemia (ALL), an aggressive leukemia subtype with an event-free survival of 19-45% in infant-ALL (<1 year of age) and 50-59% in childhood-ALL cases. Targeting enhancer activity may provide a novel therapeutic strategy for this hard-to-treat disease. Transcription factors are key effectors of enhancer activity. They recognise and bind to specific DNA sequences at enhancers, driving the expression of target genes from these distal loci. MYB is an important hematopoietic transcription factor that is commonly overexpressed in KMT2A-rearranged leukemias and is required for leukemia survival, where loss of MYB has been shown to abrogate the leukemia phenotype. However, the precise mechanism by which MYB contributes to enhancer function, and the stages of enhancer activation at which it acts are unknown. We employed a well-characterized chromatin anchoring system to interrogate MYB function and determine whether it can recruit specific enhancer activities de novo. By fusing the MYB transactivation domain (MYB TA) to the Tet repressor (TetR) DNA binding domain, MYB TA is directed to an array of Tet operator ( TetO) sequences incorporated into a neutral genomic region devoid of native gene regulatory elements. We found that MYB TA alone is sufficient to recruit key enhancer-associated proteins, including P300, Brd4 and Mediator to the TetO locus. We observed MYB TA-dependent deposition of H3K27ac, increased chromatin accessibility and transcription not just at the TetO locus, but remarkably, also at regions more than 50kb away. Using the chromosome conformation capture (3C) technique NG Capture-C, we observed MYB TA-induced DNA interactions between TetO and the distal sites of transcription, as well as more long-range interactions up to 400kb distal to where MYB TA is bound. Together, these observations suggest that MYB TA binding alone is sufficient to establish an enhancer-like regulatory element, which activates distant cryptic promoters. This activity is dependent on the continued presence of MYB TA, as disruption of MYB TA binding results in a rapid loss of transcription and acquired chromatin features, implicating MYB in enhancer maintenance as well as initiation. To confirm our observations in a disease context, we tested the requirement for MYB at endogenous enhancers in SEM cells, a model of KMT2A-rearranged ALL. By stably introducing a degradation tag (FKBP12 F36V) to the endogenous MYB gene, we were able to rapidly induce protein degradation and observe the immediate effects on chromatin features. Using transient transcriptome sequencing, we noted significant transcriptional downregulation of known MYB target genes, including BCL2, MYC, LMO4 and CDK6. We also observed decreases in enhancer H3K27ac levels and eRNA transcription, most marked at MYB-bound loci. To investigate the requirement for MYB in maintaining physical proximity between enhancer and promoter, we performed the base pair-resolution 3C technique Micro-Capture-C (MCC). MCC revealed highly punctate interactions between enhancer elements and their target promoters. At BCL2 and MYC, MYB degradation results in a reduction in interactions with the promoter at discrete regions within the enhancers. This argues for highly localised MYB-dependent activity at enhancers, driving interaction with and activation of distal target promoters. Overall, our data argue for a role for MYB in establishing enhancers de novo, indicating that MYB overexpression may directly result in onco-enhancer activity at key genes in vivo. At a subset of MYB-bound enhancers, the continued presence of MYB is absolutely required to maintain enhancer activity and oncogene upregulation. Further exploration of associated co-factors that mediate MYB-dependent enhancer activity may identify novel targets for therapeutic exploitation.

DHT and Insulin Upregulate Secretion of the Soluble Decoy Receptor of IL-33 From Decidualized Endometrial Stromal Cells.

Interleukin 33 (IL-33) signaling regulates most of the key processes of pregnancy, including decidualization, trophoblast proliferation and invasion, vascular remodeling, and placental growth. Accordingly, dysregulation of IL-33, its membrane-bound receptor (ST2L, transducer of IL-33 signaling), and its soluble decoy receptor (sST2, inhibitor of IL-33 signaling) has been linked to a wide range of adverse pregnancy outcomes that are common in women with obesity and polycystic ovary syndrome, that is, conditions associated with hyperandrogenism, insulin resistance, and compensatory hyperinsulinemia. To reveal if androgens and insulin might modulate uteroplacental IL-33 signaling, we investigated the effect of dihydrotestosterone (DHT) and/or insulin on the expression of ST2L and sST2 (along with the activity of their promoter regions), IL-33 and sIL1RAP (heterodimerization partner of sST2), during in vitro decidualization of endometrial stromal cells from 9 healthy women. DHT and insulin markedly upregulated sST2 secretion, in addition to the upregulation of its messenger RNA (mRNA) expression, while the proximal ST2 promoter, from which the sST2 transcript originates, was upregulated by insulin, and in a synergistic manner by DHT and insulin combination treatment. On the other hand, sIL1RAP was slightly downregulated by insulin and IL-33 mRNA expression was not affected by any of the hormones, while ST2L mRNA expression and transcription from its promoter region (distal ST2 promoter) could not be detected or showed a negligibly low level. We hypothesize that high levels of androgens and insulin might lead to subfertility and pregnancy complications, at least partially, through the sST2-dependent downregulation of uteroplacental IL-33 signaling.

Differential usage of DNA modifications in neurons, astrocytes, and microglia

BackgroundCellular identity is determined partly by cell type-specific epigenomic profiles that regulate gene expression. In neuroscience, there is a pressing need to isolate and characterize the epigenomes of specific CNS cell types in health and disease. In this study, we developed an in vivo tagging mouse model (Camk2a-NuTRAP) for paired isolation of neuronal DNA and RNA without cell sorting and then used this model to assess epigenomic regulation, DNA modifications in particular, of gene expression between neurons and glia.ResultsAfter validating the cell-specificity of the Camk2a-NuTRAP model, we performed TRAP-RNA-Seq and INTACT-whole genome oxidative bisulfite sequencing (WGoxBS) to assess the neuronal translatome and epigenome in the hippocampus of young mice (4 months old). WGoxBS findings were validated with enzymatic methyl-Seq (EM-Seq) and nanopore sequencing. Comparing neuronal data to microglial and astrocytic data from NuTRAP models, microglia had the highest global mCG levels followed by astrocytes and then neurons, with the opposite pattern observed for hmCG and mCH. Differentially modified regions between cell types were predominantly found within gene bodies and distal intergenic regions, rather than proximal promoters. Across cell types there was a negative correlation between DNA modifications (mCG, mCH, hmCG) and gene expression at proximal promoters. In contrast, a negative correlation of gene body mCG and a positive relationship between distal promoter and gene body hmCG with gene expression was observed. Furthermore, we identified a neuron-specific inverse relationship between mCH and gene expression across promoter and gene body regions.ConclusionsNeurons, astrocytes, and microglia demonstrate different genome-wide levels of mCG, hmCG, and mCH that are reproducible across analytical methods. However, modification-gene expression relationships are conserved across cell types. Enrichment of differential modifications across cell types in gene bodies and distal regulatory elements, but not proximal promoters, highlights epigenomic patterning in these regions as potentially greater determinants of cell identity. These findings also demonstrate the importance of differentiating between mC and hmC in neuroepigenomic analyses, as up to 30% of what is conventionally interpreted as mCG can be hmCG, which often has a different relationship to gene expression than mCG.

Mapping Mammalian 3D Genome Interactions with Micro-C-XL.

Three-dimensional (3D) chromosome organization is a major factor in genome regulation and cell-type specification. For example, cis-regulatory elements, known as enhancers, are thought to regulate the activity of distal promoters via interaction in 3D space. Genome-wide chromosome conformation capture (3C)-technologies, such as Hi-C, have transformed our understanding of how genomes are organized in cells. The current understanding of 3D genome organization is limited by the resolution with which the topological organization of chromosomes in 3D space can be resolved. Micro-C-XL measures chromosome folding with resolution at the level of the nucleosome, the basic unit of chromatin, by utilizing micrococcal nuclease (MNase) to fragment genomes during the chromosome conformation capture protocol. This results in an improved signal-to-noise ratio in the measurements, thus facilitating the better detection of insulation sites and chromosome loops compared to other genome-wide 3D technologies. A visually supported, detailed, step-by-step protocol for preparing high-quality Micro-C-XL samples from mammalian cells is presented in this article.

Differential Chromatin Organization between t(4;14) and Non-t(4;14) Multiple Myeloma Driven By the Histone Methyltransferase NSD2

Background: The incidence of multiple myeloma (MM) with t(4;14) translocation is ~15% of newly diagnosed MM (NDMM), and is associated with poor prognosis and high-risk disease. Recent analysis of t(4;14) genomics indicate that the location of the chromosome 4 (chr4) breakpoint can significantly segregate patient outcome regardless of therapeutic intervention. Translocation of IGH regulatory elements from chr14, drive the expression of NSD2 (chr4), a histone (H3) methyltransferase responsible for di-methylation of H3 lysine 36 (H3K36me2), a mark of active chromatin. Depending upon the location of the chr4 breakpoint, different isoforms of NSD2 are expressed, including full-length (NSD2-FL) and various N-terminus truncations (NSD2-TR). The mechanisms of how the location of the chr4 breakpoint and the functionalities of NSD2 isoforms may contribute to a higher risk biology in t(4;14) MM patients is not known. Here we compare and contrast the 3D chromosomal organization of t(4;14)-NSD2-FL, t(4;14)-NSD2-TR and non-t(4;14) in pre-clinical MM models. Results: Comparative analysis of chromosome conformation capture (Hi-C) indicates that there is a decrease in medium-to-long range chromosomal looping interactions accompanied by an increase in short-range interactions from t(4;14) vs. non-t(4;14) models (Figure A). Accordingly, relative quantification of chromosomal contact probabilities of intra-chromosomal (cis) contacts increase in frequency from non-t(4;14) (77%) vs. t(4;14)-NSD2-FL (87%) vs. t(4;14)-NSD2-TR (89%). Evaluation of active/inactive (A/B) 3D genomic compartments between non-t(4;14) and t(4;14) models suggest significant compartment switches (~30% genome; combined A-to-B and B-to-A). Analysis of gene neighborhoods through topologically associated domains (TADs) suggests that there are decreases in contacts at the TAD boundary regions in t(4;14)-NSD2-TR vs. t(4;14)-NSD2-FL vs. non-t(4;14) models. TAD sizes are also generally smaller in size, from t(4;14)-NSD2-TR vs. t(4;14)-NSD2-FL vs. non-t(4;14). Complimentary analysis of activity by contact (ABC) to assess the potential influence of proximal and distal enhancers on gene activity confirms there is a decrease in ABC scores when comparing t(4;14) vs. non-t(4;14). Assessment of global chromatin accessibility (ATAC-Seq) indicate significant consensus open peaks (~135k), differentially accessible (DA) peaks (~2.6k), and differential closed peaks (~8.8k) in t(4;14) vs. non-t(4;14). Principal component (PC) analysis demonstrated a unique grouping of t(4;14) models compared to a very heterogenous distribution of non-t(4;14) models. DA regions in t(4;14) were found to be primarily more enriched in both intronic and intergenic regions, and within proximal and distal promoters. Analysis of histone marks was used to assess differences in global active enhancers (H3K4me1, H3K27ac), repressive chromatin (H3K9me3, H3K27me3), and active transcription (H3K36me2/3, H3K4me3) between non-t(4;14) and t(4;14). Active marks were shown to be increased, and repressive marks decreased in t(4;14) vs. non-t(4;14) models that also correlated with DA peaks. Enhancer ranking suggests the frequency of potentially poised superenhancers increase in t(4;14) vs. non-t(4;14). Differential expression (DE) analysis between t(4;14) and non-t(4;14) was highly correlative with DA peaks (Figure B). In addition, distinctive transcriptional signatures between non-t(4;14), t(4;14)-NSD2-FL, and t(4;14)-NSD2-TR were identified. Conservative (DA+DE) pathway enrichment suggests multiple and distinct biological processes including cell adhesion through adherens junctions, cell-cycle regulation through mitotic checkpoints, chromosomal reorganization through depolymerization of the nuclear lamina, and epithelial to mesenchymal transition. Conclusions: Systematic chromosomal re-arrangement and accessibility driven by differential epigenetic deposition of histone regulatory marks, differentiates t(4;14) from non-t(4;14) MM biology. These analyses suggest that given the nature of the translocation, and additionally the location of the chr4 breakpoint, result in major 3D chromosomal structural rearrangement, mediating biological programs that may explain the differences between successful and poor interventional therapies for this subset of MM disease.

LncRNA MACC1-AS1 associates with DDX5 to modulate MACC1 transcription in breast cancer cells

MACC1 is a master oncogene involved in multiple aspects of cancer metastasis in a broad variety of tumors. However, the molecular mechanism by which MACC1 transcription is regulated remains unclear. Here, we show that in breast cancer cells, lncRNA MACC1-AS1 serves as a cis-factor to up-regulate MACC1 transcription and this regulation increases the cell proliferation potential. Mechanistically, MACC1-AS1 forms a complex with DEAD-Box helicase 5 (DDX5) and simultaneously interacts with the distal region of the MACC1 promoter. The interaction allows its associated DDX5 to spatially contact the MACC1 core promoter and shift from MACC1-AS1 to the core promoter. Moreover, binding of DDX5 to the core promoter results in local recruitment of the transcription factor SP-1, thus enhancing MACC1 transcription. Our findings reveal a molecular mechanism by which MACC1-AS1 cis-regulates MACC1 transcription by interacting with the distal promoter region and delivering DDX5 to the core-promoter of the gene.

Transcriptional regulation of the mouse EphA4, Ephrin-B2 and Ephrin-A3 genes by the circadian clock machinery

ABSTRACT Circadian rhythms originate from molecular feedback loops. In mammals, the transcription factors CLOCK and BMAL1 act on regulatory elements (i.e. E-boxes) to shape biological functions in a rhythmic manner. The EPHA4 receptor and its ligands Ephrins (EFN) are cell adhesion molecules regulating neurotransmission and neuronal morphology. Previous studies showed the presence of E-boxes in the genes of EphA4 and specific Ephrins, and that EphA4 knockout mice have an altered circadian rhythm of locomotor activity. We thus hypothesized that the core clock machinery regulates the gene expression of EphA4, EfnB2 and EfnA3. CLOCK and BMAL1 (or NPAS2 and BMAL2) were found to have transcriptional activity on distal and proximal regions of EphA4, EfnB2 and EfnA3 putative promoters. A constitutively active form of glycogen synthase kinase 3β (GSK3β; a negative regulator of CLOCK and BMAL1) blocked the transcriptional induction. Mutating the E-boxes of EphA4 distal promoter sequence reduced transcriptional induction. EPHA4 and EFNB2 protein levels did not show circadian variations in the mouse suprachiasmatic nucleus or prefrontal cortex. The findings uncover that core circadian transcription factors can regulate the gene expression of elements of the Eph/Ephrin system, which might contribute to circadian rhythmicity in biological processes in the brain or peripheral tissues.

Evaluating the potential of (epi)genotype-by-low pass nanopore sequencing in dairy cattle: a study on direct genomic value and methylation analysis

BackgroundGenotype-by-sequencing has been proposed as an alternative to SNP genotyping arrays in genomic selection to obtain a high density of markers along the genome. It requires a low sequencing depth to be cost effective, which may increase the error at the genotype assigment. Third generation nanopore sequencing technology offers low cost sequencing and the possibility to detect genome methylation, which provides added value to genotype-by-sequencing. The aim of this study was to evaluate the performance of genotype-by-low pass nanopore sequencing for estimating the direct genomic value in dairy cattle, and the possibility to obtain methylation marks simultaneously.ResultsLatest nanopore chemistry (LSK14 and Q20) achieved a modal base calling accuracy of 99.55%, whereas previous kit (LSK109) achieved slightly lower accuracy (99.1%). The direct genomic value accuracy from genotype-by-low pass sequencing ranged between 0.79 and 0.99, depending on the trait (milk, fat or protein yield), with a sequencing depth as low as 2 × and using the latest chemistry (LSK114). Lower sequencing depth led to biased estimates, yet with high rank correlations. The LSK109 and Q20 achieved lower accuracies (0.57–0.93). More than one million high reliable methylated sites were obtained, even at low sequencing depth, located mainly in distal intergenic (87%) and promoter (5%) regions.ConclusionsThis study showed that the latest nanopore technology in useful in a LowPass sequencing framework to estimate direct genomic values with high reliability. It may provide advantages in populations with no available SNP chip, or when a large density of markers with a wide range of allele frequencies is needed. In addition, low pass sequencing provided nucleotide methylation status of > 1 million nucleotides at ≥ 10 × , which is an added value for epigenetic studies.

Elevated sHLA-G plasma levels post chemotherapy combined with ILT-2 rs10416697C allele status of the sHLA-G-related receptor predict poorest disease outcome in early triple-negative breast cancer patients.

Triple negative breast cancer (TNBC) shows an aggressive growing and spreading behavior and has limited treatment options, often leading to inferior disease outcome. Therefore, surrogate markers are urgently needed to identify patients at high risk of recurrence and more importantly, to identify additional therapeutic targets enabling further treatment options. Based on the key role of the non-classical human leukocyte antigen G (HLA-G) and its related receptor immunoglobulin-like transcript receptor-2 (ILT-2) in immune evasion mechanisms of tumors, members of this ligand-receptor axis appear to be promising tool for both, defining risk groups and potential therapeutic targets. To follow this, sHLA-G levels before and after chemotherapy (CT), HLA-G 3' UTR haplotypes, and allele variations rs10416697 at the distal gene promoter region of ILT-2 were defined in healthy female controls and early TNBC patients. The results obtained were associated with clinical status, presence of circulating tumor cell (CTC) subtypes, and disease outcome of patients in terms of progression-free or overall survival. sHLA-G plasma levels were increased in TNBC patients post-CT compared to levels of patients pre-CT or controls. High post-CT sHLA-G levels were associated with the development of distant metastases, the presence of ERCC1 or PIK3CA-CTC subtypes post-CT, and poorer disease outcome in uni- or multivariate analysis. HLA-G 3' UTR genotypes did not influence disease outcome but ILT-2 rs10416697C allele was associated with AURKA-positive CTC and with adverse disease outcome by uni- and multivariate analysis. The prognostic value of the combined risk factors (high sHLA-G levels post-CT and ILT-2 rs10416697C allele carrier status) was an even better independent indicator for disease outcome in TNBC than the lymph nodal status pre-CT. This combination allowed the identification of patients with high risk of early progression/death with positive nodal status pre-CT or with non-pathological complete therapy response. The results of this study highlight for the first time that the combination of high levels of sHLA-G post-CT with ILT-2 rs10416697C allele receptor status is a promising tool for the risk assessment of TNBC patients and support the concept to use HLA-G/ILT-2 ligand-receptor axis as therapeutic targets.

Protein tyrosine phosphatase nonreceptor type 2 exerts antimetastatic functions in pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is a highly invasive tumor with a dismal prognosis. Recent studies have demonstrated PTPN2 (protein tyrosine phosphatase nonreceptor type 2) as a potential target for cancer therapy. However, the functions of PTPN2 in PDAC progression remain poorly understood. In this study, we found PTPN2 expression was downregulated in PDAC tissues, and decreased PTPN2 expression was associated with unfavorable prognosis. Functional studies indicated that PTPN2 knockdown promoted the migration and invasion abilities of PDAC cells in vitro, and the liver metastasis in vivo through epithelial-mesenchymal transition process. Mechanistically, MMP-1 was identified as a downstream target of PTPN2 via RNA-seq data and was responsible for the enhanced metastasis of PDAC cells upon PTPN2 knockdown. Moreover, according to chromatin immunoprecipitation and electrophoretic mobility shift assay, PTPN2 depletion transcriptionally activated MMP-1 via regulating the interaction of p-STAT3 with its distal promoter. This study, for the first time, demonstrated that PTPN2 inhibited PDAC metastasis, and presented a novel PTPN2/p-STAT3/MMP-1 axis in PDAC progression.

Prevalent use and evolution of exonic regulatory sequences in the human genome

AbstractIt has long been known that exons can serve as cis‐regulatory sequences, such as enhancers. However, the prevalence of such dual‐use of exons and how they evolve remain elusive. Based on our recently predicted, highly accurate large sets of cis‐regulatory module candidates (CRMCs) and non‐CRMCs in the human genome, we find that exonic transcription factor binding sites (TFBSs) occupy at least a third of the total exon lengths, and 96.7% of genes have exonic TFBSs. Both A/T and C/G in exonic TFBSs are more likely under evolutionary constraints than those in non‐CRMC exons. Exonic TFBSs in codons tend to encode loops rather than more critical helices and strands in protein structures, while exonic TFBSs in untranslated regions (UTRs) tend to avoid positions where known UTR‐related functions are located. Moreover, active exonic TFBSs tend to be in close physical proximity to distal promoters whose genes have elevated transcription levels. These results suggest that exonic TFBSs might be more prevalent than originally thought and likely in dual‐use. We proposed a parsimonious model that well explains the observed evolutionary behaviors of exonic TFBS as well as how a stretch of codons evolve into a TFBS.Key points There are more exonic regulatory sequences in the human genome than originally thought. Exonic transcription factor binding sites are more likely under negative selection or positive selection than counterpart nonregulatory sequences. Exonic transcription factor binding sites tend to be located in genome sequences that encode less critical loops in protein structures, or in less critical parts in 5′ and 3′ untranslated regions.

NFATc1 induction by an intronic enhancer restricts NKT γδ cell formation

In thymus, the ablation of Tcell receptor (TCR)-activated transcription factor NFATc1 or its inducible isoforms during the double-negative (DN) stages of thymocyte development leads to a marked increase in γδ thymocytes whereas the development of αβ thymocytes remains mostly unaffected. These γδ thymocytes are characterized by the upregulation of the promyelocytic leukemia zinc-fingerfactor (PLZF), the "master regulator" of natural killer T (NKT) cell development, and the acquisition of an NKT γδ cell phenotype with higher cell survival rates. The suppressive function of NFATc1 in NKT γδ cell formation critically depends on the remote enhancer E2, which is essential for the inducible expression of NFATc1 directed by its distal promoter P1. Thus, the enhancer deciphers a strong γδ TCR signal into the expression of inducible NFATc1 isoforms resulting in high levels of NFATc1 protein that are essential to control the numbers of NKT γδ cells.