The gonadotropin treatment of infertile men may improve spermatogenesis and lead to sperm cell production, however, only a small fraction of treated patients positively responds to such therapy. To identify individual treatment prognostic biomarkers associated with responsiveness to gonadotropins, we compared the gene expression profiles of testicular oligobiopsies from 3 patients with non-obstructive azoospermia (NOA) who positively responded to therapy with a combination of human chorionic gonadotropin and recombinant follicle-stimulating hormone (hCG/rFSH) to those of 3 non-responders. We used Affymetrix Human Gene 1.0 ST microarrays. The results of the microarray evaluation were validated by the qPCR technique while gene variants of the HLA-DQB1 (major histocompatibility complex, class II, DQ beta 1) were subsequently sequenced. In our microarrays, we have identified most significantly 5 transcripts with different expression levels in responders versus non-responders groups. Our interest has been primarily focused on the transcript associated with the HLA-DQB1 gene. Because the expression of this gene was up-regulated in the non-responding patients and only patients with heterozygotic alleles of HLA-DQB1 turned out to be positive to gonadotropin therapy, we suggest that this gene may be a biomarker of potential significance for the gonadotropin treatment of male infertility. We also compared the testicular gene expression profile in one individual before and after gonadotropin treatment. In the re-biopsied sample, we have identified over 600 genes that showed differences in testicular expression; some of these genes are critical for spermiogenesis. Thus, we documented that the applied gonadotropins successfully stimulated the spermatogenetic wave in patients with NOA.