Six novel single-nucleotide polymorphisms in SPAG11 gene and their association with sperm quality traits in Chinese Holstein bulls

Sperm-Borne Mitochondrial Activity Influenced by Season and Age of Holstein Bulls

Sperm mitochondria are vital organelles for energy production and pre- and post-fertilization sperm functions. The potential influence of the age of the bull and season on the sperm-borne mitochondrial copy number and the transcription activity has not yet been investigated. Therefore, the expression patterns of all protein-coding mitochondrial genes were identified throughout the year along with mitochondrial copy numbers in young and old bulls' spermatozoa. For that, high-quality semen samples (n = 32) with more than 80% quality for the morphological parameters, from young (n = 4, aged 18-24 months old) and old (n = 4, aged 40-54 months old) Holstein bulls, were collected during the four seasons (n = 4 samples each animal/season). The DNA and RNA were isolated from sperm cells and subjected to the DNA copy number and expression analyses using qPCR. Furthermore, an in silico analysis using gene ontology online tools for the abundantly expressed genes was utilized. The data were statistically analyzed using Prism10 software. There was a significant reduction in the mitochondria copy number of young bulls' spermatozoa compared to their old counterparts during the summer (29 ± 3 vs. 51 ± 6, p < 0.001) and winter (27 ± 3 vs. 43 ± 7, p < 0.01) seasons. However, sperm-borne mitochondrial protein-coding genes were transcriptionally higher in young bulls throughout the year. Within the same group of bulls, unlike the old bulls, there was a significant (p < 0.05) induction in the transcription activity accompanied by a significant (p < 0.05) reduction in the mitochondrial copy numbers in the summer (29 ± 3) and winter (27 ± 3) compared to the spring (42 ± 9) and autumn (36 ± 5) seasons in young bulls. Additionally, the pathway enrichment of the top six expressed genes differed between age groups and seasons. In conclusion, under the same quality of semen, the early stages of age are associated with mitochondrial biogenesis and transcription activity dysregulation in a season-dependent manner.

Determinant genetic markers of semen quality in livestock

The reproductive efficiency of livestock is crucial for agricultural productivity and economic sustainability. One critical factor in successful fertilization and the viability of offspring is the quality of semen. Poor semen quality, especially in frozen-thawed semen used in artificial insemination (AI) have been shown to influence conception outcomes, resulting a negative impact on livestock production. Recent advancements in genetic research have identified specific markers linked to semen quality traits in various livestock species, such as cattle, sheep, goats, pigs, buffalo, and equines. These genetic markers are essential in screening males for breeding suitability, which in turn enhances selective breeding programs. Understanding these markers is crucial for improving reproductive performance and increasing productivity in livestock populations. This review offers a comprehensive overview of the genetic markers associated with semen quality in key livestock. It explores the underlying genetic mechanisms and their practical implications in animal breeding and management. The review underscores the importance of integrating genetic insights into breeding strategies to optimize reproductive efficiency and ensure the sustainable development of livestock industries.

Candidate genes and proteins regulating bull semen quality: a review

Poor semen profile reflected by suboptimum fertility statistics is a concern in bulls reared for breeding purpose. A critical review of research on candidate genes and proteins associated with semen quality traits will be useful to understand the progress of molecular marker development for bull semen quality traits. Here, we have tabulated and classified candidate genes and proteins associated with bull semen quality based on a literature survey. A total of 175 candidate genes are associated with semen quality traits in various breeds of cattle. Several studies using candidate gene approach have identified 26 genes carrying a total of 44 single nucleotide polymorphisms. Furthermore, nine genome-wide association studies (GWASes) have identified 150 candidate genes using bovine single nucleotide polymorphisms (SNP) chips. Three genes, namely membrane-associated ring-CH-type finger 1 (MARCH1), platelet-derived growth factor receptor beta, and phosphodiesterase type 1, were identified commonly in two GWASes, which, especially MARCH1, are required to explore their regulatory roles in bull semen quality in in-depth studies. With the advancement of high-throughput-omic technologies, more candidate genes associated with bull semen quality may be identified in the future. Therefore, the functional significance of candidate genes and proteins need to be delved further into future investigations to augment bull semen quality.

Beta-defensins as marker for male fertility: a comprehensive review†

