Genome-wide association study for reproduction traits in Colombian Creole Blanco Orejinegro cattle

The profitability of the beef cattle production system relies heavily on reproductive traits. Unfortunately, certain traits, such as age at first calving (AFC), calving interval (CI), and gestation length (GL), can pose challenges in traditional breeding programs because of their low heritability (0.01-0.12) and sex-limited characteristics. Another important aspect is the conservation of the genetic resources of animals adapted to the Colombian regions, which implies the preservation and rational use of the creole breeds in the country market. Therefore, this study aimed to identify genomic regions in the creole cattle breed Blanco Orejinegro (BON) that influence the reproductive traits in females. The dataset comprised 439 animals and 118,116 single-nucleotide polymorphisms' (SNPs) markers. The GS3 program was used to identify the SNP effects employing the BAYES Cπ methodology. The number of SNPs with effect for AFC was 25, 1527 for CI, and 23 for GL. Some of the genes found associated with reproductive and growth traits as well as immune response and environmental adaptation ECE1, EPH, EPHB2, SMARCAL1, IGFBP5, IGFBP2, FCGRT, EGFR, MUL1, PINK1, STPG1, CNGB1, TGFB1, OXTR, IL22RA1, MYOM3, OXTR, CNR2, HIVEP3, CTPS1, CXCL8, FCGRT, MREG, TMEM169, PECR, and MC1R. Our results evidenced a high contribution of the genetic architecture of the Colombian creole cattle breed Blanco Orejinegro that may impact should be included in implementing genetic improvement and conservation programs.

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.

Sperm chromatin integrity and DNA methylation in Norwegian Red bulls of contrasting fertility

In this study, the complexity of chromatin integrity was investigated in frozen-thawed semen samples from 37 sires with contrasting fertility, expressed as 56-day non-return rates (NR56). Protamine deficiency, thiols, and disulfide bonds were assessed and compared with previously published data for DNA fragmentation index (DFI) and high DNA stainability (HDS). In addition, in vitro embryo development and sperm DNA methylation were assessed using semen samples from 16 of these bulls. The percentages of DFI and HDS were negatively associated with NR56 and cleavage rate and positively associated with sperm protamine deficiency (p < 0.05). Significant differences in cleavage and blastocyst rates were observed between bulls of high and low NR56. However, once fertilization occurred, further development into blastocysts was not associated with NR56. The differential methylation analysis showed that spermatozoa from bulls of low NR56 were hypermethylated compared to bulls of high NR56. Pathway analysis showed that genes annotated to differentially methylated cytosines could participate in different biological pathways and have important biological roles related to bull fertility. In conclusion, sperm cells from Norwegian Red bulls of inferior fertility have less compact chromatin structure, higher levels of DNA damage, and are hypermethylated compared with bulls of superior fertility.

Identification of new biomarkers in immune microenvironment of testicular germ cell tumour

To seek novel prognostic biomarkers for testicular germ cell tumour (TGCT) and investigate the tumour immune microenvironment, we identified critical differentially expressed genes (DEGs) by overlapping GSE1818 dataset from Gene Expression Omnibus (GEO). Protein-protein interaction (PPI) network was used to investigate key modules and hub genes. Functional enrichment analysis was performed to investigate the underlying molecular functions of the DEGs in TGCT development and progression. The following survival analysis based on The Cancer Genome Atlas (TCGA) TGCT dataset indicated that AKAP4, SPA17 and TNP1 are correlated with TGCT prognosis. Immunohistochemistry and quantitative real-time polymerase chain reaction verified the down-regulation of the 3 hub genes in TGCT. Gene set enrichment analysis was conducted to further explore the role of the 3 hub genes in TGCT respectively. In addition, TGCT samples had high infiltration of CD8+ T cells, M0 and M1 macrophage cells, and resting myeloid dendritic cells in immune microenvironment. We also constructed the microRNA-gene regulatory networks to identify the key upstream microRNAs in TGCT. In conclusion, our findings indicated that AKAP4, SPA17 and TNP1 are promising biomarkers of TGCT. AKAP4 and TNP1 might regulate immune cells infiltration in immune microenvironment.

Fertility Relevance Probability Analysis Shortlists Genetic Markers for Male Fertility Impairment

