Chromatin remodeling mediated by the FOXA1/A2 transcription factors activates CFTR expression in intestinal epithelial cells

The impact of genomic distance on enhancer-promoter interactions at the CFTR locus

We identified and characterized multiple cell-type selective enhancers of the CFTR gene promoter in previous work and demonstrated active looping of these elements to the promoter. Here we address the impact of genomic spacing on these enhancer:promoter interactions and on CFTR gene expression. Using CRISPR/Cas9, we generated clonal cell lines with deletions between the -35 kb airway enhancer and the CFTR promoter in the 16HBE14o- airway cell line, or between the intron 1 (185 + 10 kb) intestinal enhancer and the promoter in the Caco2 intestinal cell line. The effect of these deletions on CFTR transcript abundance, as well as the 3D looping structure of the locus was investigated in triplicate clones of each modification. Our results indicate that both small and larger deletions upstream of the promoter can perturb CFTR expression and -35 kb enhancer:promoter interactions in the airway cells, though the larger deletions are more impactful. In contrast, the small intronic deletions have no effect on CFTR expression and intron 1 enhancer:promoter interactions in the intestinal cells, whereas larger deletions do. Clonal variation following a specific CFTR modification is a confounding factor particularly in 16HBE14o- cells.

Peak Scores Significantly Depend on the Relationships between Contextual Signals in ChIP-Seq Peaks

Chromatin immunoprecipitation followed by massively parallel DNA sequencing (ChIP-seq) is a central genome-wide method for in vivo analyses of DNA-protein interactions in various cellular conditions. Numerous studies have demonstrated the complex contextual organization of ChIP-seq peak sequences and the presence of binding sites for transcription factors in them. We assessed the dependence of the ChIP-seq peak score on the presence of different contextual signals in the peak sequences by analyzing these sequences from several ChIP-seq experiments using our fully enumerative GPU-based de novo motif discovery method, Argo_CUDA. Analysis revealed sets of significant IUPAC motifs corresponding to the binding sites of the target and partner transcription factors. For these ChIP-seq experiments, multiple regression models were constructed, demonstrating a significant dependence of the peak scores on the presence in the peak sequences of not only highly significant target motifs but also less significant motifs corresponding to the binding sites of the partner transcription factors. A significant correlation was shown between the presence of the target motifs FOXA2 and the partner motifs HNF4G, which found experimental confirmation in the scientific literature, demonstrating the important contribution of the partner transcription factors to the binding of the target transcription factor to DNA and, consequently, their important contribution to the peak score.

FoxA1/2-dependent epigenomic reprogramming drives lineage switching in lung adenocarcinoma

The ability of cancer cells to alter their identity is essential for tumor survival and progression. Loss of the pulmonary lineage specifier NKX2-1 within KRAS-driven lung adenocarcinoma (LUAD) enhances tumor progression and results in a pulmonary-to-gastric lineage switch that is dependent upon the activity of pioneer factors FoxA1 and FoxA2; however, the underlying mechanism remains largely unknown. Here, we show that FoxA1/2 reprogram the epigenetic landscape of NKX2-1-negative LUAD to facilitate a gastric identity. After Nkx2-1 deletion, FoxA1/2 mediate demethylation of gastric-defining genes through recruitment of TET3, an enzyme that induces DNA demethylation. H3K27ac ChIP-seq and HiChIP show that FoxA1/2 also control the activity of regulatory elements and their 3D interactions at gastric loci. Furthermore, oncogenic KRAS is required for the FoxA1/2-dependent epigenetic reprogramming. This work demonstrates the role of FoxA1/2 in rewiring the methylation and histone landscape and cis-regulatory dynamics of NKX2-1-negative LUAD to drive cancer cell lineage switching.

The chromatin network helps prevent cancer-associated mutagenesis at transcription-replication conflicts

Genome instability is a feature of cancer cells, transcription being an important source of DNA damage. This is in large part associated with R-loops, which hamper replication, especially at head-on transcription-replication conflicts (TRCs). Here we show that TRCs trigger a DNA Damage Response (DDR) involving the chromatin network to prevent genome instability. Depletion of the key chromatin factors INO80, SMARCA5 and MTA2 results in TRCs, fork stalling and R-loop-mediated DNA damage which mostly accumulates at S/G2, while histone H3 Ser10 phosphorylation, a mark of chromatin compaction, is enriched at TRCs. Strikingly, TRC regions show increased mutagenesis in cancer cells with signatures of homologous recombination deficiency, transcription-coupled nucleotide excision repair (TC-NER) and of the AID/APOBEC cytidine deaminases, being predominant at head-on collisions. Thus, our results support that the chromatin network prevents R-loops and TRCs from genomic instability and mutagenic signatures frequently associated with cancer.

Tissue-Specific Regulation of CFTR Gene Expression

More than 2000 variations are described within the CFTR (Cystic Fibrosis Transmembrane Regulator) gene and related to large clinical issues from cystic fibrosis to mono-organ diseases. Although these CFTR-associated diseases have been well documented, a large phenotype spectrum is observed and correlations between phenotypes and genotypes are still not well established. To address this issue, we present several regulatory elements that can modulate CFTR gene expression in a tissue-specific manner. Among them, cis-regulatory elements act through chromatin loopings and take part in three-dimensional structured organization. With tissue-specific transcription factors, they form chromatin modules and can regulate gene expression. Alterations of specific regulations can impact and modulate disease expressions. Understanding all those mechanisms highlights the need to expand research outside the gene to enhance our knowledge.

Cross-talk between enhancers, structural elements and activating transcription factors maintains the 3D architecture and expression of the CFTR gene

Robust protocols to examine 3D chromatin structure have greatly advanced knowledge of gene regulatory mechanisms. Here we focus on the cystic fibrosis transmembrane conductance regulator (CFTR) gene, which provides a paradigm for validating models of gene regulation built upon genome-wide analysis. We examine the mechanisms by which multiple cis-regulatory elements (CREs) at the CFTR gene coordinate its expression in intestinal epithelial cells. Using CRISPR/Cas9 to remove CREs, individually and in tandem, followed by assays of gene expression and higher-order chromatin structure (4C-seq), we reveal the cross-talk and dependency of two cell-specific intronic enhancers. The results suggest a mechanism whereby the locus responds when CREs are lost, which may involve activating transcription factors such as FOXA2. Also, by removing the 5' topologically-associating domain (TAD) boundary, we illustrate its impact on CFTR gene expression and architecture. These data suggest a multi-layered regulatory hierarchy that is highly sensitive to perturbations.

Krüppel-Like Factor 5 Regulates CFTR Expression Through Repression by Maintaining Chromatin Architecture Coupled with Direct Enhancer Activation

