A promoter-dependent upstream activator augments CFTR expression in diverse epithelial cell types.

The cystic fibrosis transmembrane conductance regulator (CFTR) gene encodes an anion-selective channel found in epithelial cell membranes. Mutations in CFTR cause cystic fibrosis (CF), an inherited disorder that impairs epithelial function in multiple organs. Most men with CF are infertile due to loss of intact genital ducts. Here we investigated a novel epididymis-selective cis-regulatory element (CRE), located within a peak of open chromatin at -9.5 kb 5' to the CFTR gene promoter. Activation of the -9.5 kb CRE alone by CRISPRa had no impact on CFTR gene expression. However, CRISPRa co-activation of the -9.5 kb CRE and the CFTR gene promoter in epididymis cells significantly augmented CFTR mRNA and protein expression when compared to promoter activation alone. This increase was accompanied by enhanced chromatin accessibility at both sites. Furthermore, the combined CRISPRa strategy activated CFTR expression in other epithelial cells that lack open chromatin at the -9.5 kb site and in which the locus is normally inactive. However, the -9.5 kb CRE does not function as a classical enhancer of the CFTR promoter in transient reporter gene assays. These data provide a novel mechanism for activating/augmenting CFTR expression, which may have therapeutic utility for mutations that perturb CFTR transcription.

Functional genomics of the human epididymis: Further characterization of efferent ducts and model systems by single-cell RNA sequencing analysis.

BACKGROUND: The human epididymis is poorly studied due to the lack of availability of tissue samples. Our understanding of its structure and function depends upon anatomical and histological observations of archived material. OBJECTIVES: Here we used single-cell RNA sequencing (scRNA-seq) technologies to elucidate the identity of cells within the human efferent ducts (EDs) and compared them to caput epididymis cells. We also compared the cellularity of primary tissues with those of 2D and 3D (organoid) culture models used for functional studies. MATERIALS AND METHODS: Human epididymis tissue was dissected to separate different anatomical regions and digested to release single cells for processing on the 10X Genomics Chromium platform. Primary human epididymis epithelial (HEE) cells and HEE organoids were grown as described previously and subjected to scRNA-seq. scRNA-seq data were processed by standard bioinformatics pipelines and used for comparative analysis. RESULTS: We define the cell types in the EDs which include specialized epithelial cells, connective tissue stromal cells, vascular endothelial cells, smooth muscle cells, and immune cells, but lack basal cells that are seen in the caput epididymis. Furthermore, we identify a sub-population of epithelial cells which have marker genes found in the bladder and urothelium. Comparative genomics analysis of the 2D and 3D culture models shows cellular identities adapted to the culture environment while still maintaining similarity to the primary tissue. DISCUSSION: Our data suggest that EDs are lined with a transitional epithelium, which like the urothelium is able to stretch and contract depending on luminal volume. This is consistent with its primary role in seminal fluid resorption and sperm concentration. Moreover, we describe the cellularity of models to study the human epididymis epithelium in vitro. CONCLUSION: Single-cell RNA-seq data from the human epididymis make a valuable contribution to our understanding of this highly specialized organ.

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.

Early developmental phenotypes in the cystic fibrosis sheep model.

Highly effective modulator therapies for cystic fibrosis (CF) make it a treatable condition for many people. However, although CF respiratory illness occurs after birth, other organ systems particularly in the digestive tract are damaged before birth. We use an ovine model of CF to investigate the in utero origins of CF disease since the sheep closely mirrors critical aspects of human development. Wildtype (WT) and CFTR -/- sheep tissues were collected at 50, 65, 80, 100, and 120 days of gestation and term (147 days) and used for histological, electrophysiological, and molecular analysis. Histological abnormalities are evident in CFTR-/- -/-  animals by 80 days of gestation, equivalent to 21 weeks in humans. Acinar and ductal dilation, mucus obstruction, and fibrosis are observed in the pancreas; biliary fibrosis, cholestasis, and gallbladder hypoplasia in the liver; and intestinal meconium obstruction, as seen at birth in all large animal models of CF. Concurrently, cystic fibrosis transmembrane conductance regulator (CFTR)-dependent short circuit current is present in WT tracheal epithelium by 80 days gestation and is absent from CFTR -/- tissues. Transcriptomic profiles of tracheal tissues confirm the early expression of CFTR and suggest that its loss does not globally impair tracheal differentiation.