Impairment of male fertility is one of the major public health issues worldwide. Nevertheless, genetic causes of male sub- and infertility can often only be suspected due to the lack of reliable and easy-to-use routine tests. Yet, the development of a marker panel is complicated by the large quantity of potentially predictive markers. Actually, hundreds or even thousands of genes could have fertility relevance. Thus, a systematic method enabling a selection of the most predictive markers out of the many candidates is required. As a criterion for marker selection, we derived a gene-specific score, which we refer to as fertility relevance probability (FRP). For this purpose, we first categorized 2,753 testis-expressed genes as either candidate markers or non-candidates, according to phenotypes in male knockout mice. In a parallel approach, 2,502 genes were classified as candidate markers or non-candidates based on phenotypes in men. Subsequently, we conducted logistic regression analyses with evolutionary rates of genes (dN/dS), transcription levels in testis relative to other organs, and connectivity of the encoded proteins in a protein-protein interaction network as covariates. In confirmation of the procedure, FRP values showed the expected pattern, thus being overall higher in genes with known relevance for fertility than in their counterparts without corresponding evidence. In addition, higher FRP values corresponded with an increased dysregulation of protein abundance in spermatozoa of 37 men with normal and 38 men with impaired fertility. Present analyses resulted in a ranking of genes according to their probable predictive power as candidate markers for male fertility impairment. Thus, AKAP4, TNP1, DAZL, BRDT, DMRT1, SPO11, ZPBP, HORMAD1, and SMC1B are prime candidates toward a marker panel for male fertility impairment. Additional candidate markers are DDX4, SHCBP1L, CCDC155, ODF1, DMRTB1, ASZ1, BOLL, FKBP6, SLC25A31, PRSS21, and RNF17. FRP inference additionally provides clues for potential new markers, thereunder TEX37 and POU4F2. The results of our logistic regression analyses are freely available at the PreFer Genes website (https://prefer-genes.uni-mainz.de/).

The roles of cirRNA in the development of germ cells

Circular RNA (CircRNA), a type of endogenous non-coding RNAs (ncRNAs), is generally generated from precursor mRNA (pre-mRNA) by canonical splicing and head-to-tail back splicing. The structure without a polyA tail renders circRNA highly insensitive to ribonuclease. Simultaneously, the distribution of circRNAs is tissue and developmental stage-specific. There are five potential biological functions of circRNAs: 1) promote transcription of their parental genes; 2) function as a miRNA sponge; 3) RNA binding protein (RBP) sponge; 4) encode protein; 5) act as an mRNA trap. Recently, circRNA has attracted attention because studies have shown that circRNAs are associated with follicular development, ovarian senescence, spermatogenesis, and germ cell development process, suggesting that circRNAs may function in germ cells regulation. The investigation of circRNAs in germ cells will provide an excellent opportunity to understand its potential molecular basis, and potentially improving reproduction status in human. In this article, the relationship between circRNA and germ cell development will be discussed.

Pig SOX9: Expression profiles of Sertoli cell (SCs) and a functional 18 bp indel affecting testis weight

Sex determining region Y-box 9 (SOX9), an important member of the SRY- type HMGbox (SOX) gene family, plays an important role in the regulation of mammalian reproduction, including sex differentiation during the embryonic development stage and spermatogenesis after birth. To explore the roles of polymorphism and expression of the SOX9 gene in the development of testes, we analyzed the indel of SOX9 in pigs and the corresponding expression level of the SOX9 gene in 7-day and 5-month-old porcine Sertoli cells. Results revealed that the DD haplotype of SOX9 gene as well as the ID genotype were significantly associated with larger testicular weight, while the II haplotype was closely related to the smaller testicular weight. More importantly, the SOX9 gene expression of ID genotyped group was significantly higher than that in II genotyped group. Our results first revealed that the indel polymorphism and expression of SOX9 were significantly associated with pig reproduction traits indicating the critical roles of SOX9 gene in testes development. The study provides a new clue for understanding the regulation of animal reproductive activities.

Functional haplotypes of ARID4A affect promoter activity and semen quality of bulls