Single cell RNA-sequencing has accurately identified cell types within the human airway that express the Cystic Fibrosis Transmembrane Conductance regulator (CFTR) gene. Low abundance CFTR transcripts are seen in many secretory cells, while high levels are restricted to rare pulmonary ionocytes. Here we focus on the mechanisms coordinating basal CFTR expression in the secretory compartment. Cell-selective regulation of CFTR is achieved within its invariant topologically associating domain by the recruitment of cis-regulatory elements (CREs). CRE activity is coordinated by cell-type-selective transcription factors. One such factor, Krüppel-Like Factor 5 (KLF5), profoundly represses CFTR transcript and protein in primary human airway epithelial cells and airway cell lines. Here we reveal the mechanism of action of KLF5 upon the CFTR gene. We find that depletion or ablation of KLF5 from airway epithelial cells changes higher order chromatin structure at the CFTR locus. Critical looping interactions that are required for normal gene expression are altered, the H3K27ac active chromatin mark is redistributed, and CTCF occupancy is modified. However, mutation of a single KLF5 binding site within a pivotal airway cell CRE abolishes CFTR expression. Hence, KLF5 has both direct activating and indirect repressive effects, which together coordinate CFTR expression in the airway.

Epigenome editing of the CFTR-locus for treatment of cystic fibrosis

BACKGROUND

Mechanisms governing the diversity of CFTR gene expression throughout the body are complex. Multiple intronic and distal regulatory elements are responsible for regulating differential CFTR expression across tissues.

METHODS

Drawing on published data, 18 high-priority genomic regions were identified and interrogated for CFTR-enhancer function using CRISPR/dCas9-based epigenome editing tools. Each region was evaluated by dCas9p300 and dCas9KRAB for its ability to enhance or repress CFTR expression, respectively.

RESULTS

Multiple genomic regions were tested for enhancer activity using CRISPR/dCas9 epigenome editing. dCas9p300 mediates a significant increase in CFTR mRNA levels when targeted to the promoter and a region 44 kb upstream of the transcriptional start site in a CFTR-low expressing cell line. Multiple gRNAs targeting the promoter induced a robust increase in CFTR protein levels. In contrast, dCas9KRAB-mediated repression is much more robust with 10 of the 18 evaluated genomic regions inducing CFTR protein knockdown. To evaluate the therapeutic efficacy of modulating CFTR gene regulation, dCas9p300 was used to induce elevated levels of CFTR from the endogenous locus in ΔF508/ΔF508 human bronchial epithelial cells. Ussing chamber studies demonstrated a synergistic increase in ion transport in response to CRISPR-induced expression of ΔF508 CFTR mRNA along with VX809 treatment.

CONCLUSIONS

CRISPR/dCas9-based epigenome-editing provides a previously unexplored tool for interrogating CFTR enhancer function. Here, we demonstrate that therapeutic interventions that increase the expression of CFTR may improve the efficacy of CFTR modulators. A better understanding CFTR regulatory mechanisms could uncover novel therapeutic interventions for the development of cystic fibrosis therapies.

DNA Methylation-Mediated Low Expression of CFTR Stimulates the Progression of Lung Adenocarcinoma

In recent years, the mortality rate of lung adenocarcinoma (LUAD) is persistently increasing, which has already caused a huge impact on human living standards. Hence, there is an urgent need to probe the molecular mechanism of LUAD progression, so as to disclose prognostic and diagnostic markers for patients with LUAD. Methylation 450 K data and mRNA expression data of LUAD were obtained via bioinformatics analysis to screen methylation-driven genes. The expression of the target gene was detected through qRT-PCR, while the methylation level was evaluated via methylation-specific PCR (MSP). The impact of the gene on cell proliferation, migration, invasion, apoptosis and cell cycle was measured through CCK-8, wound healing, Transwell invasion assay, and flow cytometry. CFTR was defined by bioinformatics analysis as the target gene for this study. qRT-PCR revealed that CFTR was lowly expressed in LUAD cells. MSP displayed that the CFTR promoter region in LUAD cells was hypermethylated, and demethylation could pronouncedly increase the level of CFTR mRNA in LUAD cells. Cell biological functional experiments exhibited that CFTR hindered cell proliferation, migration, and invasion, fostered cell apoptosis of LUAD, and blocked the cell cycle in G2-M phase. CFTR was hypermethylated in LUAD, which mediated the low expression of CFTR in LUAD to stimulate the progression of LUAD.

OTX2 regulates CFTR expression during endoderm differentiation and occupies 3' cis-regulatory elements

BACKGROUND

Cell-specific and developmental mechanisms contribute to expression of the cystic fibrosis transmembrane conductance regulator (CFTR) gene; however, its developmental regulation is poorly understood. Here we use human induced pluripotent stem cells differentiated into pseudostratified airway epithelial cells to study these mechanisms.

RESULTS

Changes in gene expression and open chromatin profiles were investigated by RNA-seq and ATAC-seq, and revealed that alterations in CFTR expression are associated with differences in stage-specific open chromatin. Additionally, two novel open chromatin regions, at +19.6 kb and +22.6 kb 3' to the CFTR translational stop signal, were observed in definitive endoderm (DE) cells, prior to an increase in CFTR expression in anterior foregut endoderm (AFE) cells. Chromatin studies in DE and AFE cells revealed enrichment of active enhancer marks and occupancy of OTX2 at these sites in DE cells. Loss of OTX2 in DE cells alters histone modifications across the CFTR locus and results in a 2.5-fold to 5-fold increase in CFTR expression. However, deletion of the +22.6 kb site alone does not affect CFTR expression in DE or AFE cells.

CONCLUSIONS

These results suggest that a network of interacting cis-regulatory elements recruit OTX2 to the locus to impact CFTR expression during early endoderm differentiation.

CFTR Cooperative Cis-Regulatory Elements in Intestinal Cells

About 8% of the human genome is covered with candidate cis-regulatory elements (cCREs). Disruptions of CREs, described as "cis-ruptions" have been identified as being involved in various genetic diseases. Thanks to the development of chromatin conformation study techniques, several long-range cystic fibrosis transmembrane conductance regulator (CFTR) regulatory elements were identified, but the regulatory mechanisms of the CFTR gene have yet to be fully elucidated. The aim of this work is to improve our knowledge of the CFTR gene regulation, and to identity factors that could impact the CFTR gene expression, and potentially account for the variability of the clinical presentation of cystic fibrosis as well as CFTR-related disorders. Here, we apply the robust GWAS3D score to determine which of the CFTR introns could be involved in gene regulation. This approach highlights four particular CFTR introns of interest. Using reporter gene constructs in intestinal cells, we show that two new introns display strong cooperative effects in intestinal cells. Chromatin immunoprecipitation analyses further demonstrate fixation of transcription factors network. These results provide new insights into our understanding of the CFTR gene regulation and allow us to suggest a 3D CFTR locus structure in intestinal cells. A better understand of regulation mechanisms of the CFTR gene could elucidate cases of patients where the phenotype is not yet explained by the genotype. This would thus help in better diagnosis and therefore better management. These cis-acting regions may be a therapeutic challenge that could lead to the development of specific molecules capable of modulating gene expression in the future.

