|
|
||||||||
Center for Biomedical Research, Population Council, 1230 York Avenue, New York, New York 10021, USA
1 Department of Zoology, University of Hong Kong, Hong Kong, China
(Requests for offprints should be addressed to C Y Cheng; Email: y-cheng{at}popcbr.rockefeller.edu)
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
) that work in concert with proteases, protease inhibitors, kinases, and phosphatases in the microenvironment of the seminiferous epithelium behind the basement membrane. The collective interactions of these molecules determine the timely opening and closing of the bloodtestis barrier (BTB) to facilitate preleptotene spermatocyte migration across the barrier and germ cell movement during the epithelial cycle (Mruk & Cheng 2004b). Additionally, supporting evidence for this local regulatory mechanism is derived from germ cell transplantation studies showing that rat spermatogonia injected into seminiferous tubules in mice differentiate independent of the surrounding milieu (for reviews, see McLean et al. 2001 and Matzuk 2004). Adjudin [1-(2,4-dichlorobenzyl)-1H-indazole-3-carbohy-drazide], formerly known as AF-2364, is a molecule that mediates adherens junction disruption at the Sertoligerm cell interface (Mruk & Cheng 2004b). It is a derivative of lonidamine having a more potent aspermatogenic activity but less toxicity versus lonidamine (Cheng et al. 2001, Grima et al. 2001). When administered to adult rats by gavage, Adjudin exerts its effects primarily at the Sertoligerm cell interface, causing germ cell sloughing, in particular elongating/elongate/round spermatids and spermatocytes from the epithelium without perturbing adhesion between spermatogonia and Sertoli cells; as such, its effects are reversible (Mruk & Cheng 2004b, Cheng et al. 2005). Based on these initial observations, Adjudin has been used to develop an in vivo model to characterize cellcell interactions and junction dynamics pertinent to spermatogenesis (Siu et al. 2003b, 2005, Lee & Cheng 2005, Xia & Cheng 2005, Yan & Cheng 2005, Xia et al. 2006). For instance, the integrin/focal adhesion kinase (FAK)/phosphatidylinositol 3-kinase (PI-3) kinase/extracellular signal regulated kinase (ERK) signaling pathway was shown to regulate Sertoligerm cell adherens junction (AJ) dynamics, particularly the apical ectoplasmic specialization (ES), using Adjudin-treated rat testes (Siu et al. 2003b, 2005). Furthermore, a loss of proteinprotein interactions between integral membrane proteins and their adaptors at the AJ (e.g. apical and basal ES) to facilitate AJ and BTB restructuring were also detected using this model (Lee et al. 2005, Xia & Cheng 2005). This signal pathway activation and the loss of proteinadaptor interactions at the AJ were also demonstrated during spermatid loss from the epithelium, which was induced by suppressing intratesticular androgen level using testosterone and estrogen implants in adult rats (Wong et al. 2005, Xia et al. 2005b, Zhang et al. 2005). Collectively, these data clearly illustrate that the Adjudin model is a valuable tool to identify signaling pathways pertinent to AJ dynamics and possibly the regulatory mechanisms pertinent to spermatogenesis. Since DNA microarray technique has been widely used to unravel global transcriptional changes (for a review, see Stoughton 2005), we sought to identify these potential regulators of junction remodeling pertinent to spermatogenesis using expression microarray. In this report, we describe findings based on the use of Affymetrix Genechips (rat genome) that contain ~30 000 probe sets to characterize the expression profile in rat testes following treatment with Adjudin at the time of AJ restructuring. The genes and the signaling conduits identified by microarray could provide a framework to further probe the biological processes of junction restructuring pertinent to spermatogenesis.
| Materials and Methods |
|---|
|
|
|---|
Male SpragueDawley rats (~300 g b.w.) were purchased from Charles River Laboratories (Kingston, MA, USA). The use of animals in this study was approved by The Rockefeller University Animal Care and Use Committee, with protocol numbers 03017 and 06018. Rats were treated with a single dose of Adjudin at time 0 by gavage at 50 mg/kg body weight (b.w.) as described earlier (Cheng et al. 2001, Grima et al. 2001). For microarray analysis, rats (n = 3 for time 0 and 96 h; n = 2 for 8 h) were killed by CO2 asphyxiation. These two time points were selected since by 8 h, even though some elongate and elongating spermatids had begun to exfoliate from the seminiferous epithelium, most germ cells remained attached to Sertoli cells in the seminiferous epithelium. Even for those that were detached from the epithelium, they were still found in the tubule lumen. As such, this is the time that extensive AJ restructuring was occurring at the Sertoligerm cell interface, but germ cells remained in the samples being analyzed by microarray. By day 4, virtually all elongating/elongate spermatids were depleted from the epithelium and most round spermatids and spermatocytes had left the epithelium and tubule lumen as well, suggesting the induced AJ restructuring had subsided as shown in earlier studies (Chen et al. 2003, Mruk & Cheng 2004a,b). Testes were frozen immediately in liquid nitrogen upon termination and stored at 80 °C until use for RNA extraction. Rat Genome 230 2.0 GeneChips were purchased from Affymetrix, Inc. (Santa Clara, CA, USA) and used for this study.
Total RNA preparation
About 0.1 g testis from each sample was used for total RNA isolation using TRIZOL (Invitrogen). RNA was then purified with an RNeasy cleanup kit (Qiagen). Total RNA quality was assessed with an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA).
cRNA labeling, hybridization, and array scanning
cDNA synthesis and cRNA labeling by in vitro transcription (IVT) were performed using GeneChip One-Cycle Target Labeling and Control Reagents (Affymetrix, P/N 900493). Briefly, total RNA (~5 µg) prepared as mentioned above was first reverse transcribed using a T7-oligo (dT) promoter primer in the first-strand cDNA synthesis reaction. Following RNase H-mediated second-strand cDNA synthesis, the double-stranded cDNA was purified and served as a template in the subsequent in vitro transcription reaction. The IVT reaction was carried out in the presence of T7 RNA polymerase and a biotinylated nucleotide analog/ribonucleotide mix for cRNA amplification and biotinylation following the protocols provided by the manufacturer. The biotinylated cRNA targets were then purified using RNeasy spin columns, and fragmented at 94 °C for 35 min in a fragmentation buffer (40 mM Tris-acetate (pH 8.1) containing 100 mM potassium acetate and 30 mM magnesium acetate), which were used for hybridization with GeneChips. The integrity of the labeled cRNA was again assessed by Agilent 2100 Bioanalyzer before hybridization. The GeneChip Fluidics Station 450 was used for array hybridization and processing and the Affymetrix GeneChip Scanner 3000 for array scanning. cRNA labeling, hybridization, and array scanning were performed at The Rockefeller University Genomics Resource Center (New york, NY, USA). The experiment was repeated twice using different samples.
Microarray data analysis
Scanned images of the microarray were first visually examined to rule out physical anomalies or excessive high background. Data were initially accessed using the Microarray Suite 5.0 (MAS) software (Affymetrix) and all probe sets on each array were scaled to a target signal intensity of 250. This prescaling assured data uniformity in subsequent cross-chip comparison. MAS assessment also reported the background noise, the number of probe sets that are present, marginal, or absent in the samples out of a total of 31 099 probe sets on Rat Genome 230 2.0 GeneChip based on signals. Data were then exported and analyzed using the GeneSpring software (Version 7.2, Agilent). All samples were grouped accordingly into time 0 (control), 8-h and 4-day post-Adjudin treatment. Data were normalized in three steps: step 1, data transformation to set measurements from <0.01 to 0.01. Step 2, per chip normalization to the median or 50th percentile intensity of the probe set in each chip to remove non-biological variations such as differences in the experiment conditions. Step 3, per gene normalization to normalize all samples to time 0 (the median value for each probe set among the three control samples was used). Data were further filtered using the filtering on flags feature (namely signal detection, present, marginal or absent; at least present or marginal in one of the sample after Adjudin treatment) in Genespring to identify genes that are expressed in one of the experimental conditions under study (e.g. Adjudin treatment) and expression raw value (at least 50) to remove unreliable genes. Thereafter, genes that were expressed differentially and were statistically significantly different from control (normal rat testes, P< 0.05) were determined by one-way ANOVA for multiple comparisons between samples in control (n = 3) and treatment groups (n = 2 for 8 h or n = 3 for 4-day post-Adjudin treatment) with multiple test correction (default option) and multiple test correction set to none, using the Agilent GeneSpring software. Genes with significant changes in their steady-state mRNA levels in comparison with control were scored and grouped. Hierarchical clustering analysis (unsupervised) was performed in the GeneSpring software using default settings. Gene ontology and signaling pathway analysis was performed using Onto Express and Pathway Express programs available from Wayne State University (Draghici et al. 2003a,b, Khatri et al. 2004, 2005). The Affymetrix probe set ID corresponding to the specific gene at the GeneBank database with fold changes in the steady state mRNA levels at 8 h or at day 4 were imported into the programs and the top-ranking pathways were identified.
Quantitative real-time RT-PCR and data analysis
To validate relative expression levels of selected genes from microarray studies, quantitative real-time RT-PCR (qPCR) was performed using an Applied Biosystems Prism 7700 Sequence Detection System with SYBR Green PCR Master Mix (Applied Biosystems, Inc., Foster City, CA, USA). Primers for qPCR were designed using either the Oligonucleotide Properties Calculator at www.basic.northwestern.edu/biotools/oligocalc.html or the Primer Express (Version 2.0) from Applied Biosystems, Inc. and were compared with the existing database at GenBank using the basic local alignment search tool (BLAST) to ensure specificity (see Table 1
). Briefly, total RNA from different samples were prepared using TRIZOL (Invitrogen) and treated with DNaseI (Invitrogen) to remove contaminating genomic DNA before reverse transcription. Equal amount of total RNA was then reverse transcribed using Moloney murine leukemia virus reverse transcriptase (M-MLV RT, Promega) with 2 µg RNA in a 25 µl reaction. Both primers (5 pmol each) and various templates (serially diluted control testis cDNA for generating the standard curve, and 1 µl Reverse transcription product from all samples) were mixed with 2 x SYBR Green PCR Master Mix in a 25 µl final reaction volume. The thermal cycler conditions were 10 min at 95 °C, followed by two-step (15 s at 95 °C and 1 min at 60 °C) PCR for 40 cycles, and a slow heating up to 95 °C for dissociation analysis at the end to ensure the purity of the PCR product. Amplification data were analyzed with an Applied Biosystems Prism Sequence Detection Software (Version 2.1). The cycle threshold (CT) value corresponding to the PCR cycle number at which fluorescence emission in real time reached a threshold above the baseline emission was determined. The CT value for each sample was then exported and analyzed in Microsoft Excel. A combination of standard curve and comparative CT methods were applied for the analysis (User Bulletin #2 ABI PRISM 7700 Sequence Detection System, Applied Biosystems). For standard curve method, CT for Stat3 and S-16 were plotted as CT value against logarithm of relative input RNA amount respectively, and then the input RNA amount for each sample was calculated based on its CT value in reference to the corresponding standard curve. Relative expression of Stat3 was normalized to the S-16 control for each sample. For comparative CT method, the difference of CT between the target gene and S-16 for each sample was calculated (
CT) and then this difference for other samples was compared with that of a particular sample (
CT). The expression level in any sample relative to that particular sample was given by 2
CT. All reactions were prepared in triplicate and three independent sets of samples were used in each experiment. Results of qPCR were analyzed by one-way ANOVA for multiple comparisons using the JMP IN statistical analysis software package (Version 4, SAS, Inc., Cary, NC, USA).
|
Immunohistochemistry (IHC) for selected genes was performed essentially as previously described (Xia & Cheng 2005). Antibodies used in this study were either raised in this laboratory (rabbit anti-testin (Cheng & Bardin 1987), anti-cathepsin L (Lee et al. 1986),
2-MG (Cheng et al. 1990)) or purchased from Zymed-Invitrogen (Carlsbad, CA, USA) (rabbit anti-claudin 3 polyclonal antibody). Immunoblot analysis was performed as described earlier (Xia & Cheng 2005). Testis lysates from samples were prepared and resolved by SDS-PAGE and probed against a panel of antibodies.
Assessing changes in the steady-state mRNA levels of target genes in testes taking into account cellular composition alteration in the seminiferous epithelium following Adjudin treatment a mathematical model
Adjudin treatment induces massive germ cell depletion from the seminiferous epithelium, with initial loss of elongating/elongate spermatids, to be followed by round spermatids and spermatocytes (Chen et al. 2003; Fig. 1AC
). As such, the cellular composition in the epithelium changes with time following treatment, particularly between days 1 and 4, during germ cell sloughing. For instance, if a protein is exclusively expressed by Sertoli cells, such as occludin, its induction that was detected by RT-PCR and/or immunoblotting following Adjudin treatment versus controls by day 4 may simply reflect an increase in cellular RNAs or proteins contributed by Sertoli cells in the samples being analyzed because of declining germ cell numbers as a result of germ cell loss. This thus makes it difficult to assess if induction of a target protein is in fact related to junction restructuring or the result of cellular composition alteration. Herein, we provide a mathematical model to relate changes in cellular composition and an alteration of expression level of a target protein/gene so as to provide a threshold by which one can accept or reject data generated in this Adjudin or similar animal models. Because non-germ cell types in the testis (e.g. Sertoli, Leydig, and myoid cells) do not have any significant changes after Adjudin treatment (Cheng et al. 2001, 2005, Grima et al. 2001), we collectively assigned these cell types as somatic cells and categorized them under the Sertoli cell item as an entity. We also did not discriminate germ cells at different developmental stages and categorized them under the germ cell as an entity. This model was being used to screen target proteins/genes that are pertinent to junction restructuring by accepting their induction in the microarray experiments as physiologically relevant.
|
Control (normal) testis status (A) is : p + q = 1
Testis status after Adjudin treatment (B) is : p' + q' = m
where p, p' are the corresponding Sertoli cell fractions before and after Adjudin treatment, and q, q' are the corresponding germ cell fractions before and after Adjudin treatment (assuming p = p' since Sertoli cells were not affected by Adjudin; and 0<p
m<1, 0<q'
m, q>q').
For a target gene of interest, such as Z, its expression level in A (total) is:
![]() |
and Z expression level in B (total) is:
![]() |
where
,
' are the corresponding expression coefficients for gene Z in Sertoli cells before and after Adjudin treatment in A or B, and ß, ß' are the corresponding expression coefficients for gene Z in germ cells and after Adjudin treatment in A or B with 0
,
', ß, ß'. The expression coefficient of a target gene is the total number of copies of gene Z mRNA in Sertoli (
) or germ (ß) cells. For instance, in normal testes occludin is expressed in Sertoli cells only. Therefore,
occludin is a positive value and ßoccludin = 0 (note: the exact value for
occludin is not required for the following calculations, see below).
From the microarray experiments, changes in the expression level of a target gene Z for testis status B versus control status A were determined as x. Because the same amount of total RNA has used in treatment versus control samples, i.e. from the same unit weight of testis tissue (assuming the same unit weight of tissue from different samples contain the same amount of total RNA), the following formula was used:
![]() |
We consider the following three possibilities regarding target gene/protein changes during Adjudin-induced germ cell loss from the seminiferous epithelium:
1. If gene Z is only expressed in Sertoli cells: ß = 0, ß' = 0. Thus, ZB/ZA =
'/
= mx. For example, 4 days after Adjudin treatment, m (weight of testis versus control) = 0.73 (see Fig. 1Ad
), then if
'/
>1, x should be >1/0.73
x
1.4. This indicates only those genes with
1.4-fold of change are considered to be induced in the testis after Adjudin treatment. We set the threshold x
2, thus
'/
= mx
0.73 x 2
1.5. For genes listed in Table 4
, the threshold is actually x
4 by day 4, meaning
'/
0.73 x 4
3.
|
= 0,
' = 0
Thus,
(since m<1).
In this case, the actual fold changes in gene induction in germ cells are greater than the reported fold of changes for x.
3. If gene Z is expressed in both Sertoli and germ cells:
![]() |
Herein,
![]() |
where min [(p
')/(p
), (q'ß')/(qß)] is the minimal value of the two and max[(p
')/(p
), (q'ß')/(qß)] is the maximal value of the two (this is because p
', q'ß', p
and qß are all positive values. When a, b, c, d>0, (a/c)
(b/d), then (a/c)
(a+ b)/(c+ d)
(b/d). Therefore, either (p
')/(p
)
mx or (q'ß')/(qß)
mx.
When we select the threshold for x such that mx
1 (at day 4, x
2 and mx
1.46), then either p
'/p
mx>1 (
'/
1, i.e. induction in Sertoli cells) or q'ß'/qß
mx>1 (
ß'/ß
q/q'>1, i.e. induction in germ cells). On the other hand, the possibility in (3) can be treated as a mixture of (1) and (2), and on average, ZB/ZA should be between the values in (1) and (2) as well. Whichever the possibility, the threshold x
2 translates into a minimal 1.5-fold induction in either germ or Sertoli cells. Table 4
lists those target genes with fourfold of induction (x
4) by day 4 showing Adjudin treatment, which translates into a minimal threefold of induction in either germ or Sertoli cells.
We selected target gene Myd116 for the sake of discussion, which was shown to be induced 3.5-fold at 8 h and 1.7-fold at day 4 (see Table 3
). At 8 h, 3.5 x 0.9 = 3.15 and at day 4, 1.7 x 0.73 = 1.241. These results indicate that (1) if this gene is exclusively expressed in Sertoli cells (case 1), there was a 3.15-fold induction at 8 h (significant) and only 1.241-fold induction at day 4 (not likely to be significant); (2) if this gene is exclusively expressed in germ cells (case 2), there was more than 3.5-fold induction at 8 h (significant) and more than 1.7-fold induction at day 4 (not certain); (3) if this gene is expressed both in Sertoli and germ cells (case 3), there was an overall 3.15- to 3.5-fold induction at 8 h (significant) and 1.241- to 1.7-fold induction at day 4 (not certain). Collectively, this gene is significantly induced at 8 h and not likely to be induced significantly at day 4 under all possibilities. Thus, the twofold threshold is a reliable filter to accept or reject any data to assess if a target gene has been induced. Nonetheless, additional cautious steps have to be taken since these data should be validated by techniques of IHC, fluorescent microscopy, and immunoblottings.
|
| Results |
|---|
|
|
|---|
The original data sets generated in this report have been deposited at the Gene Expression Omnibus (GEO) data repository website: http://www.ncbi.nih.gov/geo/ with accession number: GSE5131 [NCBI GEO] .
Adjudin-induced morphological changes in adult rat testes
The rationale of selecting the 8-h and 4-day time points post-Adjudin treatment for their comparison with control (normal testes) was based on an earlier report that investigated the kinetics of germ cell loss from the epithelium (Chen et al. 2003). By 8 h after Adjudin treatment, while more than half of the seminiferous tubules examined (n = 600 from four testes with ~125 tubules randomly scored per testis) displayed signs of germ cell loss with
10 elongating/elongate spermatids found in tubule lumen, most elongating spermatids, round spermatids, and spermatocytes still remained attached to Sertoli cells in the epithelium (see Fig. 1Ab
vs Fig. 1Aa
). This observation also illustrates that even though there was extensive restructuring at the Sertoligerm cell interface and depleting elongating/elongate spermatids were found in the lumen, the cellular composition as well as RNAs contributed by Sertoli and germ cells in the samples being analyzed by microarrays were similar to control testes. This argument is supported by the finding that only a ~10% decline, which was not statistically significant, in testicular weight was detected by 8 h (see Fig. 1Ad
). By 4 days after Adjudin treatment, about 80% of the tubules that were scored were shown to have the number of round spermatid layers reduced by twofold or more with virtually all elongating/elongate spermatids depleted from the epithelium, and most of the tubules were only populated with round spermatids, spermatocytes, and spermatogonia along with Sertoli cells (Chen et al. 2003; Fig. 1Ac
vs Fig. 1Aa
, and b). These results were consistent with a significant loss of testicular weight, by ~27% at day 4 (Fig. 1Ad
). Due to the changes in cellular composition in the tubules as a result of progressive germ cell loss, we sought to correct the expression levels of genes detected by microarray studies by normalizing the results by the testis weight loss since total RNA from whole testes was used for microarray analysis. In our analysis, a threshold was set at twofold, which translated into ~1.5-fold of net changes (2 x 73%) for 4-day samples and ~1.8-fold of net changes (2 x 90%) for 8-h samples, so as to set the basis to accept that a target gene had been induced following Adjudin treatment. A detailed quantitative mathematical analysis of gene expression/protein level changes in concert with cellular composition alteration in the seminiferous epithelium is described in Materials and Methods. Based on this analysis, a twofold induction of a target gene/protein expression has exceeded the threshold that could be simply accounted for by cellular composition alterations in the epithelium by 8 h 4-day post Adjudin treatment, which would translate into, at the maximum, a 1.5-fold change for a target gene. As such, the focus of this study was placed on genes that were induced by twofold or more in particular by day 4. As expected, many genes appeared to be down-regulated in 4-day samples when most germ cells were depleted from the epithelium, illustrating that these genes were exclusively or abundantly expressed by germ cells. As such, these data (see Supplementary Table 1
and raw data deposited at GEO website), also offer the means to identify germ cell-specific genes.
Assessment of microarray data
The Affymetrix GeneChip uses external RNA controls, also known as spikes, or spike-in controls, to assess the inter-microarray experimental variations (Bakel & Holstege 2004, Stoughton 2005). Statistical analysis was performed by comparing data obtained from normal (control) rat testes versus treatment groups by one-way ANOVA for multiple comparisons with multiple test correction (default option), and multiple test correction was set to none using the GeneSpring software (see Materials and Methods). The signals for these controls, as well as for housekeeping genes (e.g. glyceraldehyde-3-phosphate-dehydrogenase (GAPDH) and ß-actin) were consistent between chips (i.e. between different samples) in a single experiment, and were also consistent between samples from different experiments. For testis samples from control rats (n = 3) and 8 h (n = 2) after Adjudin treatment, about 55% of 31 099 probe sets were positively detected (i.e. present), about 500 (~1.6%) were absent and 44% probe sets were detected as marginal, (Fig. 1Ae
). However, samples from rats on day 4 (4 days, n = 3) after Adjudin treatment have shown an increase in the number of genes present (~59%), and a reduction in the number of genes marginal (~40%), whereas the number of genes that were absent remained unchanged at ~1.6% (Fig. 1Ae
). These results indicated that ~4% of genes in the rat genome or ~1240 probe sets became detectable following a longer period of Adjudin treatment (by day 4) and when the epithelium was devoid of virtually all elongating/elongate spermatids, and ~4050% round spermatids and ~1020% spermatocytes (see Fig. 1A
). Indeed, an examination of all probe sets by plotting the normalized expression levels at 8 h or at day 4 versus their corresponding raw expression levels or signal intensities in control samples revealed a more dispersed distribution for day 4 than 8 h (Fig. 1Af
vs g). In particular, there were more probe sets at day 4 than at 8 h that fell above or below the twofold threshold lines with signal intensities >100, indicating more genes were differentially regulated at the later time point after Adjudin treatment. We filtered out 1466 transcripts which were significantly regulated after Adjudin treatment (control, 8 h and day 4) by one-way ANOVA. There were many more transcripts that were up- or down-regulated by at least twofold in 4-day samples than in 8 h samples (Fig. 1B
), consistent with the data shown in Fig. 1Af and g
. There were 79 transcripts showing more than twofold change in expression; within this group, 31 had more than twofold change detectable only at 8 h but not at day 4, 48 transcripts were regulated by more than twofold both at 8 h and at day 4 (Fig. 1B
). However, 436 transcripts showed more than twofold change only at day 4 but not at 8 h (Fig. 1B
). Selected genes belonging to these three categories (
twofold change at 8 h only;
twofold change at day 4 only; and
twofold change at both time points) are summarized in Tables 3
5![]()
. Table 4
summarizes the partial list of 436 transcripts with more than fourfold induction. The full list can be found in Supplementary Table 1
(see supplementary data in the online version of the Journal of Endocrinology at http://joe.endocrinology-journals.org/content/vol192/issue3/).
|
|
Genes with more than twofold change in their expression levels were tabulated and grouped by their general functions (Tables 3
5![]()
). Among these are cytokines, transcription factors, kinases, phosphatases, genes that regulate ion homeostasis and transport, proteases, protease inhibitors, cell adhesion molecules, and signaling molecules. These tables listed the top-ranked genes from classifications defined in Fig. 1B
and these top-ranked genes fell into the aforementioned functional groups, indicating these genes are likely involved in the event of germ cell depletion from the seminiferous epithelium induced by Adjudin.
In the cytokine group, interleukin (IL) and chemokines are most common in early and late time points after Adjudin treatment. IL-1
and Cxcl10 showed strong induction by 8 h, but their expression level receded by day 4 (Table 3
). Cxcl12 and Il17re were only induced by day 4 (Table 4
). Cxcl1, Ccl2, and Cx3cl1 showed consistent up-regulation at both time points (Table 5
).
A number of transcription factors were significantly induced after Adjudin treatment. Some were induced mostly at 8 h, which included Egr1 and Egr2, Bhlhb2, and Jun (Table 3
). At day 4, Nupr1 and Pawr became significantly up-regulated versus control (Table 4
). Transcription factors, Myc, Stat3, Atf3, and c-fos, also showed strong induction at 8 h and at day 4 (Table 5
). Transcription factors that were induced during Adjudin-induced germ cell loss from the testis as shown in Tables 3
5![]()
were summarized in Table 6
along with their known function in the testis, illustrating these may be potential regulators of junction restructuring events at the Sertoligerm cell interface, pertinent to spermatogenesis.
|
We selected a number of target genes to validate their expression changes detected by microarray using other methods as follows: (1) transcription factors: Stat3, Jun, Bhlhb2, and c-fos (Fig. 2A and B
); (2) cytokines and related molecules: IL-1
, Cxcl10, and Tnfrsf1a (Fig. 2C
); and (3) proteases and protease inhibitors: Prss11 and Serping1 (Fig. 2D
). The fold changes determined by qPCR were consistent with microarray data for all nine genes selected for analysis (Fig. 2
vs Tables 3
5![]()
). For instance, both Prss11 and Serping1 showed relatively little change at 8 h but were induced by about four- to fivefold at day 4 (Table 4
) in microarray analysis, and results from qPCR showed that at 8 h, the steady-state mRNA levels of these two genes were not significantly different from controls but were induced by approximately fivefold at day 4 (Fig. 2D
). Since these nine genes were randomly selected from Tables 3
5![]()
for analysis, it is logical to conclude that the fold change detected for other genes listed in these tables is reliable. Furthermore, the changes in the steady-state mRNA levels of cytokines reported in Tables 3
5![]()
, such as Cxcl1, were confirmed using cytokine antibody array (data not shown).
|
2-Macroglobulin (
2-MG) was found to be induced by 1.1-fold at 8 h and threefold at day 4 (Supplementary Table 1
|
We performed unsupervised hierarchical clustering of the eight samples (Ctrl (n = 3), 8 h (n = 2), and day 4 (n = 3)) and the 1466 transcripts using the GeneSpring software and the results were summarized in Fig. 4
. The branches of the dendrogram were removed from this presentation, but the main subtrees were illustrated by a line showing five different domains (Figure 4ae
), which were again classified into two larger subtrees of i and ii. This clustering analysis grouped similarly regulated genes. Two representative regions in the dendrogram were selected and examined in detail in Fig. 4B and C
, which showed the induced expression of these genes at day 4 but not at 8 h, and reduced expression at 8 h but mildly increased expression at day 4 respectively.
|
We next used Pathway Express (PE) to identify the signaling pathways that were activated during Adjudin-induced germ cell loss from the epithelium (Draghici et al. 2003a,b, Khatri et al. 2004, 2005; Table 7
). PE is a tool in the Onto-Tools ensemble (Wayne State University, see http://vortex.cs.wayne.edu (Draghici et al. 2003a,b, Khatri et al. 2004, 2005)) that was designed to identify the signaling pathways that contain the corresponding target genes. A list of genes with their fold changes were first submitted to the system, which then performed a search and built a list of all associated pathways. The system calculated an impact factor to rank the pathways based on the importance of a gene in the pathway (indicating how much it affects the downstream signaling), the strength or fold of change of the gene, and percentage of genes in the pathway that were included in this list of genes (Khatri et al. 2004, 2005). The probe set ID number with the corresponding fold changes at 8 h or at day 4 were tabulated (all 1466 probe sets) and imported into the program and the top signaling pathways were identified. At 8 h, the top-ranking pathways identified based on the impact factor were phosphatidylinositol (PI) signaling (7/43: input genes/pathway genes), focal adhesion (13/109), cytokinecytokine receptor interactions (16/110), MAPK signaling pathway (21/153), tight junction (10/185), apoptosis (11/66), and calcium signaling pathway (13/161; Table 7
). At day 4, the top-ranking pathways were phosphatidylinositol signaling, focal adhesion, cytokinecytokine receptor interactions, tight junction, MAPK signaling pathway, apoptosis, and calcium signaling pathway (Table 7
). The numbers of input and pathway genes for 8 h and day 4 were the same. Collectively, these results were consistent with recently published findings that these are putative signaling pathways that regulate Sertoligerm cell adhesion (Lui et al. 2003a, Siu et al. 2003b, 2005, Wong et al. 2005, Xia & Cheng 2005, Zhang et al. 2005) and probably suggest that these are the most important signaling events taking place in the testis after the rats were treated with Adjudin.
|
| Discussion |
|---|
|
|
|---|
This microarray study has identified genes that are differentially regulated during Adjudin-induced AJ restructuring at the Sertoligerm cell interface prior to the massive physical dislodging of germ cells, at 8 h vs day 4, when virtually all elongating/elongate spermatids and most round spermatids and spermatocytes were depleted from the epithelium, and compared with normal testes. Among these genes, we summarize those with most drastic induction in the steady-state mRNA level versus controls and categorize them into functional groups (Tables 3
5![]()
). Although only the expression data were obtained using this genome-wide gene profiling approach without considering translational and post-translational regulations, these genes nonetheless are potentially key players in the junction restructuring events at spermatogenesis since recent studies using specific inhibitors and the other models involving steroid implants to suppress endogenous testosterone level in the testis have clearly demonstrated that data obtained from this model is relevant to the in vivo physiology of Sertoligerm cell interactions, particularly cell adhesion regulation (Lui et al. 2003a, Siu et al. 2003b, 2005, Siu & Cheng 2004, Lee & Cheng 2005, Lee et al. 2005, Wong et al. 2005, Xia & Cheng 2005). The identified genes include those that are involved in signaling and/or regulatory functions pertinent to the junction restructuring: (1) cellcell communication, including cytokine-mediated signaling, ion transport and homeostasis; (2) intracellular signaling mediated via kinases, phosphatase, and adaptors; (3) structurally important cell adhesion and cytoskeleton-related proteins; (4) proteases and protease inhibitors that affect cell matrix, cellcell and intracellular interactions; and (5) transcription factors. Except for the last group and ion transport/homeostasis proteins, genes belonging to all other groups have been investigated and shown to be involved in regulating junction restructuring in earlier studies (Mruk & Cheng 2004b). We had since proposed a junction restructuring theory to describe the mechanisms of germ cell movement in the seminiferous epithelium (Mruk & Cheng 2004b). This theory was formulated based on studies of in vivo (e.g. Adjudin model, cadmium model), and in vitro (e.g. Sertoligerm cell co-cultures) which have also identified several putative signaling pathways that are known to regulate Sertoligerm cell adhesion (Lui et al. 2003a, Siu et al. 2003b, Lee et al. 2005), which are also confirmed in this microarray study. Perhaps the most important of all, results of this gene profiling study have unraveled two important lines of research pertinent to Sertoligerm cell adhesion regulation. First, additional mechanistic pathways that can be potentially crucial to Sertoligerm cell adhesion have been identified. For instance, this gene profiling study has demonstrated that PI-3 kinase, FAK, and MAPK are the leading signaling molecules that are involved in the Adjudin-induced Sertoligerm cell AJ restructuring, confirming results of earlier studies (Lui et al. 2003a, Siu et al. 2003b, 2005, Wong et al. 2005, Xia & Cheng 2005, Zhang et al. 2005, Xia et al. 2006). However, cytokinecytokine receptor interactions, besides TGF-ß3, TNF
, and their corresponding receptors, IL-1
and its receptors, have been identified to be a crucial mechanism to AJ dynamic regulation. This suggests that an array of cytokines produced by Sertoli and germ cells into the microenvironment at the seminiferous epithelium such as the blood testis barrior, are working in concert to activate their receptors in tandem to regulate junction restructuring events pertinent to spermatogenesis. Furthermore, the calcium signaling pathway and apoptosis regulatory genes are both activated, which should be carefully evaluated in future studies. Secondly, a set of transcription factors have now been identified, some of which can potentially regulate TJ and/or AJ dynamics, which should be carefully evaluated in future studies. It will also be crucial in future studies to examine if the downstream action of PI-3 kinase, FAK, and MAPK are mediated via one or several of these transcription factors.
Does the Adjudin model reveal physiologically relevant information on the biology of spermatogenesis?
Besides building expression profiles of thousands of genes in the testis, results from this microarray study could be used to identify genes pertinent to other testicular function. Indeed, many transcripts were found to be significantly induced after multiple rounds of filtering. We did not explore the potential functions of these new genes in this study but they are summarized in Supplementary Table 2
(http://joe.endocrinology-journals.org/content/vol192/issue3/). At present, a number of expression profile studies in the testis are available in the literature (Schlecht et al. 2004, Small et al. 2005, Zhou et al. 2005) pertinent to development, spermatogenesis, and/or fertility regulation. The present study, however, uses an exogenic compound and examines its effect in the testicular gene expression during the time course of germ cell loss with extensive restructuring at the Sertoligerm cell interface. An important question remains to be addressed: does this drug model reveal the molecular targets pertinent to spermatogenesis?
The approach of using the Adjudin model to identify molecular targets pertinent to junction remodeling during spermatogenesis may be criticized that the changes in the expression of multiple genes are the results of drug toxicity. However, recently completed toxicity studies have shown that the dose used in the present study that was effective to induce germ cell loss had no toxic effects in rats or mice (Cheng et al. 2005, Mruk et al. 2006). Moreover, Adjudin-induced germ cell depletion is reversible (Cheng et al. 2001, Grima et al. 2001). The Adjudin effects on a number of genes seemingly support the notion that this model offers a snapshot of expression alterations that represent their changes during spermatogenesis. For instance, cathepsin L becomes highly expressed at stages VIIVIII (Erikson-Lawrence et al. 1991, Anway et al. 2004), which coincides with spermiation and BTB restructuring, and is also induced during Adjudin-induced germ cell loss from the epithelium, mimicking the cellular events of spermiation. TGF-ß3, which was shown to disrupt both Sertoligerm cell adhesion and BTB integrity and was known to be stage-specifically regulated with the highest expression at stages VIVIII (Xia & Cheng 2005, Xia et al. 2006), is also induced after Adjudin treatment at the time of germ cell loss, possibly because of AJ restructuring and this observation is consistent with the results of an earlier report (Xia & Cheng 2005).
Thus, we hypothesize that the Adjudin-induced restructuring at the Sertoligerm cell interface mimics the anchoring junction restructuring events in the epithelium during spermatogenesis. This postulate was supported by recent findings that the signaling molecules crucial to Sertoligerm cell AJ dynamics as identified using the Adjudin model, such as p-FAK, p-c-Src, rMTMR2, and pERK (Li et al. 2000, Siu et al. 2003b, 2005, Xia & Cheng 2005) have also been found in another in vivo model using androgen/estrogen implants to suppress the intratesticular androgen level which is known to perturb Sertolispermatid adhesion in adult testes (ODonnell et al. 2000, Wong et al. 2005, Xia et al. 2005b, Zhang et al. 2005). This is plausibly due to the fact that there are only a limited number of signaling pathways and molecular mechanisms that control a specific phenotype. For instance, if classified by the types of cell surface receptors involved, the number of signaling pathways in human cells can be as small as 16, and as many as 200 only if counting all the variants of receptors (Fishman & Porter 2005). Normal spermatogenesis and the Adjudin-induced germ cell depletion model share some phenotypic features, such as junction restructuring, and it is not unusual that they utilize similar signals and/or regulatory mechanisms. This hypothesis should be carefully evaluated in future studies. Likewise, results using Pathway Express to rank the top signaling pathways that were activated by Adjudin in testes also support this postulate. Lipid and protein kinases, such as PI-3 kinase, are important transducers to relay signals from matrix to cells, regulating cytoskeletal remodeling, and cell migration (Carragher & Frame 2004, Webb et al. 2004). It is not surprising to find that phosphatidylinositol signaling and focal adhesion kinase signaling were the top two pathways. Indeed in studies using the Adjudin model, Sertoligerm cell co-cultures, and the androgen-suppression model have shown that integrin ß1, PI-3 kinase, and its downstream protein kinases are crucial regulators of restructuring at the apical ES (Siu et al. 2003b, 2005, Wong et al. 2005). Additionally, MAPKs have also been identified in earlier studies as important regulators of junction dynamics in the testis (Lui et al. 2001, 2003b,Lui et al. c, Wong et al. 2004, Wong & Yan Cheng 2005, Xia & Cheng 2005, Xia et al. 2005a, 2006). Taken collectively, these results provide strong evidence that the Adjudin model is useful to study the junction restructuring events pertinent to spermatogenesis.
Molecular signatures for the regulation of junction restructuring during spermatogenesis
It is of interest to note that junction restructuring and migration of different cell types share similar underlying regulatory mechanisms, since these events are fundamentally conserved across species and cell types regarding the molecules that are involved in these cellular processes. For instance, results from this study have pinpointed the important players in these events in the testis (Tables 3
5![]()
). Unsurprisingly, these results were in agreement with other genomic-scale profiling studies which have also identified genes belonging to kinases, proteases, cell adhesion molecules, and cytokine signalings are enriched in migratory cells (Borghese et al. 2006, Wang et al. 2006). We herein discuss target gene groups identified in this Adjudin model, which may be crucial to junction remodeling during germ cell migration.
Cytokines and cytokine receptor signalings
Besides TGF-ß3 (Xia & Cheng 2005), Adjudin indeed activated the expression of several cytokines and chemokines in adult rat testes. ILs and chemokines are typically involved in inflammatory responses and they are also crucial to T-cell migration (Petersen et al. 2006). Other studies have shown that cytokines are important players of immunological responses in the testis after stress or injury (Hedger & Meinhardt 2003). However, as no macrophage infiltration in the testis was observed after Adjudin treatment, it is logical to conclude that the inflammation response to Adjudin is probably a negligible possibility. The activation of multiple cytokines and chemokines at the time of germ cell loss implies that these paracrine/autocrine factors are involved in the signaling cascades that lead to the restructuring at the Sertoligerm cell interface. For instance, IL-1
is a pro-inflammatory cytokine but in the testis it is also secreted by Sertoli cells in the absence of local inflammation (Sultana et al. 2004). An IL-1
isoform was found to be secreted bidirectionally to the interstitium and intratubular compartment (Sultana et al. 2004). As a paracrine factor, it regulates Leydig cell steroidogenesis with age-dependent effects (Sultana et al. 2004). By 8 h after Adjudin treatment, a drastic induction (8.6-fold) of IL-1
was detected, which strongly suggests that its early response signifies its significance in AJ restructuring, which should be investigated in future studies.
Cx3cl1 (or fractalkine) is also expressed in Leydig cells, Sertoli cells, and spermatogonia and it is believed to be regulating cellcell interactions in and out of the seminiferous tubules and other epithelia (Rossi et al. 1998, Habasque et al. 2003). In the Adjudin model, a 3.2- and 5.3-fold increase in expression at 8 h and at day 4 respectively, were detected. Vgf (VGF nerve growth factor inducible) was found to be highly induced after Adjudin treatment (3- and 16-fold increase at 8 h and at day 4 respectively), although its basal expression level in the rat testis was relatively low (raw signal intensity at ~30). There is no literature on Vgf in the testis, since Vgf is a neuron-secreted polypeptide, particularly abundant in adult hypothalamus. Vgf/ mice were found to be small, hypermetabolic, hyperactive, and infertile for both males and females (Hahm et al. 1999), which is likely the result of a reduced production of pituitary gonadotropin hormones (Hahm et al. 1999), which in turn, affected spermatogenesis in the testis. In short, these cytokines should be examined closely to determine if they play any role in Sertoligerm cell AJ restructuring during spermatogenesis.
Proteases and protease inhibitors
Cathepsin L and tissue-type plasminogen activator (tPA) are known to be expressed stage specifically in the testis (Vihko et al. 1989, Chung et al. 1998). Cathepsin L and urokinase-type plasminogen activator (uPA) are also induced in the cadmium model (Wong et al. 2004), supporting the notion that they take part in the regulation of junction restructuring. Similarly, tissue inhibitor of metalloproteases 1 (TIMP-1) was also found to regulate BTB and AJ dynamics (Mruk et al. 2003, Siu et al. 2003a), but its activation in the Adjudin model has not been reported until now (see Table 4
). The function of other proteases and protease inhibitors, such as Prss11 and Serping1, in these events are not known. For instance, no reports are found in the literature regarding the expression of Spin2c or Serping1 by the testis, but they are highly induced after Adjudin treatment (Table 4
). It will be of interest to examine the regulation of these proteins in the testis during the epithelial cycle and their role in AJ restructuring.
Transcription factors
Transcription factors that were induced by Adjudin in the testis were summarized in Table 6
. For example, early growth responses 1 and 2 (Egr1 and Egr2) were greatly induced 8 h after Adjudin treatment (Table 3
). They belong to the Egr family of zinc-finger <