Bovine male fertility in animals has a direct impact on the productivity of dairy herds. The epididymal sperm maturations involve extensive sperm surface modifications to gain the fertilizing ability, especially by absorptions of the plethora of biomolecules, including glycoprotein beta-defensins (BDs), enzymes, organic ions, protein, and phospholipids. Defensins are broad-range nonspecific antimicrobial peptides that exhibit strong relations with innate and adaptive immunity, but their roles in male fertility are relatively recently identified. In the course of evolution, BD genes give rise to different clusters with specific functions, especially reproductive functions, by undergoing duplications and nonsynonymous mutations. BD polymorphisms have been reported with milk compositions, disease resistance, and antimicrobial activities. However, in recent decades, the link of BD polymorphisms with fertility has emerged as an appealing improvement of reproductive performance such as sperm motility, membrane integrity, cervical mucus penetration, evading of uterus immunosurveillance, oviduct cell attachment, and egg recognition. The reproductive-specific glycosylated BD class-A BDs (CA-BDs) have shown age- and sex-specific expressions in male reproductive organs, signifying their physiological pleiotropism, especially in the sperm maturation and sperm transport in the female reproductive tract. By considering adult male reproductive organ-specific BD expressions, importance in sperm functionalities, and bioinformatic analysis, we have selected two bovine BBD126 and BBD129 genes as novel potential biomarkers of bovine male fertility. Despite the importance of BDs, however, genomic characterization of most BD genes across most livestock and nonmodel organisms remains predictive/incomplete. The current review discusses our understanding of BD pleiotropic functions, polymorphism, and genomic structural attributes concerning the fertilizability of the male gamete in dairy animals.

Analysis of amplification and association polymorphisms in the bovine beta-defensin 129 (BBD129) gene revealed its function in bull fertility

β-defensins are adsorbable on the sperm surface in the male reproductive tract (MRT) and enhance sperm functional characteristics. The beta-defensin 129 (DEFB129) antimicrobial peptide is involved in sperm maturation, motility, and fertilization. However, its role in bovine fertility has not been well investigated. This study examines the relationship between the bovine BBD129 gene and Bos indicus x Bos taurus bull fertility. The complete coding sequence of BBD129 mRNA was identified by RNA Ligase Mediated-Rapid Amplification of cDNA End (RLM-RACE) and Sanger sequencing methodologies. It consisted of 582 nucleotides (nts) including 5' untranslated region (UTR) (46nts) and 3'UTR (23nts). It conserves all beta-defensin-like features. The expression level of BBD129 was checked by RT-qPCR and maximal expression was detected in the corpus-epididymis region compared to other parts of MRT. Polymorphism in BBD129 was also confirmed by Sanger sequencing of 254 clones from 5 high fertile (HF) and 6 low fertile (LF) bulls at two positions, 169 T > G and 329A > G, which change the S57A and N110S in the protein sequence respectively. These two mutations give rise to four types of BBD129 haplotypes. The non-mutated TA-BBD129 (169 T/329A) haplotype was substantially more prevalent among high-fertile bulls (P < 0.005), while the double-site mutated GG-BBD129 (169 T > G/329A > G) haplotype was significantly more prevalent among low-fertile bulls (P < 0.005). The in silico analysis confirmed that the polymorphism in BBD129 results in changes in mRNA secondary structure, protein conformations, protein stability, extracellular-surface availability, post-translational modifications (O-glycosylation and phosphorylation), and affects antibacterial and immunomodulatory capabilities. In conclusion, the mRNA expression of BBD129 in the MRT indicates its region-specific dynamics in sperm maturation. BBD129 polymorphisms were identified as the deciding elements accountable for the changed proteins with impaired functionality, contributing to cross-bred bulls' poor fertility.

Sperm Transcriptome Analysis Accurately Reveals Male Fertility Potential in Livestock

Nowadays, selection of superior male candidates in livestock as a source of frozen semen based on sperm quality at the cellular level is not considered accurate enough for predicting the potential of male fertility. Sperm transcriptome analysis approaches, such as messenger RNA levels, have been shown to correlate with fertility rates. Using this technology in livestock growth has become the principal method, which can be widely applied to predict male fertility potential in the livestock industry through the analysis of the sperm transcriptome. It provides the gene expression to validate the function of sperm in spermatogenesis, fertilization, and embryo development, as the parameters of male fertility. This review proposes a transcriptomic analysis approach as a high-throughput method to predict the fertility potential of livestock more accurately in the future.

Novel module and hub genes of distinctive breast cancer associated fibroblasts identified by weighted gene co-expression network analysis

BACKGROUND

As abundant and heterogeneous stromal cells in tumor microenvironment, carcinoma-associated fibroblasts (CAFs) are critically involved in cancer progression.

