Stellate cells are in utero markers of pancreatic disease in cystic fibrosis.

BACKGROUND: Pancreatic fibrosis is an early diagnostic feature of the common inherited disorder cystic fibrosis (CF). Many people with CF (pwCF) are pancreatic insufficient from birth and the replacement of acinar tissue with cystic lesions and fibrosis is a progressive phenotype that may later lead to diabetes. Little is known about the initiating events in the fibrotic process though it may be a sequela of inflammation in the pancreatic ducts resulting from loss of CFTR impairing normal fluid secretion. Here we use a sheep model of CF (CFTR-/-) to examine the evolution of pancreatic disease through gestation. METHODS: Fetal pancreas was collected at six time points from 50-days of gestation through to term, which is equivalent to ~ 13 weeks to term in human. RNA was extracted from tissue for bulk RNA-seq and single cells were prepared from 80-day, 120-day and term samples for scRNA-seq. Data were validated by immunochemistry. RESULTS: Transcriptomic evidence from bulk RNA-seq showed alterations in the CFTR-/- pancreas by 65-days of gestation, which are accompanied by marked pathological changes by 80-days of gestation. These include a fibrotic response, confirmed by immunostaining for COL1A1, αSMA and SPARC, together with acinar loss. Moreover, using scRNA-seq we identify a unique cell population that is significantly overrepresented in the CFTR-/- animals at 80- and 120-days gestation, as are stellate cells at term. CONCLUSION: The transcriptomic changes and cellular imbalance that we observe likely have pivotal roles in the evolution of CF pancreatic disease and may provide therapeutic opportunities to delay or prevent pancreatic destruction in CF.

Transcriptomic analysis of lung development in wildtype and CFTR-/- sheep suggests an early inflammatory signature in the CF distal lung.

The precise molecular events initiating human lung disease are often poorly characterized. Investigating prenatal events that may underlie lung disease in later life is challenging in man, but insights from the well-characterized sheep model of lung development are valuable. Here, we determine the transcriptomic signature of lung development in wild-type sheep (WT) and use a sheep model of cystic fibrosis (CF) to characterize disease associated changes in gene expression through the pseudoglandular, canalicular, saccular, and alveolar stages of lung growth and differentiation. Using gene ontology process enrichment analysis of differentially expressed genes at each developmental time point, we define changes in biological processes (BP) in proximal and distal lung from WT or CF animals. We also compare divergent BP in WT and CF animals at each time point. Next, we establish the developmental profile of key genes encoding components of ion transport and innate immunity that are pivotal in CF lung disease and validate transcriptomic data by RT-qPCR. Consistent with the known pro-inflammatory phenotype of the CF lung after birth, we observe upregulation of inflammatory response processes in the CF sheep distal lung during the saccular stage of prenatal development. These data suggest early commencement of therapeutic regimens may be beneficial.

Krüppel-like factor 5 regulates wound repair and the innate immune response in human airway epithelial cells.

A complex network of transcription factors regulates genes involved in establishing and maintaining key biological properties of the human airway epithelium. However, detailed knowledge of the contributing factors is incomplete. Here we characterize the role of Krüppel-like factor 5 (KLF5), in controlling essential pathways of epithelial cell identity and function in the human lung. RNA-seq following siRNA-mediated depletion of KLF5 in the Calu-3 lung epithelial cell line identified significant enrichment of genes encoding chemokines and cytokines involved in the proinflammatory response and also components of the junctional complexes mediating cell adhesion. To determine direct gene targets of KLF5, we defined the cistrome of KLF5 using ChIP-seq in both Calu-3 and 16HBE14o- lung epithelial cell lines. Occupancy site concordance analysis revealed that KLF5 colocalized with the active histone modification H3K27ac and also with binding sites for the transcription factor CCAAT enhancer-binding protein beta (C/EBPß). Depletion of KLF5 increased both the expression and secretion of cytokines including IL-1ß, a response that was enhanced following exposure to Pseudomonas aeruginosa lipopolysaccharide. Calu-3 cells exhibited faster rates of repair after KLF5 depletion compared with negative controls in wound scratch assays. Similarly, CRISPR-mediated KLF5-null 16HBE14o- cells also showed enhanced wound closure. These data reveal a pivotal role for KLF5 in coordinating epithelial functions relevant to human lung disease.

BACH1, the master regulator of oxidative stress, has a dual effect on CFTR expression.

