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

Genome-Wide Association Study of Reproductive Traits in Large White Pigs

(1) Background: Reproductive performance is crucial for the pork industry's success. The Large White pig is central to this, yet the genetic factors influencing its reproductive traits are not well understood, highlighting the need for further research. (2) Methods: This study utilized Genome-Wide Association Studies to explore the genetic basis of reproductive traits in the Large White pig. We collected data from 2237 Large White sows across four breeding herds in southern China, focusing on eight reproductive traits. Statistical analyses included principal component analysis, linkage disequilibrium analysis, and univariate linear mixed models to identify significant single-nucleotide polymorphisms and candidate genes. (3) Results: Forty-five significantly related SNPs and 17 potential candidate genes associated with litter traits were identified. Individuals with the TT genotype at SNP rs341909772 showed an increase of 1.24 in the number of piglets born alive (NBA) and 1.25 in the number of healthy births (NHBs) compared with those with the CC genotype. (4) Conclusions: The SNPs and genes identified in this study offer insights into the genetics of reproductive traits in the Large White pig, potentially guiding the development of breeding strategies to improve litter size.

Identification of Copy Number Variations and Selection Signatures in Wannan Spotted Pigs by Whole Genome Sequencing Data: A Preliminary Study

Copy number variation (CNV) is an important structural variation used to elucidate complex economic traits. In this study, we sequenced 25 Wannan spotted pigs (WSPs) to detect their CNVs and identify their selection signatures compared with those of 10 Asian wild boars. A total of 14,161 CNVs were detected in the WSPs, accounting for 0.72% of the porcine genome. The fixation index (Fst) was used to identify the selection signatures, and 195 CNVs with the top 1% of the Fst value were selected. Eighty genes were identified in the selected CNV regions. Functional GO and KEGG analyses revealed that the genes within these selected CNVs are associated with key traits such as reproduction (GAL3ST1 and SETD2), fatty acid composition (PRKG1, ACACA, ACSL3, UGT8), immune system (LYZ), ear size (WIF1), and feed efficiency (VIPR2). The findings of this study contribute novel insights into the genetic CNVs underlying WSP characteristics and provide essential information for the protection and utilization of WSP populations.

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.

Review: Genetic mutations affecting bull fertility

Cattle are a well-suited "model organism" to study the genetic underpinnings of variation in male reproductive performance. The adoption of artificial insemination and genomic prediction in many cattle breeds provide access to microarray-derived genotypes and repeated measurements for semen quality and insemination success in several thousand bulls. Similar-sized mapping cohorts with phenotypes for male fertility are not available for most other species precluding powerful association testing. The repeated measurements of the artificial insemination bulls' semen quality enable the differentiation between transient and biologically relevant trait fluctuations, and thus, are an ideal source of phenotypes for variance components estimation and genome-wide association testing. Genome-wide case-control association testing involving bulls with either aberrant sperm quality or low insemination success revealed several causal recessive loss-of-function alleles underpinning monogenic reproductive disorders. These variants are routinely monitored with customised genotyping arrays in the male selection candidates to avoid the use of subfertile or infertile bulls for artificial insemination and natural service. Genome-wide association studies with quantitative measurements of semen quality and insemination success revealed quantitative trait loci for male fertility, but the underlying causal variants remain largely unknown. Moreover, these loci explain only a small part of the heritability of male fertility. Integrating genome-wide association studies with gene expression and other omics data from male reproductive tissues is required for the fine-mapping of candidate causal variants underlying variation in male reproductive performance in cattle.

Genome-wide association study to identify genomic regions and positional candidate genes associated with male fertility in beef cattle

Fertility plays a key role in the success of calf production, but there is evidence that reproductive efficiency in beef cattle has decreased during the past half-century worldwide. Therefore, identifying animals with superior fertility could significantly impact cow-calf production efficiency. The objective of this research was to identify candidate regions affecting bull fertility in beef cattle and positional candidate genes annotated within these regions. A GWAS using a weighted single-step genomic BLUP approach was performed on 265 crossbred beef bulls to identify markers associated with scrotal circumference (SC) and sperm motility (SM). Eight windows containing 32 positional candidate genes and five windows containing 28 positional candidate genes explained more than 1% of the genetic variance for SC and SM, respectively. These windows were selected to perform gene annotation, QTL enrichment, and functional analyses. Functional candidate gene prioritization analysis revealed 14 prioritized candidate genes for SC of which MAP3K1 and VIP were previously found to play roles in male fertility. A different set of 14 prioritized genes were identified for SM and five were previously identified as regulators of male fertility (SOD2, TCP1, PACRG, SPEF2, PRLR). Significant enrichment results were identified for fertility and body conformation QTLs within the candidate windows. Gene ontology enrichment analysis including biological processes, molecular functions, and cellular components revealed significant GO terms associated with male fertility. The identification of these regions contributes to a better understanding of fertility associated traits and facilitates the discovery of positional candidate genes for future investigation of causal mutations and their implications.

