Alternative 3′ UTR polyadenylation of Bzw1 transcripts display differential translation efficiency and tissue-specific expression

BZW1 is a prognostic and immunological biomarker in pancreatic adenocarcinoma

Pancreatic adenocarcinoma is the most common malignant tumor of the digestive system and is called the "king of cancer" because it has been labeled with high malignancy, rapid progression, poor survival, and poor prognosis. Previously, it was reported that the basic leucine zipper and W2 domains 1 (BZW1) is involved in the progression of many tumors. However, its research in digestive system tumors such as pancreatic cancer is rarely studied. To explore potential biomarkers related to survival and prognosis of pancreatic cancer and provide a new targeted therapy for it. We first analyzed the mRNA and protein expression of BZW1 in pancreatic cancer. We then explored the correlation of BZW1 with survival prognosis and immune infiltration in pancreatic cancer patients. Finally, we explored BZW1-related gene enrichment analysis, including protein-protein interaction networks, gene ontology functional enrichment analysis, and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis. The mRNA and protein expression of the BZW1 gene in pancreatic cancer tissues were higher than those in adjacent normal tissues, and pancreatic cancer patients with high BZW1 expression had a poor prognosis. In addition, the expression of BZW1 was positively or negatively correlated with different immune cells of pancreatic cancer, such as CD4 + T lymphocytes, CD8 + T lymphocytes, B cells, macrophages, neutrophils, etc. Correlation enrichment analysis showed that we obtained 50 available experimentally determined BZW1-binding proteins and 100 targeted genes related to BZW1, and the intersection genes were eukaryotic translation termination factor 1 and Guanine nucleotide binding protein, alpha inhibiting activity polypeptide 3. Moreover, there was a positive correlation between BZW1 and eukaryotic translation termination factor 1 and Guanine nucleotide binding protein, alpha inhibiting activity polypeptide 3 genes in pancreatic cancer. Gene ontology enrichment analysis showed BZW1 was mainly related to biological processes such as "mRNA processing," "RNA splicing," "regulation of translational initiation," and "activation of innate immune response." The results of Kyoto Encyclopedia of Genes and Genomes pathway analysis further indicated that BZW1 may be involved in pancreatic carcinogenesis through the "spliceosome" and "ribosome." The BZW1 gene may be a potential immunotherapy target and a promising prognostic marker for pancreatic cancer.

Activation of BZW1 by CEBPB in macrophages promotes eIF2α phosphorylation-mediated metabolic reprogramming and endoplasmic reticulum stress in MRL/lpr lupus-prone mice

BACKGROUND

Lupus nephritis (LN) is associated with significant mortality and morbidity, while effective therapeutics and biomarkers are limited since the pathogenesis is complex. This study investigated the roles of the CEBPB/BZW1/eIF2α axis in metabolic reprogramming and endoplasmic reticulum stress in LN.

METHOD

The differentially expressed genes in LN were screened using bioinformatics tools. The expression of CEBPB in the renal tissue of patients with LN and its correlation with the levels of creatinine and urinary protein were analyzed. We used adenoviral vectors to construct LN mice with knockdown CEBPB using MRL/lpr lupus-prone mice and analyzed the physiological and autoimmune indices in mice. Chromatin immunoprecipitation quantitative polymerase chain reaction (ChIP-qPCR) and dual-luciferase reporter assays were conducted to explore the regulation of BZW1 by CEBPB, followed by glycolytic flux analysis, glucose uptake, and enzyme-linked immunosorbent assay (ELISA). Finally, the role of eIF2α phosphorylation by BZW1 in bone marrow-derived macrophages (BMDM) was explored using eIF2α phosphorylation and endoplasmic reticulum stress inhibitors.

RESULTS

CEBPB was significantly increased in renal tissues of patients with LN and positively correlated with creatinine and urine protein levels in patients. Downregulation of CEBPB alleviated the autoimmune response and the development of nephritis in LN mice. Transcriptional activation of BZW1 by CEBPB-mediated glucose metabolic reprogramming in macrophages, and upregulation of BZW1 reversed the mitigating effect of CEBPB knockdown on LN. Regulation of eIF2α phosphorylation levels by BZW1 promoted endoplasmic reticulum stress-amplified inflammatory responses in BMDM.

CONCLUSION