The cystic fibrosis transmembrane conductance regulator (CFTR) gene lies within a topologically associated domain (TAD) in which multiple cis-regulatory elements (CREs) and transcription factors (TFs) regulate its cell-specific expression. The CREs are recruited to the gene promoter by a looping mechanism that depends upon both architectural proteins and specific TFs. An siRNA screen to identify TFs coordinating CFTR expression in airway epithelial cells suggested an activating role for BTB domain and CNC homolog 1 (BACH1). BACH1 is a ubiquitous master regulator of the cellular response to oxidative stress. Here, we show that BACH1 may have a dual effect on CFTR expression by direct occupancy of CREs at physiological oxygen (∼8%), while indirectly modulating expression under conditions of oxidative stress. Hence BACH1, can activate or repress the same gene, to fine tune expression in response to environmental cues such as cell stress. Furthermore, our 4C-seq data suggest that BACH1 can also directly regulate CFTR gene expression by modulating locus architecture through occupancy at known enhancers and structural elements, and depletion of BACH1 alters the higher order chromatin structure.

Looping of upstream cis-regulatory elements is required for CFTR expression in human airway epithelial cells.

The CFTR gene lies within an invariant topologically associated domain (TAD) demarcated by CTCF and cohesin, but shows cell-type specific control mechanisms utilizing different cis-regulatory elements (CRE) within the TAD. Within the respiratory epithelium, more than one cell type expresses CFTR and the molecular mechanisms controlling its transcription are likely divergent between them. Here, we determine how two extragenic CREs that are prominent in epithelial cells in the lung, regulate expression of the gene. We showed earlier that these CREs, located at -44 and -35 kb upstream of the promoter, have strong cell-type-selective enhancer function. They are also responsive to inflammatory mediators and to oxidative stress, consistent with a key role in CF lung disease. Here, we use CRISPR/Cas9 technology to remove these CREs from the endogenous locus in human bronchial epithelial cells. Loss of either site extinguished CFTR expression and abolished long-range interactions between these sites and the gene promoter, suggesting non-redundant enhancers. The deletions also greatly reduced promoter interactions with the 5' TAD boundary. We show substantial recruitment of RNAPII to the -35 kb element and identify CEBPß as a key activator of airway expression of CFTR, likely through occupancy at this CRE and the gene promoter.

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.

A functional genomics approach to investigate the differentiation of iPSCs into lung epithelium at air-liquid interface.

The availability of robust protocols to differentiate induced pluripotent stem cells (iPSCs) into many human cell lineages has transformed research into the origins of human disease. The efficacy of differentiating iPSCs into specific cellular models is influenced by many factors including both intrinsic and extrinsic features. Among the most challenging models is the generation of human bronchial epithelium at air-liquid interface (HBE-ALI), which is the gold standard for many studies of respiratory diseases including cystic fibrosis. Here, we perform open chromatin mapping by ATAC-seq and transcriptomics by RNA-seq in parallel, to define the functional genomics of key stages of the iPSC to HBE-ALI differentiation. Within open chromatin peaks, the overrepresented motifs include the architectural protein CTCF at all stages, while motifs for the FOXA pioneer and GATA factor families are seen more often at early stages, and those regulating key airway epithelial functions, such as EHF, are limited to later stages. The RNA-seq data illustrate dynamic pathways during the iPSC to HBE-ALI differentiation, and also the marked functional divergence of different iPSC lines at the ALI stages of differentiation. Moreover, a comparison of iPSC-derived and lung donor-derived HBE-ALI cultures reveals substantial differences between these models.

Functional genomics analysis of human colon organoids identifies key transcription factors.

Organoids are a valuable three-dimensional (3D) model to study the differentiated functions of the human intestinal epithelium. They are a particularly powerful tool to measure epithelial transport processes in health and disease. Though biological assays such as organoid swelling and intraluminal pH measurements are well established, their underlying functional genomics are not well characterized. Here we combine genome-wide analysis of open chromatin by ATAC-Seq with transcriptome mapping by RNA-Seq to define the genomic signature of human intestinal organoids (HIOs). These data provide an important tool for investigating key physiological and biochemical processes in the intestinal epithelium. We next compared the transcriptome and open chromatin profiles of HIOs with equivalent data sets from the Caco2 colorectal carcinoma line, which is an important two-dimensional (2D) model of the intestinal epithelium. Our results define common features of the intestinal epithelium in HIO and Caco2 and further illustrate the cancer-associated program of the cell line. Generation of Caco2 cysts enabled interrogation of the molecular divergence of the 2D and 3D cultures. Overrepresented motif analysis of open chromatin peaks identified caudal type homeobox 2 (CDX2) as a key activating transcription factor in HIO, but not in monolayer cultures of Caco2. However, the CDX2 motif becomes overrepresented in open chromatin from Caco2 cysts, reinforcing the importance of this factor in intestinal epithelial differentiation and function. Intersection of the HIO and Caco2 transcriptomes further showed functional overlap in pathways of ion transport and tight junction integrity, among others. These data contribute to understanding human intestinal organoid biology.

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.