Genome-Wide Association Study and Fine Mapping Reveals Candidate Genes for Birth Weight of Yorkshire and Landrace Pigs

Birth weight of pigs is an important economic factor in the livestock industry. The identification of the genes and variants that underlie birth weight is of great importance. In this study, we integrated two genotyping methods, single nucleotide polymorphism (SNP) chip analysis and restriction site associated DNA sequencing (RAD-seq) to genotype genome-wide SNPs. In total, 45,175 and 139,634 SNPs were detected with the SNP chip and RAD-seq, respectively. The genome-wide association study (GWAS) of the combined SNP panels identified two significant loci located at chr1: 97,745,041 and chr4: 112,031,589, that explained 6.36% and 4.25% of the phenotypic variance respectively. To reduce interval containing causal variants, we imputed sequence-level SNPs in the GWAS identified regions and fine-mapped the causative variants into two narrower genomic intervals: a ∼100 kb interval containing 71 SNPs and a broader ∼870 kb interval with 432 SNPs. This fine-mapping highlighted four promising candidate genes, SKOR2, SMAD2, VAV3, and NTNG1. Additionally, the functional genes, SLC25A24, PRMT6 and STXBP3, are also located near the fine-mapping region. These results suggest that these candidate genes may have contribute substantially to the birth weight of pigs.

Weighted Single-Step Genome-Wide Association Study of Semen Traits in Holstein Bulls of China

Efficient production of high-quality semen is a crucial trait in the dairy cattle breeding due to the widespread use of artificial insemination. However, the genetic architecture (e.g., distributions of causal variants and their corresponding effects) underlying such semen quality traits remains unclear. In this study, we performed genome-wide association studies to identify genes associated with five semen quality traits in Chinese Holstein population, including ejaculate volume, progressive sperm motility, sperm concentration, number of sperm, and number of progressive motile sperm. Our dataset consisted of 2,218 Holstein bulls in China with full pedigree information, representing 12 artificial insemination centers, with 1,508 genotyped using the Illumina BovineSNP50 BeadChip. We used a weighted single-step genome-wide association method with 10 adjacent Single nucleotide polymorphisms (SNPs) as sliding windows, which can make use of individuals without genotypes. We considered the top 10 genomic regions in terms of their explained genomic variants as candidate window regions for each trait. In total, we detected 36 window regions related to one or multiple semen traits across 19 chromosomes. Promising candidate genes of PSMB5, PRMT5, ACTB, PDE3A, NPC1, FSCN1, NR5A2, IQCG, LHX8, and DMRT1 were identified in these window regions for these five semen traits. Our findings provided a solid basis for further research into genetic mechanisms underlying semen quality traits, which may contribute to their accurate genomic prediction in Chinese Holstein population.

Validation of 46 loci associated with female fertility traits in cattle

Background

Subfertility is one challenge facing the dairy industry as the average Holstein heifer conception rate (HCR), the proportion of heifers that conceive and maintain a pregnancy per breeding, is estimated at 55–60%. Of the loci associated with HCR, few have been validated in an independent cattle population, limiting their usefulness for selection or furthering our understanding of the mechanisms involved in successful pregnancy. Therefore, the objectives here were to identify loci associated with HCR: 1) to the first artificial insemination (AI) service (HCR1), 2) to repeated AI services required for a heifer to conceive (TBRD) and 3) to validate loci previously associated with fertility. Breeding and health records from 3359 Holstein heifers were obtained after heifers were bred by AI at observed estrus, with pregnancy determined at day 35 via palpation. Heifer DNA was genotyped using the Illumina BovineHD BeadChip, and genome-wide association analyses (GWAA) were performed with additive, dominant and recessive models using the Efficient Mixed Model Association eXpedited (EMMAX) method with a relationship matrix for two phenotypes. The HCR1 GWAA compared heifers that were pregnant after the first AI service (n = 497) to heifers that were open following the first AI service (n = 405), which included those that never conceived. The TBRD GWAA compared only those heifers which did conceive, across variable numbers of AI service (n = 712). Comparison of loci previously associated with fertility, HCR1 or TBRD were considered the same locus for validation when in linkage disequilibrium (D’ > 0.7).

