Exploiting RNA-sequencing data from the porcine testes to identify the key genes involved in spermatogenesis in Large White pigs

Comprehensive analysis of differentially expressed mRNAs, circRNAs, and miRNAs and their ceRNA network in the testis of cattle-yak, yak, and cattle

Cattle-yak is a hybrid offspring resulting from the crossbreeding of yak and cattle, and it exhibits substantial heterosis in production performance. However, male sterility in cattle-yak remains a concern. Reports suggest that noncoding RNAs are involved in the regulation of spermatogenesis. Therefore, in this study, we comprehensively compared testicular transcription profiles among cattle, yak, and cattle-yak. Numerous differentially expressed genes (DEGs), differentially expressed circRNAs (DECs), and differentially expressed miRNAs (DEMs) were identified in the intersection of two comparison groups, namely cattle versus cattle-yak and yak versus cattle-yak, with the number of DEGs, DECs, and DEMs being 4968, 360, and 59, respectively. The DEGs in cattle-yaks, cattle, and yaks were mainly associated with spermatogenesis, male gamete generation, and sexual reproduction. Concurrently, GO and KEGG analyses indicated that DEC host genes and DEM source genes were involved in the regulation of spermatogenesis. The construction of a potential competing endogenous RNA network revealed that some differentially expressed noncoding RNAs may be involved in regulating the expression of genes related to testicular spermatogenesis, including miR-423-5p, miR-449b, miR-34b/c, and miR-15b, as well as previously unreported miR-6123 and miR-1306, along with various miRNA-circRNA interaction pairs. This study serves as a valuable reference for further investigations into the mechanisms underlying male sterility in cattle-yaks.

Complete characterization of the yak testicular development using accurate full-length transcriptome sequencing

Alternative splicing is a prevalent phenomenon in testicular tissues. Due to the low assembly accuracy of short-read RNA sequencing technology in analyzing post-transcriptional regulatory events, full-length (FL) transcript sequencing is highly demanded to accurately determine FL splicing variants. In this study, we performed FL transcriptome sequencing of testicular tissues from 0.5, 1.5, 2.5, and 4-year-old yaks and 4-year-old cattle-yaks using Oxford Nanopore Technologies. The obtained sequencing data were predicted to have 47,185 open reading frames (ORFs), including 26,630 complete ORFs, detected 7645 fusion transcripts, 15,355 alternative splicing events, 25,798 simple sequence repeats, 7628 transcription factors, and 35,503 long non-coding RNAs. A total of 40,038 novel transcripts were obtained from the sequencing data, and the proportion was almost close to the number of known transcripts identified. Structural analysis and functional annotation of these novel transcripts resulted in the successful annotation of 9568 transcripts, with the highest and lowest annotation numbers in the Nr and KOG databases, respectively. Weighted gene co-expression network analysis revealed the key regulatory pathways and hub genes at various stages of yak testicular development. Our findings enhance our comprehension of transcriptome complexity, contribute to genome annotation refinement, and provide foundational data for further investigations into male sterility in cattle-yaks.

SPAG6 c.900 T>C affects boar semen quality and blood-testis barrier function by creating a new splice acceptor site

Sperm associated antigen 6 (SPAG6) acts as a scaffolding protein in the center of the flagellar axoneme and has an impact on the maturation of the motility of mammalian sperm flagella and the maintenance of sperm structure. In our previous research, SPAG6 c.900 T>C in exon 7 and exon 7 skipped transcript was identified by analyzing RNA-seq data of testicular tissues from 60 day (sexually immature) and 180 day (sexually mature) Large White boars. Herein, we found porcine SPAG6 c.900 T>C to be associated with semen quality traits in Duroc, Large White and Landrace pigs. SPAG6 c.900 C can generate a new splice acceptor site, inhibit the occurrence of SPAG6 exon 7 skipping to a certain extent, thereby promote the growth of Sertoli cells and maintain the normal blood-testis barrier function. This study provides new insights into the molecular regulation of spermatogenesis and a new genetic marker for the improvement of semen quality in pigs.

BCL2-associated athanogene 6 exon24 contributes to testosterone synthesis and male fertility in mammals

Objectives

