Two-stage genome-wide association study identifies integrin beta 5 as having potential role in bull fertility

Genome-Wide Association Studies for Body Conformation Traits in Korean Holstein Population

The objective of this study was to identify quantitative trait loci (QTL) and nearby candidate genes that influence body conformation traits. Phenotypic data for 24 body conformation traits were collected from a population of 2329 Korean Holstein cattle, and all animals were genotyped using the 50 K Illumina bovine SNP chip. A total of 24 genome-wide significant SNPs associated with 24 body conformation traits were identified by genome-wide association analysis. The selection of the most promising candidate genes was based on gene ontology (GO) terms and the previously identified functions that influence various body conformation traits as determined in our study. These genes include KCNA1, RYBP, PTH1R, TMIE, and GNAI3 for body traits; ANGPT1 for rump traits; MALRD1, INHBA, and HOXA13 for feet and leg traits; and CDK1, RHOBTB1, and SLC17A1 for udder traits, respectively. These findings contribute to our understanding of the genetic basis of body conformation traits in this population and pave the way for future breeding strategies aimed at enhancing desirable traits in dairy cattle.

Genome-wide association study reveals candidate markers related to field fertility and semen quality traits in Holstein-Friesian bulls

In vitro assessment of bull semen quality is routinely used in bull semen processing centres in order to ensure that semen destined to be used in the field has passed minimum standards. Despite these stringent quality control checks, individual bulls that pass the quality control checks can still vary in field fertility by up to 25%. A genome-wide association study was undertaken to determine genetic markers associated with prefreeze and post-thaw bull sperm quality traits as well as field fertility. Genome-wide association analysis was performed using a single nucleotide polymorphism (SNP) regression mixed linear model in WOMBAT. Genes within a 250 Kb span of a suggestive (P ≤ 1 × 10^-5) SNP were considered as candidate genes. One SNP was associated with adjusted pregnancy rate, and 21 SNPs were associated across the seven semen quality traits (P ≤ 1 × 10^-5). Functional candidate genes include SIPA1L2 which was associated with adjusted pregnancy rate. This encodes a Rap GTPase-activating protein involved in Rap1 signalling pathway and was previously found to play a role in the process of sperm differentiation. Gene ontology (GO) analysis also identified significantly enriched biological processes involved protein tyrosine kinase activity including genes such as DYRK1A, TEC and TXK that were associated with sperm motility prior to freezing. Another candidate gene associated with post-thaw sperm motility was FHDC1 which coordinates actin filament and microtubule dynamics. The induced 11 GO terms in the ejaculates rejected after freezing trait were related to ATPase, phosphatase and hydrolase activity. These results reveal novel specific genomic regions and candidate genes associated with economically important phenotypes such as field fertility and semen quality traits.

Identifying pleiotropic variants and candidate genes for fertility and reproduction traits in Holstein cattle via association studies based on imputed whole-genome sequence genotypes

Background

Genetic progress for fertility and reproduction traits in dairy cattle has been limited due to the low heritability of most indicator traits. Moreover, most of the quantitative trait loci (QTL) and candidate genes associated with these traits remain unknown. In this study, we used 5.6 million imputed DNA sequence variants (single nucleotide polymorphisms, SNPs) for genome-wide association studies (GWAS) of 18 fertility and reproduction traits in Holstein cattle. Aiming to identify pleiotropic variants and increase detection power, multiple-trait analyses were performed using a method to efficiently combine the estimated SNP effects of single-trait GWAS based on a chi-square statistic.

Results

There were 87, 72, and 84 significant SNPs identified for heifer, cow, and sire traits, respectively, which showed a wide and distinct distribution across the genome, suggesting that they have relatively distinct polygenic nature. The biological functions of immune response and fatty acid metabolism were significantly enriched for the 184 and 124 positional candidate genes identified for heifer and cow traits, respectively. No known biological function was significantly enriched for the 147 positional candidate genes found for sire traits. The most important chromosomes that had three or more significant QTL identified are BTA22 and BTA23 for heifer traits, BTA8 and BTA17 for cow traits, and BTA4, BTA7, BTA17, BTA22, BTA25, and BTA28 for sire traits. Several novel and biologically important positional candidate genes were strongly suggested for heifer (SOD2, WTAP, DLEC1, PFKFB4, TRIM27, HECW1, DNAH17, and ADAM3A), cow (ANXA1, PCSK5, SPESP1, and JMJD1C), and sire (ELMO1, CFAP70, SOX30, DGCR8, SEPTIN14, PAPOLB, JMJD1C, and NELL2) traits.

Conclusions

These findings contribute to better understand the underlying biological mechanisms of fertility and reproduction traits measured in heifers, cows, and sires, which may contribute to improve genomic evaluation for these traits in dairy cattle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08555-z.

Identification of Candidate Genes Regulating Carcass Depth and Hind Leg Circumference in Simmental Beef Cattle Using Illumina Bovine Beadchip and Next-Generation Sequencing Analyses

Genome-wide association studies are a robust means of identifying candidate genes that regulate economically important traits in farm animals. The aim of this study is to identify single-nucleotide polymorphisms (SNPs) and candidate genes potentially related to carcass depth and hind leg circumference in Simmental beef cattle. We performed Illumina Bovine HD Beadchip (~670 k SNPs) and next-generation sequencing (~12 million imputed SNPs) analyses of data from 1252 beef cattle, to which we applied a linear mixed model. Using a statistical threshold (p = 0.05/number of SNPs identified) and adopting a false discovery rate (FDR), we identified many putative SNPs on different bovine chromosomes. We identified 12 candidate genes potentially annotated with the markers identified, including CDKAL1 and E2F3, related to myogenesis and skeletal muscle development. The identification of such genes in Simmental beef cattle will help breeders to understand and improve related traits, such as meat yield.

Role of Integrins in Sperm Activation and Fertilization

In mammals, integrins are heterodimeric transmembrane glycoproteins that represent a large group of cell adhesion receptors involved in cell-cell, cell-extracellular matrix, and cell-pathogen interactions. Integrin receptors are an important part of signalization pathways and have an ability to transmit signals into and out of cells and participate in cell activation. In addition to somatic cells, integrins have also been detected on germ cells and are known to play a crucial role in complex gamete-specific physiological events, resulting in sperm-oocyte fusion. The main aim of this review is to summarize the current knowledge on integrins in reproduction and deliver novel perspectives and graphical interpretations presenting integrin subunits localization and their dynamic relocation during sperm maturation in comparison to the oocyte. A significant part of this review is devoted to discussing the existing view of the role of integrins during sperm migration through the female reproductive tract; oviductal reservoir formation; sperm maturation processes ensuing capacitation and the acrosome reaction, and their direct and indirect involvement in gamete membrane adhesion and fusion leading to fertilization.

Comprehensive Analysis of miRNAs and Target mRNAs between Immature and Mature Testis Tissue in Chinese Red Steppes Cattle

Simple Summary