Chromatin regulatory dynamics of early human small intestinal development using a directed differentiation model

The establishment of the small intestinal (SI) lineage during human embryogenesis ensures functional integrity of the intestine after birth. The chromatin dynamics that drive SI lineage formation and regional patterning in humans are essentially unknown. To fill this knowledge void, we apply a cutting-edge genomic technology to a state-of-the-art human model of early SI development. Specifically, we leverage chromatin run-on sequencing (ChRO-seq) to define the landscape of active promoters, enhancers and gene bodies across distinct stages of directed differentiation of human pluripotent stem cells into SI spheroids with regional specification. Through comprehensive ChRO-seq analysis we identify candidate stage-specific chromatin activity states, novel markers and enhancer hotspots during the directed differentiation. Moreover, we propose a detailed transcriptional network associated with SI lineage formation or regional patterning. Our ChRO-seq analyses uncover a previously undescribed pattern of enhancer activity and transcription at HOX gene loci underlying SI regional patterning. We also validated this unique HOX dynamics by the analysis of single cell RNA-seq data from human fetal SI. Overall, the results lead to a new proposed working model for the regulatory underpinnings of human SI development, thereby adding a novel dimension to the literature that has relied almost exclusively on non-human models.

Hepatocyte nuclear factor-4α regulates expression of the serotonin transporter in intestinal epithelial cells

The serotonin transporter (SERT) functions to regulate the availability of serotonin (5-HT) in the brain and intestine. An intestine-specific mRNA variant arising from a unique transcription start site and alternative promoter in the SERT gene has been identified (iSERT; spanning exon 1C). A decrease in SERT is implicated in several gut disorders, including inflammatory bowel diseases (IBD). However, little is known about mechanisms regulating the iSERT variant, and a clearer understanding is warranted for targeting SERT for the treatment of gut disorders. The current studies examined the expression of iSERT across different human intestinal regions and investigated its regulation by HNF4α (hepatic nuclear factor-4α), a transcription factor important for diverse cellular functions. iSERT mRNA abundance was highest in the human ileum and Caco-2 cell line. iSERT mRNA expression was downregulated by loss of HNF4α (but not HNF1α, HNF1β, or FOXA1) in Caco-2 cells. Overexpression of HNF4α increased iSERT mRNA concomitant with an increase in SERT protein. Progressive promoter deletion and site-directed mutagenesis revealed that the HNF4α response element spans nucleotides -1,163 to -1150 relative to the translation start site. SERT mRNA levels in the intestine were drastically reduced in the intestine-specific HNF4α-knockout mice relative to HNF4αFL/FL mice. Both HNF4α and SERT mRNA levels were also downregulated in mouse model of ileitis (SAMP) compared with AKR control mice. These results establish the transcriptional regulation of iSERT at the gut-specific internal promoter (hSERTp2) and have identified HNF4α as a critical modulator of basal SERT expression in the intestine.

The FOXA1 transcriptional network coordinates key functions of primary human airway epithelial cells

The differentiated functions of the human airway epithelium are coordinated by a complex network of transcription factors. These include the pioneer factors Forkhead box A1 and A2 (FOXA1 and FOXA2), which are well studied in several tissues, but their role in airway epithelial cells is poorly characterized. Here, we define the cistrome of FOXA1 and FOXA2 in primary human bronchial epithelial (HBE) cells by chromatin immunoprecipitation with deep-sequencing (ChIP-seq). Next, siRNA-mediated depletion of each factor is used to investigate their transcriptome by RNA-seq. We found that, as predicted from their DNA-binding motifs, genome-wide occupancy of the two factors showed substantial overlap; however, their global impact on gene expression differed. FOXA1 is an abundant transcript in HBE cells, while FOXA2 is expressed at low levels, and both these factors likely exhibit autoregulation and cross-regulation. FOXA1 regulated loci are involved in cell adhesion and the maintenance of epithelial cell identity, particularly through repression of genes associated with epithelial to mesenchymal transition (EMT). FOXA1 also directly targets other transcription factors with a known role in the airway epithelium such as SAM-pointed domain-containing Ets-like factor (SPDEF). The intersection of the cistrome and transcriptome for FOXA1 revealed enrichment of genes involved in epithelial development and tissue morphogenesis. Moreover, depletion of FOXA1 was shown to reduce the transepithelial resistance of HBE cells, confirming the role of this factor in maintaining epithelial barrier integrity.

Genetic variation in CFTR and modifier loci may modulate cystic fibrosis disease severity

In patients with cystic fibrosis (CF), genetic variants within and outside the CFTR locus contribute to the variability of the disease severity. CFTR transcription is tightly regulated by cis-regulatory elements (CREs) that control the three-dimensional structure of the locus, chromatin accessibility and transcription factor recruitment. Variants within these CREs may contribute to the pathophysiology and to the phenotypic heterogeneity by altering CFTR transcript abundance. In addition to the CREs, variants outside the CFTR locus, namely "modifiers genes", may also be associated with the clinical variability. This review addresses variants at the CFTR locus itself and CFTR CREs, together with the outcomes of the latest modifier gene studies with respect to the different CF phenotypes.

TCam-2 Cells Deficient for SOX2 and FOXA2 Are Blocked in Differentiation and Maintain a Seminoma-Like Cell Fate In Vivo

Testicular germ cell tumors (GCTs) are very common in young men and can be stratified into seminomas and non-seminomas. While seminomas share a similar gene expression and epigenetic profile with primordial germ cells, the stem cell population of the non-seminomas, the embryonal carcinoma (EC), resembles malignant embryonic stem cells. Thus, ECs are able to differentiate into cells of all three germ layers (teratomas) and even extra-embryonic-tissue-like cells (yolk-sac tumor, choriocarcinoma). In the last years, we demonstrated that the cellular microenvironment considerably influences the plasticity of seminomas (TCam-2 cells). Upon a microenvironment-triggered inhibition of the BMP signaling pathway in vivo (murine flank or brain), seminomatous TCam-2 cells reprogram to an EC-like cell fate. We identified SOX2 as a key factor activated upon BMP inhibition mediating the reprogramming process by regulating pluripotency, reprogramming and epigenetic factors. Indeed, CRISPR/Cas9 SOX2-deleted TCam-2 cells were able to maintain a seminoma-cell fate in vivo for about six weeks, but after six weeks in vivo still small sub-populations initiated differentiation. Closer analyses of these differentiated clusters suggested that the pioneer factor FOXA2 might be the driving force behind this induction of differentiation, since many FOXA2 interacting genes and differentiation factors like AFP, EOMES, CDX1, ALB, HAND1, DKK, DLK1, MSX1 and PITX2 were upregulated. In this study, we generated TCam-2 cells double-deficient for SOX2 and FOXA2 using the CRISPR/Cas9 technique and xenografted those cells into the flank of nude mice. Upon loss of SOX2 and FOXA2, TCam-2 maintained a seminoma cell fate for at least twelve weeks, demonstrating that both factors are key players in the reprogramming to an EC-like cell fate. Therefore, our study adds an important piece to the puzzle of GCT development and plasticity, providing interesting insights in what can be expected in a patient, when GCT cells are confronted with different microenvironments.

