Circular Chromosome Conformation Capture Research Articles

A distal enhancer of GATA3 regulates Th2 differentiation and allergic inflammation

Asthma is a widespread airway disorder where GATA3-dependent Type-2 helper T (Th2) cells and group 2 innate lymphoid cells (ILC2s) play vital roles. Asthma-associated single nucleotide polymorphisms (SNPs) are enriched in a region located 926-970 kb downstream from GATA3 in the 10p14 (hG900). However, it is unknown how hG900 affects the pathogenesis of allergic airway inflammation. To investigate the roles of the asthma-associated GATA3 enhancer region in experimental allergic airway inflammation, we first examined the correlation between GATA3 expression and the activation of the hG900 region was analyzed by flow cytometry and ChIP-qPCR. We found that The activation of enhancers in the hG900 region was strongly correlated to the levels of GATA3 in human peripheral T cell subsets. We next generated mice lacking the mG900 region (mG900KO mice) were generated by the CRISPR-Cas9 system, and the development and function of helper T cells and ILCs in mG900KO mice were analyzed in steady-state conditions and allergic airway inflammation induced by papain or house dust mite (HDM). The deletion of the mG900 did not affect the development of lymphocytes in steady-state conditions or allergic airway inflammation induced by papain. However, mG900KO mice exhibited reduced allergic inflammation and Th2 differentiation in the HDM-induced allergic airway inflammation. The analysis of the chromatin conformation around Gata3 by circular chromosome conformation capture coupled to high-throughput sequencing (4C-seq) revealed that the mG900 region interacted with the transcription start site of Gata3 with an influencing chromatin conformation in Th2 cells. These findings indicate that the mG900 region plays a pivotal role in Th2 differentiation and thus enhances allergic airway inflammation.

A novel cis-regulatory element regulates αD and αA-globin gene expression in chicken erythroid cells.

Cis-regulatory elements (CREs) play crucial roles in regulating gene expression during erythroid cell differentiation. Genome-wide erythroid-specific CREs have not been characterized in chicken erythroid cells, which is an organism model used to study epigenetic regulation during erythropoiesis. Analysis of public genome-wide accessibility (ATAC-seq) maps, along with transcription factor (TF) motif analysis, CTCF, and RNA Pol II occupancy, as well as transcriptome analysis in fibroblasts and erythroid HD3 cells, were used to characterize erythroid-specific CREs. An α-globin CRE was identified, and its regulatory activity was validated in vitro and in vivo by luciferase activity and genome-editing assays in HD3 cells, respectively. Additionally, circular chromosome conformation capture (UMI-4C) assays were used to distinguish its role in structuring the α-globin domain in erythroid chicken cells. Erythroid-specific CREs displayed occupancy by erythroid TF binding motifs, CTCF, and RNA Pol II, as well as an association with genes involved in hematopoiesis and cell differentiation. An α-globin CRE, referred to as CRE-2, was identified as exhibiting enhancer activity over αD and αA genes in vitro and in vivo. Induction of terminal erythroid differentiation showed that α-globin CRE-2 is required for the induction of αD and αA. Analysis of TF binding motifs at α-globin CRE-2 shows apparent regulation mediated by GATA-1, YY1, and CTCF binding. Our findings demonstrate that cell-specific CREs constitute a key mechanism that contributes to the fine-tuning gene regulation of erythroid cell differentiation and provide insights into the annotation and characterization of CREs in chicken cells.

Three-Dimensional Chromosomal Landscape Revealing miR-146a Dysfunctional Enhancer in Lupus and Establishing a CRISPR-Mediated Approach to Inhibit the Interferon Pathway.