MicroRNAs are small molecules that can regulate the relative abundance of their target genes by binding to the 3′ untranslated region of the target genes at the post-transcriptional level to affect various biological processes, such as biosynthesis, fat metabolism and proliferation, apoptosis, and cell differentiation. Fertility is one of the most important economic traits in livestock production. Bulls require the continuous production of high-quality spermatozoa in abundance. The quality of semen is an exceptionally important factor affecting the fertilization rate of the dairy cow and is also associated with the increasing conception rate in the process of artificial insemination. Therefore, accurately predicting fertility potential for a semen sample from donor bull for artificial insemination is crucial for consistently high reproductive efficiency. The present study performed a genome-wide sequencing analysis of microRNAs and mRNAs between immature and mature testes of Chinese Red Steppes. These results provide novel candidate microRNAs and functional genes related to bull reproduction traits and the networks between microRNAs and target genes, which will provide a useful genetic mechanism and epigenetic information for marker-assisted selection of bulls with excellent sperm quality in the future.

Abstract

This study aims to screen potential regulators and regulate fecundity networks between microRNAs (miRNAs) and target genes. The bovine testes of immature and mature Chinese Red Steppes were performed by genome-wide analysis of mRNAs and miRNAs. Compared with testicular tissues of newborns, 6051 upregulated genes and 7104 downregulated genes in adult cattle were identified as differentially expressed genes (DEGs). The DEGs were significantly enriched in 808 GO terms (p < 0.05) including male gonad development, male genitalia development, spermatogenesis, and sperm motility. Moreover, DEGs were also significantly enriched in 105 KEGG pathways (p < 0.05), including cGMP-PKG signaling pathway and calcium signaling pathway. To explore the expression of miRNA-regulated gene expression, 896 differentially expressed target genes negatively regulated with the expression levels of 31 differentially expressed miRNAs (DERs) were predicted and analyzed, and a network-integrated analysis was constructed. Furthermore, real-time PCR was performed to verify the expression levels of DEGs and DERs. Our results identified novel candidate DEGs and DERs correlated with male reproduction and intricate regulating networks between miRNAs and genes, which will be valuable for future genetic and epigenetic studies of sperm development and maturity, as well as providing valuable insights into the molecular mechanisms of male fertility and spermatogenesis in cattle.

Sperm Functional Genome Associated With Bull Fertility

Bull fertility is an important economic trait in sustainable cattle production, as infertile or subfertile bulls give rise to large economic losses. Current methods to assess bull fertility are tedious and not totally accurate. The massive collection of functional data analyses, including genomics, proteomics, metabolomics, transcriptomics, and epigenomics, helps researchers generate extensive knowledge to better understand the unraveling physiological mechanisms underlying subpar male fertility. This review focuses on the sperm phenomes of the functional genome and epigenome that are associated with bull fertility. Findings from multiple sources were integrated to generate new knowledge that is transferable to applied andrology. Diverse methods encompassing analyses of molecular and cellular dynamics in the fertility-associated molecules and conventional sperm parameters can be considered an effective approach to determine bull fertility for efficient and sustainable cattle production. In addition to gene expression information, we also provide methodological information, which is important for the rigor and reliability of the studies. Fertility is a complex trait influenced by several factors and has low heritability, although heritability of scrotal circumference is high and that it is a known fertility maker. There is a need for new knowledge on the expression levels and functions of sperm RNA, proteins, and metabolites. The new knowledge can shed light on additional fertility markers that can be used in combination with scrotal circumference to predict the fertility of breeding bulls. This review provides a comprehensive review of sperm functional characteristics or phenotypes associated with bull fertility.

GALLO: An R package for genomic annotation and integration of multiple data sources in livestock for positional candidate loci

BACKGROUND

The development of high-throughput sequencing and genotyping methodologies has enabled the identification of thousands of genomic regions associated with several complex traits. The integration of multiple sources of biological information is a crucial step required to better understand patterns regulating the development of these traits.

FINDINGS

Genomic Annotation in Livestock for positional candidate LOci (GALLO) is an R package developed for the accurate annotation of genes and quantitative trait loci (QTLs) located in regions identified in common genomic analyses performed in livestock, such as genome-wide association studies and transcriptomics using RNA sequencing. Moreover, GALLO allows the graphical visualization of gene and QTL annotation results, data comparison among different grouping factors (e.g., methods, breeds, tissues, statistical models, studies), and QTL enrichment in different livestock species such as cattle, pigs, sheep, and chickens.

CONCLUSIONS

Consequently, GALLO is a useful package for annotation, identification of hidden patterns across datasets, and data mining previously reported associations, as well as the efficient examination of the genetic architecture of complex traits in livestock.

Bovine seminal plasma osteopontin: Structural modelling, recombinant expression and its relationship with semen quality

Osteopontin (OPN) is a multifunctional phosphoprotein that has been linked to fertility in bulls. However, the exact mechanism by which OPN contributes to fertilisation is yet unknown. The biotechnological use of OPN in bovine reproduction is promising but some gaps remain unfilled. The present work aimed: (a) to verify whether the seminal plasma OPN is associated with seminal traits and a standard breeding soundness exam; (b) to predict OPN interactions with integrins, CD44 and glycosaminoglycans through molecular docking; and (c) to develop a protocol for recombinant expression of OPN from vesicular gland cDNA. Ejaculates from top ranked bulls had higher amounts of seminal plasma OPN in comparison with bulls classified as questionable (p < .01). The structural modelling and molecular docking predictions indicated that bovine OPN binds to heparin disaccharide, hyaluronic acid and hyaluronan. In addition, docking studies described the binding complexes of OPN with CD44 and the integrin heterodimers α5β1 and αVβ3. Finally, expression of rOPN-6His was successfully obtained after 3 hr of induction with 0.5 mM IPTG at 37°C and a denaturing purification protocol resulted in efficiently purified recombinant OPN. The present results contribute to the development of biotechnological uses of OPN as a biomarker in bovine reproduction.

Marker-assisted selection vis-à-vis bull fertility: coming full circle-a review

Bull fertility is considered an indispensable trait, as far as farm economics is concerned since it is the successful conception in a cow that provides calf crop, along with the ensuing lactation. This ensures sustainability of a dairy farm. Traditionally, bull fertility did not receive much attention by the farm managers and breeding animals were solely evaluated based on phenotypic predictors, namely, sire conception rate and seminal parameters in bull. With the advent of the molecular era in animal breeding, attempts were made to unravel the genetic complexity of bull fertility by the identification of genetic markers related to the trait. Marker-Assisted Selection (MAS) is a methodology that aims at utilizing the genetic information at markers and selecting improved populations for important traits. Traditionally, MAS was pursued using a candidate gene approach for identifying markers related to genes that are already known to have a physiological function related to the trait but this approach had certain shortcomings like stringent criteria for significance testing. Now, with the availability of genome-wide data, the number of markers identified and variance explained in relation to bull fertility has gone up. So, this presents a unique opportunity to revisit MAS by selection based on the information of a large number of genome-wide markers and thus, improving the accuracy of selection.

Adhesion molecules in gamete transport, fertilization, early embryonic development, and implantation-role in establishing a pregnancy in cattle: A review

