Identification of novel candidate genes for age at first calving in Nellore cows using a SNP chip specifically developed for Bos taurus indicus cattle

Puberty classifications in beef heifers are moderately to highly heritable and associated with candidate genes related to cyclicity and timing of puberty

Introduction: Pubertal attainment is critical to reproductive longevity in heifers. Previously, four heifer pubertal classifications were identified according to attainment of blood plasma progesterone concentrations > 1 ng/ml: 1) Early; 2) Typical; 3) Start-Stop; and 4) Non-Cycling. Early and Typical heifers initiated and maintained cyclicity, Start-Stop started and then stopped cyclicity and Non-Cycling never initiated cyclicity. Start-Stop heifers segregated into Start-Stop-Discontinuous (SSD) or Start-Stop-Start (SSS), with SSD having similar phenotypes to Non-Cycling and SSS to Typical heifers. We hypothesized that these pubertal classifications are heritable, and loci associated with pubertal classifications could be identified by genome wide association studies (GWAS). Methods: Heifers (n = 532; 2017 - 2022) genotyped on the Illumina Bovine SNP50 v2 or GGP Bovine 100K SNP panels were used for variant component estimation and GWAS. Heritability was estimated using a univariate Bayesian animal model. Results: When considering pubertal classifications: Early, Typical, SSS, SSD, and Non-Cycling, pubertal class was moderately heritable (0.38 ± 0.08). However, when heifers who initiated and maintained cyclicity were compared to those that did not cycle (Early+Typical vs. SSD+Non-Cycling) heritability was greater (0.59 ± 0.19). A GWAS did not identify single nucleotide polymorphisms (SNPs) significantly associated with pubertal classifications, indicating puberty is a polygenic trait. A candidate gene approach was used, which fitted SNPs within or nearby a set of 71 candidate genes previously associated with puberty, PCOS, cyclicity, regulation of hormone secretion, signal transduction, and methylation. Eight genes/regions were associated with pubertal classifications, and twenty-two genes/regions were associated with whether puberty was attained during the trial. Additionally, whole genome sequencing (WGS) data on 33 heifers were aligned to the reference genome (ARS-UCD1.2) to identify variants in FSHR, a gene critical to pubertal attainment. Fisher's exact test determined if FSHR SNPs segregated by pubertal classification. Two FSHR SNPs that were not on the bovine SNP panel were selected for additional genotyping and analysis, and one was associated with pubertal classifications and whether they cycled during the trial. Discussion: In summary, these pubertal classifications are moderately to highly heritable and polygenic. Consequently, genomic tools to inform selection/management of replacement heifers would be useful if informed by SNPs associated with cyclicity and early pubertal attainment.

Genome-wide association study identifies genomic regions associated with key reproductive traits in Korean Hanwoo cows

BACKGROUND

Conducting genome-wide association studies (GWAS) for reproductive traits in Hanwoo cattle, including age at first calving (AFC), calving interval (CI), gestation length (GL), and number of artificial inseminations per conception (NAIPC), is of paramount significance. These analyses provided a thorough exploration of the genetic basis of these traits, facilitating the identification of key markers for targeted trait improvement. Breeders can optimize their selection strategies, leading to more efficient and sustainable breeding programs, by incorporating genetic insights. This impact extends beyond individual traits and contributes to the overall productivity and profitability of the Hanwoo beef cattle industry. Ultimately, GWAS is essential in ensuring the long-term genetic resilience and adaptability of Hanwoo cattle populations. The primary goal of this study was to identify significant single nucleotide polymorphisms (SNPs) or quantitative trait loci (QTLs) associated with the studied reproductive traits and subsequently map the underlying genes that hold promise for trait improvement.

RESULTS

A genome-wide association study of reproductive traits identified 68 significant single nucleotide polymorphisms (SNPs) distributed across 29 Bos taurus autosomes (BTA). Among them, BTA14 exhibited the highest number of identified SNPs (25), whereas BTA6, BTA7, BTA8, BTA10, BTA13, BTA17, and BTA20 exhibited 8, 5, 5, 3, 8, 2, and 12 significant SNPs, respectively. Annotation of candidate genes within a 500 kb region surrounding the significant SNPs led to the identification of ten candidate genes relevant to age at first calving. These genes were: FANCG, UNC13B, TESK1, TLN1, and CREB3 on BTA8; FAM110B, UBXN2B, SDCBP, and TOX on BTA14; and MAP3K1 on BTA20. Additionally, APBA3, TCF12, and ZFR2, located on BTA7 and BTA10, were associated with the calving interval; PAX1, SGCD, and HAND1, located on BTA7 and BTA13, were linked to gestation length; and RBM47, UBE2K, and GPX8, located on BTA6 and BTA20, were linked to the number of artificial inseminations per conception in Hanwoo cows.