The diminished expression of microRNA-146a (miR-146a) in systemic lupus erythematosus (SLE) contributes to the aberrant activation of the interferon pathway. Despite its significance, the underlying mechanism driving this reduced expression remains elusive. Considering the integral role of enhancers in steering gene expression, our study seeks to pinpoint the SLE-affected enhancers responsible for modulating miR-146a expression. Additionally, we aim to elucidate the mechanisms by which these enhancers influence the contribution of miR-146a to the activation of the interferon pathway. Circular chromosome conformation capture sequencing and epigenomic profiles were applied to identify candidate enhancers of miR-146a. CRISPR activation was performed to screen functional enhancers. Differential analysis of chromatin accessibility was used to identify SLE-dysregulated enhancers, and the mechanism underlying enhancer dysfunction was investigated by analyzing transcription factor binding. The therapeutic value of a lupus-related enhancer was further evaluated by targeting it in the peripheral blood mononuclear cells (PBMCs) of patients with SLE through a CRISPR activation approach. We identified shared and cell-specific enhancers of miR-146a in distinct immune cells. An enhancer 32.5 kb downstream of miR-146a possesses less accessibility in SLE, and its chromatin openness was negatively correlated with SLE disease activity. Moreover, CCAAT/enhancer binding protein α, a down-regulated transcription factor in patients with SLE, binds to the 32.5-kb enhancer and induces the epigenomic change of this locus. Furthermore, CRISPR-based activation of this enhancer in SLE PBMCs could inhibit the activity of interferon pathway. Our work defines a promising target for SLE intervention. We adopted integrative approaches to define cell-specific and functional enhancers of the SLE critical gene and investigated the mechanism underlying its dysregulation mediated by a lupus-related enhancer.

Identification of enhancers responsible for the coordinated expression of myosin heavy chain isoforms in skeletal muscle

BackgroundSkeletal muscles consist of fibers of differing contractility and metabolic properties, which are primarily determined by the content of myosin heavy chain (MYH) isoforms (MYH7, MYH2, MYH1, and MYH4). The regulation of Myh genes transcription depends on three-dimensional chromatin conformation interaction, but the mechanistic details remain to be determined.ResultsIn this study, we characterized the interaction profiles of Myh genes using 4C-seq (circular chromosome conformation capture coupled to high-throughput sequencing). The interaction profile of Myh genes changed between fast quadriceps and slow soleus muscles. Combining chromatin immunoprecipitation-sequencing (ChIP-seq) and transposase accessible chromatin with high-throughput sequencing (ATAC-seq), we found that a 38 kb intergenic region interacting simultaneously with fast Myh genes promoters controlled the coordinated expression of fast Myh genes. We also identified four active enhancers of Myh7, and revealed that binding of MYOG and MYOD increased the activity of Myh7 enhancers.ConclusionsThis study provides new insight into the chromatin interactions that regulate Myh genes expression.

MYC overexpression leads to increased chromatin interactions at super-enhancers and MYC binding sites.

The MYC oncogene encodes for the MYC protein and is frequently dysregulated across multiple cancer cell types, making it an attractive target for cancer therapy. MYC overexpression leads to MYC binding at active enhancers, resulting in a global transcriptional amplification of active genes. Because super-enhancers are frequently dysregulated in cancer, we hypothesized that MYC preferentially invades into super-enhancers and alters the cancer genome organization. To that end, we performed ChIP-seq, RNA-seq, circular chromosome conformation capture (4C-seq), and Spike-in Quantitative Hi-C (SIQHiC) on the U2OS osteosarcoma cell line with tetracycline-inducible MYC. MYC overexpression in U2OS cells modulated histone acetylation and increased MYC binding at super-enhancers. SIQHiC analysis revealed increased global chromatin contact frequency, particularly at chromatin interactions connecting MYC binding sites at promoters and enhancers. Immunofluorescence staining showed that MYC molecules formed punctate foci at these transcriptionally active domains after MYC overexpression. These results demonstrate the accumulation of overexpressed MYC at promoter–enhancer hubs and suggest that MYC invades into enhancers through spatial proximity. At the same time, the increased protein–protein interactions may strengthen these chromatin interactions to increase chromatin contact frequency. CTCF siRNA knockdown in MYC-overexpressed U2OS cells demonstrated that removal of architectural proteins can disperse MYC and abrogate the increase in chromatin contacts. By elucidating the chromatin landscape of MYC-driven cancers, we can potentially target MYC-associated chromatin interactions for cancer therapy.