BCL2‐associated athanogene 6 (BAG6) plays critical roles in spermatogenesis by maintaining testicular cell survival. Our previous data showed porcine BAG6 exon24‐skipped transcript is highly expressed in immature testes compared with mature testes. The objective of this study is to reveal the functional significance of BAG6 exon24 in mammalian spermatogenesis.

Materials and Methods

CRISPR/Cas9 system was used to generate Bag6 exon24 knockout mice. Testes and cauda epididymal sperm were collected from mice. TMT proteomics analysis was used to discover the protein differences induced by Bag6 exon24 deletion. Testosterone enanthate was injected into mice to generate a high‐testosterone mice model. H&E staining, qRT‐PCR, western blotting, vector/siRNA transfection, immunofluorescence, immunoprecipitation, transmission electron microscopy, TUNEL and ELISA were performed to investigate the phenotypes and molecular basis.

Results

Bag6 exon24 knockout mice show sub‐fertility along with partially impaired blood‐testis barrier, increased apoptotic testicular cell rate and abnormal sperm morphology. Endoplasmic reticulum stress occurs in Bag6 exon24‐deficient testes and sterol regulatory element‐binding transcription factor 2 is activated; as a result, cytochrome P450 family 51 subfamily A member 1 expression is up‐regulated, which causes a high serum testosterone level. Additionally, serine/arginine‐rich splicing factor 1 down‐regulates BAG6 exon24‐skipped transcripts in porcine Sertoli cells by binding to 35–51 nt on BAG6 exon24 via its N‐terminal RNA‐recognition domain.

Conclusions

Our findings reveal the critical roles of BAG6 exon24 in testosterone biosynthesis and male fertility, which provides new insights into the regulation of spermatogenesis and pathogenesis of subfertility in mammals.

Identification of Key Genes Related to Postnatal Testicular Development Based on Transcriptomic Data of Testis in Hu Sheep

The selection of testis size can improve the reproductive capacity of livestock used for artificial insemination and has been considered as an important strategy for accelerating the breeding process. Although much work has been done to investigate the mechanisms of testis development in various species, there is little information available in regard to the differences in transcriptomic profiling of sheep testes at different developmental stages. In this work, we aimed to identify differentially expressed genes (DEGs) by RNA-Seq in sheep during different growth stages, including 0 month old (infant, M0), 3 months old (puberty, M3), 6 months old (sexual maturity, M6) and 12 months old (body maturity, M12). A total of 4,606 (2,381 up and 2,225 down), 7,500 (4,368 up and 3,132 down), 15 (8 up and seven down) DEGs were identified in M3_vs_M0, M6_vs_M3, and M12_vs_M6 comparison, respectively. Of which, a number of genes were continuously up-regulated and down-regulated with testicular development, including ODF3, ZPBP1, PKDREJ, MYBL1, PDGFA, IGF1, LH, INSL3, VIM, AMH, INHBA, COL1A1, COL1A2, and INHA. GO analysis illustrated that DEGs were mainly involved in testis development and spermatogenesis. KEGG analysis identified several important pathways and verified several reproduction-associated DEGs such as COL1A1, COL1A2, PDGFA, and IGF1. In addition, two gene modules highly associated with testis development and core genes with testis size were identified using weighted gene co-expression network analysis (WGCNA), including hub genes positively associated with testis size such as RANBP9, DNAH17, SPATA4, CIB4 and SPEM1, and those negatively associated with testis size such as CD81, CSK, PDGFA, VIM, and INHBA. This study comprehensively identified key genes related to sheep testicular development, which may provide potential insights for understanding male fertility and better guide in animal breeding.

Circular RNA Expression and Regulation Profiling in Testicular Tissues of Immature and Mature Wandong Cattle (Bos taurus)