CONCLUSIONS

The findings of this study enhance our knowledge of the genetic factors that influence reproductive traits in Hanwoo cattle populations and provide a foundation for future breeding strategies focused on improving desirable traits in beef cattle. This research offers new evidence and insights into the genetic variants and genome regions associated with reproductive traits and contributes valuable information to guide future efforts in cattle breeding.

Identification of quantitative trait loci and associated candidate genes for pregnancy success in Angus-Brahman crossbred heifers

BACKGROUND

In beef cattle, more than 50% of the energy input to produce a unit of beef is consumed by the female that produced the calf. Development of genomic tools to identify females with high genetic merit for reproductive function could increase the profitability and sustainability of beef production.

RESULTS

Genome-wide association studies (GWAS) were performed using a single-step genomic best linear unbiased prediction approach on pregnancy outcome traits from a population of Angus-Brahman crossbred heifers. Furthermore, a validation GWAS was performed using data from another farm. Heifers were genotyped with the Bovine GGP F250 array that contains 221,077 SNPs. In the discovery population, heifers were bred in winter breeding seasons involving a single round of timed artificial insemination (AI) followed by natural mating for 3 months. Two phenotypes were analyzed: pregnancy outcome to first-service AI (PAI; n = 1,481) and pregnancy status at the end of the breeding season (PEBS; n = 1,725). The heritability was estimated as 0.149 and 0.122 for PAI and PEBS, respectively. In the PAI model, one quantitative trait locus (QTL), located between 52.3 and 52.5 Mb on BTA7, explained about 3% of the genetic variation, in a region containing a cluster of γ-protocadherin genes and SLC25A2. Other QTLs explaining between 0.5% and 1% of the genetic variation were found on BTA12 and 25. In the PEBS model, a large QTL on BTA7 was synonymous with the QTL for PAI, with minor QTLs located on BTA5, 9, 10, 11, 19, and 20. The validation population for pregnancy status at the end of the breeding season were Angus-Brahman crossbred heifers bred by natural mating. In concordance with the discovery population, the large QTL on BTA7 and QTLs on BTA10 and 12 were identified.

CONCLUSIONS

In summary, QTLs and candidate SNPs identified were associated with pregnancy outcomes in beef heifers, including a large QTL associated with a group of protocadherin genes. Confirmation of these associations with larger populations could lead to the development of genomic predictions of reproductive function in beef cattle. Moreover, additional research is warranted to study the function of candidate genes associated with QTLs.

Comparative expression analysis of water buffalo (Bubalus bubalis) to identify genes associated with economically important traits

The milk, meat, skins, and draft power of domestic water buffalo (Bubalus bubalis) provide substantial contributions to the global agricultural economy. The world's water buffalo population is primarily found in Asia, and the buffalo supports more people per capita than any other livestock species. For evaluating the workflow, output rate, and completeness of transcriptome assemblies within and between reference-free (RF) de novo transcriptome and reference-based (RB) datasets, abundant bioinformatics studies have been carried out to date. However, comprehensive documentation of the degree of consistency and variability of the data produced by comparing gene expression levels using these two separate techniques is lacking. In the present study, we assessed the variations in the number of differentially expressed genes (DEGs) attained with RF and RB approaches. In light of this, we conducted a study to identify, annotate, and analyze the genes associated with four economically important traits of buffalo, viz., milk volume, age at first calving, post-partum cyclicity, and feed conversion efficiency. A total of 14,201 and 279 DEGs were identified in RF and RB assemblies. Gene ontology (GO) terms associated with the identified genes were allocated to traits under study. Identified genes improve the knowledge of the underlying mechanism of trait expression in water buffalo which may support improved breeding plans for higher productivity. The empirical findings of this study using RNA-seq data-based assembly may improve the understanding of genetic diversity in relation to buffalo productivity and provide important contributions to answer biological issues regarding the transcriptome of non-model organisms.

Major locus on ECA18 influences effectiveness of GonaCon vaccine in feral horses

Contraceptive vaccines are used to reduce birth rates in wild and feral animal populations. While the immunocontraceptive GonaCon-Equine has proven effective in reducing fertility among female feral horses, there is individual variation in the duration of infertility following treatment. To identify genetic factors influencing the effectiveness of GonaCon-Equine, we conducted a genome-wide association study of 88 mares from a feral population genotyped using the Illumina GGP Equine 70k SNP array. Contraceptive treatment schedules and long-term foaling rates have been recorded for each individual. We used mixed linear models to control for relatedness among mares. We found a significant association (p < 5 ×10^-8) with a locus on equine chromosome 18. The most likely candidate genes in this region are STAT1 and STAT4, which are both involved in immune system function. Variation in STAT function could affect the immune response to the vaccine, leading to variation in contraceptive efficacy. Additional SNPs reaching a less stringent threshold of significance (p < 5 ×10^-6) were located on other chromosomes near known immune system genes, supporting the hypothesis that variation in immunocontraceptive efficacy can be attributed to genetic variation in immune response rather than fertility genes.