Transcriptional activation of BZW1 by CEBPB promoted phosphorylation of eIF2α to promote macrophage glycolysis and endoplasmic reticulum stress in the development of LN.

Expression, oncological and immunological characterizations of BZW1/2 in pancreatic adenocarcinoma

Background: Despite the progress in early diagnosis and treatment, prognosis of pancreatic adenocarcinoma (PAAD) is still poor. Basic leucine zipper and W2 domain-containing protein 1 (BZW1) and protein 2 (BZW2) are attached to the basic leucine zipper (bZIP) superfamily. Recently, BZW1 was identified as an important role in glycolysis of PAAD. However, the comprehensive reports about BZW1/2 in PAAD are not sufficient.

Methods: RNA-seq data in the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases were retrospectively analyzed. We explored the expression of BZW1/2 in PAAD tissues and the associations between BZW1/2 and prognosis. In addition, the potential roles of BZW1/2 in tumor microenvironment (TME) of PAAD were analyzed. Finally, clinicopathological data of 49 patients with PAAD in our institution were collected. Immunohistochemistry was used to determine the expression of BZW1/2 in PAAD samples.

Results: BZW1 and BZW2 were upregulated in PAAD tissues compared to normal tissues (p < 0.05). The expression of BZW1/2 were not significantly correlated with gender, grade and stage of PAAD (p > 0.05). High expression of BZW2 was an independent predictor for poor prognosis of PAAD (HR 1.834, 95%CI 1.303–2.581, p = 0.001). And a nomogram to predict overall survival (OS) of PAAD was established with a C-index of 0.685. BZW1 and BZW2 expression were positively associated with T cell mediated immune response to tumor cell and Th2 cells in xCell database. Tumor Immune Single-Cell Hub (TISCH) analyses indicated that BZW1 and BZW2 were mainly expressed in B cells and malignant cells. External cohort furtherly validated that high expression of BZW1 and BZW2 were predictors for poor prognosis of PAAD.

Conclusion: We found that BZW1 and BZW2 are highly expressed in malignant cells and B cells in the TME of PAAD. BZW2 is an independent predictor for OS of PAAD. BZW1 and BZW2 expression are positively associated with T cell mediated immune response to tumor cell and Th2 cells in PAAD.

BZW1 Facilitates Glycolysis and Promotes Tumor Growth in Pancreatic Ductal Adenocarcinoma Through Potentiating eIF2α Phosphorylation

BACKGROUND & AIMS

Pancreatic ductal adenocarcinoma (PDAC) is characterized by severe metabolic stress due to fibrosis and poor vascularization. BZW1 is an eIF5-mimic protein involved in tumorigenesis and progression. The aim of this study was to investigate the role of BZW1 in metabolic stress resistance in PDAC.

METHODS

BZW1 expression was evaluated in human PDAC tissue microarray and PDAC cells. Glycolysis regulation of BZW1 and its correlation with glycolysis-related genes was analyzed. Tumor growth, cell proliferation, and apoptosis were evaluated in mice xenograft tumors and patient-derived organoids.

RESULTS

The results of bioinformatic screening identified that BZW1 was 1 of the top 3 genes favorable for tumor progression in PDAC. The analysis of our cohort confirmed that BZW1 was overexpressed in human PDAC tissues compared with nontumor tissues, and its abnormal expression was correlated with large tumor size and poor prognosis. BZW1 promoted cell proliferation and inhibited apoptosis in both mouse xenograft models and PDAC-derived organoids via facilitating glycolysis in the oxygen-glucose-deprivation condition. Mechanically, BZW1 served as an adaptor for PKR-like endoplasmic reticulum (ER) kinase (PERK), facilitated the phosphorylation of eIF2α, promoted internal ribosome entry site-dependent translation of HIF1α and c-Myc, and thereby boosted the Warburg effect. In organoid-based xenografts with high BZW1 levels, both the PERK/eIF2α phosphorylation inhibitor GSK2606414 and ISRIB significantly suppressed tumor growth and prolonged animal survival.

CONCLUSIONS

BZW1 is a key molecule in the internal ribosome entry site-dependent translation of HIF1α/c-Myc and plays crucial roles in the glycolysis of PDAC. BZW1 might serve as a therapeutic target for patients with pancreatic cancer.

Dynamic Variations of 3'UTR Length Reprogram the mRNA Regulatory Landscape