The AT-rich interaction domain 4 A (ARID4A) has an important role in regulating Sertoli cell function and male fertility. Its molecular mechanisms, however, remain largely unknown. In this study, two single nucleotide polymorphisms (SNPs) (g.53 G > T, ss 1966531596, and g.826 G > A, rs 210809648) were identified in the promoter region of ARID4A in 215 Chinese Holstein bulls using polymerase chain reaction (PCR)-restriction fragment length polymorphism and created restriction site-PCR. Results revealed that bulls with g.53 G > T-GG and g.826 G > A-G G genotype exhibited higher sperm deformity rate than those with g.53 G > T-TT and g.826 G > A-AA genotype (P <  0.01). Furthermore, three haplotypes (H1 (GG), H3 (TG), H4 (TA)) and six haplotype combinations (H1H1, H1H3, H1H4, H3H3, H3H4, H4H4) were obtained. The bulls with H4H4 exhibited lower sperm deformity rate than those with H1H1 and H1H3 (P <  0.05). In addition, results of bioinformatics analysis revealed that ARID4A has two promoters and that two SNPs of ARID4A are located in transcription factor binding sites. Compared with g.53 G > T-G and g.826 G > A-G allele, there was a greater fluorescence intensity in g.53 G > T-T and g.826 G > A-A allele by transient transfection in MLTC-1 cells and the luciferase report assay. qRT-PCR indicated the ARID4A expression was greater in bull spermatozoa with H4H4 haplotype combination than those with H1H1 haplotype combination (P < 0.05). Results of the present study indicate that g.53 G > T and g.826 G > A are functional mutations that are involved in regulation of ARID4A gene expression by affecting promoter activity and thus semen quality of Chinese Holstein bulls.

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.

Review: Sperm-oocyte interactions and their implications for bull fertility, with emphasis on the ubiquitin-proteasome system

Fertilization is an intricate cascade of events that irreversibly alter the participating male and female gamete and ultimately lead to the union of paternal and maternal genomes in the zygote. Fertilization starts with sperm capacitation within the oviductal sperm reservoir, followed by gamete recognition, sperm-zona pellucida interactions and sperm-oolemma adhesion and fusion, followed by sperm incorporation, oocyte activation, pronuclear development and embryo cleavage. At fertilization, bull spermatozoon loses its acrosome and plasma membrane components and contributes chromosomes, centriole, perinuclear theca proteins and regulatory RNAs to the zygote. While also incorporated in oocyte cytoplasm, structures of the sperm tail, including mitochondrial sheath, axoneme, fibrous sheath and outer dense fibers are degraded and recycled. The ability of some of these sperm contributed components to give rise to functional zygotic structures and properly induce embryonic development may vary between bulls, bearing on their reproductive performance, and on the fitness, health, fertility and production traits of their offspring. Proper functioning, recycling and remodeling of gamete structures at fertilization is aided by the ubiquitin-proteasome system (UPS), the universal substrate-specific protein recycling pathway present in bovine and other mammalian oocytes and spermatozoa. This review is focused on the aspects of UPS relevant to bovine fertilization and bull fertility.

MiR-27a-5p Increases Steer Fat Deposition Partly by Targeting Calcium-sensing Receptor (CASR)

Castration increases fat deposition, improving beef quality in cattle. Here, the steer group exhibited a significantly higher intramuscular fat (IMF) content than the bull group. To determine the potential roles of microRNAs (miRNAs) in castration-induced fat deposition, differential expression patterns of miRNA in liver tissue were investigated in bulls and steers. A total of 7,827,294 clean reads were obtained from the bull liver library, and 8,312,483 were obtained from the steer liver library; 452 conserved bovine miRNAs and 20 novel miRNAs were identified. The results showed that the expression profiles of miRNA in liver tissue were changed by castration, and 12 miRNAs that were differentially expressed between bulls and steers were identified. Their target genes were majorly involved in the metabolic, PI3K-Akt, and MAPK signaling pathways. Furthermore, six differentially expressed miRNAs were validated by quantitative real-time PCR, and luciferase reporter assays verified that calcium-sensing receptor (CASR) was the direct target of miR-27a-5p. Meantime, we found that the expression level of CASR was significantly higher in steers than in bulls, and revealed that CASR gene silencing in bovine hepatocytes significantly inhibited triacylglycerol (TAG) accumulation and reduced secretion of very low density lipoprotein (VLDL). These results obtained in the liver indicate that miR-27a-5p may increase fat deposition partly by targeting CASR in steers.

Genetic effects of PDGFRB and MARCH1 identified in GWAS revealing strong associations with semen production traits in Chinese Holstein bulls

Background

Using a genome-wide association study strategy, our previous study discovered 19 significant single-nucleotide polymorphisms (SNPs) related to semen production traits in Chinese Holstein bulls. Among them, three SNPs were within or close to the phosphodiesterase 3A (PDE3A), membrane associated ring-CH-type finger 1 (MARCH1) and platelet derived growth factor receptor beta (PDGFRB) genes. The present study was designed with the objectives of identifying genetic polymorphism of the PDE3A, PDGFRB and MARCH1 genes and their effects on semen production traits in a Holstein bull population.