Results

The HCR1 GWAA identified 116, 187 and 28 loci associated (P < 5 × 10^− 8) in additive, dominant and recessive models, respectively. The TBRD GWAA identified 235, 362, and 69 QTL associated (P < 5 × 10^− 8) with additive, dominant and recessive models, respectively. Loci previously associated with fertility were in linkage disequilibrium with 22 loci shared with HCR1 and TBRD, 5 HCR1 and 19 TBRD loci.

Conclusions

Loci associated with HCR1 and TBRD that have been identified and validated can be used to improve HCR through genomic selection, and to better understand possible mechanisms associated with subfertility.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5935-3) contains supplementary material, which is available to authorized users.

Genetic Variation in Pan Species Is Shaped by Demographic History and Harbors Lineage-Specific Functions

Chimpanzees (Pan troglodytes) and bonobos (Pan paniscus) are the closest living relatives of humans, but the two species show distinct behavioral and physiological differences, particularly regarding female reproduction. Despite their recent rapid decline, the demographic histories of the two species have been different during the past 1–2 Myr, likely having an impact on their genomic diversity. Here, we analyze the inferred functional consequences of genetic variation across 69 individuals, making use of the most complete data set of genomes in the Pan clade to date. We test to which extent the demographic history influences the efficacy of purifying selection in these species. We find that small historical effective population sizes (N_e) correlate not only with low levels of genetic diversity but also with a larger number of deleterious alleles in homozygosity and an increased proportion of deleterious changes at low frequencies. To investigate the putative genetic basis for phenotypic differences between chimpanzees and bonobos, we exploit the catalog of putatively deleterious protein-coding changes in each lineage. We show that bonobo-specific nonsynonymous changes are enriched in genes related to age at menarche in humans, suggesting that the prominent physiological differences in the female reproductive system between chimpanzees and bonobos might be explained, in part, by putatively adaptive changes on the bonobo lineage.

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.

Identification of genome regions determining semen quality in Holstein-Friesian bulls using information theory

Use of information theory can be an alternative statistical approach to detect genome regions and candidate genes that are associated with livestock traits. The aim of this study was to verify the validity of the SNPs effects on some semen quality variables of bulls using entropy analysis. Records from 288 Holstein-Friesian bulls from one AI station were included. The following semen quality variables were analyzed: CASA kinematic variables of sperm (total motility, average path velocity, straight line velocity, curvilinear velocity, amplitude of lateral head displacement, beat cross frequency, straightness, linearity), sperm membrane integrity (plazmolema, mitochondrial function), sperm ATP content. Molecular data included 48,192 SNPs. After filtering (call rate = 0.95 and MAF = 0.05), 34,794 SNPs were included in the entropy analysis. The entropy and conditional entropy were estimated for each SNP. Conditional entropy quantifies the remaining uncertainty about values of the variable with the knowledge of SNP. The most informative SNPs for each variable were determined. The computations were performed using the R statistical package. A majority of the loci had relatively small contributions. The most informative SNPs for all variables were mainly located on chromosomes: 3, 4, 5 and 16. The results from the study indicate that important genome regions and candidate genes that determine semen quality variables in bulls are located on a number of chromosomes. Some detected clusters of SNPs were located in RNA (U6 and 5S_rRNA) for all the variables for which analysis occurred. Associations between PARK2 as well GALNT13 genes and some semen characteristics were also detected.

Annotating long intergenic non-coding RNAs under artificial selection during chicken domestication

Background

Numerous biological functions of long intergenic non-coding RNAs (lincRNAs) have been identified. However, the contribution of lincRNAs to the domestication process has remained elusive. Following domestication from their wild ancestors, animals display substantial changes in many phenotypic traits. Therefore, it is possible that diverse molecular drivers play important roles in this process.

Results

We analyzed 821 transcriptomes in this study and annotated 4754 lincRNA genes in the chicken genome. Our population genomic analysis indicates that 419 lincRNAs potentially evolved during artificial selection related to the domestication of chicken, while a comparative transcriptomic analysis identified 68 lincRNAs that were differentially expressed under different conditions. We also found 47 lincRNAs linked to special phenotypes.

Conclusions

