HNF1 regulates critical processes in the human epididymis epithelium

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.

HNF1A induces glioblastoma by upregulating EPS8 and activating PI3K/AKT signaling pathway

Despite the exact biological role of HNF1 homolog A (HNF1A) in the regulatory mechanism of glioblastoma (GBM), the molecular mechanism, especially the downstream regulation as a transcription factor, remains to be further elucidated. Immunohistochemistry was used to detect the expression and clinical relevance of HNF1A in GBM patients. CCK8, TUNEL, and subcutaneous tumor formation in nude mice were used to evaluate the effect of HNF1A on GBM in vitro and in vivo. The correction between HNF1A and epidermal growth factor receptor pathway substrate 8 (EPS8) was illustrated by bioinformatics analysis and luciferase assay. Further mechanism was explored that the transcription factor HNF1A regulated the expression of EPS8 and downstream signaling pathways by directly binding to the promoter region of EPS8. Our comprehensive analysis of clinical samples in this study showed that upregulated expression of HNF1A was associated with poor survival in GBM patients. Further, we found that knockdown of HNF1A markedly suppressed the malignant phenotype of GBM cells in vivo and in vitro as well as promoted apoptosis of tumor cells, which was reversed by upregulation of HNF1A. Mechanistically, HNF1A could significantly activate PI3K/AKT signaling pathway by specifically binding to the promoter regions of EPS8. Moreover, overexpression of EPS8 was able to reverse the apoptosis of tumor cells caused by HNF1A knockdown, thereby exacerbating the GBM progression. Correctively, our study has clarified the explicit mechanism by which HNF1A promotes GBM malignancy and provides a new therapeutic target for further clinical application.

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.

Dual SMAD inhibition enhances the longevity of human epididymis epithelial cells

Primary human epididymis epithelial (HEE) cells are valuable reagents for functional studies on the human epididymis. We used them previously to determine the transcriptional networks that establish cell identity along the length of the epididymis from caput, corpus, and cauda. These studies on HEE cells and organoids derived from them revealed important cellular properties. However, similar to other primary cells, HEE cells undergo replicative senescence and de-differentiation in culture. A cocktail of small molecules was shown elsewhere to extend longevity of epithelial basal cells. The components included transforming growth factor β (TGF-β)/bone morphogenetic protein (BMP) antagonists, WNT agonist, and Rho-associated and coiled-coil containing protein kinase (ROCK) inhibitor (ROCKi), which together prevented the senescence-related upregulation of TGF-β signaling pathway members. Here, we treat HEE cells with the same cocktail and observed enhanced replicative potential and prolonged expression of markers of HEE differentiation. This treatment expands the differentiated HEE cell population available from individual epididymis tissue samples that can be used for molecular, cellular, and functional studies.

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.

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.

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.

Transcriptome analysis of turkey (Meleagris gallopavo) reproductive tract revealed key pathways regulating spermatogenesis and post-testicular sperm maturation

The application of transcriptomics to the study of the reproductive tract in male turkeys can significantly increase our current knowledge regarding the specifics of bird reproduction. To characterize the complex transcriptomic changes that occur in the testis, epididymis, and ductus deferens, deep sequencing of male turkey RNA samples (n = 6) was performed, using Illumina RNA-Seq. The obtained sequence reads were mapped to the turkey genome, and relative expression values were calculated to analyze differentially expressed genes (DEGs). Statistical analysis revealed 1,682; 2,150; and 340 DEGs in testis/epididymis, testis/ductus deferens, and epididymis/ductus deferens comparisons, respectively. The expression of selected genes was validated using quantitative real-time reverse transcriptase-polymerase chain reaction. Bioinformatics analysis revealed several potential candidate genes involved in spermatogenesis, spermiogenesis and flagellum formation in the testis, and in post-testicular sperm maturation in the epididymis and ductus deferens. In the testis, genes were linked with the mitotic proliferation of spermatogonia and the meiotic division of spermatocytes. Histone ubiquitination and protamine phosphorylation were shown to be regulatory mechanisms for nuclear condensation during spermiogenesis. The characterization of testicular transcripts allowed a better understanding of acrosome formation and development and flagellum formation, including axoneme structures and functions. Spermatozoa motility during post-testicular maturation was linked to the development of flagellar actin filaments and biochemical processes, including Ca2+ influx and protein phosphorylation/dephosphorylation. Spermatozoa quality appeared to be controlled by apoptosis and antioxidant systems in the epididymis and ductus deferens. Finally, genes associated with reproductive system development and morphogenesis were identified. To the best of our knowledge, this is the first genome-wide functional investigation of genes associated with tissue-specific processes in turkey reproductive tract. A catalog of genes worthy of further studies to understand the avian reproductive physiology and regulation was provided.

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.

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.