Cell-cell adhesion molecules have critically important roles in the early events of reproduction including gamete transport, sperm-oocyte interaction, embryonic development, and implantation. Major adhesion molecules involved in reproduction include cadherins, integrins, and disintegrin and metalloprotease domain-containing (ADAM) proteins. ADAMs on the surface of sperm adhere to integrins on the oocyte in the initial stages of sperm-oocyte interaction and fusion. Cadherins act in early embryos to organize the inner cell mass and trophectoderm. The trophoblast and uterine endometrial epithelium variously express cadherins, integrins, trophinin, and selectin, which achieve apposition and attachment between the elongating conceptus and uterine epithelium before implantation. An overview of the major cell-cell adhesion molecules is presented and this is followed by examples of how adhesion molecules help shape early reproductive events. The argument is made that a deeper understanding of adhesion molecules and reproduction will inform new strategies that improve embryo survival and increase the efficiency of natural mating and assisted breeding in cattle.

Transcriptome profile and association study revealed STAT3 gene as a potential quality marker of bovine gametes

The present study was aimed to investigate differences in molecular signatures in oocytes derived from Holstein-Friesian heifers with different genetic merit for fertility, euthanized during day 0 or day 12 of the estrous cycle. Moreover, association between single nucleotide polymorphisms (SNPs) of ODC1 and STAT3 genes and bull fertility traits was investigated. The gene expression patterns were analyzed using cDNA array and validated with quantitative real-time polymerase chain reaction (PCR). The result revealed that several genes have shown not only to be regulated by fertility merit but also by the day of oocyte recovery during the estrous cycle. The STAT3 gene was found to be upregulated in oocytes recovered from animals with high fertility merit at both day 0 and day 12. Some other genes like PTTG1, ODC1 and TUBA1C were downregulated at day 0 and upregulated at day 12 in high, compared with low, fertility merit recovered oocytes. In contrast, the transcript abundance of TPM3 was upregulated at day 0 and downregulated at day 12 in high, compared with low, fertility merit recovered oocytes. In addition, ODC1 and STAT3 were found to be associated (P < 0.05) with sperm quality traits as well as flow cytometry parameters. Therefore, the expression of several candidate genes including ODC1 and STAT3 was related to the genetic merit of the cow. In addition polymorphisms in these two genes were found to be associated with bull semen quality.

A targeted genotyping approach to enhance the identification of variants for lactation persistency in dairy cows

Lactation persistency (LP), defined as the ability of a cow to maintain milk production at a high level after milk peak, is an important phenotype for the dairy industry. In this study, we used a targeted genotyping approach to scan for potentially functional single nucleotide polymorphisms (SNPs) within 57 potential candidate genes derived from our previous genome wide association study on LP and from the literature. A total of 175,490 SNPs were annotated within 10-kb flanking regions of the selected candidate genes. After applying several filtering steps, a total of 105 SNPs were retained for genotyping using target genotyping arrays. SNP association analyses were performed in 1,231 Holstein cows with 69 polymorphic SNPs using the univariate liner mixed model with polygenic effects using DMU package. Six SNPs including rs43770847, rs208794152, and rs208332214 in ADRM1; rs209443540 in C5orf34; rs378943586 in DDX11; and rs385640152 in GHR were suggestively significantly associated with LP based on additive effects and associations with 4 of them (rs43770847, rs208794152, rs208332214, and rs209443540) were based on dominance effects at P < 0.05. However, none of the associations remained significant at false discovery rate adjusted P (FDR) < 0.05. The additive variances explained by each suggestively significantly associated SNP ranged from 0.15% (rs43770847 in ADRM1) to 5.69% (rs209443540 in C5orf34), suggesting that these SNPs might be used in genetic selection for enhanced LP. The percentage of phenotypic variance explained by dominance effect ranged from 0.24% to 1.35% which suggests that genetic selection for enhanced LP might be more efficient by inclusion of dominance effects. Overall, this study identified several potentially functional variants that might be useful for selection programs for higher LP. Finally, a combination of identification of potentially functional variants followed by targeted genotyping and association analysis is a cost-effective approach for increasing the power of genetic association studies.

Selection for bull fertility: a review

Fertility is a critically important factor in cattle production because it directly relates to the ability to produce the offspring necessary to offset costs in production systems. Female fertility has received much attention and has been enhanced through assisted reproductive technologies, as well as genetic selection; however, improving bull fertility has been largely ignored. Improvements in bull reproductive performance are necessary to optimize the efficiency of cattle production. Selection and management to improve bull fertility not only have the potential to increase conception rates but also have the capacity to improve other economically relevant production traits. Bull fertility has reportedly been genetically correlated with traits such as average daily gain, heifer pregnancy, and calving interval. Published studies show that bull fertility traits are low to moderately heritable, indicating that improvements in bull fertility can be realized through selection. Although female fertility has continued to progress according to increasing conception rates, the reported correlation between male and female fertility is low, indicating that male fertility cannot be improved by selection for female fertility. Correlations between several bull fertility traits, such as concentration, number of spermatozoa, motility, and number of spermatozoa abnormalities, vary among studies. Using male fertility traits in selection indices would provide producers with more advanced selection tools. The objective of this review was to discuss current beef bull fertility measurements and to discuss the future of genetic evaluation of beef bull fertility and potential genetic improvement strategies.

Identification of muscle-specific candidate genes in Simmental beef cattle using imputed next generation sequencing

Genome-wide association studies (GWAS) have commonly been used to identify candidate genes that control economically important traits in livestock. Our objective was to detect potential candidate genes associated mainly with muscle development traits related to dimension of hindquarter in cattle. A next generation sequencing (NGS) dataset to imputed to 12 million single nucleotide polymorphisms (SNPs) (from 1252 Simmental beef cattle) were used to search for genes affecting hindquarter traits using a linear, mixed model approach. We also used haplotype and linkage disequilibrium blocks to further support our identifications. We identified 202 significant SNPs in the bovine BTA4 chromosome region associated with width of hind leg, based on a stringent statistical threshold (p = 0.05/ effective number of SNPs identified). After exploring the region around these SNPs, we found candidate genes that were potentially related to the associated markers. More importantly, we identified a region of approximately 280 Kb on the BTA4 chromosome that harbored several muscle specific candidate genes, genes to be in a potential region for muscle development. However, we also found candidate gene SLC13A1 on BTA4, which seems to be associated with bone disorders (such as chondrodysplasia) in Simmental beef cattle.

A Nutrigenomics Approach Using RNA Sequencing Technology to Study Nutrient-Gene Interactions in Agricultural Animals

Thorough understanding of animal gene expression driven by dietary nutrients can be regarded as a bottom line of advanced animal nutrition research. Nutrigenomics (including transcriptomics) studies the effects of dietary nutrients on cellular gene expression and, ultimately, phenotypic changes in living organisms. Transcriptomics can be applied to investigate animal tissue transcriptomes at a defined nutritional state, which can provide a holistic view of intracellular RNA expression. As a novel transcriptomics approach, RNA sequencing (RNA-Seq) technology can monitor all gene expressions simultaneously in response to dietary intervention. The principle and history of RNA-Seq are briefly reviewed, and its 3 principal steps are described in this article. Application of RNA-Seq in different areas of animal nutrition research is summarized. Lastly, the application of RNA-Seq in swine science and nutrition is also reviewed. In short, RNA-Seq holds significant potential to be employed for better understanding the nutrient-gene interactions in agricultural animals.

Expression dynamics of Integrin Subunit Beta 5 in bovine gametes and embryos imply functions in male fertility and early embryonic development