Wandong cattle are an autochthonous Chinese breed used extensively for beef production. The breed tolerates extreme weather conditions and raw feed and is resistant to tick-borne diseases. However, the genetic basis of testis development and sperm production as well as breeding management is not well established in local cattle. Therefore, improving the reproductive efficiency of bulls via genetic selection is crucial as a single bull can breed thousands of cows through artificial insemination (AI). Testis development and spermatogenesis are regulated by hundreds of genes and transcriptomes. However, circular RNAs (circRNAs) are the key players in many biological developmental processes that have not been methodically described and compared between immature and mature stages in Bovine testes. In this study, we performed total RNA-seq and comprehensively analyzed the circRNA expression profiling of the testis samples of six bulls at 3 years and 3 months of developmental age. In total, 17,013 circRNAs were identified, of which 681 circRNAs (p-adjust < 0.05) were differentially expressed (DE). Among these DE circRNAs, 579 were upregulated and 103 were downregulated in calf and bull testes. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses revealed that the identified target genes were classified into three broad functional categories, including biological process, cellular component, and molecular function, and were enriched in the lysine degradation, cell cycle, and cell adhesion molecule pathways. The binding interactions between DE circRNAs and microRNAs (miRNAs) were subsequently constructed using bioinformatics approaches. The source genes ATM, CCNA1, GSK3B, KMT2C, KMT2E, NSD2, SUCLG2, QKI, HOMER1, and SNAP91 were found to be actively associated with bull sexual maturity and spermatogenesis. In addition, a real-time quantitative polymerase chain reaction (RT-qPCR) analysis showed a strong correlation with the sequencing data. Moreover, the developed model of Bovine testes in the current study provides a suitable framework for understanding the mechanism of circRNAs in the development of testes and spermatogenesis.

Selection signatures in two oldest Russian native cattle breeds revealed using high-density single nucleotide polymorphism analysis

Native cattle breeds can carry specific signatures of selection reflecting their adaptation to the local environmental conditions and response to the breeding strategy used. In this study, we comprehensively analysed high-density single nucleotide polymorphism (SNP) genotypes to characterise the population structure and detect the selection signatures in Russian native Yaroslavl and Kholmogor dairy cattle breeds, which have been little influenced by introgression with transboundary breeds. Fifty-six samples of pedigree-recorded purebred animals, originating from different breeding farms and representing different sire lines, of the two studied breeds were genotyped using a genome-wide bovine genotyping array (Bovine HD BeadChip). Three statistical analyses—calculation of fixation index (F_ST) for each SNP for the comparison of the pairs of breeds, hapFLK analysis, and estimation of the runs of homozygosity (ROH) islands shared in more than 50% of animals—were combined for detecting the selection signatures in the genome of the studied cattle breeds. We confirmed nine and six known regions under putative selection in the genomes of Yaroslavl and Kholmogor cattle, respectively; the flanking positions of most of these regions were elucidated. Only two of the selected regions (localised on BTA 14 at 24.4–25.1 Mbp and on BTA 16 at 42.5–43.5 Mb) overlapped in Yaroslavl, Kholmogor and Holstein breeds. In addition, we detected three novel selection sweeps in the genome of Yaroslavl (BTA 4 at 4.74–5.36 Mbp, BTA 15 at 17.80–18.77 Mbp, and BTA 17 at 45.59–45.61 Mbp) and Kholmogor breeds (BTA 12 at 82.40–81.69 Mbp, BTA 15 at 16.04–16.62 Mbp, and BTA 18 at 0.19–1.46 Mbp) by using at least two of the above-mentioned methods. We expanded the list of candidate genes associated with the selected genomic regions and performed their functional annotation. We discussed the possible involvement of the identified candidate genes in artificial selection in connection with the origin and development of the breeds. Our findings on the Yaroslavl and Kholmogor breeds obtained using high-density SNP genotyping and three different statistical methods allowed the detection of novel putative genomic regions and candidate genes that might be under selection. These results might be useful for the sustainable development and conservation of these two oldest Russian native cattle breeds.

The Function of Pre-mRNA Alternative Splicing in Mammal Spermatogenesis