Circular Chromosome Conformation Capture Sequencing (4C-Seq ) in Primary Adherent Cells.

The three-dimensional structure of the genome is highly organized and is an important aspect of gene regulation. Chromatin interactions can be identified using chromosome conformation capture-based techniques, which rely on proximity ligation. Of these techniques, circular chromosome conformation capture sequencing (4C-seq) is used to identify all chromatin interactions occurring with a single chromosomal location (one versus all). Here we describe a 4C-seq protocol that has been optimized for primary adherent cells, for which the first digestion step is inefficient using standard 4C-seq protocols. It can, however, also be applied to other cell or tissue types. This protocol utilizes a standard DNA library preparation method using a commercial kit, and includes a description of the data processing steps.

Detection of Allele-Specific 3D Chromatin Interactions Using High-Resolution In-Nucleus 4C-seq.

Chromosome conformation capture techniques are a set of methods used to determine 3D genome organization through the capture and identification of physical contacts between pairs of genomic loci. Among them, 4C-seq (circular chromosome conformation capture coupled to high-throughput sequencing) allows for the identification and quantification of the sequences interacting with a preselected locus of interest. 4C-seq has been widely used in the literature, mainly to study chromatin loops between enhancers and promoters or between CTCF binding sites and to identify chromatin domain boundaries. As 3D-contacts may be established in an allele-specific manner, we describe an up-to-date allele-specific 4C-seq protocol, starting from the selection of allele-specific viewpoints to Illumina sequencing. This protocol has mainly been optimized for cultured mammalian cells, but can be adapted for other cell types with relatively minor changes in initial steps.

Chromosomal Translocations in NK-Cell Lymphomas Originate from Inter-Chromosomal Contacts of Active rDNA Clusters Possessing Hot Spots of DSBs.

Simple SummaryThere are nine DSB hot spots located in the non-transcribed spacer of human rDNA units. Circular chromosome conformation capture data indicate that the rDNA clusters often shape contact with a specific set of chromosomal regions containing genes controlling differentiation and cancer, and often possessing the DSB hot spots. The data suggest a mechanism for rDNA-mediated translocation, and some of them could lead to tumorigenesis. Here, we searched for translocations in which rDNA clusters are involved. WGS data of normal T cells and NK-cell lymphomas from the same individuals were used. We revealed numerous translocations in which rDNA units are involved. The sites of these translocations in normal T cells and in the lymphomas were mostly different, but occurred at about the same frequency in both cell types. We conclude that oncogenic translocations lead to dysregulation of a specific set of genes controlling development.Endogenous hot spots of DNA double-strand breaks (DSBs) are tightly linked with transcription patterns and cancer. There are nine hot spots of DSBs (denoted Pleiades) in human rDNA units that are located exclusively inside the intergenic spacer (IGS). Profiles of Pleiades coincide with the profiles of γ-H2AX, suggesting a high level of in vivo breakage inside rDNA genes. The data were confirmed by microscopic observation of the largest γ-H2AX foci inside nucleoli in interphase chromosomes. Circular chromosome conformation capture (4C) data indicate that the rDNA units often make contact with a specific set of chromosomal regions containing genes that are involved in differentiation and cancer. Interestingly, these regions also often possess hot spots of DSBs that provide the potential for Robertsonian and oncogenic translocations. In this study, we searched for translocations in which rDNA clusters are involved. The whole genome sequence (WGS) data of normal T cells and NK-cell lymphomas from the same individuals revealed numerous translocations in which Pleiades were involved. The sites of these translocations in normal T cells and in the lymphomas were mostly different, although there were also some common sites. The genes at translocations in normal cells and in lymphomas are associated with predominantly non-overlapping lists of genes that are depleted with silenced genes. Our data indicate that rDNA-mediated translocations occur at about the same frequency in the normal T cells and NK-lymphoma cells but differ at particular sites that correspond to open chromatin. We conclude that oncogenic translocations lead to dysregulation of a specific set of genes controlling development. In normal T cells and in NK cells, there are hot spots of translocations at sites possessing strong H3K27ac marks. The data indicate that Pleiades are involved in rDNA-mediated translocation.