An organoid model to assay the role of CFTR in the human epididymis epithelium.

Organoid cultures derived from primary human tissues facilitate the study of disease processes and the development of new therapeutics. Most men with cystic fibrosis (CF) are infertile due to defects in the epididymis and vas deferens; however, the causative mechanisms are still unclear. We used human epididymis epithelial cell (HEE) organoids and polarized HEE cell cultures to assay the CF transmembrane conductance regulator (CFTR) in the human epididymis. 3D HEE organoids and polarized 2D HEE cell cultures on membrane inserts were established from human caput epididymis. Single-cell RNA sequencing (scRNA-seq) was performed to map cell type-specific gene expression in the organoids. Using forskolin (FSK) to activate CFTR and inhibitor CFTRinh172 to block its activity, we assessed how CFTR contributes to organoid swelling and epithelial barrier function. The scRNA-seq data showed key caput epididymis cell types present in HEE organoid cultures. FSK at 10 µM induced HEE organoid swelling by 20% at 16 h, while 5 and 10 µM CFTRinh172 treatment significantly reduced HEE organoid size. In transepithelial resistance (TER) measurements, FSK reduced TER, while inhibition of CFTR increased TER; also, depletion of CFTR with specific siRNAs significantly increased TER. FSK treatment significantly increased the flux of 4-kDa but not 70-kDa dextran, suggesting activation of CFTR mainly enhances transcellular diffusion. We have demonstrated that CFTR contributes to the maintenance of HEE cell TER and that cultured HEE organoids are a useful model to investigate human epididymis function. These results facilitate progress in elucidating how CFTR-dependent cellular processes impair fertility in CF.

Coordinate regulation of ELF5 and EHF at the chr11p13 CF modifier region.

E74-like factor 5 (ELF5) and ETS-homologous factor (EHF) are epithelial selective ETS family transcription factors (TFs) encoded by genes at chr11p13, a region associated with cystic fibrosis (CF) lung disease severity. EHF controls many key processes in lung epithelial function so its regulatory mechanisms are important. Using CRISPR/Cas9 technology, we removed three key cis-regulatory elements (CREs) from the chr11p13 region and also activated multiple open chromatin sites with CRISPRa in airway epithelial cells. Deletion of the CREs caused subtle changes in chromatin architecture and site-specific increases in EHF and ELF5. CRISPRa had most effect on ELF5 transcription. ELF5 levels are low in airway cells but higher in LNCaP (prostate) and T47D (breast) cancer cells. ATAC-seq in these lines revealed novel peaks of open chromatin at the 5' end of chr11p13 associated with an expressed ELF5 gene. Furthermore, 4C-seq assays identified direct interactions between the active ELF5 promoter and sites within the EHF locus, suggesting coordinate regulation between these TFs. ChIP-seq for ELF5 in T47D cells revealed ELF5 occupancy within EHF introns 1 and 6, and siRNA-mediated depletion of ELF5 enhanced EHF expression. These results define a new role for ELF5 in lung epithelial biology.

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.

Regulatory dynamics of 11p13 suggest a role for EHF in modifying CF lung disease severity.

Mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene cause cystic fibrosis (CF), but are not good predictors of lung phenotype. Genome-wide association studies (GWAS) previously identified additional genomic sites associated with CF lung disease severity. One of these, at chromosome 11p13, is an intergenic region between Ets homologous factor (EHF) and Apaf-1 interacting protein (APIP). Our goal was to determine the functional significance of this region, which being intergenic is probably regulatory. To identify cis-acting elements, we used DNase-seq and H3K4me1 and H3K27Ac ChIP-seq to map open and active chromatin respectively, in lung epithelial cells. Two elements showed strong enhancer activity for the promoters of EHF and the 5' adjacent gene E47 like ETS transcription factor 5 (ELF5) in reporter gene assays. No enhancers of the APIP promoter were found. Circular chromosome conformation capture (4C-seq) identified direct physical interactions of elements within 11p13. This confirmed the enhancer-promoter associations, identified additional interacting elements and defined topologically associating domain (TAD) boundaries, enriched for CCCTC-binding factor (CTCF). No strong interactions were observed with the APIP promoter, which lies outside the main TAD encompassing the GWAS signal. These results focus attention on the role of EHF in modifying CF lung disease severity.

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.

