Aberrant gene expression by Sertoli cells in infertile men with Sertoli cell-only syndrome

Abstract

Sertoli cell-only (SCO) syndrome is a severe form of human male infertility seemingly characterized by the lack all spermatogenic cells. However, tubules of some SCO testes contain small patches of active spermatogenesis and thus spermatogonial stem cells. We hypothesized that these stem cells cannot replicate and seed spermatogenesis in barren areas of tubule because as-of-yet unrecognized deficits in Sertoli cell gene expression disable most stem cell niches. Performing the first thorough comparison of the transcriptomes of human testes exhibiting complete spermatogenesis with the transcriptomes of testes with SCO syndrome, we defined transcripts that are both predominantly expressed by Sertoli cells and expressed at aberrant levels in SCO testes. Some of these transcripts encode proteins required for the proper assembly of adherent and gap junctions at sites of contact with other cells, including spermatogonial stem cells (SSCs). Other transcripts encode GDNF, FGF8 and BMP4, known regulators of mouse SSCs. Thus, most SCO Sertoli cells can neither organize junctions at normal sites of cell-cell contact nor stimulate SSCs with adequate levels of growth factors. We propose that the critical deficits in Sertoli cell gene expression we have identified contribute to the inability of spermatogonial stem cells within small patches of spermatogenesis in some SCO testes to seed spermatogenesis to adjacent areas of tubule that are barren of spermatogenesis. Furthermore, we predict that one or more of these deficits in gene expression are primary causes of human SCO syndrome.

About the Authors:

Darius A. Paduch

Roles Conceptualization, Investigation, Writing – review & editing

Affiliations Department of Urology, Weill Cornell Medical College, New York, NY, United States of America, Department of Molecular Biology and Genetics, Cornell University, Ithaca, NY, United States of America

Stephanie Hilz

Roles Formal analysis

Affiliations Department of Neurological Surgery, University of California, San Francisco, California, United States of America, Genomic Analysis and Sequencing Core Facility, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland, United States of America

Andrew Grimson

Roles Formal analysis, Writing – review & editing

Affiliation: Department of Neurological Surgery, University of California, San Francisco, California, United States of America

Peter N. Schlegel

Roles Funding acquisition, Writing – review & editing

Affiliation: Department of Urology, Weill Cornell Medical College, New York, NY, United States of America

Anne E. Jedlicka

Roles Formal analysis

Affiliation: Department of Biochemistry and Molecular Biology, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland, United States of America

William W. Wright

Roles Formal analysis, Methodology, Writing – original draft

* E-mail: [email protected]

Affiliation: Consulting Research Services, Inc, North Bergen, N.J., United States of America

ORCID logo http://orcid.org/0000-0002-2413-8886

Introduction

Infertility is a problem that besets approximately 15% of couples of childbearing age, with men being the sole cause of the couple’s infertility approximately one-third of the time [1]. A severe form of human male infertility is characterized histologically by the apparent lack of all spermatogenic cells in almost all seminiferous tubules, a condition called Sertoli cell-only (SCO) syndrome. While tubules without germ cells are markedly atrophied, small segments of seminiferous tubules within some SCO testes are dilated, exhibit active spermatogenesis and produce sperm [2]. As the cellular foundation of spermatogenesis is the spermatogonial stem cell (SSC), the question arises as to why self-renewing replication of the SSCs in these segments does not generate new stem cells that can colonize adjacent germ cell-deficient areas of tubule and thereby seed spermatogenesis. Normally, SSCs have this capacity [3, 4]. We propose that one of the reasons that this does not occur is that the somatic cells associated with those barren areas of tubule do not form a functional stem cell niche. Thus, identifying deficits in gene expression by those somatic cells may provide new insights into the mechanistic basis for this form of male infertility and serve as a potential first step toward developing new therapies for men whose apparently SCO testes contain some SSCs. This study focuses primarily on Sertoli cells, the only somatic cells in direct physical contact with SSCs, and an important source of growth factors that regulate the replication and differentiation of SSCs and their immediate progeny, transit amplifying progenitor spermatogonia [5, 6]. Thus, Sertoli cells within a fertile testis are both the architects of the stem cell niche and the source of important growth factors that act within it. It follows that aberrant gene expression by Sertoli cells may disable some stem cell niches, and, thus, be one cause of Sertoli cell-only syndrome.

The analyses presented herein are an important step towards the full characterization of aberrant gene expression by Sertoli cells within human SCO testes. We collected testicular biopsies from patients who were undergoing micro-testicular sperm extraction, and then performed the first thorough comparison of the transcriptomes of testis with complete spermatogenesis (i.e. spermatogenesis was either quantitatively or qualitatively normal) with the transcriptomes of testes with SCO syndrome. We then focused on transcripts expressed primarily by Sertoli cells in a fertile human testis. We refer to these as Sertoli cell signature transcripts.

Because the transcriptome of normal human Sertoli cells has yet to be analyzed directly, we used a two-step bioinformatic approach to define human Sertoli cell signature transcripts. This approach was founded on the following considerations. First, Sertoli cells of all mammalian species fulfill the same essential functions, which are determined in large part by their transcriptomes. Second, the transcriptomes of adult rat Sertoli cells and other testicular cell types are known and thus, rat Sertoli cell signatures transcripts can be identified directly [7, 8]. Third, that rat genes and their human orthologues predominantly encode proteins with the same fundamental functions. Thus, we anticipated that the compendium of rat and human Sertoli cell signature transcripts would be largely congruent. Fourth, the transcriptomes of a human Leydig and spermatogenic cells and numbers of these cells in an adult testis have been defined [9–12]. Thus, Leydig and spermatogenic cell transcriptomes could be used to estimate the contribution of each cell type to the overall testis content of a given transcript, thereby allowing us to winnow and refine our list of human Sertoli cell signature transcripts.

A goal of our studies was to identify signature transcripts whose aberrant expression might explain important characteristics of human SCO testes and their dysfunctional niches. The first characteristic is that the SCO Sertoli cells fail to form highly organized and dynamic junctions with each other and with spermatogenic cells [13, 14]. In a fertile testis, junctions between Sertoli cells create the blood-testis barrier and thus the intratubular environment required by spermatogenesis [15]. Sertoli cells also form junctions with SSCs and other spermatogenic cells, and these junctions provide the positional information needed for the stem cells to undergo the asymmetric divisions that produce one new stem cell and one progenitor spermatogonium [16, 17]. The disorganization of cell-cell junctions in SCO testes implies that such positional information is diminished.

A second characteristic of SCO testes is the inability of its few SSCs to expand in number, fill empty niches and seed new areas of tubule with spermatogenesis. We hypothesized that one reason that SSC expansion does not occur is inadequate production of essential paracrine regulators of these stem cells by Sertoli cells. Mouse Sertoli cells produce a number of such factors, including glial cell line-derived neurotrophic factor (GDNF), which stimulate replication and migration of these stem cells and suppresses their differentiation [6, 18]. We previously demonstrated that Sertoli cells in human SCO testes produce markedly less GDNF than their counterparts in fertile testes [19]. However, as multiple paracrine factors act in concert to regulate SSCs, we sought to determine if Sertoli cells in human SCO testes mis-express other important paracrine regulators of these stem cells.

Our analyses provide the first detailed comparison of the transcriptomes of human testes with complete spermatogenesis to transcriptomes of testes with SCO syndrome and identify deficits in gene expression by Sertoli cells that potentially explain the two characteristics of the human SCO testis noted above. Importantly, these explanations are strengthened by published studies on rodents that document the importance of specific genes to the functions of niches and the SSCs that they shelter.

Materials and methods

Human subjects

Testis biopsies were obtained from 11 patients, who were undergoing micro-testicular sperm extraction as part of treatment for infertility [2, 20]. All patients gave written informed consent for the collection of the biopsies, and The Institutional Review Board of Weill Medical College of Cornell University approved both the collection of tissues and their analyses. Table 1 describes the age, weight and height of each patient as well as their peripheral blood cytogenetics (determination of karyotype and a screen for Y chromosome microdeletions), success of sperm retrieval and fertilization by ICSI, serum FSH levels, and testicular histology. Each biopsy was separated into two pieces; the first was used for histological analysis; the second (~15–40 mg) was immediately snap frozen in liquid nitrogen and subsequently used for transcriptional analysis [20]. Tissue for histological analysis was fixed in 4% buffered paraformaldehyde, embedded in paraffin and 5-micron hematoxylin and eosin-stained stained histological sections prepared. Three sections were obtained from three different levels of each biopsy. Dr. Paduch and a member of the Weill-Cornell Department of Pathology independently evaluated the histology of each biopsy, which they characterized as: normal, no substantial abnormalities; exhibiting mild or moderate hypospermatogenesis, spermatogenesis is qualitatively normal, but there are reduced numbers of germ cells in tubule cross-sections; maturation arrest, a failure of spermatogenesis during meiosis or spermiogenesis; or SCO, exhibiting Sertoli cell-only syndrome. S1 Fig. shows representative micrographs of a normal testis and of testes exhibiting either mild or moderate hypospermatogenesis. We have previously published a representative photomicrograph of a section of an SCO testis, which appears devoid of all spermatogenic cells [19]. In our current study, data from testes with normal spermatogenesis or with hypospermatogenesis are combined under a single category, complete spermatogenesis (CS). When there was significant variation in the histology of tubules within a biopsy, we estimated the percentages of tubules exhibiting each of the four characteristics described above. Dr. Paduch and the pathologist concurred on the histology of the 11 testis biopsies that are the subject of this study.

[Figure omitted. See PDF.]

Table 1. Description of donors.

https://doi.org/10.1371/journal.pone.0216586.t001

Transcriptional analysis of testis biopsies

Total RNA was isolated from frozen biopsies by use of the miRCURY RNA Isolation Kit (EXIQON, Wolbern, MA). Frozen tissue was homogenized in lysis buffer with TissueRuptor (Quigen, Valencia, CA), and total RNA isolated following manufacturer’s instructions. RNA purity and integrity were evaluated with an Agilent Bioanalyzer 2100 (Agilent Technologies, Lexington, MA).