siRNA-mediated knockdown of sperm-associated antigen 11a (Spag11a) mRNA in epididymal primary epithelial cells affects proliferation: a transcriptome analyses

Differential expression of a variety of proteins in the four major regions of the epididymis contributes to maturation of spermatozoa and region-specific cellular functions as well. Proliferation of epithelial cells of the epididymis is highly controlled and thus is one of the major reasons for the nonoccurrence of cancers in this organ system. The molecular mechanisms and the contribution of region-specific genes in epithelial cell proliferation are not yet fully understood. In this study, for the first time, we analyzed the role of sperm-associated antigen 11a (Spag11a), a caput-specific beta-defensin-like antimicrobial gene in governing epididymal cell proliferation and global gene expression. siRNA-mediated knockdown of Spag11a mRNA in epididymal primary epithelial cells resulted in increased cell proliferation. Out of the 68,842 genes analyzed, 4182 genes were differentially expressed (2154 upregulated and 2028 downregulated). A variety of genes that participate in different cellular processes and pathways were differentially regulated. Genes that are important for epithelial cell proliferation were found to be differentially regulated and these changes were confirmed by real-time PCR. Overexpression of Spag11a in immortalized rat caput epididymal cells resulted in decreased proliferation capacity. Results of this study indicate that Spag11a plays a crucial role in governing epididymal epithelial cell proliferation.

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.

Transcriptional networks in the human epididymis

BACKGROUND

The epithelial lining of the human epididymis is critical for sperm maturation. This process requires distinct specialized functions in the head, body, and tail of the duct. These region-specific properties are maintained by distinct gene expression profiles which are governed by transcription factor networks, non-coding RNAs, and other factors.

MATERIALS AND METHODS

We used genome-wide protocols including DNase-seq, RNA-seq and ChIP-seq to characterize open (active) chromatin, the transcriptome and occupancy of specific transcription factors (TFs) respectively, in caput, corpus, and cauda segments of adult human epididymis tissue and primary human epididymis epithelial (HEE) cell cultures derived from them. RNA-seq following TF depletion or activation, combined with gene ontology analysis also determined TF targets.

RESULTS

Among regional differentially expressed transcripts were epithelial-selective transcription factors (TFs), microRNAs, and antiviral response genes. Caput-enriched TFs included hepatocyte nuclear factor 1 (HNF1) and the androgen receptor (AR), both of which were also predicted to occupy cis-regulatory elements identified as open chromatin in HEE cells. HNF1 targets were identified genome-wide using ChIP-seq, in HEE cells. Next, siRNA-mediated depletion of HNF1 revealed a pivotal role for this TF in coordinating epithelial water and solute transport in caput epithelium. The importance of AR in HEE cells was shown by AR ChIP-seq, and by RNA-seq after synthetic androgen (R1881) treatment. AR has a distinct transcriptional program in the HEE cells and likely recruits different co-factors (RUNX1 and CEBPβ) in comparison to those used in prostate epithelium.

DISCUSSION AND CONCLUSION

Our data identify many transcription factors that regulate the development and differentiation of HEE cells. Moreover, a comparison between immature and adult HEE cells showed key TFs in the transition to fully differentiated function of this epithelium. These data may help identify new targets to treat male infertility and have the potential to open new avenues for male contraception.

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.