Alternative pre-mRNA splicing plays important roles in co-transcriptional and post-transcriptional regulation of gene expression functioned during many developmental processes, such as spermatogenesis. The studies focusing on alternative splicing on spermatogenesis supported the notion that the development of testis is regulated by a higher level of alternative splicing than other tissues. Here, we aim to review the mechanisms underlying alternative splicing, particularly the splicing variants functioned in the process of spermatogenesis and the male infertility. There are five points regarding the alternative splicing including ⅰ) a brief introduction of alternative pre-mRNA splicing; ⅱ) the alternative splicing events in spermatogenesis-associated genes enriched in different stages of spermatogenesis; ⅲ) the mechanisms of alternative splicing regulation, such as splicing factors and m⁶A demethylation; ⅳ) the splice site recognition and alternative splicing, including the production and degradation of abnormal transcripts caused by gene variations and nonsense-mediated mRNA decay, respectively; ⅴ) abnormal alternative splicing correlated with male infertility. Taking together, this review highlights the impacts of alternative splicing and splicing variants in mammal spermatogenesis and provides new insights of the potential application of the alternative splicing into the therapy of male infertility.

Weighted single-step GWAS and gene network analysis reveal new candidate genes for semen traits in pigs

Background

In recent years, there has been increased interest in the study of the molecular processes that affect semen traits. In this study, our aim was to identify quantitative trait loci (QTL) regions associated with four semen traits (motility, progressive motility, number of sperm cells per ejaculate and total morphological defects) in two commercial pig lines (L1: Large White type and L2: Landrace type). Since the number of animals with both phenotypes and genotypes was relatively small in our dataset, we conducted a weighted single-step genome-wide association study, which also allows unequal variances for single nucleotide polymorphisms. In addition, our aim was also to identify candidate genes within QTL regions that explained the highest proportions of genetic variance. Subsequently, we performed gene network analyses to investigate the biological processes shared by genes that were identified for the same semen traits across lines.

Results

We identified QTL regions that explained up to 10.8% of the genetic variance of the semen traits on 12 chromosomes in L1 and 11 chromosomes in L2. Sixteen QTL regions in L1 and six QTL regions in L2 were associated with two or more traits within the population. Candidate genes SCN8A, PTGS2, PLA2G4A, DNAI2, IQCG and LOC102167830 were identified in L1 and NME5, AZIN2, SPATA7, METTL3 and HPGDS in L2. No regions overlapped between these two lines. However, the gene network analysis for progressive motility revealed two genes in L1 (PLA2G4A and PTGS2) and one gene in L2 (HPGDS) that were involved in two biological processes i.e. eicosanoid biosynthesis and arachidonic acid metabolism. PTGS2 and HPGDS were also involved in the cyclooxygenase pathway.

Conclusions

We identified several QTL regions associated with semen traits in two pig lines, which confirms the assumption of a complex genetic determinism for these traits. A large part of the genetic variance of the semen traits under study was explained by different genes in the two evaluated lines. Nevertheless, the gene network analysis revealed candidate genes that are involved in shared biological pathways that occur in mammalian testes, in both lines.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0412-z) contains supplementary material, which is available to authorized users.

Identification and characterization of circular RNAs in Qinchuan cattle testis

Circular RNA (circRNA) is a new class of non-coding RNA that has recently attracted researchers' interest. Studies have demonstrated that circRNA can function as microRNA sponges or competing endogenous RNAs. Although circRNA has been explored in some species and tissues, the genetic basis of testis development and spermatogenesis in cattle remains unknown. We performed ribo-depleted total RNA-Seq to detect circRNA expression profiles of neonatal (one week old) and adult (4 years old) Qinchuan cattle testes. We obtained 91 112 596 and 80 485 868 clean reads and detected 21 753 circRNAs. A total of 4248 circRNAs were significantly differentially expressed between neonatal and adult cattle testes. Among these circRNAs, 2225 were upregulated, and 2023 were downregulated in adult cattle testis. Genomic feature, length distribution and other characteristics of the circRNAs in cattle testis were studied. Moreover, Gene Ontology and KEGG pathway analyses were performed for source genes of circRNAs. These source genes were mainly involved in tight junction, adherens junction, TGFβ signalling pathway and reproduction, such as PIWIL1, DPY19L2, SLC26A8, IFT81, SMC1B, IQCG and TTLL5. CircRNA expression patterns were validated by RT-qPCR. Our discoveries provide a solid foundation for the identification and characterization of key circRNAs involved in testis development or spermatogenesis.

Identification and characterization of circular RNAs in Qinchuan cattle testis

