|Home | About | Journals | Submit | Contact Us | Français|
Y-chromosomal microdeletions (YCMD) are the major genetic cause of male infertility. To date, it is not known which global changes are induced by the presence of AZFc or AZFb + c deletions in the human testicular transcriptome. We investigated this question by microarray analysis in which we had to eliminate the ‘germ cell effect’, i.e., the dominating effect of germ cell transcripts due to the quantitative difference in germ cell composition in samples with/without YCMD. This problem was tackled by selecting 26 samples from an initial cohort of 34 samples by their homogeneity in respect to cellular composition as obtained from gene expression clustering. This way, the ‘germ cell effect’ was minimized, and a distinct ‘deletion effect’ became more apparent. Several hundred genes are influenced by YCMD as shown on the three different phenotypes hypospermatogenesis, meiotic arrest, and Sertoli-cell only syndrome. We validated on an independent cohort of samples five genes by quantitative real-time PCR that are expressed in germ cells or the somatic compartment and which are exclusively altered by the presence of YCMD. We conclude that the deletion of Y-chromosomal genes has a significant effect on spermatogenesis by modulating the transcriptional network of the germ cell and somatic compartment.
The online version of this article (doi:10.1007/s00439-010-0865-9) contains supplementary material, which is available to authorized users.
Y-chromosomal microdeletions (YCMD) are the major genetic cause for primary spermatogenic failure. 5–10% of all infertile men can be attributed to a YCMD genotype (Ferlin et al. 2006) although there is also a frequency of 1% in men without fertility problems (Kent-First et al. 1999). Three recurrently deleted regions are implicated in YCMD, namely the azoospermia factor regions AZFa, AZFb, and AZFc located in Yq11. Recent sequencing and mapping of the Y-chromosome has led to the understanding that most of the YCMD are the consequence of de novo homologous rearrangements of repeated sequence blocks at the time of recombination (Skaletsky et al. 2003). No single genes within these regions can be linked to spermatogenic failure as each of the deletions is accompanied by the removal of several genes (Krausz et al. 2006), and deletion patterns are so variable in detail that they can hardly be related to a specific phenotype. Although reports exist of an azoospermic phenotype correlating to mutations in the USP9Y gene of the AZFa locus (Sun et al. 1999), this is still in question (Luddi et al. 2009). Linking YCMD to a specific spermatogenic phenotype is further complicated by the observation of partial YCMD (Hucklenbroich et al. 2005) and the development of at least 153 haplogroups during human evolution, so that the functional contribution of the AZF region to spermatogenesis may vary significantly in men with different haplotypes (Vogt 2005).
A minimum of 14 protein-encoding Y gene families localize to AZF regions (Vogt et al. 2008), namely two for AZFa, seven for AZFb, four for AZFc, and one for AZFb + c. The most common deletions pertain to the AZFc (60%), AZFb (16%), and AZFb + c (14%) regions (Cram et al. 2006), but there is evidence for a potential bias between different populations (Simoni et al. 2008). Several studies indicate that the expression profile of genes residing within the AZFc/AZFb regions is specific to and essential for individual developmental stages of human spermatogenesis (Kleiman et al. 2007; Lardone et al. 2007; Shinka et al. 2004; Sato et al. 2006; Stouffs et al. 2004). In this report, we interrogated human testicular biopsies from men with hypospermatogenesis (HYS), meiotic arrest (MA), and Sertoli-cell-only syndrome (SCO) by global gene expression profiling, comparing phenotypes presenting no or AZFc/AZFb + c deletions. This approach has the advantage that not only transcripts of the Y-linked genes are under investigation, but one may discover alterations in gene expression that are up-/downstream of the AZF genes and therefore regulatory networks in which AZF genes play a putative role. Unveiling possible effects of AZF deletions on the testicular transcriptome is feasible since only expression data in respect to single AZF genes exist to date, and gene expression profiling of human spermatogenesis has already given vital insight to stage-specific differences in cases of spermatogenic impairment (Lian et al. 2009; Okada et al. 2008; Ellis et al. 2007; Spiess et al. 2007). A major obstacle that had to be overcome was to identify transcriptional regulation resulting from the ‘deletion effect’ and not from the dominating effect of altered germ cell composition as a consequence of YCMD (Ferras et al. 2004), herein termed as ‘germ cell effect’. This problem was minimized by (1) reducing an initial cohort of 34 samples with and without YCMD to a subset of 26 samples that featured similar germ cell composition based on their selection by gene expression profiling, (2) enrichment of differentially expressed genes for the somatic compartment by an intersection approach, and finally, (3) validation of the absence of any residual ‘germ cell effect’ in these samples by querying independent datasets of gene expression in human spermatogenesis.
Testicular biopsies were obtained from patients presenting at the Department of Andrology, University Hospital Hamburg-Eppendorf, Germany, between August 2004 and August 2006. Informed consent and Ethic Committee Approval was obtained (OB/X/2000), and the study was conducted in accordance with the guidelines of the ‘Helsinki Declaration’. Tissues used for the microarray analysis were taken independently for therapeutic testicular sperm extraction (TESE), RNA extraction, and diagnostic purposes as described previously (Jezek et al. 1998). In total, we selected seven samples with hypospermatogenesis (4 controls and 3 AZFc), 14 samples with meiotic arrest (10 controls, 3 AZFc, and 1 AZFbc) and 13 samples with Sertoli-cell-only syndrome (7 controls, 5 AZFc, and 1 AZFbc) from a cohort of 700 biopsies. Selection criteria were based on homogeneous histology within the same group, in which the histology was extrapolated from examinations of parallel biopsies from the same testis as part of the routine diagnostic work-up. As a final selection step, results from a trial-TESE were taken into account to exclude any discrepancies between the morphological classification of the first biopsy and the spermatogenic activity in a second biopsy from another area of the same testis (compare Feig et al. 2007).
All samples used in this work were screened for AZFa, AZFb, or AZFc deletions according to the ‘Best practice guidelines for molecular diagnosis of y-chromosomal microdeletions’ (Simoni et al. 2004). Amplicons of the multiplex PCR were analyzed by electrophoresis on MetaPhor™ agarose gels (3%, Cambrex Corporation) in TAE buffer by running for 16 h at 25 V and subsequent ethidium bromide staining.
Briefly, at least four fragments per testis were processed for histological analysis, post-surgical trial-TESE, and cryopreservation. An additional sample of rice grain size (approximately 30 mg) was submerged in 2 ml RNAlater® (Ambion, Austin, TX, USA) for RNA extraction and subsequent gene expression profiling. Total RNA was extracted in RNApure (Peqlab, Erlangen, Germany) and re-purified on RNeasy columns (Qiagen, Hilden, Germany) according to the manufacturers’ protocols. Purity and integrity (28S/18S ratio) were assessed by loading 200 ng aliquots onto RNA 6000 nano assay chips using an Agilent Bioanalyzer 2100 (Agilent Technologies, Palo Alto, CA, USA). Only samples with an RNA integrity number higher than 7.5 (RIN, Agilent software) were included in the analyses. Amplification and labeling of cRNA was done according to the manufacturer’s protocol (GE Healthcare, Piscataway, NJ, USA). 10 μg of cRNA were hybridized to a Codelink™ Human 20 K Bioarray (GE Healthcare) containing 19,881 gene-specific probes. Arrays were stringency washed, stained with Cy5™-Streptavidin (GE Healthcare) and washed again according to the manufacturer’s instructions. Slides were dried by centrifugation and immediately scanned on a 428™ Array Scanner (Affymetrix) using Jaguar 2.0 Software. Images were analyzed and quantified using the CodeLink™ Expression Analysis Software v4.1 (GE Healthcare) and rigorously quality controlled for hybridization artefacts.
Genes were filtered that had significant expression values across all samples and average expression higher than the 5th percentile of all experiments. The experiments were normalized by quantile normalization (R package affy, http://www.bioconductor.org). In compliance with the MIAME guidelines (Brazma et al. 2001), raw and processed data files of all samples were deposited in the Gene Expression Omnibus (GSE21613).
The top 200 genes based on their variance across samples or genes that were obtained by statistical filtering were subjected to hierarchical clustering analysis (Manhattan distance, average metric). Additionally, the complete dataset without prior filtering was interrogated by principle component analysis and three-dimensional display of the first three components to visualize the main structure in the data. All clustering was conducted with the freely available TIGR MeV v4.5 software (http://www.tm4.org/mev.html). After biopsy sample selection, the top 200 genes based on variance were analyzed in respect to their cluster stability by bootstrap hierarchical clustering (Manhattan distance, average metric).
To select genes whose expressions are altered by AZFb/AZFb + c deletions, genes were filtered by t-tests between HYS/HYS-AZFc, MA/MA-AZFc, and MA/MA-AZFb + c. We used an equal-variance t-test as homogeneity of variance for the different groups, which was found for 91, 90, and 93% of all 19,227 genes from the complete data set, respectively. The obtained p values were subjected to a Benjamini–Hochberg multiple testing correction (Benjamini and Hochberg 1995). An additional t-test was conducted between HYS and SCO subtypes (without deletions) to filter a set of genes that define transcriptional changes throughout the complete human spermatogenesis (‘germ cell effect’), and an intersection between this list and those obtained from the AZFc/AZFb + c deletion analysis was made. Furthermore, to enrich for genes exhibiting a ‘deletion effect’, this overlap was removed from the initial AZFc/AZFb + c lists. Finally, to validate the enrichment of AZFc/AZFb + c affected genes with reduced ‘germ cell effect’, these genes were matched (either by identical probe ID’s or HUGO gene abbreviations) with expression data from two different independent microarray analyses of spermatogenic stages using a completely different cohort of samples (GSE4797, NCBI Gene Expression Omnibus; Feig et al. 2007; E-TABM-234, EMBL ArrayExpress; Spiess et al. 2007).
The complete analysis above was done within the freely available R statistical environment (http://www.r-project.org). Over-representation analysis of functional groups using gene ontology (GO) terms was done using the HUGO gene abbreviations and the DAVID 2008 server (http://david.abcc.ncifcrf.gov/).
Reverse transcription from total RNA was carried out with the SuperScript® III First-Strand Synthesis SuperMix for qRT-PCR as described by the manufacturer (Invitrogen).
Quantitative RT-PCR (qRT-PCR) was performed on a LightCycler™ (Roche, Basel, Switzerland) using the SYBR®Premix ExTaq™ (Takara Bio Inc.) with working procedures according to the manufacturer’s instructions, but with primer concentrations of 0.4 μM in a total volume of 20 μl. Primer sequences and amplicon sizes are available in Supplemental Data 4. Cycling conditions were: 30 s 95°C initial denaturation, 30–40 cycles of 10 s at 95°C, 20 s at 60–64°C annealing temperature, and 30 s at 72°C. Transcript levels were normalized to ribosomal protein S27 transcript, whose expression was found to be highly stable over the different samples by microarray expression values. Fold differences were calculated by use of the five-parameter sigmoidal model (Ritz and Spiess 2008). Products were analyzed in respect to correct amplicon size on an Agilent Bioanalyzer (Model 2100; Agilent Technologies, Palo Alto, CA, USA). We validated the most interesting candidates by quantitative real-time PCR on a completely independent test set of 13 samples.
A major issue to resolve in this work was the discovery and selection of genes that are affected by AZFc/AZFb + c deletions in a background of transcriptional changes that are dominated by the depletion of specific germ cell types in the different spermatogenic stages when the deletions are present. To eliminate the dominating effect of germ cell-specific transcript levels within the comparison of two (or more) biopsies, it is an indispensable prerequisite to start with samples that present a homogeneous germ cell phenotype. Figure 1 exemplifies this problem using samples from our initial set. When comparing biopsies presenting continuously hypospermatogenesis with those that have an AZFc deletion and a ‘mixed atrophy’ phenotype with predominantly Sertoli-cell-only tubules (upper panel), any statistical filtering method will select genes specific to all germ cells that differ between the two states, i.e., spermatogonia, primary/secondary spermatocytes, and round/elongated spermatids. Furthermore, any clustering method will result in groups defining the subtypes of spermatogenic failure. In contrast, if one employs samples with similar germ cell composition (lower panel), the dominating ‘germ cell effect’ is minimized, and a ‘deletion effect’ might be more apparent.
Our initial set of 34 biopsies (Supplemental Data 1, ‘initial set’) had been classified on the basis of germ cell types that are present in the seminiferous tubules. We grouped the top 200 variant genes by hierarchical clustering (Fig. 2, upper panel). Additional information was acquired about the relative quantity of germ cell composition found in the histological examination of the biopsies (lower panel), illustrating the paradigm shift from the most differentiated germ cell type to that which is most frequent: the separation of groups in the dendrogram correlates not so much with initial histological classification but significantly more with the dominant germ cell type present. For instance, sample SCO_10_c clusters within the MA samples, which tallies with the observation of tubules containing spermatocytes (MA) being the most frequent observation in this sample. A comparison of the dendrogram with the histology for all samples clearly shows (and without exception) that all samples cluster according to their quantitatively prevailing germ cell type. This is an important observation as many of the initial samples with YCMD cluster in groups distinct from their most differentiated germ cell type. To complement these results, we employed principle component analysis (PCA) on the 34 samples using all genes. As PCA is a data reduction technique, one can visualize the first three main variances (principle components) within the gene expression data by three-dimensional display (Supplemental Data 2). Corroborating the dominant effect of relative germ cell composition on the grouping outcome, the main variance of gene expression within the 34 samples is again defined by the most frequent germ cell type (spread along the first principle component PCA1). A similar observation had been made by Ellis et al. (2007) using a mixture of arrest types with variable etiology. The second and third principle components (PCA2 and PCA3) separate a mixture of germ cell composition and YCMD status. It must be noted, however, that PCA1-3 only explained 65% of the total variance found in the data, and therefore other, yet undescribed effects seem likely. We subjected the top 200 variant genes to GO (gene ontology) term analysis and found a highly significant over-representation of terms such as spermatogenesis (p = 2E−14), mitosis (p = 5E−6) and meiosis (p = 1E−3). These results strengthen once again the fact that germ cell transcription is the dominating effect found in our initial sample set. Furthermore, the findings above indicate that the clustering of gene expression profiles from testicular biopsies is extremely sensitive to the relative germ cell composition, a property that enables us to select homogeneous samples from the clustering outcome.
Based on the grouping of samples in Figs. 2 and and33 we selected 26 samples for further analysis (Supplemental Data 1, ‘selected set’). We clustered the top 200 variant genes from the selected subset of 26 samples by bootstrap hierarchical clustering to acquire some additional information about cluster stabilities (Fig. 3). Similar to the previous grouping structure of the initial set of 34 samples, we obtained three main clusters defining the different pathological conditions HYS, MA, and SCO as the most variant genes correlate also here with the germ cell state. In contrast to the initial clustering, samples with YCMD now grouped with the corresponding germ cell phenotypes without deletion and this with relatively high cluster stability. The procedure therefore resulted in a subset in which the difference in germ cell composition between samples of the same arrest type but with absence/presence of YCMD was minimized.
We conducted t-tests with subsequent Benjamini–Hochberg multiple testing correction on the MA and SCO phenotypes of this subset between normal/AZFc samples. For the HYS group, we were only able to calculate a ratio between the averaged normal samples and one sample with AZFc due to the fact that possible replicates had been sorted out by the sample selection procedure. We filtered 214 genes with fourfold up/down-regulation for the HYS/HYS-AZFc group, 495 genes with pcorr < 0.05 for the MA/MA-AZFc group, 22 genes with pcorr < 0.05 for the SCO/SCO-AZFc group, and 18 genes with praw < 0.05 for the SCO/SCO-AZFb + c group (Supplemental Data 1, ‘HYS AZFc’, ‘MA AZFc’, ‘SCO AZFc’, and ‘SCO AZFbc’). This way, we obtained lists of differential gene expression between normal samples and those with either AZFc or AZFbc deletions specific to each of the three different spermatogenic failure types and putatively enriched for the ‘deletion effect’. In terms of gene numbers, the largest difference could be attributed to the MA phenotype with 495 genes that account for roughly 3% of genes that were under investigation. A further interest was to obtain a set of genes that show differential expression irrespective of the spermatogenic state. Interesting to note here is that for both HYS and MA, DAZ2 was found to be the transcript with the highest ratio of downregulation (~20-fold), which demonstrates that our gene expression data mirrors the depletion of DAZ transcripts as a consequence of the AZFc deletion. However, the probes used on our array platform were not able to discriminate between the four isoforms of DAZ (DAZ1-4). Furthermore, no overlap was observed between all three failure types, suggesting that the majority of affected genes correlate with specific germ cell expression.
Gene ontology analysis of differential genes within HYS and MA types in the presence of an AZFc deletion (SCO failed to be analyzed due to the small gene number) indicates functional differences in dependence of the phenotype (Table 1). Nevertheless, terms such as ‘spermatogenesis’, ‘sexual reproduction’, and ‘meiosis’, although having been decreased by the selection procedure by 2–10 orders of magnitude in p values, could still be attributed to the fact that our selection procedure was not optimal. We added another step of minimizing the ‘germ cell effect’ by eliminating all genes by an intersection approach from the obtained lists that are differentially expressed throughout human spermatogenesis, based on 2,568 genes obtained from t-tests between HYS and SCO normal samples, in analogy to Feig et al. (2007) (Supplemental Data 1, ‘Germ cell genes’). Intersections are given in Fig. 4. We removed these intersections from the initial HYS and MA gene lists, and reanalyzed in respect to GO term over-representation (Table 1, terms in bold). This resulted in the elimination of the term ‘sexual reproduction’ from the HYS group and a significant reduction of genes that correlate with the terms ‘spermatogenesis’ and ‘meiosis’ as shown on the increase of p values at 4–8 orders of magnitude. Several other GO terms such as ‘negative regulation of glucocorticoid synthesis’ and ‘killer cell-mediated cytotoxicity’ appeared after elimination of intersecting genes. An overlap analysis of the reduced sets resulted in seven genes that are up-/downregulated in the presence of AZFc in HYS as well as MA phenotype (Supplemental Data 1, ‘HYS AZFc vs. MA AZFc’).
In a last step, we queried the ‘germ cell-depleted’ gene lists on two gene expression datasets of human spermatogenesis that were obtained from a completely different cohort of samples as the analysis of microarray data has the inherent problem of ‘overtraining’ and thus can lack generality (Zervakis et al. 2009). The remaining 178 genes of the reduced HYS dataset (Supplemental Data 1, ‘HYS wo GSG’) showed no significant trend of expression changes throughout spermatogenesis (Supplemental Data 3, upper panel). The analysis of the reduced MA dataset of 232 genes (Supplemental Data 1, ‘MA wo GSG’) revealed a small increase from SCO to MA that was not statistically significant on the averaged values (Supplemental Data 3, lower panel; Codelink platform: 3.8-fold, p = 0.27; Affymetrix platform: 2.23-fold; p = 0.08), which may be interpreted as a remnant and non-removable ‘germ cell effect’ although other interpretations are possible (see “Discussion”).
A transcriptional deletion map for Y-chromosomal genes residing within the deleted loci was built, demonstrating that the deletions were mirrored almost completely on the transcriptional level, and that our microarray analysis was sensitive enough to reveal the absence of transcripts as a consequence of the deleted genes (Fig. 5a, b). However, the general impression of Y-chromosomal transcript levels was relatively heterogeneous in the different samples. We also noticed a complete absence of any signal for NLGN4Y transcripts in sample SCO_13_bc in contrast to the second sample with AZFb + c deletion, MA_14_bc (Supplemental Data 1, ‘SCO AZFbc’). To confirm absence of this gene, which is located proximal to the classic AZFb boundary defined by Vogt et al. (1996), we mapped a possible deletion in this region by using gene-specific primers (Supplemental Data 4) for NLGN4Y, its proximal genes VCY and TMSB4Y, and two STS markers proximal and distal of the NLGN4Y gene (sY95 and sY98). For this sample, the complete chromosomal region for AZFb + c reaching proximally to sY95 was deleted (Fig. 5c–e), and therefore, NLGN4Y transcripts are undetectable due to the presence of a relatively rare massive AZFb + c deletion (Silber et al. 1998), also referred to as ‘cytogenetic abnormality’ (Krausz et al. 1999).
To validate our microarray data, we conducted quantitative real-time PCR on selected genes using 13 independent samples with/without deletions. All transcripts were normalized to the S27a transcript, which was shown in recent microarray investigations (compare GEO Database Accessions GSE4797, GSE6023, GSE18997) and also this study (GSE21613) to be highly stable expressed in testicular tissue with different states of spermatogenic impairment. It was therefore selected as a good normalization gene to prevent wrong interpretations of expression changes due to differences in the amount of input cDNA in the real-time PCR reactions (Supplemental Data 5A + B). We confirmed a strong downregulation for HIV-1 Rev-binding protein (HRB; 7.5-fold array, 2.8-fold qPCR) in samples with hypospermatogenesis (HYS_7_c). For this sample, there were several other transcripts with more than tenfold upregulation in the presence of an AZFc deletion (marked in red in Supplemental Data 1), but we failed to verify several of these by qPCR (data not shown). Here, the results for genes with a very high upregulation must be taken critically as we cannot identify hybridization artefacts when replicates are lacking. In contrast to this, other genes with interesting functions in respect to AZF deletions that had prominent down-/upregulation in samples with meiotic arrest or Sertoli-cell-only syndrome (meiotic arrest: breast cancer antigen NY-BR-1, apolipoprotein E; Sertoli-cell only: ring finger protein 126; bradykinin receptor B1) could be validated in respect to the direction of expression changes by qPCR, most probably because replicates were available. The ratios obtained from microarray data and qPCR analysis differed slightly, a common finding that may be attributed to the effect of ‘ratio compression’ for microarray data (Wang et al. 2006).
It is necessary to dissect deletion-induced disturbances of transcriptional networks from those changes that are a result of relative transcript levels when the germ cell composition is different between non-deleted and deleted samples. We and others have found that the phenotype of ‘mixed atrophy’, i.e., a high heterogeneity in tubular germ cell composition within the same biopsy, is a common finding in patients with YCMD (Sigg and Hedinger 1981; Luetjens et al. 2002; Ferras et al. 2004). This quantitative difference in cell type-specific transcripts is significantly shielding and obscuring any deletion-induced effects for downstream analysis with sensitive methods such as microarray analysis, leading to gene expression patterns that reflect only cellular composition of the samples (Okada et al. 2008; Ellis et al. 2007; Feig et al. 2007). This classification, which is based exclusively on the dominating cell type, diverges from any classification using conventional histological descriptions. After enrichment of the ‘deletion effect’, we obtained a relatively high number of differential genes when YCMD are present. If using these as an indicator of “deletion effect strength”, we obtained for the different testicular phenotypes the order SCO < HYS < MA. This could be the consequence of decreased SCO differential expression pertaining mostly to the somatic compartment of the testis as germ cells are lacking, and meiosis being most strongly perturbed as haploid gamete formation is severely modulated by DAZ (Kee et al. 2009; Geoffroy-Siraudin et al. 2007).
The magnitude of transcriptional changes was the highest for those genes that reside within the respective deleted region, i.e., DAZ for AZFc deletions and some of the AZFb region located genes (SMCY, EIF1AY, CYorf14, NLGN4Y). This anticipated result tallies with the transcriptional changes observed in YCMD samples based on PCR methodology (Kleiman et al. 2007; Lardone et al. 2007). Nevertheless, several hundred other genes were down-/upregulated in the presence of YCMD, although the magnitude of transcriptional changes was mainly in a region of 2–5-fold, indicating a transcriptional system with a possible high redundancy that would result in the high heterogeneity of spermatogenic arrests in the presence of YCMD as frequently observed (Foresta et al. 2001). The precise mechanisms for the relative extensive amount of genes with altered transcription in the presence of YCMD remain to be resolved. However, the disturbance of the testicular transcriptional network by YCMD is most likely a consequence of the lacking RNA-/DNA-binding functionality of the deleted genes which is important for any transcriptional network. For instance, DAZ-family members have been shown to bind to specific transcripts in vivo and regulate pre-meiotic translation of important factors such as the mouse VASA-homolog MVH (Reynolds et al. 2005), which in turn plays a major role in the processing of PIWI-interacting RNAs (Kuramochi-Miyagawa et al. 2010). RNA binding has also been observed for the splicing-associated RBMY (Elliott et al. 2000), while other genes residing within the AZFb locus exhibit transcription factor activity (HSFY; Tessari et al. 2004) or are involved in ribosome assembly (VCY; Zou et al. 2003). Hence, impairment of the testicular transcriptional network by altering the function of transcription factors through modulation of their mRNAs or their translation could play a putative role.
From our set of real-time PCR-validated genes, HRB (HIV-1 Rev-binding protein) could be confirmed to be downregulated approximately threefold in the presence of an AZFc deletion. This finding is interesting in respect to HRB’s vital function for acrosome formation (Kang-Decker et al. 2001). Due to the lack of replicates in this group, these results should be viewed critically. For the phenotype of meiotic arrest in the presence of AZFc deletions, we confirmed ANKRD30A (breast cancer antigen NY-BR-1) to be downregulated approximately threefold, while APOE (apolipoprotein E) was upregulated approximately twofold. While the former is expressed almost exclusively in primary spermatocytes (Theurillat et al. 2008), the latter is largely expressed in the Leydig cells (Law et al. 1997), suggesting modulations of the interstitial compartment in the presence of AZFc deletions. For ring finger protein 126 (RNF126), which was validated to be approximately twofold upregulated in SCO samples with AZFc, no function in testis is known to date. In case of AZFb + c deletions with SCO phenotype, we validated BDKRB1 (bradykinin receptor B1) as one of the few transcripts with increased transcript levels (approximately fourfold). As BDKRB1, at least in the rat model, elevates intracellular calcium levels in peritubular cells after bradykinin stimulation (Wennemuth et al. 2003), involvement of the lamina propria is also likely.
In this study, we did not include other factors such as the patients’ haplotypes, first, because we wanted to identify transcriptional changes independent of haplotype background, and second, because there are indications that YCMD formation is more or less a stochastic event (Quintana-Murci et al. 2001). Nevertheless, future work should investigate the effect of haplotype background on global testicular expression profiles. Another drawback of our procedure is that it might mask interesting germ cell genes that are influenced by YCMD, and that a significantly higher number of microarray experiments might be necessary to reveal further or more subtle effects of YCMD-induced expression changes.
The obtained list of several hundred differential genes may serve as a starting point for a more detailed analysis of YCMD-induced testicular changes, either in germ cells or the somatic compartment. Moreover, the identified transcriptional networks that are disturbed might be pharmaceutically targeted to reconstitute impaired spermatogenesis in men with YCMD offering potential for male infertility treatment.
Note added During the review process of this paper, we became aware of an investigation with a similar topic that was published by Gatta et al. (2010). In this study, transcriptional changes due to Y-chromosomal microdeletions were investigated by microarray analysis of testicular biopsies without pre-selection based on clustering to minimize the quantitative ‘germ cell effect’. The results in respect to differential genes from this and our study differ substantially.
Below is the link to the electronic supplementary material.
Financial support for this work was in part from grant Sp721/1-3 of the German Research Foundation to ANS and HCO. We would like to thank Prof. Christiane Kirchhoff for fruitful discussions and Dr. Caroline Feig for technical assistance.
Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.