Genetic and Genomic Analyses of Service Sire Effect on Female Reproductive Traits in Holstein Cattle

Association of Genes TRH, PRL and PRLR with Milk Performance, Reproductive Traits and Heat Stress Response in Dairy Cattle

In our previous study, we found that changes in plasma prolactin (PRL) concentration were significantly associated with heat stress in dairy cows, and that PRL plays an important role in milk performance. Microarray sequencing revealed that thyrotropin releasing hormone (TRH) and prolactin receptor (PRLR), two genes important for PRL expression or function, may affect milk performance, reproduction, and heat stress response in dairy cattle. In this study, we further validated the genetic effects of the three genes in Chinese Holsteins. The potential variants within the three genes were first detected in 70 Chinese Holstein bulls and then screened in 1152 Chinese Holstein cows using the KASP (Kompetitive allele-specific PCR) method. In total, 42 variants were identified. Further, 13 SNPs were retained for KASP genotyping, including 8 in TRH, 3 in PRL, and 2 in PRLR. Using SNP-based association analyses, the multiple significant (p < 0.05) associations of these 13 SNPs with milk performance, reproduction, and heat stress response traits were found in the Holstein population. Furthermore, linkage disequilibrium analysis found a haplotype block in each of the TRH and PRL genes. Haplotype-based association analyses showed that haplotype blocks were also significantly (p < 0.05) associated with milk performance, reproduction, and heat stress response traits. Collectively, our results identified the genetic associations of TRH, PRL, and PRLR with milk performance, reproduction, and heat stress response traits in dairy cows, and found the important roles of SNP g.55888602A/C and g.55885455A/G in TRH in all traits, providing important molecular markers for genetic selection of high-yielding dairy cows.

Weighted single step GWAS reveals genomic regions associated with economic traits in Murrah buffaloes

The objective of the present study was to identify genomic regions influencing economic traits in Murrah buffaloes using weighted single step Genome Wide Association Analysis (WssGWAS). Data on 2000 animals, out of which 120 were genotyped using a double digest Restriction site Associated DNA (ddRAD) sequencing approach. The phenotypic data were collected from NDRI, India, on growth traits, viz., body weight at 6M (month), 12M, 18M and 24M, production traits like 305D (day) milk yield, lactation length (LL) and dry period (DP) and reproduction traits like age at first calving (AFC), calving interval (CI) and first service period (FSP). The biallelic genotypic data consisted of 49353 markers post-quality check. The heritability estimates were moderate to high, low to moderate, low for growth, production, reproduction traits, respectively. Important genomic regions explaining more than 0.5% of the total additive genetic variance explained by 30 adjacent SNPs were selected for further analysis of candidate genes. In this study, 105 genomic regions were associated with growth, 35 genomic regions with production and 42 window regions with reproduction traits. Different candidate genes were identified in these genomic regions, of which important are OSBPL8, NAP1L1 for growth, CNTNAP2 for production and ILDR2, TADA1 and POGK for reproduction traits.

Effect of Dam Body Conformations on Birth Traits of Calves in Chinese Holsteins

The objective of this study was to explore the effect of dam body conformations on birth traits including stillbirth, dystocia, gestation length and birth weight of Chinese Holstein calves and to provide a reference for improving cow reproductive performance. We collected phenotype data on 20 conformation traits of Chinese Holstein cows and analyzed the impact of dam conformation trait linear scores on stillbirth, dystocia, gestation length and calf birth weight. The feet angle, set of rear legs, fore udder attachment and rear attachment height traits of the dairy cows significantly affected the risk of stillbirth. The risk of dystocia decreases with the increase in stature and pin width. The bone quality of dairy cows had a significant positive correlation with gestation length. There was a significant positive correlation between fore udder attachment and calf weight at birth. The birth weight of calves from cows with high body conformation traits was significantly higher than that of calves with a low composite score. These results suggest that improving the body conformation traits, especially the selection of mammary system and body shape total score, will be beneficial in improving the reproductive performance of dairy cows.

Comparative analyses of copy number variations between swamp buffaloes and river buffaloes