From Pioneer to Repressor: Bimodal foxd3 Activity Dynamically Remodels Neural Crest Regulatory Landscape In Vivo

The neural crest (NC) is a transient embryonic stem cell-like population characterized by its multipotency and broad developmental potential. Here, we perform NC-specific transcriptional and epigenomic profiling of foxd3-mutant cells in vivo to define the gene regulatory circuits controlling NC specification. Together with global binding analysis obtained by foxd3 biotin-ChIP and single cell profiles of foxd3-expressing premigratory NC, our analysis shows that, during early steps of NC formation, foxd3 acts globally as a pioneer factor to prime the onset of genes regulating NC specification and migration by re-arranging the chromatin landscape, opening cis-regulatory elements and reshuffling nucleosomes. Strikingly, foxd3 then gradually switches from an activator to its well-described role as a transcriptional repressor and potentially uses differential partners for each role. Taken together, these results demonstrate that foxd3 acts bimodally in the neural crest as a switch from “permissive” to “repressive” nucleosome and chromatin organization to maintain multipotency and define cell fates.

•

FoxD3 primes neural crest specification by modulating distal enhancers

•

FoxD3 represses a number of neural crest migration and differentiation genes

•

In neural crest, FoxD3 acts to switch chromatin from “permissive” to “repressive”

•

Distinctive gene regulatory mechanisms underlie the bimodal action of FoxD3

Cell-Selective Regulation of CFTR Gene Expression: Relevance to Gene Editing Therapeutics

The cystic fibrosis transmembrane conductance regulator (CFTR) gene is an attractive target for gene editing approaches, which may yield novel therapeutic approaches for genetic diseases such as cystic fibrosis (CF). However, for gene editing to be effective, aspects of the three-dimensional (3D) structure and cis-regulatory elements governing the dynamic expression of CFTR need to be considered. In this review, we focus on the higher order chromatin organization required for normal CFTR locus function, together with the complex mechanisms controlling expression of the gene in different cell types impaired by CF pathology. Across all cells, the CFTR locus is organized into an invariant topologically associated domain (TAD) established by the architectural proteins CCCTC-binding factor (CTCF) and cohesin complex. Additional insulator elements within the TAD also recruit these factors. Although the CFTR promoter is required for basal levels of expression, cis-regulatory elements (CREs) in intergenic and intronic regions are crucial for cell-specific and temporal coordination of CFTR transcription. These CREs are recruited to the promoter through chromatin looping mechanisms and enhance cell-type-specific expression. These features of the CFTR locus should be considered when designing gene-editing approaches, since failure to recognize their importance may disrupt gene expression and reduce the efficacy of therapies.

FOXA1 reprograms the TGF-β-stimulated transcriptional program from a metastasis promoter to a tumor suppressor in nasopharyngeal carcinoma

Nasopharyngeal carcinoma (NPC) is a unique subtype of head and neck squamous carcinoma that is notorious for its high metastatic potential. In this study, we reported that FOXA1 protein was decreased in NPC cells. Loss of FOXA1 is associated with lymph node metastasis and poor prognosis. Silencing FOXA1 in NP69 and C666-1 NPC cells accelerated cell proliferation and migration, while re-expression of FOXA1 has opposite effects. Microarray and RNA-seq analysis revealed that re-expression of FOXA1 in NPC cells reprogrammed the TGF-β-stimulated transcription program, which is characterized by promotion of TGF-β-inducible tumor-suppressive targets but repression of TGF-β-inducible oncogenes expression in NPC cells, leading to restoration of NPC cell sensitivity to TGF-β's growth-inhibitory effect. BAMBI, a TGF-β responsive tumor suppressor, was induced by FOXA1 in NPC cells. FOXA1 binding on the BAMBI gene facilitated SMAD2/3 binding to the BAMBI promoter via increasing BAMBI associated H3K4me1 and H3K27ac modification. Enforced expression of BAMBI in NPC cells suppressed cell proliferation and invasiveness. Our data suggested that FOXA1 is a master factor in controlling the TGF-β-stimulated transcriptome and a regulator of TGF-β biological functions in NPC oncogenesis.

Screening for Regulatory Variants in 460 kb Encompassing the CFTR Locus in Cystic Fibrosis Patients

It is estimated that up to 5% of cystic fibrosis transmembrane conductance regulator (CFTR) pathogenic alleles are unidentified. Some of these errors may lie in noncoding regions of the locus and affect gene expression. To identify regulatory element variants in the CFTR locus, SureSelect targeted enrichment of 460 kb encompassing the gene was optimized to deep sequence genomic DNA from 80 CF patients with an unequivocal clinical diagnosis but only one or no CFTR-coding region pathogenic variants. Bioinformatics tools were used to identify sequence variants and predict their impact, which were then assayed in transient reporter gene luciferase assays. The effect of five variants in the CFTR promoter and four in an intestinal enhancer of the gene were assayed in relevant cell lines. The initial analysis of sequence data revealed previously known CF-causing variants, validating the robustness of the SureSelect design, and showed that 85 of 160 CF alleles were undefined. Of a total 1737 variants revealed across the extended 460-kb CFTR locus, 51 map to known CFTR cis-regulatory elements, and many of these are predicted to alter transcription factor occupancy. Four promoter variants and all those in the intestinal enhancer significantly repress reporter gene activity. These data suggest that CFTR regulatory elements may harbor novel CF disease-causing variants that warrant further investigation, both for genetic screening protocols and functional assays.

The pioneer factor activity of c-Myb involves recruitment of p300 and induction of histone acetylation followed by acetylation-induced chromatin dissociation

Background

The concept of pioneer transcription factors is emerging as an essential part of the epigenetic regulation, taking place during cell development and differentiation. However, the precise molecular mechanism underlying pioneer factor activity remains poorly understood. We recently reported that the transcription factor c-Myb acts as a pioneer factor in haematopoiesis, and a point mutation in its DNA binding domain, D152V, is able to abrogate this function.

Results

Here, we show that specific histone modifications, including H3K27ac, prevent binding of c-Myb to histone tails, representing a novel effect of histone modifications, namely restricting binding of a pioneer factor to chromatin. Furthermore, we have taken advantage of the pioneer-defect D152V mutant to investigate mechanisms of c-Myb’s pioneer factor activity. We show that c-Myb-dependent transcriptional activation of a gene in inaccessible chromatin relies on c-Myb binding to histones, as well as on c-Myb interacting with the histone acetyltransferase p300. ChIP assays show that both wild type and the D152V mutant of c-Myb bind to a selected target gene at its promoter and enhancer, but only wild-type c-Myb causes opening and activation of the locus. Enhancement of histone acetylation restores activation of the same gene in the absence of c-Myb, suggesting that facilitating histone acetylation is a crucial part of the pioneer factor function of c-Myb.