This paper concerns 3′-untranslated regions (3′UTRs) of mRNAs, which are non-coding regulatory platforms that control stability, fate and the correct spatiotemporal translation of mRNAs. Many mRNAs have polymorphic 3′UTR regions. Controlling 3′UTR length and sequence facilitates the regulation of the accessibility of functional effectors (RNA binding proteins, miRNAs or other ncRNAs) to 3′UTR functional boxes and motifs and the establishment of different regulatory landscapes for mRNA function. In this context, shortening of 3′UTRs would loosen miRNA or protein-based mechanisms of mRNA degradation, while 3′UTR lengthening would strengthen accessibility to these effectors. Alterations in the mechanisms regulating 3′UTR length would result in widespread deregulation of gene expression that could eventually lead to diseases likely linked to the loss (or acquisition) of specific miRNA binding sites. Here, we will review the mechanisms that control 3′UTR length dynamics and their alterations in human disorders. We will discuss, from a mechanistic point of view centered on the molecular machineries involved, the generation of 3′UTR variability by the use of alternative polyadenylation and cleavage sites, of mutually exclusive terminal alternative exons (exon skipping) as well as by the process of exonization of Alu cassettes to generate new 3′UTRs with differential functional features.

Seeking a Role for Translational Control by Alternative Polyadenylation in Saccharomyces cerevisiae

Alternative polyadenylation (APA) represents an important mechanism for regulating isoform-specific translation efficiency, stability, and localisation. Though some progress has been made in understanding its consequences in metazoans, the role of APA in the model organism Saccharomyces cerevisiae remains a relative mystery because, despite abundant studies on the translational state of mRNA, none differentiate mRNA isoforms’ alternative 3′-end. This review discusses the implications of alternative polyadenylation in S. cerevisiae using other organisms to draw inferences. Given the foundational role that research in this yeast has played in the discovery of the mechanisms of cleavage and polyadenylation and in the drivers of APA, it is surprising that such an inference is required. However, because advances in ribosome profiling are insensitive to APA, how it impacts translation is still unclear. To bridge the gap between widespread observed APA and the discovery of any functional consequence, we also provide a review of the experimental techniques used to uncover the functional importance of 3′ UTR isoforms on translation.

Alternative Polyadenylation of Mammalian Transcripts Is Generally Deleterious, Not Adaptive

Alternative polyadenylation (APA) produces from the same gene multiple mature RNAs with varying 3' ends. Although APA is commonly believed to generate beneficial functional diversity and be adaptive, we hypothesize that most genes have one optimal polyadenylation site and that APA is caused largely by deleterious polyadenylation errors. The error hypothesis, but not the adaptive hypothesis, predicts that, as the expression level of a gene increases, its polyadenylation diversity declines, relative use of the major (presumably optimal) polyadenylation site increases, and that of each minor (presumably nonoptimal) site decreases. It further predicts that the number of polyadenylation signals per gene is smaller than the random expectation and that polyadenylation signals for major but not minor sites are under purifying selection. All of these predictions are confirmed in mammals, suggesting that numerous defective RNAs are produced in normal cells, many phenotypic variations at the molecular level are nonadaptive, and cellular life is noisier than is appreciated.

RAX2: a genome-wide detection method of condition-associated transcription variation

Most mammalian genes have mRNA variants due to alternative promoter usage, alternative splicing, and alternative cleavage and polyadenylation. Expression of alternative RNA isoforms has been found to be associated with tumorigenesis, proliferation and differentiation. Detection of condition-associated transcription variation requires association methods. Traditional association methods such as Pearson chi-square test and Fisher Exact test are single test methods and do not work on count data with replicates. Although the Cochran Mantel Haenszel (CMH) approach can handle replicated count data, our simulations showed that multiple CMH tests still had very low power. To identify condition-associated variation of transcription, we here proposed a ranking analysis of chi-squares (RAX2) for large-scale association analysis. RAX2 is a nonparametric method and has accurate and conservative estimation of FDR profile. Simulations demonstrated that RAX2 performs well in finding condition-associated transcription variants. We applied RAX2 to primary T-cell transcriptomic data and identified 1610 (16.3%) tags associated in transcription with immune stimulation at FDR < 0.05. Most of these tags also had differential expression. Analysis of two and three tags within genes revealed that under immune stimulation short RNA isoforms were preferably used.