Integrins have been shown to act as signalling receptors, and they primarily recognise extracellular matrix ligands on the oocyte surface. However, their possible roles in oocyte activation and embryo development are not clearly understood. The objectives of this study were to evaluate expression of Integrin Subunit Beta 5 (ITGβ5) in bovine sperm, oocytes, and early embryos and to ascertain the evolutionary conservation of ITGβ5. To accomplish these objectives, we used western blotting to study expression levels of ITGβ5 protein in sperm and RT-qPCR to determine expression levels of ITGβ5 transcripts in oocytes and embryos. We have also used bioinformatic analysis to determine the evolutionary conservation of the ITGβ5 protein among various species. Western blotting showed that ITGβ5 protein was detectable in bull sperm. Moreover, results of RT-qPCR showed that levels of ITGβ5 were significantly higher in the two-cell embryos, followed by the 8-16-cell embryos. However, no significant difference in expression levels were noted for the morula and blastocyst stages as compared to MII oocytes. Bioinformatic analysis revealed that ITGβ5 is conserved among various species. We conclude that expression of ITGβ5 in bovine gametes and embryos implies an important role in fertilisation and embryogenesis.

Genetic mechanisms underlying spermatic and testicular traits within and among cattle breeds: systematic review and prioritization of GWAS results

Reduced bull fertility imposes economic losses in bovine herds. Specifically, testicular and spermatic traits are important indicators of reproductive efficiency. Several genome-wide association studies (GWAS) have identified genomic regions associated with these fertility traits. The aims of this study were as follows: 1) to perform a systematic review of GWAS results for spermatic and testicular traits in cattle and 2) to identify key functional candidate genes for these traits. The identification of functional candidate genes was performed using a systems biology approach, where genes shared between traits and studies were evaluated by a guilt by association gene prioritization (GUILDify and ToppGene software) in order to identify the best functional candidates. These candidate genes were integrated and analyzed in order to identify overlapping patterns among traits and breeds. Results showed that GWAS for testicular-related traits have been developed for beef breeds only, whereas the majority of GWAS for spermatic-related traits were conducted using dairy breeds. When comparing traits measured within the same study, the highest number of genes shared between different traits was observed, indicating a high impact of the population genetic structure and environmental effects. Several chromosomal regions were enriched for functional candidate genes associated with fertility traits. Moreover, multiple functional candidate genes were enriched for markers in a species-specific basis, taurine (Bos taurus) or indicine (Bos indicus). For the different candidate regions identified in the GWAS in the literature, functional candidate genes were detected as follows: B. Taurus chromosome X (BTX) (TEX11, IRAK, CDK16, ATP7A, ATRX, HDAC6, FMR1, L1CAM, MECP2, etc.), BTA17 (TRPV4 and DYNLL1), and BTA14 (MOS, FABP5, ZFPM2). These genes are responsible for regulating important metabolic pathways or biological processes associated with fertility, such as progression of spermatogenesis, control of ciliary activity, development of Sertoli cells, DNA integrity in spermatozoa, and homeostasis of testicular cells. This study represents the first systematic review on male fertility traits in cattle using a system biology approach to identify key candidate genes for these traits.

Review: Genomics of bull fertility

Fertility is one of the most economically important traits in both beef and dairy cattle production; however, only female fertility is typically subjected to selection. Male and female fertility have only a small positive genetic correlation which is likely due to the existence of a relatively small number of genetic variants within each breed that cause embryonic and developmental losses. Genomic tools have been developed that allow the identification of lethal recessive loci based upon marker haplotypes. Selection against haplotypes harbouring lethal alleles in conjunction with selection to improve female fertility will result in an improvement in male fertility. Genomic selection has resulted in a two to fourfold increase in the rate of genetic improvement of most dairy traits in US Holstein cattle, including female fertility. Considering the rapidly increasing rate of adoption of high-throughput single nucleotide polymorphism genotyping in both the US dairy and beef industries, genomic selection should be the most effective of all currently available approaches to improve male fertility. However, male fertility phenotypes are not routinely recorded in natural service mating systems and when artificial insemination is used, semen doses may be titrated to lower post-thaw progressively motile sperm numbers for high-merit and high-demand bulls. Standardization of sperm dosages across bull studs for semen distributed from young bulls would allow the capture of sire conception rate phenotypes for young bulls that could be used to generate predictions of genetic merit for male fertility in both males and females. These data would allow genomic selection to be implemented for male fertility in addition to female fertility within the US dairy industry. While the rate of use of artificial insemination is much lower within the US beef industry, the adoption of sexed semen in the dairy industry has allowed dairy herds to select cows from which heifer replacements are produced and cows that are used to produce terminal crossbred bull calves sired by beef breed bulls. Capture of sire conception rate phenotypes in dairy herds utilizing sexed semen will contribute data enabling genomic selection for male fertility in beef cattle breeds. As the commercial sector of the beef industry increasingly adopts fixed-time artificial insemination, sire conception rate phenotypes can be captured to facilitate the development of estimates of genetic merit for male fertility within US beef breeds.

Whole-genome scan reveals significant non-additive effects for sire conception rate in Holstein cattle

Background

Service sire has a considerable impact on reproductive success in dairy cattle. Most gene mapping studies for bull fertility have focused on additive effects, while non-additive effects have been largely ignored. The main goal of this study was to assess the relevance of non-additive effects on Sire Conception Rate (SCR) in Holstein dairy cattle. The analysis included 7.5 k Holstein bulls with both SCR records and 57.8 k single nucleotide polymorphism (SNP) markers spanning the entire genome.

Results

The importance of non-additive effects was evaluated using an efficient two-step mixed model-based approach. Four genomic regions located on chromosomes BTA8, BTA9, BTA13 and BTA17 showed marked dominance and/or recessive effects. Most of these regions harbor genes, such as ADAM28, DNAJA1, TBC1D20, SPO11, PIWIL3 and TMEM119, that are directly implicated in testis development, male germ line maintenance, and sperm maturation.

Conclusions

This study provides further evidence for the relevance of non-additive effects in fitness-related traits, such as male fertility. In addition, these findings may point out new strategies for improving service sire fertility in dairy cattle via marker-assisted selection.

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.

Identification of genomic variants causing sperm abnormalities and reduced male fertility

Whole genome sequencing has identified millions of bovine genetic variants; however, there is currently little understanding about which variants affect male fertility. It is imperative that we begin to link detrimental genetic variants to sperm phenotypes via the analysis of semen samples and measurement of fertility for bulls with alternate genotypes. Artificial insemination (AI) bulls provide a useful model system because of extensive fertility records, measured as sire conception rates (SCR). Genetic variants with moderate to large effects on fertility can be identified by sequencing the genomes of fertile and subfertile or infertile sires identified with high or low SCR as adult AI bulls or yearling bulls that failed Breeding Soundness Evaluation. Variants enriched in frequency in the sequences of subfertile/infertile bulls, particularly those likely to result in the loss of protein function or predicted to be severely deleterious to genes involved in sperm protein structure and function, semen quality or sperm morphology can be designed onto genotyping assays for validation of their effects on fertility. High throughput conventional and image-based flow cytometry, proteomics and cell imaging can be used to establish the functional effects of variants on sperm phenotypes. Integrating the genetic, fertility and sperm phenotype data will accelerate biomarker discovery and validation, improve routine semen testing in bull studs and identify new targets for cost-efficient AI dose optimization approaches such as semen nanopurification. This will maximize semen output from genetically superior sires and will increase the fertility of cattle. Better understanding of the relationships between male genotype and sperm phenotype may also yield new diagnostic tools and treatments for human male and idiopathic infertility.