As an important source of genomic variation, copy number variation (CNV) contributes to environmental adaptation in worldwide buffaloes. Despite this importance, CNV divergence between swamp buffaloes and river buffaloes has not been studied previously. Here, we report 21 152 CNV regions (CNVRs) in 141 buffaloes of 20 breeds detected through multiple CNV calling strategies. Only 248 CNVRs were shared between river buffalo and swamp buffalo, reflecting great variation of CNVRs between the two subspecies. Population structure analysis based on CNVs successfully separated the two buffalo subspecies. We further assessed CNV divergence by calculating F_ST for genome-wide CNVs. Totally, we identified 110 significantly divergent CNV segments and 44 putatively selected genes between river buffaloes and swamp buffaloes. In particular, LALBA, a key gene controlling milk production in cattle, presented a highly differentiated CNV in the promoter region, which makes it a strong functional candidate gene for differences between swamp buffaloes and river buffaloes in traits related to milk production. Our study provides useful information of CNVs in buffaloes, which may help explain the genetic differences between the two subspecies.

Genetic parameters for various semen production and quality traits and indicators of male and female reproductive performance in Nellore cattle

BACKGROUND

Given the economic relevance of fertility and reproductive traits for the beef cattle industry, investigating their genetic background and developing effective breeding strategies are paramount. Considering their late and sex-dependent phenotypic expression, genomic information can contribute to speed up the rates of genetic progress per year. In this context, the main objectives of this study were to estimate variance components and genetic parameters, including heritability and genetic correlations, for fertility, female precocity, and semen production and quality (andrological attributes) traits in Nellore cattle incorporating genomic information.

RESULTS

The heritability estimates of semen quality traits were low-to-moderate, while moderate-to-high estimates were observed for semen morphological traits. The heritability of semen defects ranged from low (0.04 for minor semen defects) to moderate (0.30 for total semen defects). For seminal aspect (SMN_ASPC) and bull reproductive fitness (BULL_FIT), low (0.19) and high (0.69) heritabilities were observed, respectively. The heritability estimates for female reproductive traits ranged from 0.16 to 0.39 for rebreeding of precocious females (REBA) and probability of pregnancy at 14 months (PP14), respectively. Semen quality traits were highly genetically correlated among themselves. Moderate-to-high genetic correlations were observed between the ability to remain productive in the herd until four years of age (stayability; STAY) and the other reproductive traits, indicating that selection for female reproductive performance will indirectly contribute to increasing fertility rates. High genetic correlations between BULL_FIT and female reproductive traits related to precocity (REBA and PP14) and STAY were observed. The genetic correlations between semen quality and spermatic morphology with female reproductive traits ranged from -0.22 (REBA and scrotal circumference) to 0.48 (REBA and sperm vigor). In addition, the genetic correlations between REBA with semen quality traits ranged from -0.23 to 0.48, and with the spermatic morphology traits it ranged from -0.22 to 0.19.

CONCLUSIONS

All male and female fertility and reproduction traits evaluated are heritable and can be improved through direct genetic or genomic selection. Selection for better sperm quality will positively influence the fertility and precocity of Nellore females. The findings of this study will serve as background information for designing breeding programs for genetically improving semen production and quality and reproductive performance in Nellore cattle.

Genomic prediction of service sire effect on female reproductive performance in Holstein cattle: A comparison between different methods, validation population and marker densities

Reproductive traits of dairy cattle are bound to the actual efficiency of farm operation, which therefore show great economic importance. Among them, some traits were deemed to be simultaneously affected by service sire and mating cow. Service sires are proved to play an important role in reproduction process of cows. However, limited study explored the genetic effect of service sire (GESS), let alone the genomic prediction of this effect. In the present study, 2244 genotyped bulls together with phenotypic records were used to predict the GESS on conception rate, 56-day non-return rate, calving ease, stillbirth and gestation length. The feasibilities of multi-step genomic best linear unbiased predictor (msGBLUP) and single-step genomic best linear unbiased predictor (ssGBLUP) were investigated under different scenarios, that is, different marker densities and validation population. The predictive accuracies and unbiasedness for GESS ranged from 0.159 to 0.647 and from 0.202 to 2.018, respectively, when validated on young bulls, while the accuracies and unbiasedness ranged from 0.409 to 0.802 and 0.333 to 1.146 when validated on random split data sets. It is feasible to predict GESS on reproductive traits by using a linear mixed model and genomic data, and high-density marker panel had limited contribution to the prediction. This research investigated the potential factors that influence the genomic prediction of GESS on reproductive traits and indicated the possibility of genomic selection on GESS, both in ideal and practical circumstances.