Genetics and Genomics of SOST: Functional Analysis of Variants and Genomic Regulation in Osteoblasts.

SOST encodes the sclerostin protein, which acts as a key extracellular inhibitor of the canonical Wnt pathway in bone, playing a crucial role in skeletal development and bone homeostasis. The objective of this work was to assess the functionality of two variants previously identified (the rare variant rs570754792 and the missense variant p.Val10Ile) and to investigate the physical interactors of the SOST proximal promoter region in bone cells. Through a promoter luciferase reporter assay we show that the minor allele of rs570754792, a variant located in the extended TATA box motif, displays a significant decrease in promoter activity. Likewise, through western blot studies of extracellular and intracellular sclerostin, we observe a reduced expression of the p.Val10Ile mutant protein. Finally, using a circular chromosome conformation capture assay (4C-seq) in 3 bone cell types (MSC, hFOB, Saos-2), we have detected physical interactions between the SOST proximal promoter and the ECR5 enhancer, several additional enhancers located between EVT4 and MEOX1 and a distant region containing exon 18 of DHX8. In conclusion, SOST presents functional regulatory and missense variants that affect its expression and displays physical contacts with far reaching genomic sequences, which may play a role in its regulation within bone cells.

Regulatory CDH4 Genetic Variants Associate With Risk to Develop Capecitabine-Induced Hand-Foot Syndrome.

Capecitabine-induced hand-foot syndrome (CiHFS) is a common dermatological adverse reaction affecting around 30% of patients with capecitabine-treated cancer, and the main cause of dose reductions and chemotherapy delays. To identify novel genetic factors associated with CiHFS in patients with cancer, we carried out an extreme-phenotype genomewide association study in 166 patients with breast and colorectal capecitabine-treated cancer with replication in a second cohort of 85 patients. We discovered and replicated a cluster of four highly correlated single-nucleotide polymorphisms associated with susceptibility to CiHFS at 20q13.33 locus (top hit=rs6129058, hazard ratio=2.40, 95% confidence interval=1.78-3.20; P=1.2×10-8 ). Using circular chromosome conformation capture sequencing, we identified a chromatin contact between the locus containing the risk alleles and the promoter of CDH4, located 90 kilobases away. The risk haplotype was associated with decreased levels of CDH4 mRNA and the protein it encodes, R-cadherin (RCAD), which mainly localizes in the granular layer of the epidermis. In human keratinocytes, CDH4 downregulation resulted in reduced expression of involucrin, a protein of the cornified envelope, an essential structure for skin barrier function. Immunohistochemical analyses revealed that skin from patients with severe CiHFS exhibited low levels of RCAD and involucrin before capecitabine treatment. Our results uncover a novel mechanism underlying individual genetic susceptibility to CiHFS with implications for clinically relevant risk prediction.

4C-Seq: Interrogating Chromatin Looping with Circular Chromosome Conformation Capture.

Chromosome conformation capture and its variants have allowed chromatin topology to be interrogated at a superior resolution and throughput than by microscopic methods. Among the method derivatives, 4C-seq (circular chromosome conformation capture, coupled to high-throughput sequencing) is a versatile, cost-effective means of assessing all chromatin interactions with a specific genomic region of interest, making it particularly suitable for interrogating chromatin looping events. We present the principles and procedures for designing and implementing successful 4C-seq experiments.

TP63, SOX2, and KLF5 Establish a Core Regulatory Circuitry That Controls Epigenetic and Transcription Patterns in Esophageal Squamous Cell Carcinoma Cell Lines

