mRNA-Seq of testis and liver tissues reveals a testis-specific gene and alternative splicing associated with hybrid male sterility in dzo

Alternative splicing (AS) plays an important role in the co-transcription and post-transcriptional regulation of gene expression during mammalian spermatogenesis. The dzo is the male F1 offspring of an interspecific hybrid between a domestic bull (Bos taurus male) and a yak (Bos grunniens female) which exhibits male sterility. This study aimed to identify the testis-specific genes and AS associated with hybrid male sterility in dzo. The iDEP90 program and rMATS software were used to identify the differentially expressed genes (DEG) and differential alternative splicing genes (DSG) based on RNA-seq data from the liver (n = 9) and testis (n = 6) tissues of domestic cattle, yak, and dzo. Splicing factors (SF) were obtained from the AmiGO2 and the NCBI databases, and Pearson correlation analysis was performed on the differentially expressed SFs and DSGs. We focused on the testis-specific DEGs and DSGs between dzo and cattle and yak. Among the top 3,000 genes with the most significant variations between these 15 samples, a large number of genes showed testis-specific expression involved with spermatogenesis. Cluster analysis showed that the expression levels of these testis-specific genes were dysregulated during mitosis with a burst downregulation during the pachynema spermatocyte stage. The occurrence of AS events in the testis was about 2.5 fold greater than in the liver, with exon skipping being the major AS event (81.89% to 82.73%). A total of 74 DSGs were specifically expressed in the testis and were significantly enriched during meiosis I, synapsis, and in the piRNA biosynthesis pathways. Notably, STAG3 and DDX4 were of the exon skipping type, and DMC1 was a mutually exclusive exon. A total of 36 SFs were significantly different in dzo testis, compared with cattle and yak. DDX4, SUGP1, and EFTUD2 were potential SFs leading to abnormal AS of testis-specific genes in dzo. These results show that AS of testis-specific genes can affect synapsis and the piRNA biosynthetic processes in dzo, which may be important factors associated with hybrid male sterility in dzo. The dzo is the interspecific hybrid male offspring of a domestic cattle (Bos taurus) and a yak (Bos grunniens), and displays arrested spermatogenesis at meiosis that causes male sterility. This study discovered differential alternative splicing of testis-specific genes and revealed the potential splicing factors associated with hybrid male sterility in dzo, highlighting the important role of alternative splicing in their meiotic arrest. The interspecific hybrid offspring of a domestic bull (Bos taurus) and a yak (Bos grunniens) display heterosis in meat and milk production. The hybrid offspring are particularly adaptable to the harsh environments of the Qinghai-Tibet Plateau. However, the male F1 to F3 offspring of this interspecies hybrid are infertile, and spermatogenesis is arrested at meiosis preventing the prolonged utilization of the benefits of heterosis. This study aimed to identify the testis-specific genes and alternative splicing (AS) associated with hybrid male sterility using RNA-Seq data from the liver and testis tissues of domestic cattle, yaks, and their F1 offspring (dzo). The expression of the testis-specific genes became disordered during mitosis and meiosis in dzo. Their testis-specific genes with AS events were enriched during synapsis and in the piRNA biosynthetic processes. In addition, we identified the potential splicing factors associated with abnormal testis-specific AS gene expression in dzo. These results reveal the important role of AS in the meiotic arrest of dzo.