Non-stranded RNAseq libraries were prepared (TruSeq, Illumina; 15 cycles of amplification) and sequenced on an Illumina HiSeq 2500 (50 bp; single read; average library depth was 17 million reads). Reads were aligned to genes and numbers of transcripts quantified using the Burrows-Wheller Alignment tool and HTSeq [21, 22]. The Pre-normalized data were filtered using a threshold of > 0.4 CPMs (counts per million mapped reads) in at least 2 of the libraries prepared from testes with complete spermatogenesis (https://support.bioconductor.org/p/60205/) and purged of genes that did not encode mRNAs. edgeR [23] was used to calculate normalized gene expression from raw counts of transcripts for human genes (iGenomes hg19 gene annotation file). CPM values for each library then were normalized, log2 transformed, and significant differences in expression between the two sets of testes defined as a false discovery rate (FDR) < 0.05.

To confirm that the threshold of 0.4 CPM noted above was a good choice for a cutoff in our whole transcriptome differential expression analysis, we examined the relationship between transcript expression and logFC (MA plot), as well as FDR p-value (volcano plot), and observed that there was no over-abundance of low-expression genes with apparent outlier logFC values that were being identified as differentially expressed. R script used for analysis of the sequencing data is available online (https://github.com/SRHilz/SCO).

As somatic cells are enriched in SCO testes, we initially evaluated expression of well-characterized transcripts after data were normalized for beta actin mRNA. Two tailed T-tests were used to compare beta-actin-normalized CPM values in CS and SCO testes. These comparisons were performed with GraphPad Prism 6.0 (La Jolla, CA) and significant differences defined as p≤0.005.

Defining signature transcripts of human Sertoli cells

The very limited size of the testis biopsies precluded isolating Sertoli cells prior to sequencing. To circumvent this limitation, we first assessed the entire testis transcriptome of each sample, and then focused our analyses on transcripts that were predominantly expressed by human Sertoli cells and thus were signature transcripts for this cell-type. We developed a two-step screen to define such transcripts, which is described in detail S1 Text. Supplementary Methods. The first step was based on both the knowledge that Sertoli cells of all mammalian species fulfill the same essential functions and on our expectation that orthologous genes encode analogous functions. Thus, we used the published transcriptomes of rat Sertoli cells, Leydig cells, spermatogonia, pachytene spermatocytes and round spermatids to identify transcripts that were expressed at least 4-fold higher by rat Sertoli cells than by any of the other 4 cell types [8, 24]. We then identified orthologues of these transcripts in human testes.

The second step estimated the level of expression of orthologous transcripts by human Sertoli cells and identified those transcripts that met the criteria for signature transcripts (described below). To generate this estimate, we combined the human testis transcriptome defined by our analyses with published transcriptomes of human Leydig cells and spermatogenic cells and the relative numbers of each cell type, including Sertoli cells, in a human testis [9–12]. We calculated the contribution of Leydig cells and of each germ cell type to the amount of each transcript in the total testis transcriptome. Then we estimated the expression of a transcript by Sertoli cells with complete spermatogenesis by subtracting the total amount of a given transcript expressed by Leydig cells plus spermatogenic cells from the total testis transcriptome. We defined signature transcripts of Sertoli as meeting the following criteria: (a) Sertoli cells are the source of at least 60% of that transcript in the testis, (b) the expression of the transcript by Sertoli cells was at least twice that of Leydig cells, and (c) expression by Sertoli cells was greater than expression by any germ cell. Prior to comparing expression of Sertoli cells signature transcripts in testes with complete spermatogenesis with testes with SCO syndrome, CPM values for the signature transcripts in each sample were normalized to the same median value and data statistically reanalyzed as described above for the total testis transcriptomes

Finally, as signature transcripts are essential to Sertoli cell differentiated function and as these functions are important to fertility of all male mammals, we tested the prediction that most human Sertoli cell signature transcripts were also expressed by mature mouse Sertoli cells, which have been thoroughly studied by single-cell RNA sequencing [25]. To this end, we searched for orthologues of the human Sertoli cell signature transcripts in the mature mouse Sertoli cell transcriptome.

Identification of FGF8 as a product of human Sertoli cells and measurement of FGF8 mRNA by real time PCR

Our analyses identified FGF8 as a product of Sertoli cells that was expressed at remarkably reduced levels in the SCO testis. As FGF8 is a regulator of mouse SSCs [26] we used FACs analysis, to confirm that this growth factor is expressed primarily by Sertoli cells in the human testis, and real time RT-PCR, to corroborate that FGF8 mRNA expression is marked reduced in SCO testes. Details of these two procedures are described in Supplemental methods.

Results

Comparison of the transcriptomes of testes with complete spermatogenesis with the transcriptomes of testes with SCO syndrome

Table 1 describes the donors of the biopsies, and Table 2 describes the histology of the biopsies. Testes from these patients were assigned to one of two groups based on histology. Spermatogenesis was complete in four of the biopsies (C1 to C4). Sertoli cell-only syndrome was evident in all tubules of 6 of the seven other testes. The histology of the seventh SCO testis, SCO11, was mixed. Eighty percent of the tubules were SCO and 20% exhibited severe hypospermatogenesis. The patients in the two groups were similar in age, weight and height, but patients with SCO syndrome had elevated serum FSH levels. All patients but one had a normal karyotype and a normal Y chromosome. The donor of sample SCO11 was diagnosed with Kleinfelter syndrome. Sperm were retrieved from one or more segments of dilated tubules and fertilization by intracytoplasmic sperm injection was achieved for all patients. Dilated tubule segments in testes with SCO syndrome could not be collected for RNA seq analysis, since all were required for sperm retrieval and subsequent ICSI.

[Figure omitted. See PDF.]

Table 2. Histology of the biopsies obtained from donors.

https://doi.org/10.1371/journal.pone.0216586.t002

RNA was isolated from each biopsy and the transcriptomes determined by RNAseq. These data have been uploaded to the NCBI database of Genotypes and Phenotypes (Study Name: NICHD/SCO_Syndrome; Accession Number phs001777.vi.p1). A total of 17,208 distinct mRNAs were identified in testes, and of those 4,748 and 5,274 transcripts had significantly higher or lower expression, respectively, in SCO testes. These differentially expressed transcripts sorted into 10 clusters based on expression levels and subsequent Gene Ontology (GO) enrichment analysis (Fig 1A). As anticipated, transcripts with GO terms “Spermatogenesis”.and “Reproduction” were substantially less abundant in SCO testes. However, genes with GO terms “G-protein coupled receptor activity” were enriched in the transcriptomes of SCO testes. A likely explanation for this enrichment was that the somatic cells in the testis express many G-protein coupled receptors and somatic cells are significantly enriched in the SCO testis. In fact, many transcripts expressed primarily by somatic cells, including Sertoli cells, were enriched approximately 3-fold in SCO testes; so, too, was beta actin mRNA. Therefore, to initially compare expression of characterized- somatic cell and germ cell transcripts in the two patient populations, we normalized CPMs (counts per million map reads) for each transcript to CPMs of beta actin mRNA. Fig 1B shows data for vimentin (VIM), smooth muscle actin (ACTA2) and inhibin alpha (INHA), which are considered markers of Sertoli cells, peritubular myoid cells and Leydig cells, respectively [19, 27, 28]. Normalized expressions of these three transcripts were similar in human testes with complete spermatogenesis and in testes exhibiting SCO syndrome. Testes in both groups also expressed similar levels of transcripts that encode receptors for the three hormones essential to male fertility, FSH, testosterone and LH (Fig 1B). Thus, some genes essential for male fertility are expressed normally by somatic cells in SCO testes.

[Figure omitted. See PDF.]

Fig 1. Comparison of the transcriptomes of testes with complete spermatogenesis with the transcriptomes of testes with SCO syndrome.

(A) K means cluster analysis comparing the transcriptomes of testes with complete spermatogenesis with testes with SCO syndrome. Genes were clustered under 5 Gene Ontology terms; one or two representative terms per cluster are listed here. (B) The abundance in of a few well-characterized transcripts considered markers of three testicular somatic cells: Sertoli cells (VIM), peritubular myoid cells (ACTA2), and Leydig cells (INHA). Also shown are transcripts encoding receptors for three hormones required for normal testis function: FSHR, AR, and LHCGR. (C) Abundance of transcripts expressed by human SSCs (GFRA1, RET, ZBTB16, UCHL1, ENO2), by spermatocytes (SYCP2, SYCP3) and by round spermatids (PRM1, PRM2). Data in panels B and C are presented as CPM values in the total testis transcriptome divided by the CPM value for beta actin in the same sample. Asterisks over a pair of bars identify significant differences between testes with complete spermatogenesis and SCO testes (p<0.005).

https://doi.org/10.1371/journal.pone.0216586.g001

We previously reported that the expression of a marker of human SSCs, GFRA1 mRNA, was significantly reduced to 10% of its expression in testes with complete spermatogenesis, while expression of a marker of spermatocytes and spermatids, DDX4 mRNA was reduced to 0.1% [19]. To confirm and extend those results we compared beta actin-normalized expression of 9 different germ cell transcripts in the two groups of testes (Fig 1C). Consistent with results of the study noted above, expression of GFRA1 mRNA in SCO testes was 8% of its expression in testes with complete spermatogenesis. Expressions of four other markers of SSCs (Ret, ZBTB16, UCLH1 and ENO2) were also reduced to 19% to 30% of the levels in testes with complete spermatogenesis. Transcripts expressed by spermatocytes (SYCP2 and SYCP3) or by spermatids (PRM1 and PRM2) were reduced to an even greater extent, approximately 1% and 0.1%, respectively, of what was measured in testes with complete spermatogenesis. The concentrations of germ cell-expressed transcripts in individual biopsies were consistent with their histology (S2 Fig). All of these transcripts were markedly higher in biopsies with complete spermatogenesis than in those exhibiting SCO syndrome. Furthermore, the concentrations of these transcripts were higher in testes with normal spermatogenesis (C1, C2) than in testes with hypospermatogenesis (C3, C4) (S2 Fig). Finally, our analyses demonstrated that in SCO testes there is a much greater decrease in the concentrations of the four transcripts expressed by meiotic or haploid germ cells than in the concentrations of the five of transcripts expressed by SSCs (Fig 1 and S2 Fig). This suggests that SSCs are present in some apparently barren tubules in SCO testes, but produce few differentiated progeny.

Identification and categorization of human Sertoli cell signature transcripts

In order to define Sertoli cell signature transcripts, we next subjected published data on the transcriptomes of isolated testicular cells to the two-step screen that is described in Methods. The first step used the published transcriptomes of adult rat Sertoli cells, Leydig cells and spermatogenic cells to identify 539 transcripts that were expressed at least 4-fold higher by Sertoli cells than any other cell types in the rat testis [8, 24] (S1 Table). An additional 84 transcripts were added to the tentative Sertoli cell list, even though they were not detected by the arrays, but are known to encode components of the blood-testis barrier or are growth factors for receptors expressed by human SSCs [29, 30]. The second step estimated the expression of these 623 transcripts by human Sertoli cells and compared these estimates to the expression of those transcripts by human Leydig cells and/or by spermatogenic cells [9–12]. We then used three criteria to define a Sertoli cell signature transcript: (a) Sertoli cells were the source of at least 60% of that transcript in the testis, (b) the expression of the transcript by Sertoli cells was at least twice that of Leydig cells, and (c) expression by Sertoli cells was greater than expression by any germ cell. Two hundred forty three transcripts met these criteria. The gene symbols for these transcripts and the calculations that identified them as Sertoli cell signature transcripts are presented in S2 Table.

As discussed in Methods, we predicted that given the important and conserved functions of Sertoli cells, the transcriptomes of mouse Sertoli cells should contain orthologues of most human Sertoli cell signature transcripts. It was therefore of note that 211 of the 243 human Sertoli cells signature transcripts were detected when mature mouse Sertoli cells were subjected to single cell sequencing (S3 Table) [25]. There is precedence for a small difference in the transcriptomes of human and rodent Sertoli cells. In rodent testes Sertoli cells express the sex hormone binding globulin gene [24, 31]; however, in the human testis this gene is not expressed by Sertoli cells but rather by spermatogenic cells [32]. Nonetheless, after we made an overall comparison of the expression of Sertoli cell signature transcripts in testes with complete spermatogenesis with testes with SCO syndrome, we focused our analyses on the human signature transcripts that are also expressed by mature mouse Sertoli cells.

Our next step was to identify Sertoli cell signature transcripts whose expressions differed between testes with complete spermatogenesis and testes with SCO syndrome. Expression of these transcripts were normalized to the same median CPM value and compared statistically (S3 Table). Principal component analyses demonstrated that overall, the expression of signature Sertoli cell transcripts differed between the two groups of testes (Fig 2A). Seventy-five and 58 Sertoli cell signature transcripts were significantly more or less abundant, respectively, in SCO testes (Fig 2B). It is noteworthy that four transcripts (WT1, DHH, DKK3 and PRND) that were expressed at higher levels in human SCO testes are essential for fertility of male mice [33–36] (S3 Fig). The remaining 111 signature transcripts were expressed at similar levels in testes with complete spermatogenesis and testes with SCO syndrome. Gene ontology analysis identified some of the less abundant transcripts as associated with cell development or signaling, while more abundant transcripts were associated with actin filaments and non-membrane bound organelles. Thus, this rigorous transcriptional analysis identified subsets of transcripts that could be analyzed further, considering the known characteristics of the SCO testis and its dysfunctional spermatogonial stem cell niche.

[Figure omitted. See PDF.]

Fig 2. The expression of Sertoli cell signature transcripts in testes with complete spermatogenesis and in testes with SCO syndrome.

(A) PCA analysis of the expression of Sertoli cell signature transcripts in individual testes. Data were obtained from 4 testes with complete spermatogenesis (C1-4, green symbols) and seven SCO testes (SCO5, SCO7-12, red symbols) (B). A comparison of the mean level of expression of each Sertoli cell signature transcript in testes with complete spermatogenesis and in testes with SCO syndrome. Transcripts that were expressed in SCO testes at significantly higher levels (FDR<0.05), the same level or significantly lower levels are identified by blue, white and black dots, respectively.

https://doi.org/10.1371/journal.pone.0216586.g002

Expression of genes encoding cell surface adhesion proteins and the gap junction protein, Connexin 26

Sertoli cells not only provide the architectural framework for the stem cell niche but also for the entire seminiferous epithelium. In normal testes this epithelium is divided into two compartments, a basal compartment in which spermatogonia and early preleptotene spermatocytes reside and an adluminal compartment that contains more mature spermatogenic cells. The blood-testis barrier separates these two compartments, which are formed by two different types of junctions on apposing Sertoli cell surfaces: adhesive junctions (tight junctions, ectoplasmic specializations (ES)) and gap junctions (communicating junctions), which allow passage of small second messenger molecules, such as cAMP and Ca+2, between cells [37]. The same junctions also mediate physical and physiological interactions between Sertoli cells and germ cells, including spermatogonia [38]. It is, therefore, of note that human SCO testes lack a functional blood-testis barrier and that some proteins that mediate adhesion or communication between two adjacent Sertoli cells and between Sertoli cells and spermatogenic cells are either mislocated or undetectable in human SCO testes [13, 14].

As in other species, human Sertoli cells transcriptionally express multiple integral membrane proteins that form tight junctions and ectoplasmic specializations (ES) [24]. Sertoli cell signature transcripts encode claudins or junction adhesion molecules, adhesive proteins in tight junctions, and cadherins or nectins, the adhesive proteins in ectoplasmic specializations (S2 Table). Fig 3A presents data for the two most abundant transcripts encoding adhesive proteins in tight junctions, Claudin 11 (CLDN 11) and junction adhesion molecule 3 (JAM3), and the most abundant transcripts for the adhesive proteins in ectoplasmic specializations, nectin cell adhesion molecule 2 (NECTIN 2) and cadherin 2 (CDH2). CLDN11 and Nectin 2 were expressed at higher levels in SCO testes; Jam3 and CHD2 were expressed at similar levels in testes with complete spermatogenesis and testes with SCO syndrome. In fact, 8 of the 13 transcripts that encode other adhesive proteins were either expressed at the same or higher levels in SCO testes (S2 Table). Thus, failure of adhesive junction formation is not due to a lack of expression of transcripts that encode the integral membrane adhesive proteins at cell-cell junctions.

[Figure omitted. See PDF.]

Fig 3. Expression of Sertoli cell signature transcripts encoding proteins that form adhesive or gap junctions or that organize and polarize the Sertoli cell plasma membrane.

(A) CLDN11 and JAM3 are tight junction adhesive proteins. Nectin2 and CDH2 are adhesive proteins of ectoplasmic specializations. GJB2 encodes the gap junction protein, Connexin 26. (B) Transcripts that encode proteins that organize or polarize the plasma membrane. CRB2 is a cell polarity component protein that defines the site of separation of the apical and basal-lateral plasma membrane domains. TJP3 is an adaptor protein that links the adhesive proteins of a tight junction to the actin cytoskeleton. SPTB, SPTPN4 and SPTPN5 are scaffolding proteins that concentrate beneath the plasma membrane, seed formation of membrane microdomains, and link those domains to the actin cytoskeleton. All data in this figure are expressed as mean + SEM of CPMs. The presence of an asterisk over a pair of bars indicates a significant difference between testes with complete spermatogenesis and testes with SCO syndrome (FDR≤0.05).

https://doi.org/10.1371/journal.pone.0216586.g003

We note, however, that not all genes for proteins that mediate cell-cell interactions were expressed at normal or elevated levels in SCO testes. The expression of GJB2, which encodes the gap junction protein connexin 26, was markedly reduced in SCO testes (Fig 3A).

Reduced transcriptional expression of proteins that polarize and organize the plasma membrane

Although CLDN11 expression is increased in SCO testes, the protein is not concentrated at the presumptive sites of the blood-testis barrier [39, 40]. Rather, most is distributed across the entire Sertoli cell plasma membrane, indicating that these cells lack normal apical-basal polarity and/or membrane microdomains that are required to concentrate adhesive proteins at sites of cell-cell contact [41]. One probable cause for loss of polarity in these cells is inadequate expression of Crumbs (Crb) polarity complex proteins that mark the boundary between the apical and basal-lateral plasma membrane and adaptor proteins (TJPs) that stabilize the tight junctions assembled at that site [41]. Sertoli cell signature transcripts encode CRB2 and TJP3, and their expression in SCO testes were reduced to 16% and 1%, respectively, of their expression in testes with complete spermatogenesis (Fig 3B). We recognize that other polarity and adaptor proteins, CRB3, TJP1 and TJP2, have been implicated in formation of the blood-testis barrier by rodent Sertoli cells [42, 43]. While we did not define these 3 proteins as Sertoli cell signature transcripts, it should be noted that when the abundance of their transcripts in total testis transcriptomes was normalized to beta actin mRNA, expressions of all three were significantly reduced in SCO testes (S4 Fig). Thus, markedly reduced expression of CRB2 and TJP3 in SCO testes is not compensated by increased expression of other polarity and adaptor proteins.

The localization of adhesive and gap junctions within sites in the plasma membrane that abut other cells requires a higher level of membrane organization than is achieved by separation of the apical from the basal-lateral plasma membrane. Proteins in the spectrin family of actin-binding proteins play a key role in the formation of plasma membrane microdomains and the concentration of junctional molecules in those domains [44, 45]. Three Sertoli cell signature transcripts encode the spectrins, SPTB, SPTBN4, and SPTBN5, which were reduced in expression from 23% to 60% in SCO testes (Fig 3B). These decreases and those for transcripts encoding CRB2 and TJP3 lead us to conclude that in SCO testes Sertoli cells can form neither a barrier that separates the apical and basal-lateral plasma membrane nor the membrane microdomains that normally concentrate cell adhesion complexes at sites of cell-cell contact.

Reduced expression of genes that act in pathways, which regulate vesicular trafficking

An essential requirement for male fertility is the dynamic assembly and disassembly of junctions between adjacent Sertoli cells and between these somatic cells and spermatogenic cells [46, 47]. This dynamism requires vesicular trafficking of membrane proteins in and out of the plasma membrane [48]. This trafficking is regulated by a pathway that begins with the binding of RAS GTPases to catalytic subunits of phosphoinositide 3-kinase (PI3K) [49]. The resulting activation of PI3K enzyme activity stimulates RAB GTPases, the immediate regulators of exocytic and endocytic vesicular trafficking [15, 50]. Sertoli cells within SCO testes expressed markedly reduced levels transcripts that encode key enzymes in this pathway. Expressions of RRAS2 and HRAS mRNAs were reduced to 40% of what we measured in testes with complete spermatogenesis; expression of the catalytic subunit of PI3K, PI3CA, was reduced by 50% (Fig 4A). In contrast, expression of two PI3K regulatory subunits, PI3KR1 and PI3KR2, were increased or similar to expression in testes with complete spermatogenesis (Fig 4B). Multiple RABS encoded by Sertoli cell signature transcripts were also significantly reduced in SCO testes (Fig 4C and S3 Table). For example, SCO testes expression of RAB20, a regulator of endocytic vesicular trafficking, was reduced by 65% [51], while expression of the regulators of exocytic vesicles, RAB3D and RAB40B, were reduced by 60% and 70%, respectively [52, 53]. Thus, Sertoli cells within SCO testes exhibit significant deficits in expression of multiple genes that act at different points in the pathway that regulates the vesicular trafficking of junctional proteins to and from the plasma membrane.

[Figure omitted. See PDF.]

Fig 4. Expression of Sertoli cell signature transcripts encoding regulators of vesicular trafficking.

(A) Expression of RRAS2 and HRAS and the catalytic subunit of phosphoinositide 3-kinase (PI3K), PIK3CA. Binding of RRAS2 and HRAS to PI3KCA stimulates its catalytic activity. (B) Expression of regulatory subunits of PI3K. (C) Small GTPase family members that regulate vesicular trafficking. RAB20 regulates endocytic trafficking; RAB40B and RAB3D regulate exocytic trafficking. Activities of these small GTPases are stimulated by PI3K. Data are expressed as mean + SEM of CPMs for each transcript. The presence of an asterisk over a pair of bars indicates a significant difference between testes with complete spermatogenesis and testes with SCO syndrome (FDR<0.05).

https://doi.org/10.1371/journal.pone.0216586.g004

In summary, we conclude that there are at least three reasons that Sertoli cells in SCO testes cells do not form normal cell-cell junctions with each other and with spermatogenic cells, including SSCs: First, these somatic cells are unable to form a boundary that is required for the separation of the apical from the basal-lateral plasma membrane. Thus, the plasma membrane is not properly polarized. Second, these cells do not form the proper membrane microdomains into which adhesive proteins and communicating junctions are normally concentrated. Third, there is inefficient vesicular trafficking of proteins that form adhesive and communicating junctions, thereby suppressing the formation and turnover of junctions that tether Sertoli cells to each other or to any neighboring spermatogenic cells.

Deficient expression of three growth factors that regulate replication and differentiation of SSCs

There is considerable evidence from studies of mice that multiple growth factors regulate the size and function of the spermatogonial stem cell pool [5]. Human SSCs undoubtedly share this requirement for they express receptors for all of the growth factors that are required by mouse SSCs [30]. Therefore, it follows that inadequate stimulation by growth factors may be a significant contributor to SCO syndrome.

Our analyses identified 19 Sertoli cell signature transcripts that encode growth factors or cytokines that are ligands for those receptors (S2 Table). Five factors, GDNF, FGF8, BMP4, Wnt4 and Wnt5A, have defined effects on mouse SSCs and/or progenitor spermatogonia [3, 6, 54–57]. Expressions of 3 of these growth factors were significantly reduced in SCO testes (Fig 5). Consistent with our earlier report [19], expression of GDNF in SCO testes was only 20% of expression in testes with complete spermatogenesis. FGF8 mRNA expression was also markedly reduced in SCO testes. Our analyses indicated that there was no compensation for this marked decrease in FGF8 mRNA by increased transcriptional expression of another member of the fibroblast growth factor family, FGF2, which binds the same receptor as FGF8 [29]. While FGF2 is not encoded by a Sertoli cell signature transcript, it is worth noting that when the abundance of FGF2 mRNA in the total testis transcriptome was normalized to beta actin mRNA, expression of FGF2 mRNA was significantly reduced in SCO testes (S5 Fig).

[Figure omitted. See PDF.]

Fig 5. Expression of Sertoli cell signature transcripts that encode growth factors that are known regulators of SSCs and/or progenitor spermatogonia.

Data are expressed as mean + SEM of CPMs. An asterisk over a pair of bars indicates a significant difference (FDR≤0.05).

https://doi.org/10.1371/journal.pone.0216586.g005

BMP4, WNT4 and WNT5A have distinct effects on mouse SSCs; BMP4 and WNT4 promote SSC differentiation, while WNT5A promotes stem cell maintenance [55, 57–59].

Since human SSCs express receptors for these three growth factors [30], it is probable that they have similar effects on human SSCs. Our analyses demonstrated that in SCO testes BMP4 mRNA expression was reduced to 30% of expression in testes with complete spermatogenesis. In contrast, the two WNT family members were expressed at similar levels (Fig 5). So, too, was an important regulator of the final steps in spermatogonial differentiation, kit ligand (KITLG) [60] (Fig 5).

CSF1 and CXCL12 are two additional growth factors that are known to regulate mouse SSCs [61, 62]. While they were not identified as Sertoli cell signature transcripts, it should be noted that examination of the total testis transcriptomes revealed that beta actin-normalized expression of CSF1 and CXCL12 mRNA was similar in normal and SCO testes (S5 Fig).

While our screening process identified FGF8 as a Sertoli cell signature transcript, studies of mice concluded that spermatogenic cells and not Sertoli cells are this growth factor’s source [56]. Given the almost complete lack of expression of FGF8 in SCO testes and its well-characterized effect on mouse SSCs, we performed additional experiments to prove that FGF8 was a product of human Sertoli cells, and we examined whether FGF8 protein content of SCO Sertoli cells was significantly diminished. FACs analysis demonstrated that all cells that express FGF8 also express SOX9, a specific marker of Sertoli cells [63] (Fig 6A). Furthermore, while in human testes with complete spermatogenesis FGF8 expression per Sertoli cell exhibited a normal distribution, FGF8 protein expression in SCO testes was bimodal with approximately 50% of Sertoli cells expressing abnormally low levels (Fig 6B). However, those low levels of FGF8 were substantially greater than background in this assay (Fig 6C). Finally, quantitative PCR analysis confirmed that FGF8 mRNA expression was markedly reduced in SCO testes (Fig 6D). Thus, the SSCs that are present in some human SCO testes are deprived of two growth factors (GDNF and FGF8) that promote”stemness” and one factor (BMP4) that stimulates differentiation. This suggests that the regulation of the balance between self-renewing stem cell replication and stem cell differentiation is askew.

[Figure omitted. See PDF.]

Fig 6. Confirmation that FGF8 is a product of human Sertoli cells and verification that expressions of both FGF8 protein and mRNA are reduced in SCO testes.

Testes with complete spermatogenesis (A) and with SCO syndrome (B) were dispersed into single cells and then incubated with antibodies for the Sertoli cell marker SOX9 and for FGF8. Negative control cells (C) were not incubated with antibodies. Cells were then analyzed FACs. These data are representative of three independent analyses. (D) Quantification by real-time PCR of levels of FGF8 mRNA in biopsies of human testes with complete spermatogenesis and in testes with SCO syndrome. Data (mean + SEM; n = 3) were normalized for expression of beta actin in the same sample (p<0.001).

https://doi.org/10.1371/journal.pone.0216586.g006

Discussion

There is considerable evidence that a fertile testis is characterized by precise and complex interactions between testicular somatic cells and spermatogenic cells [24, 64]. One result of these interactions is the formation of a functional niche for spermatogonial stem cells by Sertoli and other somatic cells. It follows that aberrant gene expression by Sertoli cells might decrease niche function and contribute to the dearth and dysfunction of SSCs in infertile men with SCO syndrome. Furthermore, niche dysfunction might also explain why the few SSCs in tubule segments with active spermatogenesis do not replicate and fill empty niches within SCO testes. This first thorough comparison of the transcriptome of testes with complete spermatogenesis with the transcriptome of testes with SCO syndrome was motivated by the possibility that defining the transcriptional basis of niche dysfunction would be an important first step to development of new therapies for men with SCO syndrome whose testes contain small segments of tubules with active spermatogenesis. We have focused our analyses on 244 transcripts, which our two-step screen identified as human Sertoli cell signature transcripts. Two hundred eleven of these are also expressed by mature mouse Sertoli cells. Such signature transcripts are a foundation of these cells’ critical but varied functions, and one such function is the creation of the niche for SSCs. Our analyses demonstrated that 69% of these transcripts were expressed at normal or elevated levels in SCO testes. However, 31% were expressed at significantly lower levels and their decreased expression provides probable explanations for two important characteristics of SCO testes. As summarized in Fig 7, we conclude that the failure of Sertoli cells in SCO testes to form normal cell-cell junctions results from their inability to assemble and stabilize those junctions at normal sites of cell-cell contact. We also conclude that Sertoli cells within SCO testes express abnormally low levels of GDNF, FGF8 and BMP4, thereby depriving SSCs of growth factors that have been well characterized as regulators of mouse SSCs.

[Figure omitted. See PDF.]

Fig 7. Summary: Deficits in Sertoli cell gene expression that are associated with specific aspects of human SCO syndrome.

(A) Human testis with complete spermatogenesis. Polarity complex proteins (red) and tight junctions (gray) polarize the plasma membrane thereby generating apical and basal-lateral domains (identified by up and down arrows). Additional tight junctions, ectoplasmic specialization (yellow) and gap junctions (blue and green) are concentrated at sites of cell-cell contact. Vesicular trafficking (VT; tan circles) of adhesive and gap junction proteins facilitates the assembly and disassembly of junctions, controlled by a pathway that includes Ras GTPases, PI3K and RAB small GTPases (purple shading). Sertoli cells express normal amounts of GDNF, FGF8 and BMP4, growth factors that regulate replication and differentiation of SSCs. (B) The infertile SCO testis. Inadequate expression of polarity and adaptor proteins renders Sertoli cells incapable of properly polarizing their plasma membranes. There is reduced, inefficient vesicular trafficking of proteins that form tight junctions, ectoplasmic specializations and gap junctions, stranding these proteins in vesicles within the cytoplasm. Finally, there is reduced expression of GDNF, FGF8 and BMP4. (C) The effects of GDNF, FGF8 and BMP4 on SSCs and progenitor spermatogonia. GDNF and FGF8 promote self-renewing replication of SSCs, while BMP4 stimulates differentiation of SSCs into progenitor spermatogonia and progenitors into fully differentiated B spermatogonia.

https://doi.org/10.1371/journal.pone.0216586.g007

Failure of Sertoli cells to form normal cell-cell junctions

The failure of Sertoli cells in SCO testes to form normal cell-cell junctions is not due to inadequate expression of transcripts that encode the adhesive proteins within those junctions. Our studies corroborated earlier reports that these somatic cells express elevated levels of the tight junction protein claudin 11 in SCO testes [65]. A new finding from our analyses, however, is that these Sertoli cells express normal or elevated levels of transcripts encoding many other major integral membrane adhesive proteins. Moreover, our analyses identified deficits in the expression of genes that polarize the plasma membrane or concentrate adhesive junctions in the microdomains of plasma membranes that abut other cells. Furthermore, our data identified markedly reduced expression of genes in a pathway that regulates vesicular trafficking of junctional proteins in and out of the plasma membrane.

In SCO testes much of the tight junction adhesive molecule, claudin 11, is dispersed across the entire Sertoli cell plasma membrane rather than at the presumptive sites of the blood-testis barrier [40]. This dispersal indicates that the cells have failed to segregate their plasma membranes into apical and basal-lateral domains. Our data identified deficits in expression of two genes that are the likely causes of this failure, the Crumbs (Crb) polarity complex component CRB2 and the adaptor protein, TJP3. Both normally play important roles in membrane polarization. The insertion of CRB2 into the plasma membrane marks the site at which a barrier between the apical and basal-lateral plasma membrane is constructed [41]. This insertion also recruits tight junction adhesive proteins, such as CLDN11 and JAM3, which are then anchored to the actin cytoskeleton by adaptor proteins, such as TJP3. In the absence of Crb proteins, tight junction formation remains incomplete and lacks barrier function [41]. Without adaptor proteins tight junction adhesive proteins are not anchored at the presumptive site of the barrier [66]. The expression of CRB2 and TJP3 was significantly reduced in SCO testes, and the lack of compensatory increases in expression of other CRBs or adaptor proteins makes it highly probable that Sertoli cells nether efficiently mark the site in the plasma membrane for barrier formation nor anchor tight junctions at that site. Unmoored, tight junction adhesive proteins would, thus, be dispersed across the entire Sertoli cell surface. Studies of human SCO testes have reported such dispersion [40]. This lack of membrane polarity would also place fewer constraints on the distribution of all other junctional molecules within the plasma membrane.

Our analyses indicated that the effects of diminished plasma membrane polarity on the formation of cell-cell junctions are likely further exacerbated by a significantly reduced expression of three members of the spectrin-family of actin-binding proteins, SPTB, SBTBN4 and SPTBN5. Spectrins play a key role assembling lipid rafts within a plasma membrane and in concentrating multiprotein adhesion and signaling complexes within those rafts [44, 45, 67, 68]. Thus, reduced expression of these three spectrins synergizes with inadequate expression of CRB2 and TJP3 to disrupt the ability of Sertoli cells to form adhesive junctions with each other and with spermatogenic cells, including SSCs.

Reduced expression of genes in pathways that regulate vesicular trafficking

In a normal testis, interactions between adjacent Sertoli cells and between Sertoli cells and germ cells, including spermatogonia, are dynamic. Movement of preleptotene spermatocytes into the adluminal compartment of the seminiferous tubule requires simultaneous elimination of Sertoli-Sertoli junctions above the germ cells and assembly of new junctions below [47, 69]. Concurrently, junctions that hold elongate spermatids to Sertoli cells are disassembled in preparation for spermiation [70]. Vesicular trafficking is responsible for assembling and disassembling cell-cell junctions, and phosphoinositide 3-kinase (PI3K) and the pathway in which it acts are key traffic regulators [50]. As binding of HRAS and RRAS-2 to PI3K’s catalytic subunits stimulates PI3K activity [49], the markedly reduced expression of HRAS and RRAS-2 and of PI3CA must greatly diminish Sertoli cell PI3K enzyme activity in Sertoli cells within SCO testes. We suggest that coupling this reduction with significantly lower expression of some RAB GTPases decreases the efficiency of vesicular trafficking to and from the plasma membrane, leaving some vesicles and the junctional proteins they contain stalled in the cytoplasm. This suggestion is supported by the report that in SCO testes some claudin 11 is located in the cytoplasm of Sertoli cells [71]. It follows that inefficient trafficking further reduces the ability of these Sertoli cells to assemble adhesive and communicating junctions at sites of cell-cell contact.

Connexin 26

An unexpected finding from our analyses was that Sertoli cells in SCO testes expressed markedly reduced levels of connexin 26 (GJB2). In rodents connexin 26 is concentrated at sites where Sertoli cells and spermatogonia abut [37]. It is, therefore, of interest that rodent Sertoli cells communicate efficiently via gap junctions with spermatogonia, and spermatocytes [72, 73]. While GJB2 is also expressed by human spermatogonia (S1 Table), the reduction in GJB2 expression to 2% of what is observed in testes with complete spermatogenesis suggests that there is reduced gap junctional communication between Sertoli cells and the few SSCs that exist in an SCO testis.

The importance of cell adhesion to SSCs

Since Sertoli cells form adhesive junctions with spermatogonia [40], one must consider the possible consequence of loss of cell adhesion to SSCs. First, this loss will reduce signaling by focal adhesion kinases within stem cells [74]. Second, decreased adhesive strength deprives SSCs of the positional information they need to divide asymmetrically. (Studies of pig SSCs document such divisions [16]). Asymmetric division is a common strategy that adult stem cells use to balance generation of new stem cells with the production of differentiating progeny [16, 75, 76]. Thus, the loss of positional information in SCO testes reduces the probability of asymmetric division by SSCs, thereby threatening the maintenance of the stem cell pool.

Deficient expression of growth factors that regulate SSCs

Replication and differentiation of SSCs are regulated by multiple growth factors, and our analyses demonstrated that GDNF, FGF8 and BMP4 were expressed at abnormally low levels in SCO testes. All three are well-characterized regulators of mouse SSCs and progenitor spermatogonia [6, 56, 58, 59]. Sertoli cells are the only major source of GDNF in the human testis, and between the current and previous studies [19], we have evaluated GDNF expression in 21 separate human SCO testes. Here and in our previous study, the average GDNF mRNA level in SCO testes was reduced to approximately 20% of normal, and all SCO testes expressed lower levels of GDNF than testis with complete spermatogenesis. Studies in mice suggest that as little as a 50% reduction in GDNF expression is sufficient to cause gradual development of SCO syndrome [77]. Thus, it is becoming increasingly clear that inadequate GDNF expression by Sertoli cells is a fundamental characteristic, as well as one probable cause, of human SCO syndrome.

Our studies led to important new conclusions: human Sertoli cells are the predominant testicular source of FGF8, and expression of this growth factor is less than 10% of normal in SCO testes. Additionally, our data indicate that there is no compensatory increase by SCO testes in FGF2 expression, which binds the same receptors as FGF8 [29]. The probable consequence of this 95% decrease in FGF8 expression is evident in studies of mice. FGF8 suppresses the early steps in differentiation of SSCs and early progenitor spermatogonia, thereby promoting the maintenance of a normal stem cell pool [56]. Thus, studies of mouse SSCs predict that the greatly diminished expression of both GDNF and FGF8 in SCO testes deprives these stem cells of stimuli that both promotes self-renewing replication and suppresses differentiation of these stem cells. Furthermore, loss of adequate GDNF stimulation deprives the few remaining SSCs of chemotactic stimuli that would normally drive these stem cells to migrate to empty niches [18, 78]. BMP4 is the third growth factor that we identified as expressed at substantially lower levels in SCO testes. Studies of mice indicate that this growth factor stimulates stem cells to form progenitor spermatogonia and progenitors to express Kit, the receptor for Kit ligand, a growth factor that drives the final steps in differentiation of human progenitor spermatogonia into Type B spermatogonia [3, 58]. Thus, while Kit ligand expression was at normal levels in SCO testes, its effect on progenitors is probably dampened by the lack of BMP4-driven Kit expression by these cells. Thus, the reduced expression of BMP4 provides one explanation for our observation that there is a substantially greater decrease in expression of markers for spermatocytes and spermatids than for SSCs in SCO testes.

Next steps to a full understanding of Sertoli cell-only syndrome

This study described herein is but a first step towards a full understanding of the molecular basis of Sertoli cell-only syndrome. It is critical that we attain this understanding. First, Sertoli cell-only syndrome is one of the severest forms of male infertility, which is a major medical and public health problem [1]. Second, small segments of tubules in SCO testes produce sperm, and thus must contain SSCs [2]. Third, there is a pressing need to develop new therapies for human male infertility. Fourth, identifying aberrations in gene expression by SCO Sertoli cells might lead to new approaches to correct or compensate for niche dysfunction, thereby increasing the chances of SCO infertile men becoming biological fathers. The data presented described herein are a necessary first step to defining those aberrations. However, as with all first steps, this one has its limitations. Because the transcriptome of normal adult human Sertoli cells has not been characterized directly, we relied on a rigorous two-step bioinformatic analytical approach to define human Sertoli cell signature transcripts. As we previously noted, all mammalian Sertoli cells perform the same essential functions and it is highly probable that Sertoli cells of different species express orthologous genes in order to perform these functions. Nonetheless, a next important step will be to thoroughly describe the entire transcriptome of highly purified human Sertoli cells that have been rapidly isolated from fertile testes. A prerequisite for this step is the development of the isolation protocol. Rapidity is necessary to minimize changes in these cells’ transcriptomes. A second limitation of the current study is that the small mass and immediate freezing of the biopsies (15–40 mg) precluded the isolation of Sertoli cells prior to the evaluation of their transcriptomes. We anticipate that this second limitation will be circumvented by single cell transcriptome analysis of cells in freshly obtained biopsies. However, clear and compelling scientific rationale is required before such a study can be proposed and conducted. Our analyses described herein build such rationale. Finally, it is likely that a comparison of the transcriptomes of dialated, germ cell-filled tubules and atretic tubules from the same SCO testes would provide significant insights into the transcriptional foundation of SCO syndrome. However, as dialated tubules are potential sources of sperm, and as the clinical objective of the micro-Tese procedure is to help infertile men conceive a child, it is necessary that all dialated tubules be used for sperm retrieval.

We also acknowledge that signals from spermatogenic cells regulate a subset of the genes expressed by Sertoli cells. Thus, some of the aberrations in Sertoli cell gene expression noted herein could be secondary to loss of such regulatory signals [7]. However, we doubt that lack of germ cells is the primary cause of the aberrant gene expression by Sertoli cells in SCO testes. Consider the following: While GDNF mRNA expression is markedly reduced in SCO testes, its expression is markedly increased in germ cell-deficient mouse testes [79]. Additionally, human Sertoli cells in SCO testes do not form a functional blood-testis barrier [13, 14], and our data indicate that they do not express adequate levels of genes required for polarization of their plasma membrane or efficient transportation of cell adhesion proteins to the plasma membrane. In contrast, Sertoli cells in germ cell deficient rat testes form a functional blood-testis barrier [80]. It follows that those rat Sertoli cells must express adequate levels of the genes that are essential for polarization of the Sertoli cell plasma membrane and assembly of junctions at sites of normal cell-contact. Finally, others and we have proven that both stimulatory and inhibitory signals from spermatogenic cells regulate gene expression by Sertoli cells. This complex regulation is clearly evident in the stage-specific changes in rat Sertoli cell expression of the cathepsin L (CTSL) gene [7]. Stimulatory and inhibitor signals from germ cells act via enhancers and repressors in the CTSL gene promoter to regulate stage-specific transcription. Loss of germ cells does not substantially change total testis content of cathepsin L mRNA; rather all Sertoli cells express similar levels of the transcript. Other genes show a similar response to loss of germ cells [7]. Thus, we propose that the aberrant expression of some and possibly many genes by Sertoli cells in SCO testis is a characteristic of the syndrome and is not secondary to the loss of spermatogenic cells.

Conclusion

The goal of our analyses was to identify changes in gene expression by human Sertoli cells that might explain two characteristics of the dysfunctional stem cell niche in SCO testes and thereby potentially identify some of the mechanistic underpinnings of this infertile human phenotype (Fig 7). The inability of these Sertoli cells to form normal cell-cell junctions can be explained by their inadequate expression of genes required to properly polarize and organize the plasma membrane and by inefficient cycling of proteins that form adhesive and communicating junctions in and out of the plasma membrane. Sertoli cells in SCO testes also express abnormally low levels of three growth factors GDNF, FGF8 and BMP4, which are important regulators of the replication and differentiation of mouse SSCs and progenitor spermatogonia. We propose that throughout most of an SCO testis, GDNF and FGF8 levels are inadequate to stimulate both the replication of SSCs and their migration into germ cell-deficient areas of seminiferous tubules. Low levels of BMP4 may deprive stimulus that normally drives differentiation of SSCs and progenitor spermatogonia. Thus, our data indicate that SCO syndrome is a complex condition associated with multiple deficits in Sertoli cell gene expression. We acknowledge that further studies are needed to gain a more clear and precise understanding of the contribution of Sertoli cells to SCO syndrome. We hope that such an understanding might foster the development of new therapies to stimulate the expansion of the few spermatogonial stem cells that are present some SCO testes, thereby increasing the possibility that an infertile man can become a father.

Supporting information

[Figure omitted. See PDF.]

S1 Text. Supplementary methods.

https://doi.org/10.1371/journal.pone.0216586.s001

(DOCX)

S1 Fig. The histologies of human testes with normal spermatogenesis, with mild hypospermatogenesis and moderate hypospermatogenesis.

The testes with these three histologies are all referred to in the manuscript as exhibiting complete spermatogenesis. Photomicrographs were captured using 10X and 40x microscope lenses.

https://doi.org/10.1371/journal.pone.0216586.s002

(TIF)

S2 Fig. Abundance of transcripts characteristic of SSCs (GFRA1), spermatocytes (SYCP2, SYCP3) and spermatids (PRM1, PRM2) in each of the 4 testes with complete spermatogenesis (C1-C4) and in each of the 7 testes with SCO Syndrome (SCO5, SCO5-12).

Data are from total testis transcriptomes and are expressed as CPM for each of the above transcripts divided by CPM for beta actin.

https://doi.org/10.1371/journal.pone.0216586.s003

(TIF)

S3 Fig. Increased abundance in SCO testes of four Sertoli cell signature transcripts essential for spermatogenesis.

Data (mean+ SEM) are expressed CPM. An asterisk over a pair of bars indicates a significant difference between testes with complete spermatogenesis and SCO testes (FDR≤ 0.05).

https://doi.org/10.1371/journal.pone.0216586.s004

(TIF)

S4 Fig. Abundance of transcripts encoding the cell polarity protein, CRB2 and the adapter proteins TJP1 and TJP2 in the total testis transcriptomes of normal and SCO testes.

Data (mean + CPM) are expressed as CPM in the total testis transcriptome divided by CPM of ACTB in the same sample. Asterisks over a pair of bars indicate that normalized expression of a transcript differs between testes with complete spermatogenesis and SCO testes (p≤0.005).

https://doi.org/10.1371/journal.pone.0216586.s005

(TIF)

S5 Fig. Abundance in the total testis transcriptomes of normal and SCO testes of FGF2, CSF1 and CXCL12, three growth factors that have been demonstrated in studies of mice to regulate SSCs or progenitor spermatogonia.

Data (mean + SEM) are expressed as CPM in the total testis transcriptome divided by CPM of ACTB in the same sample. An asterisk over a pair of bars indicates a significant difference between testes with complete spermatogenesis and SCO testes (p≤0.005).

https://doi.org/10.1371/journal.pone.0216586.s006

(TIF)

S1 Table. Identification of transcripts expressed at least 4-fold higher by rat Sertoli cells than by rat Leydig cells, pachytene spermatocytes, round spermatids and spermatogonia.

https://doi.org/10.1371/journal.pone.0216586.s007

(XLSX)

S2 Table. Definition of human Sertoli cell signature transcripts.

https://doi.org/10.1371/journal.pone.0216586.s008

(XLSX)

S3 Table. Expression of Sertoli cell signature transcripts in testis with complete spermatogenesis and in testes with Sertoli cell-only syndrome.

https://doi.org/10.1371/journal.pone.0216586.s009

(XLS)

Acknowledgments

We thank Alex Bolyakov and Anna Mielnik for excellent technical support. We thank Qianyi Ma and Sue Hammoud for providing us with the list of transcripts expressed by mature mouse Sertoli cells. We thank Pamela Wright for expert assistance with editing of the manuscript and with preparation of figures.

Citation: Paduch DA, Hilz S, Grimson A, Schlegel PN, Jedlicka AE, Wright WW (2019) Aberrant gene expression by Sertoli cells in infertile men with Sertoli cell-only syndrome. PLoS ONE 14(5): e0216586. https://doi.org/10.1371/journal.pone.0216586

References

1. Irvine DS. Epidemiology and aetiology of male infertility. Hum Reprod. 1998;13 Suppl 1:33–44. Epub 1998/07/15. pmid:9663768.

2. Schlegel PN, Palermo GD, Goldstein M, Menendez S, Zaninovic N, Veeck LL, et al. Testicular sperm extraction with intracytoplasmic sperm injection for nonobstructive azoospermia. Urology. 1997;49(3):435–40. Epub 1997/03/01. pmid:9123710.

3. Nagano M, Ryu BY, Brinster CJ, Avarbock MR, Brinster RL. Maintenance of mouse male germ line stem cells in vitro. Biol Reprod. 2003;68(6):2207–14. Epub 2003/02/28. pmid:12606373.

4. Ogawa T, Ohmura M, Yumura Y, Sawada H, Kubota Y. Expansion of murine spermatogonial stem cells through serial transplantation. Biol Reprod. 2003;68(1):316–22. Epub 2002/12/21. pmid:12493728.

5. Oatley JM, Brinster RL. The germline stem cell niche unit in mammalian testes. Physiol Rev. 2012;92(2):577–95. Epub 2012/04/27. pmid:22535892; PubMed Central PMCID: PMC3970841.

6. Parker N, Falk H, Singh D, Fidaleo A, Smith B, Lopez MS, et al. Responses to glial cell line-derived neurotrophic factor change in mice as spermatogonial stem cells form progenitor spermatogonia which replicate and give rise to more differentiated progeny. Biol Reprod. 2014;91(4):92. Epub 2014/08/29. pmid:25165119; PubMed Central PMCID: PMC4435032.

7. Wright WW. Stage-specific gene expression by Sertoli Cells. In: Griswold M, editor. Sertoli Cell Biology. Amsterdam, Netherlands: Elsevier; 2015. p. 273–306.

8. Stanley EL, Johnston DS, Fan J, Papadopoulos V, Chen H, Ge RS, et al. Stem Leydig cell differentiation: gene expression during development of the adult rat population of Leydig cells. Biol Reprod. 2011;85(6):1161–6. Epub 2011/08/13. pmid:21832170; PubMed Central PMCID: PMC3223250.

9. Jan SZ, Vormer TL, Jongejan A, Roling MD, Silber SJ, de Rooij DG, et al. Unraveling transcriptome dynamics in human spermatogenesis. Development. 2017;144(20):3659–73. Epub 2017/09/25. pmid:28935708; PubMed Central PMCID: PMCPMC5675447.

10. Lottrup G, Belling K, Leffers H, Nielsen JE, Dalgaard MD, Juul A, et al. Comparison of global gene expression profiles of microdissected human foetal Leydig cells with their normal and hyperplastic adult equivalents. Mol Hum Reprod. 2017;23(5):339–54. Epub 2017/03/24. pmid:28333300.

11. Paniagua R, Codesal J, Nistal M, Rodriguez MC, Santamaria L. Quantification of cell types throughout the cycle of the human seminiferous epithelium and their DNA content. A new approach to the spermatogonial stem cell in man. Anat Embryol (Berl). 1987;176(2):225–30. Epub 1987/01/01. pmid:3619074.

12. Petersen PM, Seieroe K, Pakkenberg B. The total number of Leydig and Sertoli cells in the testes of men across various age groups—a stereological study. J Anat. 2015;226(2):175–9. Epub 2014/12/30. pmid:25545958; PubMed Central PMCID: PMCPMC4304573.

13. Camatini M, Franchi E, Decurtis I. Permeability to lanthanum of blood testis barrier in human germinal aplasia. Anat Rec. 1981;200(3):293–7. Epub 1981/07/01. pmid:7270929.

14. Matsuo Y, Nomata K, Eguchi J, Aoki D, Hayashi T, Hishikawa Y, et al. Immunohistochemical analysis of connexin43 expression in infertile human testes. Acta Histochem Cytochem. 2007;40(3):69–75. Epub 2007/07/27. pmid:17653298; PubMed Central PMCID: PMCPMC1931485.

15. Mruk DD, Cheng CY. The Mammalian Blood-Testis Barrier: Its Biology and Regulation. Endocr Rev. 2015;36(5):564–91. Epub 2015/09/12. pmid:26357922; PubMed Central PMCID: PMC4591527.

16. Luo J, Megee S, Dobrinski I. Asymmetric distribution of UCH-L1 in spermatogonia is associated with maintenance and differentiation of spermatogonial stem cells. J Cell Physiol. 2009;220(2):460–8. Epub 2009/04/24. pmid:19388011; PubMed Central PMCID: PMCPMC2732714.

17. Spradling A, Fuller MT, Braun RE, Yoshida S. Germline stem cells. Cold Spring Harb Perspect Biol. 2011;3(11):a002642. Epub 2011/07/28. pmid:21791699; PubMed Central PMCID: PMCPMC3220357.

18. Dovere L, Fera S, Grasso M, Lamberti D, Gargioli C, Muciaccia B, et al. The niche-derived glial cell line-derived neurotrophic factor (GDNF) induces migration of mouse spermatogonial stem/progenitor cells. PLoS One. 2013;8(4):e59431. Epub 2013/04/25. pmid:23613711; PubMed Central PMCID: PMCPMC3632550.

19. Singh D, Paduch DA, Schlegel PN, Orwig KE, Mielnik A, Bolyakov A, et al. The production of glial cell line-derived neurotrophic factor by human sertoli cells is substantially reduced in sertoli cell-only testes. Hum Reprod. 2017;32(5):1108–17. Epub 2017/04/04. pmid:28369535.

20. Schlegel PN. Causes of azoospermia and their management. Reprod Fertil Dev. 2004;16(5):561–72. Epub 2004/09/16. pmid:15367371.

21. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–60. Epub 2009/05/20. pmid:19451168; PubMed Central PMCID: PMCPMC2705234.

22. Anders S, Pyl PT, Huber W. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9. Epub 2014/09/28. pmid:25260700; PubMed Central PMCID: PMCPMC4287950.

23. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40. Epub 2009/11/17. pmid:19910308; PubMed Central PMCID: PMCPMC2796818.

24. Johnston DS, Wright WW, Dicandeloro P, Wilson E, Kopf GS, Jelinsky SA. Stage-specific gene expression is a fundamental characteristic of rat spermatogenic cells and Sertoli cells. Proc Natl Acad Sci U S A. 2008;105(24):8315–20. Epub 2008/06/12. pmid:18544648; PubMed Central PMCID: PMC2448834.

25. Green CD, Ma Q, Manske GL, Shami AN, Zheng X, Marini S, et al. A Comprehensive Roadmap of Murine Spermatogenesis Defined by Single-Cell RNA-Seq. Dev Cell. 2018;46(5):651–67 e10. Epub 2018/08/28. pmid:30146481.

26. Hasegawa K, Saga Y. Retinoic acid signaling in Sertoli cells regulates organization of the blood-testis barrier through cyclical changes in gene expression. Development. 2012;139(23):4347–55. Epub 2012/10/26. pmid:23095883.

27. Cowan G, Childs AJ, Anderson RA, Saunders PT. Establishment of long-term monolayer cultures of somatic cells from human fetal testes and expansion of peritubular myoid cells in the presence of androgen. Reproduction. 2010;139(4):749–57. Epub 2010/01/22. pmid:20089665.

28. Marchetti C, Hamdane M, Mitchell V, Mayo K, Devisme L, Rigot JM, et al. Immunolocalization of inhibin and activin alpha and betaB subunits and expression of corresponding messenger RNAs in the human adult testis. Biol Reprod. 2003;68(1):230–5. Epub 2002/12/21. pmid:12493718.

29. Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, et al. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45(D1):D362–D8. Epub 2016/12/08. pmid:27924014; PubMed Central PMCID: PMCPMC5210637.

30. Guo J, Grow EJ, Yi C, Mlcochova H, Maher GJ, Lindskog C, et al. Chromatin and Single-Cell RNA-Seq Profiling Reveal Dynamic Signaling and Metabolic Transitions during Human Spermatogonial Stem Cell Development. Cell Stem Cell. 2017;21(4):533–46 e6. Epub 2017/10/07. pmid:28985528; PubMed Central PMCID: PMCPMC5832720.

31. Wang YM, Sullivan PM, Petrusz P, Yarbrough W, Joseph DR. The androgen-binding protein gene is expressed in CD1 mouse testis. Mol Cell Endocrinol. 1989;63(1–2):85–92. Epub 1989/05/01. pmid:2753230.

32. Selva DM, Hammond GL. Human sex hormone-binding globulin is expressed in testicular germ cells and not in sertoli cells. Horm Metab Res. 2006;38(4):230–5. Epub 2006/05/16. pmid:16700003.

33. Wang XN, Li ZS, Ren Y, Jiang T, Wang YQ, Chen M, et al. The Wilms tumor gene, Wt1, is critical for mouse spermatogenesis via regulation of sertoli cell polarity and is associated with non-obstructive azoospermia in humans. PLoS Genet. 2013;9(8):e1003645. Epub 2013/08/13. pmid:23935527; PubMed Central PMCID: PMCPMC3731222.

34. 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. Epub 2000/11/25. pmid:11090455.

35. Das DS, Wadhwa N, Kunj N, Sarda K, Pradhan BS, Majumdar SS. Dickkopf homolog 3 (DKK3) plays a crucial role upstream of WNT/beta-CATENIN signaling for Sertoli cell mediated regulation of spermatogenesis. PLoS One. 2013;8(5):e63603. Epub 2013/05/15. pmid:23667645; PubMed Central PMCID: PMCPMC3647036.

36. Behrens A, Genoud N, Naumann H, Rulicke T, Janett F, Heppner FL, et al. Absence of the prion protein homologue Doppel causes male sterility. Embo J. 2002;21(14):3652–8. Epub 2002/07/12. pmid:12110578; PubMed Central PMCID: PMCPMC125402.

37. Kidder GM, Cyr DG. Roles of connexins in testis development and spermatogenesis. Semin Cell Dev Biol. 2016;50:22–30. Epub 2016/01/19. pmid:26780117.

38. Cheng CY, Mruk DD. Biochemistry of Sertoli cell/gern cell junctions, germ cell transport, and spermiation in the seminiferous epithelium. In: Griswold MD, editor. Sertoli Cell Biology, Second Edition. Amsterdam, Netherlands: Elsevier; 2015. p. 333–83.

39. Fink C, Weigel R, Fink L, Wilhelm J, Kliesch S, Zeiler M, et al. Claudin-11 is over-expressed and dislocated from the blood-testis barrier in Sertoli cells associated with testicular intraepithelial neoplasia in men. Histochem Cell Biol. 2009;131(6):755–64. Epub 2009/02/26. pmid:19241088.

40. Stammler A, Luftner BU, Kliesch S, Weidner W, Bergmann M, Middendorff R, et al. Highly Conserved Testicular Localization of Claudin-11 in Normal and Impaired Spermatogenesis. PLoS One. 2016;11(8):e0160349. Epub 2016/08/04. pmid:27486954; PubMed Central PMCID: PMC4972306.

41. Rodriguez-Boulan E, Macara IG. Organization and execution of the epithelial polarity programme. Nat Rev Mol Cell Biol. 2014;15(4):225–42. Epub 2014/03/22. pmid:24651541; PubMed Central PMCID: PMCPMC4211427.

42. Gao Y, Lui WY, Lee WM, Cheng CY. Polarity protein Crumbs homolog-3 (CRB3) regulates ectoplasmic specialization dynamics through its action on F-actin organization in Sertoli cells. Sci Rep. 2016;6:28589. Epub 2016/07/01. pmid:27358069; PubMed Central PMCID: PMCPMC4928075.

43. Gruber M, Mathew LK, Runge AC, Garcia JA, Simon MC. EPAS1 Is Required for Spermatogenesis in the Postnatal Mouse Testis. Biol Reprod. 2010;82(6):1227–36. Epub 2010/02/26. pmid:20181618; PubMed Central PMCID: PMC2874504.

44. Bennett V, Baines AJ. Spectrin and ankyrin-based pathways: metazoan inventions for integrating cells into tissues. Physiol Rev. 2001;81(3):1353–92. Epub 2001/06/28. pmid:11427698.

45. Machnicka B, Czogalla A, Hryniewicz-Jankowska A, Boguslawska DM, Grochowalska R, Heger E, et al. Spectrins: a structural platform for stabilization and activation of membrane channels, receptors and transporters. Biochim Biophys Acta. 2014;1838(2):620–34. Epub 2013/05/16. pmid:23673272.

46. Siu MK, Wong CH, Lee WM, Cheng CY. Sertoli-germ cell anchoring junction dynamics in the testis are regulated by an interplay of lipid and protein kinases. J Biol Chem. 2005;280(26):25029–47. Epub 2005/05/05. pmid:15870075.

47. Smith BE, Braun RE. Germ cell migration across Sertoli cell tight junctions. Science. 2012;338(6108):798–802. Epub 2012/09/22. pmid:22997133; PubMed Central PMCID: PMCPMC3694388.

48. Xiao X, Mruk DD, Wong EW, Lee WM, Han D, Wong CK, et al. Differential effects of c-Src and c-Yes on the endocytic vesicle-mediated trafficking events at the Sertoli cell blood-testis barrier: an in vitro study. Am J Physiol Endocrinol Metab. 2014;307(7):E553–62. Epub 2014/08/15. pmid:25117412; PubMed Central PMCID: PMC4187029.

49. Nakhaeizadeh H, Amin E, Nakhaei-Rad S, Dvorsky R, Ahmadian MR. The RAS-Effector Interface: Isoform-Specific Differences in the Effector Binding Regions. PLoS One. 2016;11(12):e0167145. Epub 2016/12/10. pmid:27936046; PubMed Central PMCID: PMCPMC5147862.

50. Bhattacharya S, McElhanon KE, Gushchina LV, Weisleder N. Role of phosphatidylinositol-4,5-bisphosphate 3-kinase signaling in vesicular trafficking. Life Sci. 2016;167:39–45. Epub 2016/10/21. pmid:27760304; PubMed Central PMCID: PMCPMC5204450.

51. Egami Y, Araki N. Rab20 regulates phagosome maturation in RAW264 macrophages during Fc gamma receptor-mediated phagocytosis. PLoS One. 2012;7(4):e35663. Epub 2012/05/01. pmid:22545127; PubMed Central PMCID: PMCPMC3335809.

52. Jacob A, Jing J, Lee J, Schedin P, Gilbert SM, Peden AA, et al. Rab40b regulates trafficking of MMP2 and MMP9 during invadopodia formation and invasion of breast cancer cells. J Cell Sci. 2013;126(Pt 20):4647–58. Epub 2013/08/02. pmid:23902685; PubMed Central PMCID: PMCPMC3795337.

53. Millar AL, Pavios NJ, Xu J, Zheng MH. Rab3D: a regulator of exocytosis in non-neuronal cells. Histol Histopathol. 2002;17(3):929–36. Epub 2002/08/10. pmid:12168804.

54. Boyer A, Hermo L, Paquet M, Robaire B, Boerboom D. Seminiferous tubule degeneration and infertility in mice with sustained activation of WNT/CTNNB1 signaling in sertoli cells. Biol Reprod. 2008;79(3):475–85. Epub 2008/05/16. pmid:18480464.

55. Boyer A, Yeh JR, Zhang X, Paquet M, Gaudin A, Nagano MC, et al. CTNNB1 signaling in sertoli cells downregulates spermatogonial stem cell activity via WNT4. PLoS One. 2012;7(1):e29764. Epub 2012/01/19. pmid:22253774; PubMed Central PMCID: PMCPMC3257228.

56. Hasegawa K, Saga Y. FGF8-FGFR1 signaling acts as a niche factor for maintaining undifferentiated spermatogonia in the mouse. Biol Reprod. 2014;91(6):145. Epub 2014/11/02. pmid:25359900.

57. Yeh JR, Zhang X, Nagano MC. Wnt5a is a cell-extrinsic factor that supports self-renewal of mouse spermatogonial stem cells. J Cell Sci. 2011;124(Pt 14):2357–66. Epub 2011/06/23. pmid:21693582.

58. Pellegrini M, Grimaldi P, Rossi P, Geremia R, Dolci S. Developmental expression of BMP4/ALK3/SMAD5 signaling pathway in the mouse testis: a potential role of BMP4 in spermatogonia differentiation. J Cell Sci. 2003;116(Pt 16):3363–72. Epub 2003/07/15. pmid:12857787.

59. Yang Y, Feng Y, Feng X, Liao S, Wang X, Gan H, et al. BMP4 Cooperates with Retinoic Acid to Induce the Expression of Differentiation Markers in Cultured Mouse Spermatogonia. Stem Cells Int. 2016;2016:9536192. Epub 2016/11/01. pmid:27795714; PubMed Central PMCID: PMC5067322.

60. Ohta H, Yomogida K, Dohmae K, Nishimune Y. Regulation of proliferation and differentiation in spermatogonial stem cells: the role of c-kit and its ligand SCF. Development. 2000;127(10):2125–31. Epub 2000/04/19. pmid:10769236.

61. Oatley JM, Oatley MJ, Avarbock MR, Tobias JW, Brinster RL. Colony stimulating factor 1 is an extrinsic stimulator of mouse spermatogonial stem cell self-renewal. Development. 2009;136(7):1191–9. Epub 2009/03/10. pmid:19270176; PubMed Central PMCID: PMC2685936.

62. Yang QE, Kim D, Kaucher A, Oatley MJ, Oatley JM. CXCL12-CXCR4 signaling is required for the maintenance of mouse spermatogonial stem cells. J Cell Sci. 2013;126(Pt 4):1009–20. Epub 2012/12/15. pmid:23239029; PubMed Central PMCID: PMCPMC4074255.

63. Chui K, Trivedi A, Cheng CY, Cherbavaz DB, Dazin PF, Huynh AL, et al. Characterization and Functionality of Proliferative Human Sertoli Cells. Cell Transplant. 2010. Epub 2010/11/09. pmid:21054948.

64. Johnston DS, Olivas E, DiCandeloro P, Wright WW. Stage-specific changes in GDNF expression by rat Sertoli cells: a possible regulator of the replication and differentiation of stem spermatogonia. Biol Reprod. 2011;85(4):763–9. Epub 2011/06/10. pmid:21653894; PubMed Central PMCID: PMC3184291.

65. Chiba K, Kondo Y, Yamaguchi K, Miyake H, Fujisawa M. Inhibition of claudin-11 and occludin expression in rat Sertoli cells by mono-(2-ethylhexyl) phthalate through p44/42 mitogen-activated protein kinase pathway. J Androl. 2012;33(3):368–74. Epub 2011/08/27. pmid:21868750.

66. Hartsock A, Nelson WJ. Adherens and tight junctions: structure, function and connections to the actin cytoskeleton. Biochim Biophys Acta. 2008;1778(3):660–9. Epub 2007/09/15. pmid:17854762; PubMed Central PMCID: PMCPMC2682436.

67. Nusrat A, Parkos CA, Verkade P, Foley CS, Liang TW, Innis-Whitehouse W, et al. Tight junctions are membrane microdomains. J Cell Sci. 2000;113 (Pt 10):1771–81. Epub 2000/04/19. pmid:10769208.

68. Pradhan D, Lombardo CR, Roe S, Rimm DL, Morrow JS. alpha -Catenin binds directly to spectrin and facilitates spectrin-membrane assembly in vivo. J Biol Chem. 2001;276(6):4175–81. Epub 2000/11/09. pmid:11069925.

69. Cheng CY, Mruk DD. Regulation of blood-testis barrier dynamics by focal adhesion kinase (FAK): an unexpected turn of events. Cell Cycle. 2009;8(21):3493–9. Epub 2009/10/14. pmid:19823026; PubMed Central PMCID: PMCPMC2804914.

70. Cheng CY, Mruk DD. An intracellular trafficking pathway in the seminiferous epithelium regulating spermatogenesis: a biochemical and molecular perspective. Crit Rev Biochem Mol Biol. 2009;44(5):245–63. Epub 2009/07/23. pmid:19622063; PubMed Central PMCID: PMCPMC2758293.

71. Chiba K, Yamaguchi K, Ando M, Miyake H, Fujisawa M. Expression pattern of testicular claudin-11 in infertile men. Urology. 2012;80(5):1161 e13–7. Epub 2012/09/07. pmid:22951003.

72. Risley MS, Tan IP, Farrell J. Gap junctions with varied permeability properties establish cell-type specific communication pathways in the rat seminiferous epithelium. Biol Reprod. 2002;67(3):945–52. Epub 2002/08/24. pmid:12193406.

73. Decrouy X, Gasc JM, Pointis G, Segretain D. Functional characterization of Cx43 based gap junctions during spermatogenesis. J Cell Physiol. 2004;200(1):146–54. Epub 2004/05/12. pmid:15137067.

74. Vitillo L, Kimber SJ. Integrin and FAK Regulation of Human Pluripotent Stem Cells. Curr Stem Cell Rep. 2017;3(4):358–65. Epub 2017/11/28. pmid:29177133; PubMed Central PMCID: PMC5683053.

75. Fuchs E, Tumbar T, Guasch G. Socializing with the neighbors: stem cells and their niche. Cell. 2004;116(6):769–78. Epub 2004/03/24. pmid:15035980.

76. Morrison SJ, Kimble J. Asymmetric and symmetric stem-cell divisions in development and cancer. Nature. 2006;441(7097):1068–74. Epub 2006/07/01. pmid:16810241.

77. Meng X, Lindahl M, Hyvonen ME, Parvinen M, de Rooij DG, Hess MW, et al. Regulation of cell fate decision of undifferentiated spermatogonia by GDNF. Science. 2000;287(5457):1489–93. pmid:10688798.

78. Huleihel M, Fadlon E, Abuelhija A, Piltcher Haber E, Lunenfeld E. Glial cell line-derived neurotrophic factor (GDNF) induced migration of spermatogonial cells in vitro via MEK and NF-kB pathways. Differentiation. 2013;86(1–2):38–47. Epub 2013/08/14. pmid:23939027.

79. Tadokoro Y, Yomogida K, Ohta H, Tohda A, Nishimune Y. Homeostatic regulation of germinal stem cell proliferation by the GDNF/FSH pathway. Mech Dev. 2002;113(1):29–39. Epub 2002/03/20. pmid:11900972.

80. Tindall DJ, Vitale R, Means AR. Androgen binding protein as a biochemical marker of formation of the blood-testis barrier. Endocrinology. 1975;97(3):636–48. Epub 1975/09/01. pmid:1175511.