METHODS

To identify co-expression module and hub genes of distinctive breast CAFs, weighted gene co-expression network analysis (WGCNA) was conducted based on the expression array results of CAFs from seven chemo-sensitive breast cancer (BC) patients and seven chemo-resistant ones before neo-adjuvant chemotherapy.

RESULTS

A total of 4916 genes were included in WGCNA, and 12 modules were determined. Module-trait assay showed that the blue module (cor = 0.97, P < 0.001) was associated with CAF-related chemo-resistance, which was enriched mainly as "inflammatory response", "interferon-gamma-mediated signaling" and "NIK/NF-kappaB signaling" pathways. Moreover, CXCL8, CXCL10, CXCL11, PLSCR1, RIPK2 and USP18 were found to be potentially associated with chemo-resistance related to CAFs and prognosis of BC.

CONCLUSIONS

Our current data offered valuable insights into the molecular mechanisms of distinctive breast CAFs, which was beneficial for revealing how chemo-resistance of BC was initiated.

Detection of bovine viral diarrhea virus genotype 1 in aerosol by a real time RT-PCR assay

BACKGROUND

As a pestivirus of the Flaviviridae family, bovine viral diarrhea virus (BVDV), has imposed a large burden on animal husbandry worldwide, and such virus can be transmitted mainly through direct contact with other infected animals and probably via aerosols. In the present study, we aimed to develop a real-time RT-PCR method for detection of BVDV-1 in aerosol samples.

METHODS

A pair of primers specific for highly conserved regions of the BVDV-1 5'-UTR was designed. The standard curve and sensitivity of the developed assay were assessed based on 10-fold serial dilutions of RNA molecular standard. The specificity of the assay was evaluated with other pestiviruses and infectious bovine viruses. The clinical performance was examined by testing 169 aerosol samples.

RESULTS

The results showed that a good linear relationship existed between the standard curve and the concentration of template. The lowest detection limit was 5.2 RNA molecules per reaction. This assay was specific for detection of BVDV-1, and no amplification was found for other pestiviruses such as classical swine fever virus (CSFV), border disease virus (BDV), and common infectious bovine viruses, including BVDV-2, infectious bovine rhinotracheitis virus (IBRV), bovine parainfluenza virus type 3 (BPIV-3), bovine respiratory syncytial virus (BRSV), bovine ephemeral fever virus (BEFV) and bovine coronavirus (BcoV). The assay was highly reproducible with low variation coefficient values (CVs) for intra-assay and inter-assay. A total of 169 aerosol samples collected from six dairy herds were tested using this method. The results showed that the positive detection rate of BVDV-1 was 17.2% (29/169), which was significantly higher compared with the conventional RT-PCR. Additionally, the positive samples (n = 29) detected by real-time RT-PCR were verified by BVDV RPA-LFD, and a concordance rate of 100% was obtained between them.

CONCLUSIONS

Taken together, we developed a real-time RT-PCR assay for quantitative analysis of BVDV-1 in aerosol samples, and our finding provided valuable insights into the risk on aerosol transmission of BVDV-1.

Gene Polymorphisms in Boar Spermatozoa and Their Associations with Post-Thaw Semen Quality

Genetic markers have been used to assess the freezability of semen. With the advancement in molecular genetic techniques, it is possible to assess the relationships between sperm functions and gene polymorphisms. In this study, variant calling analysis of RNA-Seq datasets was used to identify single nucleotide polymorphisms (SNPs) in boar spermatozoa and to explore the associations between SNPs and post-thaw semen quality. Assessment of post-thaw sperm quality characteristics showed that 21 boars were considered as having good semen freezability (GSF), while 19 boars were classified as having poor semen freezability (PSF). Variant calling demonstrated that most of the polymorphisms (67%) detected in boar spermatozoa were at the 3’-untranslated regions (3’-UTRs). Analysis of SNP abundance in various functional gene categories showed that gene ontology (GO) terms were related to response to stress, motility, metabolism, reproduction, and embryo development. Genomic DNA was isolated from sperm samples of 40 boars. Forty SNPs were selected and genotyped, and several SNPs were significantly associated with motility and membrane integrity of frozen-thawed (FT) spermatozoa. Polymorphism in SCLT1 gene was associated with significantly higher motility and plasma membrane integrity of FT spermatozoa from boars of the GSF group compared with those of the PSF group. Likewise, polymorphisms in MAP3K20, MS4A2, and ROBO1 genes were significantly associated with reduced cryo-induced lipid peroxidation and DNA damage of FT spermatozoa from boars of the GSF group. Candidate genes with significant SNP associations, including APPL1, PLBD1, FBXO16, EML5, RAB3C, OXSR1,PRICKLE1, and MAP3K20 genes, represent potential markers for post-thaw semen quality, and they might be relevant for future improvement in the selection procedure of boars for cryopreservation. The findings of this study provide evidence indicating that polymorphisms in genes expressed in spermatozoa could be considered as factors associated with post-thaw semen quality.