Conclusions

We suggest a pioneer factor model in which c-Myb binds to regions of closed chromatin and then recruits histone acetyltransferases. By binding to histones, c-Myb facilitates histone acetylation, acting as a cofactor for p300 at c-Myb bound sites. The resulting H3K27ac leads to chromatin opening and detachment of c-Myb from the acetylated chromatin. We propose that the latter phenomenon, acetylation-induced chromatin dissociation, represents a mechanism for controlling the dynamics of pioneer factor binding to chromatin.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0208-y) contains supplementary material, which is available to authorized users.

A transcription factor network represses CFTR gene expression in airway epithelial cells

Mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene cause the inherited disorder cystic fibrosis (CF). Lung disease is the major cause of CF morbidity, though CFTR expression levels are substantially lower in the airway epithelium than in pancreatic duct and intestinal epithelia, which also show compromised function in CF. Recently developed small molecule therapeutics for CF are highly successful for one specific CFTR mutation and have a positive impact on others. However, the low abundance of CFTR transcripts in the airway limits the opportunity for drugs to correct the defective substrate. Elucidation of the transcriptional mechanisms for the CFTR locus has largely focused on intragenic and intergenic tissue-specific enhancers and their activating trans-factors. Here, we investigate whether the low CFTR levels in the airway epithelium result from the recruitment of repressive proteins directly to the locus. Using an siRNA screen to deplete ∼1500 transcription factors (TFs) and associated regulatory proteins in Calu-3 lung epithelial cells, we identified nearly 40 factors that upon depletion elevated CFTR mRNA levels more than 2-fold. A subset of these TFs was validated in primary human bronchial epithelial cells. Among the strongest repressors of airway expression of CFTR were Krüppel-like factor 5 and Ets homologous factor, both of which have pivotal roles in the airway epithelium. Depletion of these factors, which are both recruited to an airway-selective cis-regulatory element at -35 kb from the CFTR promoter, improved CFTR production and function, thus defining novel therapeutic targets for enhancement of CFTR.

SOX2 is essential for in vivo reprogramming of seminoma-like TCam-2 cells to an embryonal carcinoma-like fate

Type II germ cell cancers (GCC) are divided into seminomas, which are highly similar to primordial germ cells and embryonal carcinomas (EC), often described as malignant counterparts to embryonic stem cells.

Previously, we demonstrated that the development of GCCs is a highly plastic process and strongly influenced by the microenvironment. While orthotopic transplantation into the testis promotes seminomatous growth of the seminoma-like cell line TCam-2, ectopic xenotransplantation into the flank initiates reprogramming into an EC-like fate.

During this reprogramming, BMP signaling is inhibited, leading to induction of NODAL signaling, upregulation of pluripotency factors and downregulation of seminoma markers, like SOX17. The pluripotency factor and EC-marker SOX2 is strongly induced.

Here, we adressed the molecular role of SOX2 in this reprogramming. Using CRISPR/Cas9-mediated genome-editing, we established SOX2-deficient TCam-2 cells. Xenografting of SOX2-deficient cells into the flank of nude mice resulted in maintenance of a seminoma-like fate, indicated by the histology and expression of OCT3/4, SOX17, TFAP2C, PRDM1 and PRAME. In SOX2-deficient cells, BMP signaling is inhibited, but NODAL signaling is not activated. Thus, SOX2 appears to be downstream of BMP signaling but upstream of NODAL activation. So, SOX2 is an essential factor in acquiring an EC-like cell fate from seminomas.

A small population of differentiated cells was identified resembling a mixed non-seminoma. Analyses of these cells revealed downregulation of the pluripotency and seminoma markers OCT3/4, SOX17, PRDM1 and TFAP2C. In contrast, the pioneer factor FOXA2 and its target genes were upregulated, suggesting that FOXA2 might play an important role in induction of non-seminomatous differentiation.

Mechanistic Insights Into the Interaction Between Transcription Factors and Epigenetic Modifications and the Contribution to the Development of Obesity

Objective: The development of obesity is inseparable from genetic and epigenetic factors, and transcription factors (TFs) play an essential role in these two mechanisms. This review analyzes the interaction of TFs with epigenetic modifications and the epigenetic mechanisms underlying peroxisome proliferator-activated receptor (PPAR)γ, an important transcription factor, in the development of obesity. Methods: We describe the relationship between TFs and different epigenetic modifications and illustrate the several mechanisms described. Next, we summarize the epigenetic mechanisms of PPARs, an important class of transcription factors involved in obesity, that induce obesity with different triggering factors. Finally, we discuss the mechanisms of epigenetic modification of PPAR-related ligands in lipid metabolism and propose future avenues of research. Results: TFs participate in epigenetic modifications in different forms, causing changes in gene expression. The interactions between the different epigenetic modifications and PPARs are important biological developments that affect fat tissue differentiation, lipogenesis, and lipid metabolism, thereby inducing or inhibiting the development of obesity. We then highlight the need for more research to understand the role of epigenetic modifications and PPARs. Conclusions: Epigenetic mechanisms involved in the regulation of PPARs may be excellent therapeutic targets for obesity treatment. However, there is a need for a deeper understanding of how PPARs and other obesity-related transcription factors interact with epigenetic modifications.

ABCMdb reloaded: updates on mutations in ATP binding cassette proteins

ABC (ATP-Binding Cassette) proteins with altered function are responsible for numerous human diseases. To aid the selection of positions and amino acids for ABC structure/function studies we have generated a database, ABCMdb (Gyimesi et al., ABCMdb: a database for the comparative analysis of protein mutations in ABC transporters, and a potential framework for a general application. Hum Mutat 2012; 33:1547–1556.), with interactive tools. The database has been populated with mentions of mutations extracted from full text papers, alignments and structural models. In the new version of the database we aimed to collect the effect of mutations from databases including ClinVar. Because of the low number of available data, even in the case of the widely studied disease-causing ABC proteins, we also included the possible effects of mutations based on SNAP2 and PROVEAN predictions. To aid the interpretation of variations in non-coding regions, the database was supplemented with related DNA level information. Our results emphasize the importance of in silico predictions because of the sparse information available on variants and suggest that mutations at analogous positions in homologous ABC proteins have a strong predictive power for the effects of mutations. Our improved ABCMdb advances the design of both experimental studies and meta-analyses in order to understand drug interactions of ABC proteins and the effects of mutations on functional expression.

Database URL:http://abcm2.hegelab.org

A Prostate Cancer Risk Element Functions as a Repressive Loop that Regulates HOXA13