A novel transcriptional network for the androgen receptor in human epididymis epithelial cells.

STUDY QUESTION: What is the transcriptional network governed by the androgen receptor (AR) in human epididymis epithelial (HEE) cells from the caput region and if the network is tissue-specific, how is this achieved? SUMMARY ANSWER: About 200 genes are differentially expressed in the caput HEE cells after AR activation; the AR transcriptional network is tissue-specific and may be mediated in part by distinct AR co-factors including CAAT-enhancer binding protein beta (CEBPB) and runt-related transcription factor 1 (RUNX1). WHAT IS KNOWN ALREADY: Little is known about the AR transcriptional program genome wide in HEE cells, nor its co-factors in those cells. AR has been best studied in the prostate gland epithelium and prostate cancer cell lines, due to the important role of this factor in prostate cancer. However AR-associated differentially expressed genes (DEGs) and AR co-factors have not yet been compared between human epididymis and prostate epithelial cells. STUDY DESIGN, SIZE, DURATION: Caput HEE cells from two donors were exposed to the synthetic androgen R1881 at 1 nM for 12-16 h after 72 h of hormone starvation. PARTICIPANTS/MATERIALS, SETTING, METHODS: Chromatin was prepared from R1881-treated and vehicle control HEE cells. AR-associated chromatin was purified by chromatin immunoprecipitation (ChIP) and AR occupancy genome wide was revealed by deep sequencing (ChIP-seq). Two independent biological replicates were performed. Total RNA was prepared from R1881 and control-treated HEE cells and gene expression profiles were documented by RNA-seq. The interaction of the potential novel AR co-factors CEBPB and RUNX1, identified through in-silico motif analysis of AR ChIP-seq data, was examined by ChIP-qPCR after siRNA-mediated depletion of each co-factor individually or simultaneously. MAIN RESULTS AND THE ROLE OF CHANCE: The results identify about 200 genes that are differentially expressed (DEGs) in HEE cells after AR activation. Some of these DEGs show occupancy of AR at their promoters or cis-regulatory elements suggesting direct regulation. However, there is little overlap in AR-associated DEGs between HEE and prostate epithelial cells. Inspection of over-represented motifs in AR ChIP-seq peaks identified CEBPB and RUNX1 as potential co-factors, with no evidence for FOXA1, which is an important co-factor in the prostate epithelium. CEBPB and RUNX1 ChIP-seq in HEE cells showed that both these factors often occupied AR-binding sites, though rarely simultaneously. Further analysis at a single AR-regulated locus (FK506-binding protein 5, FKPB5) suggests that CEBPB may be a co-activator. These data suggest a novel AR transcriptional network governs differentiated functions of the human epididymis epithelium. LARGE SCALE DATA: AR ChIP-seq and RNA-seq data are deposited at GEO: GSE109063. LIMITATIONS, REASONS FOR CAUTION: There is substantial donor-to-donor variation in primary HEE cells cultures. We applied stringent statistical tests with a false discovery rate (FDR) of 0.1% for ChIP-seq and standard pipelines for RNA-seq so it is possible that we have missed some AR-regulated genes that are important in caput epididymis function. WIDER IMPLICATIONS OF THE FINDINGS: Our data suggest that a novel AR transcriptional network governs differentiated functions of the human epididymis epithelium. Since this cell layer has a critical role in normal sperm maturation, the results are of broader significance in understanding the mechanisms underlying the maintenance of fertility in men. STUDY FUNDING/COMPETING INTERESTS: This work was funded by the National Institutes of Health, Eunice Kennedy Shriver National Institute of Child Health and Development: R01 HD068901 (PI: Harris). The authors have no competing interests to declare.

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.

Expression profiles of human epididymis epithelial cells reveal the functional diversity of caput, corpus and cauda regions.