Cell Cycle Regulation in the Plant Response to Stress

As sessile organisms, plants face a variety of environmental challenges. Their reproduction and survival depend on their ability to adapt to these stressors, which include water, heat stress, high salinity, and pathogen infection. Failure to adapt to these stressors results in programmed cell death and decreased viability, as well as reduced productivity in the case of crop plants. The growth and development of plants are maintained by meiosis and mitosis as well as endoreduplication, during which DNA replicates without cytokinesis, leading to polyploidy. As in other eukaryotes, the cell cycle in plants consists of four stages (G1, S, G2, and M) with two major check points, namely, the G1/S check point and G2/M check point, that ensure normal cell division. Progression through these checkpoints involves the activity of cyclin-dependent kinases and their regulatory subunits known as cyclins. In order for plants to survive, cell cycle control must be balanced with adaption to dynamic environmental conditions. In this review, we summarize recent advances in our understanding of cell cycle regulation in plants, with a focus on the molecular interactions of cell cycle machinery in the context of stress tolerance.

Identification of modules and hub genes associated with platinum-based chemotherapy resistance and treatment response in ovarian cancer by weighted gene co-expression network analysis

High-grade serous ovarian carcinoma (HGSOC) is the most prevalent and malignant ovarian tumor.To identify co-expression modules and hub genes correlated with platinum-based chemotherapy resistant and sensitive HGSOC, we performed weighted gene co-expression network analysis (WGCNA) on microarray data of HGSOC with 12 resistant samples and 16 sensitive samples of GSE51373 dataset.A total of 5122 genes were included in WGCNA, and 16 modules were identified. Module-trait analysis identified that the module salmon (cor = 0.50), magenta (cor = 0.49), and black (cor = 0.45) were discovered associated with chemotherapy resistant, and the significance for these platinum-resistant modules were validated in the GSE63885 dataset. Given that the black module was validated to be the most related one, hub genes of this module, alcohol dehydrogenase 1B, cadherin 11, and vestigial like family member 3were revealed to be expressional related with platinum resistance, and could serve as prognostic markers for ovarian cancer.Our analysis might provide insight for molecular mechanisms of platinum-based chemotherapy resistance and treatment response in ovarian cancer.

Phylostratigraphic analysis of gene co-expression network reveals the evolution of functional modules for ovarian cancer

Ovarian cancer (OV) is an extremely lethal disease. However, the evolutionary machineries of OV are still largely unknown. Here, we used a method that combines phylostratigraphy information with gene co-expression networks to extensively study the evolutionary compositions of OV. The present co-expression network construction yielded 18,549 nodes and 114,985 edges based on 307 OV expression samples obtained from the Genome Data Analysis Centers database. A total of 20 modules were identified as OV related clusters. The human genome sequences were divided into 19 phylostrata (PS), the majority (67.45%) of OV genes was already present in the eukaryotic ancestor. There were two strong peaks of the emergence of OV genes screened by hypergeometric test: the evolution of the multicellular metazoan organisms (PS5 and PS6, P value = 0.002) and the emergence of bony fish (PS11 and PS12, P value = 0.009). Hence, the origin of OV is far earlier than its emergence. The integrated analysis of the topology of OV modules and the phylogenetic data revealed an evolutionary pattern of OV in human, namely, OV modules have arisen step by step during the evolution of the respective lineages. New genes have evolved and become locked into a pathway, where more and more biological pathways are fixed into OV modules by recruiting new genes during human evolution.

Targeted inhibition of histone deacetylase 6 in inflammatory diseases

Targeting epigenetic modification of gene expression represents a promising new approach under investigation for the treatment of inflammatory diseases. Accumulating evidence suggests that epigenetic mechanisms, such as histone modification, play a crucial role in a number of inflammatory diseases, including rheumatoid arthritis, asthma, and contact hypersensitivity. Consistent with this role, histone deacetylase (HDAC) inhibitors have shown efficacy in the treatment of inflammatory diseases. In particular, selective inhibitors of HDAC6, a cytoplasmic member of the HDAC family that contains two deacetylase domains, are under investigation as a potential treatment strategy for inflammatory diseases due to their ability to regulate inflammatory cells and cytokines. Here, we review recent findings highlighting the critical roles of HDAC6 in a variety of inflammatory diseases, and discuss the therapeutic potential of HDAC6 inhibitors in these settings.