Candidate genes for male and female reproductive traits in Canchim beef cattle

Background

Beef cattle breeding programs in Brazil have placed greater emphasis on the genomic study of reproductive traits of males and females due to their economic importance. In this study, genome-wide associations were assessed for scrotal circumference at 210 d of age, scrotal circumference at 420 d of age, age at first calving, and age at second calving, in Canchim beef cattle. Data quality control was conducted resulting in 672,778 SNPs and 392 animals.

Results

Associated SNPs were observed for scrotal circumference at 420 d of age (435 SNPs), followed by scrotal circumference at 210 d of age (12 SNPs), age at first calving (six SNPs), and age at second calving (four SNPs). We investigated whether significant SNPs were within genic or surrounding regions. Biological processes of genes were associated with immune system, multicellular organismal process, response to stimulus, apoptotic process, cellular component organization or biogenesis, biological adhesion, and reproduction.

Conclusions

Few associations were observed for scrotal circumference at 210 d of age, age at first calving, and age at second calving, reinforcing their polygenic inheritance and the complexity of understanding the genetic architecture of reproductive traits. Finding many associations for scrotal circumference at 420 d of age in various regions of the Canchim genome also reveals the difficulty of targeting specific candidate genes that could act on fertility; nonetheless, the high linkage disequilibrium between loci herein estimated could aid to overcome this issue. Therefore, all relevant information about genomic regions influencing reproductive traits may contribute to target candidate genes for further investigation of causal mutations and aid in future genomic studies in Canchim cattle to improve the breeding program.

Electronic supplementary material

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

Oligo-dT selected spermatozoal transcript profiles differ among higher and lower fertility dairy sires

Spermatozoal messenger RNA (mRNA) has the potential as a molecular marker for sire fertility because this population can reflect gene expression that occurred during spermatogenesis and may have a functional role in early embryonic development. The goal of this study was to compare the oligo-dT selected spermatozoal transcript profiles of higher fertility (Conception Rate (CR) 1.8-3.5) and lower fertility (CR -2.9 to -0.4) sires using Ribonucleic Acid Sequencing (RNA-Seq). A total of 3227 transcripts and 5366 transcripts were identified in the higher and lower fertility populations, respectively. While common transcripts between the two populations were identified (2422 transcripts), several transcripts were also unique to the fertility populations including 805 transcripts that were unique to the higher fertility population and 2944 transcripts that were unique to the lower fertility population. From gene ontological analysis, the transcripts unique to each fertility population differed in Biological Processes (BP), including enrichment of regulatory transcripts for growth and protein kinase activity in the higher fertility bulls. Biological variation in transcript presence among individual sires was also found. Of the candidate fertility spermatozoal transcripts chosen from the RNA-Seq population analysis reported here and previous publications, COX7C was negatively correlated with sire fertility. Using high-throughput sequencing, candidate spermatozoal transcripts were identified for further study as potential markers for sire fertility.

Unravelling the genomic architecture of bull fertility in Holstein cattle

BACKGROUND

Fertility is considered an important economic trait in dairy cattle. Most studies have investigated cow fertility while bull fertility has received much less consideration. The main objective of this study was to perform a comprehensive genomic analysis in order to unravel the genomic architecture underlying sire fertility in Holstein dairy cattle. The analysis included the application of alternative genome-wide association mapping approaches and the subsequent use of diverse gene set enrichment tools.

RESULTS

The association analyses identified at least eight genomic regions strongly associated with bull fertility. Most of these regions harbor genes, such as KAT8, CKB, TDRD9 and IGF1R, with functions related to sperm biology, including sperm development, motility and sperm-egg interaction. Moreover, the gene set analyses revealed many significant functional terms, including fertilization, sperm motility, calcium channel regulation, and SNARE proteins. Most of these terms are directly implicated in sperm physiology and male fertility.

CONCLUSIONS

This study contributes to the identification of genetic variants and biological processes underlying sire fertility. These findings can provide opportunities for improving bull fertility via marker-assisted selection.

Pathway-Based Genome-Wide Association Studies for Two Meat Production Traits in Simmental Cattle

Most single nucleotide polymorphisms (SNPs) detected by genome-wide association studies (GWAS), explain only a small fraction of phenotypic variation. Pathway-based GWAS were proposed to improve the proportion of genes for some human complex traits that could be explained by enriching a mass of SNPs within genetic groups. However, few attempts have been made to describe the quantitative traits in domestic animals. In this study, we used a dataset with approximately 7,700,000 SNPs from 807 Simmental cattle and analyzed live weight and longissimus muscle area using a modified pathway-based GWAS method to orthogonalise the highly linked SNPs within each gene using principal component analysis (PCA). As a result, of the 262 biological pathways of cattle collected from the KEGG database, the gamma aminobutyric acid (GABA)ergic synapse pathway and the non-alcoholic fatty liver disease (NAFLD) pathway were significantly associated with the two traits analyzed. The GABAergic synapse pathway was biologically applicable to the traits analyzed because of its roles in feed intake and weight gain. The proposed method had high statistical power and a low false discovery rate, compared to those of the smallest P-value and SNP set enrichment analysis methods.

Differential histone modification status of spermatozoa in relation to fertility of buffalo bulls

In this study genome-wide di-methylated H3K4 (H3K4me2) and tri-methylated H3K27 (H3K27me3) modification profiles were analyzed in spermatozoa of buffalo bulls having wide fertility differences. The custom designed 4 × 180 K buffalo (Bubalus bubalis) ChIP-on-chip array was fabricated by employing array-based sequential hybridization using bovine and buffalo genomic DNA for comparative hybridization. The buffalo specific array developed had 177,440 features assembled from Coding sequences, Promoter and CpG regions comprising 2967 unique genes. A total of 84 genes for H3K4me2 and 80 genes for H3K27me3 were found differentially enriched in mature sperm of high and sub-fertile buffalo bulls. Gene Ontology analysis of these genes revealed their association with different cellular functions and biological processes. Genes identified as differentially enriched between high and sub-fertile bulls were found to be involved in the processes of germ cell development, spermatogenesis and embryonic development. This study presents the first genome-wide H3K4me2 and H3K27me3 profiling of buffalo bull sperm. Results provide a list of specific genes which could be made responsible for differential bull fertility.

Genomic prediction and genome-wide association study for dagginess and host internal parasite resistance in New Zealand sheep

Background

Dagginess (faecal soiling of the perineum region) and host nematode parasite resistance are important animal welfare traits in New Zealand sheep. Genomic prediction (GP) estimates the genetic merit, as a molecular breeding value (mBV), for each trait based on many SNPs. The additional information the mBV provides (as determined by its accuracy) has led to its incorporation into breeding schemes. Some GP methods give SNP effects, which provide additional information to identify genome-wide associations (GWAS) for a trait of interest. Here we report results from a GP and GWAS study for dagginess and host nematode parasite resistance in a New Zealand sheep industry resource.