We investigated the transcriptome of esophageal squamous cell carcinoma (ESCC) cells, activity of gene regulatory (enhancer and promoter regions), and the effects of blocking epigenetic regulatory proteins. We performed chromatin immunoprecipitation sequencing with antibodies against H3K4me1, H3K4me3, and H3K27ac and an assay for transposase-accessible chromatin to map the enhancer regions and accessible chromatin in 8 ESCC cell lines. We used the CRC_Mapper algorithm to identify core regulatory circuitry transcription factors in ESCC cell lines, and determined genome occupancy profiles for 3 of these factors. In ESCC cell lines, expression of transcription factors was knocked down with small hairpin RNAs, promoter and enhancer regions were disrupted by CRISPR/Cas9 genome editing, or bromodomains and extraterminal (BET) family proteins and histone deacetylases (HDACs) were inhibited with ARV-771 and romidepsin, respectively. ESCC cell lines were then analyzed by whole-transcriptome sequencing, immunoprecipitation, immunoblots, immunohistochemistry, and viability assays. Interactions between distal enhancers and promoters were identified and verified with circular chromosome conformation capture sequencing. NOD-SCID mice were given injections of modified ESCC cells, some mice where given injections of HDAC or BET inhibitors, and growth of xenograft tumors was measured. We identified super-enhancer-regulated circuits and transcription factors TP63, SOX2, and KLF5 as core regulatory factors in ESCC cells. Super-enhancer regulation of ALDH3A1 mediated by core regulatory factors was required for ESCC viability. We observed direct interactions between the promoter region of TP63 and functional enhancers, mediated by the core regulatory circuitry transcription factors. Deletion of enhancer regions from ESCC cells decreased expression of the core regulatory circuitry transcription factors and reduced cell viability; these same results were observed with knockdown of each core regulatory circuitry transcription factor. Incubation of ESCC cells with BET and HDAC disrupted the core regulatory circuitry program and the epigenetic modifications observed in these cells; mice given injections of HDAC or BET inhibitors developed smaller xenograft tumors from the ESCC cell lines. Xenograft tumors grew more slowly in mice given the combination of ARV-771 and romidepsin than mice given either agent alone. In epigenetic and transcriptional analyses of ESCC cell lines, we found the transcription factors TP63, SOX2, and KLF5 to be part of a core regulatory network that determines chromatin accessibility, epigenetic modifications, and gene expression patterns in these cells. A combination of epigenetic inhibitors slowed growth of xenograft tumors derived from ESCC cells in mice.

The Position and Complex Genomic Architecture of Plant T-DNA Insertions Revealed by 4SEE.

The integration of T-DNA in plant genomes is widely used for basic research and agriculture. The high heterogeneity in the number of integration events per genome, their configuration, and their impact on genome integrity highlight the critical need to detect the genomic locations of T-DNA insertions and their associated chromosomal rearrangements, and the great challenge in doing so. Here, we present 4SEE, a circular chromosome conformation capture (4C)-based method for robust, rapid, and cost-efficient detection of the entire scope of T-DNA locations. Moreover, by measuring the chromosomal architecture of the plant genome flanking the T-DNA insertions, 4SEE outlines their associated complex chromosomal aberrations. Applying 4SEE to a collection of confirmed T-DNA lines revealed previously unmapped T-DNA insertions and chromosomal rearrangements such as inversions and translocations. Uncovering such events in a feasible, robust, and cost-effective manner by 4SEE in any plant of interest has implications for accurate annotation and phenotypic characterization of T-DNA insertion mutants and transgene expression in basic science applications as well as for plant biotechnology.

Epigenetic specifications of host chromosome docking sites for latent Epstein-Barr virus

Epstein-Barr virus (EBV) genomes persist in latently infected cells as extrachromosomal episomes that attach to host chromosomes through the tethering functions of EBNA1, a viral encoded sequence-specific DNA binding protein. Here we employ circular chromosome conformation capture (4C) analysis to identify genome-wide associations between EBV episomes and host chromosomes. We find that EBV episomes in Burkitt’s lymphoma cells preferentially associate with cellular genomic sites containing EBNA1 binding sites enriched with B-cell factors EBF1 and RBP-jK, the repressive histone mark H3K9me3, and AT-rich flanking sequence. These attachment sites correspond to transcriptionally silenced genes with GO enrichment for neuronal function and protein kinase A pathways. Depletion of EBNA1 leads to a transcriptional de-repression of silenced genes and reduction in H3K9me3. EBV attachment sites in lymphoblastoid cells with different latency type show different correlations, suggesting that host chromosome attachment sites are functionally linked to latency type gene expression programs.