Therapeutic targeting of cellular stress responses in cancer

Similar to bacteria, yeast, and other organisms that have evolved pathways to respond to environmental stresses, cancer cells develop mechanisms that increase genetic diversity to facilitate adaptation to a variety of stressful conditions, including hypoxia, nutrient deprivation, exposure to DNA-damaging agents, and immune responses. To survive, cancer cells trigger mechanisms that drive genomic instability and mutation, alter gene expression programs, and reprogram the metabolic pathways to evade growth inhibition signaling and immune surveillance. A deeper understanding of the molecular mechanisms that underlie the pathways used by cancer cells to overcome stresses will allow us to develop more efficacious strategies for cancer therapy. Herein, we overview several key stresses imposed on cancer cells, including oxidative, metabolic, mechanical, and genotoxic, and discuss the mechanisms that drive cancer cell responses. The therapeutic implications of these responses are also considered, as these factors pave the way for the targeting of stress adaption pathways in order to slow cancer progression and block resistance to therapy.

Survival Mechanisms to Selective Pressures and Implications

Organisms have evolved a spectrum of strategies that facilitate survival in the face of adverse environmental conditions. In order to make full use of the unfavorable resources of nature, human beings usually impose selective pressures to breed phenotypic traits that can survive in adverse environments. Animals are frequently under attack by biotic stress, such as bacterial and viral infections, while plants are more often subjected to abiotic stress, including high salinity, drought, and cold. In response to these diverse stresses, animals and plants initiate wide-ranging changes in gene expression by altering regulation of transcriptional and post-transcriptional activities. Recent studies have identified a number of key responsive components that promote survival of animals and plants in response to biotic and abiotic stresses. Importantly, with recent developments in genome-editing technology based on the CRISPR/Cas9 system, manipulation of genetic elements to generate stress-resistant animals and plants has become both feasible and cost-effective. Herein, we review important mechanisms that govern the response of organisms to biotic and abiotic stresses with the aim of applying our understanding to the agriculture and animal husbandry industries.

Molecular mechanisms underlying stress response and adaptation

Environmental stresses are ubiquitous and unavoidable to all living things. Organisms respond and adapt to stresses through defined regulatory mechanisms that drive changes in gene expression, organismal morphology, or physiology. Immune responses illustrate adaptation to bacterial and viral biotic stresses in animals. Dysregulation of the genotoxic stress response system is frequently associated with various types of human cancer. With respect to plants, especially halophytes, complicated systems have been developed to allow for plant growth in high salt environments. In addition, drought, waterlogging, and low temperatures represent other common plant stresses. In this review, we summarize representative examples of organismal response and adaptation to various stresses. We also discuss the molecular mechanisms underlying the above phenomena with a focus on the improvement of organismal tolerance to unfavorable environments.

Molecular Characterization and Comparison of Phospholipase C zeta (PLCZ1) Gene Between Swamp (Bubalus carabanensis) and Riverine (Bubalus bubalis) Buffaloes: Its Implications and Future Perspectives

Phospholipase C zeta, a novel sperm-specific protein which is widely known to induce oocyte activation following fertilization, had already been characterized in various mammalian species, but not in water buffaloes thus far. The present study was conducted to initially characterize and compare the sequences of PLCZ1 gene of swamp and riverine buffaloes. Semen samples were collected; total RNA was extracted and reverse-transcribed. PLCZ1 cDNA was then amplified, and submitted for sequencing. Buffalo PLCZ1 gene yielded a sequence of 1905 base pair nucleotides translated into 634 bp amino acids. In general, the buffalo PLCZ1 gene was found to have high sequence identity with cattle and other domestic species. Similarly, significant residues and motifs in PLCZ1 gene sequence are found conserved in water buffaloes. However, there are variations in sequences identified between types of water buffaloes that may play a role in species-specific differences in terms of gene and protein expression, physiological mechanisms, and biological functions. The molecular information on buffalo PLCZ1 gene is highly valuable in subsequent works such as correlation studies on the identified gene variations with semen quality and fertility, and the development of biomarkers for bull fertility.

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.