Results

A total of 20 SNPs were detected and genotyped in 730 bulls. Association analyses using de-regressed estimated breeding values of each semen production trait revealed four statistically significant SNPs for one or more semen production traits (P < 0.05): one SNP was located downstream of PDGFRB and three SNPs were located in the promoter of MARCH1. Interestingly, for MARCH1, haplotype-based analysis revealed significant associations of haplotypes with semen volume per ejaculate. Furthermore, high expression of the MARCH1 gene was observed in sperm cells. One SNP (rs43445726) in the regulatory region of MARCH1 had a significant effect on gene expression.

Conclusion

Our study demonstrated the significant associations of genetic variants of the PDGFRB and MARCH1 genes with semen production traits. The identified SNPs may serve as genetic markers to optimize breeding programs for semen production traits in Holstein bull populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0527-1) contains supplementary material, which is available to authorized users.

Functional SNPs of INCENP Affect Semen Quality by Alternative Splicing Mode and Binding Affinity with the Target Bta-miR-378 in Chinese Holstein Bulls

Inner centromere protein (INCENP) plays an important role in mitosis and meiosis as the main member of chromosomal passenger protein complex (CPC). To investigate the functional markers of the INCENP gene associated with semen quality, the single nucleotide polymorphisms (SNPs) g.19970 A>G and g.34078 T>G were identified and analyzed. The new splice variant INCENP-TV is characterized by the deletion of exon 12. The g.19970 A>G in the exonic splicing enhancer (ESE) motif region results in an aberrant splice variant by constructing two minigene expression vectors using the pSPL3 exon capturing vector and transfecting vectors into MLTC-1 cells. INCENP-TV was more highly expressed than INCENP-reference in adult bull testes. The g.34078 T>G located in the binding region of bta-miR-378 could affect the expression of INCENP, which was verified by luciferase assay. To analyze comprehensively the correlation of SNPs with sperm quality, haplotype combinations constructed by g.19970 A>G and g.34078 T>G, as well as g.-692 C>T and g.-556 G>T reported in our previous studies, were analyzed. The bulls with H1H12 and H2H2 exhibited a higher ejaculate volume than those with H2H10 and H9H12, respectively (P < 0.05). Bulls with H11H11 and H2H10 exhibited higher initial sperm motility than those with H2H2 (P < 0.05). The expression levels of INCENP in bulls with H1H12 and H11H11 were significantly higher than those in bulls with H9H12 (P < 0.05), as determined by qRT-PCR. Findings suggest that g.19970 A>G and g.34078 T>G in INCENP both of which appear to change the molecular and biological characteristics of the mRNA transcribed from the locus may serve as a biomarkers of male bovine fertility by affecting alternative splicing mode and binding affinity with the target bta-miR-378.

Association of polymorphisms at the microRNA binding site of the caprine KITLG 3'-UTR with litter size

This study identified three novel single nucleotide polymorphisms (SNPs) (c.1389C > T, c.1457A > C and c.1520G > A) in the caprine KITLG 3'-UTR through DNA sequencing. The three SNP loci were closely linked in Guanzhong dairy (GD) goats. Two alleles of the c.1457A > C SNP introduced two miRNA sites (chi-miR-204-5p and chi-miR-211). Individuals with combined genotype TT-CC-AA had a higher litter size compared with those with combined genotypes CC-AA-GG, TC-CC-GA and CC-AC-GG (P < 0.05). Luciferase assays showed that chi-miR-204-5p and chi-miR-211 suppressed luciferase expression in the presence of allele 1457A compared with negative control (NC) and allele 1457C (P < 0.05). Western blot revealed that KITLG significantly decreased in the granulosa cells (GCs) of genotype AA compared with that in the GCs of genotype CC and NC (P < 0.05). The KITLG mRNA levels of the CC-AA-GG carriers significantly decreased compared with those of the TT-CC-AA, TC-CC-GA and CC-AC-GG carriers. In addition, cell proliferation was reduced in haplotype C-A-G GCs compared with that in haplotype T-C-A GCs. These results suggest that SNPs c.1389C > T, c.1457A > C and c.1520G > A account for differences in the litter size of GD goats because chi-miR-204-5p and chi-miR-211 could change the expression levels of the KITLG gene and reduce GC proliferation.