Weighted single-step GWAS and RNA sequencing reveals key candidate genes associated with physiological indicators of heat stress in Holstein cattle

Background

The study of molecular processes regulating heat stress response in dairy cattle is paramount for developing mitigation strategies to improve heat tolerance and animal welfare. Therefore, we aimed to identify quantitative trait loci (QTL) regions associated with three physiological indicators of heat stress response in Holstein cattle, including rectal temperature (RT), respiration rate score (RS), and drooling score (DS). We estimated genetic parameters for all three traits. Subsequently, a weighted single-step genome-wide association study (WssGWAS) was performed based on 3200 genotypes, 151,486 phenotypic records, and 38,101 animals in the pedigree file. The candidate genes located within the identified QTL regions were further investigated through RNA sequencing (RNA-seq) analyses of blood samples for four cows collected in April (non-heat stress group) and four cows collected in July (heat stress group).

Results

The heritability estimates for RT, RS, and DS were 0.06, 0.04, and 0.03, respectively. Fourteen, 19, and 20 genomic regions explained 2.94%, 3.74%, and 4.01% of the total additive genetic variance of RT, RS, and DS, respectively. Most of these genomic regions are located in the Bos taurus autosome (BTA) BTA3, BTA6, BTA8, BTA12, BTA14, BTA21, and BTA24. No genomic regions overlapped between the three indicators of heat stress, indicating the polygenic nature of heat tolerance and the complementary mechanisms involved in heat stress response. For the RNA-seq analyses, 2627 genes were significantly upregulated and 369 downregulated in the heat stress group in comparison to the control group. When integrating the WssGWAS, RNA-seq results, and existing literature, the key candidate genes associated with physiological indicators of heat stress in Holstein cattle are: PMAIP1, SBK1, TMEM33, GATB, CHORDC1, RTN4IP1, and BTBD7.

Conclusions

Physiological indicators of heat stress are heritable and can be improved through direct selection. Fifty-three QTL regions associated with heat stress indicators confirm the polygenic nature and complex genetic determinism of heat tolerance in dairy cattle. The identified candidate genes will contribute for optimizing genomic evaluation models by assigning higher weights to genetic markers located in these regions as well as to the design of SNP panels containing polymorphisms located within these candidate genes.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00748-6.

Identification of Genomic Regions and Candidate Genes for Litter Traits in French Large White Pigs Using Genome-Wide Association Studies

The reproductive traits of sows are one of the important economic traits in pig production, and their performance directly affects the economic benefits of the entire pig industry. In this study, a total of 895 French Large White pigs were genotyped by GeneSeek Porcine 50K SNP Beadchip and four phenotypic traits of 1407 pigs were recorded, including total number born (TNB), number born alive (NBA), number healthy piglets (NHP) and litter weight born alive (LWB). To identify genomic regions and genes for these traits, we used two approaches: a single-locus genome-wide association study (GWAS) and a single-step GWAS (ssGWAS). Overall, a total of five SNPs and 36 genomic regions were identified by single-locus GWAS and ssGWAS, respectively. Notably, fourof all five significant SNPs were located in 10.72–11.06 Mb on chromosome 7, were also identified by ssGWAS. These regions explained the highest or second highest genetic variance in the TNB, NBA and NHP traits and harbor the protein coding gene ENSSSCG00000042180. In addition, several candidate genes associated with litter traits were identified, including JARID2, PDIA6, FLRT2 and DICER1. Overall, these novel results reflect the polygenic genetic architecture of the litter traits and provide a theoretical reference for the following implementation of molecular breeding.