Prostate cancer (PCa) is the leading cancer among men in the United States, with genetic factors contributing to ∼42% of the susceptibility to PCa. We analyzed a PCa risk region located at 7p15.2 to gain insight into the mechanisms by which this noncoding region may affect gene regulation and contribute to PCa risk. We performed Hi-C analysis and demonstrated that this region has long-range interactions with the HOXA locus, located ∼873 kb away. Using the CRISPR/Cas9 system, we deleted a 4-kb region encompassing several PCa risk-associated SNPs and performed RNA-seq to investigate transcriptomic changes in prostate cells lacking the regulatory element. Our results suggest that the risk element affects the expression of HOXA13 and HOTTIP, but not other genes in the HOXA locus, via a repressive loop. Forced expression of HOXA13 was performed to gain further insight into the mechanisms by which this risk element affects PCa risk.

FOXA transcriptional factor modulates insect susceptibility to Bacillus thuringiensis Cry1Ac toxin by regulating the expression of toxin-receptor ABCC2 and ABCC3 genes

Cry toxins produced by Bacillus thuringiensis (Bt) are insecticidal proteins widely used in insect control. Recently, it was shown that ATP-binding cassette transporter proteins (ABC) such as ABCC2, ABCC3, ABCG1 and ABCA2 are implicated in the insecticidal action of Cry toxins as putative receptors. However, the transcriptional regulators involved in the expression of ABC transporter genes remain unknown. Sequence analysis of promoter regions of ABCC2 gene from Helicoverpa armigera and ABCC3 gene from Spodoptera litura Sl-HP cultured cells, revealed the potential participation of Forkhead box protein A (FOXA), a transcription factor that regulates the expression of genes through remodeling chromatin. To determine if FOXA was involved in regulating expression of ABCC2 and ABCC3 genes, the expression of FOXA, ABCC2 and ABCC3 was compared in Sl-HP cells that are sensitive to Cry1Ac toxin with those in S. frugiperda Sf9 cells that are not sensitive to the toxin. Expression levels of those genes were significantly higher in Sl-HP than in Sf9 cells. Transient expression of FOXA in Sf9 cells activated ABCC2 and ABCC3 transcription, which directly correlated with enhanced Cry1Ac-susceptibility in these cells. Silencing of FOXA gene expression by RNAi in H. armigera larvae resulted in a decreased expression of ABCC2 and ABCC3 without affecting expression of other Cry toxin receptor genes such as alkaline phosphatase, aminopeptidase or cadherin. Silencing of FOXA gene expression also resulted in a Cry1Ac-tolerant phenotype since lower mortality and higher pupation rate were observed in diet containing Cry1Ac protoxin in comparison with the control group. These results demonstrate that FOXA up-regulates expression of the Cry1Ac-toxin receptor ABCC2 and ABCC3 genes, and that lower FOXA expression correlates with tolerance to Cry toxin in cell lines and in lepidopteran larvae.

Ets homologous factor (EHF) has critical roles in epithelial dysfunction in airway disease

The airway epithelium forms a barrier between the internal and external environments. Epithelial dysfunction is critical in the pathology of many respiratory diseases, including cystic fibrosis. Ets homologous factor (EHF) is a key member of the transcription factor network that regulates gene expression in the airway epithelium in response to endogenous and exogenous stimuli. EHF, which has altered expression in inflammatory states, maps to the 5' end of an intergenic region on Chr11p13 that is implicated as a modifier of cystic fibrosis airway disease. Here we determine the functions of EHF in primary human bronchial epithelial (HBE) cells and relevant airway cell lines. Using EHF ChIP followed by deep sequencing (ChIP-seq) and RNA sequencing after EHF depletion, we show that EHF targets in HBE cells are enriched for genes involved in inflammation and wound repair. Furthermore, changes in gene expression impact cell phenotype because EHF depletion alters epithelial secretion of a neutrophil chemokine and slows wound closure in HBE cells. EHF activates expression of the SAM pointed domain-containing ETS transcription factor, which contributes to goblet cell hyperplasia. Our data reveal a critical role for EHF in regulating epithelial function in lung disease.

Feature-based classification of human transcription factors into hypothetical sub-classes related to regulatory function

Background

Transcription factors are key proteins in the regulation of gene transcription. An important step in this process is the opening of chromatin in order to make genomic regions available for transcription. Data on DNase I hypersensitivity has previously been used to label a subset of transcription factors as Pioneers, Settlers and Migrants to describe their potential role in this process. These labels represent an interesting hypothesis on gene regulation and possibly a useful approach for data analysis, and therefore we wanted to expand the set of labeled transcription factors to include as many known factors as possible. We have used a well-annotated dataset of 1175 transcription factors as input to supervised machine learning methods, using the subset with previously assigned labels as training set. We then used the final classifier to label the additional transcription factors according to their potential role as Pioneers, Settlers and Migrants. The full set of labeled transcription factors was used to investigate associated properties and functions of each class, including an analysis of interaction data for transcription factors based on DNA co-binding and protein-protein interactions. We also used the assigned labels to analyze a previously published set of gene lists associated with a time course experiment on cell differentiation.

Results

The analysis showed that the classification of transcription factors with respect to their potential role in chromatin opening largely was determined by how they bind to DNA. Each subclass of transcription factors was enriched for properties that seemed to characterize the subclass relative to its role in gene regulation, with very general functions for Pioneers, whereas Migrants to a larger extent were associated with specific processes. Further analysis showed that the expanded classification is a useful resource for analyzing other datasets on transcription factors with respect to their potential role in gene regulation. The analysis of transcription factor interaction data showed complementary differences between the subclasses, where transcription factors labeled as Pioneers often interact with other transcription factors through DNA co-binding, whereas Migrants to a larger extent use protein-protein interactions. The analysis of time course data on cell differentiation indicated a shift in the regulatory program associated with Pioneer-like transcription factors during differentiation.

Conclusions

The expanded classification is an interesting resource for analyzing data on gene regulation, as illustrated here on transcription factor interaction data and data from a time course experiment. The potential regulatory function of transcription factors seems largely to be determined by how they bind DNA, but is also influenced by how they interact with each other through cooperativity and protein-protein interactions.

Electronic supplementary material

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

Differential contribution of cis-regulatory elements to higher order chromatin structure and expression of the CFTR locus

Higher order chromatin structure establishes domains that organize the genome and coordinate gene expression. However, the molecular mechanisms controlling transcription of individual loci within a topological domain (TAD) are not fully understood. The cystic fibrosis transmembrane conductance regulator (CFTR) gene provides a paradigm for investigating these mechanisms.CFTR occupies a TAD bordered by CTCF/cohesin binding sites within which are cell-type-selective cis-regulatory elements for the locus. We showed previously that intronic and extragenic enhancers, when occupied by specific transcription factors, are recruited to the CFTR promoter by a looping mechanism to drive gene expression. Here we use a combination of CRISPR/Cas9 editing of cis-regulatory elements and siRNA-mediated depletion of architectural proteins to determine the relative contribution of structural elements and enhancers to the higher order structure and expression of the CFTR locus. We found the boundaries of the CFTRTAD are conserved among diverse cell types and are dependent on CTCF and cohesin complex. Removal of an upstream CTCF-binding insulator alters the interaction profile, but has little effect on CFTR expression. Within the TAD, intronic enhancers recruit cell-type selective transcription factors and deletion of a pivotal enhancer element dramatically decreases CFTR expression, but has minor effect on its 3D structure.