Circular RNA (circRNA) is a new class of non-coding RNA that has recently attracted researchers' interest. Studies have demonstrated that circRNA can function as microRNA sponges or competing endogenous RNAs. Although circRNA has been explored in some species and tissues, the genetic basis of testis development and spermatogenesis in cattle remains unknown. We performed ribo-depleted total RNA-Seq to detect circRNA expression profiles of neonatal (one week old) and adult (4 years old) Qinchuan cattle testes. We obtained 91 112 596 and 80 485 868 clean reads and detected 21 753 circRNAs. A total of 4248 circRNAs were significantly differentially expressed between neonatal and adult cattle testes. Among these circRNAs, 2225 were upregulated, and 2023 were downregulated in adult cattle testis. Genomic feature, length distribution and other characteristics of the circRNAs in cattle testis were studied. Moreover, Gene Ontology and KEGG pathway analyses were performed for source genes of circRNAs. These source genes were mainly involved in tight junction, adherens junction, TGFβ signalling pathway and reproduction, such as PIWIL1, DPY19L2, SLC26A8, IFT81, SMC1B, IQCG and TTLL5. CircRNA expression patterns were validated by RT-qPCR. Our discoveries provide a solid foundation for the identification and characterization of key circRNAs involved in testis development or spermatogenesis.

Identification and characterisation of microRNAs and Piwi-interacting RNAs in cockerels' spermatozoa by Solexa sequencing

1. There has been substantial research focused on the roles of microRNAs (miRNAs) and Piwi-interacting RNAs (piRNAs) derived from mammalian spermatozoa; however, comparatively little is known about the role of spermatozoa-derived miRNAs and piRNAs within breeding cockerels' spermatozoa. 2. A small RNA library of cockerels' spermatozoa was constructed using Illumina high-throughput sequencing technology. Unique sequences with lengths of 18-26 nucleotides were mapped to miRBase 21.0 and unique sequences with lengths of 25-37 nucleotides were mapped to a piRNA database. A total of 1311 miRNAs and 2448 potential piRNAs were identified. Based on stem-loop qRT-PCR, 8 miRNAs were validated. 3. Potential target genes of the abundant miRNAs were predicted, and further Kyoto Encyclopedia of Genes and Genomes database (KEGG) and Gene Ontology (GO) analyses were performed, which revealed that some candidate miRNAs were involved in the spermatogenesis process, spermatozoa epigenetic programming and further embryonic development. 5. GO and KEGG analyses based on mapping genes of expressed piRNAs were performed, which revealed that spermatozoal piRNAs could play important regulatory roles in embryonic development of offspring. 6. The search for endogenous spermatozoa miRNAs and piRNAs will contribute to a preliminary database for functional and molecular mechanistic studies in embryonic development and spermatozoa epigenetic programming.

RNA sequencing reveals candidate genes and polymorphisms related to sperm DNA integrity in testis tissue from boars

Background

Sperm DNA is protected against fragmentation by a high degree of chromatin packaging. It has been demonstrated that proper chromatin packaging is important for boar fertility outcome. However, little is known about the molecular mechanisms underlying differences in sperm DNA fragmentation. Knowledge of sequence variation influencing this sperm parameter could be beneficial in selecting the best artificial insemination (AI) boars for commercial production. The aim of this study was to identify genes differentially expressed in testis tissue of Norwegian Landrace and Duroc boars, with high and low sperm DNA fragmentation index (DFI), using transcriptome sequencing.

Results

Altogether, 308 and 374 genes were found to display significant differences in expression level between high and low DFI in Landrace and Duroc boars, respectively. Of these genes, 71 were differentially expressed in both breeds. Gene ontology analysis revealed that significant terms in common for the two breeds included extracellular matrix, extracellular region and calcium ion binding. Moreover, different metabolic processes were enriched in Landrace and Duroc, whereas immune response terms were common in Landrace only. Variant detection identified putative polymorphisms in some of the differentially expressed genes. Validation showed that predicted high impact variants in RAMP2, GIMAP6 and three uncharacterized genes are particularly interesting for sperm DNA fragmentation in boars.

Conclusions