Analysis of runs of homozygosity of cattle living in different climate zones

Aim of this study was to analyze the distribution and characteristics of runs of homozygosity in Bos taurus taurus and Bos taurus indicus breeds, as well as their crosses, farmed all around the world. With this aim in view, we used single-nucleotide polymorphisms (SNP) genotypes for 3,263 cattle belonging to 204 different breeds. After quality control, 23,311 SNPs were retained for the analysis. Animals were divided into seven different groups: 1) continental taurus, 2) temperate taurus, 3) temperate indicus, 4) temperate composite, 5) tropical taurus, 6) tropical indicus, and 7) tropical composite. The climatic zones were created according to the latitude of the breeds' country of origin: i) continental, latitude ≥ 45°; ii) temperate, 45°< Latitude >23.26°; iii) tropics, latitude ≤ 23.26°. Runs of homozygosity were computed as 15 SNPs spanning in at least 2 Mb; number of ROH per animal (nROH), average ROH length (meanMb), and ROH-based inbreeding coefficients (FROH) were also computed. Temperate indicus showed the largest nROH, whereas Temperate taurus the lowest value. Moreover, the largest meanMb was observed for Temperate taurus, whereas the lowest value for Tropics indicus. Temperate indicus breeds showed the largest FROH values. Genes mapped in the identified ROH were reported to be associated with the environmental adaptation, disease resistance, coat color determinism, and production traits. Results of the present study confirmed that runs of homozygosity could be used to identify genomic signatures due to both artificial and natural selection.

Gene expression in the amygdala and hippocampus of cyclic and acyclic gilts

Age at first estrus is the earliest phenotypic indicator of future reproductive success of gilts. Prebreeding anestrus is a major reason for reproductive failure leading to culling of replacement gilts. The two types of prebreeding anestrus are delay in attaining puberty (prepubertal anestrus, PPA) and silent ovulation (behavioral anestrus, BA). Neural tissues such as amygdala and hippocampus play a major role in regulating sexual behavior, social interactions, and receptivity to males. Differences in gene expression in the amygdala and hippocampus of gilts were analyzed in three comparisons: 1) PPA cases and cyclic controls at follicular phase of estrous cycle, 2) BA cases and cyclic controls at luteal phase of estrous cycle, and 3) gilts at different stages of the ovarian cycle (cyclic gilts at follicular phase and luteal phase of estrous cycle) to gain functional understanding of how these rarely studied tissues may differ between pubertal phenotypes and different stages of the estrous cycle of gilts. Differentially expressed genes (DEG) between PPA and BA cases and their respective cyclic controls were involved in neurological and behavioral disorders as well as nervous system functions that could directly or indirectly involved in development of behaviors related to estrus. The comparison between cyclic follicular and luteal phase control gilts identified the greatest number of DEG in the hippocampus and amygdala. These DEG were involved in adult neurogenesis and neural synapse (e.g., GABAergic, dopamine, cholinergic), suggesting that these tissues undergo structural changes and synaptic plasticity in gilts. This is the first report to demonstrate that the stage of estrous cycle is associated with dynamic changes in gene expression within porcine hippocampus and amygdala and indicates a role of gonadal steroids in regulating their biology.

The Role of microRNAs in the Mammary Gland Development, Health, and Function of Cattle, Goats, and Sheep

Milk is an integral and therefore complex structural element of mammalian nutrition. Therefore, it is simple to conclude that lactation, the process of producing milk, is as complex as the mammary gland, the organ responsible for this biochemical activity. Nutrition, genetics, epigenetics, disease pathogens, climatic conditions, and other environmental variables all impact breast productivity. In the last decade, the number of studies devoted to epigenetics has increased dramatically. Reports are increasingly describing the direct participation of microRNAs (miRNAs), small noncoding RNAs that regulate gene expression post-transcriptionally, in the regulation of mammary gland development and function. This paper presents a summary of the current state of knowledge about the roles of miRNAs in mammary gland development, health, and functions, particularly during lactation. The significance of miRNAs in signaling pathways, cellular proliferation, and the lipid metabolism in agricultural ruminants, which are crucial in light of their role in the nutrition of humans as consumers of dairy products, is discussed.