MicroRNA-99a and 100 mediated upregulation of FOXA1 in bladder cancer

Urothelial cell carcinoma of the bladder (UCC) is a common disease often characterized by FGFR3 dysregulation. Whilst upregulation of this oncogene occurs most frequently in low-grade non-invasive tumors, recent data reveal increased FGFR3 expression characterizes a common sub-type of invasive UCC sharing molecular similarities with breast cancer. These similarities include upregulation of the FOXA1 transcription factor and reduced expression of microRNAs-99a/100. We have previously identified direct regulation of FGFR3 by these two microRNAs and now search for further targets. Using a microarray meta-database we find potential FOXA1 regulation by microRNAs-99a/100. We confirm direct targeting of the FOXA1 3′UTR by microRNAs-99a/100 and also potential indirect regulation through microRNA-485-5p/SOX5/JUN-D/FOXL1 and microRNA-486/FOXO1a. In 292 benign and malignant urothelial samples, we find an inverse correlation between the expression of FOXA1 and microRNAs-99a/100 (r=−0.33 to −0.43, p<0.05). As for FGFR3 in UCC, tumors with high FOXA1 expression have lower rates of progression than those with low expression (Log rank p=0.009). Using global gene expression and CpG methylation profiling we find genotypic consequences of FOXA1 upregulation in UCC. Genetic changes are associated with regional hypomethylation, occur near FOXA1 binding sites, and mirror gene expression changes previously reported in FGFR3 mutant-UCC. These include gene silencing through aberrant hypermethylation (e.g. IGFBP3) and affect genes characterizing breast cancer sub-types (e.g. ERBB2). In conclusion, we have identified microRNAs-99a/100 mediate a direct relationship between FGFR3 and FOXA1 and potentially facilitate cross talk between these pathways in UCC.

Analysis of long-range interactions in primary human cells identifies cooperative CFTR regulatory elements

A mechanism by which control DNA elements regulate transcription over large linear genomic distances is by achieving close physical proximity with genes, and looping of the intervening chromatin paths. Alterations of such regulatory 'chromatin looping' systems are likely to play a critical role in human genetic disease at large. Here, we studied the spatial organization of a ≈790 kb locus encompassing the cystic fibrosis transmembrane conductance regulator (CFTR) gene. Dysregulation of CFTR is responsible for cystic fibrosis, which is the most common lethal genetic disorder in Caucasian populations. CFTR is a relatively large gene of 189 kb with a rather complex tissue-specific and temporal expression profile. We used chromatin conformation at the CFTR locus to identify new DNA sequences that regulate its transcription. By comparing 5C chromatin interaction maps of the CFTR locus in expressing and non-expressing human primary cells, we identified several new contact points between the CFTR promoter and its surroundings, in addition to regions featuring previously described regulatory elements. We demonstrate that two of these novel interacting regions cooperatively increase CFTR expression, and suggest that the new enhancer elements located on either side of the gene are brought together through chromatin looping via CTCF.

Chromatin Dynamics in the Regulation of CFTR Expression

The contribution of chromatin dynamics to the regulation of human disease-associated loci such as the cystic fibrosis transmembrane conductance regulator (CFTR) gene has been the focus of intensive experimentation for many years. Recent technological advances in the analysis of transcriptional mechanisms across the entire human genome have greatly facilitated these studies. In this review we describe the complex machinery of tissue-specific regulation of CFTR expression, and put earlier observations in context by incorporating them into datasets generated by the most recent genomics methods. Though the gene promoter is required for CFTR expression, cell-type specific regulatory elements are located elsewhere in the gene and in flanking intergenic regions. Probably within its own topological domain established by the architectural proteins CTCF and cohesin, the CFTR locus utilizes chromatin dynamics to remodel nucleosomes, recruit cell-selective transcription factors, and activate intronic enhancers. These cis-acting elements are then brought to the gene promoter by chromatin looping mechanisms, which establish long-range interactions across the locus. Despite its complexity, the CFTR locus provides a paradigm for elucidating the critical role of chromatin dynamics in the transcription of individual human genes.

FOXA2 regulates a network of genes involved in critical functions of human intestinal epithelial cells

The forkhead box A (FOXA) family of pioneer transcription factors is critical for the development of many endoderm-derived tissues. Their importance in regulating biological processes in the lung and liver is extensively characterized, though much less is known about their role in intestine. Here we investigate the contribution of FOXA2 to coordinating intestinal epithelial cell function using postconfluent Caco2 cells, differentiated into an enterocyte-like model. FOXA2 binding sites genome-wide were determined by ChIP-seq and direct targets of the factor were validated by ChIP-qPCR and siRNA-mediated depletion of FOXA1/2 followed by RT-qPCR. Peaks of FOXA2 occupancy were frequent at loci contributing to gene ontology pathways of regulation of cell migration, cell motion, and plasma membrane function. Depletion of both FOXA1 and FOXA2 led to a significant reduction in the expression of multiple transmembrane proteins including ion channels and transporters, which form a network that is essential for maintaining normal ion and solute transport. One of the targets was the adenosine A2B receptor, and reduced receptor mRNA levels were associated with a functional decrease in intracellular cyclic AMP. We also observed that 30% of FOXA2 binding sites contained a GATA motif and that FOXA1/A2 depletion reduced GATA-4, but not GATA-6 protein levels. These data show that FOXA2 plays a pivotal role in regulating intestinal epithelial cell function. Moreover, that the FOXA and GATA families of transcription factors may work cooperatively to regulate gene expression genome-wide in the intestinal epithelium.

Oxidative stress regulates CFTR gene expression in human airway epithelial cells through a distal antioxidant response element