We identified differentially expressed genes between groups of boars with high and low sperm DFI, and functional annotation of these genes point towards important biochemical pathways. Moreover, variant detection identified putative polymorphisms in the differentially expressed genes. Our results provide valuable insights into the molecular network underlying DFI in pigs.

Electronic supplementary material

The online version of this article (10.1186/s12917-017-1279-x) contains supplementary material, which is available to authorized users.

miR-638 Inhibits immature Sertoli cell growth by indirectly inactivating PI3K/AKT pathway via SPAG1 gene

Numerous studies have demonstrated that microRNAs (miRNAs) play important roles in cell growth, apoptosis and spermatogenesis. Our previous study showed that miR-638 was differentially expressed in sexually immature and mature testes of Large White boars. Here we reported that sperm-associated antigen 1 (SPAG1) was a direct target gene of miR-638. Moreover, miR-638 inhibited cell proliferation and cell cycle, and promoted apoptosis of porcine immature Sertoli cells. Key genes including phosphorylated phosphatidylinositide 3-kinases (p-PI3K) and phosphorylated serine/ threonine kinase (p-AKT) in PI3K/AKT pathway as well as cell cycle factors including c-MYC, cyclin-D1 (CCND1), cyclin-E1 (CCNE1) and cyclin-dependent kinase 4 (CDK4) were all significantly down-regulated after overexpression of miR-638 or RNAi of SPAG1. Notably, mRNA levels of SRY-related HMG-box 2 (SOX2) and POU domain, class 5, transcription factor 1 (POU5F1) essential for spermatogonia proliferation were significantly suppressed in SPAG1 siRNA- transfected ST cells. This study suggests that miR-638 regulates immature Sertoli cell growth and apoptosis by targeting SPAG1 gene which can indirectly inactivate PI3K/AKT pathway, and plays roles in pig spermatogenesis.

miRNA and piRNA expression profiles of breeder cock testes detected by next-generation sequencing

miRNAs are small non-coding regulatory RNAs that play key roles in diverse biological processes. In this study, we used the Solexa sequencing technique to profile miRNAs in breeder cock testes to illustrate their functions. A total of 663 co-expressed miRNAs and 3,180 co-expressed piRNAs were detected in three libraries. Based on Mir-X^™ miRNA qRT-PCR, three miRNAs representing low, medium and high expression levels according to the sequencing results were selected randomly to validate the miRNAs' expression profiles. Results suggested that the miRNA expression profiles data could represent actual miRNA expression levels. Moreover, target genes prediction of the co-expressed miRNAs and further Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses were performed, which revealed that some candidate miRNAs were involved in the regulation of the spermatogenesis process, spermatozoa function and testicular metabolism. In conclusion, we provided a useful resource for further elucidation of the miRNAs' regulatory role in spermatogenesis, contributing to a preliminary database for functional and molecular mechanistic studies in testicular metabolism, spermatogenesis and other testes functions.

Genome-wide differential expression of genes and small RNAs in testis of two different porcine breeds and at two different ages

Some documented evidences proved small RNAs (sRNA) and targeted genes are involved in mammalian testicular development and spermatogenesis. However, the detailed molecular regulation mechanisms of them remain largely unknown so far. In this study, we obtained a total of 10,716 mRNAs, 67 miRNAs and 16,953 piRNAs which were differentially expressed between LC and LW pig breeds or between the two sexual maturity stages. Of which, we identified 16 miRNAs and 28 targeted genes possibly related to spermatogenesis; 14 miRNA and 18 targeted genes probably associated with cell adhesion related testis development. We also annotated 579 piRNAs which could potentially regulate cell death, nucleosome organization and other basic biology process, which implied that those piRNAs might be involved in sexual maturation difference. The integrated network analysis results suggested that some differentially expressed genes were involved in spermatogenesis through the ECM-receptor interaction, focal adhesion, Wnt and PI3K-Akt signaling pathways, some particular miRNAs have the negative regulation roles and some special piRNAs have the positive and negative regulation roles in testicular development. Our data provide novel insights into the molecular expression and regulation similarities and diversities of spermatogenesis and testicular development in different pig breeds at different stages of sexual maturity.