Splicing-related single nucleotide polymorphism of RAB, member of RAS oncogene family like 2B (RABL2B) jeopardises semen quality in Chinese Holstein bulls

RAB, member of RAS oncogene family like 2B (RABL2B) is a member of a poorly characterised clade of the RAS GTPase superfamily, which plays an essential role in male fertility, sperm intraflagellar transport and tail assembly. In the present study, we identified a novel RABL2B splice variant in bovine testis and spermatozoa. This splice variant, designated RABL2B-TV, is characterised by exon 2 skipping. Moreover, a single nucleotide polymorphism (SNP), namely c.125G>A, was found within the exonic splicing enhancer (ESE) motif, indicating that the SNP caused the production of the RABL2B-TV aberrant splice variant. This was demonstrated by constructing a pSPL3 exon capturing vector with different genotypes and transfecting these vectors into murine Leydig tumour cell line (MLTC-1) cells. Expression of the RABL2B-TV transcript was lower in semen from high- versus low-performance bulls. Association analysis showed that sperm deformity rate was significantly lower in Chinese Holstein bulls with the GG or GA genotype than in bulls with the AA genotype (P < 0.05). In addition, initial sperm motility was significantly higher in individuals with the GG or GA genotype than in individuals with the AA genotype (P < 0.05). The findings of the present study suggest that the difference in semen quality in bulls with different RABL2B genotypes is generated via an alternative splicing mechanism caused by a functional SNP within the ESE motif.

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.

A g.-1256 A>C in the promoter region of CAPN1 is associated with semen quality traits in Chinese Holstein bulls

The micromolar calcium-activated neutral protease gene (CAPN1) is a physiological candidate gene for sperm motility. However, the molecular mechanisms involved in regulating the expression of the CAPN1 gene in bulls remain unknown. In this study, we investigated the expression pattern of CAPN1 in testis, epididymis, and sperm at the RNA and protein levels by qRT-PCR, western blot, immunohistochemistry, and immunofluorescence assay. Results revealed that the expression of CAPN1 levels was higher in the sperm head compared with that in other tissues. Moreover, we identified a novel single-nucleotide polymorphism (g.-1256 A>C, ss 1917715340) in the noncanonical core promoter of the CAPN1 gene between base g.-1306 and g.-1012. Additionally, we observed greater sperm motility in bulls with the genotype CC than in those with the genotype AA (P<0.01), indicating that different genotypes were associated with the bovine semen trait. Furthermore, a higher fluorescence intensity of the C allele than that of the A allele at g. -1256 A>C was revealed by transient transfection in MLTC-1 cells and luciferase report assay. Finally, CAPN1 was highly expressed in the spermatozoa with the CC genotype compared with that with the AA genotype by qRT-PCR. This study is the first report on genetic variant g.-1256 A>C in the promoter region of CAPN1 gene association with the semen quality of Chinese Holstein bulls by influencing its expression. g.-1256 A>C can be a functional molecular marker in cattle breeding.

The g.-165 T>C Rather than Methylation Is Associated with Semen Motility in Chinese Holstein Bulls by Regulating the Transcriptional Activity of the HIBADH Gene

The 3-hydroxyisobutyrate dehydrogenase (HIBADH) is regarded as a human sperm-motility marker. However, the molecular mechanisms involved in the regulation of expression of the HIBADH gene in bulls remain largely unknown. HIBADH was detected in the testis, epididymis, and sperm via reverse transcription polymerase chain reaction and Western blot analysis. It is also expressed in the seminiferous epithelium, spermatids, and the entire epididymis, as detected by immunohistochemistry. Furthermore, HIBADH was expressed in the neck-piece and mid-piece of bull spermatids, as shown in the immunofluorescence assay. Using serially truncated bovine HIBADH promoters and luciferase constructs, we discovered an 878 bp (-703 bp to +175 bp) fragment that constitutes the core promoter region. One SNP g.-165 T>C of HIBADH was identified and genotyped in 307 Chinese Holstein bulls. Correlation analysis revealed that bulls with the TT genotype had higher initial sperm motility than those with the CC genotype (P < 0.05). Furthermore, the T- or C-containing loci (designated as pGL3-T and pGL3-C) were transiently transfected into MLTC-1 to test the effect of SNP on HIBADH expression. The luciferase reporter assay showed that the pGL3-T genotype exhibited 58% higher transcriptional activity than the pGL3-C genotype (P < 0.05). The bisulfite sequencing analysis revealed that the methylation pattern of the core promoter presented hypomethylation in the ejaculated semen in high-motility and low-motility bulls. The results demonstrated for the first time that the g.-165 T>C rather than methylation in the 5'-flanking region could affect the bovine sperm motility through the regulation of HIBADH gene transcriptional activity.

