- Research article
- Open Access
Transfection of Sertoli cells with androgen receptor alters gene expression without androgen stimulation
BMC Molecular Biology volume 16, Article number: 23 (2015)
Androgens play an important role for the development of male fertility and gained interest as growth and survival factors for certain types of cancer. Androgens act via the androgen receptor (AR/Ar), which is involved in various cell biological processes such as sex differentiation. To study the functional mechanisms of androgen action, cell culture systems and AR-transfected cell lines are needed. Transfection of AR into cell lines and subsequent gene expression analysis after androgen treatment is well established to investigate the molecular biology of target cells. However, it remains unclear how the transfection with AR itself can modulate the gene expression even without androgen stimulation. Therefore, we transfected Ar-deficient rat Sertoli cells 93RS2 by electroporation using a full length human AR.
Transfection success was confirmed by Western Blotting, immunofluorescence and RT-PCR. AR transfection-related gene expression alterations were detected with microarray-based genome-wide expression profiling of transfected and non-transfected 93RS2 cells without androgen stimulation. Microarray analysis revealed 672 differentially regulated genes with 200 up- and 472 down-regulated genes. These genes could be assigned to four major biological categories (development, hormone response, immune response and metabolism). Microarray results were confirmed by quantitative RT-PCR analysis for 22 candidate genes.
We conclude from our data, that the transfection of Ar-deficient Sertoli cells with AR has a measurable effect on gene expression even without androgen stimulation and cause Sertoli cell damage. Studies using AR-transfected cells, subsequently stimulated, should consider alterations in AR-dependent gene expression as off-target effects of the AR transfection itself.
Androgens play a pivotal role for the development of the male phenotype, the initiation and maintenance of spermatogenesis and therefore male fertility (for review see ). The action of the most important androgens testosterone (T) and dihydrotestosterone (DHT) is mediated by the androgen receptor (AR/Ar). It is a ligand-activated transcriptional factor belonging to the nuclear receptor superfamily. The AR/Ar gene is located on the X chromosome and consists of eight exons, coding for the N-terminal transcription regulation domain, the DNA binding domain (DBD) in the middle of the protein and the C-terminal ligand binding domain (LBD). The DBD as well as the LBD are highly conserved throughout species (for review see ). Bound to its ligand, the androgen-AR complex is translocated into the nucleus, binds to the DNA (androgen responsive elements, AREs) and is able to activate or repress gene expression by recruiting co-activators or co-repressors (for review see ). The activity of steroid hormone receptors is also regulated by post-transcriptional modifications. In case of AR/Ar, a great variety of these modifications has been described, i.e. phosphorylation, acetylation, ubiquitination and also methylation (for review see ).
The AR/Ar is expressed in all tissues except the spleen (for review see ). In the testis, it is expressed in interstitial Leydig cells and endothelial cells, as well as in peritubular myoid cells and tubular Sertoli cells , for review see . Since germ cells do not express AR/Ar, the androgen action has to be mediated towards the germ cells by Sertoli cells. These somatic cells have been described as branched cells surrounding all germ cell stages [6, 7]. As was shown by Willems et al. , a selective ablation of Ar in mouse Sertoli cells (SCARKO) leads to a disturbed Sertoli cell maturation including a delayed and defective establishment of the blood-testis barrier. Moreover, no meiotic germ cells were observed in SCARKO mice, showing the importance of a functional AR/Ar on Sertoli cell biology and for the development of germ cells.
To examine the role of the AR/Ar in different biological processes such as cell growth and survival as well as AR/Ar-dependent gene expression, cell culture systems are needed. Therefore, administration of T and/or the more efficient metabolite DHT has widely been used to investigate the effect of androgens and AR/Ar, respectively, in diverse cultured cells such as human breast cancer cells, adrenocortical carcinoma cells, murine skeletal muscle cells or liver carcinoma cells [9–12]. Additionally, AR/Ar-deficient cell lines were used, e.g. AR-deficient MCF-7 breast cancer cells, to examine the effect on estrogen administration in a system lacking AR . Szelei et al.  transfected AR-deficient MCF-7 breast cancer cells with human AR and showed an inhibition of proliferation. Also prostate cancer cells devoid of AR were transfected with human AR and showed a decreased proliferation rate . The question is, whether the transfection procedure itself might have led to an altered expression of AR/Ar-dependent and AR/Ar-independent genes. Xiao et al.  demonstrated equal concerns in Amh-Cre-transfected mouse, where an increase of oxidative stress and lipid peroxidation in Sertoli cells was detected even without stimulation.
Beside “classical” androgen-dependent prostate cancer development, also androgen-independent signalling pathways gained increasing interest as shown recently by Li et al. . The authors describe a persistent transcriptional activity in castration-resistant prostate cancer cell lines in the absence of androgens. This transcriptional activity was mediated by a truncated AR protein lacking the LBD. This raises the question, whether an androgen-independent AR/Ar action is always important in cell biology and which genes might be expressed or repressed by AR/Ar presence alone.
For this purpose, we transfected rat Sertoli cells which have been shown to be deficient of Ar with full length human AR DNA. After transfection, we performed genome-wide microarray analysis and compared the gene expression pattern with non-transfected Sertoli cells to identify a possible “intrinsic” activity of AR/Ar without androgen administration. We found significantly altered gene expression in transfected compared with non-transfected cells, possibly influencing Sertoli cell function.
Transfection of 93RS2 cells with the human AR
Success of transfection with full length human AR CDS was validated by immunofluorescence (IF, Fig. 2a), Western Blot (Fig. 2b) and RT-PCR (Fig. 2c). As the commercially available human AR was introduced in a GFP-coupled vector system, we used a rabbit anti-GFP antibody for IF experiments in transfected cells whereas non-transfected cells were used as internal negative control. Using PAGE, we were able to show the CAG repeat length of 17 to be stable throughout different settings (Fig. 2d).
Microarray analysis revealed an altered gene expression in transfected 93RS2 cells
Microarray analysis revealed 672 significantly regulated genes (p < 0.01 and fold change (FC) >2.0). Of these, 200 genes showed higher gene expression values, whereas 472 revealed a lower gene expression in 93RShAR17 cells compared with non-transfected cells.
Hierarchical clustering of the 672 significantly regulated genes shows two clusters clearly differentiating between transfected and non-transfected cells (Fig. 3). Three biological replicates have been tested and show a homogeneous expression pattern, indicating high reproducibility of microarray results. An overview of the ten highest regulated genes for down- and up-regulation is given in Table 1. Complete array data may be found following the link provided .
“Development”, “Hormone response” and “Immune response” are the predominant functions of the differently regulated genes
Of 370 annotated down-regulated genes, 330 could be assigned to DAVID functional categories, and 124 out of 142 annotated up-regulated genes, respectively.
An overview of the functional categories that have been inferred with DAVID is given in Table 2. Down-regulation is predominant in “Cell development/Cell contact”, “Response to hormone stimulus” and “Nucleotide catabolic process”, whereas regulation is evenly distributed in “Immune response”. The highest score values are achieved by four significantly overrepresented gene ontology (GO) categories clustered under “Biological adhesion”, whereas the highest number of genes is assigned to 15 GO categories grouped as a cluster named “Epithelium development”. More than half of the functional assigned groups belong to cell development and cell contact while 25 % of the functionally assigned genes are related to immune response. 36 genes can be attributed to “Hormone stimulus” and a minority of 12 genes contributes to “Nucleotide catabolic process”.
Upstream regulation analysis identified more activation than de-activation
Upstream regulation analysis with IPA is based on gene expression patterns and predicts activation or deactivation of regulators of the differentially regulated genes. The results show that more upstream regulators are predicted to be activated (n = 51) than inhibited (n = 20).
These predictions are based on 220 genes from which 95 contributed to activation as well as to deactivation. The proportion of overall down- and up-regulation is mirrored in these genes with more down-regulation in inhibition as well as in activation (Tables 3, 4, 5, 6). The majority of deactivated upstream regulators (8 out of 20) are classified as transcription regulators. Activation is mainly predicted for cytokines (14 out of 51).
Validation of microarray data by RT-qPCR
For validation of microarray results we performed RT-qPCR for 22 candidate genes, showing different ranges of regulation (up, down). Among the chosen genes, some are mainly associated with development and are known Sertoli cell markers, such as Dhh , Gja1 , Inhbb , and Tf . Other genes are markers for differentiation and proliferation (e.g. Bambi and Tgfb1i1 ) and some are involved in apoptosis, such as Myc and Tnfrsf1a . Results from RT-qPCR were mostly consistent with data from microarray analysis (Fig. 4). Relative gene expression was lower in transfected compared to non-transfected Sertoli cells in 13 of 22 cases. Gene expression of Cdkn1a, Egr1, Fst, Gja1, Myc, Pmepa1, Ptsg2, Rarg and Tnfrsf1a was higher in 93RShAR17 cells compared to Ar-deficient 93RS2 cells. In the latter case, it has to be mentioned, that differences of the means did not reach significance in four genes, due to high variability of Cq.
To study the effects of androgens and AR/Ar on diverse cell culture systems and the relevance for cell biology, cell culture experiments were conducted in different human cell lines (e.g. breast cancer cells, adrenocortical carcinoma cells, murine skeletal muscle cells or liver carcinoma cells [9–12]). Also AR-deficient cell lines have been used, either transfected with AR [14, 15] or without . Both groups working with transfected cell lines performed their experiments using either not stimulated  or mock-transfected cells  as negative controls. Moreover, Jacobsen et al.  showed, that transfection of MCF-7 breast cancer cells lead to severe differences in gene expression levels in distinct genes, depending on the transfection reagent used. Interestingly, transfection with a vector encoding for a reporter gene and a vector without insert, respectively, revealed no differences in gene expression. This implies, that the transfection procedure itself might alter gene expression in these cells. Therefore, we performed gene expression analysis with AR-transfected rat Sertoli cells using non-transfected cells as controls to show “intrinsic” gene expression alterations due to the transfection procedure. As electroporation has been shown to be superior with respect to cell viability and also transfection efficiency compared to chemical transfection using lipofectamine , we applied this technique to introduce the AR. Cell viability was not influenced by electroporation, but whole genome microarray analysis showed an altered gene expression. Surprisingly, more genes have been down-regulated than up-regulated comparing transfected and non-transfected cells. We selected 22 genes showing an altered expression pattern and confirmed microarray results with RT-qPCR analysis. In the following, we will discuss in more depth interesting genes and pathways, respectively.
Among the down-regulated genes, many are involved in metabolic processes, as for example in iron transport and metabolism (cytochrome b reductase 1 (Cybrd1), FC = −107; transferrin (Tf), FC = −6,898; six-transmembrane epithelial antigen of the prostate 2 (Steap2), FC = −2.3). Iron is relevant for Sertoli cells in two different aspects: as supervisors of germ cell development, Sertoli cells provide iron which is needed for DNA synthesis and cell growth of germ cells that undergo multiple mitotic divisions . On the other hand, Sertoli cells avoid toxic environmental conditions that might be given at elevated concentrations of insoluble ferric iron (Fe3+). Therefore, Sertoli cells secrete transferrin, the product of the Tf gene , which may be used as a marker for Sertoli cell function and differentiation  as it creates an environment low in free iron that impedes bacterial survival in a process called iron withholding. The protein level of Tf decreases in inflammation. The lower expression of Tf gene, which was confirmed in RT-qPCR (Fig. 4), could be interpreted as a sign of severe disturbance and inflammation of cells.
The latter is reflected by the high proportion of upstream regulators related to immune response (= cytokines and members of the MAP kinase signalling pathway) that are predicted to be activated (Table 3) and the presence of multiple immune response-related genes on top of the list in up-regulation (Table 1). “Immune response” is the second huge cluster of altered genes in our study, represented by e.g. prostaglandin-endoperoxidase synthase 2 (Ptgs2, FC = 3.558) also known as cyclooxygenase 2 (Cox2). An increase in Cox2 expression was observed by Matzkin et al.  in Leydig cells of infertile men showing either hypospermatogenesis, Sertoli cell only syndrome or maturational arrest. By increased numbers of testicular macrophages, levels of interleukin 1β (Il-1β) are increased and activates Ptgs2, the key enzyme in prostaglandin synthesis culminating in inflammation. The expression of Tf, Ptgs2 and interleukins is coupled in Sertoli cells; as shown by Yamaguchi et al. , incubation with cisplatin lead to an increase in Ptgs2 and a decrease in Tf expression in Sertoli cell cultures, similar to our study. Additionally, an analysis of upstream regulation using IPA revealed a high number of key players in inflammation to be activated showing congruently high FCs for Ccl5, Irf7, and Ifnb1. This might on the one side be due to the transfection procedure itself and/or reflect inflammatory processes in the cells due to increased cell damage. Remarkably, an influence of molecular biological techniques on gene expression and immune response has been observed also in regard to short-interfering RNAs (siRNAs). Sledz et al. reported an induction of interferon β levels in a human glioblastoma cell line which was transfected with siRNAs as a non-specific side effect additionally to the silencing of the target gene lamin  .
Not only metabolism and immune response gene expression seem to be altered in transfected Sertoli cells, but also cell cycle and development genes (desert hedge hog (Dhh) FC = −2.032; fibroblast growth factor receptor 2 (Fgfr2) FC = −8.239; follistatin (Fst) FC = 2.162; inhibin beta b (Inhbb) FC = −3.126). Dhh is involved in various areas of embryonic development, including testicular cord formation. It is expressed in mouse Sertoli cell precursors during mid- to late gestation  and also is important for germ cell development after puberty in mouse  and rat testis . In the mouse, a lack of Dhh results in a severe impairment of spermatogenesis due to a lack of spermatogonial development beyond primary spermatocytes . Fgfr2 is a known differentiation factor in prenatal Sertoli cells as it is concomitantly expressed with Sry and is essential for subsequent expression of anti-muellerian hormone (Amh) and Sox9 . IPA analysis of upstream regulation predicted an inhibition of transcription factor Sox9 with a z-score of −2.2 (Table 6). Moreover, lack of Fgfr2 might cause a partial XY sex reversal, as loss of Fgfr2 leads to an up-regulation of Follistatin (Fst), a female somatic cell marker , which was confirmed by microarray and RT-qPCR. A down-regulation of the Sertoli cell marker Inhbb (for review see ) also points to a decreased Sertoli cell function and a severe disturbance of spermatogenesis in the rat . Figure 5 shows the association of Inhbb, Fst, Dhh, Pmepa1, Fgfr2, Ptgs2, Tf and Myc as especially interesting genes on known pathways as predicted by IPA.
A disturbance of Sertoli cell function is also visible in gene expression alteration concerning the functional cluster “Cell adherence” or “Cell adhesion” (Fig. 6, e.g. collagen type IV alpha (Col4a1) FC = −12.503; gap junction protein 1 (Gja1) FC = −1.188). Cell adhesion and formation of tight junctions between Sertoli cells generating the blood-testis barrier is one of the most important features of Sertoli cell maturation and function (for review see ) as it is a prerequisite for intact spermatogenesis. Also cell-to-cell contact and communication seem to be disturbed in transfected cells as indicated by the down-regulation of Gja1, also known as connexin 43 (for review see ).
Our results indicate a severe disturbance of Sertoli cell metabolism, function and cell biology concerning immune status and generation of blood-testis barrier, caused by the transfection procedure even without androgen stimulation of cells. The alterations in gene expression levels might either be related to the transfection procedure itself and/or to the insertion of human AR into Ar-free rat Sertoli cells. A microarray analysis with mock-transfected Sertoli cell line would be needed to distinguish both possibilities. We consider the altered gene expression to be caused by AR insertion, as many of the altered genes were identified as AR and Sertoli cell specific. In either case, incubation of transfected cell lines with testosterone or dihydrotestosterone might lead to false-positive or false-negative results and additionally, also non-genomic pathways including AR/Ar action may be altered by transfection procedures. Therefore, suitable negative controls are needed for stimulation experiments with T or DHT, i.e. non-transfected cells as appropriate negative controls. Gene expression has to be normalized by these non-transfected cells to avoid false-positive or false-negative results regarding gene regulation.
Cell lines and culture conditions, human positive control tissue
We determined the expression of Ar in different Sertoli cell lines by RT-PCR. For this study, we used four existing immortalized Sertoli cell lines from either mouse (WL3, SK-11) [40, 41] or rat testis (93RS2)  kindly provided by our collaborators. Additionally, SCIT-C8 cells were generated from immortalized Sertoli cells from rat testis as described by Konrad et al. . We did not conduct any animal research in our study and therefore ethics approval was not required. Total RNA of these cells was isolated by peqGold Total RNA Kit (Peqlab, Erlangen, Germany), set to a concentration of 200 ng/µl and genomic DNA was digested by RNase-Free DNase Set (Qiagen, Hilden, Germany). Reverse transcription was performed with Omniscript RT Kit (Qiagen). The mastermix was prepared as follows: 2 µl Buffer RT (10×), 2 µl dNTP mix (5 mM each), 0,7 µl RNAse inhibitor (20 units/µl, Invitrogen via LifeTechnologies, Carlsbad, CA, USA), 0,2 µl Oligo-dT primer (10 µM, Qiagen) and 1 µl Omniscript Reverse Transcriptase were mixed and RNase free water was added to a final volume of 10 µl. To test genomic DNA digestion success, we omitted reverse transcriptase and replaced it by RNase free water for one reaction. After addition of 1 µl RNA (200 µg/µl), we incubated the reaction mix for 1 h at 37 °C. cDNA not directly used for further experiments was stored at −20 °C. Amplification of Ar was achieved with a matching primer pair for murine and rat Ar obtained from Eurofins MWG Operon (Huntsville, AL, USA) as can be seen in Table 7 and Taq PCR Master Mix Kit (Qiagen). Mastermix was prepared as follows: 34 µl Taq PCR Master Mix, forward and reverse primer (2,5 µl each) and RNase free water as well as 1 µl cDNA were mixed to a final volume of 50 µl. Amplification was performed with 1× 94 °C for 4 min, 35× (94 °C for 40 s, 60 °C for 45 s, 72 °C for 90 s) and 1× 72 °C for 5 min.
As the prepubertal rat Sertoli cell line 93RS2 proved to be devoid of Ar, we chose this cell line for further experiments. The cells were maintained in a 5 % CO2 atmosphere at 34 °C. The standard culture media consists of DMEM high glucose mixed 1:1 with Ham’s F-12 media plus 100 units/ml penicillin, 0,1 mg/ml streptomycin, 10 % FBS-Gold (total protein 3.0–4.5 g/dl), and 1 % ITS (1000 mg/l Insulin, 550 mg/l Transferrin, 0.68 mg/l Selenin). Unless otherwise stated, cell culture media were purchased from Invitrogen (via Life Technologies, Carlsbad, CA, USA).
Ethics, consent and permissions
For positive control used in RT-PCR and Western Blotting, we used testis homogenate from a patient showing normal spermatogenesis attending the andological clinic in Münster for re-fertilization surgery. After written informed consent, biopsies were taken under general anesthesia. The reported study has been approved by the Ethics committee of the Medical Faculty of the Justus Liebig University Giessen (decision 75/00 and 56/05).
AR transfection in 93RS2 Sertoli cells
We introduced a commercial available full length human AR (OriGene, Rockville, MD, USA), containing 17 CAG triplets, into the expression vector pcDNA 6.2 C-EmGFP (Invitrogen) after amplification of AR using GC-Rich PCR System (Roche, Basel, Switzerland), according to manufacturer’s instructions. Transfection of 93RS2 cells was performed using the microporation system MP-100 (Peqlab). Cells were detached by Trypsin (PAA, Piscataway, NY, USA) and subsequently adjusted to 200,000 cells per well in a 6-well-plate. After re-suspending the cells in the provided buffer, plasmid DNA was added. We used a current strength of 1150 V for 20 ms with two pulses.
Validation of transfection success in 93RS2 by immunofluorescence, RT-PCR and Western Blotting
24 h after transfection, transfected cells (93RShAR17) were fixed in 6-well-plates with 4 % paraformaldehyde for 20 min at room temperature, washed three times with PBS and permeabilized with 0.1 % Triton x-100. After transferring the cells to a 12-well-plate and washing with PBS, unspecific binding sites were blocked with 3 % BSA (bovine serum albumin, Carl Roth GmbH + Co.KG, Karlsruhe, Germany) in TBST (Tris-Buffered Saline and Tween 20, Carl Roth) and washed again with PBS. The rabbit anti-GFP antibody (ab290, Abcam, Cambrigde, UK) was added in a dilution of 1:200. After incubation for 3 h and washing with PBS, goat anti-rabbit Alexa 488 antibody (Invitrogen) was added in a dilution of 1:200. After a final incubation for 1 h in the dark, cells were washed and embedded with Vectashield mounting medium with DAPI (H-1200, Vector Laboratories, Dossenheim, Germany). Transfection efficiency was evaluated using a fluorescence microscope (AxioPhot, Zeiss, Oberkochen, Germany). Western Blot analysis to prove antibody specificity and AR protein expression in transfected Sertoli cells was performed as described elsewhere . Shortly, proteins extracted from cell lysates of transfected 93RS2 cells and human testis tissue were submitted to protein extraction using TRI Reagent® RNA Isolation Reagent (Sigma-Aldrich, St. Louis, MO, USA) according to Chomczynski . Proteins were run on a 3–8 % Tris–acetate gel (Life Technologies, Carlsbad, CA, USA) for 75 min at 150 V and blotted on nitrocellulose membrane for 75 min at 30 V. A polyclonal rabbit anti-human AR antibody (sc-816, Santa Cruz Biotechnology Inc., Dallas, TX, USA) in a 1:500 dilution and a biotinylated goat anti-rabbit antibody (E0432, Dako, Glostrup, Denmark) in a 1:1000 dilution were used. As weight marker, we used HiMark™ Pre-Stained Protein Standard (Life Technologies). Signal detection was performed by incubating the membrane with Vectastain Elite ABC Standard Kit (Vector Laboratories, Inc., Burlingame, CA, USA) and TrueBlue™ Peroxidase Substrate (KPL, Gaithersburg, MD, USA). A negative control was performed by omitting the primary antibody.
To detect AR mRNA in transfected cells, we performed RT-PCR (primers may be seen in Table 7) as described earlier with minor changes concerning the cycling conditions: 1× 94 °C for 4 min, 35× (94 °C for 45 s, 55 °C for 45 s, 72 °C for 90 s) and 72 °C for 5 min resulting in a 591 bp amplicon. The CAG repeat length was confirmed using RT-PCR with subsequent high resolution polyacrylamide gel electrophoresis (PAGE) as described recently .
RNA isolation for microarray analysis
Total RNA of transfected 93RShAR17 cells as well as of non-transfected 93RS2 cells (using three technical replicates (N1-N3) each) was extracted using the peqGold total RNA kit (Peqlab) following manufacturer’s instructions. The amount of RNA was measured on a BioPhotometer (Eppendorf, Hamburg, Germany) as follows: 93RS2 N1 2200 ng/µl, N2 2130 ng/µl and N3 1920 ng/µl and 93RS2hAR17 N1 990 ng/µl, N2 1150 ng/µl and N3 1065 ng/µl (each replicate with a total volume of 15 µl). RNA was stored after extraction until use at −80 °C and transported in liquid nitrogen. The quality of total RNA was checked on a 1 % agarose gel stained with ethidium bromide (Sigma-Aldrich) as well as on Agilent 2100 Bioanalyzer using Eukaryote Total RNA Nano Assay (Agilent Technologies, Santa Clara, CA, USA). For this purpose, RNA was diluted to a concentration of 300 ng/µl. Only high quality RNA samples were used for microarray analysis.
cRNA synthesis and hybridization
Extracted RNA was transcribed into biotinylated cRNA using MessageAmp™ II-Biotin Enhanced Kit (LifeTechnologies). Biotinylated cRNA again was quality checked on Agilent 2100 Bioanalyzer as stated above followed by cRNA fragmentation and finally hybridization on CodeLink Rat Whole Genome using the CodeLink Expression Assay Kit (GE Healthcare, Chalfont St. Giles, Buckinghamshire, UK). For this, 10 µg cRNA was diluted with nuclease-free water to final volume of 20 µl and mixed with 5 µl fragmentation buffer (taken from CodeLink iExpress iAmplify cRNA Prep & Hyb Kit, GE Healthcare) and fragmented at 94 °C for 20 min and subsequent cooling to 0 °C on ice. Hybridization solution was prepared by mixing hybridization buffer component A and B (taken from CodeLink iExpress iAmplify cRNA Prep & Hyb Kit), nuclease-free water and 25 µl fragmented cRNA. Denaturation of cRNA was performed at 90 °C for 5 min with subsequent cooling on ice. Hybridization reaction was carried out at 37 °C for 18 h. Subsequent washing was performed with 0.75 × TNT (1 M Tris–HCl, 5 M NaCl and 20 % Tween 20) buffer. Bioarrays were stained with Cy5™-streptavadin (GE Healthcare) and scanned using the GenePix® 4000 B scanner and the GenePix Pro 4.0 Software (Axon Instruments, Arlington, USA). Scan resolution was set to 5 microns. A total of 2 × 3 = 6 array images were subjected to data analysis. Spot signals of CodeLink bioarrays were quantified using the CodeLink System Software 126.96.36.199312 which generated local background corrected raw as well as median centred intra-slide normalized data.
Quality control of microarray data
The genes represented by probe sets were annotated using the biocLite package (BioConductor) with the library “rwgcod.db” for CodeLink Rat Whole Genome arrays. The intra-slide normalized data containing 35129 rows and 6 columns (200 k values) were processed by an automated workflow that includes omission of control genes (n = 1280), removal of genes with poor QC (n = 1300 values, 0.6 %) or negative sign (n = 1603 values, 0.8 %), removal of probe sets with too high proportion (≥50 %) of missing values per group (n = 203 probe sets, 0.5 %) or with not any group having at least 50 % of values flagged as “G = good” and 50 % values above threshold (n = 7177 probe sets, 21.2 %), removal of outliers (expression values deviating more than fourfold from the group median, n = 427 values, 0.3 %). A total of 26452 probe sets remained after quality control with 1257 probe sets (=4.7 %) containing 1235 missing values (=0.8 %).
Remaining missing values were imputed by probabilistic principal component analysis (PPCA) using the R-package pca Methods. Imputed dataset was quantile normalized using the R-package limma , and logarithm for the base 2 was calculated.
Differential gene expression
Students t test was applied and a false discovery rate (FDR) ≤0.01 was set for the significance level with an absolute fold change (FC) ≥2 between transfected and non-transfected cells.
Functional gene analysis: overrepresentation analysis
Enriched functional gene ontology (GO) categories within the differentially regulated genes were determined using DAVID version 6.7 [47, 48]. Functional annotation clustering as well as an enrichment score was calculated for each cluster.
Upstream regulation analysis
To identify the regulators responsible for the observed gene expression profiles, we performed prediction analysis for activation or inhibition of upstream regulators using the Ingenuity® Pathway Analyzer and the Ingenuity®Knowledge Base (IPA, Qiagen). Prediction is given as a z-score with >2 for activated and <2 for inactivated upstream regulators.
Validation of microarray results by quantitative RT-PCR (RT-qPCR)
For validation of microarray data, we performed RT-qPCR with 93RShAR17 and non-transfected 93RS2 cells for 22 genes (Table 7) that have been shown to be significantly altered in microarray analysis. All primer pairs obtained from MWG Operon have been validated in standard RT-PCR using rat testis as positive control. For this pupose, total RNA from rat testis was extracted using TRI Reagent® RNA Isolation Reagent (Sigma-Aldrich) according to Chomczynski . Genomic DNA was digested by using DNase I (Roche). For this, 6,65 µl RNA (200 ng/µl) were incubated with 1 µl MgCl2 (25 mM, Thermo Fisher Scientific), 1 µl DNase Buffer (Roche), 0,25 µl RNase inhibitor (40 units/µl, Thermo Fisher Scientific) and 1 µl DNase I for 25 min at 37 °C in a thermocycler. After a enzyme heat inactivation for 5 min at 75 °C, RNA was immediately reversely transcribed into cDNA. For this, 1,5 µl DNase-treated RNA was mixed with 1 µl 10x PCR Gold Buffer, 2 µl MgCl2 (25 mM), 1 µl dNTP mix (each 2,5 mM), 0,5 µl random hexamer primer (50 mM), 0,5 µl RNase inhibitor (20 units/µl), 0,5 µl MultiScribe® Reverse Transcriptase (50 units/µl) and RNase free water to a final volume of 9 µl. All reagents were obtained from Thermo Fisher Scientific. For -RT control, reverse transcriptase was replaced by the same amount of RNase free water. Incubation was performed as follows: 8 min at 21 °C, 15 min at 42 °C and 5 min at 99 °C. cDNA was stored at −20 °C until use. For primer validation in standard RT-PCR, 1 µl cDNA was mixed with 2,5 µl 10× PCR Gold Buffer, 2 µl MgCl2 (25 mM), dNTP mix (each 2.5 mM), 1 µl forward and reverse primer, respectively (each 10 pM), 0.125 µl AmpliTaq Gold® DNA Polymerase (5 units/µl) and RNase free water to a final volume of 25 µl. Cycling conditions were: 1× 94 °C for 9 min, 35× (94 °C for 45 s, 60 °C for 45 s, 72 °C for 45 s) and 72 °C for 5 min. Length of the resulting amplicons was checked in an agarose gel electrophoresis as described earlier. For RT-qPCR dilution series we used rat Rplp and Ubc as internal reference genes and performed triple determination in a decreasing 10- fold dilution series (undil., 1:10, 1:100). RT-qPCR efficiency (E) has been calculated using Bio-Rad CFX Manager version 3.1 (Bio-Rad) from the standard curve’s slope and may be seen in Table 7. Reference genes have been determined by using a TaqMan® Array Rat Endogenous Control Plate (96-well, 32 reference genes pre-plated, Applied Biosystems via Thermo Fisher Scientific, Waltham, MA, USA).
For RT-qPCR, total RNA from transfected and non-transfected cells was extracted using peqGold Total RNA Kit (PEQlab) and reversely transcribed into cDNA as described above. As technical replicates we used cell pellets from three independent passages and for each specimen, double determination was performed using 1 µl of cDNA, 4 µl EvaGreen mastermix (no Rox) (Bio&Sell, Feucht, Germany), 0.6 µl forward and reverse primer each and 12.8 µl sterile aqua bidest to a final volume of 20 µl. RT-qPCR conditions were 1× 95 °C for 15 min, 40× (95 °C for 15 s, 60 °C for 30 s, 72 °C for 20 s) followed by melt curve analysis (1× 95 °C for 10 s, 65 °C to 95 °C, increment 0.5 °C for 5 s) on a CFX96 RealTime cycler (Bio-Rad Laboratories, Hercules, CA, USA). Relative gene expression was calculated by the 2−ΔΔCq method, using Rplp and Ubc as internal reference genes. Expression levels represent x fold higher expression in the transfected than in the non-transfected cells (set as “1”). For statistical analysis, differences of the mean were assessed by ANOVA analysis. P-values of p ≤ 0.05 are set as statistically significant. The Cq values for all transcripts may be seen in Additional file 1: Table S1.
Availability of supporting data
Complete microarray data may be found on GEO Accession Viewer database  with accession number GSE57653. Single Sertoli cell line data may be found under accession numbers GSM1385418 (Sertoli Cell Line noAR_1), GSM1385419 (Sertoli Cell Line noAR_2), GSM1386001 (Sertoli Cell Line noAR_3), GSM1385420 (Sertoli Cell Line AR17_1), GSM1385421 (Sertoli Cell Line AR17_2), GSM1385422 (Sertoli Cell Line AR17_3). Raw data of RT-qPCR experiments can be seen in Additional file 1: Table S1.
- 93RS2, SCIT-C8:
rat Sertoli cell lines
with human AR transfected 93RS cells
androgen insensitivity syndrome
- Amh :
- AR/Ar :
androgen responsive element
- Bambi :
BMP and activin membrane-bound inhibitor
- Ccl5 :
chemokine (C–C motif) ligand 5
- Ckdn1a :
cyclin-dependent kinase inhibitor 1a
coding DNA sequence
- Col4a1 :
collagen type IV alpha 1
- Cox2 :
- Cybrd1 :
cytochrome b reductase 1
DNA binding domain
- Dhh :
- Egr1 :
early growth response 1
fetal bovine serum
false discovery rate
- Fgfr2 :
fibroblast growth factor receptor 2
- Fst :
- Gja1 :
gap junction protein alpha 1
- Ifnb1 :
interferon beta 1, fibroblast
- Ifr7 :
interferon regulatory factor 7
- Il-1β :
- Inhbb :
inhibin beta B
Ingenuity® Pathway Analyzer
ligand binding domain
- Myc :
no template control
polyacrylamide gel electrophoresis
phosphate buffered saline
- Pmepa1 :
prostate transmembrane protein, androgen induced 1
probabilistic principal component analysis
- Ptsg2 :
prostaglandin-endoperoxide synthase 2
- Rarg :
retinoid acid receptor, gamma
- Rplp :
ribosomal protein l 16
reverse transcription polymerase chain reaction
WL3 mouse Sertoli cell lines
- Smad 1 :
SMAD familiy member 1
- Sox9 :
- Sry :
sex determining reagion on Y chromosome
- Steap2 :
six-transmembrane epithelial antigen of the prostate 2
Tris-Buffered Saline and Tween 20
- Tf :
- Tgfb1i1 :
transforming growth factor beta 1 induced transcript 1
- Tnfrsf1a :
tumor necrosis factor receptor superfamily, member 1a
- Ubc :
Walker WH. Molecular mechanisms of testosterone action in spermatogenesis. Steroids. 2009;74(7):602–7.
Gelmann EP. Molecular biology of the androgen receptor. J Clin Oncol. 2002;20(13):3001–15.
Yong EL, Loy CJ, Sim KS. Androgen receptor gene and male infertility. Hum Reprod Update. 2003;9(1):1–7.
Coffey K, Robson CN. Regulation of the androgen receptor by post-translational modifications. J Endocrinol. 2012;215(2):221–37.
Bergh A, Damber JE. Immunohistochemical demonstration of androgen receptors on testicular blood vessels. Int J Androl. 1992;15(5):425–34.
Sertoli E. Dell’ esistenza di particolari cellule ramificante nei cunalicoli seminiferi del testicoli umano. Morgagni. 1865;7:31–40.
Lee H, Chang C. Recent advances in androgen receptor action. Cell Mol Life Sci. 2003;60(8):1613–22.
Willems A, Batlouni SR, Esnal A, Swinnen JV, Saunders PTK, Sharpe RM, et al. Selective ablation of the androgen receptor in mouse sertoli cells affects sertoli cell maturation, barrier formation and cytoskeletal development. PLoS One. 2010;5(11):e14168.
Ortmann J, Prifti S, Bohlmann MK, Rehberger-Schneider S, Strowitzki T, Rabe T. Testosterone and 5 alpha-dihydrotestosterone inhibit in vitro growth of human breast cancer cell lines. Gynecol Endocrinol. 2002;16(2):113–20.
Jaroenporn S, Furuta C, Nagaoka K, Watanabe G, Taya K. Comparative effects of prolactin versus ACTH, estradiol, progesterone, testosterone, and dihydrotestosterone on cortisol release and proliferation of the adrenocortical carcinoma cell line H295R. Endocrinology. 2008;33(2):205–9.
Pronsato L, Boland R, Milanesi L. Testosterone exerts antiapoptotic effects against H2O2 in C2C12 skeletal muscle cells through the apoptotic intrinsic pathway. J Endocrinol. 2012;212(3):371–81.
Chen G, Li S, Dong X, Bai Y, Chen A, Yang S, et al. Investigation of testosterone, androstenone, and estradiol metabolism in HepG2 cells and primary culture pig hepatocytes and their effects on 17βHSD7 gene expression. PLoS One. 2012;7(12):e52255.
Yeh S, Hu Y, Wang P, Xie C, Xu Q, Tsai M, et al. Abnormal mammary gland development and growth retardation in female mice and MCF7 breast cancer cells lacking androgen receptor. J Exp Med. 2003;198(12):1899–908.
Szelei J, Jimenez J, Soto AM, Luizzi MF, Sonnenschein C. Androgen-induced inhibition of proliferation in human breast cancer MCF7 cells transfected with androgen receptor. Endocrinology. 1997;138(4):1406–12.
Yuan S, Trachtenberg J, Mills GB, Brown TJ, Xu F, Keating A. Androgen-induced inhibition of cell proliferation in an androgen-insensitive prostate cancer cell line (PC-3) transfected with a human androgen receptor complementary DNA. Cancer Res. 1993;53(6):1304–11.
Xiao Y, Karnati S, Qian G, Nenicu A, Fan W, Tchatalbachev S, et al. Cre-mediated stress affects sirtuin expression levels, peroxisome biogenesis and metabolism, antioxidant and proinflammatory signaling pathways. PLoS One. 2012;7(7):e41097.
Li Y, Chan SC, Brand LJ, Hwang TH, Silverstein KAT, Dehm SM. Androgen receptor splice variants mediate enzalutamide resistance in castration-resistant prostate cancer cell lines. Cancer Res. 2013;73(2):483–9.
Jiang C, Hall SJ, Boekelheide K. Development and characterization of a prepubertal rat Sertoli cell line, 93RS2. J Androl. 1997;18(4):393–9.
GEO Accession viewer. http://www.ncbi.nlm.nih.gov/gds/?term=GSE57653 Accession 15 May 2014.
Bitgood MJ, Shen L, McMahon AP. Sertoli cell signaling by Desert hedgehog regulates the male germline. Curr Biol. 1996;6(3):298–304.
Weider K, Bergmann M, Brehm R. Connexin 43: its regulatory role in testicular junction dynamics and spermatogenesis. Histol Histopathol. 2011;26(10):1343–52.
Sharpe RM, McKinnell C, Kivlin C, Fisher JS. Proliferation and functional maturation of Sertoli cells, and their relevance to disorders of testis function in adulthood. Reproduction. 2003;125(6):769–84.
Skinner MK, Schlitz SM, Anthony CT. Regulation of Sertoli cell differentiated function: testicular transferrin and androgen-binding protein expression. Endocrinology. 1989;124(6):3015–24.
Barakat B, Itman C, Mendis SH, Loveland KL. Activins and inhibins in mammalian testis development: new models, new insights. Mol Cell Endocrinol. 2012;359(1–2):66–77.
Ashe PC, Berry MD. Apoptotic signaling cascades. Prog Neuropsychopharmacol Biol Psychiatry. 2003;27(2):199–214.
Jacobsen L, Calvin S, Lobenhofer E. Transcriptional effects of transfection: the potential for misinterpretation of gene expression data generated from transiently transfected cells. Biotechniques. 2009;47(1):617–24.
Li F, Yamaguchi K, Okada K, Matsushita K, Enatsu N, Chiba K, et al. Efficient transfection of DNA into primarily cultured rat sertoli cells by electroporation. Biol Reprod. 2013;88(3):61.
Leichtmann-Bardoogo Y, Cohen LA, Weiss A, Marohn B, Schubert S, Meinhardt A, et al. Compartmentalization and regulation of iron metabolism proteins protect male germ cells from iron overload. Am J Physiol Endocrinol Metab. 2012;302(12):E1519–30.
Skinner MK, Griswold MD. Sertoli cells synthesize and secrete transferrin-like protein. J Biol Chem. 1980;255(20):9523–5.
Matzkin ME, Mayerhofer A, Rossi SP, Gonzalez B, Gonzalez CR, Gonzalez-Calvar SI, et al. Cyclooxygenase-2 in testes of infertile men: evidence for the induction of prostaglandin synthesis by interleukin-1β. Fertil Steril. 2010;94(5):1933–6.
Yamaguchi K, Ishikawa T, Kondo Y, Fujisawa M. Cisplatin regulates Sertoli cell expression of transferrin and interleukins. Mol Cell Endocrinol. 2008;283(1–2):68–75.
Sledz CA, Holko M, de Veer MJ, Silverman RH, Williams BR. Activation of the interferon system by short-interfering RNAs. Nat Cell Biol. 2003;5(9):834–9.
Bitgood MJ, McMahon AP. Hedgehog and Bmp genes are coexpressed at many diverse sites of cell-cell interaction in the mouse embryo. Dev Biol. 1995;172(1):126–38.
Clark AM, Garland KK, Russell LD. Desert hedgehog (Dhh) gene is required in the mouse testis for formation of adult-type Leydig cells and normal development of peritubular cells and seminiferous tubules. Biol Reprod. 2000;63(6):1825–38.
Mäkelä J, Saario V, Bourguiba-Hachemi S, Nurmio M, Jahnukainen K, Parvinen M, et al. Hedgehog signalling promotes germ cell survival in the rat testis. Reproduction. 2011;142(5):711–21.
Schmahl J, Kim Y, Colvin JS, Ornitz DM, Capel B. Fgf9 induces proliferation and nuclear localization of FGFR2 in Sertoli precursors during male sex determination. Development. 2004;131(15):3627–36.
Bagheri-Fam S, Sim H, Bernard P, Jayakody I, Taketo MM, Scherer G, et al. Loss of Fgfr2 leads to partial XY sex reversal. Dev Biol. 2008;314(1):71–83.
Valeri C, Schteingart HF, Rey RA. The prepubertal testis: biomarkers and functions. Curr Opin Endocrinol Diabetes Obes. 2013;20(3):224–33.
Pfaff T, Rhodes J, Bergmann M, Weinbauer GF. Inhibin B as a marker of sertoli cell damage and spermatogenic disturbance in the rat. Birth Defects Res B Dev Reprod Toxicol. 2013;98(1):91–103.
Mueller SO, Korach KS. Immortalized testis cell lines from estrogen receptor (ER) alpha knock-out and wild-type mice expressing functional ERalpha or ERbeta. J Androl. 2001;22(4):652–64.
Sneddon SF, Walther N, Saunders PTK. Expression of androgen and estrogen receptors in sertoli cells: studies using the mouse SK11 cell line. Endocrinology. 2005;146(12):5304–12.
Konrad L, Munir Keilani M, Cordes A, Völck-Badouin E, Laible L, Albrecht M, et al. Rat Sertoli cells express epithelial but also mesenchymal genes after immortalization with SV40. Biochim Biophys Acta. 2005;1722(1):6–14.
Weider K, Bergmann M, Giese S, Guillou F, Failing K, Brehm R. Altered differentiation and clustering of Sertoli cells in transgenic mice showing a Sertoli cell specific knockout of the connexin 43 gene. Differentiation. 2011;82(1):38–49.
Chomczynski P. A reagent for the single-step simultaneous isolation of RNA, DNA and proteins from cell and tissue samples. BioTechniques. 1993;15(3):532–4.
Fietz D, Geyer J, Kliesch S, Gromoll J, Bergmann M. Evaluation of CAG repeat length of androgen receptor expressing cells in human testes showing different pictures of spermatogenic impairment. Histochem Cell Biol. 2011;136(6):689–97.
Bolstad BM, Irizarry RA, Astrand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003;19(2):185–93.
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, et al. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 2003;4(5):P3.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
DF and MM drafted the manuscript and performed quantitative RT-PCR, microarray analysis and statistical analysis, respectively. DL performed screening for cell lines, transfection, and immunofluorescence staining. KL was responsible for study design and supervision of the transfection procedure; furthermore he was involved in cell line acquisition and cell culture experiments. GJ supervised transfection procedure as well and generously provided the AR-GFP construct. KS provided human testis material for positive control purposes. Both GJ and KS critically revised the manuscript. CT and HH participated in the design of microarray experiments and HH also helped to draft the manuscript. MB was responsible for study design, supervision and critically revised the manuscript. All authors read and approved the final manuscript.
The authors want to thank J. Vogelsberg and D. Zoltan for their skilful technical assistance. Cell lines were generously provided by Profs. Korach, Boekelheide and Gromoll. The presented study was funded by German Research Foundation (DFG KFO181, BE1061/7-1).
The authors declare that they have no competing interests.
D. Fietz, M. Markmann, D. Lang contributed equally to the study
About this article
Cite this article
Fietz, D., Markmann, M., Lang, D. et al. Transfection of Sertoli cells with androgen receptor alters gene expression without androgen stimulation. BMC Molecular Biol 16, 23 (2015). https://doi.org/10.1186/s12867-015-0051-7
- Gene expression analysis
- Androgen receptor
- Sertoli cells