Enhancement of airway epithelial cell differentiation by pulmonary endothelial cell co-culture.

Cross-talk between lung epithelial cells and their microenvironment has an important physiological role in development. Using an in vitro model of differentiation of human induced pluripotent stem cells (iPSCs) to air-liquid interface (ALI)-cultured lung epithelial cells, we investigated the contribution of the microenvironment to maintenance of the lung progenitor cell state. Our protocol modeled in vivo cell-to-matrix and cell-to-cell interactions. These included growth of iPSCs on inserts coated with different basement membrane proteins (collagen I, IV, fibronectin, heparan sulfate or Matrigel plus collagen IV) and co-culture with human pulmonary microvascular endothelial cells (HPMECs). Marker gene expression was measured by RT-qPCR and protein expression and localization was confirmed by immunocytochemistry. The results showed that iPSCs grown on collagen IV had the highest success rate in terms of differentiation to robust ALI-cultured lung epithelial cells, followed by fibronectin, collagen I and heparan sulfate. Coating with Matrigel mixed with collagen IV further increased the success rate to > 97 %. Co-culture of iPSCs with HPMECs enhanced the expression of key airway lineage markers (NKX2.1, KRT5, TP63, MUC5AC, MUC16, FOXJ1, CFTR and SCGB1A1) during ALI culture. Cross-talk between iPSCs and their microenvironment during cell differentiation had a significant effect on lung epithelial cell differentiation in these 3D in vitro models. Both matrix proteins and endothelial cells play critical roles in the differentiation of lung progenitor cells.

Cell function and identity revealed by comparative scRNA-seq analysis in human nasal, bronchial and epididymis epithelia.

The evolutionary relationship of cells within tissues having a similar function but located in different anatomical sites is of considerable biological interest. The development of single-cell RNA sequencing (scRNA-seq) protocols has greatly enhanced opportunities to address this topic. Here we focus on cells in the epithelium which lines two regions of the human respiratory tract and the male genital ducts to delineate the shared, differentiated functions of the different cell populations. Transcriptomic data were used to assess the gene expression profiles of human bronchial, nasal, and epididymal epithelium (HBE, HNE, and HEE). Bulk RNA-seq showed many shared genes expressed in cells from the nasal and bronchial epithelium and highlighted their divergence from the epididymal epithelium. ScRNA-seq in HBE and HNE cells demonstrated overlapping gene expression patterns within basal and secretory cell populations. Moreover, the distribution of cell types was altered in HNE cells derived from donors with cystic fibrosis (CF) when compared to cells from healthy donors. Next, the HBE and HNE datasets were merged and confirmed intersection of cell type gene expression profiles from the two sites. However, secretory and ciliated cells were the most abundant types in the HBE samples, while more basal cells were seen in the HNE populations. We then merged single-cell data from the epididymis to determine if overlapping functions of these cells corresponded to those in the airway. Of note, only the pulmonary ionocytes/epididymis clear cells showed a strongly conserved identity, which was confirmed by imputation in bulk RNA-seq datasets from the same cells.

Sheep models of F508del and G542X cystic fibrosis mutations show cellular responses to human therapeutics.