Genome-wide association study for sperm concentration in Holstein-Friesian bulls

The aim of the study was to screen the entire bull genome to identify markers and candidate genes underlying sperm concentration. The analysed data set originates from a population of 877 Polish Holstein-Friesian bulls. Based on sperm concentration value, two extreme groups of bulls were created: Low (L, n = 126) and High (H, n = 140). Each bull was genotyped using the Illumina BovineSNP50 BeadChip. Genome-wide association analysis was performed with the use of GoldenHelix SVS7 software. An additive model with a Cohran-Armitage test, Correlation/Trend adjusted by a Bonferroni test, was used to estimate the effect of SNP marker for sperm concentration. Thirteen markers reached genome-wide significance. The most significant SNPs were located on chromosome 3 (rs109154964 and rs108965556), 14 (rs41621145) and 18 (rs41615539), in the close vicinity of protein arginine methyltransferase 6 (PRMT6), Sel1 repeat containing 1 (SELRC1), triple QxxK/R motif containing (TRIQK) and zinc finger homeobox 3 (ZFHX3) genes, respectively. For three other candidate genes located close to significant markers (within a distance of ca 1 Mb), namely histone deacetylase 9 (HDAC9), an inhibitor of DNA binding 2 (ID2) and glutathione S-transferase theta 1 (GSTT1), their potential role in the production of male germ cells was confirmed in earlier studies. Six additional candidate genes (Vav3, GSTM1, CDK5, NOS3, PDP1 and GAL3ST1) were suspected of being significantly associated with sperm concentration or semen biochemistry. Our results indicate the genetic complexity of sperm concentration but also open the possibility for finding causal polymorphism useful in marker-assisted selection.

Association between an alternative promoter polymorphism and sperm deformity rate is due to modulation of the expression of KATNAL1 transcripts in Chinese Holstein bulls

Katanin p60 subunit A-like 1 (KATNAL1) is an ATPase that regulates Sertoli cell microtubule dynamics and sperm retention. We evaluated one novel splice variant and characterized the promoter and a functional single nucleotide polymorphism (SNP) of the bovine KATNAL1 gene to explore its expression pattern, possible regulatory mechanism and relationship with semen traits in Chinese Holstein bulls. A novel splice variant, KATNAL1 transcript variant 2 (KATNAL1-TV2) of the retained 68 bp in intron 2, was identified by RT-PCR and compared with KATNAL1 transcript variant 1 (KATNAL1-TV1, NM 001192918.1) in various tissues. Bioinformatics analyses predicted that KATNAL1 transcription was regulated by two promoters: P1 in KATNAL1-TV1 and P2 in KATNAL1-TV2. Results of qRT-PCR revealed that KATNAL1-TV1 had higher expression than did KATNAL1-TV2 in testes of adult bulls (P < 0.05). Promoter luciferase activity analysis suggested that the core sequences of P1 and P2 were mapped to the region of c.-575˜c.-180 and c.163-40˜c.333+59 respectively. One novel SNP (c.163-210T>C, ss836312085) located in intron 1 was found using sequence alignment. The SNP in P2 resulted in the presence of the DeltaE binding site, improving its basal promoter activity (P < 0.05); and we observed a greater sperm deformity rate in bulls with the genotype CC than in those with the genotype TT (P < 0.05), which indicated that different genotypes were associated with the bovine semen traits. Bioinformatics analysis of the KATNAL1 protein sequence predicted that the loss of the MIT domain in the KATNAL1-TV2 transcript resulted in protein dysfunction. These findings help us to understand that a functional SNP in P2 and subsequent triggering of expression diversity of KATNAL1 transcripts are likely to play an important role with regard to semen traits in bull breeding programs.

Genome-wide association study for poor sperm motility in Holstein-Friesian bulls