Results

Genomic prediction analysis was performed using 50K SNP chip data and parent average-removed, de-regressed BVs for five traits, from a resource of 8705 pedigree recorded animals. The five traits were dag score at three and eight months (DAG3, DAG8) and nematode faecal egg count in summer (FEC1), autumn (FEC2) and as an adult (AFEC). The resource consisted of Romney, Coopworth, Perendale, Texel and various breed crosses (designated: CompRCP, CompRCPT and CompCRP). The pure breeds, apart from Texel, plus CompRCP were used to develop the GP. The resulting SNP effects were used to identify genetic regions associated with dagginess and parasite resistance. Accuracies of the weighted correlation between mBV and true BV ranged between −0.07 (Texel) and 0.56 (Coopworth) for DAG3 and DAG8. For FEC1, FEC2 and AFEC accuracies ranged between −0.22 (CompRCPT) and 0.69 (Coopworth). The weighted average individual accuracy (calculated from theory) ranges were 0.13 (Texel) to 0.52 (Coopworth) and 0.11 (Texel) to 0.55 (Coopworth) respectively, for dagginess and parasite traits. There was one SNP for DAG8 that reached Bonferroni significance threshold (P < 1 × 10⁻⁶) on OAR15, the same two SNPs for each of the parasite traits (OAR26) and none for DAG3. A notable peak was also observed on OAR7 for all the parasite traits, however, it did not reach the Bonferroni significance threshold.

Conclusions

This study presents the first results of a GWAS on dagginess and faecal egg count traits in New Zealand sheep. The results suggest that there are quantitative trait loci on OAR 15 for dagginess and on OAR26 and seven for faecal egg count.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2148-2) contains supplementary material, which is available to authorized users.

Identification of putative drug targets for human sperm-egg interaction defect using protein network approach

Background

Sperm-egg interaction defect is a significant cause of in-vitro fertilization failure for infertile cases. Numerous molecular interactions in the form of protein-protein interactions mediate the sperm-egg membrane interaction process. Recent studies have demonstrated that in addition to experimental techniques, computational methods, namely protein interaction network approach, can address protein-protein interactions between human sperm and egg. Up to now, no drugs have been detected to treat sperm-egg interaction disorder, and the initial step in drug discovery research is finding out essential proteins or drug targets for a biological process. The main purpose of this study is to identify putative drug targets for human sperm-egg interaction deficiency and consider if the detected essential proteins are targets for any known drugs using protein-protein interaction network and ingenuity pathway analysis.

Results

We have created human sperm-egg protein interaction networks with high confidence, including 106 nodes and 415 interactions. Through topological analysis of the network with calculation of some metrics, such as connectivity and betweenness centrality, we have identified 13 essential proteins as putative drug targets. The potential drug targets are from integrins, fibronectins, epidermal growth factor receptors, collagens and tetraspanins protein families. We evaluated these targets by ingenuity pathway analysis, and the known drugs for the targets have been detected, and the possible effective role of the drugs on sperm-egg interaction defect has been considered. These results showed that the drugs ocriplasmin (Jetrea©), gefitinib (Iressa©), erlotinib hydrochloride (Tarceva©), clingitide, cetuximab (Erbitux©) and panitumumab (Vectibix©) are possible candidates for efficacy testing for the treatment of sperm-egg interaction deficiency. Further experimental validation can be carried out to confirm these results.

Conclusion

We have identified the first potential list of drug targets for human sperm-egg interaction defect using the protein interaction network approach. The essential proteins or potential drug targets were found using topological analysis of the protein network. These putative targets are promising for further experimental validation. These study results, if validated, may develop drug discovery techniques for sperm-egg interaction defect and also improve assisted reproductive technologies to avoid in-vitro fertilization failure.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0186-7) contains supplementary material, which is available to authorized users.

Differential proteomic profile of spermatogenic and Sertoli cells from peri-pubertal testes of three different bovine breeds

Sub-fertility is one of the most common problems observed in crossbred males, but the etiology remains unknown in most of the cases. Although proteomic differences in the spermatozoa and seminal plasma between breeds have been investigated, the possible differences at the sperm precursor cells and supporting/nourishing cells have not been studied. The present study reports the differential proteomic profile of spermatogenic and Sertoli cells in crossbred and purebred bulls. Testis was removed by unilateral castration of 12 peri-pubertal bulls (10 months age), four each from crossbred (Holstein Friesian × Tharparkar), exotic purebred [Holstein Friesian (HF)] and indigenous purebred [Tharparkar (TP)] bulls. Spermatogenic and Sertoli cells were isolated and subjected to proteomic analysis. Protein extracts from the Sertoli and spermatogenic cells of each breed were analyzed with 2-dimensional difference gel electrophoresis (2D-DIGE) and analyzed with Decyder™ software. Compared to HF, 26 protein spots were over expressed and 14 protein spots were under expressed in spermatogenic cells of crossbred bulls. Similarly, 7 protein spots were over expressed and 15 protein spots were under expressed in the spermatogenic cells of TP bulls compared to that of crossbred bulls. Out of 12 selected protein spots identified through mass spectrometry, Phosphatidyl ethanolamine binding protein was found to be over expressed in the spermatogenic cells of crossbred bulls compared to TP bulls. The protein, gamma actin was found to be over expressed in the Sertoli cells of HF bulls, whereas Speedy Protein-A was found to be over expressed in Sertoli cells of crossbred bulls. It may be concluded that certain proteomic level differences exist in sperm precursor cells and nourishing cells between breeds, which might be associated with differences in the fertility among these breeds.

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.

Genome wide association studies for body conformation traits in the Chinese Holstein cattle population

BACKGROUND

Genome-wide association study (GWAS) is a powerful tool for revealing the genetic basis of quantitative traits. However, studies using GWAS for conformation traits of cattle is comparatively less. This study aims to use GWAS to find the candidates genes for body conformation traits.

RESULTS

The Illumina BovineSNP50 BeadChip was used to identify single nucleotide polymorphisms (SNPs) that are associated with body conformation traits. A least absolute shrinkage and selection operator (LASSO) was applied to detect multiple SNPs simultaneously for 29 body conformation traits with 1,314 Chinese Holstein cattle and 52,166 SNPs. Totally, 59 genome-wide significant SNPs associated with 26 conformation traits were detected by genome-wide association analysis; five SNPs were within previously reported QTL regions (Animal Quantitative Trait Loci (QTL) database) and 11 were very close to the reported SNPs. Twenty-two SNPs were located within annotated gene regions, while the remainder were 0.6-826 kb away from known genes. Some of the genes had clear biological functions related to conformation traits. By combining information about the previously reported QTL regions and the biological functions of the genes, we identified DARC, GAS1, MTPN, HTR2A, ZNF521, PDIA6, and TMEM130 as the most promising candidate genes for capacity and body depth, chest width, foot angle, angularity, rear leg side view, teat length, and animal size traits, respectively. We also found four SNPs that affected four pairs of traits, and the genetic correlation between each pair of traits ranged from 0.35 to 0.86, suggesting that these SNPs may have a pleiotropic effect on each pair of traits.

CONCLUSIONS