4See: A Flexible Browser to Explore 4C Data.

It is established that transcription of many metazoan genes is regulated by distal regulatory sequences beyond the promoter. Enhancers have been identified at up to megabase distances from their regulated genes, and/or proximal to or within the introns of unregulated genes. The unambiguous identification of the target genes of newly identified regulatory elements can thus be challenging. Well-studied enhancers have been found to come into direct physical proximity with regulated genes, presumably by the formation of chromatin loops. Chromosome conformation capture (3C) derivatives that assess the frequency of proximity between different genetic elements is thus a popular method for exploring gene regulation by distal regulatory elements. For studies of chromatin loops and promoter-enhancer communication, 4C (circular chromosome conformation capture) is one of the methods of choice, optimizing cost (required sequencing depth), throughput, and resolution. For ease of visual inspection of 4C data we present 4See, a versatile and user-friendly browser. 4See allows 4C profiles from the same bait to be flexibly plotted together, allowing biological replicates to either be compared, or pooled for comparisons between different cell types or experimental conditions. 4C profiles can be integrated with gene tracks, linear epigenomic profiles, and annotated regions of interest, such as called significant interactions, allowing rapid data exploration with limited computational resources or bioinformatics expertise.

Drosophila rDNA Genes Shape the Stable Contacts with the Tlk Gene at the Expression Area of Small RNAs and Affect on Looped Domains inside the Gene

In experiments on mouse and human cells it was demonstrated that rDNA plays an important role in epigenetic regulation of many genes. To identify and study rDNA-contacting genes in Drosophila we used the 4С (circular chromosome conformation capture) approach. We detected very stable contacts of rDNA genes within a 5-kb region inside the Tlk gene residing in X chromosome. This 5-kb region corresponds to small RNAs. After heat shock treatment both the amount of contacts, and the expression level of the gene were increased. Tlk and Rala are genes that share the same short bidirectional promoter but exhibit different expression levels. Around the region of rDNA contacts inside the Tlk gene, looped domains were formed. We conclude that rDNA contact-dependent epigenetic regulation is guided by small RNAs and that the contacts are involved in rearrangements of the looped domains.

Contact Sites of rDNA Clusters with FANK1 Gene Correspond to Repressed Chromatin

rDNA genes play an important role in epigenetic regulation and in differentiation of eukaryotic cells. Using the 4C (circular chromosome conformation capture) approach and model HEK293T cells, we analyzed the rDNA-contacting gene FANK1, using anchor located inside rDNA genes. At the 5' end of the FANK1 gene we detected frequent contacts with rDNA clusters. The contact sites coincide with the border where chromatin state changes and nucleosome positioning. The adjacent genes DHX32, BCCIP and UROS are located in the active chromatin and are transcribed, but do not contact with rDNA genes, while FANK1 gene is silenced, and is located in repressed chromatin. Heat shock treatment dramatically changes the pattern of rDNA contacts in the region and induces about 4-fold increase in activation of the FANK1 gene. We conclude that rDNA contacts may be involved in repression of the FANK1 gene.

A non-coding genetic variant associated with abdominal aortic aneurysm alters ERG gene regulation.