The aim of the study was to screen the whole bull genome to identify markers and candidate genes underlying poor sperm motility. The analyzed data set originates from the Polish Holstein-Friesian bull population and consists of 41 Case and 279 Control bulls (selected from 1581 bulls). The most distinguishing trait of case group was very poor sperm motility (average 25.61%) when compared to control samples (average 72.95%). Each bull was genotyped using the Illumina BovineSNP50 BeadChip. Genome-wide association analysis was performed with the use of GoldenHelix SVS7 software. An additive model with a Cohran-Armitage test, Correlation/Trend adjusted by Bonferroni test were used to estimate the effect of Single Nucleotide Polymorphism (SNP) marker for poor sperm motility. Markers (n=34) reached genome-wide significance. The most significant SNP were located on chromosome 24 (rs110876480), 5 (rs110827324 and rs29011704), and 1 (rs110596818), in the close vicinity of melanocortin 4 receptor (MC4R), PDZ domain containing ring finger 4 (PDZRN4) and ethanolamine kinase 1 (ETNK1), olfactory receptor 5K3-like (LOC785875) genes, respectively. For five other candidate genes located close to significant markers (in distance of ca. 1 Mb), namely alkaline phosphatase, liver/bone/kidney (ALPL), tripartite motif containing 36 (TRIM36), 3-hydroxyisobutyrate dehygrogenase (HIBADH), kelch-like 1 (KLHL1), protein kinase C, beta (PRKCB), their potential role in sperm motility was confirmed in the earlier studies. Five additional candidate genes, cystic fibrosis transmembrane conductance regulator (CFTR), insulin-like growth factor 1 receptor (IGF1R), steroid-5-alpha-reductase, alpha polypeptide 2 (SRD5A2), cation channel, sperm associated 1 (CATSPER1) calpain 1 (mu/I) large subunit (CAPN1) were suggested to be significantly associated with sperm motility or semen biochemistry. Results of the present study indicate there is a genetic complexity of poor sperm motility but also indicate there might be a causal polymorphism useful in marker-assisted selection. Identifying genomic regions associated with poor sperm motility may be very important for early recognition of a young sire as unsuitable for effective semen production in artificial insemination centers.

A novel subclass of bovine β-defensins links reproduction and immunology

β-defensins are effector molecules of the innate immune system, found in many diverse species. Their presence in invertebrates as well as vertebrates suggests highly conserved functional roles. Most β-defensins are believed to act as antimicrobial agents at epithelial surfaces, although additional functions have also been described, including immune regulatory activity, wound repair and a role in coat-colour determination. High expression of β-defensins have been found in testis and epididymidal epithelium as well as in the seminal fluid of humans, macaque, rat, mouse and cow. Human and macaque β-defensins have recently been shown to affect sperm motility while a mutation in β-defensin 126 is associated with reduced fertility in men. Genetic variation in bovine defensin genes may explain the increased incidence of low fertility in cattle. Here, we present a summary of the known functions of β-defensins as well as their emerging role in reproduction and their potential to improve fertility in cattle.

PLCz functional haplotypes modulating promoter transcriptional activity are associated with semen quality traits in Chinese Holstein bulls

The sperm-specific phospholipase C zeta (PLCz) is a candidate sperm-borne oocyte-activating factor that triggers a characteristic series of physiological stimuli via cytoplasmic Ca²⁺ oscillations during fertilization. The molecular mechanisms involved in the regulation of PLCz gene expression remain largely unknown. To explore the genetic variations in the 5′-flanking region of the PLCz gene and their common haplotypes in Chinese Holstein bulls, as well as to determine whether these variations affect bovine semen quality traits and transcriptional activity, DNA samples were collected from Chinese Holstein bulls and sequenced for the identification of genetic variants in the 5′-flanking region of PLCz. Two genetic variants were identified, and their haplotypic profiles were constructed. The two novel genetic variations (g. −456 G>A and g. +65 T>C) were genotyped in 424 normal Chinese Holstein bulls. Bioinformatics analysis revealed that both loci are in transcription factor binding sites of the core promoter region. The association studies revealed that the two genetic variations and their haplotype combinations significantly affected semen quality traits. Using serially truncated constructs of the bovine PLCz promoters and the luciferase reporter, we found that a 726 bp (−641 nt to +112 nt) fragment constitutes the core promoter region. Furthermore, four haplotypes, H1H1 (GTGT), H2H2 (GCGC), H3H3 (ATAT), and H4H4 (ACAC), were significantly associated with semen quality traits and successfully transfected into MLTC-1 cell lines. The luciferase reporter assay showed that the different haplotypes exhibited distinct promoter activities. Maximal promoter activity was demonstrated by the H2H2 haplotypes, as compared with the other haplotypes. To the best of our knowledge, this study is the first report on genetic variants and their respective haplotypes in the 5′-flanking region of PLCz gene that can influence the semen quality of Chinese Holstein bulls as well as contribute to the transcriptional activity of the PLCz promoter.