DDX3X, the X homologue of AZFa gene DDX3Y, expresses a complex pattern of transcript variants only in the male germ line

DDX3X, the functional X homologue of the major AZFa gene, DDX3Y, belongs to the highly conserved PL10-subfamily of DEAD-box RNA helicase genes which are functionally conserved from yeast to man. They are mainly involved in cell cycle control and translation initiation control of gene transcripts with long 5'UTR extensions containing complex secondary structures. Interestingly, in humans both gene copies were found to be expressed at different phases of human spermatogenesis. Whereas DDX3Y transcripts are translated only in premeiotic male germ cells, the DDX3X protein is expressed only in postmeiotic spermatids. In this study, we found that the major class of DDX3X transcripts in human testis become activated first after meiosis and at a specific core promoter not active in somatic tissues and not present upstream of the DDX3Y homologue. Two alternative 5'UTR transcript lengths are subsequently produced by an additional testis-specific 5'UTR splicing event. Both transcripts are mainly processed for polyadenylation in their proximal 3'UTR. A minor transcript class starting at the same male germ line-specific core promoter produces primary transcripts with an extremely long 3'UTR (∼ 17 kb), which is subsequently spliced at distinct sites resulting in six short 3'UTR splice variants (I-VI). Comparative analyses of the DDX3X transcripts in mouse and primates revealed that this complex pattern of male germ line-specific transcript variants first evolved in primates. Our data thus suggest complex translational control mechanism(s) for the human DDX3X gene locus functioning only in the male germ line and resulting in expression of its protein only in the postmeiotic spermatids.

A dynamic interplay between alternative polyadenylation and microRNA regulation: implications for cancer (Review)

Alternative polyadenylation and microRNA regulation are both mechanisms of post-transcriptional regulation of gene expression. Alternative polyadenylation often results in mRNA isoforms with the same coding sequence but different lengths of 3' UTRs, while microRNAs regulate gene expression by binding to specific mRNA 3' UTRs. In this sense, different isoforms of an mRNA may be differentially regulated by microRNAs, sometimes resulting in cellular proliferation and this mechanism is being speculated on as a potential cause for cancer development.

Maternal inheritance of the Gnas cluster mutation Ex1A-T affects size, implicating NESP55 in growth

Genes subjected to genomic imprinting are often associated with prenatal and postnatal growth. Furthermore, it has been observed that maternally silenced/paternally expressed genes tend to favour offspring growth, whilst paternally silenced/maternally expressed genes will restrict growth. One imprinted cluster in which this has been shown to hold true is the Gnas cluster; of the three proteins expressed from this cluster, two, Gsα and XLαs, have been found to affect postnatal growth in a number of different mouse models. The remaining protein in this cluster, NESP55, has not yet been shown to be involved in growth. We previously described a new mutation, Ex1A-T, which upon paternal transmission resulted in postnatal growth retardation due to loss of imprinting of Gsα and loss of expression of the paternally expressed XLαs. Here we describe maternal inheritance of Ex1A-T which gives rise to a small but highly significant overgrowth phenotype which we attribute to reduction of maternally expressed NESP55.

Tissue-specific alternative polyadenylation at the imprinted gene Mest regulates allelic usage at Copg2

The gene Mest (also known as Peg1) is regulated by genomic imprinting in the mouse and only the paternal allele is active for transcription. MEST is similarly imprinted in humans, where it is a candidate for the growth retardation Silver-Russell syndrome. The MEST protein belongs to an ancient family of hydrolases but its function is still unknown. It is highly conserved in vertebrates although imprinted expression is only observed in marsupials and eutherians, thus a recent evolutionary event. Here we describe the identification of new imprinted RNA products at the Mest locus, longer variants of the RNA, called MestXL, transcribed >10 kb into the downstream antisense gene Copg2. During development MestXL is produced exclusively in the developing central nervous system (CNS) by alternative polyadenylation. Copg2 is biallelically expressed in the embryo except in MestXL-expressing tissues, where we observed preferential expression from the maternal allele. To analyze the function of the MestXL transcripts in Copg2 regulation, we studied the effects of a targeted allele at Mest introducing a truncation in the mRNA. We show that both the formation of the MestXL isoforms and the allelic bias at Copg2 are lost in the CNS of mutants embryos. Our results propose a new mechanism to regulate allelic usage in the mammalian genome, via tissue-specific alternative polyadenylation and transcriptional interference in sense-antisense pairs at imprinted loci.