STUDY HYPOTHESIS: Region-specific transcriptional profiling of tissues and cultured epithelial cells from the human epididymis will predict functional specialization along the duct. STUDY FINDING: We identified the molecular signature driving functions of the caput, corpus and cauda epithelium, and determined how these differ to establish the regional differentiation of the organ. WHAT IS KNOWN ALREADY: The epithelium lining the human male genital ducts has a critical role in fertility. In particular, it controls the luminal environment in the epididymis, which is required for normal sperm maturation and reproductive competence. Studies in many animal species have largely informed our understanding of the molecular basis of epididymis function. However, there are substantial differences between species. STUDY DESIGN, SAMPLES/MATERIALS, METHODS: Using RNA sequencing on biological replicates, we described gene expression profiles for tissue from each region of the epididymis and cultured epithelial cells derived from these regions. Bioinformatic tools were then utilized to identify differentially expressed genes (DEGs) between tissues and cells from the caput, corpus and cauda. MAIN RESULTS AND THE ROLE OF CHANCE: The data showed that the caput is functionally divergent from the corpus and cauda, which have very similar transcriptomes. Interrogation of DEGs using gene ontology process enrichment analyses showed that processes of ion transport, response to hormone stimulus and urogenital tract development are more evident in the caput, while defense response processes are more important in the corpus/cauda. Consistent with these regional differences in epididymis function, we observed differential expression of transcription factors in the caput and corpus/cauda. LIMITATIONS, REASONS FOR CAUTION: Cultured caput, corpus and cauda cells may not faithfully represent the same cells in the intact organ, due to loss of hormonal signals from the testis and communication from other cell types. WIDER IMPLICATIONS OF THE FINDINGS: Our data provide a molecular characterization that will facilitate advances in understanding human epididymis epithelium biology in health and disease. They may also reveal the mechanisms coordinating epididymis luminal environment and sperm maturation. LARGE SCALE DATA: Data deposited at http://www.ncbi.nlm.nih.gov/geo/GSE72986. STUDY FUNDING AND COMPETING INTERESTS: This work was supported by the National Institutes of Health: R01HD068901 (PI: A.H.). The authors declare no conflict of interest.

Ets homologous factor regulates pathways controlling response to injury in airway epithelial cells.

Ets homologous factor (EHF) is an Ets family transcription factor expressed in many epithelial cell types including those lining the respiratory system. Disruption of the airway epithelium is central to many lung diseases, and a network of transcription factors coordinates its normal function. EHF can act as a transcriptional activator or a repressor, though its targets in lung epithelial cells are largely uncharacterized. Chromatin immunoprecipitation followed by deep sequencing (ChIP-seq), showed that the majority of EHF binding sites in lung epithelial cells are intergenic or intronic and coincide with putative enhancers, marked by specific histone modifications. EHF occupies many genomic sites that are close to genes involved in intercellular and cell-matrix adhesion. RNA-seq after EHF depletion or overexpression showed significant alterations in the expression of genes involved in response to wounding. EHF knockdown also targeted genes in pathways of epithelial development and differentiation and locomotory behavior. These changes in gene expression coincided with alterations in cellular phenotype including slowed wound closure and increased transepithelial resistance. Our data suggest that EHF regulates gene pathways critical for epithelial response to injury, including those involved in maintenance of barrier function, inflammation and efficient wound repair.

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.

Nucleosome mapping across the CFTR locus identifies novel regulatory factors.

Nucleosome positioning on the chromatin strand plays a critical role in regulating accessibility of DNA to transcription factors and chromatin modifying enzymes. Hence, detailed information on nucleosome depletion or movement at cis-acting regulatory elements has the potential to identify predicted binding sites for trans-acting factors. Using a novel method based on enrichment of mononucleosomal DNA by bacterial artificial chromosome hybridization, we mapped nucleosome positions by deep sequencing across 250 kb, encompassing the cystic fibrosis transmembrane conductance regulator (CFTR) gene. CFTR shows tight tissue-specific regulation of expression, which is largely determined by cis-regulatory elements that lie outside the gene promoter. Although multiple elements are known, the repertoire of transcription factors that interact with these sites to activate or repress CFTR expression remains incomplete. Here, we show that specific nucleosome depletion corresponds to well-characterized binding sites for known trans-acting factors, including hepatocyte nuclear factor 1, Forkhead box A1 and CCCTC-binding factor. Moreover, the cell-type selective nucleosome positioning is effective in predicting binding sites for novel interacting factors, such as BAF155. Finally, we identify transcription factor binding sites that are overrepresented in regions where nucleosomes are depleted in a cell-specific manner. This approach recognizes the glucocorticoid receptor as a novel trans-acting factor that regulates CFTR expression in vivo.