Identification of key genes associated with spermatogenesis arrest in fox hybrids using weighted gene co-expression network analysis

Abstract The silver fox and the blue fox represent different genera, but produce viable offspring. Although these hybrids show obvious heterosis, they are completely sterile due to spermatogenic arrest at the early stages of spermatogenesis, especially mitosis and meiosis I; the hybrids produce few spermatogonia and primary spermatocytes, and no secondary spermatocytes. Although the mechanisms of spermatogenic arrest have been well investigated, transcriptomic differences between hybrid and the pure-species testes have not clarified. In the present study, we used RNA sequencing (RNA-Seq) to generate testicular transcriptomic profiles for silver foxes, blue foxes, and reciprocal hybrids during the pre-breeding period and the breeding season. In total, 1,344,022 transcripts (≥200 bp) were generated; 1,057,724 genes were obtained; and 33,423 genes were shown to have fragments per kilobase of transcript per million mapped reads (FPKM) > 0.3. To identify the hub genes associated with spermatogenesis arrest, weighted gene co-expression network analysis (WGCNA) was used. Nine modules were explored. Genes in only a single module were consistently downregulated in the hybrids as compared to the pure species; these genes were significantly associated with fox hybrid male infertility. Six of the genes in this module (CATSPERD, DMRTC2, RNF17, NME5, SPEF2, SPINK2) also play key roles in mitosis and meiosis during spermatogenesis. Therefore, these six genes might be associated with fox hybrid male infertility.