Mechanisms and consequences of alternative polyadenylation

Alternative polyadenylation (APA) is emerging as a widespread mechanism used to control gene expression. Like alternative splicing, usage of alternative poly(A) sites allows a single gene to encode multiple mRNA transcripts. In some cases, this changes the mRNA coding potential; in other cases, the code remains unchanged but the 3' UTR length is altered, influencing the fate of mRNAs in several ways, for example, by altering the availability of RNA binding protein sites and microRNA binding sites. The mechanisms governing both global and gene-specific APA are only starting to be deciphered. Here we review what is known about these mechanisms and the functional consequences of alternative polyadenylation.

Molecular characterization of two isoforms of ZFAND3 cDNA from the Japanese quail and the leopard gecko, and different expression patterns between testis and ovary

Zing finger AN1-type domain 3 (ZFAND3), also known as testis expressed sequence 27 (Tex27), is a gene found in the mouse testis, but its physiological function is unknown. We identified the full-length sequences of two isoforms (short and long) of ZFAND3 cDNA from Japanese quail and leopard gecko. This is the first cloning of avian and reptilian ZFAND3 cDNA. The two isoforms are generated by alternative polyadenylation in the 3'UTR and have the same ORF sequences encoding identical proteins. There were highly conserved regions in the 3'UTR of the long form near the polyadenylation sites from mammals to amphibians, suggesting that the features for determining the stability of mRNA or translation efficiency differ between isoforms. The deduced amino acid sequence of ZFAND3 has two putative zinc finger domains, an A20-like zinc finger domain at the N-terminal and an AN1-like zinc finger domain at the C-terminal. Sequence analysis revealed an additional exon in the genomic structures of the avian and reptilian ZFAND3 genes which is not present in mammals, amphibians, or fish, and this exon produces additional amino acid residues in the A20-like zinc finger domain. Expression analysis in Japanese quail revealed that the expression level of ZFAND3 mRNA was high in not only the testis but also the ovary, and ZFAND3 mRNA was expressed in both spermatides of the testis and oocytes of the ovary. While the short form mRNA was mainly expressed in the testis, the expression level of the long form mRNA was high in the ovary. These results suggest that ZFAND3 has physiological functions related to germ cell maturation and regulatory mechanisms that differ between the testis and ovary.

Coordinated post-transcriptional regulation of Hsp70.3 gene expression by microRNA and alternative polyadenylation

Heat shock protein 70 (Hsp70) is well documented to possess general cytoprotective properties in protecting the cell against stressful and noxious stimuli. We have recently shown that expression of the stress-inducible Hsp70.3 gene in the myocardium in response to ischemic preconditioning is NF-κB-dependent and necessary for the resulting late phase cardioprotection against a subsequent ischemia/reperfusion injury. Here we show that the Hsp70.3 gene product is subject to post-transcriptional regulation through parallel regulatory processes involving microRNAs and alternative polyadenylation of the mRNA transcript. First, we show that cardiac ischemic preconditioning of the in vivo mouse heart results in decreased levels of two Hsp70.3-targeting microRNAs: miR-378* and miR-711. Furthermore, an ischemic or heat shock stimulus induces alternative polyadenylation of the expressed Hsp70.3 transcript that results in the accumulation of transcripts with a shortened 3'-UTR. This shortening of the 3'-UTR results in the loss of the binding site for the suppressive miR-378* and thus renders the alternatively polyadenylated transcript insusceptible to miR-378*-mediated suppression. Results also suggest that the alternative polyadenylation-mediated shortening of the Hsp70.3 3'-UTR relieves translational suppression observed in the long 3'-UTR variant, allowing for a more robust increase in protein expression. These results demonstrate alternative polyadenylation of Hsp70.3 in parallel with ischemic or heat shock-induced up-regulation of mRNA levels and implicate the importance of this process in post-transcriptional control of Hsp70.3 expression.

RNA polymerase II kinetics in polo polyadenylation signal selection