Cystic fibrosis transmembrane conductance regulator gene (CFTR) expression in human airway epithelial cells involves the recruitment of distal cis-regulatory elements, which are associated with airway-selective DNase hypersensitive sites at -44 kb and -35 kb from the gene. The -35-kb site encompasses an enhancer that is regulated by the immune mediators interferon regulatory factor 1 and 2 and by nuclear factor Y. Here we investigate the -44-kb element, which also has enhancer activity in vitro in airway epithelial cells but is inactive in intestinal epithelial cells. This site contains an antioxidant response element (ARE) that plays a critical role in its function in airway cell lines and primary human bronchial epithelial cells. The natural antioxidant sulforaphane (SFN) induces nuclear translocation of nuclear factor, erythroid 2-like 2 (Nrf2), a transcription factor that regulates genes with AREs in their promoters, many of which are involved in response to injury. Under normal conditions, the -44-kb ARE is occupied by the repressor BTB and CNC homology 1, basic leucine zipper transcription factor (Bach1), and v-Maf avian musculoaponeurotic fibrosarcoma oncogene homolog K (MafK) heterodimers. After 2 hours of SFN treatment, Nrf2 displaces these repressive factors and activates CFTR expression. Site-directed mutagenesis shows that both the ARE and an adjacent NF-κB binding site are required for activation of the -44-kb element in airway epithelial cells. Moreover, this element is functionally linked to the -35-kb enhancer in modulating CFTR expression in response to environmental stresses in the airway.

FOXA2 attenuates the epithelial to mesenchymal transition by regulating the transcription of E-cadherin and ZEB2 in human breast cancer

The Forkhead Box A2 (FOXA2) transcription factor is required for embryonic development and for normal functions of multiple adult tissues, in which the maintained expression of FOXA2 is usually related to preventing the progression of malignant transformation. In this study, we found that FOXA2 prevented the epithelial to mesenchymal transition (EMT) in human breast cancer. We observed a strong correlation between the expression levels of FOXA2 and the epithelial phenotype. Knockdown of FOXA2 promoted the mesenchymal phenotype, whereas stable overexpression of FOXA2 attenuated EMT in breast cancer cells. FOXA2 was found to endogenously bind to and stimulate the promoter of E-cadherin that is crucial for epithelial phenotype of the tumor cells. Meanwhile, FOXA2 prevented EMT of breast cancer cells by repressing the expression of EMT-related transcription factor ZEB2 through recruiting a transcriptional corepressor TLE3 to the ZEB2 promoter. The stable overexpression of FOXA2 abolished metastasis of breast cancer cells in vivo. This study confirmed that FOXA2 inhibited EMT in breast cancer cells by regulating the transcription of EMT-related genes such as E-cadherin and ZEB2.

A property-based analysis of human transcription factors

Background

Transcription factors are essential proteins for regulating gene expression. This regulation depends upon specific features of the transcription factors, including how they interact with DNA, how they interact with each other, and how they are post-translationally modified. Reliable information about key properties associated with transcription factors will therefore be useful for data analysis, in particular of data from high-throughput experiments.

Results

We have used an existing list of 1978 human proteins described as transcription factors to make a well-annotated data set, which includes information on Pfam domains, DNA-binding domains, post-translational modifications and protein–protein interactions. We have then used this data set for enrichment analysis. We have investigated correlations within this set of features, and between the features and more general protein properties. We have also used the data set to analyze previously published gene lists associated with cell differentiation, cancer, and tissue distribution.

Conclusions

The study shows that well-annotated feature list for transcription factors is a useful resource for extensive data analysis; both of transcription factor properties in general and of properties associated with specific processes. However, the study also shows that such analyses are easily biased by incomplete coverage in experimental data, and by how gene sets are defined.

Electronic supplementary material

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

Transcription factors and miRNAs that regulate fetal to adult CFTR expression change are new targets for cystic fibrosis

The CFTR gene displays a tightly regulated tissue-specific and temporal expression. Mutations in this gene cause cystic fibrosis (CF). In this study we wanted to identify trans-regulatory elements responsible for CFTR differential expression in fetal and adult lung, and to determine the importance of inhibitory motifs in the CFTR-3'UTR with the aim of developing new tools for the correction of disease-causing mutations within CFTR. We show that lung development-specific transcription factors (FOXA, C/EBP) and microRNAs (miR-101, miR-145, miR-384) regulate the switch from strong fetal to very low CFTR expression after birth. By using miRNome profiling and gene reporter assays, we found that miR-101 and miR-145 are specifically upregulated in adult lung and that miR-101 directly acts on its cognate site in the CFTR-3'UTR in combination with an overlapping AU-rich element. We then designed miRNA-binding blocker oligonucleotides (MBBOs) to prevent binding of several miRNAs to the CFTR-3'UTR and tested them in primary human nasal epithelial cells from healthy individuals and CF patients carrying the p.Phe508del CFTR mutation. These MBBOs rescued CFTR channel activity by increasing CFTR mRNA and protein levels. Our data offer new understanding of the control of the CFTR gene regulation and new putative correctors for cystic fibrosis.

Architectural proteins CTCF and cohesin have distinct roles in modulating the higher order structure and expression of the CFTR locus

Higher order chromatin structures across the genome are maintained in part by the architectural proteins CCCTC binding factor (CTCF) and the cohesin complex, which co-localize at many sites across the genome. Here, we examine the role of these proteins in mediating chromatin structure at the cystic fibrosis transmembrane conductance regulator (CFTR) gene. CFTR encompasses nearly 200 kb flanked by CTCF-binding enhancer-blocking insulator elements and is regulated by cell-type-specific intronic enhancers, which loop to the promoter in the active locus. SiRNA-mediated depletion of CTCF or the cohesin component, RAD21, showed that these two factors have distinct roles in regulating the higher order organization of CFTR. CTCF mediates the interactions between CTCF/cohesin binding sites, some of which have enhancer-blocking insulator activity. Cohesin shares this tethering role, but in addition stabilizes interactions between the promoter and cis-acting intronic elements including enhancers, which are also dependent on the forkhead box A1/A2 (FOXA1/A2) transcription factors (TFs). Disruption of the three-dimensional structure of the CFTR gene by depletion of CTCF or RAD21 increases gene expression, which is accompanied by alterations in histone modifications and TF occupancy across the locus, and causes internalization of the gene from the nuclear periphery.

A balance between activating and repressive histone modifications regulates cystic fibrosis transmembrane conductance regulator (CFTR) expression in vivo

The genetic mechanisms that regulate CFTR, the gene responsible for cystic fibrosis, have been widely investigated in cultured cells. However, mechanisms responsible for tissue-specific and time-specific expression are not completely elucidated in vivo. Through the survey of public databases, we found that the promoter of CFTR was associated with bivalent chromatin in human embryonic stem (ES) cells. In this work, we analyzed fetal (at different stages of pregnancy) and adult tissues and showed that, in digestive and lung tissues, which expressed CFTR, H3K4me3 was maintained in the promoter. Histone acetylation was high in the promoter and in two intronic enhancers, especially in fetal tissues. In contrast, in blood cells, which did not express CFTR, the bivalent chromatin was resolved (the promoter was labeled by the silencing mark H3K27me3). Cis-regulatory sequences were associated with lowly acetylated histones. We also provide evidence that the tissue-specific expression of CFTR is not regulated by dynamic changes of DNA methylation in the promoter. Overall, this work shows that a balance between activating and repressive histone modifications in the promoter and intronic enhancers results in the fine regulation of CFTR expression during development, thereby ensuring tissue specificity.