Our study provides a comprehensive view of the genome-wide landscape of lincRNAs in chicken. This will promote a better understanding of the roles of lincRNAs in domestication, and the genetic mechanisms associated with the artificial selection of domestic animals.

Electronic supplementary material

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

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.

Genome-wide association study for sperm membrane integrity in frozen-thawed semen of Holstein-Friesian bulls

The aim of the study was to screen the entire bull genome to identify SNP markers and propose candidate genes potentially involved in the variation of sperm membrane integrity in Holstein-Friesian bulls. Two hundred eighty eight bulls kept in one AI center were included in the study. Each bull was genotyped for 54.001 Single Nucleotide Polymorpisms (SNP) by the Illumina BovineSNP50 BeadChip. Commercial straws of frozen-thawed semen were used for the evaluation of sperm plasma membrane integrity (SYBR-14/PI staining) and sperm mitochondrial function (JC1/PI staining). An additive model for Linear Regression Analysis was applied to estimate the effect of SNP marker for sperm membrane integrity (by the use of GoldenHelix SVS7 software). Five significant markers (encompassing 2,2 MB region located on chromosome 6) for SYBR-14/PI were found. Among them one marker-rs41570391 passed Bonferroni correction test. Within approximately 3 Mb genomic region including significant markers three candidate genes: SGMS2 (Sphingomyelin Synthase 2), TET2 (Methylcytosine dioxygenase 2) and GSTCD genes (Gluthatione S-transferase C terminal domain) were proposed as potentially involved in sperm membrane integrity in frozen-thawed semen of Holstein-Friesian bulls.

Second messengers, steroids and signaling cascades: Crosstalk in sperm development and function

Signaling cascades control numerous aspects of sperm physiology, ranging from creation to fertilization. Novel aspects of several kinases and their influence on sperm development will be discussed in the first section and cover proliferation, chromatin remodeling and morphology. In the second section, protein kinases (A, B and C) that affect sperm function and their regulation by second messengers, cyclic-AMP and phosphoinositides, as well as steroids will be featured. Key areas of integration will be presented on the topics of sperm motility, capacitation, acrosome reaction and fertilization.

Variations in Antioxidant Genes and Male Infertility

Oxidative stress and reactive oxygen species (ROS) are generated from both endogenous and environmental resources, which in turn may cause defective spermatogenesis and male infertility. Antioxidant genes, which include catalase (CAT), glutathione peroxidase (GPX), glutathione S-transferase (GST), nitric oxide synthase (NOS), nuclear factor erythroid 2-related factor 2 (NRF2), and superoxide dismutase (SOD), play important roles in spermatogenesis and normal sperm function. In this review, we discuss the association between variations in major antioxidant genes and male infertility. Numerous studies have suggested that genetic disruption or functional polymorphisms in these antioxidant genes are associated with a higher risk for male infertility, which include low sperm quality, oligoasthenoteratozoospermia, oligozoospermia, and subfertility. The synergistic effects of environmental ROS and functional polymorphisms on antioxidant genes that result in male infertility have also been reported. Therefore, variants in antioxidant genes, which independently or synergistically occur with environmental ROS, affect spermatogenesis and contribute to the occurrence of male infertility. Large cohort and multiple center-based population studies to identify new antioxidant genetic variants that increase susceptibility to male infertility as well as validate its potential as genetic markers for diagnosis and risk assessment for male infertility for precise clinical approaches are warranted.

Missense mutation in glutathione-S-transferase M1 gene is associated with sperm motility and ATP content in frozen-thawed semen of Holstein-Friesian bulls

Glutathione-S-transferase genes (GSTs) encode enzymes that are involved in detoxification and neutralization of reactive oxygen species (ROS) in male reproductive system and play protective role during spermatogenesis. The aim of the study was to evaluate whether C/G missense mutation (rs135955605) within glutathione-S-transferase M1 (GSTM1) gene is associated with selected parameters of frozen-thawed semen in 309 Holstein-Friesian bulls. Single nucleotide substitution C/G was identified by amplification of GSTM1 gene fragment followed be digestion with restriction enzyme DdeI. Bulls with GG genotype were the most frequent (67.96%), in comparison to CC (2.59%) and GC (29.45%). Significant associations were found between GSTM1 genotypes and ATP content and total sperm motility. Bulls with GG genotype had the highest values for both traits. Rare variant C of GSTM1 was associated with significant decrease of sperm motility and ATP content. Our results demonstrate that C/G missense mutation within GSTM1 gene is involved in bull sperm quality.