Abdominal aortic aneurysm (AAA) is a major cause of sudden death in the elderly. While AAA has some overlapping genetic and environmental risk factors with atherosclerosis, there are substantial differences, and AAA-specific medication is lacking. A recent meta-analysis of genome-wide association studies has identified four novel single-nucleotide polymorphisms (SNPs) specifically associated with AAA. Here, we investigated the gene regulatory function for one of four non-coding SNPs associated with AAA, rs2836411, which is located in an intron of the ERG gene. Rs2836411 resides within a >70 kb super-enhancer that has high levels of H3K27ac and H3K4me1 in vascular endothelial and haematopoietic cell types. Enhancer luciferase assays in cell lines showed that the risk allele significantly alters enhancer activity. The risk allele also correlates with reduced ERG expression in aortic and other vascular tissues. To identify whether rs2836411 directly contacts the promoters of ERG and/or of genes further away, we performed allele-specific circular chromosome conformation capture sequencing. In vascular endothelial cells, which express ERG, the SNP region interacts highly within the super-enhancer, while in vascular smooth muscle cells, which do not express ERG, the interactions are distributed across a wider region that includes neighbouring genes. Furthermore, the risk allele has fewer interactions within the super-enhancer compared to the protective allele. In conclusion, our results indicate that rs2836411 likely affects ERG expression by altering enhancer activity and changing local chromatin interactions. ERG is involved in vascular development, angiogenesis, and inflammation in atherosclerosis; therefore mechanistically, rs2836411 could contribute to AAA by modulating ERG levels.

Male-Specific Long Noncoding RNA TTTY15 Inhibits Non-Small Cell Lung Cancer Proliferation and Metastasis via TBX4.

Gender affects cancer susceptibility. Currently, there are only a few studies on Y chromosome-linked long noncoding RNAs (lncRNAs), and the potential association between lncRNAs and cancers in males has not been fully elucidated. Here, we examined the expression of testis-specific transcript Y-linked 15 (TTTY15) in 37 males with non-small cell lung cancer (NSCLC), and performed circular chromosome conformation capture with next-generation sequencing to determine the genomic interaction regions of the TTTY15 gene. Our results showed that the expression levels of TTTY15 were lower in NSCLC tissues. Lower TTTY15 expression levels were associated with Tumor-Node-Metastasis (TNM) stage. A TTTY15 knockdown promoted malignant transformation of NSCLC cells. Based on the bioinformatics analysis of circular chromosome conformation capture data, we found that T-box transcription factor 4 (TBX4) may be a potential target gene of TTTY15. The RNA immunoprecipitation and chromatin immunoprecipitation results showed that TTTY15 may interact with DNA (cytosine-5)-methyltransferase 3A (DNMT3A), and the TTTY15 knockdown increased the binding of DNMT3A to the TBX4 promoter. We concluded that low TTTY15 expression correlates with worse prognosis among patients with NSCLC. TTTY15 promotes TBX4 expression via DNMT3A-mediated regulation. The identification of lncRNAs encoded by male-specific genes may help to identify potential targets for NSCLC therapy.

Identification of a novel enhancer of CEBPE essential for granulocytic differentiation

AbstractCCAAT/enhancer binding protein ε (CEBPE) is an essential transcription factor for granulocytic differentiation. Mutations of CEBPE occur in individuals with neutrophil-specific granule deficiency (SGD), which is characterized by defects in neutrophil maturation. Cebpe-knockout mice also exhibit defects in terminal differentiation of granulocytes, a phenotype reminiscent of SGD. Analysis of DNase I hypersensitive sites sequencing data revealed an open chromatin region 6 kb downstream of the transcriptional start site of Cebpe in murine myeloid cells. We identified an interaction between this +6-kb region and the core promoter of Cebpe using circular chromosome conformation capture sequencing (4C-seq). To understand the role of this putative enhancer in transcriptional regulation of Cebpe, we targeted it using catalytically inactive Cas9 fused to Krüppel-associated box (KRAB) domain and observed a significant downregulation of transcript and protein levels of CEBPE in cells expressing guide RNA targeting the +6-kb region. To further investigate the role of this novel enhancer further in myelopoiesis, we generated mice with deletion of this region using CRISPR/Cas9 technology. Germline deletion of the +6-kb enhancer resulted in reduced levels of CEBPE and its target genes and caused a severe block in granulocytic differentiation. We also identified binding of CEBPA and CEBPE to the +6-kb enhancer, which suggests their role in regulating the expression of Cebpe. In summary, we have identified a novel enhancer crucial for regulating expression of Cebpe and required for normal granulocytic differentiation.