Cystic Fibrosis (CF) is a genetic disease caused by mutations in the CF transmembrane conductance regulator (CFTR) gene. The F508del and G542X are the most common mutations found in US patients, accounting for 86.4% and 4.6% of all mutations, respectively. The F508del causes deletion of the phenylalanine residue at position 508 and is associated with impaired CFTR protein folding. The G542X is a nonsense mutation that introduces a stop codon into the mRNA, thus preventing normal CFTR protein synthesis. Here, we describe the generation of CFTRF508del / F508del and CFTRG542X / G542X lambs using CRISPR/Cas9 and somatic cell nuclear transfer (SCNT). First, we introduced either F508del or G542X mutations into sheep fetal fibroblasts that were subsequently used as nuclear donors for SCNT. The newborn CF lambs develop pathology similar to CFTR -/- sheep and CF patients. Moreover, tracheal epithelial cells from the CFTRF508del / F508del lambs responded to a human CFTR (hCFTR) potentiator and correctors, and those from CFTRG542X / G542X lambs showed modest restoration of CFTR function following inhibition of nonsense-mediated decay (NMD) and aminoglycoside antibiotic treatments. Thus, the phenotype and electrophysiology of these novel models represent an important advance for testing new CF therapeutics and gene therapy to improve the health of patients with this life-limiting disorder.

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.

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.

The Bromodomain Containing 8 (BRD8) transcriptional network in human lung epithelial cells.

Mechanisms regulating gene expression in the airway epithelium underlie its response to the environment. A network of transcription factors (TFs) and architectural proteins, modulate chromatin accessibility and recruit activating or repressive signals. Bromodomain-containing proteins function as TFs or by engaging methyltransferase or acetyltransferase activity to induce chromatin modifications. Here we investigate the role of Bromodomain Containing 8 (BRD8) in coordinating lung epithelial function. Sites of BRD8 occupancy genome-wide were mapped in human lung epithelial cell lines (Calu-3 and 16HBE14o-). CCCTC-Binding Factor (CTCF) was identified as a predicted co-factor of BRD8, based upon motif over-representation under BRD8 ChIP-seq peaks. Following siRNA-mediated depletion of BRD8, differentially expressed genes with nearby peaks of BRD8 occupancy were subject to gene ontology process enrichment analysis. BRD8 targets are enriched for genes involved in the innate immune response and the cell cycle. Depletion of BRD8 increased the secretion of the antimicrobial peptide beta-defensin 1 and multiple chemokines, and reduced cell proliferation.

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.

An atlas of human proximal epididymis reveals cell-specific functions and distinct roles for CFTR.

Spermatozoa released from the testis are unable to fertilize an egg without a coordinated process of maturation in the lumen of the epididymis. Relatively little is known about the molecular events that integrate this critical progression along the male genital ducts in man. Here, we use single cell RNA-sequencing to construct an atlas of the human proximal epididymis. We find that the CFTR, which is pivotal in normal epididymis fluid transport, is most abundant in surface epithelial cells in the efferent ducts and in rare clear cells in the caput epididymis, suggesting region-specific functional properties. We reveal transcriptional signatures for multiple cell clusters, which identify the individual roles of principal, apical, narrow, basal, clear, halo, and stromal cells in the epididymis. A marked cell type-specific distribution of function is seen along the duct with local specialization of individual cell types integrating processes of sperm maturation.

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.

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.

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 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.

Region-specific microRNA signatures in the human epididymis.

The epithelium of the human epididymis maintains an appropriate luminal environment for sperm maturation that is essential for male fertility. Regional expression of small noncoding RNAs such as microRNAs contributes to segment-specific gene expression and differentiated functions. MicroRNA profiles were reported in human epididymal tissues but not specifically in the epithelial cells derived from those regions. Here, we reveal miRNA signatures of primary cultures of caput, corpus, and cauda epididymis epithelial cells and of the tissues from which they were derived. We identify 324 epithelial cell-derived microRNAs and 259 tissue-derived microRNAs in the epididymis, some of which displayed regionalized expression patterns in cells and/or tissues. Caput cell-enriched miRNAs included miR-573 and miR-155. Cauda cell-enriched miRNAs included miR-1204 and miR-770. Next, we determined the gene ontology pathways associated with in silico predicted target genes of the differentially expressed miRNAs. The effect of androgen receptor stimulation on miRNA expression was also investigated. These data show novel epithelial cell-derived miRNAs that may regulate the expression of important gene networks that are responsible for the regionalized gene expression and function of the epididymis.