A total of 59 significant SNPs associated with 26 conformation traits were identified in the Chinese Holstein population. Six promising candidate genes were suggested, and four SNPs showed genetic correlation for four pairs of traits.

Polymorphism of growth hormone gene and its association with seminal and sexual behavioral traits in crossbred cattle

The decline in the male reproductive ability in terms of sexual behavior and seminal traits might lead to nonavailability of required number of bulls in a progeny testing program. The present study was conducted in 493 crossbred cattle (Bos taurus × Bos indicus) bulls to study polymorphisms of growth hormone (GH) gene and its association with seminal and sexual behavioral characteristics. A 428-base pair fragment of GH gene spanning over the fourth exon, fourth intron, and fifth exon was amplified and digested with AluI restriction enzyme. Bulls were found to be polymorphic, with two variants, LL and LV, and higher genotypic frequency for LL being 0.88. Twelve nucleotide changes and a single nonsynonymous substitution of Leucine by Valine were observed from GH1 (L) to GH2 (V). Statistical analysis revealed that the genotype of the GH gene had a significant effect on libido score, reaction time, Flehmen response, requirement of mounting stimulus, sperm mass activity, number of semen doses per collection, individual fresh sperm motility, postthaw sperm motility, acrosome integrity, hypo-osmotic swelling test, live and dead count, total morphological abnormality, and head abnormality of sperm in crossbred bulls. Growth hormone gene might be considered a candidate gene for seminal and sexual behavioral traits in crossbred cattle.

Discovery of single nucleotide polymorphisms in candidate genes associated with fertility and production traits in Holstein cattle

Background

Identification of single nucleotide polymorphisms (SNPs) for specific genes involved in reproduction might improve reliability of genomic estimates for these low-heritability traits. Semen from 550 Holstein bulls of high (≥ 1.7; n = 288) or low (≤ −2; n = 262) daughter pregnancy rate (DPR) was genotyped for 434 candidate SNPs using the Sequenom MassARRAY® system. Three types of SNPs were evaluated: SNPs previously reported to be associated with reproductive traits or physically close to genetic markers for reproduction, SNPs in genes that are well known to be involved in reproductive processes, and SNPs in genes that are differentially expressed between physiological conditions in a variety of tissues associated in reproductive function. Eleven reproduction and production traits were analyzed.

Results

A total of 40 SNPs were associated (P < 0.05) with DPR. Among these were genes involved in the endocrine system, cell signaling, immune function and inhibition of apoptosis. A total of 10 genes were regulated by estradiol. In addition, 22 SNPs were associated with heifer conception rate, 33 with cow conception rate, 36 with productive life, 34 with net merit, 23 with milk yield, 19 with fat yield, 13 with fat percent, 19 with protein yield, 22 with protein percent, and 13 with somatic cell score. The allele substitution effect for SNPs associated with heifer conception rate, cow conception rate, productive life and net merit were in the same direction as for DPR. Allele substitution effects for several SNPs associated with production traits were in the opposite direction as DPR. Nonetheless, there were 29 SNPs associated with DPR that were not negatively associated with production traits.

Conclusion

SNPs in a total of 40 genes associated with DPR were identified as well as SNPs for other traits. It might be feasible to include these SNPs into genomic tests of reproduction and other traits. The genes associated with DPR are likely to be important for understanding the physiology of reproduction. Given the large number of SNPs associated with DPR that were not negatively associated with production traits, it should be possible to select for DPR without compromising production.

Genome-wide association study for fatty acid composition in Japanese Black cattle

Fatty acid composition is one of the important traits in beef. The aim of this study was to identify candidate genomic regions for fatty acid composition by genome-wide association study with 50 K single nucleotide polymorphism (SNP) array in Japanese Black cattle. A total of 461 individuals and 40 657 SNPs were used in this study. We applied genome-wide rapid association using mixed model and regression (GRAMMAR) and genomic control approaches to estimate the associations between genotypes and fatty acid composition. In addition, two SNPs in fatty acid synthase (FASN) (T1952A) and stearoyl-CoA desaturase (SCD) (V293A) genes were also genotyped. Association analysis revealed that 30 significant SNPs for several fatty acids (C14:0, C14:1, C16:1 and C18:1) were located in the BTA19 FASN gene located within this region but the FASN mutation had no significant effect on any traits. We also detected one significant SNP for C18:1 on BTA23 and two SNPs for C16:0 on BTA25. The region around 17 Mb on BTA26 harbored two significant SNPs for C14:1 and SNP in SCD in this region showed the strongest association with C14:1. This study demonstrated novel candidate regions in BTA19, 23 and 25 for fatty acid composition.

Stallion sperm transcriptome comprises functionally coherent coding and regulatory RNAs as revealed by microarray analysis and RNA-seq

Mature mammalian sperm contain a complex population of RNAs some of which might regulate spermatogenesis while others probably play a role in fertilization and early development. Due to this limited knowledge, the biological functions of sperm RNAs remain enigmatic. Here we report the first characterization of the global transcriptome of the sperm of fertile stallions. The findings improved understanding of the biological significance of sperm RNAs which in turn will allow the discovery of sperm-based biomarkers for stallion fertility. The stallion sperm transcriptome was interrogated by analyzing sperm and testes RNA on a 21,000-element equine whole-genome oligoarray and by RNA-seq. Microarray analysis revealed 6,761 transcripts in the sperm, of which 165 were sperm-enriched, and 155 were differentially expressed between the sperm and testes. Next, 70 million raw reads were generated by RNA-seq of which 50% could be aligned with the horse reference genome. A total of 19,257 sequence tags were mapped to all horse chromosomes and the mitochondrial genome. The highest density of mapped transcripts was in gene-rich ECA11, 12 and 13, and the lowest in gene-poor ECA9 and X; 7 gene transcripts originated from ECAY. Structural annotation aligned sperm transcripts with 4,504 known horse and/or human genes, rRNAs and 82 miRNAs, whereas 13,354 sequence tags remained anonymous. The data were aligned with selected equine gene models to identify additional exons and splice variants. Gene Ontology annotations showed that sperm transcripts were associated with molecular processes (chemoattractant-activated signal transduction, ion transport) and cellular components (membranes and vesicles) related to known sperm functions at fertilization, while some messenger and micro RNAs might be critical for early development. The findings suggest that the rich repertoire of coding and non-coding RNAs in stallion sperm is not a random remnant from spermatogenesis in testes but a selectively retained and functionally coherent collection of RNAs.

Application of quantum dot nanoparticles for potential non-invasive bio-imaging of mammalian spermatozoa

Background

Various obstacles are encountered by mammalian spermatozoa during their journey through the female genital tract, and only few or none will reach the site of fertilization. Currently, there are limited technical approaches for non-invasive investigation of spermatozoa migration after insemination. As the knowledge surrounding sperm behavior throughout the female genital tract still remains elusive, the recent development of self-illuminating quantum dot nanoparticles may present a potential means for real-time in vitro and in vivo monitoring of spermatozoa.

Results