Regulated alternative polyadenylation is an important feature of gene expression, but how gene transcription rate affects this process remains to be investigated. polo is a cell-cycle gene that uses two poly(A) signals in the 3' untranslated region (UTR) to produce alternative messenger RNAs that differ in their 3'UTR length. Using a mutant Drosophila strain that has a lower transcriptional elongation rate, we show that transcription kinetics can determine alternative poly(A) site selection. The physiological consequences of incorrect polo poly(A) site choice are of vital importance; transgenic flies lacking the distal poly(A) signal cannot produce the longer transcript and die at the pupa stage due to a failure in the proliferation of the precursor cells of the abdomen, the histoblasts. This is due to the low translation efficiency of the shorter transcript produced by proximal poly(A) site usage. Our results show that correct polo poly(A) site selection functions to provide the correct levels of protein expression necessary for histoblast proliferation, and that the kinetics of RNA polymerase II have an important role in the mechanism of alternative polyadenylation.

Alternative polyadenylation: a twist on mRNA 3' end formation

Regulation of gene expression by RNA processing mechanisms is now understood to be an important level of control in mammalian cells. Regulation at the level of RNA transcription, splicing, polyadenylation, nucleo-cytoplasmic transport, and translation into polypeptides has been well-studied. Alternative RNA processing events, such as alternative splicing, also have been recognized as key contributors to the complexity of mammalian gene expression. Pre-messenger RNAs (pre-mRNAs) may be polyadenylated in several different ways due to more than one polyadenylation signal, allowing a single gene to encode multiple mRNA transcripts. However, alternative polyadenylation has only recently taken the field as a major player in gene regulation. This review summarizes what is currently known about alternative polyadenylation. It covers results from bioinformatics, as well as those from investigations of viral and tissue-specific studies and, importantly, will set the stage for what is yet to come.

Alternative polyadenylation produces two major transcripts of Alix

The mammalian adaptor protein Alix participates in multiple cellular processes. Since mouse Alix cDNA detects two distinct transcripts of approximately 3.5 and approximately 7.0 kb in various mouse tissues, it is possible that there exist isoforms of Alix protein that perform varied biological functions. In this study, we first demonstrate that four different anti-Alix monoclonal antibodies immunoblot the single Alix protein in nine different mouse tissues. We then show that the two transcripts of 3.2 and 6.4 kb are widely expressed in various human tissues and cell lines. These two transcripts are generated from the same Alix gene localizing at 3p22.3 via alternative polyadenylation, thus containing an identical open reading frame. However, the 3.2-kb transcript is much more active in translation than the 6.4-kb transcript in a randomly selected cell line. These results eliminate the possibility that the two transcript variants encode different isoforms of Alix protein and suggest that alternative polyadenylation is one of the mechanisms controlling Alix protein expression.

Expression profile of AZF genes in testicular biopsies of azoospermic men

BACKGROUND

The Y-chromosome AZF regions include genes whose functions and specific roles in spermatogenesis have not been fully clarified. This study investigated the expression of several AZF (USP9Y, DDX3Y/DDX3Yt1, EIF1AY and PRY) and USP9X transcripts in testicular biopsies of 89 azoospermic men who had been classified by histology and cytology assessments.

METHODS

Expression was analysed by RT-PCR, and some biopsies were evaluated by multiplex RT-PCR. Quantitative PCR was performed in some biopsies to determine the ratio of the testis-specific transcript DDX3Yt1 to the total DDX3Y transcription.

RESULTS

The expression of USP9Y, USP9X and DDX3Y was found in all the specimens tested, whereas DDX3Yt1 expression was diminished or undetectable in several biopsies with impaired spermatogenesis. EIF1AY was detected in all except two of the specimens. Noteworthy, PRY expression was detected mainly in biopsies with germ cells, and this association was significant (P < 0.001). An identical expression profile was obtained by either single or multiplex RT-PCR.

CONCLUSIONS

These findings suggest that PRY is usually expressed in germ cells, whereas the other transcripts are also expressed in testicular somatic cells. The absence of EIF1AY expression might sporadically contribute to azoospermia. The decreased ratio of DDX3Yt1/DDX3Y transcript in impaired spermatogenesis suggests that the DDX3Yt1 transcript is under-expressed in impaired spermatogenesis. The findings contribute to the search and selection of the most valuable gene markers potentially useful as additional tools for predicting complete spermatogenesis by multiplex expression analysis.