Transcriptome analysis to identify the downstream genes of androgen receptor in dermal papilla cells
BMC Genomic Data volume 23, Article number: 2 (2022)
Testosterone signaling mediates various diseases, such as androgenetic alopecia and prostate cancer. Testosterone signaling is mediated by the androgen receptor (AR). In this study, we fortuitously found that primary and immortalized dermal papilla cells suppressed AR expression, although dermal papilla cells express AR in vivo. To analyze the AR signaling pathway, we exogenously introduced the AR gene via a retrovirus into immortalized dermal papilla cells and comprehensively compared their expression profiles with and without AR expression.
Whole-transcriptome profiling revealed that the focal adhesion pathway was mainly affected by the activation of AR signaling. In particular, we found that caveolin-1 gene expression was downregulated in AR-expressing cells, suggesting that caveolin-1 is controlled by AR.
Our whole transcriptome data is critical resources for discovery of new therapeutic targets for testosterone-related diseases.
The comprehensive gene expression profiling were obtained by RNA-Seq analysis about AR negative and AR positive dermal papilla cells.
The bioinformatics analysis suggested that caveolin-1 and EGF receptors are the downstream of AR signaling.
Our study showed the combination of pinpoint mutant cells and global transcriptome is effective to identify the downstream genes.
The software and its versions were used for the data analysis.
FastQC, version 0.11.8 was used for quality check of raw FASTQ sequencing file. https://www.bioinformatics.babraham.ac.uk/projects/fastqc/
PRINSEQ, version 0.20.4 was used for the removal of low quality reads. http://prinseq.sourceforge.net/
PEAT, version 1.2 was used for the removal of the adaptor sequence. https://github.com/jhhung/PEAT
STAR, version 2.6.1 was used for the mapping. https://github.com/alexdobin/STAR
featureCount, SUBREAD, release 1.6.5 was used for the expression counting. http://subread.sourceforge.net
R package, version 4.0.3, was used for the downstream analysis. https://www.r-project.org
TCC-GUI, tool for the downstream analysis. https://github.com/swsoyee/TCC-GUI
Testosterone is a hormone that controls cell growth or sexual differentiation of the reproductive organs. This hormone is known to be associated with various diseases, such as androgenetic alopecia (AGA) and prostate cancer. Approximately 30% of males are reported to be affected by AGA . The mechanism of AGA progression is the activation of AR signaling and secretion of Dickkopf-related protein 1 (DKK1) or tumor growth factor-β2 (TGF-β2), which suppress the growth of hair matrix cells . Testosterone signaling is mainly mediated by the androgen receptor (AR) . The activation of AR signaling starts with the binding of ligands, such as testosterone or dihydroxy-testosterone, which has a high affinity for AR. Ligand-activated AR forms a dimer on the cellular membrane and translocates from the cytoplasm to the nucleus. After nuclear translocation, AR forms a transcriptional complex, including coactivators and RNA polymerase, resulting in the transcriptional activation of downstream genes . In prostate cancer, the growth of cancer cells strongly depends on the activation of testosterone . Therefore, the identification of a chemical that can inhibit testosterone signaling may provide a strong candidate for the treatment of prostate cancer and/or AGA.
The hair growth is controlled by the growth signal from human follicle dermal papilla cells (HFDPCs). Interestingly, Kwack et al. showed that primary HFDPCs cause a decrease in AR expression even after passage 6. In agreement with this phenomenon, we could not measure the expression level of AR, even in early passages of primary HFDPCs . Although the detailed mechanism is not clear, the expression level of AR dramatically decreases after sequential passages of the DPCs. The immortalized cells established from corresponding primary cells were also negative for AR expression . These situations led us to hypothesize that we could reconstitute the AR signaling pathway if AR is exogenously introduced. In fact, we introduced an expression cassette of AR with a hemagglutinin (HA) tag using a retrovirus in our previous study . Interestingly, the introduction of exogenous AR caused elevated expression of DKK1, indicating that the AR signaling pathway remains intact in HFDPCs even after immortalization . Therefore, we concluded that we successfully obtained genetic mutants of HFDPCs that were negative or positive for AR expression. The comparison of expression profiles between AR-positive and AR-negative immortalized DPCs would allow us to identify the downstream genes of AR. In this study, we comprehensively compared the expression profiles of immortalized HFDPCs with or without AR expression to identify the downstream gene network of AR signaling, which will contribute to finding out of new target molecules to androgen related diseases.
Cell culture and RNA extraction
We previously reported the establishment of immortalized DPCs with the expression of R24C mutant cyclin-dependent kinase 4 (CDK4), cyclin D1, and telomerase reverse transcriptase (TERT) via lentiviral gene transfer. The immortalized human follicle dermal papilla cells were named as “HFDPC_K4DT” from the last characters of the introduced genes (CDK4, cyclin D1, TERT). We previously confirmed that AR expression in established immortalized cells (HFDPC_K4DT) is undetectable, almost identical to that of fibroblasts. To reconstitute the AR signaling pathway, we introduced an expression cassette of AR with an HA tag through the retrovirus expression system. Total RNA was extracted from immortalized HFDPCs using the K4DT method and AR-expressing immortalized HFDPCs. To detect the reproducibility of the expression counts, we extracted RNA from three biological replicates.
We checked the quality of the raw sequencing reads using FastQC. Data were obtained using paired-end sequencing. We analyzed the data from six samples of HFDPC-K4DT and AR-expressing HFDPC-K4DT. Raw sequencing reads were processed with PEAT program to remove the adaptor sequence. We removed the low-quality reads using PRINSEQ, and the remaining reads were mapped to the human reference genome (GRC38, NCBI) with the STAR mapping program. The expression counts for each sample were obtained using the featureCounts program. We first extracted genes which have more than 3000 counts on any sample. To normalize expression counts, we used the DEGES normalization method . We set the parameters of TCC-GUI as default setting [7, 8] The differentially expressed genes (DEGs) were determined using the TCC-GUI program developed by Dr. Koji Kadota (Tokyo University, Tokyo, Japan; https://infinityloop.shinyapps.io/TCC-GUI/). For the determination of differentially expressed genes, we used Q value (FDR) less than 0.01.
Downstream pathway analysis
After mapping all sequencing reads, the expression count of all genes was determined (33,122 genes). Based on the results of DEGs in TCC-GUI, we extracted 1196 genes classified as DEGs. The list of DEGs was processed using the DAVID pathway analysis tool. The expression levels of genes in the listed pathways were visualized using the heatmap function of TCC-GUI.
Detection of expression level of caveolin-1 and EGFR with qPCR
The total RNAs of HFDPC-K4DT and AR expressing K4DT cells were extracted with NucleoSpin RNA (Takara Bio, Shiga, Japan). The cells were treated with basal medium containing 50 nM DHT or without DHT (no treatement) for 8 h. Total cDNAs were obatained with PrimeScript RT reagent Kit with gDNA Eraser with random primer method. The cDNA were used as the template samples of qRT-PCR with THUNDERBIRD SYBR qPCR Mix (TOYOBO, Osaka, Japan) under the reaction condition recommended from the manufacture. The sequences of the detection primers were listed in below. Primers for caveolin 1, TF1162; 5′-CCCGCAGCCTGGGAGTGCCCTGA-3′ and TF1163, 5′-GCTTGTAGATGTTGCCCTGTTCCCGGAT-3′. Primers for EGFP, TF1164, 5′-ATGATGCAAATAAAACCGGACTGAAGGA-3′, TF1165, 5′-CTGCACCCCAGCAGCTCCCATTG-3′. Primers for GAPDH, TF999, 5′- GAGGTGCACCACCAACTGCTTAGC-3′ and TF1000, 5′-TCGGCATGGACTGTGGTCATGAG-3′. The detections with qRT-PCR were carried out with Thermal Cycler Dice Real Time System II (Takara Bio, Shiga, Japan) under the relative quantitation with GAPDH.
Biological background of immortalized HFDPCs and AR-expressing immortalized cells
In our previous study, we exogenously introduced an AR expression cassette through retrovirus gene transfer. Although the expression level of AR was undetectable in parent cells, the AR-expressing immortalized HFDPCs showed an intense signal with the expected molecular weight in western blotting (see Fig. 1 of our previous publication, Fukuda et al., 2020). Furthermore, we detected the nuclear localization of AR even without ligand stimulation, which may explain auto-dimerization based on the force expression system (see Fig. 2 of our previous publication, in Fukuda et al., 2020). Furthermore, we detected activation of DKK1 (a major downstream gene of AR) expression in AR-expressing immortalized HFDPCs, which indicated that the AR signaling pathway was reconstituted by exogenous AR introduction.
Whole-gene transcriptome analysis of parent HFDPC-K4DT and AR-expressing HFDPC-K4DT
To comprehensively compare the expression patterns of whole genes, we carried out RNA-Seq analysis using an Illumina Hiseq X sequencing machine (Illumina, San Diego, CA, USA) and a 150-bp paired end. The sequencing workflow is shown in Fig. 1A. After the removal of the adaptor, we initiated the mapping process. The number of obtained sequence reads was at least 22 M, indicating that the read number was sufficient for quantitation . To evaluate the reproducibility of the data, we carried out RNA-Seq reactions with three biological replicates. We evaluated the quality of the read data using the FASTQC program (Fig. S1 and S2). The results of FASTQC indicated that the average of almost all sequencing data was mapped within the green area, suggesting that the read data was reliable. We next mapped sequencing reads using the STAR program and human reference genome (GRCh38). The mapping ratio and read number are shown in Fig. 1C. The mapping ratio of the samples was more than 95%, indicating that our mapping method was suitable for detecting gene expression. The mapping ratios were 95.9% (HFDPC_K4DT1), 96.2% (HFDPC_K4DT2), 96.3% (HFDPC_K4DT3), 96.4% (HFDPC_K4DT_AR1), 96.3% (HFDPC_K4DT_AR2), and 95.1% (HFDPC_K4DT_AR3).
The complete list of expression counts of parent HFDPC-K4DT and AR-expressing HFDPC-K4DT is provided in Figshare (https://figshare.com/articles/dataset/HFDPC_K4DT_HFDPC_K4DT_AR/13567343). The sequencing data were submitted to the DDBJ database under Bioproject Submission ID PRJDB10909. We first filtered genes at least 3000 counts on any sample. The number of genes remained after the filtration was 1537 genes. Next, we input the expression counts of the whole genome into TCC-GUI. First, we analyzed the correlation plots of expression profiles, as shown in Fig. 1B. The biological replicates of parent HFDPC-K4DT and AR-expressing HFDPC-K4DT formed unique clusters, indicating that the sequencing results were reproducible. The triplicated data also formed unique clusters in three-dimensional PCA (Fig. 1).
Pathway analysis results
We further analyzed the DEGs with a FDR based Q value of less than 0.01. In total, we narrowed down the DEGs to 1196 as candidate genes. The list of DEGs was submitted to the pathway analysis tool DAVID. The most significant pathway determined in DAVID was the focal adhesion pathway (61 hits in the annotation list), followed by Proteoglycans in cancer (51 hits in the annotation list). Based on the pathway analysis results, we compared 61 genes listed in focal adhesion using heatmap analysis (Fig. 2). Furthermore, the expression levels of Proteoglycans in cancer related genes (51 genes) are shown in Fig. 6. We also mapped the DEGs in the Kyoto Encyclopedia of Genes and Genomes (KEGG), which showed more than a 2-fold increase or 0.5-fold decrease (at least 3000 counts in any sample) based on bar plots of expression counts and decision criteria [10,11,12]. Mapping of the focal adhesion pathway indicated that collagen-related molecules (FN1, COL1A1, COL27A1, COL4A1, COL4A2, and COL5A3) were either downregulated or upregulated in AR-expressing cells. Furthermore, the expression level of Caveolin 1 was downregulated and EGFR was upregulated in AR expressing HFDPC-K4DT cell (Fig. 5).
In Proteoglycans in cancer related pathway, the expression level of Twist was elevated which is one of the transcriptional factor in AR expressing cell.
Validation of RNA-Seq results with qPCR analysis
We furthermore detected the expression levels of Caveolin 1 downregulated and EGFR with qPCR analysis. As the first evidence, we detected expression of these two genes in identical RNAs which used for RNA-Seq. The downregulation of Caveolin 1 and upregulation of EGFR in AR expressing cell were reproduced with qPCR analysis (Fig. 3). Furthermore, we carried out treatment of 50 nM of dihydrotestosterone (DHT) with and without AR expressing cells. Under the intact condition, downregulation of Caveolin 1 and upregulation of EGFR in AR expressing cell were also reproduced (Fig. 4A and B, Left side). However, expression of Caveolin 1 was strongly suppressed after the treatment of DHT even in parent K4DT cell. The expression levek of Caveolin 1 was elevated in AR expressing cells. These data indicate that 50 nM DHT treatment causes various types of expression change, which independent to AR signaling.
In this study, we comprehensively compared the expression profiles of parent immortalized human HFDPCs with K4DT and AR-expressing offspring cells. The genetic background of these cell lines was identical except for the expression status of AR. The comparison of global expression profiles allowed us to identify the downstream genes involved in AR signaling.
We already detected that HA-labelled androgen receptor get into the nuclei without treatment of testosterone in AR expressing cell in our previous publication. The nuclear localization of exogenously introduced HA-labelled androgen receptor without ligand treatment can be possibly explained by following two reasons; first possibility is the auto-dimerization. Due to the high level of expression of HA-labelled AR, receptor might cause auto-dimerization, which is critical process for the activation of AR signaling. The second possibility is the testosterone concentration within the serum in cell culture medium. Since the testosterone is quite sensitive as the hormone receptor, we need to perform serum withdrawn to exactly detect the nuclear translocation. We furthermore carried out the exposure of 50 nM DHT to parent K4DT cell and AR expressing cell. The exposure of DHT causes the elevation of Caveolin 1 even in the AR negative parent cell. The results of DHT exposure in parent cell showed that DHT treatment causes changes in signaling pathway, which independent from AR manner. We need to pay attention for the use of the ligand. If we use synthetic antrogen recptor ligand, such as R1881, the effect of ligand treatment might show different response. However, although the synthetic ligand might be specific from the view point of androgen signaling pathway, the biological explanation would be difficult when it compared with natural ligands such as DHT.
We found that caveolin-1 is downregulated in AR-expressing cells. In support of this result, caveolin-1 is reported to be controlled by AR signaling . Furthermore, caveolin-1 has been identified as a malignant marker of prostate cancer, and it controls the survival ratio of prostate cancer cells [14, 15]. In addition, in a mouse study, caveolin-1 expression was suppressed after testosterone treatment. The expression level of caveolin-1 increased in castrated male mice . In the mouse genome, the two binding sites were identified within intron 2, indicating that caveolin-1 expression is controlled by AR. Although the association of caveolin-1 and AR has been suggested in previous studies, the connection between these two molecules has not yet been fully elucidated. Our expression profiling data suggest that the strength of the AR signaling pathway is controlled by caveolin-1. We also identified the expression level of EGFR is upregulated in AR expressing cell. Identification of EGFR as the downstream suggest the cross-talk of AR and EGF signals. These data indicate the existence of a gene network under the control of the AR nuclear receptor. Improved understanding of AR-related networks may contribute to the discovery of new therapeutic targets for AR-related diseases, such as prostate cancer or AGA.
Our whole transcriptome data is critical resources for discovery of new therapeutic targets for testosterone-related diseases.
Availability of data and materials
The datasets generated during and/or analysed during the current study are available in the Figshare and DNA data base of Japan (DDBJ repository, [https://figshare.com/articles/dataset/HFDPC_K4DT_HFDPC_K4DT_AR/13567343] and Bioproject Submission ID PRJDB10909, https://www.ncbi.nlm.nih.gov/bioproject/686284.
Cyclin-dependent kinase 4
Dickkopf-related protein 1
False discovery rate
Graphical user interface
Human follicle dermal papilla cells
Principle component analysis
Tag count comparison
Telomerase reverse transcriptase
Tumor growth factor-β2
Krupa Shankar DS, Chakravarthi M, Shilpakar R. Male androgenetic alopecia: population-based study in 1,005 subjects. Int J Trichol. Wolters Kluwer -- Medknow Publications. 2009:131–3. https://doi.org/10.4103/0974-7753.58556.
Kwack MH, Sung YK, Chung EJ, Im SU, Ahn JS, Kim MK, et al. Dihydrotestosterone-inducible dickkopf 1 from balding dermal papilla cells causes apoptosis in follicular keratinocytes. J Invest Dermatol. 2008;128:262–9. https://doi.org/10.1038/sj.jid.5700999.
Pihlajamaa P, Sahu B, Jänne OA. Determinants of receptor- and tissue-specific actions in androgen signaling. Endocr Rev. 2015;36:357–84. https://doi.org/10.1210/er.2015-1034.
Yassin A, AlRumaihi K, Alzubaidi R, Alkadhi S, Al AA. Testosterone, testosterone therapy and prostate cancer. Aging Male. 2019;22:219–27. https://doi.org/10.1080/13685538.2018.1524456.
Fukuda T, Takahashi K, Takase S, Orimoto A, Eitsuka T, Nakagawa K, et al. Human derived immortalized dermal papilla cells with a constant expression of testosterone receptor. Front cell. Dev Biol. 2020;8. https://doi.org/10.3389/fcell.2020.00157.
Kadota K, Nishiyama T, Shimizu K. A normalization strategy for comparing tag count data. Algorithms Mol Biol. 2012;7:5.
Sun J, Nishiyama T, Shimizu K, Kadota K. TCC: an R package for comparing tag count data with robust normalization strategies. BMC Bioinformatics. 2013;14:219.
Su W, Sun J, Shimizu K, Kadota K. TCC-GUI: a shiny-based application for differential expression analysis of RNA-Seq count data. BMC Res Notes. 2019;12:133.
Fukuda T, Doi K, Donai K, Takahashi K, Kobayashi H, Hirano T, et al. Data descriptor: global transcriptome analysis of pig induced pluripotent stem cells derived from six and four reprogramming factors. Sci Data. 2019;6. https://doi.org/10.1038/sdata.2019.34.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;11:1947-51. https://doi.org/10.1002/pro.3715.
Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: Integrating viruses and cellular organisms. Nucleic Acids Res. 2021:49.
Shuang-Gang H, Guang-Xin Yao YS. [androgen/androgen receptor directly regulates the expression of Caveolin-1 in mouse epididymides] - PubMed. Zhonghua Nan Ke Xue. 2013;10:867–72 https://pubmed.ncbi.nlm.nih.gov/24218937/. Accessed 16 Jan 2021.
Ayala G, Morello M, Frolov A, You S, Li R, Rosati F, et al. Loss of caveolin-1 in prostate cancer stroma correlates with reduced relapse-free survival and is functionally relevant to tumour progression. J Pathol. 2013;231:77–87. https://doi.org/10.1002/path.4217.
Bryant KG, Camacho J, Jasmin JF, Wang C, Addya S, Casimiro MC, et al. Caveolin-1 overexpression enhances androgen-dependent growth and proliferation in the mouse prostate. Int J Biochem Cell Biol. 2011;43:1318–29. https://doi.org/10.1016/j.biocel.2011.04.019.
We thank Dr. Koji Kadata (Tokyo University) for the technical help for the gene expression analysis. We also thank Dr. Taku Ozaki and Dr. Tetsuro Yamashita (Iwate University) for the mentoring of students in our laboratory.
This work was supported in part by a basic managing budget of Iwate University. The founder has no role in the study design, decision to publish or data production.
Ethics approval and consent to participate
This study used the human derived primary cell. However, the original primary cell were obtained from PromoCell (Heidelberg, Germany) through the local distributor (Takara Bio, Shiga, Japan). Although we asked the necessity of Ethics approval of human derived samples to comity of Iwate university, the final conclusion of comity was “This study do not require the approval of ethics, since the sample is commercially distributed cells”. Based on this final conclusion, we did not submit the paperwork to the approval comity.
Consent for publication
All authors do not have any conflict of interest and competing interest to disclose.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Furuya, K., Fujibayashi, S., Wu, T. et al. Transcriptome analysis to identify the downstream genes of androgen receptor in dermal papilla cells. BMC Genom Data 23, 2 (2022). https://doi.org/10.1186/s12863-021-01018-6
- Androgen receptor
- Dermal papilla cells