Here, we show the ability of boar spermatozoa to harmlessly interact and incorporate bioluminescent resonance energy transfer-conjugated quantum dot (BRET-QD) nanoparticles. The confocal microscope revealed in situ fluorescence of BRET-QD in the entire spermatozoon, while the ultra-structural analysis using the transmission electron microscope indicated BRET-QD localization on the sperm plasma membrane and intracellular compartment. In controlled-in vitro assays, bioluminescent imaging demonstrated that spermatozoa incubated with BRET-QD and luciferase substrate (coelenterazine) emit light (photons/sec) above the background, which confirmed the in situ fluorescence imaging. Most importantly, sperm motility, viability, and fertilizing potential were not affected by the BRET-QD incorporation when used at an appropriated ratio.

Conclusions

Our results demonstrate that pig spermatozoa can incorporate BRET-QD nanoparticles without affecting their motility and capacity to interact with the oocyte when used at an appropriated balance. We anticipate that our study will enable in-depth exploration of the male components of in vivo migration, fertilization, and embryonic development at the molecular level using this novel approach.

Allelic gene expression imbalance of bovine IGF2, LEP and CCL2 genes in liver, kidney and pituitary

Allelic expression imbalance (AEI) is an important genetic factor being the cause of differences in phenotypic traits that can be heritable. Studying AEI can be useful in searching for factors that modulate gene expression and help to understand molecular mechanisms underlying phenotypic changes. Although it was commonly recognized in many species and we know many genes show allelic expression imbalance, this phenomena was not studied on a larger scale in cattle. Using the pyrosequencing method we analyzed a set of 29 bovine genes in order to find those that have preferential allelic expression. The study was conducted in three tissues: liver, pituitary and kindey. Out of the studied group of genes 3 of them—LEP (leptin), IGF2 (insulin-like growth factor 2), CCL2 (chemokine C–C motif ligand 2) showed allelic expression imbalance.

Genome-wide association study identifies candidate markers for bull fertility in Holstein dairy cattle

The decline in the reproductive efficiency of dairy cattle has become a challenging problem worldwide. Female fertility is now taken into account in breeding goals while generally less attention is given to male fertility. The objective of this study was to perform a genome-wide association study in Holstein bulls to identify genetic variants significantly related to sire conception rate (SCR), a new phenotypic evaluation of bull fertility. The analysis included 1755 sires with SCR data and 38,650 single nucleotide polymorphisms (SNPs) spanning the entire bovine genome. Associations between SNPs and SCR were analyzed using a mixed linear model that included a random polygenic effect and SNP genotype either as a linear covariate or as a categorical variable. A multiple testing correction approach was used to account for the correlation between SNPs because of linkage disequilibrium. After genome-wide correction, eight SNPs showed significant association with SCR. Some of these SNPs are located close to or in the middle of genes with functions related to male fertility, such as the sperm acrosome reaction, chromatin remodeling during the spermatogenesis, and the meiotic process during male germ cell maturation. Some SNPs showed marked dominance effects, which provide more evidence for the relevance of non-additive effects in traits closely related to fitness such as fertility. The results could contribute to the identification of genes and pathways associated with male fertility in dairy cattle.

Genome-Wide SNP Discovery from Transcriptome of Four Common Carp Strains

Background

Single nucleotide polymorphisms (SNPs) have been used as genetic marker for genome-wide association studies in many species. Gene-associated SNPs could offer sufficient coverage in trait related research and further more could themselves be causative SNPs for traits. Common carp (Cyprinus carpio) is one of the most important aquaculture species in the world accounting for nearly 14% of freshwater aquaculture production. There are various strains of common carp with different economic traits, however, the genetic mechanism underlying the different traits have not been elucidated yet. In this project, we identified a large number of gene-associated SNPs from four strains of common carp using next-generation sequencing.

Results

Transcriptome sequencing of four strains of common carp (mirror carp, purse red carp, Xingguo red carp, Yellow River carp) was performed with Solexa HiSeq2000 platform. De novo assembled transcriptome was used as reference for alignments, and SNP calling was done through BWA and SAMtools. A total of 712,042 Intra-strain SNPs were discovered in four strains, of which 483,276 SNPs for mirror carp, 486,629 SNPs for purse red carp, 478,028 SNPs for Xingguo red carp and 488,281 SNPs for Yellow River carp were discovered, respectively. Besides, 53,893 inter-SNPs were identified. Strain-specific SNPs of four strains were 53,938, 53,866, 48,701, 40,131 in mirror carp, purse red carp, Xingguo red carp and Yellow River carp, respectively. GO and KEGG pathway analysis were done to reveal strain-specific genes affected by strain-specific non-synonymous SNPs. Validation of selected SNPs revealed that 48% percent of SNPs (12 of 25) were tested to be true SNPs.

Conclusions

Transcriptome analysis of common carp using RNA-Seq is a cost-effective way of generating numerous reads for SNP discovery. After validation of identified SNPs, these data will provide a solid base for SNP array designing and genome-wide association studies.

Progress of genome wide association study in domestic animals

Domestic animals are invaluable resources for study of the molecular architecture of complex traits. Although the mapping of quantitative trait loci (QTL) responsible for economically important traits in domestic animals has achieved remarkable results in recent decades, not all of the genetic variation in the complex traits has been captured because of the low density of markers used in QTL mapping studies. The genome wide association study (GWAS), which utilizes high-density single-nucleotide polymorphism (SNP), provides a new way to tackle this issue. Encouraging achievements in dissection of the genetic mechanisms of complex diseases in humans have resulted from the use of GWAS. At present, GWAS has been applied to the field of domestic animal breeding and genetics, and some advances have been made. Many genes or markers that affect economic traits of interest in domestic animals have been identified. In this review, advances in the use of GWAS in domestic animals are described.

Genome-wide associations for milk production and somatic cell score in Holstein-Friesian cattle in Ireland

BACKGROUND

Contemporary dairy breeding goals have broadened to include, along with milk production traits, a number of non-production-related traits in an effort to improve the overall functionality of the dairy cow. Increased indirect selection for resistance to mastitis, one of the most important production-related diseases in the dairy sector, via selection for reduced somatic cell count has been part of these broadened goals. A number of genome-wide association studies have identified genetic variants associated with milk production traits and mastitis resistance, however the majority of these studies have been based on animals which were predominantly kept in confinement and fed a concentrate-based diet (i.e. high-input production systems). This genome-wide association study aims to detect associations using genotypic and phenotypic data from Irish Holstein-Friesian cattle fed predominantly grazed grass in a pasture-based production system (low-input).

RESULTS

Significant associations were detected for milk yield, fat yield, protein yield, fat percentage, protein percentage and somatic cell score using separate single-locus, frequentist and multi-locus, Bayesian approaches. These associations were detected using two separate populations of Holstein-Friesian sires and cows. In total, 1,529 and 37 associations were detected in the sires using a single SNP regression and a Bayesian method, respectively. There were 103 associations in common between the sires and cows across all the traits. As well as detecting associations within known QTL regions, a number of novel associations were detected; the most notable of these was a region of chromosome 13 associated with milk yield in the population of Holstein-Friesian sires.

CONCLUSIONS

A total of 276 of novel SNPs were detected in the sires using a single SNP regression approach. Although obvious candidate genes may not be initially forthcoming, this study provides a preliminary framework upon which to identify the causal mechanisms underlying the various milk production traits and somatic cell score. Consequently this will deepen our understanding of how these traits are expressed.